Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Re-examining the silicon self-interstitial charge states and defect levels: A density functional theory and bounds analysis study

The self-interstitial atom (SIA) is one of two fundamental point defects in bulk Si. Isolated Si SIAs are extremely difficult to observe experimentally. Even at very low temperatures, they anneal before typical experiments can be performed. Given the challenges associated with experimental characterization, accurate theoretical calculations provide valuable information necessary to elucidate the properties of these defects. Previous studies have applied Kohn–Sham density functional theory (DFT) to the Si SIA, using either the local density approximation or the generalized gradient approximation to the exchange-correlation (XC) energy. The consensus of these studies indicates that a Si SIA may exist in five charge states ranging from -2 to +2 with the defect structure being dependent on the charge state. This study aims to re-examine the existence of these charge states in light of recently derived “approximate bounds” on the defect levels obtained from finite-size supercell calculations and new DFT calculations using both semi-local and hybrid XC approximations. We conclude that only the neutral and +2 charge states are directly supported by DFT as localized charge states of the Si SIA. Within the current accuracy of DFT, our results indicate that the +1 charge state likely consists of an electron in a conduction-band-like state that is coulombically bound to a +2 SIA. Furthermore, the -1 and -2 charge states likely consist of a neutral SIA with one and two additional electrons in the conduction band, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Electronic Absorption Spectra of Neutral Perylene (C20H12), Terrylene (C30H16), and Quaterrylene (C40H20) and their Positive and Negative Ions: Ne Matrix-Isolation Spectroscopy and Time Dependent Density Functional Theory Calculations

We present a full experimental and theoretical study of an interesting series of polycyclic aromatic hydrocarbons, the oligorylenes. The absorption spectra of perylene, terrylene and quaterrylene in neutral, cationic and anionic charge states are obtained by matrix-isolation spectroscopy in Ne. The experimental spectra are dominated by a bright state that red shifts with growing molecular size. Excitation energies and state symmetry assignments are supported by calculations using time dependent density functional theory methods. These calculations also provide new insight into the observed trends in oscillator strength and excitation energy for the bright states: the oscillator strength per unit mass of carbon increases along the series.

Halasinski, Thomas M.↗
DOE OSTI2023

Density Functional Theory Guided Investigation of Ligand‐Induced Neptunyl‐Neptunyl Interactions

Abstract Actinyl‐actinyl interactions are particularly prevalent for the pentavalent neptunyl cation (Np(V)O 2 ) + where these interactions appear either as a T ‐ or D ‐shape (diamond‐shape). T ‐shaped interactions have been previously identified in high concentration Np(V) solutions containing simple anions (NO 3 − , ClO 4 − , Cl − ) whereas D ‐shaped have only been isolated in the solid‐state in the presence of carboxylate ligands. In this study, Density Functional Theory (DFT) calculations were paired with Raman spectroscopy to evaluate the formation of D ‐shaped interactions in the presence of aliphatic (R=H (formate), CH 3 (acetate), CH 2 CH 3 (propionate)) and aromatic (R=C 6 H 5 (benzoate), C 6 H 4 OH (4‐hydroxybenzoate), C 5 H 4 N (isonicotinate)) carboxylate ligands. DFT studies indicate that the Δ G to form hydrated T ‐ and D ‐shaped forms are not spontaneous but become so with the addition of the carboxylate ligands. Raman spectra of the Np(V) carboxylate solutions contained vibrational modes associated with the D ‐shaped interactions, but spectral changes observed over time indicate a dynamic system. Crystallization experiments from the Np(V) carboxylate systems confirmed the presence of D‐ shaped dimers for the aromatic carboxylates, suggesting that the choice of the anion in solution favors actinyl‐actinyl interactions even at low concentrations (≤20 mM) of Np(V).

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

Broken-Symmetry Density Functional Theory Analysis of the Ω Intermediate in Radical S -Adenosyl- l -methionine Enzymes: Evidence for a Near-Attack Conformer over an Organometallic Species

Radical S-adenosyl-L-methionine (SAM) enzymes are found in all domains of life and catalyze a wide range of biochemical reactions. Recently, an organometallic intermediate, Ω, has been experimentally implicated in the 5'-deoxyadenosyl radical generation mechanism of the radical SAM superfamily. In this work, we employ broken-symmetry density functional theory to evaluate several structural models of Ω. The results show that the calculated hyperfine coupling constants (HFCCs) for the proposed organometallic structure of Ω are inconsistent with the experiment. In contrast, a near-attack conformer of SAM bound to the catalytic [4Fe–4S] cluster, in which the distance between the unique iron and SAM sulfur is ~3 Å, yields HFCCs that are all within 1 MHz of the experimental values. Furthermore, these results clarify the structure of the ubiquitous Ω intermediate and suggest a paradigm shift reversal regarding the mechanism of SAM cleavage by members of the radical SAM superfamily.

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

Microscopic force driving the photoinduced ultrafast phase transition: Time-dependent density functional theory simulations of IrTe 2

Photoinduced phase transitions can have complex and intriguing behaviors more than material ground-state dynamics. Understanding the underlying mechanism can help us to design new ways to manipulate the materials. A variety of mechanisms has been proposed to explain the photoinduced phase transitions of IrTe 2 , but a consensus has yet to be reached. Here, we study the photo-induced phase transitions of IrTe 2 by performing the real-time time-dependent density functional theory (rt-TDDFT) simulations in combination with the occupation constrained DFT method. We reveal that the microscopic driving force for the photo-induced phase transitions arises from the tendency to lower the energy levels occupied by the excited carriers, through the increase or decrease of the associated atomic pair distances, depending on whether the newly occupied states are antibonding or bonding states, respectively. The geometric constraints between different bonds represented by the Poisson ratio can bring together different tendencies from different atomic pairs, thus forming a complex intriguing dynamic picture depending on the intensity of the excitation. Furthermore, we also find that phonons don't play a primary role, but can assist the phase transition. These findings imply that one can control the structural phase transitions by selectively exciting photocarriers into designated atomic states using appropriate photon sources.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

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 OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

Polaron-induced metal-to-insulator transition in vanadium oxides from density functional theory calculations

Vanadium oxides have been extensively studied as phase-change memory units in artificial synapses for neuromorphic computing due to their metal-insulator transitions (MIT) at or near room temperature. Recently, injection of charge carriers into vanadium oxides, e.g., via optically via a heterostructure, has been proposed as an alternative switching mechanism and also potentially as a means to tune the MIT temperature. In this study, we explore the formation of small polarons in the low temperature (LT) insulating phases for V 3 O 5 ,VO 2 , and V 2 O 3 , and the barriers to their migration using density functional theory calculations. We find that V 3 O 5 exhibits very low hole and electron polaron migration barriers (<100 meV) compared to V 2 O 3 and VO 2 , leading to much higher estimated polaronic conductivity. We also link the relative migration barriers to the amount of distortion that has to travel when the polaron migrate from one site to another. Polarons in V 3 O 5 also have smaller binding energies to vanadium and oxygen vacancy defects. Furthermore, these results suggest that the triggering of the MIT via injection of charge carriers are due to the formation of small polarons that can migrate rapidly through the crystal.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Revealing symmetry-broken superconducting configurations by density functional theory

A coherent theory for the superconductivity of both conventional and unconventional superconductors is currently lacking. Here we show that superconductivity arises from the formation of a symmetry-broken superconducting configuration (SCC) due to atomic perturbation of the normal conducting configuration (NCC). This electron–phonon interaction creates straight one-dimensional tunnels (SODTs) for charge density of electrons and/or holes as revealed by the calculations based on density functional theory (DFT). The SODTs act as resistance-free superhighways and are correlated to the Cooper pairs in the Bardeen–Cooper–Schrieffer (BCS) theory. The formation of SODTs implies that the electron–phonon interaction in the BCS theory can be represented by the difference in charge densities between SCC and NCC predicted by DFT. The present work highlights that in conventional superconductors, SODTs are embedded within the bulk materials and are easily destroyed by phonon vibrations, resulting in a low critical superconducting temperature (T C ). Conversely, in unconventional superconductors such as YBa 2 Cu 3 O 7 (YBCO 7 ), SODTs are protected by a layered pontoon structure with very weak bonding to the bulk materials, maintaining SODTs’ stability at higher temperatures and leading to a much higher T C . The present approach is validated for 14 conventional superconductors of 18 pure elements and MgB 2 examined in this work, including the presently predicted superconductivity in Cu, Ag, Au, Sb, Bi, and MgB 2 at 0 K and 0 GPa, and one unconventional superconductor of YBCO 7 . Our discovery indicates that DFT can be a practical tool for predicting superconductors, enabling a systematic search for new superconducting materials in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coupling computational thermodynamics with density-function-theory based calculations to design L12 precipitates in Fe Ni based alloys

Achieving a high-volume fraction of thermodynamically stable L12-type precipitates that are resistant to coarsening is of great importance for the development of low-cost FeNi based austenitic steels. With the aid of computational thermodynamics, this work designed two model alloys: Fe-37.4Ni-6.1Al-2.9Ti (FNAT) and Fe-45.2Ni-5.9Al-8.5Si (FNAS). Both alloys were designed to contain a similar amount of L1 2 precipitate in Fe-Ni matrix without forming other precipitates. Density-Function-Theory (DFT) calculation was coupled with computational thermodynamics to predict the critical radius at which the precipitates change shape from spherical to cuboidal. The calculation results suggest that critical radius for the FNAT alloy is about two orders of magnitude larger than that for the FeNiAlSi alloy. Phase stability and morphology of the L1 2 precipitates in these two alloys were experimentally investigated through X-ray diffraction, atom probe tomography, and scanning and transmission electron microscopy. The L1 2 precipitates in the FeNiAlSi system were found to be cuboidal and rod shaped, with much larger size than the spherical ones in the FeNiAlTi system, agreeing with the calculation results. This work suggested that coupling computational thermodynamics with DFT calculations can be reliably used to design L1 2 precipitates in FeNi based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Characterization of Lanthanum Monazite Surface Chemistry and Crystal Morphology through Density Functional Theory and Experimental Approaches

Monazite is a rare earth element (REE)-containing mineral that consists of (REE)PO 4 formal units and is one of the most important sources of these critical materials. The concentration of REEs from mined monazite ore often involves froth flotation, which is a beneficiation process that enhances the efficiency of downstream processing. The effectiveness of froth flotation is largely governed by the ability of collector agents to selectively bind to monazite particles. Thus, a molecular-level understanding of monazite interfacial chemistry is integral to the design of effective collector agents. To address this need, we performed density functional theory (DFT) calculations and a variety of experimental techniques to characterize La-monazite and elucidate its crystal morphology. Interestingly, we find minimal differences in the predicted morphologies of La-monazite for hydrous and anhydrous environments, which are largely dominated by low-index facets (e.g., {110}, {100}, and {010}). Indexing of synthesized La-monazite crystals via X-ray diffraction also uncovers {110} and {100} as the predominant facets. The average surface energies of 0% and 100% water coverage La-monazite crystals were predicted to be 0.87 and 0.76 J/m 2 , respectively, while calorimetry suggests values of 1.30 and 1.15 J/m 2 , respectively. The apparent discrepancies between the theoretical and experimental values are expected and attributed to defects present in physical crystals, in contrast to the perfect mineral surfaces in simulations. The difference in surface energy between the 0% and 100% water coverage morphologies predicted by theory is consistent with the value measured via calorimetry. DFT reveals a wide range of adsorption energies for water across the studied facets, but in all cases, water is predicted to strongly bind to monazite surfaces with an average adsorption energy of -92.7 kJ/mol for a La-monazite single crystal. Finally, this study provides the groundwork necessary for the rational design of froth flotation collector agents by granting molecular-level insight into the predominant facets of monazite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗