Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NLO Computations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Back-to-back inclusive dijets in DIS at small x: gluon Weizsäcker-Williams distribution at NLO

In [1], we performed the first complete computation of the back-to-back inclusive dijet cross-section in Deeply Inelastic Scattering (DIS) at small x Bj to next-to-leading order (NLO) in the Color Glass Condensate effective field theory (CGC EFT). We demonstrate here that for dijets with relative transverse momentum P ⊥ and transverse momentum imbalance q ⊥ , to leading power in q ⊥ /P ⊥ , the cross-section for longitudinally polarized photons can be fully factorized into the product of a perturbative impact factor and the non-perturbative Weizsäcker-Williams (WW) transverse momentum dependent (TMD) gluon distribution to NLO accuracy. The impact factor can further be expressed as the product of a universal soft factor which resums Sudakov double and single logs in P ⊥ /q ⊥ and a coefficient function given by a remarkably compact analytic expression. We show that in the CGC EFT the WW TMD satisfies a kinematically constrained JIMWLK renormalization group evolution in rapidity. This factorization formula is valid to all orders in Q s /q ⊥ for q ⊥ , Q s << P ⊥ , where Q s is the semi-hard saturation scale that grows with rapidity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Non-singlet quark helicity $\mathrm{PDF}$s of the nucleon from pseudo-distributions

The non-singlet helicity quark parton distribution functions (PDFs) of the nucleon are determined from lattice QCD, by jointly leveraging pseudo-distributions and the distillation spatial smearing paradigm. A Lorentz decomposition of appropriately isolated space-like matrix elements reveals pseudo-distributions that contain information on the leading-twist helicity PDFs, as well as an invariant amplitude that induces an additional z 2 contamination of the leading-twist signal. An analysis of the short-distance behavior of the space-like matrix elements using matching coefficients computed to next-to-leading order (NLO) exposes the desired PDF up to this additional z 2 contamination. Due to the non-conservation of the axial current, we elect to isolate the helicity PDFs normalized by the nucleon axial charge at the same scale μ 2 . The leading-twist helicity PDFs as well as several sources of systematic error, including higher-twist effects, discretization errors, and the aforementioned z 2 contaminating amplitude are jointly determined by characterizing the computed pseudo-distribution in a basis of Jacobi polynomials. The Akaike Information Criterion is exploited to effectively average over distinct model parameterizations and cuts on the pseudo-distribution. Encouraging agreement is observed with recent global analyses of each non-singlet quark helicity PDF, notably a rather small non-singlet anti-quark helicity PDF for all quark momentum fractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Computationally-driven discovery of second harmonic generation in EuBa 3 (B 3 O 6 ) 3 through inversion symmetry breaking

Nonlinear optical (NLO) crystals with superior properties are significant for advancing laser technologies and applications. Introducing rare earth metals to borates is a promising and effective way to modify the electronic structure of a crystal to improve its optical properties in the visible and ultraviolet range. In this work, we computationally discover inversion symmetry breaking in EuBa 3 (B 3 O 6 ) 3 , which was previously identified as centric, and demonstrate noncentrosymmetry via synthesizing single crystals for the first time by the floating zone method. We determine the correct space group to be P 6¯. The material has a large direct bandgap of 5.56 eV and is transparent down to 250 nm. The complete anisotropic linear and nonlinear optical properties were also investigated with a d 11 of ∼0.52 pm/V for optical second harmonic generation. Further, it is Type I and Type II phase matchable. This work suggests that rare earth metal borates are an excellent crystal family for exploring future deep ultraviolet (DUV) NLO crystals. It also highlights how first principles computations combined with experiments can be used to identify noncentrosymmetric materials that have been wrongly assigned to be centrosymmetric.

He, Jingyang (ORCID:0000000238863659)↗
DOE OSTI2023

MgSiP 2 : An Infrared Nonlinear Optical Crystal with a Large Non‐Resonant Phase‐Matchable Second Harmonic Coefficient and High Laser Damage Threshold

Superior infrared nonlinear optical (NLO) crystals are in urgent demand in the development of lasers and optical technologies for communications and computing. The critical challenge is to find a crystal with large non‐resonant phase‐matchable NLO coefficients and high laser damage threshold (LDTs) simultaneously, which however scale inversely. This work reports such a material, MgSiP 2, that exhibits a large second harmonic generation (SHG) coefficient of d 14 ≈ d 36 = 89 ± 5 pm V −1 at 1550 nm fundamental wavelength, surpassing the commercial NLO crystals AgGaS 2 , AgGaSe 2 , and ZnGeP 2 . First principles theory reveals the polarizability and geometric arrangement of the [SiP 4 ] tetrahedral units as the origin of this large nonlinear response. Remarkably, it also exhibits a high LDT value of 684 GW cm −2 , which is six times larger than ZnGeP 2 and three times larger than CdSiP 2 . It has a wide transparency window of 0.53–10.35 µm, allowing broadband tunability. Further, it is Type I and Type II phase‐matchable with large effective SHG coefficients of d eff,I ≈80.2 pm V −1 and d eff,II ≈73.4 pm V −1 . The outstanding properties of MgSiP 2 make it a highly attractive candidate for optical frequency conversion in the infrared.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Higgs decays at NLO in the SMEFT

The calculation of precise predictions for Higgs decays is a necessary ingredient for determining Higgs properties at the LHC and future colliders. We compute all two- and three- body Higgs decays at next-to-leading order (NLO) in both QCD and electroweak interactions using the dimension-6 Standard Model Effective Field Theory (SMEFT). Results for four-body Higgs decays that are accurate to NLO QCD/electroweak order in the SMEFT are obtained using the narrow width approximation. Our results are contained in a flexible Monte Carlo program, NEWiSH, that is publicly available and we illustrate the impact of the NLO electroweak corrections for HL-LHC, Tera-Z, and Higgstrahlung projections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Back-to-back inclusive dijets in DIS at small x : Sudakov suppression and gluon saturation at NLO

Back-to-back dijet cross-sections in deeply inelastic scattering (DIS) at small xBj are suppressed by many-body multiple scattering and screening effects arising from gluon saturation at high parton densities. They are similarly sensitive in these kinematics to large Sudakov logarithms from soft gluon radiation. Uncovering novel physics in this DIS channel therefore requires understanding the interplay of the two phenomena. In this work, we compute the small $x_\text{Bj}$ inclusive dijet DIS cross-section in back-to-back kinematics at next-to-leading order (NLO) in the Color Glass Condensate effective field theory (CGC EFT). Our result includes, for the first time, all real and virtual NLO contributions to the impact factor. These include all Sudakov double and single logarithm contributions, as well as all other finite $\mathcal{O}(α_s)$ terms that contribute at this order. We demonstrate explicitly that resummations of small x and Sudakov logarithms can be performed simultaneously in the CGC EFT. This requires that the JIMWLK kernel for small x evolution of the Weizsäcker-Williams (WW) gluon distribution satisfies a kinematic constraint imposed by lifetime ordering of successive gluon emissions; the corresponding modifications to the kernel, corresponding to resummations of large double transverse logarithms, are precisely of the type required to stabilize JIMWLK evolution beyond leading logarithmicaccuracy. We compute the azimuthal harmonics of the NLO back-to-back distributions and show their sensitivity to both the unpolarized and linearly polarized WW gluon distributions. Finally, we discuss how TMD factorization is broken by an emergent saturation scale at small $x$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The isospin-3 three-particle K -matrix at NLO in ChPT

The three-particle K-matrix, K df,3 , is a scheme-dependent quantity that parametrizes short-range three-particle interactions in the relativistic-field-theory three-particle finite-volume formalism. In this work, we compute its value for systems of three pions at maximal isospin through next-to-leading order (NLO) in Chiral Perturbation Theory (ChPT). We compare the values to existing lattice QCD results and find that the agreement between lattice QCD data and ChPT in the first two coefficients of the threshold expansion of K df,3 is significantly improved with respect to leading order once NLO effects are incorporated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

$t\overline{t}b\overline{b}$ at NLO precision in a variable flavor number scheme

Top-quark pair production in association with two $b$-jets is computed at next-to-leading order QCD precision, including effects of the $b$-quark mass, and matched to a $t\overline{t}$+jets simulation in a variable flavor number scheme. The Monte Carlo realization of this method, called fusing, consistently embeds the four-flavor calculation in a particle-level event generator. As a first phenomenological application, we present observables relevant to the data-driven estimation of irreducible backgrounds to $t\overline{t}H$ -production.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Accelerating $\mathrm{LHC}$ event generation with simplified pilot runs and fast $\mathrm{PDF}$s

Poor computing efficiency of precision event generators for LHC physics has become a bottleneck for Monte-Carlo event simulation campaigns. We provide solutions to this problem by focusing on two major components of general-purpose event generators: The PDF evaluator and the matrix-element generator. For a typical production setup in the ATLAS experiment, we show that the two can consume about 80% of the total runtime. Using NLO simulations of pp→ℓ + ℓ - +jets and pp→$t\bar{t}$+jets as an example, we demonstrate that the computing footprint of LHAPDF and SHERPA can be reduced by factors of order 10, while maintaining the formal accuracy of the event sample. The improved codes are made publicly available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Krylov spaces for truncated spectrum methodologies

We propose herein an extension of truncated spectrum methodologies, a nonperturbative numerical approach able to elucidate the low energy properties of quantum field theories. TSMs, in their various flavors, involve a division of a computational Hilbert space, H , into two parts, one part, H 1 that is “kept” for the numerical computations, and one part, H 2 , that is discarded or “truncated.” Even though H 2 is discarded, truncated spectrum methodologies will often try to incorporate the effects of H 2 in some effective way. In these terms, we propose to keep the dimension of H 1 small. We pair this choice of H 1 with a Krylov subspace iterative approach able to take into account the effects of H 2 . This iterative approach can be taken to arbitrarily high order and so offers the ability to compute quantities to arbitrary precision. In many cases it also offers the advantage of not needing an explicit UV cutoff. To compute the matrix elements that arise in the Krylov iterations, we employ a Feynman diagrammatic representation that is then evaluated with Monte Carlo techniques. Each order of the Krylov iteration is variational and is guaranteed to improve upon the previous iteration. The first Krylov iteration is akin to the next-to-leading order approach of Elias-Miró [NLO renormalization in the Hamiltonian truncation, ]. To demonstrate this approach, we focus on the ( 1 + 1 d )-dimensional ϕ 4 model and compute the bulk energy and mass gaps in both the Z 2 -broken and unbroken sectors. We estimate the critical ϕ 4 coupling in the broken phase to be g c = 0.2645 ± 0.002 . Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Developments in Performance and Portability for MadGraph5_aMC@NLO

Event generators simulate particle interactions using Monte Carlo techniques, providing the primary connection between experiment and theory in experimental high energy physics. These software packages, which are the first step in the simulation worflow of collider experiments, represent approximately 5 to 20% of the annual WLCG usage for the ATLAS and CMS experiments. With computing architectures becoming more heterogeneous, it is important to ensure that these key software frameworks can be run on future systems, large and small. In this contribution, recent progress on porting and speeding up the Madgraph5_aMC@NLO event generator on hybrid architectures, i.e. CPU with GPU accelerators, is discussed. The main focus of this work has been in the calculation of scattering amplitudes and "matrix elements", which is the computational bottleneck of an event generation application. For physics processes limited to QCD leading order, the code generation toolkit has been expanded to produce matrix element calculations using C++ vector instructions on CPUs and using CUDA for NVidia GPUs, as well as using Alpaka, Kokkos and SYCL for multiple CPU and GPU architectures. Performance is reported in terms of matrix element calculations per time on NVidia, Intel, and AMD devices. The status and outlook for the integration of this work into a production release usable by the LHC experiments, with the same functionalities and very similar user interfaces as the current Fortran version, is also described.

Valassi, Andrea↗
DOE OSTI2026

Accelerating LHC event generation with simplified pilot runs and fast PDFs

High-precision calculations are an indispensable ingredient for the success of the LHC physics programme, yet their poor computing efficiency has been a growing cause for concern, threatening to become a paralysing bottleneck in the coming years. We present solutions to eliminate the apprehension by focussing on two major components of generalpurpose Monte Carlo event generators: the evaluation of parton distribution functions, and the generation of perturbative matrix elements. We show that for the cost-driving event samples employed by the ATLAS experiment to model omnipresent, irreducible Standard Model backgrounds, such as weak boson or top-quark pair production in association with jets, these computational components dominate the overall run time by up to 80 %. We demonstrate that a reduction of the computing footprint of LHAPDF and SHERPA by factors of around 40 can be achieved for multi-leg NLO event generation.

Bothmann, Enrico [Gottingen U.]↗
DOE OSTI2025

A Demonstration of ARCANE Reweighting: Reducing the Sign Problem in the MC@NLO Generation of $e^+ e^- \rightarrow q \bar{q} + 1\, jet$ Events

Negatively weighted events, which appear in the simulation of particle collisions, significantly increase the computational requirements of collider experiments. A new technique called ARCANE reweighting has been introduced in a companion paper to tackle this problem. This paper demonstrates the technique for the next-to-leading-order generation of $e^+ e^- \rightarrow q \bar{q} + 1\, jet$ events. By redistributing the contributions of "standard" and "hard remainder" pathways in the generator that lead to the same final event, ARCANE reweighting almost completely eliminates the negative weights problem for this process. Some thoughts on implementing the technique in other scenarios are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

QED corrections to the thermal neutrino interaction rate

Abstract Motivated by precision computations of neutrino decoupling at MeV temperatures, we show how QED corrections to the thermal neutrino interaction rate can be related to the electron-positron spectral function as well as an effective$$\overline{\nu }\nu \gamma $$vertex. The spectral function is needed both in a timelike and in a spacelike domain, and for both of its physical polarization states (transverse and longitudinal with respect to spatial momentum). Incorporating an NLO evaluation of this spectral function, an estimate of the$$\overline{\nu }\nu \gamma $$vertex, and HTL resummation of scatterings mediated by soft Bose-enhancedt-channel photons, we compute the interaction rate as a function of the neutrino momentum and flavour. Effects on the – (0. . .2)% level are found, noticeably smaller than a previous estimate of a related quantity.

Physics↗
DOE OSTI2023

Matrix element corrections in the PYTHIA8 parton shower in the context of matched simulations at next-to-leading order

We discuss the role of matrix element corrections (MEC) to parton showers in the context of MC@NLO-type matchings for processes that feature unstable resonances, where MEC are liable to result in double-counting issues, and are thus generally not employed. By working with PYTHIA8 , we show that disabling all MEC is actually unnecessary in computations based on the narrow-width approximation, and we propose alternative MEC settings which, while still avoiding double counting, allow one to include hard-recoil effects in the simulations of resonance decays. We illustrate our findings by considering $t\overline{t}$ production at the LHC, and by comparing MadGraph5-aMC@NLO predictions with those of POWHEG-BOX and standalone PYTHIA8 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The three-pion K -matrix at NLO in ChPT

The three-particle K-matrix, Κ df,3 , is a scheme-dependent quantity that parametrizes short-range three-particle interactions in the relativistic-field-theory three-particle finite-volume formalism. In this work, we compute its value for systems of three pions in all isospin channels through next-to-leading order in Chiral Perturbation Theory, generalizing previous work done at maximum isospin. We obtain analytic expressions through quadratic order (or cubic order, in the case of zero isospin) in the expansion about the three-pion threshold.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Two-point correlator of twist-2 light-ray operators in N=4 SYM in BFKL approximation

We generalize local operators of the leading twist-2 of N = 4 SYM theory to the case of complex Lorentz spin j using principal series representation of sl(2,R). We give the direct computation of correlation function of two such non-local operators in the BFKL regime when j → 1. The correlator appears to have the expected conformal coordinate dependence governed by the anomalous dimension of twist-2 operator in NLO BFKL approximation predicted by Kotikov and Lipatov.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

Highly Non-Linear Optical (NLO) organic crystals and films. Electrooptical organic materials

Devices employing nonlinear optics (NLO) hold great promise for important applications in integrated optics, optical information processing and telecommunications. Properly designed organics possess outstanding optical and electrooptical properties which will substantially advance many technologies including electrooptical switching, optical amplification for communications, and parallel processing for hybrid optical computers. A brief comparison of organic and inorganic materials is given.

Mcmanus, Samuel P.↗