Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility”

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 OSTI2019

Assessing saturation physics explanations of collectivity in small collision systems with the ip-jazma model

Experimental measurements in relativistic collisions of small systems from $p + p$ to $p/d/ ^3He + A$ at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) reveal particle emission patterns that are strikingly similar to those observed in $A + A$ collisions of large nuclei. One explanation of these patterns is the formation of small droplets of quark-gluon plasma (QGP) followed by hydrodynamic evolution. A geometry engineering program was proposed to further investigate these emission patterns, and the experimental data from that program in $p + Au, d + Au, ^3 He + Au$ collisions for elliptic and triangular anisotropy coefficients $v_2$ and $v_3$ follow the pattern predicted by hydrodynamic calculations [C. Aidala et al. (PHENIX Collaboration), Nat. Phys. 15, 214 (2019)]. One alternative approach, referred to as initial-state correlations, suggests that for small systems the patterns observed in the final-state hadrons are encoded at the earliest moments of the collision and therefore require no final-state parton scattering or hydrodynamic evolution. Recently, new calculations using only initial-state correlations, in the dilute-dense approximation of gluon saturation physics, reported striking agreement with the $v_2$ patterns observed in $p/d/ ^3He + Au$ data at RHIC [M. Mace, V. V. Skokov, P. Tribedy, and R. Venugopalan, Phys. Rev. Lett. 121, 052301 (2018)]. The results reported by Mace, Skokov, Tribedy and Venugopalan (MSTV) are counterintuitive and thus we aim here to reproduce some of the basic features of these calculations. In this first investigation, we provide a description of our publicly available model, ip-jazma, and investigate its implications for saturation scales, multiplicity distributions, and eccentricities, reserving for later work the analysis of momentum spectra and azimuthal anisotropies. We find that our implementation of the saturation physics model reproduces the results of the MSTV calculation of the multiplicity distribution in $d + Au$ collisions at RHIC. However, additional aspects of studies, together with existing data, call into question some of the essential elements reported by MSTV. Resolution of these issues will require further developments of ip-jazma, in order to determine if it can replicate the qualitative agreement with the $v_2$ reported by MSTV. Both the work reported here and future studies will establish which features in the experimental data are uniquely attributable to the color glass condensate description.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Investigating the XENON1T low-energy electronic recoil excess using NEST

The search for dark matter, the missing mass of the Universe, is one of the most active fields of study within particle physics. The XENON1T experiment recently observed a 3.5⁢𝜎 excess potentially consistent with dark matter, or with solar axions. Here, we will use the Noble Element Simulation Technique (NEST) software to simulate the XENON1T detector, reproducing the excess. We utilize different detector efficiency and energy reconstruction models, but they primarily impact sub-keV energies and cannot explain the XENON1T excess. However, using NEST, we can reproduce their excess in multiple, unique ways, most easily via the addition of 31 ± 11 37 Ar decays. Furthermore, this results in new, modified background models, reducing the significance of the excess to ≤2.2⁢𝜎 at least using non-Profile Likelihood Ratio (PLR) methods. This is independent confirmation that the excess is a real effect, but potentially explicable by known physics. Many cross-checks of our 37 Ar hypothesis are presented.

particle dark matter↗
DOE OSTI2022

Dark Energy Survey Year 3 results: Cosmology from combined galaxy clustering and lensing validation on cosmological simulations

In this work, we present a validation of the Dark Energy Survey Year 3 (DES Y3) 3×2-point analysis choices by testing them on Buzzard2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak-lensing analyses. We show that the buzzard2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak-lensing statistics. We then show that our model for the 3×2-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the nonlinear matter power spectrum, nonlinear galaxy bias, and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy–galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Unified description of DGLAP, CSS, and BFKL evolution: TMD factorization bridging large and small x

This paper introduces a transverse-momentum dependent (TMD) factorization scheme designed to unify both large and small Bjorken-x regimes. We compute the next-to-leading order quantum chromodynamics (QCD) corrections to the gluon TMD operator for an unpolarized hadron within this proposed scheme. This leads to the emergence of a new TMD evolution, incorporating those in transverse momentum, rapidity, and Bjorken-x. When matched to the collinear factorization scheme, our factorization scheme faithfully reproduces the well-established Dokshitzer-Gribov-Lipatov-Altarelli-Parisi and Collins-Soper-Sterman (CSS) evolutions. Conversely, matching with high-energy factorization not only yields the Balitsky-Fadin-Kuraev-Lipatov evolution but also reveals distinctive signatures of CSS logarithms. The development of this novel TMD factorization scheme, capable of seamlessly reconciling disparate Bjorken-x regimes and faithfully reproducing established QCD evolution equations, has the potential to significantly advance our comprehension of high-energy processes and three-dimensional parton structures of hadrons. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Simulating a numerical UV completion of quartic Galileons

The Galileon theory is a prototypical effective field theory that incorporates the Vainshtein screening mechanism—a feature that arises in some extensions of general relativity, such as massive gravity. The Vainshtein effect requires that the theory contain higher order derivative interactions, which results in Galileons, and theories like them, failing to be technically well posed. While this is not a fundamental issue when the theory is correctly treated as an effective field theory, it nevertheless poses significant practical problems when numerically simulating this model. These problems can be tamed using a number of different approaches: introducing an active low-pass filter and/or constructing a UV completion at the level of the equations of motion, which controls the high momentum modes. These methods have been tested on cubic Galileon interactions, and have been shown to reproduce the correct low-energy behavior. Here we show how the numerical UV-completion method can be applied to quartic Galileon interactions, and present the first simulations of the quartic Galileon model using this technique. We demonstrate that our approach can probe physics in the regime of the effective field theory in which the quartic term dominates, while successfully reproducing the known results for cubic interactions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Sub-MeV dark sink dark matter

A dark sink uses dark-sector interactions to siphon energy from dark matter to lighter dark degrees of freedom, i.e., dark radiation. Here, we extend dark matter models containing a dark sink to sub-MeV masses. We consider a dark sink model where the dark matter is charged under a light dark photon that has kinetic mixing with the Standard Model. For sub-MeV dark matter masses, plasmon decays are the dominant mechanism for transferring energy to the dark sector. Relative to a standard freeze-in cosmology, reproducing the observed dark matter density in a dark sink structure requires an increase in the dark matter couplings to the Standard Model, and hence increased direct detection cross sections. These models provide benchmarks for current and upcoming direct detection experiments. Accounting for plasmon effects, we derive the range of possible dark matter masses and cross sections for dark sink models in the sub-MeV regime. We make the freezein code available to reproduce our benchmarks; it may be of use for other freeze-in scenarios, including those where plasmon decays to the dark matter are important.

Dark matter detectors↗
DOE OSTI2025

Phases of 2D massless QCD with qubit regularization

We investigate the possibility of reproducing the continuum physics of 2D S U ( N ) gauge theory coupled to a single flavor of massless Dirac fermion using qubit regularization. The continuum theory is described by N free fermions in the ultraviolet (UV) and a coset Wess-Zumino-Witten (WZW) model in the infrared (IR). In this work, we first explore how well these features can be reproduced using the Kogut-Susskind (KS) Hamiltonian with a finite-dimensional link Hilbert space and a generalized Hubbard coupling. We do this by analyzing the renormalization group (RG) flow diagram of the continuum theory and identifying important phases of the theory. Using strong coupling expansions, we show that our lattice model exhibits a gapped dimer phase and a spin-chain phase. Furthermore, for N = 2 , using tensor network methods, we show that there is a second-order phase transition between these two phases, which we identify as the critical surface of the continuum theory that connects the IR and UV fixed points. In the IR, we identify the critical theory at the transition as the expected S U ( 2 ) 1 WZW model. Lastly, we argue that modifications of our model may allow the study of the UV physics of free fermions. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Polarization of Λ ( Λ ¯ ) Hyperons along the Beam Direction in Au + Au Collisions at s N N = 200 GeV

The Λ ($\bar{Λ}$) hyperon polarization along the beam direction has been measured in Au+Au collisions at $\sqrt{s_{NN}}$=200 GeV, for the first time in heavy-ion collisions. The polarization dependence on the hyperons' emission angle relative to the elliptic flow plane exhibits a second harmonic sine modulation, indicating a quadrupole pattern of the vorticity component along the beam direction, expected due to elliptic flow. The polarization is found to increase in more peripheral collisions, and shows no strong transverse momentum ($p_{T}$) dependence at $p_{T}$ greater than 1 GeV/c. Furthermore, the magnitude of the signal is about 5 times smaller than those predicted by hydrodynamic and multiphase transport models; the observed phase of the emission angle dependence is also opposite to these model predictions. In contrast, the kinematic vorticity calculations in the blast-wave model tuned to reproduce particle spectra, elliptic flow, and the azimuthal dependence of the Gaussian source radii measured with the Hanbury Brown-Twiss intensity interferometry technique reproduce well the modulation phase measured in the data and capture the centrality and transverse momentum dependence of the polarization signal.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Nonlinear Magnetization Dynamics Driven by Strong Terahertz Fields

We present a comprehensive experimental and numerical study of magnetization dynamics in a thin metallic film triggered by single-cycle terahertz pulses of ~20 MV/m electric field amplitude and ~1 ps duration. The experimental dynamics is probed using the femtosecond magneto-optical Kerr effect, and it is reproduced numerically using macrospin simulations. The magnetization dynamics can be decomposed in three distinct processes: a coherent precession of the magnetization around the terahertz magnetic field, an ultrafast demagnetization that suddenly changes the anisotropy of the film, and a uniform precession around the equilibrium effective field that is relaxed on the nanosecond time scale, consistent with a Gilbert damping process. Macrospin simulations quantitatively reproduce the observed dynamics, and allow us to predict that novel nonlinear magnetization dynamics regimes can be attained with existing tabletop terahertz sources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Self-Interacting Dark Matter and the Origin of Ultradiffuse Galaxies NGC1052-DF2 and -DF4

Observations of ultradiffuse galaxies NGC 1052-DF2 and -DF4 show they may contain little dark matter, challenging our understanding of galaxy formation. Using controlled N-body simulations, we explore the possibility that their properties can be reproduced through tidal stripping from the elliptical galaxy NGC 1052, in both cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. To explain the dark matter deficiency, we find that a CDM halo must have a very low concentration so that it can lose sufficient inner mass in the tidal field. In contrast, SIDM favors a higher and more reasonable concentration as core formation enhances tidal mass loss. Final stellar distributions in our SIDM benchmarks are more diffuse than the CDM one, and hence the former are in better agreement with the data. We further show that a cored CDM halo model modified by strong baryonic feedback is unlikely to reproduce the observations. Our results indicate that SIDM is more favorable for the formation of dark-matter-deficient galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

First Application of Mass Measurements with the Rare-RI Ring Reveals the Solar r -Process Abundance Trend at A = 122 and A = 123

The Rare-RI Ring (R3) is a recently commissioned cyclotronlike storage ring mass spectrometer dedicated to mass measurements of exotic nuclei far from stability at Radioactive Isotope Beam Factory (RIBF) in RIKEN. The first application of mass measurement using the R3 mass spectrometer at RIBF is reported. Rare isotopes produced at RIBF— 127 Sn, 126 In, 125 Cd, 124 Ag, 123 Pd—were injected in R3. Masses of 126 In, 125 Cd, and 123 Pd were measured whereby the mass uncertainty of 123 Pd was improved. This is the first reported measurement with a new storage ring mass spectrometry technique realized at a heavy-ion cyclotron and employing individual injection of the preidentified rare nuclei. The latter is essential for the future mass measurements of the rarest isotopes produced at RIBF. The impact of the new 123 Pd result on the solar r-process abundances in a neutron star merger event is investigated by performing reaction network calculations of 20 trajectories with varying electron fraction Y e . It is found that the neutron capture cross section on 123 Pd increases by a factor of 2.2 and β-delayed neutron emission probability, P 1n , of 123 Rh increases by 14%. The neutron capture cross section on 122 Pd decreases by a factor of 2.6 leading to pileup of material at A = 122, thus reproducing the trend of the solar r-process abundances. The trend of the two-neutron separation energies (S 2n ) was investigated for the Pd isotopic chain. The new mass measurement with improved uncertainty excludes large changes of the S 2n value at N = 77. Finally, such large increase of the S 2n values before N = 82 was proposed as an alternative to the quenching of the N = 82 shell gap to reproduce r-process abundances in the mass region of A = 112 – 124.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

β + Gamow-Teller Strengths from Unstable 14 O via the (d, 2 He) Reaction in Inverse Kinematics

For the first time, the (d, 2 He) reaction was successfully used in inverse kinematics to extract the Gamow-Teller transition strength in the β + direction from an unstable nucleus. The new technique was made possible by the use of an active-target time-projection chamber and a magnetic spectrometer, and opens a path to addressing a range of scientific challenges, including in astrophysics and neutrino physics. Here, in this Letter, the nucleus studied was 14 O, and the Gamow-Teller transition strength to 14 N was extracted up to an excitation energy of 22 MeV. The data were compared to shell-model and state-of-the-art coupled-cluster calculations. Shell-model calculations reproduce the measured Gamow-Teller strength distribution up to about 15 MeV reasonably well, after the application of a phenomenological quenching factor. In a significant step forward to better understand this quenching, the coupled-cluster calculation reproduces the full strength distribution well without such quenching, owing to the large model space, the inclusion of strong correlations, and the coupling of the weak interaction to two nucleons through two-body currents.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Enhanced Deuteron Coalescence Probability in Jets

The transverse-momentum ($p_T$) spectra and coalescence parameters $B_2$ of (anti)deuterons are measured in p-p collisions at $\sqrt{s}=13$ $\text{TeV}$ for the first time in and out of jets. In this measurement, the direction of the leading particle with the highest $p_T$ in the event ($p^{lead}_T > 5 GeV/c$) is used as an approximation for the jet axis. The event is consequently divided into three azimuthal regions, and the jet signal is obtained as the difference between the toward region, that contains jet fragmentation products in addition to the underlying event (UE), and the transverse region, which is dominated by the UE. The coalescence parameter in the jet is found to be approximately a factor of 10 larger than that in the underlying event. This experimental observation is consistent with the coalescence picture and can be attributed to the smaller average phase-space distance between nucleons in the jet cone as compared with the underlying event. The results presented in this Letter are compared to predictions from a simple nucleon coalescence model, where the phase-space distributions of nucleons are generated using pythia8 with the Monash 2013 tuning, and to predictions from a deuteron production model based on ordinary nuclear reactions with parametrized energy-dependent cross sections tuned on data. The latter model is implemented in pythia8.3. Both models reproduce the observed large difference between in-jet and out-of-jet coalescence parameters, although the almost flat trend of the $B^{Jet}_2$ is not reproduced by the models, which instead give a decreasing trend.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Using Direct Laboratory Measurements of Electron Temperature Anisotropy to Identify the Heating Mechanism in Electron-Only Guide Field Magnetic Reconnection

Anisotropic electron heating during electron-only magnetic reconnection with a large guide magnetic field is directly measured in a laboratory plasma through in situ measurements of electron velocity distribution functions. Electron heating preferentially parallel to the magnetic field is localized to one separatrix, and anisotropies of 1.5 are measured. The mechanism for electron energization is identified as the parallel reconnection electric field because of the anisotropic nature of the heating and spatial localization. These characteristics are reproduced in a 2D particle-in-cell simulation and are also consistent with numerous magnetosheath observations. A measured increase in the perpendicular temperature along both separatrices is not reproduced by our 2D simulations. Here, this work has implications for energy partition studies in magnetosheath and laboratory reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Comparing time-dependent density functional theory with many-body perturbation theory for semiconductors: Screened range-separated hybrids and the GW plus Bethe-Salpeter approach

We present band structure and optical absorption spectra obtained from density functional theory (DFT) and linear response time-dependent DFT (TDDFT) calculations using a screened range-separated hybrid (SRSH) functional, including spin-orbit coupling, for seven prototypical semiconductors. The results are compared to those obtained from highly converged many-body perturbation theory calculations using the GW approximation and the GW plus Bethe-Salpeter equation (GW-BSE) approaches. We use a single empirical parameter for our SRSH calculations, fit such that the SRSH band gap reproduces the GW band gap at the Γ point. We then find that ground-state generalized Kohn-Sham SRSH eigenvalues accurately reproduce the band structure obtained from GW calculations, typically to within 0.1–0.2 eV, and optical absorption spectra obtained using TDDFT with the SRSH functional agree well with those of GW-BSE, with a mean deviation of 0.03 and 0.11 eV for the location of the first and second absorption peaks, respectively, at a fraction of the computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Finite-temperature simulation of anharmonicity and octahedral tilting transitions in halide perovskites

Octahedral tilting transitions are observed in most inorganic halide perovskites and play an important role in determining their functional and thermodynamic properties. Despite existing near room temperature, the cubic and tetragonal forms of halide perovskites become dynamically unstable at low temperature, making it impossible to study their thermodynamic properties with commonly used quasiharmonic models. An anharmonic vibrational Hamiltonian is constructed that accurately reproduces the low-energy portion of the potential-energy surface of the halide perovskite CsPbBr 3 . The Hamiltonian is validated using a large first-principles dataset of energies calculated within density functional theory for large-amplitude deformations of the CsPbBr 3 crystal. Monte Carlo simulations performed on the Hamiltonian reproduce the orthorhombic-tetragonal-cubic phase transitions observed in CsPbBr 3 and many other halide perovskites, demonstrating the importance of anharmonic vibrational excitations in stabilizing the tetragonal and cubic phases in these materials. Measures of local structure and octahedral tilting in the cubic and tetragonal phases, obtained from Monte Carlo simulations, confirm the connection between large anisotropic displacement factors and octahedral tilting, as observed experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronegative metal dopants improve switching variability in Al 2 ⁢O 3 resistive switching devices

Resistive random-access memories are promising for nonvolatile memory and brain-inspired computing applications. High variability and low yield of these devices are key drawbacks hindering reliable training of physical neural networks. In this paper, we show that doping an oxide electrolyte, Al 2 ⁢O 3 , with electronegative metals makes resistive switching significantly more reproducible, surpassing the reproducibility requirements for obtaining reliable hardware neuromorphic circuits. Based on density functional theory calculations, the underlying mechanism is hypothesized to be the ease of creating oxygen vacancies in the vicinity of electronegative dopants due to the capture of the associated electrons by dopant midgap states and the weakening of Al-O bonds. These oxygen vacancies and vacancy clusters also bind significantly to the dopant, thereby serving as preferential sites and building blocks in the formation of conducting paths. Throughout this work, we validate this theory experimentally by implanting different dopants over a range of electronegativities in devices made of multiple alternating layers of Al 2 ⁢O 3 and WN and find superior repeatability and yield with highly electronegative metals, Au, Pt, and Pd. These devices also exhibit a gradual SET transition, enabling multibit switching that is desirable for analog computing.

36 MATERIALS SCIENCE↗