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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2023

Atomistic theories and Simulations of Multiferroics

The goal of this award was to build a research program aimed at investigating puzzling complex phenomena in multiferroics, that are systems possessing electric and magnetic dipoles. We first emphasized the study of multiferroic nanostructures, that are 0D, 1D, 2D systems as well as superlattices and then also study multiferroics as a whole, i.e., bulks as well. In order to fully understand such complex systems, we also decided to study magnetic systems and ferroelectrics, both in their nanostructure and bulk forms. For that and in addition to use standard Density Functional Theory codes, we developed and used ab-initio tools that are able to model properties of large supercells (about a million atoms nowadays) at finite temperature and for any mechanical and boundary conditions. We also collaborated with experimentalists to fully understand the systems under study. The whole program resulted in 108 publications that are listed at the end of this report. In the following, I will put an emphasis on some publications before the last award period and also indicate a description of the last award period.

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

A deep learning interatomic potential developed for atomistic simulation of carbon materials

Interatomic potentials based on neural-network machine learning method have attracted considerable attention in recent years owing to their outstanding ability to balance the accuracy and efficiency in atomistic simulations. In this work, a neural-network potential (NNP) for carbon is generated to simulate the structural properties of various carbon structures. The potential is trained using a database consisting of crystalline and liquid structures obtained by the first-principles density functional theory (DFT) calculations. The developed potential accurately predicts the energies and forces in crystalline and liquid carbon structures, the energetic stability of defected graphene, and the structures of amorphous carbon as the function of density. As a result, the excellent accuracy and transferability of the NNP provide a promising tool for accurate atomistic simulations of various carbon materials with faster speed and much lower cost.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

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

Final Technical Report for project entitled Highly Active, Durable, and Ultra-Low PGM NSTF Thin Film ORR Catalysts and Support

In this project, the objective was to develop new oxygen reduction reaction (ORR) electrocatalysts for proton exchange membrane fuel cells (PEMFCs) which could exceed all of the Department of Energy (DOE) 2020 targets listed in DE-FOA-0001224, Subtopic 1b, Table I. The expected outcome was development of one or more electrocatalysts which are substantially improved over the current state-of-the-art in terms of overall activity, durability, and cost, and are suitable for automotive traction and stationary fuel cell applications. This report summarizes this project’s progress towards meeting the stated objectives and includes a summary of key findings and conclusions. The development activities towards new highly active and durable thin film electrocatalysts has led to the discovery of several electrocatalysts which approached or exceeded several DOE 2020 targets for activity, durability, and performance in proton exchange membrane fuel cell (membrane electrode assemblies (MEAs). In this project, electrocatalyst development focused on systematic physical and electrochemical characterization of electrocatalyst activity, durability, and performance in MEAs as a function of electrocatalyst fabrication, compositional and structural variables. The project approach was to develop relationships between the catalyst functional responses (activity, durability and performance) with catalyst physical properties and catalyst fabrication methods. Additionally, the thin film catalysts were integrated onto the unique and durable 3M Nanostructured Thin Film (NSTF) support, consisting of arrays of self-assembled organic crystalline whiskers. The catalysts were evaluated against the project targets via extensive electrochemical characterization in MEA format at 3M, advanced structural and compositional microscopy at Oak Ridge National Laboratory (ORNL), and atomic structure analysis via XAFS at Argonne National Laboratory (ANL). The extensive electrochemical and physical characterization resulted in development of several trends which correlate the catalysts’ electrochemical properties to their physical properties, including composition, structure, and method of fabrication. The development was guided by density functional theory (DFT) modeling at Purdue University and kinetic Monte Carlo (kMC) modeling at Johns Hopkins University. The catalyst simulations provided key insights into the observed experimental activity and durability trends, and additionally were utilized to assess new electrocatalyst concepts prior to or coincident with physical catalyst development. This development has led to several electrocatalyst candidates with activity and durability which approach or exceed DOE targets. One class of catalysts, based on nm-scale thin layers of Pt on Ir, met or exceeded 6 of the 6 DOE targets the project addressed.

08 HYDROGEN↗
DOE OSTI2021

Coupling of Flavonoid Initiation Sites with Monolignols Studied by Density Functional Theory

Lignin recalcitrance presents a challenge for the development of a bioeconomy that employs lignocellulosic feedstocks. The efficiency of lignin deconstruction is improved by a reduction of molecular weight, and given the discovery that flavones serve as initiation sites in lignin biosynthesis, these molecular weight reductions could potentially be achieved with plant metabolic engineering to over-express flavonoids. Upon increasing the flavonoid content in lignin, the bond strengths and properties of flavonoid-monolignol linkages become increasingly important. To that end, the current work applies density functional theory calculations to elucidate the bond dissociation enthalpies (BDEs) of flavonoid-monolignol linkages, including dimers with oxidized monolignols. Specifically, the dimer bond strengths and monomer hydrogen abstraction energies for the flavonoids tricin, chrysoeriol, luteolin, apigenin, catechin, epicatechin, epigallocatechin, and epigallocatechin gallate are calculated, when coupled to seven natural and engineered monolignols. Results indicate that 4'-O-ß linkage strengths between flavonoids and monolignols are of comparable strength to inter-monolignol ß-O-4 linkages, with average flavonoid-monolignol BDEs of 70.7 kcal/mol relative to ~69.3 kcal/mol in analogous canonical monolignols. Epigallocatechin yielded the lowest 4'-O-ß bond strength of 52.3 kcal/mol when coupled to an oxidized monolignol, while the flavones overall produced lower average BDEs, relative to the flavanols. Substituents at the 3'-C and 5'-C positions on flavonoids affected the dimer linkage strengths to a greater extent than glycosylation or substituents further from the linkage. Erythro and threo stereochemistry across the flavonoid-monolignol linkage library exhibited only small energetic differences and no pronounced correlations. Taken together, the predictions from this work support the concept that higher concentrations of flavonoid initiation sites in lignin may afford linkage properties conducive to more facile lignin depolymerization.

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

Towards atom-level understanding of metal oxide catalysts for the oxygen evolution reaction with machine learning

Abstract Green hydrogen production is crucial for a sustainable future, but current catalysts for the oxygen evolution reaction (OER) suffer from slow kinetics, despite many efforts to produce optimal designs, particularly through the calculation of descriptors for activity. In this study, we develop a dataset of density functional theory calculations of bulk and surface perovskite oxides, and adsorption energies of OER intermediates, which includes compositions up to quaternary and facets up to (555). We demonstrate that per-site properties of perovskite oxides such as Bader charge or band center can be tuned through element substitution and faceting, and develop a machine learning model that accurately predicts these properties directly from the local chemical environment. We leverage these per-site properties to identify promising perovskites with high theoretical OER activity. The identified design principles and promising materials provide a roadmap for closing the gap between current artificial catalysts and biological enzymes such as photosystem II.

Chemistry↗
DOE OSTI2020

Confirmation of the PPLB Derivative Discontinuity: Exact Chemical Potential at Finite Temperatures of a Model System

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semilocal approximations. Here, we provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

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

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations

Here we extend our recently developed sparse-stochastic fragmented exchange formalism for ground-state near-gap hybrid DFT to calculate absorption spectra within linear-response time-dependent generalized Kohn-Sham DFT (LR-GKS-TDDFT) for systems consisting of thousands of valence electrons within a grid-based/plane-wave representation. A mixed deterministic/fragmented-stochastic compression of the exchange kernel, here using long-range explicit exchange functionals, provides an efficient method for accurate optical spectra. Both real-time propagation as well as frequency-resolved Casida-equation-type approaches for spectra are presented, and the method is applied to large molecular dyes.

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

Range-separated hybrid functionals for mixed dimensional heterojunctions: Application to phthalocyanines/MoS 2

We analyze the electronic structure and level alignment of transition-metal phthalocyanine (MPc) molecules adsorbed on two-dimensional MoS 2 employing density functional theory (DFT) calculations. We develop a procedure for multi-objective optimal tuning of parameters of range-separated hybrid functionals in these mixed-dimensional systems. Using this procedure, which leads to the asymptotically-correct exchange-correlation potential between molecule and two-dimensional material, we obtain electronic structures consistent with experimental photoemission results for both energy level alignment and electronic bandgaps, representing a significant advance compared to standard DFT methods. We elucidate the MoS 2 valence resonance with the transition-metal phthalocyanine non-frontier 3d orbitals and its dependence on the transition metal atomic number. Based on our calculations, we derive parameter-free, model self-energy corrections that quantitatively accounts for the effects of the heterogeneous dielectric environment on the electronic structure of these mixed-dimensional heterojunctions.

2D materials↗
DOE OSTI2025

GPU acceleration of hybrid functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing hybrid functional calculations in generalized Kohn–Sham density functional theory. In particular, we develop a batch variant of the recently formulated Kronecker product-based linear solver for the simultaneous solution of multiple linear systems. We then develop a modular, math kernel based implementation for hybrid functionals on NVIDIA architectures, where computationally intensive operations are offloaded to the GPUs, while the remaining workload is handled by the central processing units (CPUs). Considering bulk and slab examples, we demonstrate that GPUs enable up to 8× speedup in node-hours and 80× in core-hours compared to CPU-only execution, reducing the time to solution on V100 GPUs to around 300 s for a metallic system with over 6000 electrons, and significantly reducing the computational resources required for a given wall time.

Kohn-Sham density functional theory↗
DOE OSTI2024

ReaxFF Reactive Force Field for Exploring Electronically Switchable Polarization in Zn 1– x Mg x O Ferroelectric Semiconductors

Cation misfit in traditional ferroelectric crystals offers a new material platform that can drive electronic components toward structural miniaturization and high-density integration, enabling deviation from von-Neumann architectures. Here, we explore ferroelectricity in Zn 1–x Mg x O, a nontraditional ferroelectric material with tunable properties. Using data from density-functional theory calculations, we have developed a ReaxFF reactive force field to explore the ferroelectric properties of Zn 1–x Mg x O and reveal the hysteresis behavior. We discover that ferroelectric switching can be observed at a critical thickness of 10 nm with a residual polarization of ~100 μC/cm 2 . Our analysis indicates that an increase in Mg-substitution correlates with a decrease in the coercive field. We also observe a strong temperature dependence of the coercive field in Zn 1–x Mg x O, with values decreasing as the temperature increases. Additionally, we find that the distribution of Mg atoms significantly impacts the coercive field, with a clustered distribution leading to a substantial increase. In particular, a decrease in coercive field values is observed when Mg atoms are randomly distributed, compared to uniform distribution. Furthermore, leveraging tunable hysteresis behavior offered by varied percentages and distribution of Mg-substitution provides valuable insights into the design of next-generation functional devices and will inspire further investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

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

Electron–electrophile coupled dinitrogen reduction in a cerium– meta -tetraphenolate system: a computational study

The use of lanthanide complexes for catalytic dinitrogen reduction is a new development in homogeneous catalysis. Density functional theory calculations on our recently reported cerium phenolate catalyst [K 2 Ce 2 (sol) 4 (mTP) 2 ] (mTP = {(OC 6 H 2 - 2- t Bu-4-Me) 2 CH} 2 - 1,3-C 6 H 4 ; sol = OMe 2 here; THF in the experiment) have been undertaken to elucidate the reduction, activation and silylation steps at the bound dinitrogen molecule, in the presence of the reductant, potassium metal (K 0 ) and the electrophile Me 3 SiCl (TMSCl). Out of the total of six electron reductions required to cleave the N 2 , the first two-electron reduction step was found to be highly disfavoured unless potassium cations (K + ) are included, upon which the step is rendered strongly exergonic; N–Si bond formation at the two-electron stage is predicted to be unfavourable. The three-electron-reduced N 2 -adduct is found to be at the reductive activation limit in the absence of added electrophiles, which can form N-element bonds and lower the overall charge. Added electron density beyond three-electron reduction no longer localises on N 2 , preventing formal N 2 4− formation. A pathway in which both K 0 and Me 3 SiCl work in concert was modelled, and six sequential reduction–silylation steps were calculated, showing how the N–N bond is cleaved after the third reduction, eventually releasing two equivalents of N(SiMe 3 ) 3 , and regenerating the starting complex with the highest barrier of any step being 22 kcal mol −1 . We establish alkali metal coordination and coupled electron–electrophile transfer as key factors in the design of rare-earth-mediated dinitrogen functionalisation.

Ahmad, Shahbaz [Univ. of Manchester (United Kingdo↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional

The absence of a reliable formulation of kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions which have been proposed in the literature and the appropriate coefficients are calculated via a linear regression framework. To model the nonlocal contributions, we explore two new nonlocal functionals - a functional that captures fluctuations in electronic density and a functional that incorporates gradient information. Since, the analytical functional forms of the kernels present in these nonlocal terms are not known from theory, we propose a basis function expansion to model these seemingly difficult nonlocal quantities. This allows us to easily reconstruct kernels for any system using only a few structures. The proposed method is able to learn kinetic energy densities and total kinetic energies of molecular and periodic systems, such as H 2 , LiH, LiF and a one-dimensional chain of 8 hydrogens using data from Kohn-Sham density functional theory calculations for only a few structures. For the ease of reproduction, codes used to generate the models are provided in the supporting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A reaction mechanism for plasma electrolysis of AgNO 3 forming silver nanoclusters and nanoparticles

In plasma-driven solution electrolysis (PDSE), gas-phase plasma-produced species interact with an electrolytic solution to produce, for example, nanoparticles. An atmospheric pressure plasma jet (APPJ) directed onto a liquid solution containing a metallic salt will promote reduction of metallic ions in solution, generating metallic clusters that nucleate to form nanoparticles. In this article, results from a computational investigation are discussed of a PDSE process in which a radio-frequency APPJ sustained in helium impinges on a silver nitrate solution, resulting in growth of silver nanoparticles. A reaction mechanism was developed and implemented in a global plasma chemistry model to predict nanoparticle growth. To develop the reaction mechanism, density functional theory was used to generate probable silver growth pathways up to Ag 9 . Neutral clusters larger than Ag 9 were classified as nanoparticles. Kinetic reaction rate coefficients for thermodynamically favorable growth pathways were estimated based on an existing, empirically determined base reaction mechanism for smaller Ag particle interactions. These rates were used in conjunction with diffusion-controlled reaction rate coefficients that were calculated for other Ag species. The role of anions in reduction of Agn ions in forming nanoparticles is also discussed. Oxygen containing impurities or admixtures to the helium, air entrainment into the APPJ, and dissociation of saturated water vapor above the solution can produce additional reactive oxygen species in solution, resulting in the production of anions and O 2 – in particular. For a given molarity, delivering a sufficient fluence of reducing species will produce similar nanoparticle densities and sizes for all applied power levels. In conclusion, comparisons are made to alternate models for nanoparticle formation, including charged nanoparticles and use of direct current plasmas.

36 MATERIALS SCIENCE↗