Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “APPROXIMATIONS”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Progress towards understanding ultranonlocality through the wave-vector and frequency dependence of approximate exchange-correlation kernels

In the framework of time-dependent density functional theory (TDDFT), the exact exchange-correlation (xc) kernel f xc (n, q, ω) determines the ground-state energy, excited-state energies, lifetimes, and the time-dependent linear density response of any many-electron system. The recently developed MCP07 xc kernel f xc (n, q, ω) of A. Ruzsinszky et al. [Phys. Rev. B 101, 245135 (2020)] yields excellent uniform electron gas (UEG) ground-state energies and plausible plasmon lifetimes. As MCP07 is constructed to describe f xc of the UEG, it cannot capture optical properties of real materials. To verify this claim, we follow Nazarov et al. [Phys. Rev. Lett. 102, 113001 (2009)] to construct the long-range, dynamic xc kernel, lim q→0 f xc (n, q, ω) = -α(ω)e 2 /q 2 , of a weakly inhomogeneous electron gas, using MCP07 and other common xc kernels. The strong wavevector and frequency dependence of the “ultranonlocality” coefficient α(ω) is demonstrated for a variety of simple metals and semiconductors. We examine how imposing exact constraints on an approximate kernel shapes α(ω). Comparisons to kernels derived from correlated-wavefunction calculations are drawn.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhancing T c in a composite superconductor/metal bilayer system: A dynamical cluster approximation study

It has been proposed that the superconducting transition temperature T c of an unconventional superconductor with a large pairing scale but strong phase fluctuations can be enhanced by coupling it to a metal. However, the general efficacy of this approach across different parameter regimes remains an open question. Using the dynamical cluster approximation, we study this question in a system composed of an attractive Hubbard layer in the intermediate coupling regime, where the magnitude of the attractive Coulomb interaction |U| is slightly larger than the bandwidth W, hybridized with a noninteracting metallic layer. We find that while the superconducting transition becomes more mean-field-like with increasing interlayer hopping, the superconducting transition temperature T c exhibits a nonmonotonic dependence on the strength of the hybridization t ⊥ . Furthermore, this behavior arises from a reduction of the effective pairing interaction in the correlated layer that outcompetes the growth in the intrinsic pair-field susceptibility induced by the coupling to the metallic layer. We find that the largest T c inferred here for the composite system is comparable to the maximum value currently estimated for the isolated negative-U Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Quasielastic lepton scattering and back-to-back nucleons in the short-time approximation

Understanding quasielastic electron- and neutrino-scattering from nuclei has taken on new urgency with current and planned neutrino oscillation experiments, and with electron scattering experiments measuring specific final states, such as those involving nucleon pairs in “back-to-back” configurations. Accurate many-body methods are available for calculating the response of light (A ≤ 12) nuclei to electromagnetic and weak probes, but they are computationally intensive and only applicable to the inclusive response. In the present work we introduce a novel approach, based on realistic models of nuclear interactions and currents, to evaluate the short-time (high-energy) inclusive and exclusive response of nuclei. The approach accounts reliably for crucial two-nucleon dynamics, including correlations and currents, and provides information on back-to-back nucleons observed in electron and neutrino scattering experiments. We demonstrate that in the quasielastic regime and at moderate momentum transfers both initial- and final-state correlations, and two-nucleon currents are important for a quantitatively successful description of the inclusive response and final state nucleons. Finally, the approach can be extended to include relativistic—kinematical and dynamical—effects, at least approximately in the two-nucleon sector, and to describe the response in the resonance-excitation region.

6 ≤ A ≤ 19↗
DOE OSTI2022

β — -delayed fission in the coupled quasiparticle random-phase approximation plus Hauser-Feshbach approach

Beta-delayed neutron emission and β-delayed fission (βdf) probabilities were calculated for heavy, neutronrich nuclei using the Los Alamos coupled quasiparticle random-phase approximation plus Hauser-Feshbach (QRPA + HF) approach. In this model, the compound nucleus is initially populated by β decay and is followed through subsequent statistical decays taking into account competition among neutrons, γ rays, and fission. The primary output of these calculations includes branching ratios along with neutron and γ -ray spectra. Here, we find a relatively large region of heavy nuclides where the probability of βdf is near 100%. For a subset of nuclei near the neutron drip line, delayed neutron emission and the probability of fission are both large, which leads to the possibility of multichance βdf (mc-βdf). We comment on prospective neutron-rich nuclei that could be probed by future experimental campaigns and provide a full table of branching ratios in ASCII format in the Supplemental Material for use in various applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Level scheme study of Mo 91 : Weak-coupling approximation in the N = 50 region

The level scheme of 91 Mo has been studied using the GALILEO γ-ray spectrometer at Laboratori Nazionali di Legnaro, INFN. Five transitions together with four levels are identified for the first time by the γ-γ coincidence measurement, considerably enriching the negative band in 91 Mo. Furthermore, the discrepancy in the previous parity assignments for the states higher than 3807.5 keV has been clarified with the help of new experimental results. A shell-model calculation employing the JUN45 interaction is performed in order to get a further understanding of the level structure of 91 Mo. Being close to N=50 shell closure, the weak-coupling approximation between the valence g 9/2 neutron-hole with the even-even core ( 92 Mo) is discussed in 91Mo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Quasiparticle random-phase approximation calculations for M 1 transitions with the noniterative finite-amplitude method and application to neutron radiative capture cross sections

Here, we derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method with the Hartree-Fock+Bardeen–Cooper–Schrieffer (HF+BCS) single-particle states, and calculate the magnetic dipole (M⁢1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that would correspond to the M⁢1 scissors mode observed in nuclear experiments. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the photoabsorption cross sections of even-even nuclei given by QRPA. We find that the capture cross section is enhanced due to the contribution from the low energy M⁢1 transition although the calculated capture cross section still underestimates the experimental data. This issue in the calculated capture cross section could be improved by uncertainties of the low energy E⁢1 transition neglected in our QRPA calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Landscape of nuclear deformation softness with spherical quasiparticle random-phase approximation

We investigate the stability and softness of nuclei against quadrupole, octupole, and hexadecapole deformation. By applying the spherical Skyrme-force Hartree-Fock Bardeen-Cooper-Schrieffer quasiparticle random phase approximation, we diagnose ground-state deformation when imaginary solutions are obtained, i.e., the spherical ground state collapses. We also calculate the multipole polarizability in spherical nuclei with no collapse, as a measure of softness. This numerically light and theoretically sound method is found able to capture deformation patterns across the nuclide chart. As a result, the connection between the intrinsic shape of nuclei and the dynamics of their low-lying collective states is established and the role of shell structure is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Quantum Monte Carlo calculations of electron scattering from 12 C in the short-time approximation

The short-time approximation is a method introduced to evaluate electroweak nuclear response for systems with A ≥ 12, extending the reach of first-principle many-body quantum Monte Carlo calculations. Using realistic two- and three-body nuclear interactions and consistent one- and two-body electromagnetic currents, we calculate longitudinal and transverse response densities and response functions of 12 C. Here, we compare the resulting cross sections with experimental data for electron-nucleus scattering, finding good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Transport coefficients of quasiparticle models within a new relaxation time approximation of the Boltzmann equation

We investigate the transport properties of a kinetic theory model that is tuned to describe the thermodynamic properties of QCD at zero chemical potential using a new formulation of the relaxation time approximation. In contrast to previous approaches, the latter is constructed to preserve the fundamental properties of the collision term of the Boltzmann equation for any energy dependence of the relaxation time. A novel choice of matching conditions is implemented to ensure that the background mean-field depends only on the temperature even when the system is out of equilibrium. We provide a consistent analysis of how the transport coefficients of relativistic Navier-Stokes theory vary with the energy dependence of the relaxation time. We also show that the entropy production of this theory is consistent with the second law of thermodynamics and verify that it is independent of the matching conditions employed. We used this fact to calculate the matching independent combination of transport coefficients.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Gamow-Teller Strength in Ca 48 and Ni 78 with the Charge-Exchange Subtracted Second Random-Phase Approximation

Here, we develop a fully self-consistent subtracted second random-phase approximation for charge-exchange processes with Skyrme energy-density functionals. As a first application, we study Gamow-Teller excitations in the doubly magic nucleus 48 Ca, the lightest double-β emitter that could be used in an experiment, and in 78 Ni, the single-beta-decay rate of which is known. The amount of Gamow-Teller strength below 20 or 30 MeV is considerably smaller than in other energy-density-functional calculations and agrees better with experiment in 48 Ca, as does the beta-decay rate in 78 Ni. These important results, obtained without ad hoc quenching factors, are due to the presence of two-particle–two-hole configurations. Their density progressively increases with excitation energy, leading to a long high-energy tail in the spectrum, a fact that may have implications for the computation of nuclear matrix elements for neutrinoless double-β decay in the same framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Optimal Protocols in Quantum Annealing and Quantum Approximate Optimization Algorithm Problems

Quantum annealing (QA) and the quantum approximate optimization algorithm (QAOA) are two special cases of the following control problem: apply a combination of two Hamiltonians to minimize the energy of a quantum state. Which is more effective has remained unclear. Here we analytically apply the framework of optimal control theory to show that generically, given a fixed amount of time, the optimal procedure has the pulsed (or “bang-bang”) structure of QAOA at the beginning and end but can have a smooth annealing structure in between. This is in contrast to previous works which have suggested that bang-bang (i.e., QAOA) protocols are ideal. To support this theoretical work, we carry out simulations of various transverse field Ising models, demonstrating that bang-anneal-bang protocols are more common. Futher, the general features identified here provide guideposts for the nascent experimental implementations of quantum optimization algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Particle Dynamics in a Diblock-Copolymer-Based Dodecagonal Quasicrystal and Its Periodic Approximant by X-Ray Photon Correlation Spectroscopy

Temperature-dependent X-ray photon correlation spectroscopy (XPCS) measurements are reported for a binary diblock copolymer blend that self-assembles into an aperiodic dodecagonal quasicrystal (DDQC) and a periodic Frank-Kasper σ phase approximant. The measured structural relaxation times are Bragg scattering wavevector-independent and are five times faster in the DDQC than the σ phase, with minimal temperature dependence. Further, the underlying dynamical relaxations are ascribed to differences in particle motion at the grain boundaries within each of these tetrahedrally close-packed assemblies. These results identify unprecedented particle dynamics measurements of tetrahedrally-coordinated micellar block polymers, thus expanding the application of XPCS to ordered soft materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Integer versus half-integer spin on an approximate honeycomb lattice

Recent interest in honeycomb lattice materials has focused on their potential to host quantum spin liquid (QSL) states. Variations in bond angles and spin allow a range of interesting behaviors on this lattice, from the predicted QSL ground state of the Kitaev model to exotic magnetic orders. Here we report the physical properties of two compounds with rare earths on an approximate honeycomb lattice. Here, the isostructural compounds Nd 2 S 5 Sn (J = 9/2) and Pr 2 S 5 Sn (J = 4) permit a direct comparison of half-integer versus integer spins on this lattice. We find strikingly different magnetic properties for the two compounds. Nd 2 S 5 Sn orders antiferromagnetically at T N ≈ 2.5 K and undergoes several magnetic transitions to other ordered states under applied field. Pr 2 S 5 Sn displays no magnetic ordering transition above T = 0.41 K, and may be proximate to a spin liquid state.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gaussian approximation potential for amorphous Si : H

Hydrogenation of amorphous silicon (a–Si : H) is critical for reducing defect densities, passivating midgap states and surfaces, and improving photoconductivity in silicon-based electro-optical devices. Modeling the atomic-scale structure of this material is critical to understanding these processes, which in turn is needed to describe c–Si/a–Si : H heterojunctions that are at the heart of modern solar cells with world-record efficiency. Density functional theory (DFT) studies achieve the required high accuracy but are limited to moderate system sizes of 100 atoms or so by their high computational cost. Simulations of amorphous materials have been hindered by this high cost because large structural models are required to capture the medium-range order that is characteristic of such materials. Empirical potential models are much faster, but their accuracy is not sufficient to correctly describe the frustrated local structure. Data-driven, machine-learned interatomic potentials have broken this impasse and have been highly successful in describing a variety of amorphous materials in their elemental phase. Here, we extend the Gaussian approximation potential (GAP) for silicon by incorporating the interaction with hydrogen, thereby significantly improving the degree of realism with which amorphous silicon can be modeled. We show that our Si : H GAP enables the simulation of hydrogenated silicon with an accuracy very close to DFT but with computational expense and run times reduced by several orders of magnitude for large structures. Here, we demonstrate the capabilities of the Si : H GAP by creating models of hydrogenated liquid and amorphous silicon and showing that their energies, forces, and stresses are in excellent agreement with DFT results, and their structure as captured by bond and angle distributions are in agreement with both DFT and experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

In this work, we derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of sp and sd metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r 2 SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r 2 SCAN. We compare OFR2 to the r 2 SCAN -L LL-MGGA [Mejia-Rodriguez and Trickey, Phys. Rev. B 102, 121109(R) (2020)] and show that OFR2 tends to outperform r 2 SCAN -L for the equilibrium properties of solids, but r 2 SCAN -L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and nonmetallic condensed matter, we continue to recommend SCAN and r 2 SCAN. Numerical performance is discussed in detail, and our paper provides an outlook to machine learning.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Active learning approach to simulations of strongly correlated matter with the ghost Gutzwiller approximation

Quantum embedding (QE) methods such as the ghost Gutzwiller approximation (gGA) offer a powerful approach to simulating strongly correlated systems, but come with the computational bottleneck of computing the ground state of an auxiliary embedding Hamiltonian (EH) iteratively. In this work, we introduce an active learning (AL) framework integrated within the gGA to address this challenge. The methodology is applied to the single-band Hubbard model and results in a significant reduction in the number of instances where the EH must be solved. Through a principal component analysis (PCA), we find that the EH parameters form a low-dimensional structure that is largely independent of the geometric specifics of the systems, especially in the strongly correlated regime. Our AL strategy enables us to discover this low-dimensionality structure on the fly, while leveraging it for reducing the computational cost of gGA, laying the groundwork for more efficient simulations of complex strongly correlated materials. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗