SEARCH · Engineering Papers
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.
This lecture discusses the physics implemented by Monte Carlo event generators for hadron colliders. It details the construction of parton showers and the matching of parton showers to fixed-order calculations at higher orders in perturbative QCD. It also discusses approaches to merge calculations for a varying number of jets, the interface to the underlying event and hadronization.
The rapid deployment of computing hardware different from the traditional CPU+RAM model in data centers around the world mandates a change in the design of event generators for the Large Hadron Collider, in order to provide economically and ecologically sustainable simulations for the high-luminosity era of the LHC. Parton-level event generation is one of the most computationally demanding parts of the simulation and is therefore a prime target for improvements. We present a production-ready leading-order parton-level event generation framework capable of utilizing most modern hardware and discuss its performance in the standard candle processes of vector boson and top-quark pair production with up to five additional jets.
This technical report summarizes the main accomplishments of the research performed under the U.S. Department of Energy (DOE) Office of Science, Early Career Award (ECA), grant number DE-SC0018117, from 9/1/2017 to 8/31/2022, with Prof. Dennis Perepelitsa as the single PI. The grant proposal was to investigate the quark-gluon plasma (QGP) created in ultrarelativistic collisions of large nuclei with the ATLAS experiment at the Large Hadron Collider (LHC) and to prepare for this physics with the sPHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). In particular, the focus was on QGP regions of different sizes and shapes, including the tradtional “large” collision systems such as Pb+Pb as well as small collision systems such as p+Pb, where there has been evidence of QGP “flow” but not yet the accompanying energy loss from hard-scattered partons traversing the QGP medium. Additionally, a particular tool was the use of photon+jet events, where the photon is unaffected by the QGP and thus serves as the “control” while the jet (the result of parton fragmentation) is the “probe” of the QGP. The research was successfully carried out, resulting in over ten publications of experimental measurements in peer-reviewed journals, and the significant training of junior scientists.
Within the pseudo-PDF framework, we investigate the perturbative contributions to correlators that are used to study transverse momentum dependent parton distributions (TMDs) on the lattice. Our results contain the full perturbative corrections which arise as artifacts from performing the calculation for a Euclidean separation between the parton fields, as well as the corrections which yield the evolution equations for the TMDs. This gives a path to extraction of the genuine non-perturbative features of the three-dimensional structure of the hadrons from lattice calculations.
In recent years there has been a momentous surge of interest in the determination of the parton content of nucleons using lattice QCD. Scalar field theory is a convenient toy model to understand the relation between Euclidean correlators and Parton Distribution Functions (PDFs). In this talk we review the factorization formulae that allow the extraction of PDFs from the pseudo-PDFs and quasi-PDFs computed on the lattice. The scalar theory provides a simplified setting where the main physical properties are easily analysed.
We introduce a new type of exclusive processes for a better study of generalized parton distributions (GPDs), which we refer to as single-diffractive hard exclusive processes (SDHEPs). We advocate a two-stage framework for picturing SDHEPs based on the separation of scales, which gives a clear description both kinematically and dynamically. We examine the sensitivity of the SDHEP to the parton momentum fraction x -dependence of GPDs, and demonstrate it quantitatively with two specific processes that can be readily measured at J-PARC or AMBER using a pion beam and at JLab using a photon beam, respectively. Both processes are capable of providing enhanced sensitivity to the x -dependence, overcoming the problem of shadow GPDs, and disentangling different types of GPDs with various spin asymmetries.
In this study, we review transverse momentum dependent (TMD) parton distribution functions, their application to topical issues in high-energy physics phenomenology, and their theoretical connections with QCD resummation, evolution and factorization theorems. We illustrate the use of TMDs via examples of multi-scale problems in hadronic collisions. These include transverse momentum q T spectra of Higgs and vector bosons for low q T , and azimuthal correlations in the production of multiple jets associated with heavy bosons at large jet masses. We discuss computational tools for TMDs, and present the application of a new tool, TMD LIB , to parton density fits and parameterizations.
We present the first official release of the nCTEQ nuclear parton distribution functions with errors. The main addition to the previous nCTEQ PDFs is the introduction of PDF uncertainties based on the Hessian method. Another important addition is the inclusion of pion production data from RHIC that give us a handle on constraining the gluon PDF. This contribution summarizes our results from arXiv:1509.00792 and concentrates on the comparison with other groups providing nuclear parton distributions.
We present the first official release of the nCTEQ nuclear parton distribution functions with errors. The main addition to the previous nCTEQ PDFs is the introduction of PDF uncertainties based on the Hessian method. Another important addition is the inclusion of pion production data from RHIC that give us a handle on constraining the gluon PDF. This contribution summarizes our results from arXiv:1509.00792, and concentrates on the comparison with other groups providing nuclear parton distributions.
In scattering experiments, high-virtuality partons, i.e., quarks and gluons, initiate a series of additional parton emissions to create collimated sprays of particles known as jets. This paper presents a measurement of the Lund jet plane (LJP) of high-energy jets produced in lead-lead (PbPb) collisions and compares the results to data for proton-proton (pp) collisions. The LJP is formed by iteratively declustering the constituents of a jet into consecutive emissions and recording the relative transverse momentum ($k_\mathrm{T}$) and angle of the resulting emission with respect to its emitter. The angular distributions of two different $k_\mathrm{T}$ slices of the LJP are investigated for jets with radius parameter of 0.4 and transverse momentum in the range 200$-$1000 GeV. The PbPb (pp) data were recorded by the CMS experiment in 2018 (2017) and correspond to an integrated luminosity of 1.7 nb$^{-1}$ (301 pb$^{-1}$) at a nucleon-nucleon center-of-mass energy of 5.02 TeV. The measurement was designed to test whether the earliest jet emissions are produced before the formation of the quark-gluon plasma (QGP) in PbPb collisions. Within the experimental uncertainties, no significant difference is observed between the angular distribution of high-$k_\mathrm{T}$ emissions in \pp and PbPb collisions, which is consistent with these emissions occurring early in the jet evolution, before substantial interaction with the QGP.
We introduce a technique for the next-to-leading order accurate simulation of $e^+e^-\to W^+W^-b\bar{b}$ that respects the resonant nature of the process above and near the top-quark pair production threshold. The parton-shower evolution, infrared subtraction and NLO matching account in particular for finite width effects beyond the Breit-Wigner structure considered in resonance-aware approaches. We present first phenomenological results relevant to a potential future electron-positron collider and provide a publicly available simulator based on the ALARIC parton shower and the SHERPA event generator.
We present the first lattice QCD calculations of the isovector helicity transverse momentum-dependent parton distribution function (TMDPDF) and the flavor-dependent unpolarized TMDPDFs for up and down quarks in the proton. Our computations utilize domain-wall fermion discretization with physical quark masses. Employing Coulomb-gauge-fixed quark correlation functions within the large-momentum effective theory framework, we access nonperturbative transverse quark separations 𝑏 𝑇 up to approximately 1 fm, corresponding to transverse momenta as low as 200 MeV. Based on the quasi-TMD factorization theorem, we construct renormalization-group–invariant ratios that are equal to the corresponding light-cone TMDPDF ratios. At moderate 𝑥, our results reveal that the isovector helicity and unpolarized TMDPDFs exhibit nearly identical transverse structure up to a normalization factor, and the unpolarized distributions display only mild flavor dependence. These findings not only support key trends observed in recent global analyses but also provide robust, nonperturbative constraints that can distinguish between different parameterizations. This Letter establishes a first-principles benchmark for TMDPDFs, offering valuable input for ongoing and future experimental efforts to map the proton’s three-dimensional structure.
The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.
Parton distribution functions (PDFs) at large x are poorly constrained by high-energy experimental data, but extremely important for probing physics beyond standard model at colliders. We study the calculation of PDFs at large-x through large-momentum P z expansion of the lattice quasi PDFs. Similar to deep-inelastic scattering, there are two distinct perturbative scales in the threshold limit where the matching coefficient can be factorized into a space-like jet function at scale P z |1 – y| and a pair of heavy-light Sudakov form factors at scale P z . The matching formula allows us to derive a full renormalization group resummation of large threshold logarithms, and the result is consistent with the known calculation to the next-to-next to leading order (NNLO). This paves the way for direct large-x PDFs calculations in lattice QCD. As by-products, we find that the space-like jet function is related to a time-like version calculated previously through analytic continuation, and the heavy-light Sudakov form factor, calculated here to NNLO, is a universal object appearing as well in the large momentum expansion of quasi transverse-momentum-dependent PDFs and quasi wave-function amplitudes.
We present the calculation of all non-planar master integrals that are needed to describe production of two off-shell vector bosons in collisions of two massless partons through NNLO in perturbative QCD. The integrals are computed analytically using differential equations in external kinematic variables and expressed in terms of Goncharov polylogarithms. These results provide the last missing ingredient needed for the computation of two-loop amplitudes that describe the production of two gauge bosons with different invariant masses in hadron collisions.
We present a simple parton-shower model that replaces the explicit angular ordering of the coherent branching formalism with a differentially accurate simulation of soft-gluon radiation by means of a non-trivial dependence of the splitting functions on azimuthal angles. We introduce a global kinematics mapping and provide an analytic proof that it satisfies the criteria for next-to leading logarithmic accuracy. In the new algorithm, initial and final state evolution are treated on the same footing. We provide an implementation for final-state evolution in the numerical code ALARIC and present a first comparison to experimental data.
We review recent works on the modelling of Generalised Parton Distributions within the Dyson-Schwinger formalism. We highlight how covariant computations, using the impulse approximation, allows one to fulfil most of the theoretical constraints of the GPDs. A specific attention is brought to chiral properties and especially the so-called soft pion theorem, and its link with the Axial-Vector Ward-Takahashi identity. The limitation of the impulse approximation are also explained. Beyond impulse approximation computations are reviewed in the forward case. Finally, we stress the advantages of the overlap of lightcone wave functions, and possible ways to construct covariant GPD models within this framework, in a two-body approximation
The cross section for deep inelastic lepton-proton scattering (DIS) ℓp→ℓ'p'X includes a diffractive deep inelastic (DDIS) contribution ℓp→ℓ'p'X, in which the proton remains intact with a large longitudinal momentum fraction x F greater than 0.9 and small transverse momentum. The DDIS events, which can be identified with Pomeron exchange in the t-channel, account for approximately 10% of all of the DIS events. Thus, when one measures DIS, one automatically includes the leading-twist Bjorken-scaling DDIS events as a contribution to the DIS cross section, whether or not the final-state proton p' is detected. In such events, the missing momentum fraction x p ' ~0.9 carried by the final-state proton p' in the DDIS events could be misidentified with the light-front momentum fraction carried by sea quarks or gluons in the protons' Fock structure. As we shall show in this article, the underlying QCD Pomeron-exchange amplitude which produces the DDIS events does not obey the operator product expansion nor satisfy momentum sum rules. Thus we conclude that the quark and gluon distributions measured in DIS experiments will be misidentified, unless the measurements explicitly exclude the DDIS events and that a correct determination of the parton distribution functions (PDFs) derived from the DIS data requires the explicit subtraction of the DDIS contribution from the full DIS cross section.