Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High Energy Physics”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Disentangling Centrality Bias and Final-State Effects in the Production of High-𝑝 𝑇 Neutral Pions Using Direct Photon in 𝑑+Au Collisions at $\sqrt{s_{NN}}$ = 200 GeV

Here, PHENIX presents a simultaneous measurement of the production of direct 𝛾 and 𝜋 0 in 𝑑 + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV over a 𝑝 𝑇 range of 7.5 to 18 GeV/𝑐 for different event samples selected by event activity, i.e., charged-particle multiplicity detected at forward rapidity. Direct-photon yields are used to empirically estimate the contribution of hard-scattering processes in the different event samples. Using this estimate, the average nuclear-modification factor, 𝑅$^{𝜋^0}_{dAu,EXP}$, is 0.925 ± 0.023⁢(stat) ± 0.15⁢(scale), consistent with unity for minimum-bias (MB) 𝑑+Au collisions. For event classes with low and moderate event activity, 𝑅$^{𝜋^0}_{dAu,EXP}$ is consistent with the MB value within 5% uncertainty. This result confirms that the previously observed enhancement of high-𝑝 𝑇 𝜋 0 production found in small-system collisions with low event activity is a result of a bias in interpreting event activity within the Glauber framework. In contrast, for the top 5% of events with the highest event activity, 𝑅$^{𝜋^0}_{dAu,EXP}$ is suppressed by 20% relative to the MB value with a significance of 4.5⁢𝜎, which may be due to final-state effects. This suppression corresponds to a 𝑝 𝑇 shift of 𝛿⁢𝑝 𝑇 = 0.213 ± 0.055 Gev/𝑐 at 9 Gev/𝑐.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of inclusive jet cross section and substructure in 𝑝 + 𝑝 collisions at $\sqrt{s}$ = 200 GeV

The jet cross section and jet-substructure observables in 𝑝 + 𝑝 collisions at $\sqrt{s}$ =200 GeV were measured by the PHENIX Collaboration at the Relativistic Heavy Ion Collider (RHIC). Jets are reconstructed from charged-particle tracks and electromagnetic-calorimeter clusters using the anti-𝑘 𝑡 algorithm with a jet radius of 𝑅 = 0.3 for jets with transverse momentum within 8.0 < 𝑝 𝑇 < 40.0 GeV/𝑐 and pseudorapidity |𝜂| < 0.15. Measurements include the jet cross section, as well as distributions of SoftDrop-groomed momentum fraction (𝑧 𝑔 ), charged-particle transverse momentum with respect to jet axis (𝑗 𝑇 ), and radial distributions of charged particles within jets (𝑟). Also measured was the distribution of 𝜉 = −ln⁡(𝑧), where 𝑧 is the fraction of the jet momentum carried by the charged particle. The measurements are compared to theoretical next-to and next-to-next-to-leading-order calculations, the PYTHIA and H erwig event generators, and to other existing experimental results. Indicated from these measurements is a lower particle multiplicity in jets at RHIC energies when compared to models. Also noted are implications for future jet measurements with sPHENIX at RHIC as well as at the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Detector Simulation Challenges for Future Accelerator Experiments

Detector simulation is a key component for studies on prospective future high-energy colliders, the design, optimization, testing and operation of particle physics experiments, and the analysis of the data collected to perform physics measurements. This review starts from the current state of the art technology applied to detector simulation in high-energy physics and elaborates on the evolution of software tools developed to address the challenges posed by future accelerator programs beyond the HL-LHC era, into the 2030–2050 period. New accelerator, detector, and computing technologies set the stage for an exercise in how detector simulation will serve the needs of the high-energy physics programs of the mid 21st century, and its potential impact on other research domains.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Multiplicity dependent 𝐽/𝜓 and 𝜓⁡(2⁢𝑆) production at forward and backward rapidity in 𝑝 + 𝑝 collisions at $\sqrt{𝑠}$ = 200 GeV

Recent measurements of 𝐽/𝜓 production as a function of event charged-particle multiplicity at the collision energies of both the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) show enhanced 𝐽/𝜓 production yields with increasing multiplicity. One potential explanation for this type of dependence is multiparton interactions (MPI). We present the first study of potential autocorrelations at RHIC energies and forward and backward rapidity of self-normalized 𝐽/𝜓 yields and 𝜓⁡(2⁢𝑆) to 𝐽/𝜓 ratio, as a function of self-normalized multiplicity in 𝑝 + 𝑝 collisions. In addition, detailed pythia studies tuned to RHIC energies were performed to investigate the MPI impacts. We find that the PHENIX data at RHIC are consistent with recent LHC measurements and can only be described by pythia calculations that include MPI effects. The forward and backward 𝜓⁡(2⁢𝑆) to 𝐽/𝜓 ratio is found to be less dependent on the charged-particle multiplicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of elliptic flow of 𝐽/𝜓 in $\sqrt{s_{NN}}$ = 200 GeV ⁢Au + Au collisions at forward rapidity

Here, we report the first measurement of the azimuthal anisotropy of 𝐽/𝜓 at forward rapidity (1.2 < |𝜂| < 2.2) in Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at the BNL Relativistic Heavy Ion Collider. The data were collected by the PHENIX experiment in 2014 and 2016 with integrated luminosity of 14.5 nb −1 . The second Fourier coefficient (𝑣 2 ) of the azimuthal distribution of 𝐽/𝜓 is determined as a function of the transverse momentum (𝑝𝑇) using the event-plane method. The measurements were performed for several selections of collision centrality: 0%–50%, 10%–60%, and 10%–40%. We find that in all cases the values of 𝑣 2 ⁡(𝑝 𝑇 ), which quantify the elliptic flow of 𝐽/𝜓, are consistent with zero. Within uncertainties, the results are consistent with measurements at midrapidity, indicating no significant elliptic flow of the 𝐽/𝜓 within the quark-gluon-plasma medium at collision energies of $\sqrt{s_{NN}}$ = 200 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Global tuning of hadronic interaction models with accelerator-based and astroparticle data

In high-energy and astroparticle physics, event generators play an essential role, even in the simplest data analyses. As analysis techniques become more sophisticated, e.g. based on deep neural networks, their correct description of the observed event characteristics becomes even more important. Physical processes occurring in hadronic collisions are simulated within a Monte Carlo framework. A major challenge is the modeling of hadron dynamics at low momentum transfer, which includes the initial and final phases of every hadronic collision. QCD-inspired phenomenological models used for these phases cannot guarantee completeness or correctness over the full phase space. These models usually include parameters which must be tuned to suitable experimental data. Until now, event generators have been developed and tuned mainly on the basis of data from high-energy physics experiments at accelerators. The wealth of data available from the latest generation of astroparticle experiments has not yet been fully exploited, and in many cases is not satisfactorily described. Both kinds of data sets are complementary as astroparticle experiments provide sensitivity especially to hadrons produced nearly parallel to the collision axis and cover center-of-mass energies up to several hundred TeV, well beyond those reached at colliders so far. In this report, we provide an overview of state-of-the-art event generators and their tuning, including the most relevant inputs from high-energy accelerator and astroparticle experiments. We present a road map that shows, for the first time, how the unified tuning of event generators with accelerator-based and astroparticle data can be performed.

Albrecht, J. [Ruhr U., Bochum, RAPP Ctr.; Ruhr U.,↗
DOE OSTI2023

A new operator-split Compton scattering method

In high-energy-density physics, when the temperatures are high (keV-range) and/or the densities are low, Compton scattering becomes the main energy transfer mechanism between the radiation field and the material. The interaction between the radiation field and the material is non-linear in the radiation field – due to induced processes caused by photons being bosons – and in the material temperature. Many existing codes use a Compton-Fokker-Planck (CFP) approximation – which is inaccurate for high photon energies or material temperatures – or use linear approximations within a time-step. Here, we have implemented a new Compton-heating method called OS-FBC (Operator-Split Full-Boltzmann Compton) that solves the full, nonlinear Compton-heating system with few approximations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The Backup Program of the Dark Energy Spectroscopic Instrument’s Milky Way Survey

The Milky Way Backup Program (MWBP), a survey currently underway with the Dark Energy Spectroscopic Instrument (DESI) on the Nicholas U. Mayall 4 m Telescope, works at the margins of the DESI Main surveys to obtain spectra of millions of additional stars from the Gaia catalog. Efficiently utilizing times between ∼12° and 18° twilight and poor weather conditions, the MWBP extends the range of stellar sources studied to both brighter magnitudes and lower Galactic latitude and declination than the stars studied in DESI’s Main Milky Way Survey. While the MWBP prioritizes candidate giant stars selected from the Gaia catalog (using color and parallax criteria), it also includes an unbiased sample of bright stars (i.e., 11.2 ≲ G < 16 mag) as well as fainter sources (to G ≲ 19 mag). As of 2025 March 1, the survey had obtained spectra of ∼7 million stars, approximately 1.2 million of which are included in the DESI Data Release 1. The DESI spectra cover the wavelength range from 3600 to 9800 Å at a resolution λ/Δλ varying from 2000 to 5000. The full survey, when completed, will cover an area of more than 21,000 deg 2 and include approximately 10 million Gaia sources, roughly equal to the number of stellar spectra obtained through the DESI Main Survey, while only utilizing ≈9% of all DESI observing time. This paper provides an overview of the MWBP, describing the target selection, observing strategy, and an introduction to the resulting data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Richtmyer–Meshkov instability when a shock is reflected for fluids with arbitrary equation of state

First predicted by Richtmyer in 1960 and experimentally confirmed by Meshkov in 1969, the Richtmyer–Meshkov instability (RMI) is crucial in fields such as physics, astrophysics, inertial confinement fusion and high-energy-density physics. These disciplines often deal with strong shocks moving through condensed materials or high-pressure plasmas that exhibit non-ideal equations of state (EoS), thus requiring theoretical models with realistic fluid EoS for accurate RMI simulations. Approximate formulae for asymptotic growth rates, like those proposed by Richtmyer, are helpful but rely on heuristic prescriptions for compressible materials. These prescriptions can sometimes approximate the RMI growth rate well, but their accuracy remains uncertain without exact solutions, as the fully compressible RMI growth rate is influenced by both vorticity deposited during shock refraction and multiple sonic wave refractions. This study advances previous work by presenting an analytic, fully compressible theory of RMI for reflected shocks with arbitrary EoS. It compares theoretical predictions with heuristic prescriptions using ideal gas, van der Waals gas and three-term constitutive equations for simple metals, the latter being analysed with detailed and simplified ideal-gas-like EoS. We additionally offer an alternative explicit approximate formula for the asymptotic growth rate. The comprehensive model also incorporates the effects of constant-amplitude acoustic waves at the interface, associated with the D'yakov–Kontorovich instability in shocks.

Napieralski, Mario (ORCID:0009000692344901)↗
DOE OSTI2020

How data science methods can improve the quality and efficiency of ICF and HEDP research

Data Science methods (many that are Bayesian based) are widely used in the physical sciences to estimate model parameters from experimental data, synthesize heterogeneous data, calibrate models, design experiments, and determine statistical significance of data. These methods provide a wealth of advantages over traditional analysis techniques because: 1) uncertainties are rigorously defined and propagated naturally through complex systems including covariance, 2) prior information is captured within the analysis framework (including rad-MHD and rad-hydro simulations), 3) competing models can be selected and/or ruled out using quantitative criteria, and 4) complex, heterogeneous data can be incorporated simultaneously. While these methods have been widely adopted as the gold standard in fields such as particle physics, astronomy, and biology, they have been slow to catch on in Inertial Confinement Fusion (ICF) and High Energy Density Physics (HEDP) research. Recently, several teams at LLNL, SNL, LANL, and the LLE have been exploring the use of these tools in their research and have found success. Here we propose that a concerted effort to consolidate these independent research efforts by developing and deploying common tools for use across the complex can revolutionize the way we approach data analysis, assimilation of theory and experiment, and decision making. The Bayesian formalism provides a means to accomplish this, but we are lacking certain infrastructure to make it happen on a large scale. Furthermore, once adopted, these techniques can be used to develop standards by which discoveries can be judged, similar to the so-called 5σ rule in high energy particle physics. Such standards may be used in the future to address the issue of unknown reproducibility in ICF and HED experiments caused by low shot rate and high cost per experiment. Our goals as a group are to advance the state of the art in HED measurement science by enabling: 1) better inferences from data with well-defined uncertainties, 2) better use of the data we have and continue to collect, 3) intelligent synthesis of data, 4) evaluation of the statistical significance of our data, and 5) informed decision making regarding the design of new experiments and instruments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Experimental shock Hugoniot and melt boundary temperatures of aluminum

Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.

Hartsfield, Thomas Murray [Sandia National Laborat↗
DOE OSTI2020

The First Three Seconds: a Review of Possible Expansion Histories of the Early Universe

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The Multiphysics on Advanced Platforms Project

In 2015, the Lawrence Livermore National Laboratory started development of next-generation multiphysics simulation capabilities for the National Nuclear Security Administration under the Advanced Technologies Development and Mitigation (ATDM) element of the Advanced Simulation and Computing program in collaboration with the Exascale Computing Project (ECP). A key driver for this effort across the NNSA tri-lab was the emergence of advanced high performance computing (HPC) architectures based on heterogeneous compute capabilities, including GPU based systems, as part of the national drive toward exascale computing platforms at multiple Department of Energy (DOE) facilities. Developing a multiphysics code capable of meeting the various simulation needs of the NNSA as defined by the current generation of integrated codes (or ICs), initially developed as part of the Accelerated Strategic Computing Initiative (ASCI) program beginning in 1996, and able to scale to the current 100 petaflop class pre-exascale systems, as well the forthcoming exaflop class computers, is a daunting challenge. To accomplish this ambitious goal, LLNL has embraced two key themes: use of high-order numerical methods and a modular approach to code development. The LLNL next generation effort is organized under the Multi-Physics on Advanced Platforms Project (MAPP). A foundational component of MAPP is the Axom computer science (CS) toolkit which provides infrastructure for the development of modular, performance portable, multi-physics application codes. MARBL is a next-generation application code built on the Axom base to address the modeling needs of the high energy density physics (HEDP) community for simulating high-explosive, magnetic or laser driven experiments such as inertial confinement fusion (ICF), pulsed-power magneto-hydrodynamics (MHD), equation of state (EOS) and material strength studies as part of the NNSA’s stockpile stewardship program (SSP).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

On the study of hydrodynamic instabilities in the presence of background magnetic fields in high-energy-density plasmas

Blast-wave-driven hydrodynamic instabilities are studied in the presence of a background B-field through experiments and simulations in the high-energy-density (HED) physics regime. In experiments conducted at the Laboratoire pour l’utilisation des lasers intenses (LULI), a laser-driven shock-tube platform was used to generate a hydrodynamically unstable interface with a prescribed sinusoidal surface perturbation, and short-pulse x-ray radiography was used to characterize the instability growth with and without a 10-T B-field. The LULI experiments were modeled in FLASH using resistive and ideal magnetohydrodynamics (MHD), and comparing the experiments and simulations suggests that the Spitzer model implemented in FLASH is necessary and sufficient for modeling these planar systems. These results suggest insufficient amplification of the seed B-field, due to resistive diffusion, to alter the hydrodynamic behavior. Although the ideal-MHD simulations did not represent the experiments accurately, they suggest that similar HED systems with dynamic plasma-β (=2μ0ρv2/B2) values of less than ∼100 can reduce the growth of blast-wave-driven Rayleigh–Taylor instabilities. These findings validate the resistive-MHD FLASH modeling that is being used to design future experiments for studying B-field effects in HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Search for the 𝑌⁡(2175) in the Photoproduction Cross Section Measurement of 𝛾⁢𝑝 → 𝜙⁢𝜋 + ⁢𝜋 − ⁢𝑝 at GlueX

Based on 334 pb −1 of photoproduction data collected with the GlueX detector at Jefferson Lab, we have measured for the first time the cross section of the exclusive reaction 𝛾+𝑝→𝜙⁡(1020)⁢𝜋 + ⁢𝜋 − ⁢𝑝 by reconstructing the final state 𝐾 +⁢ 𝐾 − ⁢𝜋 + ⁢𝜋 − ⁢𝑝 produced with a photon beam of energies between 8.0 and 11.6 GeV. Based on the measured differential cross section, we have performed a search for the strangeoniumlike exotic candidate 𝑌⁡(2175), recently renamed to 𝜙⁡(2170). This state has been reported by different 𝑒 + ⁢𝑒 − annihilation experiments and it is addressed here for the first time in a photoproduction experiment. We do not find evidence for this state when using the resonance parameters quoted by the Particle Data Group and provide upper limits on the photoproduction cross section. Instead, we find a structure at a mass of 𝑚⁡(𝜙⁢𝜋 + ⁢𝜋 − ) = 2.24 GeV/𝑐 2 with a statistical significance of about 5⁢𝜎. The parameters of this structure differ from those quoted by the Particle Data Group for the 𝜙⁡(2170) and are consistent with a previous observation in 𝑒 + ⁢𝑒 − annihilation. In addition, there is evidence for a second structure at 1.82 GeV/𝑐 2 .

Exotic mesons↗
DOE OSTI2025

Performance of heavy-flavour jet identification in Lorentz-boosted topologies in proton-proton collisions at √(s) = 13 TeV

Measurements in the highly Lorentz-boosted regime provoke increased interest in probing the Higgs boson properties and in searching for particles beyond the standard model at the LHC. In the CMS Collaboration, various boosted-object tagging algorithms, designed to identify hadronic jets originating from a massive particle decaying to bb̅ or cc̅, have been developed and deployed across a range of physics analyses. This paper highlights their performance on simulated events, and summarizes novel calibration techniques using proton-proton collision data collected at √(s) = 13 TeV during the 2016–2018 LHC data-taking period. Three dedicated methods are used for the calibration in multijet events, leveraging either machine learning techniques, the presence of muons within energetic boosted jets, or the reconstruction of hadronically decaying high-energy Z bosons. The calibration results, obtained through a combination of these approaches, are presented and discussed.

Pattern recognition↗
DOE OSTI2026

Cross sections of 𝜂 mesons in 𝑝 +𝑝 collisions at forward rapidity at $\sqrt{𝑠}$ = 500 GeV and central rapidity at $\sqrt{𝑠}$ = 510 GeV

We present the first measurements of the forward and midrapidity 𝜂-meson cross sections from 𝑝 + 𝑝 collisions at $\sqrt{𝑠}$ = 500 and 510 GeV, respectively. We also report the midrapidity 𝜂/𝜋 0 ratio at 510 GeV. The forward cross section is measured differentially in 𝜂-meson transverse momentum (𝑝 𝑇 ) from 1.0 to 6.5 GeV/𝑐 for pseudorapidity 3.0 < |𝜂| < 3.8. The midrapidity cross section is measured from 3.5 to 44 GeV/𝑐 for pseudorapidity |𝜂| < 0.35. Both cross sections serve as critical inputs to an updated global analysis of the 𝜂-meson fragmentation functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗