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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2021

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Differential two-particle number and momentum correlations with the AMPT, UrQMD, and EPOS models in Pb-Pb collisions at s N N = 2.76 TeV

We report studies of charge-independent and charge-dependent two-particle differential number correlation functions, R 2 (Δη,Δφ), and transverse momentum correlation functions, P 2 (Δη,Δφ), of charged particles produced in Pb–Pb collisions at the LHC center-of-mass energy √sNN= 2.76 TeV with the UrQMD, AMPT, and EPOS models. Model calculations for R2 and P2 correlation functions are presented for inclusive charged hadrons (h ± ) in selected transverse momentum ranges and with full azimuthal coverage in the pseudorapidity range |η|<1.0. We compare these calculations for the strength, shape, and particularly the width of the correlation functions with recent measurements of these observables by the ALICE collaboration. Our analysis indicates that comparative studies of R 2 and P 2 correlation functions provide valuable insight towards the understanding of particle production in Pb–Pb collisions. We find, in particular, that these models produce quantitatively different magnitudes and shapes for these correlation functions and none reproduce the results reported by the ALICE collaboration. Accounting for quantum number conservation in models, particularly charge conservation, is mandatory to reproduce the detailed measurements of number and transverse momentum correlation functions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Applications of reduced-basis methods to the nuclear single-particle spectrum

Reduced-basis methods provide a powerful framework for building efficient and accurate emulators. Although widely applied in many fields to simplify complex models, reduced-basis methods have only been recently introduced into nuclear physics. In this Letter we build an emulator to study the single-particle structure of atomic nuclei. By scaling a suitable mean-field Hamiltonian, a “universal” reduced basis is constructed capable of accurately and efficiently reproduce the entire single-particle spectrum of a variety of nuclei. Indeed, the reduced-basis model reproduces both ground- and excited-state energies as well as the associated wave functions with remarkable accuracy. Here our results bode well for more demanding applications that use Bayesian optimization to calibrate nuclear energy density functionals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Secondary neutron production from thick target deuteron breakup

Thick target deuteron breakup is a variable-energy accelerator-based source of high-energy neutrons, with applications in fundamental and applied nuclear science and engineering. However, the breakup mechanism remains poorly understood, and data on neutron yields from thick target breakup remain relatively scarce. In this work, the double-differential neutron yields from deuteron breakup have been measured on a thick beryllium target at ϵ d = 33 and 40 MeV, using both time-of-flight and activation techniques. We have also introduced a simple hybrid model for the double-differential deuteron breakup cross section, applicable in the ϵ d = 10-100 MeV energy range on light (Z ≤ 6) targets. This model features four empirical parameters that have been fit to reproduce experimental breakup measurements on beryllium targets, using the method of least squares. It was shown that these parameters extrapolate well to higher energies, and to other low-Z target materials. Here, we also include optimization of the parameters that modify the Kalbach systematics for compound and preequilibrium reactions, in order to better reproduce the experimental data for beryllium targets at large angles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Density dependence of the symmetry energy in the post–PREX-CREX era

The recently published CREX results suggest a rather peculiar picture for the density dependence of the symmetry energy. Whereas PREX favors a large neutron skin thickness in 208 Pb, thereby suggesting a stiff equation of state, CREX suggests instead a much softer equation of state. Here, this discrepancy has caused a large spur in the theoretical community, since no model has been able to simultaneously reproduce within 1σ the PREX and CREX results. Motivated by a novel correlation between a CREX observable and a combination of bulk symmetry energy parameters, we calibrate three new covariant energy density functionals that reproduce binding energies and charge radii of spherical nuclei—and also accommodate the constraints imposed by PREX and CREX. Given that these models suggest a stiff equation of state at high densities, predictions for neutron star properties are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Single-channel and single-energy partial-wave analysis with continuity improved through minimal phase constraints

Single-energy partial-wave analysis has often been applied as a way to fit data with minimal model dependence. However, remaining unconstrained, partial waves at neighboring energies will vary discontinuously because the overall amplitude phase cannot be determined through single-channel measurements. This problem can be mitigated through the use of a constraining penalty function based on an associated energy-dependent fit. However, the weight given to this constraint results in a biased fit to the data. In this paper, for the first time, we explore a constraining function which does not influence the fit to data. The constraint comes from the overall phase found in multichannel fits which, in the present study, are the Bonn-Gatchina and Jülich-Bonn multichannel analyses. The data are well reproduced and weighting of the penalty function does not influence the result. The method is applied to K⁢Λ photoproduction data and all observables can be maximally well reproduced. While the employed multichannel analyses display very different multipole amplitudes, we show that the major difference between two sets of multipoles can be related to the different overall phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Magnetic dipole γ-ray strength functions in the crossover from spherical to deformed neodymium isotopes

We calculate the magnetic dipole $\gamma$-ray strength functions in a chain of even-mass neodymium isotopes $^{144-152}$Nd in the framework of the configuration-interaction (CI) shell model. We infer the strength function by applying the maximum entropy method (MEM) to the exact imaginary-time response function calculated with the shell-model Monte Carlo (SMMC) method. The success of the MEM depends on the choice of a good strength function as a prior distribution. We investigate two choices for the prior strength function: the static path approximation (SPA) and the quasiparticle random-phase approximation (QRPA). We find that the QRPA is a better approximation at low temperatures (i.e., near the ground state), while the SPA is a better choice at finite temperatures. We identify a low-energy enhancement (LEE) in the MEM deexcitation $M1$ strength functions of the even-mass neodymium isotopes and compare with recent experimental results for the total deexcitation $\gamma$-ray strength functions. The LEE is already seen in the SPA strength function but not in the QRPA strength function, indicating the importance of large-amplitude static fluctuations around the mean field in reproducing the LEE. Our method is currently the only one which can reproduce LEE in heavy open-shell nuclei where conventional CI shell model calculations are prohibited. With the onset of deformation as number of neutrons increases along the chain of neodymium isotopes, we observe that some of the LEE strength transfers to a low-energy excitation, which we interpret as a finite-temperature ``scissors'' mode. Here, we also observe a finite-temperature spin-flip mode.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Negative components in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ of heavy-ion fusion cross sections

Negative components are observed in the representation $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ in heavy-ion fusions related to recent studies of noticeable oscillations in the high-energy region of the fusion excitation functions. A modified multiGaussian model has been developed, which includes contributions for both positive and negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ representation. The model reproduces the heavy-ion fusion excitation functions very well, including oscillations. The property of these negative components in the $\frac{𝑑^{2}⁡(𝜎⁢𝐸)}{𝑑⁡𝐸^2}$ spectrum and the reaction mechanism of fusion at energies above the Coulomb barrier are investigated. Finally, the model also reproduces well light heavy-ion fusion reactions between 12 C and 16 O at high energy; results imply that the compound nucleus channel effect may have an important contribution at the high energy region for all fusion systems.

low & intermediate energy heavy-ion reactions↗
DOE OSTI2018

Quantum Monte Carlo formalism for dynamical pions and nucleons

In most simulations of nonrelativistic nuclear systems, the wave functions found solving the many-body Schrodinger equations describe the quantum-mechanical amplitudes of the nucleonic degrees of freedom. In those simulations the pionic contributions are encoded in nuclear potentials and electroweak currents, and they determine the low-momentum behavior. In this work we present an alternative quantum Monte Carlo formalism in which both relativistic pions and nonrelativistic nucleons are explicitly included in the quantum-mechanical states of the system. We report the renormalization of the nucleon mass as a function of the momentum cutoff, a Euclidean time density correlation function that deals with the short-time nucleon diffusion, and the pion cloud density and momentum distributions. In the two-nucleon sector we show that the interaction of two static nucleons at large distances reduces to the one-pion exchange potential, and we fit the low-energy constants of the contact interactions to reproduce the binding energy of the deuteron and two neutrons in finite volumes. We show that the method can be readily applied to light-nuclei.

Madeira, Lucas↗
DOE OSTI2018

Microscopic optical potentials for calcium isotopes

We construct nucleonic microscopic optical potentials by combining the Green's function approach with the coupled-cluster method for 40 Ca and 48 Ca. For the computation of the ground states of 40 Ca and 48 Ca, we use the coupled-cluster method in the singles-and-doubles approximation, while for the A = ±1 nuclei we use particle-attached and particle-removed equation-of-motion method truncated at two-particle–one-hole and one-particle–two-hole excitations, respectively. Our calculations are based on the chiral nucleon-nucleon and three-nucleon interaction NNLO sat , which reproduces the charge radii of 40 Ca and 48 Ca, and the chiral nucleon-nucleon interaction NNLO opt . In all cases considered here, we observe that the overall form of the neutron scattering cross section is reproduced for both interactions, but the imaginary part of the potential, which reflects the loss of flux in the elastic channel, is negligible. Furthermore, the latter points to neglected many-body correlations that would appear beyond the coupled-cluster truncation level considered in this work. We show that, by artificially increasing the parameter η in the Green's function, practical results can be further improved.

74 ATOMIC AND MOLECULAR PHYSICS↗