Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “higher-order calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Second-order renormalized Hamiltonian of Yukawa theory

Using the renormalization group procedure for effective particles we calculate the effective Hamiltonians in the theory of a fermion field coupled to a scalar field via the Yukawa interaction. The theory is renormalized by the addition of counterterms. Necessary counterterms are determined by computing matrix elements of the effective Hamiltonian. All calculations are performed up to the second order in the expansion in powers of the coupling constant. Renormalized effective Hamiltonians are well-defined symmetric forms acting in the Fock space as opposed to the renormalized bare Hamiltonian, which is not well defined without regularization. We introduce computational techniques that should streamline higher-order calculations and may be of independent interest.

Ab initio calculations↗
DOE OSTI2026

Tracking discontinuities in parameter space

We develop a geometric framework in Feynman-parameter space to determine constraints on the sequential discontinuities of Feynman integrals. Our method is based on tracking the deformation of the integration contour as external kinematics are analytically continued. This procedure imposes powerful constraints on the analytic structure of Feynman integrals, providing crucial inputs for their bootstrap. We demonstrate the usefulness of this framework by applying it to integrals in dimensional regularization, with higher propagator powers, and to examples with non-uniform transcendental weight. The method is illustrated with several one- and two-loop calculations.

Differential and Algebraic Geometry↗
DOE OSTI2026

Heavy-quark contributions to the DIS structure functions F 4 and F 5 at NLO in the ACOT scheme

We compute the contributions of heavy quarks to the deep-inelastic scattering structure functions F 4 and F 5 at next-to-leading order of perturbative QCD in the ACOT scheme. Both analytic results including the details of the calculation as well as numerical results for the neutral and charged current cases are presented. Our study thus lays the groundwork for future measurements of these two structure functions in experiments such as the SHiP experiment.

deep inelastic scattering↗
DOE OSTI2025

Small radius inclusive jet production at the LHC through NNLO+NNLL

The study of hadronic jets and their substructure at hadronic colliders is crucial for improving our understanding of QCD, and searching for new physics. As such, there has been a significant effort to improve their theoretical description. In the small radius limit, inclusive jet production exhibits a universal factorization, enabling the resummation of logarithms which greatly stabilizes theoretical predictions. In this paper, we show how to combine a recently introduced framework for small-R resummation with the Stripper subtraction formalism for fragmentation, enabling next-to-next-to-leading order calculations of small-R inclusive jet production for a wide variety of processes at the LHC. We extract the two-loop constants for the jet functions, enabling for the first time next-to-next-to-leading logarithmic resummation matched to next-to-next-to-leading order perturbative calculation. We compare with CMS data for small-R jet production, and find that our results greatly improve the accuracy of the predictions at small-R, and stabilize the perturbative convergence and error estimates at larger R. Our approach is applicable to a wide class of jet substructure observables exhibiting similar factorization theorems, opening the door to an NNLO jet substructure program at the LHC.

Effective Field Theories↗
DOE OSTI2026

Radiative corrections to superallowed beta decays at $\mathcal{O}\left({\alpha}^2Z\right)$

We compute $\mathcal{O}\left({\alpha}^2Z\right)$ radiative corrections to superallowed β decays with a heavy-particle effective field theory that systematically describes the interactions of low-energy ultrasoft photons with nuclei. We calculate two-loop virtual and one-loop real-virtual amplitudes by reducing the Feynman integrals to a set of master integrals, which we solve analytically using a variety of techniques. These techniques can be applied to other phenomenologically interesting observables. The ultrasoft corrections can then be combined with contributions arising from the exchange of potential photons to obtain the complete $\mathcal{O}\left({\alpha}^2Z\right)$ correction to the decay rate, with resummation of large logarithms of the electron energy times the nuclear radius. We find that $\mathcal{O}\left({\alpha}^2Z\right)$ ultrasoft loops induce a relative correction to the decay rate that ranges from 0.7 ∙ 10 −3 in the decay of 10 C to 3.6 ∙ 10 −3 in the decay of 54 Co, and will thus impact the extraction of V ud at the permille level. We show that the inclusion of these corrections reduces the residual renormalization scale dependence of the decay rate to a negligible level, making missing ultrasoft perturbative corrections a subdominant source of theoretical uncertainty.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Analytic reconstruction with massive particles: one-loop amplitudes for $0\to \overline{q} qt\overline{t}H$

We present an analytic reconstruction of one-loop amplitudes for the process $0\to \overline{q} qt\overline{t}H$. Our calculation is a novel use of analytic reconstruction, retaining explicit covariance in the massive spin states through the massive spinor-helicity formalism. The analytic reconstruction relies on embedding the massive five-point kinematics in a fully massless eight-point phase space while still building a minimal ansatz directly in the five-point phase space. In order to obtain compact analytic expressions it is necessary to identify suitable partial fraction decompositions and extract common numerator factors, which we achieve through careful inspection of limits in which pairs of denominators vanish. We find that the resulting amplitudes are more numerically efficient than ones computed using automatic methods but that the gains are not as significant as in the massless case, at least at present. The method opens the door to applications at two-loop order, where numerical efficiency and improvements in the reconstruction methodology are more crucial, especially with regards to the number of free parameters in the ansatz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The structure of quark mass corrections in the gg → HH amplitude at high-energy

The leading and next-to-leading order QCD predictions for Higgs boson pair production at hadron colliders suffer from a significant mass renormalisation scheme uncertainty related to the choice of the top-quark mass. The functional dependence of the result on the value of the intermediate quark mass can be understood in the high-energy limit using the Method of Regions and the tools of Soft-Collinear Effective Theory. In this work, we study the origin of the sizeable logarithmic mass corrections in the gg → HH amplitudes at leading and next-to-leading power in the limit s, |t|, |u| ≫ ${m}_t^2$ ≫ ${m}_H^2$. We argue that the mass corrections follow a predictable factorised pattern that can be exploited to simplify their computation. We present results for the leading power leading logarithmic corrections, our analysis leads to a significant reduction in the theoretical uncertainty of the double Higgs production amplitudes at scattering energies ≳ 1 TeV due to the top-quark mass scheme.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

SDA: a symbolic differential algebra package in C++

Truncated Power Series Algebra (TPSA), or Differential Algebra (DA), is a well-established tool in accelerator physics, commonly used for generating high-order maps of dynamic systems, as well as in symplectic tracking, normal form analysis, verified integration, optimization, and fast multipole methods. This package is the first to perform symbolic DA computations, enabling traceability of initial condition contributions and runtime reduction for repeated DA calculations, potentially expanding DA’s applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The Complexation Properties of Self-Defensive Microgel-Modified Antimicrobial Surfaces

The complexation of cationic antimicrobials with polyanionic microgels on a biomaterial surface can render that surface self-defensive against bacteria by killing those bacteria which physically contact the antimicrobial-loaded microgels. This killing has been attributed to the contact-driven transfer of antimicrobial from a microgel to a challenging bacterium, though much remains unknown about this process. Here, in this study, we use a combination of experiments and computational modeling to identify key aspects of the complexation phenomena which influence the self-defensive properties. We synthesize poly(acrylic acid) (PAA) microgels (∼2–5 μm diameter) via membrane emulsification and electrostatically deposit them onto polycaprolactone (PCL) coupons or onto glass to form a discontinuous submonolayer. Subsequent microgel loading with colistin or with Sub5 antimicrobial peptide (AMP) causes microgel deswelling. Under physiological conditions Sub5 remains stably sequestered whereas colistin is quickly released. Coarse-grained molecular dynamics (CGMD) simulations confirm stronger Sub5/PAA complexation. CGMD calculations also indicate that Sub5 forms dimers and higher-order structures, a prediction confirmed experimentally by Small-Angle X-ray Scattering (SAXS). Supramolecular structure entropically enhances the complexation strength because of enhanced counterion release per complexation event, and this finding can help identify other antimicrobials well suited for such a nonelutive yet self-defensive strategy. CGMD simulations also show that Sub5 has a higher complexation strength with the Staphylococcus aureus membrane than it does with PAA, confirming that there is a thermodynamic driving force for antimicrobial transfer. Such self-defensive surfaces significantly reduce S. aureus colonization (over 90% reduction relative to unmodified controls) in an in vitro hematogenous contamination model and remain cyto-compatible as evidenced by mesenchymal stem cell spreading and proliferation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Field redefinitions can be nonlocal

We revisit the lore establishing the allowed space of field redefinitions and show that there are essentially no restrictions. Our conclusions hold to all orders in perturbation theory and for any dispersion relation. Field redefinitions can be nonlocal, symmetry breaking, or in certain cases have explicit dependence on spacetime. We address field redefinitions that can be resummed into the propagator, which demonstrates how to perform perturbative calculations away from the minimum in field space. Field redefinitions are used to derive higher-order Schwinger-Dyson equations, which imply multiparticle soft theorems. Non-standard field redefinitions are showcased using both relativistic and nonrelativistic examples.

effective field theories↗
DOE OSTI2025

QCD splitting functions beyond kinematical limits

In this talk, I will present a new decomposition of QCD splitting functions, carried out systematically up to second order in the strong coupling. The core idea is to separate the splitting functions into two components: scalar dipole radiator functions and pure splitting remainders. While these pieces include terms that are formally subleading in the soft or collinear limits, they turn out to be universal and well understood. Unlike conventional approaches, our construction does not rely on any soft or collinear approximations. This framework not only offers a clearer understanding of the structure of higher-order QCD corrections but also paves the way for more robust subtraction schemes in precision calculations.

Knobbe, Max [Fermilab]↗
DOE OSTI2025

First measurement of symmetric cumulants of hexagonal flow harmonics in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

Correlations between event-by-event fluctuations of anisotropic flow harmonics are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair of 5.02 TeV, as recorded by the ALICE detector at the LHC. This study presents correlations up to the hexagonal flow harmonic 𝑣 6 , which was measured for the first time. The magnitudes of these higher-order correlations are found to vary as a function of collision centrality and harmonic order. These measurements are compared to viscous hydrodynamic model calculations with EKRT initial conditions and to the iEBE-VISHNU model with T R ⁢ENTo initial conditions. The observed discrepancies between the data and the model calculations vary depending on the harmonic combinations. Due to the sensitivity of model parameters estimated with Bayesian analyses to these higher-order observables, the results presented in this work provide new and independent constraints on the initial conditions and transport properties in theoretical models used to describe the system created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
DOE OSTI2026

Les Houches study on inclusive jet production at NNLO+NNLL

Jet production at the LHC is a powerful probe of QCD, making it ideal for precision tests and determinations of QCD parameters such as parton distribution functions and the strong coupling constant. To make the most of the abundant jet production data collected at the LHC, precise calculations are required. While state-of-the-art calculations reach next-to-next-to-leading order (NNLO) QCD accuracy, a critical assessment of the remaining uncertainties arising from non-perturbative effects and missing higher orders remains crucial for correctly interpreting comparisons between theory and data. Scale variation is nearly always used to determine effects from missing higher orders. In this article, we reassess this method in the context of inclusive jet production by performing NNLO QCD calculations supplemented by small-jet-radius resummation through next-to-next-to-leading-logarithmic accuracy (NNLL). We find that NNLL resummation can have an appreciable impact on the scale uncertainty for inclusive jet cross sections, and, for some scale choices, can lead to sizeable shifts of the central cross section. We conclude that scale variations in fixed-order and resummed calculations can drastically underestimate the impact of higher orders for commonly used jet radius parameters, and that missing higher-order estimates obtained via scale variations should be considered unreliable. Our findings add further evidence to the importance of going beyond scale variations in jet and jet substructure calculations.

Lee, Kyle↗
DOE OSTI2025

Ab initio computations of the fourth-order charge density moments of 48 Ca and 208 Pb

Neutron skins of neutron-rich nuclei connect nuclei with the matter in neutron stars. High-precision measurements of nuclear charge densities to extract higher-order moments are proposed to be sensitive to neutron radii and skin thicknesses. We investigate the charge density of 48 Ca and 208 Pb, leading candidates for such studies, with ab initio nuclear structure calculations. We find strong correlations between the fourth-order charge density moment R$^{4}_{ch}$ and the charge and neutron radii, allowing us to predict R$^{4}_{ch}$ for 48 Ca and 208 Pb. We find a substantially weaker correlation between the fourth-order charge density moment and the neutron skin, limiting the ability of high-precision electron scattering to determine the neutron skin in a model-independent manner.

Ab initio nuclear theory↗
DOE OSTI2026

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2025

EFT approach to the endpoint of muon decay-in-orbit

As upcoming experiments aim to probe muon conversion with unprecedented precision, equally precise theoretical predictions are crucial to maximize discovery potential. This applies not only to the new physics signal, muon-electron conversion, but also to its only irreducible background, muon decay-in-orbit (DIO) near the endpoint. Accurate computation of higher-order corrections in bound states is a long-standing challenge due to the difficulty of systematically organizing contributions. In previous work, we developed an Effective Field Theory framework to address this issue and applied it to muon conversion. Here, we extend this approach to the DIO endpoint, a more complex problem due to the presence of a neutrino-antineutrino pair in the final state. We present the most precise prediction to date of the background spectrum relevant for future muon conversion searches, achieving next-to-leading logarithmic prime accuracy for QED corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Lattice anharmonicity effects in fluorite oxide single crystals and anomalous increase in phonon lifetime in ceria at elevated temperature

Here, we investigate the temperature dependence of the frequency and linewidth of the triply degenerate T 2g zone-centered optical phonon in flux-grown ceria and hydrothermally synthesized thoria single crystals from room temperature to 1273 K using Raman spectroscopy. Both crystals exhibit an expected increase in the phonon linewidth with temperature due to enhanced phonon–phonon scattering. However, ceria displays an anomalous linewidth reduction in the temperature range of 1023–1123 K. First-principles phonon linewidth calculations considering cubic and quartic phonon interactions within temperature-independent phonon dispersion fail to describe this anomaly. A parameterization of the temperature-dependent second-order interatomic force constants based on previously reported phonon dispersion measured at room and high temperatures predicts a deviation from the monotonic linewidth increase, albeit at temperatures lower than those observed experimentally for ceria. The qualitative agreement in the trend of temperature-dependent linewidth suggests that lattice anharmonicity-induced phonon renormalization plays a role in phonon lifetime. Specifically, a change in the overlap between softened acoustic and optical branches in the dispersion curve reduces the available phonon scattering phase space of the Raman-active mode at the zone center, leading to an increased phonon lifetime within a narrow temperature interval. These findings provide insights into higher-order anharmonic interactions in ceria and thoria, motivating further investigations into the role of anharmonicity-induced phonon renormalization on phonon lifetimes at high temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗