Search NASA⌕ Search

SEARCH · Search NASA

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

Practical overview of the ICRP Publication 145 phantoms for use with MCNP

Abstract This paper provides a practical overview of the mesh-type reference phantoms (MRCPs) published by the ICRP by providing supplementary guidance on the proper implementation of using the MRCPs in computational codes, such as Monte Carlo N-Particle. Various validation scenarios for photon and neutron exposures were also examined where published values were reproduced using the MRCPs for both external and internal exposures. Generally, the MRCPs performed well in the validation scenarios with the values reproduced being within 10% of the published values.

Hetrick, Lucas (ORCID:0000000178978007)↗
DOE OSTI2025

Comparative analysis of plasmon modes in layered Lindhard metals and strange metals

The enigmatic strange metal remains one of the central unsolved problems of 21st century science. Understanding this phase of matter requires knowledge of the momentum- and energy-resolved dynamic charge susceptibility 𝜒⁡(𝑞,𝜔), especially at finite momentum. Inelastic electron scattering (EELS), performed in either transmission or reflection geometry, is a powerful probe of 𝜒⁡(𝑞,𝜔). For the prototypical strange metal Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝑥 , transmission- and reflection EELS, and infrared (IR) spectroscopy agree at 𝑞∼0, all revealing a highly damped plasmon near 1 eV. At larger 𝑞, however, EELS results show unresolved discrepancies. Since IR data are highly reproducible, it is advantageous to use IR data to calculate what the expected EELS response should be at modest 𝑞. Building on prior momentum-resolved reflection geometry M-EELS work [Chen et al., Phys. Rev. B 109, 045108 (2024)], we extend this approach to transmission EELS for finite stacks of metallic layers, comparing a “textbook” Lindhard metal to a strange metal. In the Lindhard case, the low-𝑞 response is dominated by long-lived, standing wave plasmon modes arising from interlayer Coulomb coupling, with in-plane dispersions that resemble the well-known Fetter modes of layered metals. This behavior depends only on the geometry and the long-range nature of the Coulomb interaction and is largely insensitive to layer details. At larger 𝑞, the response reflects the microscopic properties of individual layers. For the strange metal, calculations based on IR data predict a highly damped plasmon with weak dispersion and no distinct surface mode. While our results match IR and M-EELS at low 𝑞, they do not reproduce any published EELS spectra at large 𝑞, highlighting unresolved discrepancies that demand further experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Neutrino flavor instabilities in neutron star mergers with moment transport: Slow, fast, and collisional modes

Determining where, when, and how neutrino flavor oscillations must be included in large-scale simulations of hot and dense astrophysical environments is an enduring challenge that must be tackled to obtain accurate predictions. Here, using an angular moment-based linear stability analysis framework, we examine the different kinds of flavor instabilities that can take place in the context of the postprocessing of a neutron star merger simulation, with a particular focus on the collisional flavor instability and a careful assessment of several commonly used approximations. First, neglecting anisotropies of the neutrino field, we investigate the extent to which commonly used monoenergetic growth rates reproduce the results obtained from a full multienergy treatment. Contrary to the large discrepancies found in core-collapse supernova environments, we propose a simple combination of energy-averaged estimates that reproduces the multienergy growth rates in our representative simulation snapshot. We then quantify the impact of additional physical effects, including nuclear many-body corrections, scattering opacities, and the inclusion of the vacuum term in the neutrino Hamiltonian. Finally, we include the neutrino distribution anisotropies, which allows us to explore, for the first time in a multienergy setting, the interplay between collisional, fast, and slow modes in a moment-based neutron star merger simulation. We find that, despite a dominance of the fast instability in most of the simulation volume, certain regions exhibit only a collisional instability, while others, especially at large distances, exhibit a slow instability that is largely underestimated if anisotropic effects are neglected.

neutrino oscillations↗
DOE OSTI2026

Unified Description of Cuprate Superconductors by Fractionalized Electrons Emerging from Integrated Analyses of Photoemission Spectra and Quasiparticle Interference

Electronic structure of high-temperature superconducting cuprates is studied by analyzing experimental data independently obtained from two complementary spectroscopies: one, quasiparticle interference (QPI) measured by scanning-tunneling microscopy, and the other, angle-resolved photoemission spectroscopy (ARPES). We combine these two sets of data in a unified theoretical analysis. Through explicit calculations of experimentally measurable quantities, we show that a simple two-component fermion model (TCFM) representing electron fractionalization succeeds in reproducing various detailed features of these experimental data: ARPES and QPI data are concomitantly reproduced by the TCFM in full energy and momentum spaces. The measured QPI pattern reveals a signature characteristic of the TCFM, distinct from the conventional single-component prediction, supporting the validity of the electron fractionalization in the cuprates. The integrated analysis also solves the puzzles of ARPES and QPI data that are seemingly inconsistent with each other. The overall success of the TCFM offers a comprehensive understanding of the electronic structure of the cuprates, in particular, the unoccupied side of the spectra, of which momentum-resolved structure has long been unexplored experimentally. We further predict that a characteristic QPI pattern should appear in the unoccupied high-energy part if the fractionalization is at work. We propose that integrated-spectroscopy analyses offer a promising way to explore challenging issues of strongly correlated electron systems.

Sakai, Shiro [Sophia University; RIKEN Center for ↗
DOE OSTI2023

Simulating Heisenberg interactions in the Ising model with strong drive fields

The time evolution of an Ising model with large driving fields over discrete time intervals is shown to be reproduced by an effective XXZ-Heisenberg model at leading order in the inverse field strength. For specific orientations of the drive field, the dynamics of the XXX-Heisenberg model is reproduced. Finally, these approximate equivalences, valid above a critical driving field strength set by dynamical phase transitions in the Ising model, are expected to enable quantum devices that natively evolve qubits according to the Ising model to simulate more complex systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Information Carried by Electromagnetic Radiation Launched from Accelerated Polarization Currents

We show experimentally that a continuous, linear, dielectric antenna in which a superluminal polarization-current distribution accelerates can be used to transmit a broadband signal that is reproduced in a comprehensible form at a chosen target distance and angle. The requirement for this exact correspondence between broadcast and received signals is that each moving point in the polarization-current distribution approaches the target at the speed of light at all times during its transit along the antenna. This results in a one-to-one correspondence between the time at which each point on the moving polarization current enters the antenna and the time at which all of the radiation emitted by this particular point during its transit through the antenna arrives simultaneously at the target. This has the effect of reproducing the desired time dependence of the original broadcast signal. For other observer-detector positions, the time dependence of the signal is scrambled, due to the nontrivial relationship between emission (retarded) time and reception time. This technique represents a contrast to conventional radio transmission methods; in most examples of the latter, signals are broadcast with little or no directivity, selectivity of reception being achieved through the use of narrow frequency bands. In place of this, the current paper uses a spread of frequencies to transmit information to a particular location; the signal is weaker and has a scrambled time dependence elsewhere. Here, we point out the possible relevance of this mechanism to 5G neighborhood networks and pulsar astronomy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Structure and energetics of the elbows in the Au(111) herringbone reconstruction

We study the structure of the threading edge dislocations, or “elbows,” which are an essential component of the well-known herringbone reconstruction of the (111) surface of Au. Previous work had shown that these dislocations can be stabilized by long-range elastic relaxations into the bulk. However, the validity of the harmonic spring model that had been used to estimate the energies of the dislocations is uncertain. To enable a more refined model of the dislocation energetics, we have imaged the atomic structure of these dislocations using scanning tunneling microscopy. We find that the harmonic spring model does not adequately reproduce the observed structure. We are able to reproduce the structure, however, with a two-dimensional Frenkel-Kontorova (FK) model that uses a pairwise Morse potential to describe the interactions between the top layer Au atoms on a rigid substrate. The parameters of the potential were obtained by fitting the energy of uniaxially compressed phases, or “stripes”, computed with density functional theory, as a function of surface Au density. Within this model, the formation of the threading dislocations remains unfavorable. However, the large forces on the substrate atoms near the threading-dislocation cores, render the assumption of a completely rigid substrate questionable. Indeed, if the FK parameters are modified to account for the relaxation of just one more atomic layer, threading dislocations can, in principle, become favorable, even without bulk elastic relaxations. Additional evidence for a small elbow energy is that our computed change in the Au(111) surface stress tensor caused by the (√ 3 × 22) reconstruction is considerably smaller than previous estimates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Boiling and cavitation caused by transient heat transfer in superfluid helium-4

Superfluid helium-4 (He II) has been widely utilized as a coolant in various scientific and engineering applications due to its superior heat transfer capability. An important parameter required in the design of many He II based cooling systems is the peak heat flux $q^*$, which refers to the threshold heat flux above which boiling spontaneously occurs in He II. Past experimental and numerical studies showed that $q^*$ increases when the heating time $t_h$ is reduced, which leads to an intuitive expectation that very high $q^*$ may be achievable at sufficiently small $t_h$. Knowledge on how $q^*$ actually behaves at small $t_h$ is important for applications such as laser ablation in He II. Furthermore we present a numerical study on the evolution of the thermodynamic state of the He II in front of a planar heater by solving the He II two-fluid equations of motion. For an applied heat flux, we determine the heating time beyond which the He II near the heater transits to the vapor phase. As such, a curve correlating $q^*$ and $t_h$ can be obtained, which nicely reproduces some relevant experimental data. Surprisingly, we find that there exists a critical peak heat flux $q^*_c$, above which boiling occurs nearly instantaneously regardless of $t_h$. We reveal that the boiling in this regime is essentially cavitation caused by the combined effects of the first-sound and the second-sound waves in He II. Based on this physical picture, an analytical model for $q^*_c$ is developed, which reproduces the simulated $q^*_c$ values at various He II bath temperatures and hydrostatic head pressures. This work represents a major progress in our understanding of transient heat transfer in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Hydrodynamic theory of scrambling in chaotic long-range interacting systems

The Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) equation provides a mean-field theory of out-of-time-ordered commutators in locally interacting quantum chaotic systems at high energy density. In systems with power-law interactions, the corresponding fractional-derivative FKPP equation provides an analogous mean-field theory. However, the fractional FKPP description is potentially subject to strong quantum fluctuation effects, so it is not clear a priori if it provides a suitable effective description for generic chaotic systems with power-law interactions. Here, in this work, we study this problem using a model of coupled quantum dots with interactions decaying as 1/r α , where each dot hosts N degrees of freedom. The large-N limit corresponds to the mean-field description, while quantum fluctuations contributing to the OTOC can be modeled by 1/N corrections consisting of a cutoff function and noise. Within this framework, we show that the parameters of the effective theory can be chosen to reproduce the butterfly light cone scalings previously found for N=1 and generic finite N. In order to reproduce these scalings, the fractional index μ in the FKPP equation needs to be shifted from the naïve value of μ=2⁢α–1 to a renormalized value μ=2⁢α–2. We provide supporting analytic evidence for the cutoff model and numerical confirmation for the full fractional FKPP equation with cutoff and noise.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Magic-angle twisted symmetric trilayer graphene as a topological heavy-fermion problem

Recently, Song and Bernevig [Phys. Rev. Lett. 129, 047601 (2022)] reformulated magic-angle twisted bilayer graphene as a topological heavy fermion problem, and used this reformulation to provide a deeper understanding for the correlated phases at integer fillings. Here, in this work, we generalize this heavy-fermion paradigm to magic-angle twisted symmetric trilayer graphene, and propose a low-energy f–c–d model that reformulates magic-angle twisted symmetric trilayer graphene as heavy localized f modes coupled to itinerant topological semimetalic c modes and itinerant Dirac d modes. Our f–c–d model well reproduces the single-particle band structure of magic-angle twisted symmetric trilayer graphene at low energies for displacement field $\mathcal{E}$ ϵ [0,300]⁢ meV. By performing Hartree-Fock calculations with the f–c–d model for v = 0,–1,–2 electrons per Moiré unit cell, we reproduce all the correlated ground states obtained from the previous numerical Hartree-Fock calculations with the Bistritzer-MacDonald-type model, and we find additional new correlated ground states at high displacement field. Based on the numerical results, we propose a simple rule for the ground states at high displacement fields by using the f–c–d model, and provide analytical derivation for the rule at charge neutrality. We also provide analytical symmetry arguments for the (nearly) degenerate energies of the high-$\mathcal{E}$ ground states at all the integer fillings of interest, and make experimental predictions of which charge-neutral states are stabilized in magnetic fields. Our f–c–d model provides a new perspective for understanding the correlated phenomena in magic-angle twisted symmetric trilayer graphene, suggesting that the heavy fermion paradigm of Song and Bernevig [Phys. Rev. Lett. 129, 047601 (2022)] should be the generic underpinning of correlated physics in multilayer moire graphene structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Nonlinear sigma model approach to many-body quantum chaos: Regularized and unregularized out-of-time-ordered correlators

Here, the out-of-time-ordered correlators (OTOCs) have been proposed and widely used recently as a tool to define and describe many-body quantum chaos. Here, we develop the Keldysh nonlinear sigma model technique to calculate these correlators in interacting disordered metals. In particular, we focus on the regularized and unregularized OTOCs. The calculation of the rate of OTOCs' exponential growth is reminiscent to that of the Altshuler-Aronov-Khmelnitskii dephasing rate in interacting metals, but here it involves two replicas of the system (two “worlds”). The intraworld contributions reproduce the Altshuler-Aronov-Khmelnitskii dephasing (that would correspond to a decay of the correlator), while the interworld terms provide a term of the opposite sign that exceeds dephasing. Consequently, both regularized and unregularized OTOCs grow exponentially in time, but surprisingly we find that the corresponding many-body Lyapunov exponents are different. For the regularized correlator, we reproduce an earlier perturbation theory result for the Lyapunov exponent that satisfies the Maldacena-Shenker-Stanford bound. However, the Lyapunov exponent of the unregularized correlator parametrically exceeds the bound, is not a reliable indicator of many-body quantum chaos as it contains additional contributions from elastic scattering events due to virtual processes that should not contribute to many-body chaos. These results bring up an important general question of the physical meaning of the OTOCs often used in calculations and proofs. We briefly discuss possible connections of the OTOCs to observables in quantum interference effects and level statistics via a generalization of the Bohigas-Giannoni-Schmit conjecture to many-body chaotic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Geometry and electronic structure of iridium adsorbed on graphene

We report investigation of the geometry and electronic structure of iridium atoms adsorbed onto graphene through a combined experimental and theoretical study. Ir atoms were deposited onto a flake of graphene on a Pt(111) surface and found to form clusters even at low temperatures. The areal density of the observed clusters on the graphene flake suggests the clusters are most likely pairs of Ir atoms. Theoretical ab initio density functional (DFT) calculations indicate that these Ir dimers are oriented horizontally, near neighboring "bridge" sites of the graphene lattice, as this configuration has the strongest adsorption energy of all high-symmetry configurations for the Ir dimer. A large peak in the local density of states (LDOS) at the Dirac point energy was measured via scanning tunneling spectroscopy, and this result is reproduced by a DFT calculation of the LDOS. The peak at the Dirac point energy is found to be from the Ir s and p states. Finally, the LDOS in the monomer case was also calculated, and is found to significantly differ from the experimentally determined data, further supporting the hypothesis of low-temperature clustering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Numerical investigation of spin excitations in a doped spin chain

We study the doping evolution of spin excitations in a one-dimensional (1D) Hubbard model and its downfolded spin Hamiltonians, by using exact diagonalization combined with cluster perturbation theory. In all models we observe hardening (softening) of spin excitations upon electron (hole) doping, which are reminiscent of recent experiments on two-dimensional (2D) cuprate materials. We also find that the three-site and even higher-order terms are crucial for the low-energy effective spin models to reproduce the magnetic spectra of doped Hubbard systems at a quantitative level. To interpret the numerical results, we further employ a strong coupling slave-boson mean-field theory. The mean-field theory provides an intuitive understanding of the overall compact support of dynamic spin structure factors, including the shift of zero-energy modes and change of spin excitation bandwidth with doping. Our results can serve as predictive benchmarks for future inelastic x-ray or neutron scattering experiments on doped 1D antiferromagnetic Mott insulators.

Pärschke, Ekaterina M.↗
DOE OSTI2019

Exact isovector pairing in a shell-model framework: Role of proton-neutron correlations in isobaric analog states

Here, we utilize a nuclear shell model Hamiltonian with only two adjustable parameters to generate, for the first time, exact solutions for pairing correlations for light to medium-mass nuclei, including the challenging proton-neutron pairs, while also identifying the primary physics involved. In addition to single-particle energy and Coulomb potential terms, the shell model Hamiltonian consists of an isovector T = 1 pairing interaction and an average proton-neutron isoscalar T = 0 interaction, where the T = 0 term describes the average interaction between nonpaired protons and neutrons. This Hamiltonian is exactly solvable, where, utilizing three to seven single-particle energy levels, we reproduce experimental data for 0 + state energies for isotopes with mass A = 10 through A = 62 exceptionally well including isotopes from He to Ge. Additionally, we isolate effects due to like-particle and proton-neutron pairing, provide estimates for the total and proton-neutron pairing gaps, and reproduce N (neutron) = Z (proton) irregularity. These results provide a further understanding for the key role of proton-neutron pairing correlations in nuclei, which is especially important for waiting-point nuclei on the $\mathcal{rp}$ path of nucleosynthesis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

HIJING can describe the anisotropy-scaled charge-dependent correlations at the BNL Relativistic Heavy Ion Collider

The experimentally measured charge-dependent correlations in heavy-ion collisions have been suggested as a signature of the chiral magnetic effect (CME). Early model studies could not reproduce the measurement. For example, the HIJING model yielded far smaller magnitude for the charge-dependent correlation than observed in data. This led to the conclusion that the CME had to be invoked to explain the observed correlations in heavy-ion collisions. In this paper we show that this conclusion of the CME interpretation is premature. We show that the reason that HIJING predicts a far smaller correlation than data is because the elliptic anisotropy $(v_ 2)$ parameter in HIJING is too small. When properly scaled, the HIJING model can reproduce in entirety the measured correlations. Further, we also employ the AMPT model, which has a large enough $v_2$, to demonstrate that the measured data can be easily accommodated by models without invoking the CME.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Excluded-volume model for quarkyonic matter: Three-flavor baryon-quark mixture

The single-flavor excluded-volume model based on the effective size of baryons reproduces the hard-soft density evolution of the equation of state (EoS) required by the recent studies of GW170817. This phenomenological model basically realizes the concept of quarkyonic matter which is introduced from large-$N_c$ gauge theory for dense matter. Further, enhanced nucleon interactions and dynamically generated quark degrees of freedom can reproduce the hard-soft evolution of the EoS. In this paper, we extend the excluded-volume model to a three-flavor system by considering electromagnetic charge and possible weak equilibrium in order to obtain a proper description for the hard-soft behavior of the EoS inferred from the gravitational waves observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Total absorption spectroscopy of the β decay of 101,102 Zr and 109 Tc

The β decay of 101,102 Zr and 109 Tc was studied using the technique of total absorption spectroscopy. The experiment was performed at the National Superconducting Cyclotron Laboratory using the Summing NaI(Tl) (SuN) detector in the first-ever application of total absorption spectroscopy with a fast beam produced via projectile fragmentation. The β -decay feeding intensity and Gamow-Teller transition strength distributions were extracted for these three decays. The extracted distributions were compared to three different quasiparticle random-phase approximation (QRPA) models based on different mean-field potentials. A comparison with calculations from one of the QRPA models was performed to learn about the ground-state shape of the parent nucleus. For 101 Zr and 102 Zr, calculations assuming a pure shape configuration (oblate or prolate) were not able to reproduce the extracted distributions. These results may indicate that some type of mixture between oblate and prolate shapes is necessary to reproduce the extracted distributions. For 109 Tc, a comparison of the extracted distributions with QRPA calculations suggests a dominant oblate configuration. The other two QRPA models are commonly used to provide β-decay properties in r-process network calculations. This work shows the importance of making comparisons between the experimental and theoretical β-decay distributions, rather than just half-lives and β-delayed neutron emission probabilities, as close to the r-process path as possible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Investigating high-energy proton-induced reactions on spherical nuclei: Implications for the preequilibrium exciton model

A number of accelerator-based isotope production facilities utilize $100-200$ MeV proton beams due to the high production rates enabled by high-intensity beam capabilities and the greater diversity of isotope production brought on by the long range of high-energy protons. However, nuclear reaction modeling at these energies can be challenging because of the interplay between different reaction modes and a lack of existing guiding cross section data. A Tri-lab collaboration has been formed between the Lawrence Berkeley, Los Alamos, and Brookhaven National Laboratories to address these complexities by characterizing charged-particle nuclear reactions relevant to the production of established and novel radioisotopes. In the inaugural collaboration experiments, stacked-targets of niobium foils were irradiated at the Brookhaven Linac Isotope Producer ($E_p=200$ MeV) and the Los Alamos Isotope Production Facility ($E_p=100$ MeV) to measure $^{93}$Nb(p,x) cross sections between $50-200$ MeV. The results were compared with literature data as well as the default calculations of the nuclear model codes TALYS, CoH, EMPIRE, and ALICE. The default code predictions largely failed to reproduce the measurements. Therefore, we developed a standardized procedure, which determines the reaction model parameters that best reproduce the most prominent reaction channels in a physically justifiable manner. Overall, the primary focus of the procedure was to determine the best parameterization for the pre-equilibrium two-component exciton model. This modeling study revealed a trend towards a relative decrease for internal transition rates at intermediate proton energies ($E_p=20-60$ MeV) in the current exciton model as compared to the default values. The results of this work are instrumental for the planning, execution, and analysis essential to isotope production.

43 PARTICLE ACCELERATORS↗