Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density function 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

X-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula

Here, we study ab initio approaches for calculating x-ray Thomson scattering spectra from density functional theory molecular dynamics simulations based on a modified Chihara formula that expresses the inelastic contribution in terms of the dielectric function. We study the electronic dynamic structure factor computed from the Mermin dielectric function using an ab initio electron-ion collision frequency in comparison to computations using a linear-response time-dependent density functional theory (LR-TDDFT) framework for hydrogen and beryllium and investigate the dispersion of free-free and bound-free contributions to the scattering signal. A separate treatment of these contributions, where only the free-free part follows the Mermin dispersion, shows good agreement with LR-TDDFT results for ambient-density beryllium, but breaks down for highly compressed matter where the bound states become pressure ionized. LR-TDDFT is used to reanalyze x-ray Thomson scattering experiments on beryllium demonstrating strong deviations from the plasma conditions inferred with traditional analytic models at small scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. Here, this study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

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

Symmetry-correct bonding in density functional theory calculations for delta phase Pu

The long-held conclusion that magnetic order in delta phase Pu makes the structure mechanically unstable in density functional theory calculations is incorrect if the magnetic order is allowed to be non-collinear. The non-collinear 3Q spin structure is intrinsically cubic and makes all structurally equivalent bonds equivalent in their bonding character. Applied to density functional theory calculations on delta phase Pu, the 3Q spin structure results in elastic constants and phonons with the correct symmetry. Finally, the calculated phonon dispersion agrees better with experiment than previous calculations, and the calculated thermal expansion shows the unique behavior seen experimentally: it is negative.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Density-functional theory of material design: fundamentals and applications-I

This article is part-I of a review of density-functional theory (DFT) that is the most widely used method for calculating electronic structure of materials. The accuracy and ease of numerical implementation of DFT methods has resulted in its extensive use for materials design and discovery and has thus ushered in the new field of computational material science. In this article, we start with an introduction to Schrödinger equation and methods of its solutions. After presenting exact results for some well-known systems, difficulties encountered in solving the equation for interacting electrons are described. How these difficulties are handled using the variational principle for the energy to obtain approximate solutions of the Schrödinger equation is discussed. The resulting Hartree and Hartree–Fock theories are presented along with results they give for atomic and solid-state systems. We then describe Thomas–Fermi theory and its extensions which were the initial attempts to formulate many-electron problem in terms of electronic density of a system. Having described these theories, we introduce modern DFT by discussing Hohenberg–Kohn theorems that form its foundations. We then go on to discuss Kohn–Sham (KS) formulation of DFT in its exact form. Next, local density approximation (LDA) is introduced and solutions of KS equation for some representative systems, obtained using the LDA, are presented. We end part-I of the review describing the contents of part-II.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Charge self-consistent density functional theory plus ghost rotationally invariant slave-boson theory for correlated materials

We present a charge self-consistent density functional theory combined with the ghost rotationally invariant slave-boson (DFT+gRISB) formalism for studying correlated materials. Here, this method is applied to SrVO 3 and NiO, representing prototypical correlated metals and charge-transfer insulators. For SrVO 3 , we demonstrate that DFT+gRISB yields an accurate equilibrium volume and effective mass close to experimentally observed values. Regarding NiO, DFT+gRISB enables the simultaneous description of charge-transfer and Mott-Hubbard bands, significantly enhancing the accuracy of the original DFT+RISB approach. Furthermore, the calculated equilibrium volume and spectral function reasonably agree with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Improving the precision of forces in real-space pseudopotential density functional theory

The high-order finite difference real-space pseudopotential density functional theory (DFT) approach is a valuable method for large-scale, massively parallel DFT calculations. A significant challenge in the approach is the oscillating “egg-box” error introduced by aliasing associated with a coarse grid spacing. To address this issue while minimizing computational cost, we developed a finite difference interpolation (FDI) scheme [Roller et al., J. Chem. Theory Comput. 19, 3889 (2023)] as a means of exploiting the high resolution of the pseudopotential to reduce egg-box effects systematically. Here, we show an implementation of this method in the PARSEC code and examine the practical utility of the combination of FDI with additional methods for improving force precision and/or reducing its computational cost, including orbital-based forces, compensating charges (namely, adding and subtracting a judiciously chosen charge density such that the total density is unaltered), and a modified spatial domain in which the real-space grid is defined. Using selected small molecules, as well as metallic Li, as test cases, we show that a combination of all four aspects leads to a significant reduction in computational cost while retaining a high level of precision that supports accurate structures and vibrational spectra, as well as stable and accurate molecular dynamics runs.

Chemistry↗
DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

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

Nonadiabatic Molecular Dynamics by Multiconfiguration Pair-Density Functional Theory

Herein, we present the first implementation of multiconfiguration pair-density functional theory (MC-PDFT) ab initio molecular dynamics. MC-PDFT is a multireference electronic structure method that in many cases has a similar accuracy (or even better accuracy) the complete active space second-order perturbation theory (CASPT2) at a significantly lower computational cost. In this study, we introduced MC-PDFT analytical gradients into the SHARC molecular dynamics program for ab initio, nonadiabatic molecular dynamics simulations. We verify our implementation by examining the intersystem crossing dynamics of thioformaldehyde, and we observe excellent agreement with recent CASPT2 and experimental findings. Moreover, with MC-PDFT, we could perform dynamics simulations with the 12 electron in 10 orbitals active space that was computationally too expensive for direct dynamics with CASPT2.

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

Real-space Kohn–Sham density functional theory for complex energy applications

Real-space Kohn-Sham density functional theory (real-space KS-DFT) enables large-scale electronic structure simulations that is particularly well-suited for the modern high-performance computing (HPC) architectures. This feature article reviews its theoretical foundations, highlights the algorithmic advances and recent developments, and showcases applications in complex nano systems. We aim to provide a perspective on the trajectory of real-space KS-DFT as an emerging tool for computational chemistry and materials science in the exascale era.

Zhang, Zeyi↗
DOE OSTI2022

Density functional theory investigations into the magnetic ordering of U 3 O 8

Density functional theory (DFT) has been highly successful in supporting experimental materials science; however, a correct electronic ground state is required to realize the full theoretical capacity of DFT. The uranium oxides, α–U 3 O 8 in particular, are simultaneously technologically important materials and theoretically challenging for DFT because the uranium magnetic ground state is not obvious. This is true for both experiment and theory—magnetic susceptibility measurements indicate an antiferromagnetic (AFM) ground state with transitions near 4.2 and 8.0 K, but the ordering itself is not known. Theoretical literature reports are in contradiction, with independent studies finding paramagnetic, ferromagnetic (FM), and AFM states as the lowest energy configuration. However, recent inelastic neutron scattering experiments suggested an uninvestigated magnetic structure with ordering along the [0.5 1 1] plane, motivating a theoretical reinvestigation. Using this insight, we calculated the relative energy of FM and AFM orderings along [0.5 1 1], [0.5 0 0], [0 1 0], and [0 0 1] using noncollinear DFT calculations with spin-orbital coupling. In this work, we found that the [0.5 1 1] AFM structure is lower in energy than FM or AFM orderings along the low Miller index directions. We also investigated polarization of the magnetic moment along each lattice vector and found that polarization along the out-of-plane direction is the energetically preferred orientation for the AFM structures. Additionally, we found in all calculations that moments initially pointing along the in-plane lattice vectors significantly relax until they point along the coordinate between the two crystallographically distinct uranium sites with complex noncollinear magnetic configurations. The new [0.5 1 1] AFM magnetic structure provides an additional path forward toward understanding the electronic structure of α-U 3 O 8 and lends theoretical credibility to recent neutron scattering results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗
DOE OSTI2021

Global study of separable pairing interaction in covariant density functional theory

A systematic global investigation of pairing properties based on all available experimental data on pairing indicators has been performed for the first time in the framework of covariant density functional theory. It is based on the separable pairing interaction. The optimization of the scaling factors of this interaction to experimental data clearly reveals its isospin dependence in the neutron subsystem. However, the situation is less certain in the proton subsystem since similar accuracy of the description of pairing indicators can be achieved both with isospin-dependent and mass-dependent scaling factors. The differences in the functional dependencies of scaling factors lead to uncertainties in the prediction of proton and neutron pairing properties which are especially pronounced at high isospin and could have a significant impact on some physical observables. For a given part of the nuclear chart the scaling factors for spherical nuclei are smaller than those for deformed ones; this feature exists also in nonrelativistic density functional theories. Furthermore, its origin is traced back to particle-vibration coupling in odd-A nuclei which is missing in all existing global studies of pairing. Although the present investigation is based on the NL5(E) covariant energy density functional (CEDF), its general conclusions are expected to be valid also for other CEDFs built at the Hartree level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quantum Monte Carlo and Density Functional Theory Study of Strain and Magnetism in 2D 1T-VSe 2 with Charge Density Wave States

Two-dimensional (2D) 1T-VSe 2 has prompted significant interest due to the discrepancies regarding alleged ferromagnetism (FM) at room temperature, charge density wave (CDW) states, and the interplay between the two. We employed a combined Diffusion Monte Carlo (DMC) and density functional theory (DFT) approach to accurately investigate the magnetic properties, CDW states, and their responses to strain in monolayer 1T-VSe 2 . Our calculations show the delicate competition between various phases, revealing critical insights into the relationship between their energetic and structural properties. Here, we performed classical Monte Carlo simulations informed by our DMC and DFT results and found the magnetic transition temperature (T c ) of the undistorted (non-CDW) FM phase to be 228 K and the distorted (CDW) phase to be 68 K. Additionally, we studied the response of biaxial strain on the energetic stability and magnetic properties of various phases of 2D 1T-VSe 2 and found that small amounts of strain can increase the T c , suggesting a promising route for engineering and enhancing magnetic behavior. Finally, we synthesized 1T-VSe 2 and performed Raman spectroscopy measurements, which were in close agreement with our calculated results, validating our computational approach. Our work emphasizes the role of highly accurate DMC methods in advancing the understanding of monolayer 1T-VSe 2 and provides a robust framework for future studies of 2D magnetic materials.

2D magnets↗
DOE OSTI2020

Tuned and screened range-separated hybrid density functional theory for describing electronic and optical properties of defective gallium nitride

We apply a hybrid density functional theory approach, based on a tuned and screened range-separated hybrid (SRSH) exchange-correlation functional, to describe the optoelectronic properties of defective gallium nitride (GaN). SRSH and time-dependent SRSH (TDSRSH) are tuned to produce accurate energetics for the pristine material and applied to the study of a series of point defects in bulk GaN, a blue light-emitting material that degrades in the presence of defects. We first establish the accuracy of the method by comparing the predicted quasiparticle gap and low-energy excitation spectra of (TD)SRSH and many-body perturbation theory for both pristine GaN and GaN containing a single nitrogen vacancy. Aided by the reduced computational cost of (TD)SRSH, we then report on three additional technologically relevant point defects and defect complexes in GaN: the gallium vacancy, the carbon interstitial, and the carbon-silicon complex. We compute the low-energy optical absorption spectra for these defects and show the presence of defect-centered transitions. Furthermore, by estimating the Stokes shift, we predict, in agreement with previous studies, that the carbon substitutional defect is a candidate for the detrimental yellow luminescence in GaN. Furthermore, this study indicates that TDSRSH is a promising and computationally feasible approach for quantitatively accurate, first-principles modeling of defective semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantum Chemical Simulations of CO 2 and N 2 Capture in Reline, a Prototypical Deep Eutectic Solvent

Deep eutectic solvents such as reline are an emerging class of low-cost, environmentally friendly solvents with tunable properties that are potentially applicable for the capture and separation of CO 2 . Experimental measurements showed that a reline-based membrane contactor can capture and separate CO 2 via physisorption through a dissolution process with 96.7% purity from a mixed gas containing CO 2 and N 2 (50:50% molar ratio). Here, we examine the nature of the interaction of CO 2 and N 2 with reline employing quantum chemical methods. We focus on explaining the mechanism by which CO 2 and N 2 bind to reline and the reason for the high selectivity for absorption of CO 2 compared to N 2 . We analyze the dynamics, energetics, and binding motifs for CO 2 and N 2 in reline employing density functional theory, density functional tight binding, and ab initio molecular dynamics. We also investigate the effect of reline on the vibrational spectra of CO 2 and reline. Our simulations indicate that the selective capture of CO 2 from the mixture of CO 2 and N 2 is due to the interplay between attractive electrostatic and charge polarization forces with opposing entropic effects, which shift the energetic balance and make the N 2 absorption unfavorable in reline.

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

Accelerating Finite-Temperature Kohn-Sham Density Functional Theory with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

97 MATHEMATICS AND COMPUTING↗