Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXTRAPOLATION”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Valence and magnetism in EuPd 3 S 4 and (Y, La) x Eu 1-x Pd 3 S 4

151 Eu Mössbauer spectroscopy shows that yttrium substitution in mixed-valent EuPd 3 S 4 drives the initial 50:50 mix of Eu 3+ and E 2+ towards pure Eu 2+ , whereas lanthanum substitution has the opposite effect, but only for substitution levels above 50%. Here, we find that total valence electron count and chemical pressure effects cannot account for the observed behavior; however, conserving the cell volume provides a consistent description of the changes in the Eu 2+ :Eu 3+ ratio. Remarkably, lanthanum substitution also leads to a clear transition from static mixed-valent behavior at lower temperatures to dynamic mixed-valent behavior at higher temperatures, with the onset temperature monotonically increasing with Eu content and extrapolating to a value of ~ 340 K for the pure EuPd 3 S 4 compound. Magnetic order persists at least as far as x=0.875 in both series, despite the drastic reduction in the amount of moment-carrying Eu 2+ ions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

How quantum phases on cylinders approach the two-dimensional limit

We consider the properties of $T$ = $0$ quantum phases of matter - especially superconducting and analogous spin-liquid phases - on infinite cylinders of width $L$ ⊥ and analyze the ways in which the $L$ ⊥ → ∞ (two-dimensional, 2d) limit is approached. This problem is interesting in its own right but is particularly important in the context of extrapolating accessible density matrix renormalization group (DMRG) results on model strongly interacting problems to the desired 2d limit. In conclusion, various methods for drawing firm conclusions about the quantum phases in 2d from relatively small $L$ ⊥ results are illustrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Superconducting valence bond fluid in lightly doped eight-leg $t-J$ cylinders

Superconductivity in doped quantum paramagnets has been a subject of long theoretical inquiry. In this work, we report a density matrix renormalization group study of lightly doped $t-J$ models on the finite-width square lattice (doped hole densities $δ$ = 1/12 and 1/8) with parameters for which previous studies have suggested that the undoped system in 2D is either a quantum spin liquid or a valence bond crystal. Our studies are performed on cylinders with width up to 8. Here, ground-state correlations are found to be nearly identical for the “doped quantum spin liquid” and “doped valence bond crystal.” Upon increasing the cylinder widths from 4 to 8, we observed a significant strengthening of the quasi-long-range superconducting correlations and a dramatic suppression of any “competing” charge density wave order. Extrapolating from the observed behavior of the width eight cylinders, we speculate that the system has a nodeless $d$-wave superconducting ground state in the 2D limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sound speed and Grüneisen parameter up to three terapascal in shock-compressed iron

This paper presents the first sound speed and Grüneisen parameter data for fluid iron compressed to 3 TPa (30 million atmospheres) and 20 g/cm 3 on the Hugoniot. Both the sound speed and Grüneisen parameter are derivatives of the equation of state (EOS), and thus tightly constrain the contours of the EOS surface. The sound speed data are systematically lower than expected from a simple extrapolation of previous data. The Grüneisen parameter shows a 30% drop at pressures and temperatures above the melt transition. Furthermore, while some models compare well with either the sound speed or Grüneisen parameter, none of today’s state-of-the-art models can explain both sets of data. Furthermore these new data will provide pivotal benchmarks for both future theoretical EOSs of warm dense iron and modeling planetary states and processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamical defects in a two-dimensional Wigner crystal: Self-doping and kinetic magnetism

We study the quantum dynamics of interstitials and vacancies in a two-dimensional Wigner crystal (WC) using a semi-classical instanton method that is asymptotically exact at low density, i.e., in the r s → ∞ limit. Here, the dynamics of these point defects mediates magnetism with much higher energy scales than the exchange energies of the pure WC. Via exact diagonalization of the derived effective Hamiltonians in the single-defect sectors, we find the dynamical corrections to the defect energies. The resulting expression for the interstitial (vacancy) energy extrapolates to 0 at r s = r mit ≈ 70 (r s ≈ 30), suggestive of a self-doping instability to a partially melted WC for some range of r s below r mit . We thus propose a “metallic electron crystal” phase of the two-dimensional electron gas at intermediate densities between a low density insulating WC and a high density Fermi fluid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Global polarization of ΛΛ¯ hyperons in Pb-Pb collisions at √ s NN = 2.76 and 5.02 TeV

The global polarization of the Λ and Λ ¯ hyperons is measured for Pb-Pb collisions at s N N = 2.76 and 5.02 TeV recorded with the ALICE at the Large Hadron Collider (LHC). The results are reported differentially as a function of collision centrality and hyperon's transverse momentum ( p T ) for the range of centrality 5–50%, 0.5 < p T < 5 GeV/ c , and rapidity | y | < 0.5 . The hyperon global polarization averaged for Pb-Pb collisions at s N N = 2.76 and 5.02 TeV is found to be consistent with zero, ( P H ) ( % ) ≈ 0.01 ± 0.06 (stat.) ± 0.03 (syst.) in the collision centrality range 15–50%, where the largest signal is expected. The results are compatible with expectations based on an extrapolation from measurements at lower collision energies at the Relativistic Heavy Ion Collider, hydrodynamical model calculations, and empirical estimates based on collision energy dependence of directed flow, all of which predict the global polarization values at LHC energies of the order of 0.01 % .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Measurement of the scintillation and ionization response of liquid xenon at MeV energies in the EXO-200 experiment

In this work, liquid xenon (LXe) is employed in a number of current and future detectors for rare event searches. We use the EXO-200 experimental data to measure the absolute scintillation and ionization yields generated by γ interactions from 228 Th (2615 keV), 226 Ra (1764 keV), and 60 Co (1332 keV and 1173 keV) calibration sources, over a range of electric fields. The W value that defines the recombination-independent energy scale is measured to be 11.5±0.5 (syst.) ±0.1 (stat.) eV. These data are also used to measure the recombination fluctuations in the number of electrons and photons produced by the calibration sources at the MeV scale, which deviate from extrapolations of lower-energy data. Additionally, a semiempirical model for the energy resolution of the detector is developed, which is used to constrain the recombination efficiency, i.e., the fraction of recombined electrons that result in the emission of a detectable photon. Detailed measurements of the absolute charge and light yields for MeV-scale electron recoils are important for predicting the performance of future neutrinoless double β -decay detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Analytic insights on the information content of new observables

Uncertainty quantification has emerged as a rapidly growing field in nuclear science. Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experi- ments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of limited value. Also, one often deals with theoretical constructs that involve fundamental quantities that are not accessible to experiment or observation. Particularly relevant in this context is the pressure of pure neutron matter. Furthermore, in this Rapid Communication, we develop an analytic framework to answer the question of “How can new data reduce uncertainties of current theoretical models?” Simple and insightful expressions are obtained to quantify the impact of one or two new observables on theoretical uncertainties in two critical quantities: The slope of the symmetry energy at saturation density and the pressure of pure neutron matter at twice nuclear matter saturation density.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Observation of new neutron-rich isotopes in the vicinity of Zr 110

Neutron-rich radioactive isotope (RI) beams in the vicinity of 110 Zr were produced using the in-flight fission of a 345-MeV/nucleon 238 U beam at the RIKEN Radioactive Isotope Beam Factory. Here, the RI beams were separated in flight by the large-acceptance two-stage fragment separator BigRIPS. Isotopes were clearly identified in a low-background particle-identification plot, which was obtained from the measurements of time of flight, magnetic rigidity, and energy loss using beam-line detectors placed in the BigRIPS separator and the ZeroDegree spectrometer. Nine new isotopes, 101 Br, 102 Kr, 105,106 Rb, 108 Sr, 110,111 Y, 114 Zr, and 117 Nb, are reported for the first time, with their measured production cross sections compared with the extrapolations of lighter isotopes down to the femtobarn regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Global Λ-hyperon polarization in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV

Global hyperon polarization, $\overline{P}$ H , in Au + Au collisions over a large range of collision energy, $\sqrt{S_{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While na¨ıve extrapolation of data trends suggests a large $\overline{P}$ H as the collision energy is reduced, the behavior of $\overline{P}$ H at small $\sqrt{S_{NN}}$ < 7.7 GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of Λ hyperons along the direction of global angular momentum in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV. As a result, the observation of substantial polarization of 4.91 ± 0.81 (stat.) ± 0.15 (syst.) % in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Converged ab initio calculations of heavy nuclei

We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of E 3max =28, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of 132 Sn converged to within a few MeV with respect to the E 3max truncation. In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum J rel > 9/2 to the ground-state energy is dominated by a basis-truncation artifact, which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of 132 Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Comprehensive simulation of heavy-ion collisions at nonzero baryon chemical potential

Here, we present results of hydrodynamic modeling of Au-Au collisions from $\sqrt{s_{NN}}$ = 7.7 to 200 GeV. Our simulations have three novel components. First, we use a Linear EXtrapolation of Ultratrelativistic nucleon-nucleon Scattering to nucleus-nucleus collisions (LEXUS) inspired Monte Carlo initial-state model. Second, we use a crossover equation of state at finite baryon densities without a critical point. Finally, we use departure functions derived from the quasiparticle theory of transport coefficients for hadronic matter at nonzero baryon densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Relativistic mean-field theories for neutron-star physics based on chiral effective field theory

We describe and implement a procedure for determining the couplings of a relativistic mean-field theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars which consist of neutron-rich matter at densities ranging up to several times saturation density, which allows for additional astrophysical constraints. In our approach, rather than using the RMFT to extrapolate from symmetric to neutron-rich matter and from finite-sized nuclei to uniform matter, we fit the RMFT to properties of uniform pure neutron matter obtained from chiral effective field theory. Chiral effective field theory incorporates the experimental data for nuclei in the framework of a controlled expansion for nuclear forces valid at nuclear densities and enables us to account for theoretical uncertainties when fitting the RMFT. We construct four simple RMFTs that span the uncertainties provided by chiral effective field theory for neutron matter, and are consistent with current astrophysical constraints on the equation of state. Lastly, our RMFTs can be used to model the properties of neutron-rich matter across the vast range of densities and temperatures encountered in neutron stars and their mergers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

3D structure of anisotropic flow in small collision systems at energies available at the BNL Relativistic Heavy Ion Collider

Here, we present (3 + 1)-dimensional [(3 + 1)D] dynamical simulations of asymmetric nuclear collisions at the BNL Relativistic Heavy Ion Collider (RHIC). Employing a dynamical initial state model coupled to (3 + 1)D viscous relativistic hydrodynamics, we explore the rapidity dependence of anisotropic flow in the RHIC small system scan at 200 GeV center-of-mass energy. We calibrate parameters to describe central 3 He + Au collisions and make extrapolations to d + Au and p + Au collisions. Our calculations demonstrate that approximately 50% of the v 3 (p T ) difference between the measurements by the STAR and PHENIX Collaborations can be explained by the use of reference flow vectors from different rapidity regions. This emphasizes the importance of longitudinal flow decorrelation for anisotropic flow measurements in asymmetric nuclear collisions, and the need for (3 + 1)D simulations. We also present results for the beam energy dependence of particle spectra and anisotropic flow in d + Au collisions.

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↗