Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

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

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

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

Analytic gradients for compressed multistate pair-density functional theory

Photochemical reactions often involve states that are closely coupled due to near degeneracies, for example by proximity to conical intersections. Therefore, a multistate method is used to accurately describe these states; for example, ordinary perturbation theory is replaced by quasidegenerate perturbation theory. Multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to approximate the full dynamical correlation energy of strongly correlated systems, and we recently proposed compressed multistate pair-density functional theory (CMS-PDFT) to treat closely coupled states. In the present paper, we report the implementation of analytic gradients for CMS-PDFT in both OpenMolcas and PySCF, and here we illustrate the use of these gradients by applying the method to the excited states of formaldehyde and phenol.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/$\mathrm{cm^3}$ to 16 g/$\mathrm{cm^3}$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. Here, we find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Multiferroic ground states in free standing perovskite-based nanodots: a density functional theory study

In this work, we use density functional theory to investigate the possibility of polar and multiferroic states in free-standing, perovskite-based nanodots at the atomic limit of miniaturization: single unit cells with terminations which allow centro-symmetry. We consider both A-O and B-O 2 terminations for three families of nanodots: (i) A = Ba with B = Ti, Zr, and Hf; (ii) A = Ca and Sr with B = Ti; and (iii) A = Na and K with B = Nb. We find all A–O terminated dots to be non-polar and to exhibit cubic symmetry (except for K 8 NbO 6 ), regardless of the presence of ferroelectricity in the bulk. In contrast, all the B–O 2 terminated nanodots considered relax to a non-cubic ground state. Rather surprisingly, all of these structures exhibit polar ground states (except NaNb 8 O 12 ). We propose a new structural parameter, the cluster tolerance factor (CTF), to determine whether a particular chemistry will result in a polar ground state nanodot, analogous to the Goldschmidt factor for bulk ferroelectrics. In addition, we find that all A–O terminated (except Ca 8 TiO 6 ) and all polar B–O 2 terminated nanodots are magnetic, where none show magnetism in the bulk. As with bulk systems, multiferroicity in the B–O 2 terminated dots originates from separation between spin density in peripheral B atoms and polarity primarily caused by the off-center central A atom. Our findings stress that surface termination plays a crucial role in determining whether ferroelectricity is completely suppressed in perovskite-based materials at their limit of miniaturization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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

Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package

Constrained density functional theory (CDFT) is a powerful tool for the prediction of electron transfer parameters in condensed phase simulations at a reasonable computational cost. In this work we present an extension to CDFT in the popular mixed Gaussian/plane wave electronic structure package CP2K, implementing the additional force terms arising from a constraint based on Hirshfeld charge partitioning. This improves upon the existing Becke partitioning scheme, which is prone to give unphysical atomic charges. We verify this implementation for a variety of systems: electron transfer in (H 2 O) 2 + in vacuum, electron tunnelling between oxygen vacancy centers in solid MgO, and electron self-exchange in aqueous Ru 2+ -Ru 3+ . We find good agreement with previous plane-wave CDFT results for the same systems, but at a significantly lower computational cost, and we discuss the general reliability of condensed phase CDFT calculations.

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

Computation of Auger Electron Spectra in Organic Molecules with Multiconfiguration Pair-Density Functional Theory

Efficient and accurate computation of molecular Auger electron spectra for larger systems is limited by the rapid increase in the number of doubly ionized final states as the system size grows. Here, in this work, we benchmark the application of multiconfiguration pair-density functional theory with a restricted active space (RAS) reference wave function for computing the carbon K-edge decay spectra of 20 organic molecules. Decay rates are computed within the one-center approximation. We evaluate the performance of different basis sets and on-top functionals and find that multiconfiguration pair-density functional theory achieves accuracy comparable to RAS followed by second-order perturbation theory, but at significantly lower computational cost.

Fouda, Adam E. A. [Argonne National Laboratory (AN↗
DOE OSTI2023

Non-adiabatic approximations in time-dependent density functional theory: progress and prospects

Time-dependent density functional theory continues to draw a large number of users in a wide range of fields exploring myriad applications involving electronic spectra and dynamics. Although in principle exact, the predictivity of the calculations is limited by the available approximations for the exchange-correlation functional. In particular, it is known that the exact exchange-correlation functional has memory-dependence, but in practise adiabatic approximations are used which ignore this. Here we review the development of non-adiabatic functional approximations, their impact on calculations, and challenges in developing practical and accurate memory-dependent functionals for general purposes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Importance of enforcing Hund’s rules in density functional theory calculations of rare earth magnetocrystalline anisotropy

Density functional theory (DFT) and its extensions, such as DFT+U and DFT+dynamical mean-field theory, are invaluable for studying magnetic properties in solids. However, rare-earth (R) materials remain challenging due to self-interaction errors and the lack of proper orbital polarization. We show how the orbital dependence of self-interaction error contradicts Hund’s rules and plagues magnetocrystalline anisotropy (MA) calculations, and how analyzing DFT states that respect Hund’s rules can mitigate this issue. We benchmark MA in RCo 5 , R 2 Fe 14 B, and RFe 12 , extending prior work on RMn 6 Sn 6 , achieving excellent agreement with experiments. Additionally, we illustrate a semi-analytical perturbation approach that treats crystal fields as a perturbation in the large spin-orbit coupling limit. Using Gd-4f crystal-field splitting, this method provides a microscopic understanding of MA and enables rapid screening of high-MA materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fully consistent Eulerian Monte Carlo fields method for solving probability density function transport equations in turbulence modeling

The probability density function (PDF) transport equation method is a sophisticated model for the closure of turbulent mixing and turbulent reactive flows. An efficient solution approach for solving the PDF transport equation has been vital for the method to be widely used in applications. The Eulerian Monte Carlo fields (EMCF) method has been developed to solve the PDF transport equation efficiently for decades. A recent work by Wang et al. revealed a serious issue of the EMCF method for not being fully consistent with the PDF transport equation for which the method is designed to solve. This work advances the state of the art by introducing fully consistent EMCF methods for solving the PDF transport equation. The fully consistent EMCF formulations are derived for two different PDF equation forms. The consistency of the EMCF formulations is mathematically confirmed by examining the derived moment transport equations from the EMCF formulations and from the PDF transport equation. The method of manufactured solutions is employed to further verify the consistency and convergence of the different EMCF formulations numerically. The newly introduced EMCF formulations bring the EMCF method to full consistency with the PDF transport equations for the first time.

42 ENGINEERING↗
DOE OSTI2023

Correct symmetry in density functional theory calculations of δ-Pu [Slides]

Density functional theory (DFT) calculations in the delta-phase of Pu have long relied on the antiferromagnetic (AFM) spin structure to approximate the material’s unique electronic structure. The resulting tetragonal distortion away from the correct cubic symmetry was believed to be a necessary consequence of allowing an ordered spin structure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1996

Chemistry by Way of Density Functional Theory

In this work we demonstrate that density functional theory (DFT) methods make an important contribution to understanding chemical systems and are an important additional method for the computational chemist. We report calibration calculations obtained with different functionals for the 55 G2 molecules to justify our selection of the B3LYP functional. We show that accurate geometries and vibrational frequencies obtained at the B3LYP level can be combined with traditional methods to simplify the calculation of accurate heats of formation. We illustrate the application of the B3LYP approach to a variety of chemical problems from the vibrational frequencies of polycyclic aromatic hydrocarbons to transition metal systems. We show that the B3LYP method typically performs better than the MP2 method at a significantly lower computational cost. Thus the B3LYP method allows us to extend our studies to much larger systems while maintaining a high degree of accuracy. We show that for transition metal systems, the B3LYP bond energies are typically of sufficient accuracy that they can be used to explain experimental trends and even differentiate between different experimental values. We show that for boron clusters the B3LYP energetics are not as good as for many of the other systems presented, but even in this case the B3LYP approach is able to help understand the experimental trends.

Bauschlicher, Charles W., Jr.↗
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↗
DOE OSTI2022

Sensitivity of the electronic and magnetic structures of cuprate superconductors to density functional approximations

Abstract We discuss the crystal, electronic, and magnetic structures of La 2− x Sr x CuO 4 (LSCO) for x = 0.0 and x = 0.25 employing 13 density functional approximations, representing the local, semi-local, and hybrid exchange-correlation approximations within the Perdew–Schmidt hierarchy. The meta-generalized gradient approximation (meta-GGA) class of functionals is found to perform well in capturing the key properties of LSCO, a prototypical high-temperature cuprate superconductor. In contrast, the localspin-density approximation, GGA, and the hybrid density functional fail to capture the metal-insulator transition under doping.

36 MATERIALS SCIENCE↗
DOE OSTI2020

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 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.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗