Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
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 OSTI2019

VeloxChem: A Python-driven density-functional theory program for spectroscopy simulations in high-performance computing environments

An open-source program named VeloxChem has been developed for the calculation of electronic real and complex linear response functions at the levels of Hartree–Fock and Kohn–Sham density functional theories. With an object-oriented program structure written in a Python/C++ layered fashion, VeloxChem enables time-efficient prototyping of novel scientific approaches without sacrificing computational efficiency, so that molecular systems involving up to and beyond 500 second-row atoms (or some 10,000 contracted and in part diffuse Gaussian basis functions) can be routinely addressed. In addition, VeloxChem is equipped with a polarizable embedding scheme for the treatment of the classical electrostatic interactions with an environment that in turn is modeled by atomic site charges and polarizabilities. The underlying hybrid message passing interface (MPI)/open multiprocessing (OpenMP) parallelization scheme makes VeloxChem suitable for execution in high-performance computing cluster environments, showing even slightly beyond linear scaling for the Fock matrix construction with use of up to 16,384 central processing unit (CPU) cores. An efficient—with respect to convergence rate and overall computational cost—multifrequency/gradient complex linear response equation solver enables calculations not only of conventional spectra, such as visible/ultraviolet/X-ray electronic absorption and circular dichroism spectra, but also time-resolved linear response signals as due to ultra-short weak laser pulses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A noncollinear density functional theory ansatz for the phononic and thermodynamic properties of α -Pu

Plutonium's phase diagram is host to complex structures and interactions that make the description of its ground state properties elusive. Using all-electron density functional theory calculations, we study the thermodynamic properties of α-Pu. To do this, we build on recent work in the literature by introducing a novel noncollinear magnetic ansatz for α-Pu's ground state. We show that the noncollinear ansatz improves the description of the experimental phonon density of states, heat capacity, and thermal expansion over previous results obtained with a collinear ansatz. Here, despite lacking a bond-equivalent picture of the bonding, these new results for α-Pu, along with recent results on δ-Pu, demonstrate the efficacy of noncollinear ansatzes for the description of plutonium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Radical Coupling Reactions of Hydroxystilbene Glucosides and Coniferyl Alcohol: A Density Functional Theory Study

The monolignols, p-coumaryl, coniferyl, and sinapyl alcohol, arise from the general phenylpropanoid biosynthetic pathway. Increasingly, however, authentic lignin monomers derived from outside this process are being identified and found to be fully incorporated into the lignin polymer. Among them, hydroxystilbene glucosides, which are produced through a hybrid process that combines the phenylpropanoid and acetate/malonate pathways, have been experimentally detected in the bark lignin of Norway spruce (Picea abies). Several interunit linkages have been identified and proposed to occur through homo-coupling of the hydroxystilbene glucosides and their cross-coupling with coniferyl alcohol. In the current work, the thermodynamics of these coupling modes and subsequent rearomatization reactions have been evaluated by the application of density functional theory (DFT) calculations. The objective of this paper is to determine favorable coupling and cross-coupling modes to help explain the experimental observations and attempt to predict other favorable pathways that might be further elucidated via in vitro polymerization aided by synthetic models and detailed structural studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Point defect energetics in gallium arsenide, a comprehensive density functional theory study

In materials, point defects often control or modify functional properties. To predict the performance of materials intended for application in optoelectronic devices, it is imperative to understand the properties of those point defects. For the first time, all six intrinsic defects of GaAs, a key optoelectronics material, and their charge transition levels are calculated using density functional theory with the HSE06 functional. For comparison, both PBE and r 2 SCAN calculations are also carried out. The HSE06 results are found to be in better agreement with experimental data than previous calculations. In conclusion, the importance of using the exact electron exchange present in hybrid functionals and larger supercells to accurately determine defect levels and ground state defect configurations is demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accurate frozen core approximation for all-electron density-functional theory

We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over twofold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the Periodic Table and a large-scale simulation of CsPbBr3 with 2560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below −200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.

Yu, Victor Wen-zhe↗
DOE OSTI2020

Time-dependent density-functional theory of high-order harmonic generation from noble-gas atoms driven by orthogonally polarized two-color laser fields

An orthogonally polarized two-color (OTC) field of frequency ratio 1 : 2 is known to generate both odd and even harmonic spectra in the mutually orthogonal directions of the fundamental and the second-harmonic fields, respectively. Recollision angles of an ionized electron driven back to its parent ion by an OTC field are inferred from the intensity ratio of odd and even harmonics, where the relative phase delays between the two-color fields can be used in time-resolved measurements on a subcycle time scale. Here, we present the ab initio, many-electron calculation of time-resolved recollision angles of noble-gas atoms driven by an OTC field, based on the time-dependent density-functional theory. Our results agree well with experimental data by Shafir et al. in Nat. Phys. 5, 412 (2009) of helium atom and New J. Phys. 12, 073032 (2010) of neon atom.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Discovery of single-atom alloy catalysts for CO 2 -to-methanol reaction by density functional theory calculations

The transformations of CO 2 molecules into valuable products are of increasing interest due to the negative impact of anthropogenic CO 2 emissions on global warming. The CO 2 -to-methanol hydrogenation is an economically profitable reaction of carbon fixation, but it still steps away from widespread industrialization because of the lack of efficient and selective catalysts. Recently, single-atom alloy (SAA) catalysts have been developed to work remarkably in CO 2 hydrogenation reactions. Doping isolated single atoms into metallic catalyst can dramatically alter the catalytic performance of the host. Here, we have performed a screening discovery on Ru and 6 RuX (X = Fe, Co, Ni, Cu, Ir and Pt) SAAs using density functional theory (DFT) computations. We considered 13 possible elementary reactions in 4 possible reaction pathways on Ru and all RuX surfaces. In the computed mechanisms, we found that the formation of *H 2 COOH and *HCOO intermediates plays a critical role in determining catalysts’ activities. Doping Co and Pt isolated single atoms into Ru surface can thermodynamically and kinetically facilitate these intermediates formation processes, eventually promoting the production of methanol. The combination of weak binding and enhanced charge redistribution on RuCo and RuPt surfaces give them improved catalytic activities over pure Ru. This work will ultimately facilitate the discovery and development of SAAs for CO 2 to methanol, serving as guidance to experiments and theoreticians alike.

25 ENERGY STORAGE↗
DOE OSTI2020

The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules

Abstract Maximum diversification of data is a central theme in building generalized and accurate machine learning (ML) models. In chemistry, ML has been used to develop models for predicting molecular properties, for example quantum mechanics (QM) calculated potential energy surfaces and atomic charge models. The ANI-1x and ANI-1ccx ML-based general-purpose potentials for organic molecules were developed through active learning; an automated data diversification process. Here, we describe the ANI-1x and ANI-1ccx data sets. To demonstrate data diversity, we visualize it with a dimensionality reduction scheme, and contrast against existing data sets. The ANI-1x data set contains multiple QM properties from 5 M density functional theory calculations, while the ANI-1ccx data set contains 500 k data points obtained with an accurate CCSD(T)/CBS extrapolation. Approximately 14 million CPU core-hours were expended to generate this data. Multiple QM calculated properties for the chemical elements C, H, N, and O are provided: energies, atomic forces, multipole moments, atomic charges, etc. We provide this data to the community to aid research and development of ML models for chemistry.

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

Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

In-situ Synthesis of PuCl 3 and Corresponding Raman and Density Functional Theory Vibrational Modes

The experimental and calculated Raman spectrum of PuCl 3 has been reported for the first time. PuCl 3 is a primary species found in plutonium metal refinement, specifically in pyrochemical salt processes including multicycle direct oxide reduction, metal chlorination, and electrorefining. As such, Raman signatures of PuCl 3 could serve as potential forensic indicators of material process history. A novel technique for synthesizing PuCl 3 from the in situ chlorination of plutonium metal with HCl was developed to establish these signatures. Cerium metal surrogates were utilized to ensure optimization of the plutonium experiments and to minimize personnel exposure, and all experiments were carried out in a Raman reaction chamber designed for air-tight, high vacuum environments. In situ Raman spectroscopy was employed in conjunction with density functional theory (DFT) to investigate the vibrational modes of PuCl 3 . Associated mixed oxy and hydroxyl phases are also reported. The combined Raman and DFT results have eliminated inconsistencies in Raman mode assignments for the MX 3 family of metal chlorides having P6 3 /m symmetry, and IR modes derived from the DFT calculations are additionally presented. The data observed in this study are of potential interest to nuclear forensic analyses, nuclear sample aging, nuclear energy, plutonium processing, and stockpile stewardship.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiconfiguration Pair-Density Functional Theory for Vertical Excitation Energies in Actinide Molecules

Modeling actinides with electronic structure theories is challenging because these systems present a strong ligand field and metal–ligand covalency. Here, we systematically investigate the effectiveness of pair-density functional theory (PDFT) for the calculation of vertical excitation energies in An(III), [An III Cl 6 ] 3– , and [An VI O 2 ] 2+ (An = U, Np, Pu, and Am). We compare the performance of PDFT, hybrid PDFT, and multistate PDFT with traditional active-space methods followed by perturbation theory, like multistate CASPT2, and with experimental data. Overall, multistate PDFT gives quantitative agreement with multistate CASPT2 at a significantly reduced computational cost.

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

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

In this work, we present DFT-FE 1.0, building on DFT-FE 0.6 [Comput. Phys. Commun. 246, 106853 (2020)], to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ~ 100,000 electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation—via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency—as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our implementation, which yields ~ 20× CPU-GPU speed-up by using GPU acceleration on hybrid CPU-GPU nodes. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of 80–140 seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ~ 6, 000 – 15,000 electrons.

pseudopotential↗
DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Building surrogate models of nuclear density functional theory with Gaussian processes and autoencoders

From the lightest Hydrogen isotopes up to the recently synthesized Oganesson (Z = 118), it is estimated that as many as about 8,000 atomic nuclei could exist in nature. Most of these nuclei are too short-lived to be occurring on Earth, but they play an essential role in astrophysical events such as supernova explosions or neutron star mergers that are presumed to be at the origin of most heavy elements in the Universe. Understanding the structure, reactions, and decays of nuclei across the entire chart of nuclides is an enormous challenge because of the experimental difficulties in measuring properties of interest in such fleeting objects and the theoretical and computational issues of simulating strongly-interacting quantum many-body systems. Nuclear density functional theory (DFT) is a fully microscopic theoretical framework which has the potential of providing such a quantitatively accurate description of nuclear properties for every nucleus in the chart of nuclides. Thanks to high-performance computing facilities, it has already been successfully applied to predict nuclear masses, global patterns of radioactive decay like β or γ decay, and several aspects of the nuclear fission process such as, e.g., spontaneous fission half-lives. Yet, predictive simulations of nuclear spectroscopy—the low-lying excited states and transitions between them—or of nuclear fission, or the quantification of theoretical uncertainties and their propagation to basic or applied nuclear science applications, would require several orders of magnitude more calculations than currently possible. However, most of this computational effort would be spent into generating a suitable basis of DFT wavefunctions. Such a task could potentially be considerably accelerated by borrowing tools from the field of machine learning and artificial intelligence. In this paper, we review different approaches to applying supervised and unsupervised learning techniques to nuclear DFT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Hole Polaron Transport in Bismuth Vanadate BiVO 4 from Hybrid Density Functional Theory

In this paper, we address the chemical character and the mobility of hole polarons in BiVO 4 (BVO). Two distinct structures for hole polarons in BVO have been reported in the literature to date. “Standard” hybrid density functional theory (DFT) calculations have predicted hole polarons to be (BiO 8 ) dodecahedron-centered polarons, while Hubbard U-DFT (DFT + U) calculations have predicted single oxygen atom-centered polarons. Resolving these contrasting findings is critical to validate theory in materials research. We investigated the structure of hole polarons in BVO using hybrid DFT with varying fractions (α) of the exact exchange interaction that enters hybrid functionals. For values of the fraction α, from α = 0.25 to α = 0.45, we obtained both structures, the h + (BiO 8 ) hole structure and the h + (O) hole structure. For the smaller values of α, the h + (BiO 8 ) structure is lower in energy, while for the larger values of α, the h+(O) structure is more stable. The piece-wise linearity of DFT energies E(q) with respect to a point defect charge q is maintained satisfactorily in all cases, although the deviation is the smallest for hybrid DFT with a small exact exchange fraction. For α = 0.25, we determined the diffusivity of h + (BiO 8 ) holes from the Marcus/Holstein two-state model. The activation barrier for Bi-to-Bi hops (ΔG # ~76 meV) and mobility (μ ~0.029 cm2 V –1 s –1 ) are found to be in close agreement with the published THz experimental data (ΔG # ~90 meV and μ exp ~0.02 cm2 V –1 s –1 , respectively). In h + (BiO 8 ) structures, nearly 25% of the hole charge and unpaired spin reside on the Bi atom and the rest is divided among the eight O atoms. The hole is best described as a Bi-centered hole, rather than a Bi 6s hole. In contrast, in h + (O) structures, about 75% of the unpaired spin population resides on the O atom. These findings suggest that if one wishes to resolve this theoretical dichotomy, there remains a need for further experimental characterization of hole polarons in BVO with techniques that could identify the oxidation state of Bi and O atoms in operando conditions and help validate theoretical results of structure and transport.

14 SOLAR ENERGY↗