Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional 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 OSTI2022

Stochastic Vector Techniques in Ground-State Electronic Structure

Herein we review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order many-body perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss effective techniques to reduce statistical errors. Finally, we review the use of stochastic vectors for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods.

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

Electronic Structure Modification of Rectangular Phosphorene Quantum Dots Via Edge Passivation

Phosphorene is an elemental two-dimensional material that shows great promise as a metal-free photocatalyst owing to its high carrier mobility and direct band gap. Doping and edge functionalization of phosphorene can be exploited to tune its band gap and alter other properties. Herein, the electronic properties of a series of edge-passivated monolayer phosphorene quantum dots are investigated using an array of functional groups (OH, SH, NH 2 , SCH 3 , OCN, CN, and Cl) that are known to have electron-donating or electron-withdrawing capabilities. By employing density functional theory calculations, we demonstrate that edge functionalization can perturbatively shift the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) energy levels while simultaneously preserving the geometric structure and other favorable properties of pristine phosphorene. We use the average overall electrostatic potential of these systems to successfully predict the relative positioning of the HOMO and LUMO energy levels while presenting qualitative explanations on potential jumps introduced by the dipole layers formed due to edge functionalization. Furthermore, our results suggest that selective functionalization permits considerable control over the band edge positions, which could be utilized for applications in energy and optoelectronic devices.

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

Core spectroscopy of oxazole

We have measured, analyzed and simulated the ground state valence photoelectron (PES), X-ray absorption (XAS), X-ray photoelectron (XPS) as well as normal and resonant Auger-Meitner electron (AES) spectra of oxazole at the carbon, oxygen and nitrogen K-edge in order to understand its electronic structure. Experimental data is compared to theoretical calculations performed at the coupled cluster, restricted active space perturbation theory to second order and time-dependent density functional levels of theory. We demonstrate that both N and O K-edge XAS are sensitive to the amount of dynamical electron correlation included in the theoretical description, and that for a complete description of XPS, additional orbital correlation and orbital relaxation effects need to be considered. The normal AES are dominated by a singlet excitation channel and well described by theory. Furthermore, the resonant AES, however, are more complicated. While the participator decay channels, dominating at higher kinetic energies, are well described by coupled cluster theory, spectator channels can only be described satisfactorily using a method that combines restricted active space perturbation theory to second order for the bound part and a one-center approximation for the continuum.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

True molecular conformation and structure determination by three-dimensional electron diffraction of PAH by-products potentially useful for electronic applications

The true mol­ecular conformation and the crystal structure of benzo[e]di­naphtho­[2,3-a;1',2',3',4'-ghi]fluoranthene, 7,14-di­phenyl­naphtho­[1,2,3,4-cde]bis­anthene and 7,16-di­phenyl­naphtho­[1,2,3,4-cde]heli­anthrene were determined ab initio by 3D electron diffraction. All three mol­ecules are remarkable polycyclic aromatic hydro­carbons. The mol­ecular conformation of two of these com­pounds could not be determined via classical spectroscopic methods due to the large size of the mol­ecule and the occurrence of multiple and reciprocally connected aromatic rings. The mol­ecular structure of the third mol­ecule was previously considered provisional. These com­pounds were isolated as by-products in the synthesis of similar products and were at the same time nanocrystalline and available only in very limited amounts. 3D electron diffraction data, taken from submicrometric single crystals, allowed for direct ab initio structure solution and the unbiased determination of the inter­nal mol­ecular conformation. Detailed synthetic routes and spectroscopic analyses are also discussed. Based on many-body perturbation theory simulations, benzo[e]di­naphtho­[2,3-a;1',2',3',4'-ghi]fluoranthene may be a promising candidate for triplet–triplet annihilation and 7,14-di­phenyl­naphtho­[1,2,3,4-cde]bis­anthene may be a promising candidate for inter­molecular singlet fission in the solid state.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculations of anisotropic magnetic properties using spin-orbit energy variations

We analyze several methods of obtaining the accurate relativistic total energy (TE) variations using traditional perturbation theories (PTs) and proposed coupling constant integration (CCI) methods. For this purpose, we perform benchmark calculations within the density functional theory taking the spin-orbit coupling (SOC) and its derivative as a perturbation. The TE change due to SOC addition obtained from both PTs and CCI is shown to reach the accuracy of fully self-consistent TE calculations. Similar accuracy is also obtained even for the magnetocrystalline anisotropy energy (MAE). The real advantage of the proposed methods is to use PTs and CCI methods in those electronic structure methods where accurate total energies currently cannot be obtained with required accuracy. Correspondingly, we demonstrate the applicability of suggested methods for calculations of MAE in different magnetic materials using a dynamic mean-field method. All suggested PTs and CCI methods also provide convenient site, orbital, and spin decompositions of the TE variation, creating a powerful way to analyze microscopic physics in strongly correlated materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

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

Quantum Perturbation Theory Using Tensor Cores and a Deep Neural Network

In this work, time-independent quantum response calculations are performed using Tensor cores. This is achieved by mapping density matrix perturbation theory onto the computational structure of a deep neural network. The main computational cost of each deep layer is dominated by tensor contractions, i.e., dense matrix–matrix multiplications, in mixed-precision arithmetics, which achieves close to peak performance. Quantum response calculations are demonstrated and analyzed using self-consistent charge density-functional tight-binding theory as well as coupled-perturbed Hartree–Fock theory. For linear response calculations, a novel parameter-free convergence criterion is presented that is well-suited for numerically noisy low-precision floating point operations and we demonstrate a peak performance of almost 200 Tflops using the Tensor cores of two Nvidia A100 GPUs.

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

Spin-crossover complexes: Self-interaction correction vs density correction

Complexes containing a transition metal atom with a 3d 4 –3d 7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi–Löwdin self-interaction correction. Considering that Hartree–Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree–Fock method and the Perdew–Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. Here, we find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.

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

Electronic structure and optical properties of halide double perovskites from a Wannier-localized optimally-tuned screened range-separated hybrid functional

Halide double perovskites are a chemically diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and ab initio many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here, in this study, we apply the WOT-SRSH functional to five halide double perovskites and compare the results with those obtained from other known functionals and previous GW calculations. We also use the approach as a starting point for GW calculations and we compute the band structures and optical absorption spectrum for Cs 2 AgBiBr 6 , using both time-dependent DFT and the GW-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exact two-body expansion of the many-particle wave function

Progress toward the solution of the strongly correlated electron problem has been stymied by the exponential complexity of the wave function. Previous work established an exact two-body exponential product expansion for the ground-state wave function. By developing a reduced density-matrix analog of Dalgarno-Lewis perturbation theory, we prove here that (i) the two-body exponential product expansion is rapidly and globally convergent with each operator representing an order of a renormalized perturbation theory, (ii) the energy of the expansion converges quadratically near the solution, and (iii) the expansion is exact for both ground and excited states. The two-body expansion offers a reduced parametrization of the many-particle wave function as well as the two-particle reduced density matrix with potential applications on both conventional and quantum computers for the study of strongly correlated quantum systems. In this work, we demonstrate the result with the exact solution of the contracted Schrödinger equation for the molecular chains H 4 and H 5 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Revealing symmetry-broken superconducting configurations by density functional theory

A coherent theory for the superconductivity of both conventional and unconventional superconductors is currently lacking. Here we show that superconductivity arises from the formation of a symmetry-broken superconducting configuration (SCC) due to atomic perturbation of the normal conducting configuration (NCC). This electron–phonon interaction creates straight one-dimensional tunnels (SODTs) for charge density of electrons and/or holes as revealed by the calculations based on density functional theory (DFT). The SODTs act as resistance-free superhighways and are correlated to the Cooper pairs in the Bardeen–Cooper–Schrieffer (BCS) theory. The formation of SODTs implies that the electron–phonon interaction in the BCS theory can be represented by the difference in charge densities between SCC and NCC predicted by DFT. The present work highlights that in conventional superconductors, SODTs are embedded within the bulk materials and are easily destroyed by phonon vibrations, resulting in a low critical superconducting temperature (T C ). Conversely, in unconventional superconductors such as YBa 2 Cu 3 O 7 (YBCO 7 ), SODTs are protected by a layered pontoon structure with very weak bonding to the bulk materials, maintaining SODTs’ stability at higher temperatures and leading to a much higher T C . The present approach is validated for 14 conventional superconductors of 18 pure elements and MgB 2 examined in this work, including the presently predicted superconductivity in Cu, Ag, Au, Sb, Bi, and MgB 2 at 0 K and 0 GPa, and one unconventional superconductor of YBCO 7 . Our discovery indicates that DFT can be a practical tool for predicting superconductors, enabling a systematic search for new superconducting materials in the future.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Kinetic effects of Alfven wave nonlinearity. I - Ponderomotive density fluctuations

The Vlasov theory is used to study kinetic corrections to fluid descriptions of Alfven wave nonlinearity. The method is to obtain an expression for the second-order perturbed distribution function produced by a nonlinear Alfven wave. From this distribution function a kinetically correct expression is obtained for the plasma density perturbation associated with an envelope-modulated Alfven wave. This kinetic theory result differs substantially from the fluid expression when the plasma beta is greater than about 1, and the electron and ion temperatures are approximately equal. This result is of interest because density fluctuations are an observationally accessible indicator of wave nonlinearity in solar system Alfven waves. It also will assist in the determination of properties of Alfven waves in the interstellar medium. Finally, this analysis also yields a kinetically correct expression for u, the magnetic field-aligned component of the plasma fluid velocity.

Spangler, Steven R.↗
NASA NTRS1998

C2 Fragmentation Energy of C60 Revisited: Theory Disagrees with Most Experiments

Following our earlier work on the subject, we have carried out density functional theory (DFT) and second-order Moller-Plesset perturbation theory (MP2) calculations of the dissociation energy of the reaction C60 yields C58 + C2 using polarized basis sets and geometries optimized with DFT methods. The present theoretical results support an electronic fragmentation energy D(sub e) around 10-11 eV in disagreement with most experimental results that place the dissociation energy D(sub o) (including zero point energy) around 7-8 eV. The plausible errors remaining in the theoretical calculations are unlikely to account for this big difference (2-4 eV).

Boese, A. Daniel↗
NASA NTRS1975

High zonal harmonics of rapidly rotating planets

A new perturbation expansion is derived for the structure of rotating bodies in hydrostatic equilibrium. The method uses an expansion of the density on Legendre polynomial functions of angle, and can be developed analytically in a manner analogous to the standard level-surface perturbation theory. The new theory proceeds from a prescribed pressure-density relation rather than from a prescribed density distribution, and is both simpler and more physically transparent than the level-surface approach. High zonal harmonics are shown to arise via a transfer function involving derivatives of the interior sound velocity, and via mixing of multipole density components in the outer shell of the planet. Sample calculations for polytropic sequences are presented, as well as standard gravity models for Jupiter and Saturn. Mathematical subleties of the theory are discussed in an appendix.

Hubbard, W. B.↗
DOE OSTI2023

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

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

Near-Quantitative Predictions of the First-Shell Coordination Structure of Hydrated First-Row Transition Metal Ions Using K-Edge X-ray Absorption Near-Edge Spectroscopy

Understanding the solvation structure of transition metal ions is important due to its broad impact on the kinetics, stability, and reactivity in a variety of applications in geochemistry, biochemistry, energy storage, and environmental chemistry. Using water as the common ligand, we study the X-ray absorption pre-edge and near-edge spectra at the K-edge of a near-complete series of hydrated first-row transition metal ions with d-orbital occupancy ranging from d 2 to d 10 . Starting with optimized structures that were derived from an explicit solvation treatment at the density functional theory (DFT) level, we then compute the pre-edge X-ray absorption spectra at the metal ion K-edge at the time dependent density functional theory (TDDFT) and restricted active-space second-order perturbation theory (RASPT2) levels of theory. TDDFT calculations provide accurate results for spectra that are dominated by single excitations, while significant improvements were obtained with RASPT2 calculations that correctly distinguish between singly and doubly excited states, where relevant, with quantitative accuracy compared with experiment. We analyze and assign the dierent pre-edge features for each metal ion in order to reveal the impact of the variations in the d orbital occupancy on the rst-shell coordination environment. For completeness, we also report the lowest energy ligand-field d-d transitions for all the transition metal aqua ions considered in this study using complete active space second order perturbation theory (CASPT2).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗