Search NASA⌕ Search

SEARCH · Search NASA

Results for “MBPT”

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.

18 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Accurate CI-MBPT calculation of radiative lifetimes and transition probabilities of neutral lanthanum (La I) odd states with J = 3/2

In this work, an accurate calculation of radiative lifetimes and transition probabilities of the first 21 odd states with a total angular momentum J = 3/2 was performed for La I using configuration-interaction many-body perturbation theory (CI-MBPT). A technique to obtain correct transition amplitudes was devised for pairs of states with substantial mixing that improved the theoretical accuracy. The results were compared to known experimental and theoretical data for the considered transitions. It has been shown that the transition data calculated with CI-MBPT and with corrections for pairs of mixed states are in good agreement with the latest experimental data of Den Hartog et al, 2015 J. Phys. B: At. Mol. Phys. 48, 155 001, better than other theoretical results. The transition data can find many applications in astrophysics, and the CI-MBPT formalism with mixing correction is applicable to other atoms similar to or more complex than La I.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CI-MBPT and Intensity-Based Lifetime Calculations for Th II

Lifetime calculations of Th II J = 1.5 and 2.5 odd states are performed with configuration–interaction many-body perturbation theory (CI-MBPT). For many J = 2.5 states, lifetimes are quite accurate, but two pairs of J = 2.5 odd states and many groups of J = 1.5 states are strongly mixed, making theoretical predictions unreliable. To solve this problem, a method based on intensities is used. To relate experimental intensities to lifetimes, two parameters, one an overall coefficient of proportionality for transition rates and one temperature of the Boltzmann distribution of populations, are introduced and fitted to minimize the deviation between theoretical and intensity-derived lifetimes. For strongly mixed groups of states, the averaged lifetimes obtained from averaged transition rates were used instead of individual lifetimes in the fit. Close agreement is obtained. Then intensity branching ratios are used to extract individual lifetimes for the strongly mixed states. The resulting lifetimes are compared to available directly measured lifetimes and reasonable agreement is found, considering limited accuracy of intensity measurements. The method of intensity-based lifetime calculations with fit to theoretical lifetimes is quite general and can be applied to many complex atoms where strong mixing between multiple states exists.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

CI-MBPT line strengths and atomic probabilities for some transitions of neutral iodine

Iodine is a promising alternative to expensive xenon as a propellant for electric propulsion devices. Iodine optical spectroscopy is a valuable diagnostic tool for optimization of the iodine propellant. Recently, some measurements have been conducted with cells containing iodine plasma, and iodine atomic absorption at several wavelengths (911, 906, 206 and 1315 nm) was measured with a Ti:sapphire laser. Motivated by these experiments and future applications, here we have calculated line strengths and atomic probabilities for the transitions between different types of 5s 2 5p 5 , 5s 2 5p 4 6s, 5s 2 5p 4 6p and 5s 2 5p 4 5d states of atomic iodine, using configuration-interaction many-body perturbation theory. Since experimental transition data for atomic iodine are only available for a relatively small number of transitions, we also performed calculations for the bromine atom, to further validate our approach. We found good agreement between our theory and available experimental data for both elements in most cases. Therefore, we believe that our predicted line strengths and transition probabilities are quite reliable and will be useful for iodine spectroscopy.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A new “gold standard”: Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing

A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the “gold standard” CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.

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

Relativistic Configuration-Interaction and Perturbation Theory Calculations of the Sn XV Emission Spectrum

Recently, there has been increased interest in developing advanced bright sources for lithography. Sn ions are particularly promising due to their bright emission spectrum in the required wavelength range. Cowan’s code has been used to model the emission; however, it has adjustable parameters, which limit its predictive power, and it has limited relativistic treatment. Here, we present calculations based on ab initio relativistic configuration-interaction many-body perturbation theory (CI-MBPT), with relativistic corrections included at the Dirac-Fock level and core-polarization effects with the second-order MBPT. As a proof of principle that the theory is generally applicable to other Sn ions with proper development, we focused on one ion where direct comparison with experimental observations is possible. The theory can also be used for ions of other elements to predict emissions for optimization of plasma-based bright sources.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Relativistic Configuration-Interaction and Perturbation Theory Calculations for Heavy Atoms

Heavy atoms present challenges to atomic theory calculations due to the large number of electrons and their complicated interactions. Conventional approaches such as calculations based on Cowan’s code are limited and require a large number of parameters for energy agreement. One promising approach is relativistic configuration-interaction and many-body perturbation theory (CI-MBPT) methods. We present CI-MBPT results for various atomic systems where this approach can lead to reasonable agreement: La I, La II, Th I, Th II, U I, Pu II. Among atomic properties, energies, g-factors, electric dipole moments, lifetimes, hyperfine structure constants, and isotopic shifts are discussed. While in La I and La II accuracy for transitions is better than that obtained with other methods, more work is needed for actinides.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory

Coupled cluster methods based exclusively on double excitations are comparatively “cheap” and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T 2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Here, our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.

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

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange

We report on the use of the adaptively compressed exchange (ACE) operator to accelerate many-body perturbation theory (MBPT) calculations, including G 0 W 0 and the Bethe–Salpeter equation (BSE), for hybrid density functional theory starting points. We show that by approximating the exact exchange operator with the low-rank ACE operator, substantial computational savings can be achieved with systematically controllable errors in the quasiparticle energies computed with full-frequency G 0 W 0 and the optical absorption spectra and vertical excitation energies computed by solving the BSE within density matrix perturbation theory. Our implementation makes use of the ACE-accelerated electronic Hamiltonian to carry out both G 0 W 0 and BSE without explicitly computing empty states. We show the robustness of the approach and present the computational gains obtained on both the central processing unit and graphics processing unit nodes. In conclusion, our work will facilitate the exploration and evaluation of fine-tuned hybrid starting points aimed at enhancing the accuracy of MBPT calculations without involving computationally demanding self-consistency in Hedin’s equations.

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

Ligand‐Mediated Quantum Yield Enhancement in 1‐D Silver Organothiolate Metal–Organic Chalcogenolates

X-ray free electron laser (XFEL) microcrystallography and synchrotron single-crystal crystallography are used to evaluate the role of organic substituent position on the optoelectronic properties of metal–organic chalcogenolates (MOChas). MOChas are crystalline 1D and 2D semiconducting hybrid materials that have varying optoelectronic properties depending on composition, topology, and structure. While MOChas have attracted much interest, small crystal sizes impede routine crystal structure determination. A series of constitutional isomers where the aryl thiol is functionalized by either methoxy or methyl ester are solved by small molecule serial femtosecond X-ray crystallography (smSFX) and single crystal rotational crystallography. While all the methoxy examples have a low quantum yield (0-1%), the methyl ester in the ortho position yields a high quantum yield of 22%. Here, the proximity of the oxygen atoms to the silver inorganic core correlates to a considerable enhancement of quantum yield. Four crystal structures are solved at a resolution range of 0.8–1.0 Å revealing a collapse of the 2D topology for functional groups in the 2- and 3- positions, resulting in needle-like crystals. Further analysis using density functional theory (DFT) and many-body perturbation theory (MBPT) enables the exploration of complex excitonic phenomena within easily prepared material systems.

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

Analysis of Fourth-, Fifth-, and Infinite-Order Triple Excitations in Unitary Coupled Cluster Theory

Here, in this work, we introduce a perturbative correction to the unitary coupled cluster method with single and double excitations (UCCSD) that incorporates the effects of missing triple excitations through fifth-order in many-body perturbation theory (MBPT). Referred to as UCCSD[T-5], this method is benchmarked alongside the previously developed UCCSD[T] to lend insight into the behavior of perturbative triples corrections relative to UCCSDT, which inherently provides an infinite-order treatment of triple excitations, as well as full configuration interaction (FCI). Two key findings emerge from this analysis. First, UCCSD[T] consistently yields ground-state energies in closest agreement with FCI, outperforming both UCCSD[T-5] and UCCSDT. Second, UCCSD[T-5] largely emulates the behavior of UCCSDT, suggesting a close relationship between fifth- and infinite-order triple-excitation contributions. Given the growing interest in UCC ansätze for quantum computing and the limitations of current quantum hardware, these results highlight the potential of classically computed perturbative corrections within UCC theory to capture triple-excitation effects without the additional quantum resources required by the UCCSDT ansatz.

Windom, Zachary W. [Oak Ridge National Laboratory ↗
DOE OSTI2023

First-Principles Understanding of Optical Excitations within Low-Dimensional Materials

The objective of this project is to utilize first-principles computational approaches to understand optical excitations within two-dimensional monolayers and heterostructures, with the ultimate goal of designing new materials by modifying the chemical and physical structure on the nanoscale. In 2D heterostructures, the nature and movement of the exciton is determined to a large extent by the electron-electron and electron-phonon strengths within individual monolayers, as well as the coupling between layers, all of which in turn depend on chemistry, solid-state screening properties, and quantum confinement. Controlling the nature and migration properties of excitons requires an understanding of the complex relationship between electrons, phonons, defects, disorder, and dynamics. The aim of this project is to determine, by analysis of highly accurate density functional theory (DFT) and many-body perturbation theory (MBPT) calculations, how these complexities can be decomposed into simple, tunable parameters. The specific goals of this project are to 1) understand electron conductivity, absorption and transparency of monolayer semiconductors and metals; 2) develop a theory of excitons near defects in low dimensions; and 3) understand the influence of inter-layer interactions in heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

Theoretical vibrations of carbon chains C3, C4, C5, C6, C7, C8, and C9

The MBPT (2) procedure with the 6-31g (asterisk) basis set was used to study nearly linear carbon chains. The theoretical vibrational frequencies of the molecules C3 through C9 are presented and, for C3 through C6, compared to experimental stretching frequencies and their (C-13)/(C-12) isotopomers. Predictions for C7, C8, and C9 stretching frequencies are calculated by directly scaling the theoretical frequencies with factors derived from experimental-to-theoretical ratios known for the smaller molecules.

Kurtz, Joe↗
NASA NTRS1993

An open-shell restricted Hartree-Fock perturbation theory based on symmetric spin orbitals

A new open-shell perturbation theory is formulated in terms of symmetric spin orbitals. Only one set of spatial orbitals is required, thereby reducing the number of independent coefficients in the perturbed wavefunctions. For second order, the computational cost is shown to be similar to a closed-shell calculation. This formalism is therefore more efficient than the recently developed RMP, ROMP or RMP-MBPT theories. The perturbation theory described herein was designed to have a close correspondence with our recently proposed coupled-cluster theory based on symmetric spin orbitals. The first-order wavefunction contains contributions from only doubly excited determinants. Equilibrium structures and vibrational frequencies determined from second-order perturbation theory are presented for OH, NH, CH, 02, NH2 and CH2.

Lee, Timothy J.↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Response Functions of Correlated Systems within Green's Function Theory

Why use Green's functions as the fundamental variable? Wave-function (Psi) methods are king for high-fidelity and Density-functional (Rho) methods are very efficient (Kohn-Sham). Goldilocks principle: Green's function (G) methods straddle the Rho and Psi methods, intermediate in both accuracy and efficiency. Also, when interest lies in excitations & 2-particle properties: G-methods are natural - intrinsic to the theory.

DMFT↗