Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correlation energy”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Re-examining the silicon self-interstitial charge states and defect levels: A density functional theory and bounds analysis study

The self-interstitial atom (SIA) is one of two fundamental point defects in bulk Si. Isolated Si SIAs are extremely difficult to observe experimentally. Even at very low temperatures, they anneal before typical experiments can be performed. Given the challenges associated with experimental characterization, accurate theoretical calculations provide valuable information necessary to elucidate the properties of these defects. Previous studies have applied Kohn–Sham density functional theory (DFT) to the Si SIA, using either the local density approximation or the generalized gradient approximation to the exchange-correlation (XC) energy. The consensus of these studies indicates that a Si SIA may exist in five charge states ranging from -2 to +2 with the defect structure being dependent on the charge state. This study aims to re-examine the existence of these charge states in light of recently derived “approximate bounds” on the defect levels obtained from finite-size supercell calculations and new DFT calculations using both semi-local and hybrid XC approximations. We conclude that only the neutral and +2 charge states are directly supported by DFT as localized charge states of the Si SIA. Within the current accuracy of DFT, our results indicate that the +1 charge state likely consists of an electron in a conduction-band-like state that is coulombically bound to a +2 SIA. Furthermore, the -1 and -2 charge states likely consist of a neutral SIA with one and two additional electrons in the conduction band, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Convergence acceleration of Monte Carlo many-body perturbation methods by using many control variates

Herein, the use of many control variates is proposed as a method to accelerate the second- and third-order Monte Carlo (MC) many-body perturbation (MC-MP2 and MC-MP3) calculations. A control variate is an exactly integrable function that is strongly correlated or anti-correlated with the target function to be integrated by the MC method. Evaluating both integrals as well as their covariances in the same MC run, one can effect a mutual cancellation of the statistical uncertainties and biases in the MC integrations, thereby accelerating its convergence considerably. Six and thirty-six control variates, whose integrals are known a priori, are generated for MC-MP2 and MC-MP3, respectively, by systematically replacing one or more two-electron-integral vertexes of certain configurations by zero-valued overlap-integral vertexes in their Goldstone diagrams. The variances and co- variances of these control variates are computed at a marginal cost, enhancing the overall efficiency of the MC-MP2 and MC-MP3 calculations by a factor of up to 14 and 20, respectively.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Piecewise interaction picture density matrix quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Fermionic mean-field theory as a tool for studying spin Hamiltonians

The Jordan–Wigner transformation permits one to convert spin 1/2 operators into spinless fermion ones, or vice versa. In some cases, it transforms an interacting spin Hamiltonian into a noninteracting fermionic one, which is exactly solved at the mean-field level. Even when the resulting fermionic Hamiltonian is interacting, its mean-field solution can provide surprisingly accurate energies and correlation functions. Furthermore, Jordan–Wigner is, however, only one possible means of interconverting spin and fermionic degrees of freedom. Here, we apply several such techniques to the XXZ and J 1 –J 2 Heisenberg models, as well as to the pairing or reduced Bardeen–Cooper–Schrieffer Hamiltonian, with the aim of discovering which of these mappings is most useful in applying fermionic mean-field theory to the study of spin Hamiltonians.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Intensified Tpx3Cam, a fast data-driven optical camera with nanosecond timing resolution for single photon detection in quantum applications

Here we describe a fast data-driven optical camera, Tpx3Cam, with nanosecond scale timing resolution and 80 Mpixel/sec throughput. After the addition of intensifier, the camera is single photon sensitive with quantum efficiency determined primarily by the intensifier photocathode. The single photon performance of the camera was characterized with results on the gain, timing resolution and afterpulsing reported here. The intensified camera was successfully used for measurements in a variety of applications including quantum applications. As an example of such application, which requires simultaneous detection of multiple photons, we describe registration of photon pairs from the spontaneous parametric down-conversion source in a spectrometer. We measured the photon wavelength and timing with respective precisions of 0.15 nm and 3 ns, and also demonstrated that the two photons are anti-correlated in energy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Weak structure functions in $\nu_l-N$ and $\nu_l-A$ scattering with nonperturbative and higher order perturbative QCD effects

We study the effect of various perturbative and nonperturbative QCD corrections on the free nucleon structure functions [FiNWI(x,Q2);i=1–3] and their implications in the determination of nuclear structure functions. The evaluation of the nucleon structure functions has been performed by using the MMHT 2014 parton distribution functions (PDFs) parametrization, and the target mass correction (TMC) and higher twist (HT) effects are incorporated following the works of Schienbein et al. and Dasgupta et al., respectively. These nucleon structure functions are taken as input in the determination of nuclear structure functions. The numerical calculations for the νl/ν¯l-A deep inelastic scattering (DIS) process have been performed by incorporating the nuclear medium effects like Fermi motion, binding energy, nucleon correlations, mesonic contributions, shadowing, and antishadowing in several nuclear targets such as carbon, polystyrene scintillator, iron, and lead, which are being used in MINERνA, and in argon nuclei, which is relevant for the ArgoNeuT and DUNE experiments. The differential scattering cross sections d2σAWIdxdy and (dσAWIdx/dσCHWIdx) have also been studied in the kinematic region of the MINERνA experiment. The theoretical results are compared with the recent experimental data of MINERνA and the earlier data of the NuTeV, CCFR, CDHSW, and CHORUS Collaborations. Moreover, a comparative analysis of the present results for the ratio (dσAWIdx/dσCHWIdx), and the results from the Monte Carlo (MC) generator GENIE and other phenomenological models of Bodek and Yang, and Cloet et al., has been performed in the context of the MINERνA experiment. The predictions have also been made for the ν¯l-A cross section relevant for the MINERνA experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Single-Shot Diagnosis of Electron Energy Evolution via Streaked Betatron X Rays in a Curved Laser Wakefield Accelerator

We report on an experimental observation of the streaking of betatron x rays in a curved laser wakefield accelerator. The streaking of the betatron x rays was realized by launching a laser pulse into a plasma with a transverse density gradient. By controlling the plasma density and the density gradient, we realized the steering of the laser driver, electron beam, and betatron x rays simultaneously. Furthermore, we observed an energy-angle correlation of the streaked betatron x rays and utilized it in diagnosing the electron acceleration process in a single-shot mode. Our work could also find applications in advanced control of laser beam and particle propagation. More importantly, the angular streaked betatron x ray has an intrinsic spatiotemporal correlation, which makes it a promising tool for single-shot pump-probe applications.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Transitions in the morphology and critical stresses of gliding dislocations in multiprincipal element alloys

Refractory multiprincipal element alloys (MPEAs) are promising material candidates for high-temperature, high-strength applications. However, their deformation mechanisms, particularly at the dislocation scale, are unusual compared to those of conventional alloys, making it challenging to understand the origin of their high strength. Here, using an atomistically informed phase-field dislocation dynamics model, we study transitions in the morphology and critical stresses of long gliding screw dislocations over extended distances in a refractory MPEA. The model MPEA crystal accounts for the atomic-scale fluctuations in chemical composition across the glide planes via spatially correlated lattice energies and for the differences in glide resistances between screw and edge dislocations of unit length. Here, we show that the dislocation moves in a stop-start motion, alternating between wavy morphology in free flight and nearly recovered straight screw orientation in full arrest. The periods of wavy glide are due to variable kink-pair formation and migration rates along the length, where portions with higher rates glide more quickly. The critical stress to initiate motion corresponds to the stress required to form and migrate a kink pair at the weakest region along the length of the dislocation. Heterogeneity in lattice energy leads to variability in the local stress-strain response and to a strain hardening-like response, in which the critical stress to reactivate glide increases with glide distance. Statistical assessment of hundreds of realizations of dislocations indicates that the amount of hardening directly scales with the dispersion in underlying lattice energy.

3-dimensional systems↗
DOE OSTI2024

Quantum chaos on edge

Recently, the physics of many-body quantum chaotic systems close to their ground states has come under intensified scrutiny. Such studies are motivated by the emergence of model systems exhibiting chaotic fluctuations throughout the entire spectrum [the Sachdev-Ye-Kitaev (SYK) model being a renowned representative] as well as by the physics of holographic principles, which likewise unfold close to ground states. Interpreting the edge of the spectrum as a quantum critical point, here we combine a wide range of analytical and numerical methods to the identification and comprehensive description of two different universality classes: the near edge physics of “sparse” and the near edge of “dense” chaotic systems. The distinction lies in the ratio between the number of a system's random parameters and its Hilbert space dimension, which is exponentially small or algebraically small in the sparse and dense case, respectively. Notable representatives of the two classes are generic chaotic many-body models (sparse) and invariant random matrix ensembles or chaotic gravitational systems (dense). While the two families share identical spectral correlations at energy scales comparable to the level spacing, the density of states and its fluctuations near the edge are different. Considering the SYK model as a representative of the sparse class, we apply a combination of field theory and exact diagonalization to a detailed discussion of its edge spectrum. Conversely, Jackiw-Teitelboim gravity is our reference model for the dense class, where an analysis of the gravitational path integral and random matrix theory reveal universal differences to the sparse class, whose implications for the construction of holographic principles we discuss. Published by the American Physical Society 2024

Altland, Alexander (ORCID:0000000229914805)↗
DOE OSTI2021

Artificial intelligence “sees” split electrons

Chemical bonds between atoms are stabilized by the exchange-correlation (xc) energy, a quantum-mechanical effect in which “social distancing” by electrons lowers their electrostatic repulsion energy. Kohn-Sham density functional theory (DFT) states that the electron density determines this xc energy, but the density functional must be approximated. Furthermore, this is usually done by satisfying exact constraints of the exact functional (making the approximation predictive), by fitting to data (making it interpolative), or both. Two exact constraints—the ensemble-based piecewise linear variation of the total energy with respect to fractional electron number and fractional electron z-component of spin —require hard-to-control nonlocality. On page 1385 of this issue, Kirkpatrick et al. have taken a big step toward more accurate predictions for chemistry through the machine learning of molecular data plus the fractional charge and spin constraints, expressed as data that a machine can learn.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Nuclear Data Covariances and Sampling

In this report, we present an overview of the task of producing sample nuclear data sets for the purpose of exploring the effects of nuclear data uncertainties on physics simulations. These sample data sets (or, as we sometimes call them variations) are produced based on uncertainties specified in the ENDF/B-VIII nuclear data library, although the principles discussed here apply to sampling based on similar libraries, such as JENDL. The uncertainties in ENDF and similar libraries are given in the form of mean values (a mean vector) and a covariance matrix. The components of the mean vector are what are normally thought of as the values of the nuclear data (cross sections, $\bar{v}$, the prompt fission neutron spectrum, etc.) and the covariance matrix contains the information about the uncertainties in these values and correlations between them, including correlations across different energies in a single channel and correlations across different channels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Physics-informed Estimation of the Covariance Matrix for Various Neutron Spectra

A method for estimating covariance matrices which capture the uncertainties in calculated reactor spectra has been developed. This method is based on perturbing the parameters of a physics-based analytic model fitted to a calculated spectrum. The covariance of the perturbed analytic spectra imposes energy-dependent correlations due to the physics of the neutron processes in the reactor, i.e., a fission component, a 1/E down-scatting component, and a thermal Maxwellian component. An analytic model is developed which is shown to produce good fits to several reactor environments. The covariance matrices produced via this method are then used as the prior spectrum in STAYSL least squares spectrum adjustment where it is combined with integral metrics, such as activation measurements, to produce a high-fidelity neutron spectrum characterization. It was concluded that the methodology showed agreeable results for the ACRR free-field spectrum adjustment in STAYSL resulting in a 𝜒 2 value of 2.21 (per degree of freedom), but further work is needed to describe scattering and interface regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1975

Compressor and fan wake characteristics

Approaches for developing an analytical model capable of determining the effects of rotor flow and blade parameters and turbulence properties (i.e. energy, velocity correlations, and length scale) on the rotor wake characteristics and its diffusion properties are discussed. The three-dimensional model will employ experimental measurements, instantaneous velocities, and turbulence properties at various stations downstream from a rotor. A triaxial probe and a rotating conventional probe, which is mounted on a traverse gear operated by two step motors, are to be used for these measurements. The final rotor wake model will be capable of predicting the discrete and broadband noise generated in a fan rotor and of evaluating the aerodynamic losses, efficiency and optimum spacing between a rotor and stator in turbomachinery.

Lakshminarayana, B.↗
NASA NTRS1975

Change in the high-energy radiation from the Crab

On 1973 July 23, the Cornell gamma-ray telescope had a second exposure to the Crab Nebula and pulsar via balloon flight from Palestine, Texas. On the previous flight, 1971 October 6, the signal from the pulsar had been unmistakably clear, especially at the highest energies (above 800 MeV), and an unpulsed component of similar average power had also been discernible. On the second flight, despite higher altitude, longer exposure, and improved electronics, the pulsed signal was barely detectable and the dc component not observable at all. The drop in intensity seems to increase with energy, being only a factor of 2 at 200 MeV but an order of magnitude at 1 GeV. There is some indication (not compelling) that the gamma-ray flux may have changed even within the 6-hour exposure. The first flight may have viewed the Crab in a state of enhanced high-energy emission correlated with the glitches that occurred on 1971 August 1 and October 25.

Greisen, K.↗
NASA NTRS1978

Compressor and fan wake characteristics

A triaxial probe and a rotating conventional probe, mounted on a traverse gear operated by two step motors were used to measure the mean velocities and turbulence quantities across a rotor wake at various radial locations and downstream stations. The data obtained was used in an analytical model developed to study how rotor flow and blade parameters and turbulence properties such as energy, velocity correlations, and length scale affect the rotor wake characteristics and its diffusion properties. The model, includes three dimensional attributes, can be used in predicting the discrete as well as broadband noise generated in a fan rotor, as well as in evaluating the aerodynamic losses, efficiency and optimum spacing between a rotor and stator in turbomachinery.

Reynolds, B.↗