Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory development”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Ab initio Molecular Dynamics Beyond Density Functional Theory

The funding from this award supported the development of new theoretical methods that both increase the accuracy of quantum simulations beyond that of density functional theory, as well as reduce the cost of such simulations, for example, through machine learning and algorithmic improvements. The improvements lead to, in some cases, orders of magnitude increases in simulation speed as well as orders of magnitude increase in accuracy, compared to before this work was started. This translates to both greater certainty in making predictions about chemical phenomena that are currently studied computationally, as well as laying the ground-work for the study of new areas of chemistry that have not traditionally been simulated, because of the new time-scales and accuracies that can be reached.

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

Accelerating the development of new solar absorbers by photoemission characterization coupled with density functional theory

Abstract The expectation to progress towards Terawatts production by solar technologies requires continuous development of new materials to improve efficiency and lower the cost of devices beyond what is currently available at industrial level. At the same time, the turnaround time to make the investment worthwhile is progressively shrinking. Whereas traditional absorbers have developed in a timeframe spanning decades, there is an expectation that emerging materials will be converted into industrially relevant reality in a much shorter timeframe. Thus, it becomes necessary to develop new approaches and techniques that could accelerate decision-making steps on whether further research on a material is worth pursuing or not. In this review, we will provide an overview of the photoemission characterization methods and theoretical approaches that have been developed in the past decades to accelerate the transfer of emerging solar absorbers into efficient devices.

14 SOLAR ENERGY↗
DOE OSTI2022

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

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

Understanding Twinning and Deformation in High Entropy Alloys

A combination of high strength and high ductility has been observed in multi-principal element alloys due to twin formation attributed to low stacking fault energy (SFE). In the pursuit of low SFE alloys, a key bottleneck is the lack of understanding of the composition–SFE cor- relations that would guide tailoring SFE via alloy composition. Using density functional theory (DFT), we show that dopant radius, which have been postulated as a key descriptor for SFE in dilute alloys, does not fully explain SFE trends across different host metals. Instead, charge density is a much more central descriptor. It allows us to (1) explain contrasting SFE trends in Ni and Cu host metals due to various dopants in dilute concentrations, (2) explain the large SFE variations observed in the literature even within a given alloy composition due to the nearest neighbor environments in “model” concentrated alloys, and (3) develop a machine learning model that can be used to predict SFEs in multi-elemental alloys. This model opens a possibility to use charge density as a descriptor for predicting SFE in alloys. Furthermore, a descriptor-less machine learning (ML) model based only on charge density images extracted from density functional theory (DFT) is developed to predict stacking fault energies (SFE) in concentrated alloys. The model is based on convolutional neural networks (CNNs) as one of the promising ML techniques for dealing with complex images and data. Identification of correct descriptors is a key bottleneck to develop ML models for predicting materials properties. Often, in most ML models, textbook physical descriptors such as atomic radius, valence charge and electronegativity are used as descriptors which have limitations because these properties change in concentrated alloys when multiple elements are mixed to form a solid solution. We illustrate that, within the scope of DFT, the search for descriptors can be circumvented by electronic charge density, which is the backbone of the Kohn-Sham DFT and describes the system completely. The performance of our model is demonstrated by predicting SFE of concentrated alloys with an RMSE and R2 of 6.18 mJ/m2 and 0.87, respectively, validating the accuracy of the proposed approach.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accelerating Multiscale Materials Modeling with Machine Learning

The focus of this project is to accelerate and transform the workflow of multiscale materials modeling by developing an integrated toolchain seamlessly combining DFT, SNAP, LAMMPS, (shown in Figure 1-1) and a machine-learning (ML) model that will more efficiently extract information from a smaller set of first-principles calculations. Our ML model enables us to accelerate first-principles data generation by interpolating existing high fidelity data, and extend the simulation scale by extrapolating high fidelity data (10 2 atoms) to the mesoscale (10 4 atoms). It encodes the underlying physics of atomic interactions on the microscopic scale by adapting a variety of ML techniques such as deep neural networks (DNNs), and graph neural networks (GNNs). We developed a new surrogate model for density functional theory using deep neural networks. The developed ML surrogate is demonstrated in a workflow to generate accurate band energies, total energies, and density of the 298K and 933K Aluminum systems. Furthermore, the models can be used to predict the quantities of interest for systems with more number of atoms than the training data set. We have demonstrated that the ML model can be used to compute the quantities of interest for systems with 100,000 Al atoms. When compared with 2000 Al system the new surrogate model is as accurate as DFT, but three orders of magnitude faster. We also explored optimal experimental design techniques to choose the training data and novel Graph Neural Networks to train on smaller data sets. These are promising methods that need to be explored in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optical properties of a diamond NV color center from capped embedded multiconfigurational correlated wavefunction theory

Diamond defects are among the most promising qubits. Modeling their properties through accurate quantum mechanical simulations can further their development into robust units of information. We use the recently developed capped density functional embedding theory (capped-DFET) with the multiconfigurational n-electron valence second-order perturbation theory to characterize the electronic excitation energies for different spin manifolds of the well-characterized negatively charged substitutional N defect adjacent to a vacancy (V C ) in diamond (N C V C − ). We successfully reproduce vertical excitation energies for both triplet and singlet states of N C V C − with errors < 0.1 eV. Unlike other embedding methods, capped-DFET exhibits robust predictions that are approximately independent of the embedded cluster size: it only requires a cluster to contain the defect atoms and their nearest neighbors (as small as a 40-atom capped cluster). Furthermore, our method is free from slowly converging Coulomb interactions between charged defects, and thus also only weakly dependent on supercell size.

Chemistry↗
DOE OSTI2023

Insights into the Oxygen Evolution Reaction on Graphene-Based Single-Atom Catalysts from First-Principles-Informed Microkinetic Modeling

Single-atom transition metals embedded in nitrogen-doped graphene have emerged as promising electrocatalysts due to their high activity and low material cost. These materials have been shown to catalyze a variety of electrochemical reactions, but their active sites under reaction conditions remain poorly understood. Using first-principles density functional theory calculations, we develop a pH-dependent microkinetic model to evaluate the relative performance of transition metal catalysts embedded in fourfold N-substituted double carbon vacancies in graphene for the oxygen evolution reaction. We find that reaction pathways involving intermediates co-adsorbed on the metal site are preferred on all transition metals. These pathways lead to enhancements in catalytic activity and broaden the activity peak when compared with purely thermodynamics-based predictions. Furthermore, these findings demonstrate the importance of investigating reaction pathways on graphene-based catalysts and other two-dimensional (2D) materials that involve metal active centers decorated by spectator intermediate species.

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

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Charge-density based convolutional neural networks for stacking fault energy prediction in concentrated alloys

A descriptor-less machine learning (ML) model based only on charge density images extracted from density functional theory (DFT) is developed to predict stacking fault energies (SFE) in concentrated alloys. The model is based on convolutional neural networks (CNNs) as one of the promising ML techniques for dealing with complex images and data. Identification of correct descriptors is a key bottleneck to develop ML models for predicting materials properties. Often, in most ML models, textbook physical descriptors such as atomic radius, valence charge and electronegativity are used as descriptors which have limitations because these properties change in concentrated alloys when multiple elements are mixed to form a solid solution. Here, we illustrate that, within the scope of DFT, the search for descriptors can be circumvented by electronic charge density, which is the backbone of the Kohn-Sham DFT and describes the system completely. The performance of our model is demonstrated by predicting SFE of concentrated alloys with an RMSE and R 2 of 6.18 mJ/m 2 and 0.87, respectively, validating the accuracy of the proposed approach.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Cluster model of C 12 in the density functional theory framework

We employ the constrained density functional theory to investigate cluster phenomena for the 12 C nucleus. The proton and neutron densities are generated from the placement of three 4 He nuclei (α particles) geometrically. These densities are then used in a density constrained Hartree-Fock calculation that produces an antisymmetrized state with the same densities through energy minimization. In the calculations no a priori analytic form for the single-particle states is assumed and the full energy density functional is utilized. The geometrical scan of the energy landscape provides the ground state of 12 C as an equilateral triangular configuration of three αs with molecular bond like structures. The use of the nucleon localization function provides further insight to these configurations. One can conclude that these configurations are a hybrid between a pure mean-field and a pure α particle condensate. Furthermore, this development could facilitate density functional theory based fusion calculations with a more realistic 12 C ground state.

6 ≤ A ≤ 19↗
DOE OSTI2020

Two-Dimensional Palladium Diselenide with Strong In-Plane Optical Anisotropy and High Mobility Grown by Chemical Vapor Deposition

Two-dimensional (2D) palladium diselenide (PdSe 2 ) has strong interlayer coupling and a puckered pentagonal structure, leading to remarkable layer-dependent electronic structures and highly anisotropic in-plane optical and electronic properties. However, the lack of high-quality, 2D PdSe 2 crystals grown by bottom-up approaches limits the study of their exotic properties and practical applications. In this work, chemical vapor deposition growth of highly crystalline few-layer (≥2 layers) PdSe 2 crystals on various substrates is reported. The high quality of the PdSe 2 crystals is confirmed by low-frequency Raman spectroscopy, scanning transmission electron microscopy, and electrical characterization. In addition, strong in-plane optical anisotropy is demonstrated via polarized Raman spectroscopy and second-harmonic generation maps of the PdSe 2 flakes. A theoretical model based on kinetic Wulff construction theory and density functional theory calculations is developed and described the observed evolution of “square=like” shaped PdSe 2 crystals into rhombus due to the higher nucleation barriers for stable attachment on the (1,1) and (1,-1) edges, which results in their slower growth rates. Finally, few-layer PdSe 2 field-effect transistors reveal tunable ambipolar charge carrier conduction with an electron mobility up to ≈294 cm 2 V -1 s -1 , which is comparable to that of exfoliated PdSe 2 , indicating the promise of this anisotropic 2D material for electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accelerated computation of lattice thermal conductivity using neural network interatomic potentials

We report with the development of the density functional theory (DFT) and ever-increasing computational capacity, an accurate prediction of lattice thermal conductivity based on the Boltzmann transport theory becomes computationally feasible, contributing to a fundamental understanding of thermal conductivity as well as a choice of the optimal materials for specific applications. However, steep costs in evaluating third-order force constants limit the theoretical investigation to crystals with high symmetry and few atoms in the unit cell. Currently, machine learning potentials are garnering attention as a computationally efficient high-fidelity model of DFT, and several studies have demonstrated that the lattice thermal conductivity could be computed accurately via the machine learning potentials. However, test materials were mostly crystals with high symmetries, and the applicability of machine learning potentials to a wide range of materials has yet to be demonstrated. Furthermore, establishing a standard training set that provides consistent accuracy and computational efficiencies across a wide range of materials would be useful. To address these issues, herein we compute lattice thermal conductivities at 300 K using neural network interatomic potentials. As test materials, we select 25 materials with diverse symmetries and a wide range of lattice thermal conductivities between 10 -1 and 10 3 Wm -1 K -1 . Among various choices of training sets, we find that molecular dynamics trajectories at 50–700 K consistently provide results at par with DFT for the test materials. In contrast to pure DFT approaches, the computational cost in the present approach is uniform over the test materials, yielding a speed gain of 2–10 folds. When a smaller reduced training set is used, the relative efficiency increases by up to ~50 folds without sacrificing accuracy significantly. The current work will broaden the application scope of machine learning potentials by establishing a robust framework for accurately computing lattice thermal conductivity with machine learning potentials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

One dimensional wormhole corrosion in metals

Corrosion is a ubiquitous failure mode of materials. Often, the progression of localized corrosion is accompanied by the evolution of porosity in materials previously reported to be either three-dimensional or two-dimensional. However, using new tools and analysis techniques, we have realized that a more localized form of corrosion, which we call 1D wormhole corrosion, has previously been miscategorized in some situations. Using electron tomography, we show multiple examples of this 1D and percolating morphology. To understand the origin of this mechanism in a Ni-Cr alloy corroded by molten salt, we combined energy-filtered four-dimensional scanning transmission electron microscopy and ab initio density functional theory calculations to develop a vacancy mapping method with nanometer-resolution, identifying a remarkably high vacancy concentration in the diffusion-induced grain boundary migration zone, up to 100 times the equilibrium value at the melting point. Deciphering the origins of 1D corrosion is an important step towards designing structural materials with enhanced corrosion resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-throughput dataset of impurity adsorption on common catalysts in biomass upgrading applications

Abstract An extensive dataset consisting of adsorption energies of pernicious impurities present in biomass upgrading processes on common catalysts and support materials has been generated. This work aims to inform catalyst and process development for the conversion of biomass-derived feedstocks to fuels and chemicals. A high-throughput workflow was developed to execute density functional theory calculations for a diverse set of atomic (Al, B, Ca, Cl, Fe, K, Mg, Mn, N, Na, P, S, Si, Zn) and molecular (COS, H 2 S, HCl, HCN, K 2 O, KCl, NH 3 ) species on 35 unique surfaces for transition-metal (Ag, Au, Co, Cu, Fe, Ir, Ni, Pd, Pt, Re, Rh, Ru) and metal-oxide (Al 2 O 3 , MgO, anatase-TiO 2 , rutile-TiO 2 , ZnO, ZrO 2 ) catalysts and supports. Approximately 3,000 unique adsorption geometries and corresponding adsorption energies were obtained.

09 BIOMASS FUELS↗
DOE OSTI2024

First-principles study of the magneto-Raman effect in van der Waals layered magnets

Magneto-Raman spectroscopy has been used to study spin-phonon coupling in two-dimensional (2D) magnets. Raman spectra of CrI 3 show a strong dependence on the magnetic order within a layer and between the layers. Here we carry out the first systematic theoretical investigation of the magneto-Raman effect in 2D magnets by performing density functional theory calculations and developing a generalized polarizability model. Our first-principles simulations well reproduce experimental Raman spectra of CrI 3 with different magnetic states. The model reveals how the change of spin orientation in each layer is coupled to the layer’s vibration to induce or eliminate the spin-dependent anti-symmetric off-diagonal terms in the Raman tensor for altering the selection rules. We also uncover that the correlation between phonon modes and magnetic orders is a universal phenomenon, which should exist in other phonon modes and 2D magnets. Our predictive simulations and modeling are expected to guide the research in 2D magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗
DOE OSTI2021

Determining the electrical conductivity of metals using the 2 MA Thor pulsed power driver

We present the development of a pulsed power experimental technique to infer the electrical conductivity of metals from ambient to high energy density conditions. The method is implemented on Thor, a moderate scale (1–2 MA) pulsed power driver. The electrical conductivity of copper at elevated temperature (>4000 K) and pressure (>10 GPa) is determined, and a new tabular material model is developed, guided by density functional theory, which preserves agreement with existing experimental data. Minor modifications (<10%) are found to be necessary to the previous Lee–More–Desjarlais model isotherms in the vicinity of the melt transition in order to account for observed discrepancies with the new experimental data. An analytical model for magnetic direct drive flyer acceleration and Joule heating induced vaporization based on the Tsiolkovsky “rocket equation” is presented to assess sensitivity of the method to minor changes in electrical conductivity.

47 OTHER INSTRUMENTATION↗