Search NASA⌕ Search

SEARCH · Search NASA

Results for “perturbation theory”

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

Simplified tuning of long-range corrected density functionals for use in symmetry-adapted perturbation theory

Long considered a failure, second-order symmetry-adapted perturbation theory (SAPT) based on Kohn–Sham orbitals, or SAPT0(KS), can be resurrected for semiquantitative purposes using long-range corrected density functionals whose asymptotic behavior is adjusted separately for each monomer. Here, as in other contexts, correct asymptotic behavior can be enforced via “optimal tuning” based on the ionization energy theorem of density functional theory, but the tuning procedure is tedious, expensive for large systems, and comes with a troubling dependence on system size. Here, we show that essentially identical results are obtained using a fast, convenient, and automated tuning procedure based on the size of the exchange hole. In conjunction with “extended” (X)SAPT methods that improve the description of dispersion, this procedure achieves benchmark-quality interaction energies, along with the usual SAPT energy decomposition, without the hassle of system-specific tuning.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Here we present a workflow of practical calculations of electron-phonon (e-ph) coupling with many-electron correlation effects included using the GW perturbation theory (GWPT). This workflow combines BerkeleyGW, ABINIT, and EPW software packages to enable accurate e-ph calculations at the GW self-energy level, going beyond standard calculations based on density functional theory (DFT) and density-functional perturbation theory (DFPT). This workflow begins with DFT and DFPT calculations (ABINIT) as starting point, followed by GW and GWPT calculations (BerkeleyGW) for the quasiparticle band structures and e-ph matrix elements on coarse electron k- and phonon q-grids, which are then interpolated to finer grids through Wannier interpolation (EPW) for computations of various e-ph coupling determined physical quantities such as the electron self-energies or solutions of anisotropic Eliashberg equations, among others. A gauge-recovering symmetry unfolding technique is developed to reduce the computational cost of GWPT (as well as DFPT) while fulfilling the gauge consistency requirement for Wannier interpolation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Statistical uncertainty of fission matrix eigenvalues using perturbation theory

Eigenvalue search of high dominance ratio systems may be slow to converge. The fission matrix element is defined by its element (FM){sub ij}, which are the probability for a neutron born in cell i to create a fission in cell j for a spatial mesh of n{sub i}*n{sub j} elements. Fission matrices are used in Monte Carlo criticality simulations to enhance computing speed, but also to find higher order eigenvalues. However, few studies have been made on the link between statistical uncertainties of fission matrix elements and eigenvalues uncertainties. Thus, dominance ratio statistical uncertainties remain unknown. This paper uses a new generalized perturbation theory (GPT) method to estimate sensitivities of eigenvalues to fission matrix elements and then to calculate dominance ratio uncertainties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Studying excited-state-specific perturbation theory on the Thiel set

Here, we explore the performance of a recently introduced N 5 -scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set’s inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles norm offers a relatively low-cost way to test for doubly excited character without the need to define an active space.

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

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Despite decades of practice, finite-size errors in many widely used electronic structure theories for periodic systems remain poorly understood. For periodic systems using a general Monkhorst-Pack grid, there has been no comprehensive and rigorous analysis of the finite-size error in the Hartree-Fock theory (HF) and the second order Møller-Plesset perturbation theory (MP2), which are the simplest wavefunction based method, and the simplest post-Hartree-Fock method, respectively. Such calculations can be viewed as a multi-dimensional integral discretized with certain trapezoidal rules. Due to the Coulomb singularity, the integrand has many points of discontinuity in general, and standard error analysis based on the Euler-Maclaurin formula gives overly pessimistic results. The lack of analytic understanding of finite-size errors also impedes the development of effective finite-size correction schemes. We propose a unified analysis to obtain sharp convergence rates of finite-size errors for the periodic HF and MP2 theories. Our main technical advancement is a generalization of the result of Lyness [Math. Comp. 30 (1976), pp. 1–23] for obtaining sharp convergence rates of the trapezoidal rule for a class of non-smooth integrands. Our result is applicable to three-dimensional bulk systems as well as low dimensional systems (such as nanowires and 2D materials). Our unified analysis also allows us to prove the effectiveness of the Madelung-constant correction to the Fock exchange energy, and the effectiveness of a recently proposed staggered mesh method for periodic MP2 calculations (see X. Xing, X. Li, and L. Lin [J. Chem. Theory Comput. 17 (2021), pp. 4733–4745]). In conclusion, our analysis connects the effectiveness of the staggered mesh method with integrands with removable singularities, and suggests a new staggered mesh method for reducing finite-size errors of periodic HF calculations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Perturbation theory in thermosphere dynamics

It is shown that density and pressure throughout the thermosphere can be adequately described in a logarithmic expansion that provides a sound basis for the application of perturbation theory. This expansion eliminates most of the important nonlinearities associated with density variations. On the basis of this expansion, the validity of perturbation theory can be extended to cover a large variety of atmospheric conditions in which the relative temperature amplitude is less than 0.5 and wind velocities are significantly less than the speed of sound.

Mayr, H. G.↗
NASA NTRS1977

Diagrammatic perturbation theory - N2 X1 Sigma/plus/g

The diagrammatic many-body perturbation theory is used to calculate the correlation energy of the nitrogen molecule in its electronic ground state. Using the algebraic approximation, the energy is evaluated through third order, including all many-body effects. (2/1) Pade approximants and variational upper bounds are constructed. For one of the perturbation expansions considered, the (2/1) Pade approximant leads to the recovery of 79.5 percent of the empirical correlation energy, while the variational upper bound recovers 72.0 percent. Three-body effects are examined in some detail. The relationships with previous work on N2 are discussed.

Wilson, S.↗
DOE OSTI2023

Anharmonic phonon behavior via irreducible derivatives: Self-consistent perturbation theory and molecular dynamics

Cubic phonon interactions are now regularly computed from first principles, and the quartic interactions have begun to receive more attention. Given this realistic anharmonic vibrational Hamiltonian, the classical phonon Green's function can be precisely measured using molecular dynamics, which can then be used to rigorously assess the range of validity for self-consistent diagrammatic approaches in the classical limit. Here we use the bundled irreducible derivative approach to efficiently and precisely compute the cubic and quartic phonon interactions of CaF 2 , systematically obtaining the vibrational Hamiltonian purely in terms of irreducible derivatives. We demonstrate that the 4-phonon sunset diagram has an important contribution to the optical phonon linewidths beyond T = 500 K. Reasonable results are obtained even at T = 900 K when performing self-consistency using the 4-phonon loop diagram and evaluating the 3-phonon bubble and 4-phonon sunset diagrams post-self-consistency. Further improvements are obtained by performing quasiparticle perturbation theory, where both the 4-phonon loop and the real part of the 3-phonon bubble are employed during self-consistency. Finally, our irreducible derivative approach to self-consistent perturbation theory is a robust tool for studying anharmonic phonons in both the quantum and classical regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗