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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2026

Insights on carbon dioxide adsorption in a flexible 1-D coordination polymer from in situ X-ray scattering and density functional theory

A combination of in situ small-angle X-ray scattering (SAXS) microstructure characterization over scales ranging from 1 nm to 10 µm and powder X-ray diffraction (XRD) structure characterization under various gas/pressure/temperature conditions with density functional theory (DFT) calculations provides new insights for the CO2 sorption behavior of a one-dimensional porous coordination polymer: catena-bis­(di­benzoyl­methanato)(4,4′-bi­pyridyl)­nickel(II), denoted NiDBM-Bpy. The NiDBM-Bpy chains are held together by van der Waals forces, but the structure of guest-free NiDBM-Bpy is unsolved due to a lack of suitable crystals and high-quality powder XRD patterns. Nevertheless, SAXS and powder XRD can follow microstructural and structural changes as a function of gas pressure, composition and temperature. Both mixed-gas flow and static supercritical CO2 regimes are explored experimentally. DFT calculations are used to model the structural variation associated with XRD changes and hysteresis in the sorption isotherms. XRD and DFT calculations suggest that an orthorhombic Fddd structure with two CO2 per Ni emerges following a transition from the structure with lower molar volume and symmetry that exists without CO2 present.

Allen, Andrew J.↗
DOE OSTI2020

Determining the hydration energetics on carbon-supported Ru catalysts: An adsorption calorimetry and density functional theory study

Fundamental knowledge on the energetics at the interface between a water layer and a metal catalyst is essential so as to understand the roles that water can play in the synthesis, activation and regeneration of noble metal-based catalysts. Here, we report the direct measurement of the enthalpy of water adsorption (Δh ads ) on activated carbon (C) and activated C-supported Ru nanoparticles, which are promising catalyst as applied to the hydrogenation/hydrodeoxygenation (HDO) of oxygenates (phenolics, aldehydes, etc.). Specifically, the near-zero coverage enthalpy of water adsorption on a C-supported Ru catalyst is -75.3 ± 0.4 kJ/(mol water), suggesting favorable water–metal binding. This is much more exothermic than that on C, which has an enthalpy of adsorption of -50.3 ± 1.3 kJ/(mol water). Despite the favorable initial binding, the magnitudes of enthalpies of water condensation on C and Ru-C indicate that overall, their surfaces are both hydrophobic. Moreover, the experimentally-measured near-zero coverage water adsorption enthalpy at the Ru sites is in very good agreement with our density functional theory based calculations. At low coverages, we obtain a water binding energy of -61.7 kJ/(mol water), which increases to -78.1 kJ/(mol water) at saturation. Complementary results are also obtained from a thermal analysis, which employed a thermogravimetric analysis–differential scanning calorimetry–mass spectrometry (TG-DSC-MS), a spectroscopic investigation using ex situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and a morphological evaluation with transmission electron microscopy (TEM). We point out that in carbon-supported metal catalysts, such as Ru-C, a strong hydration at near-zero coverage and relative weak water-surface interactions occurs upon saturation. Such heterogeneity is essential and crucial for catalytic hydrogenation/HDO reactions that involve balanced interactions among the water-rich reactant mixture and nonpolar organic products.

09 BIOMASS FUELS↗
DOE OSTI2023

Predicting Mössbauer Parameters of Nonheme Diiron Complexes with Density Functional Theory

Mössbauer spectroscopy provides significant insights into the electronic structure and environment of the metal centers. In this work, we investigate the electronic structures of a set of nonheme diiron complexes by evaluating two key parameters pertaining to Mössbauer spectroscopy, namely, the isomer shift (δ) and quadrupole splitting (|ΔE Q |), using different levels of density functional theory (DFT). The diiron systems investigated here span diverse oxidation states, bridging motifs, and spin coupling patterns, which present a challenging case for theoretical predictions. We demonstrate that the combination of B97-D3/def2-TZVP is an efficient approach in modeling both the δ and |ΔEQ| values with high accuracy for the representative nonheme diiron complexes. We also show that δ is accurately predicted irrespective of the choice of approximate density functional while the |ΔE Q | is sensitive to the level of theory employed. Further investigation shows that the present methodology assessed using synthetic nonheme diiron complexes could be extended to nonheme diiron enzyme active sites, featuring both ferromagnetic and antiferromagnetic coupling between the iron centers.

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

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid–amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Density-related properties from self-interaction corrected density functional theory calculations

The Perdew–Zunger self-interaction correction (PZ-SIC) removes unphysical electron self-interaction from calculations employing standard density functional approximations. Doing so improves many computed properties, bringing them into better agreement with experimental observations or with results from high-level quantum chemistry calculations. However, while PZ-SIC generally corrects in the right direction relative to the corresponding reference values, in many cases, it over-corrects. For this reason, scaled-down versions of PZ-SIC have been proposed and investigated. These approaches have mostly employed exterior scaling in which SIC correction terms are scaled in the same way at every point in space. Recently, a new local, or interior, scaling SIC method was proposed on non-empirical grounds to restore a property of the exact, but unknown, density functional that is broken in PZ-SIC. Here, the scaling at each point depends on the character of the charge density at that point. However, the local scaling can be done in various ways while still restoring the behavior of the exact functional. In this work, we compare and contrast the performance of various interior scaling approaches for addressing over-corrections of calculated molecular dipole moments and atomic polarizabilities and properties that reflect the nature of the electronic charge density.

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

Chemical bond effects in classical site density functional theory of inhomogeneous molecular liquids

Intra-molecular interactions or chemical bonds represent one of the main distinguishing characteristic of molecular fluids. Development of accurate and practical methods to treat these effects is one of the long standing problems in classical site density functional theory (SDFT). One particular instance when these issues become particularly severe is the case of classical interactions potentials with auxiliary sites or dummy atoms. In this situation current SDFT implementations, such as three-dimensional reference interaction site model (3D-RISM), lead to nonphysical results. We re-examine this issue in this work using our recent reformulation of SDFT. We put forward a simple practical solution to this problem, and illustrate its utility for the case of spherical solutes in diatomic liquids.

Chuev, Gennady N.↗
DOE OSTI2020

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO + , and the N K-edge of NO 2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. Finally, we also present general guidelines for computing core-excited states from orbital optimized DFT.

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

Electronic structure of Alq 3 and Liq using soft X-ray spectroscopy and density functional theory calculation

We show the electronic structure of aluminum tris-8-hydroxyquinoline (Alq 3 ) and 8-hydroquinolatolithium (Liq) was investigated using a combination of X-ray emission spectroscopy (XES), X-ray absorption spectroscopy (XAS), X-ray photoelectron spectroscopy (XPS), and spectral simulation with density functional theory. The chemical bonding states of Alq 3 and Liq were analyzed using XPS core-level spectra. The band gaps of Alq 3 and Liq were measured by combining the XAS and XES results. Additionally, resonant and non-resonant XES were used to measure the density of states for O sites in the molecules.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗