Search NASA⌕ Search

SEARCH · Search NASA

Results for “Many body perturbation calculations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Relativistic configuration-interaction and many-body-perturbation-theory calculations of U i hyperfine constants

Neutral uranium (U i) is a very difficult atom for theoretical calculations due to a large number (six) of valence electrons, strong valence-valence and valence-core correlations, high density of states, and relativistic effects. Configuration-interaction many-body perturbation theory (CI-MBPT) can efficiently treat valence-core correlations and relativistic effects, but because the formalism was developed for a Dirac-Hartree-Fock (DHF) starting potential that does not contain valence electrons, quite large CI space is needed to compensate for the + 6 charge of such a potential. Much more efficient is the relativistic configuration-interaction (RCI) approach, which uses a relatively accurate starting DHF potential that includes some valence electrons to make the valence-electron Hamiltonian diagonally dominated for some states. Here we report calculations of U i hyperfine constants of several low-energy states using the RCI method with the starting potential that includes four f valence electrons. With this starting potential, it is possible to use the single-configuration approximation or small basis sets to obtain quite accurate results for hyperfine-structure constants. In fact, by scaling the nuclear magnetic moment, the agreement for five levels was within 5% and a new magnetic moment can be recommended, 0.43(2). When two states are mixed, it is difficult to predict the exact mixing fractions, so in this case we proposed a method in which the mixing fraction is found by optimization of one parameter in the MBPT correction to obtain correct g factors. The resulting hyperfine constants are significantly improved. Finally, the CI-MBPT approach was also tested, with some success for low-energy states, with the limited basis set. The methods investigated here can be further developed to include more extensive data sets to improve accuracy and can be applied to other atoms and for calculations of other properties, for example, relevant to fundamental symmetry tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Mixed Stochastic-Deterministic Approach for Many-Body Perturbation Theory Calculations

Here, we present an approach for GW calculations of quasiparticle energies with quasi-quadratic scaling by approximating high-energy contributions to the Green's function in its Lehmann representation with effective stochastic vectors. The method is easy to implement without altering the GW code, converges rapidly with stochastic parameters, and treats systems of various dimensionality and screening response. Our calculations on a 5.75° twisted MoS 2 bilayer show how large-scale GW methods include geometry relaxations and electronic correlations on an equal basis in structurally nontrivial materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
DOE OSTI2022

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. Here, we present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

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

Delocalization of dark and bright excitons in flat-band materials and the optical properties of V2O5

Abstract The simplest picture of excitons in materials with atomic-like localization of electrons is that of Frenkel excitons, where electrons and holes stay close together, which is associated with a large binding energy. Here, using the example of the layered oxide V 2 O 5 , we show how localized charge-transfer excitations combine to form excitons that also have a huge binding energy but, at the same time, a large electron-hole distance, and we explain this seemingly contradictory finding. The anisotropy of the exciton delocalization is determined by the local anisotropy of the structure, whereas the exciton extends orthogonally to the chains formed by the crystal structure. Moreover, we show that the bright exciton goes together with a dark exciton of even larger binding energy and more pronounced anisotropy. These findings are obtained by combining first principles many-body perturbation theory calculations, ellipsometry experiments, and tight binding modelling, leading to very good agreement and a consistent picture. Our explanation is general and can be extended to other materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Magnetically confined surface and bulk excitons in a layered antiferromagnet

The discovery of two-dimensional van der Waals magnets has greatly expanded our ability to create and control nanoscale quantum phases. A unique capability emerges when a two-dimensional magnet is also a semiconductor that features tightly bound excitons with large oscillator strengths that fundamentally determine the optical response and are tunable with magnetic fields. Here, in this study, we report a previously unidentified type of optical excitation-a magnetic surface exciton-enabled by the antiferromagnetic spin correlations that confine excitons to the surface of CrSBr. Magnetic surface excitons exhibit stronger Coulomb attraction, leading to a higher binding energy than excitons confined in bulk layers, and profoundly alter the optical response of few-layer crystals. Distinct magnetic confinement of surface and bulk excitons is established by layer- and temperature-dependent exciton reflection spectroscopy and corroborated by ab initio many-body perturbation theory calculations. By quenching interlayer excitonic interactions, the antiferromagnetic order of CrSBr strictly confines the bound electron-hole pairs within the same layer, regardless of the total number of layers. Our work unveils unique confined excitons in a layered antiferromagnet, highlighting magnetic interactions as a vital approach for nanoscale quantum confinement, from few layers to the bulk limit.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stochastic evaluation of fourth-order many-body perturbation energies

A scalable, stochastic algorithm evaluating the fourth-order many-body perturbation (MP4) correction to energy is proposed. Three hundred Goldstone diagrams representing the MP4 correction are computer generated and then converted into algebraic formulas expressed in terms of Green’s functions in real space and imaginary time. They are evaluated by the direct (i.e., non-Markov, non-Metropolis) Monte Carlo (MC) integration accelerated by the redundant-walker and control-variate algorithms. The resulting MC-MP4 method is efficiently parallelized and is shown to display O(n 5.3 ) size-dependence of cost, which is nearly two ranks lower than the O(n 7 ) dependence of the deterministic MP4 algorithm. Furthermore, it evaluates the MP4/aug-cc-pVDZ energy for benzene, naphthalene, phenanthrene, and corannulene with the statistical uncertainty of 10 mE h (1.1% of the total basis-set correlation energy), 38 mE h (2.6%), 110 mE h (5.5%), and 280 mE h (9.0%), respectively, after about 10 9 MC steps.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗
DOE OSTI2020

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory (χ EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χ EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. Furthermore, we then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium χ EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

Here, we introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (𝜒⁢EFT ). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as 𝜒⁢EFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the 𝜒⁢EFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Transition energy measurements of the X-ray lines of neon-like europium

Here, we report the measurements of the n = 3 → n = 2 L-shell X-ray transitions of neon-like Eu 53+ . These include all seven allowed electric dipole (E1) transitions, five electric quadrupole (E2) transitions, and one magnetic quadrupole (M2) transition. The measurements were carried out at the EBIT-I electron beam ion trap facility at Livermore using an X-ray calorimeter. The measured values agree with many-body perturbation theory calculations within the experimental uncertainties. This includes the transitions that involve a 2s electron and that are, therefore, strongly affected by quantum electrodynamical effects.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Relativistic MR-MP Energy Levels for K-like Fe

Relativistic Multi-Reference Møller-Plesset (MR-MP) many-body perturbation theory calculations are reported for the ground and excited states of Fe viii. Energies are provided for states coming from the configurations $3{s}^{m}3{p}^{n}3{d}^{p}4{l}^{q}$ (l = s, p, d, and f) with $m+n+p+q=9$, p ≤ 3, and q ≤ 1. The MR-MP results were compared with the recommended observed values in the CHIANTI-v8 database and with previous ab-initio calculations. The MR-MP calculations reproduce 90% of the available observed values within 3500 cm –1 . The previous accurate ab-initio calculations reproduce only 50% of the observed values within 3500 cm –1 . Furthermore, the electronic structure of Fe viii is very complex and our MR-MP results show that accurate calculations require the inclusion of not only 3s 2 –3d 2 , 3p 2 –3d 2 , and 3s3p–3d 2 promotions but also of all electron correlations via perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Convergence acceleration of Monte Carlo many-body perturbation methods by direct sampling

In the Monte Carlo many-body perturbation (MC-MP) method, the conventional correlation-correction formula, which is a long sum of products of low-dimensional integrals, is first recast into a short sum of high-dimensional integrals over electron-pair and imaginary-time coordinates. These high-dimensional integrals are then evaluated by the Monte Carlo method with random coordinates generated by the Metropolis–Hasting algorithm according to a suitable distribution. The latter algorithm, while advantageous in its ability to sample nearly any distribution, introduces autocorrelation in sampled coordinates, which in turn increases the statistical uncertainty of the integrals and thus the computational cost. It also involves wasteful rejected moves and an initial “burn-in” step as well as displays hysteresis. Here, an algorithm is proposed that directly produces a random sequence of electron-pair coordinates for the same distribution used in the MC-MP method, which is free from autocorrelation, rejected moves, a burn-in step, or hysteresis. Furthermore, this direct-sampling algorithm is shown to accelerate second- (MC-MP2) and third-order Monte Carlo many-body perturbation (MC-MP3) calculations by up to 222% and 38%, respectively.

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

Dark-exciton energy splitting in monolayer W Se 2 : Insights from time-dependent density functional theory

Here, we present here a formalism based on time-dependent density functional theory to describe characteristics of both intra- and intervalley excitons in semiconductors, the latter of which had remained a challenge. Through the usage of an appropriate exchange-correlation kernel (nanoquanta), we trace the energy difference between the intra- and intervalley dark excitons in monolayer (1L) W⁢ Se 2 to the domination of the exchange part in the exchange-correlation energies of these states. Furthermore, our calculated transition contribution maps establish the momentum resolved weights of the electron-hole excitations in both bright and dark excitons thereby providing a comprehensive understanding of excitonic properties of 1L W⁢ Se 2 . We find that the states consist of hybridized excitations around the corresponding valleys which leads to brightening of the dark excitons, i.e., significantly decreasing their lifetime which is reflected in the PL spectrum. Using many-body perturbation theory, we calculate the phonon contribution to the energy band gap and the linewidths of the excited electrons, holes, and (bright) exciton to find that as the temperature increases the band gap significantly decreases, while the linewidths increase. Our work paves for describing the ultrafast charge dynamics of transition metal dichalcogenide within an ab initio framework.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2022

Finding predictive models for singlet fission by machine learning

Singlet fission (SF), the conversion of one singlet exciton into two triplet excitons, could significantly enhance solar cell efficiency. Molecular crystals that undergo SF are scarce. Computational exploration may accelerate the discovery of SF materials. However, many-body perturbation theory (MBPT) calculations of the excitonic properties of molecular crystals are impractical for large-scale materials screening. We use the sure-independence-screening-and-sparsifying-operator (SISSO) machine-learning algorithm to generate computationally efficient models that can predict the MBPT thermodynamic driving force for SF for a dataset of 101 polycyclic aromatic hydrocarbons (PAH101). SISSO generates models by iteratively combining physical primary features. The best models are selected by linear regression with cross-validation. The SISSO models successfully predict the SF driving force with errors below 0.2 eV. Based on the cost, accuracy, and classification performance of SISSO models, we propose a hierarchical materials screening workflow. Three potential SF candidates are found in the PAH101 set.

36 MATERIALS SCIENCE↗