Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculation”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗
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 OSTI2018

Nonequivalent Spin Exchanges of the Hexagonal Spin Lattice Affecting the Low-Temperature Magnetic Properties of RInO 3 (R = Gd, Tb, Dy): Importance of Spin–Orbit Coupling for Spin Exchanges between Rare-Earth Cations with Nonzero Orbital Moments

Here, rare-earth indium oxides RInO 3 (R = Gd, Tb, Dy) consist of spin-frustrated hexagonal spin lattices made up of rare-earth ions R 3+ , where R 3+ = Gd 3+ (f 7 , L = 0), Tb 3+ (f 8 , L = 3), and Dy 3+ (f 9 , L = 5). We carried out DFT calculations for RInO 3 , including on-site repulsion U with/without spin–orbit coupling (SOC), to explore if their low-temperature magnetic properties are related to the two nonequivalent nearest-neighbor (NN) spin exchanges of their hexagonal spin lattices. Our DFT + U + SOC calculations predict that the orbital moments of the Tb 3+ and Dy 3+ ions are smaller than their free-ion values by ~2μ B while the Tb 3+ spins have an inplane magnetic anisotropy, in agreement with the experiments. This suggests that the f orbitals of each R 3+ (R = Tb, Dy) ion are engaged, though weakly, in bonding with the surrounding ligand atoms. The magnetic properties of GdInO 3 with the zero orbital moment are adequately described by the spin exchanges extracted by DFT + U calculations. In describing the magnetic properties of TbInO 3 and DyInO 3 with nonzero orbital moments, however, the spin exchanges extracted by DFT + U + SOC calculations are necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

First-principles study of crystal structure and stability of T 1 precipitates in Al-Li-Cu alloys

Aluminum-lithium-copper alloys have a low density, high elastic modulus and high specific strength.Due to this combination of properties, alloys strengthened with the ternary (Al-Li-Cu) T 1 phase have attracted a great deal of interest especially in aerospace applications. Determining the atomic structuralinformation of the precipitate is a fundamental step in developing a basis for advanced alloy design; however, even though many experimental studies have addressed the T 1 crystal structure, it remains the subject of some controversy. Here, we use density functional theory (DFT) calculations to investigate the structure and composition of the T 1 phase by comparing the energetic stability of five previously-proposed models of the crystal structure of T 1 . The DFT formation energy of these proposed T 1 crystal structures was calculated using a special quasi-random structure (SQS) approach to describe a disordered Al-Cu sub-lattice. In conflict with the experimental phase diagram, DFT calculations of all five proposed models result in an energetically unstable T 1 phase. Here, we search for a new, lower-energy structure of T 1 using a cluster expansion approach, and find a new structural model with DFT energy that is stable (at T = 0 K), i.e., on the calculated convex hull of the Al-Li-Cu ternary system. However, this new predicted phase does not have a tie-line with Al, but the formation energy of the phase is very close to the energy required to make a tie-line with Al (Δ E = 0.013 eV/atom), which could be affected by finite temperature entropic effects (i.e., vibrational entropic stabilization).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Orbital-Free Quantum Simulation Methods for Application to Warm Dense Matter (Final Technical Report)

Predictive simulations for prediction of condensed system behavior in state conditions far from ambient is increasingly crucial to DOE priorities. Warm dense matter (WDM) is the paradigm: temperature T > 1-15 eV, pressures P to 1 Mbar or greater. Experiments under such state conditions are difficult and costly. We summarize work driven by the need and opportunity to make free-energy density functional theory (DFT) as powerful a tool for ab initio simulation of matter under such extreme conditions as ground state DFT is for ordinary matter Advancing orbital-free DFT (OF-DFT) to eliminate the Kohn-Sham (KS) scaling bottleneck in such simulations is the other priority. The concurrent challenge for both goals is the intrinsic complexity of WDM. We summarize 15 years of successes and major progress on (1) free energy exchange-correlation functionals; (2) non-interacting free energy functionals (counterpart to T=0 Kohn-Sham kinetic energy density functionals); (3) rigorous results and constraints for free-energy DFT; (4) software for free energy DFT calculations in both conventional Kohn-Sham and OF-DFT form; (5) de-orbitalization of advanced orbital-dependent ground state functionals for use in OF-DFT; (6) demonstration calculations; (7) ancillary achievements (e.g. major review articles, secondary explorations motivated by primary goals).

36 MATERIALS SCIENCE↗
DOE OSTI2018

Charge-patching method for the calculation of electronic structure of polypeptides

Theoretical study of the electronic structures of protein is a fundamental challenge in computational biochemistry due to the large size of the systems. The electronic structure of a protein is important for some of the important protein functionalities, such as photosynthesis. Here, we explored the charge-patching method to calculate the electronic structure of polypeptides. This method generates the charge densities of the systems by patching the charge motifs calculated from small prototype systems. The method was tested on a range of polypeptides, including the glycine polypeptide in 2 7 -ribbon, α-helix, 3 10 -helix, and β-strand structures. After the charge density profiles of these systems were obtained, the electronic structures of these glycine polypeptides were further calculated based on density functional theory (DFT) using a folded-spectrum method. The highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) were analyzed and compared with conventional direct DFT calculations. The charge-patching method results were found to be in good agreement with the directed DFT results.

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

Development of an interatomic potential for the W–Ta system

A physics-inspired, data-driven interatomic potential framework for multi-element system is presented. This potential is based on the generalization of the embedded atom method potential and systematically incorporates two-, three- and many-body effects. Different atomic environments are described by atom-centered Gaussian basis sets which provides sufficient flexibility to capture diverse atomic environments and allows for easy optimization of the free parameters. An interatomic potential model for the tungsten–tantalum (W–Ta) system is developed using this framework by training on data from ab initio density functional theory (DFT) calculations. The potential is thoroughly tested at various compositions by comparing bulk and defect properties from experimental and DFT data. It is shown that the potential predicts the elastic constants, defect properties, such as vacancy and interstitial formation energies, core structures of dislocations, and melting points for both pure and tantalum and tungsten with an accuracy comparable to other single element machine learning based potentials. In conclusion, the potential model is used to investigate the formation energies of ordered alloys and vacancy formation energies for chemically random W–Ta alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

In Situ Studies of Methanol Decomposition Over Cu(111) and Cu2O/Cu(111): Effects of Reactant Pressure, Surface Morphology, and Hot Spots of Active Sites

The dissociative adsorption of methanol was investigated on Cu(111) and ultrathin Cu 2 O films. We employed synchrotron-based Ambient Pressure X-ray Photoelectron Spectroscopy (AP-XPS) and Scanning Tunneling Microscopy (STM) to study the dynamics of gas–solid interactions, and calculations based on Density Functional Theory (DFT) were used to examine the reaction path. C 1s XPS spectra revealed that methanol underwent dissociative adsorption on plain Cu(111) to form methoxy (CH 3 O), formaldehyde (H 2 CO), and formate (HCOO) at a pressure range of 0.5–10 mTorr, with these species remaining on the surface after evacuation. This was accompanied by the appearance of a low coverage (~0.05 ML) of O ads in the O 1s which can be considered a highly active site for methanol activation. The high activity is apparent by a coverage of 0.8 ML of methoxy at room temperature. STM was unable to image these species at room temperature as they were highly mobile on metallic copper. In contrast, for CH 3 OH on Cu 2 O/Cu(111), STM showed clear hot spots for reaction and a complex array of adsorption structures. On the oxide substrate, there was decomposition of methanol to H 2 CO, CH 3 O, HCOO, and hydrocarbon species (CH x ) due to the subsequent interactions of methanol with lattice oxygen. Cu(111) remained entirely saturated with decomposition products under 10 mTorr of methanol (θ ≈ 0.97 ML), whereas the Cu 2 O overlayer was saturated at a much lower coverage (θ ≈ 0.30 ML). STM revealed rows and step edges of Cu 2 O decorated with decomposition products and metallic Cu islands ~5 nm in size. The difference in activity between Cu(111) and Cu 2 O/Cu(111) is attributed to the significant amount of O present on the oxide surface. Additionally, Density Functional theory (DFT) calculations described the XPS measurements well, showing a likely methanol dissociation to *CH 3 O and therefore a surface reduction. More importantly, the DFT results revealed that it was the chemisorbed oxygen on Cu 2 O/Cu(111) which oxidized the dissociated *CH 3 O to *HCOO and eventually CO 2 , while the reaction only involving upper oxygen on the Cu 2 O hexagonal ring led to the formation of H 2 CO.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting magnetic properties of single-molecule magnets from self-interaction-free density-functional theory (Final Report)

In this project we investigated electronic structure of an intermediate-sized copper-based molecule and magnetic and hyperfine properties of several small non-magnetic and magnetic molecules including transition-metal elements by applying self-interaction corrections to density functional theory (DFT). In a sufficient number of cases, DFT-calculated electronic structure and magnetic properties of single-molecule magnets qualitatively differ from corresponding experimental data. This is partly due to self-interacting electrons within the DFT formalism. Recently, an efficient method to correct the self-interactions was proposed, i.e., Fermi-Lowdin prbital (FLO) based self-interaction corrected (SIC) methodology, within DFT. Henceforth, this method is referred to as FLO-SIC method which exists in FLOSIC code. We used this FLO-SIC method for our studies of electronic structure and magnetic and hyperfine properties of small magnetic molecules and non-magnetic molecules. Our study will provide insight into predictions of magnetic properties of single-molecule magnets where self-interaction corrections play a critical role. There are two components of this project. In the first work, we studied the electronic structure of a planar mononuclear Cu-based molecule in two oxidation states, using DFT with the FLO-SIC method. We chose this system because it is small enough and it includes a transition metal element. We found that the standard FLO-SIC method takes too much compute time even for the small transition-metal molecule and so we slightly modified the method in order to expedite the process. In the dianionic state, we found that the FLO-SIC spin density agrees quantitatively with accurate quantum chemistry methods, while DFT spin density without self-interaction corrections are severely deviated from the quantum chemistry methods. We also showed that the energy gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the dianionic state is larger than that of the monoanionic state. This result is consistent with experimental data. In the second work, we investigated how the interaction between the electron spin and nuclear spin is affected by electron self-interactions within small non-magnetic and magnetic molecules. Such an interaction is called hyperfine interaction. For molecules without significant orbital angular momentum, the hyperfine interaction consists of Fermi contact and dipolar interaction terms. Since the Fermi contact term depends on electron spin density at the nuclear site, it would be highly affected by self-interaction corrections. Therefore, we calculated the hyperfine interaction for the small molecules using DFT with the slightly modified expedited FLO-SIC method which was obtained in the first work, and compared the results to experimental data and DFT calculations without self-interaction corrections. We found significant improvement of the Fermi contact term computed using the FLO-SIC method for small magnetic molecules. Overall, the first and second work provided positive outlook of application of the FLO-SIC method to magnetic molecules and systems including transition-metal elements.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Defect energetics of cubic hafnia from quantum Monte Carlo simulations

Cubic hafnia (HfO 2 ) is of great interest for a number of applications in electronics because of its high dielectric constant. However, common defects in such applications degrade the properties of hafina. We have investigated the electronic properties of oxygen vacancies and nitrogen substitution in cubic HfO 2 using first-principles calculations based on density functional theory (DFT) and many-body diffusion Monte Carlo (DMC) methods. We investigate five different charge defect states of oxygen vacancies, as well as substitutional N defects which can lead to local magnetic moments. Both DMC and DFT calculations shows that an oxygen vacancy induces strong lattice relaxations around the defect. Finally, we compare defect formation energies, charge and spin densities obtained from DMC with results obtained using DFT. While the obtained formation energies from DMC are 0.6–1.5 eV larger than those from GGA+ U , the agreement for the most important defects, neutral and positively charged oxygen vacancies, and nitrogen substitutional defect, under oxygen-poor conditions are in reasonably good agreement. Lastly, our work confirms that nitrogen can act to passivate cubic hafnia for applications in electronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Point defects in CdTe and CdSeTe alloy: a first principles investigation with DFT+U

CdTe and its alloy CdSeTe are widely used in optoelectronic devices, such as radiation detectors and solar cells, due to their superior electrical properties. However, the formation of defects and defect complexes in these materials can significantly affect their performance. As a result, understanding the defect formation and recombination processes in CdTe and CdSeTe alloy is of great importance. In recent years, density functional theory (DFT) calculations have emerged as a powerful tool for investigating the properties of defects in semiconductors. In this paper, we use DFT+U calculations to comprehensively study the properties of intrinsic defects as well as extrinsic defects induced by commonly used dopants, such as Cu and group V elements, in CdTe and CdSeTe alloy. This work provides insights into the effects of these defects on the electrical and optical properties of the material.

14 SOLAR ENERGY↗
DOE OSTI2021

Theoretical and experimental studies of compression and shear deformation behavior of Osmium to 280 GPa

Abstract The compression behavior of osmium metal was investigated up to 280 GPa (volume compression V/V o = 0.725) under nonhydrostatic conditions at ambient temperature using angle dispersive axial x-ray diffraction (A-XRD) with a diamond anvil cell (DAC). In addition, shear strength of osmium was measured to 170 GPa using radial x-ray diffraction (R-XRD) technique in DAC. Both diffraction techniques in DAC employed platinum as an internal pressure standard. Density functional theory (DFT) calculations were also performed, and the computed lattice parameters and volumes under compression are in good agreement with the experiments. DFT predicts a monotonous increase in axial ratio ( c/a ) with pressure and the structural anomalies of less than 1% in ( c/a ) ratio reported below 150 GPa were not reproduced in theoretical calculations and hydrostatic measurements. The measured value of shear strength of osmium ( τ ) approaches a limiting value of 6 GPa above a pressure of 50 GPa in contrast to theoretical predictions of 24 GPa and is likely due to imperfections in polycrystalline samples. DFT calculations also enable the studies of shear and tensile deformations. The theoretical ideal shear stress is found along the (001)[1–10] shear direction with the maximal shear stress ∼24 GPa at critical strain ∼0.13.

Lin, Chia-Min↗
DOE OSTI2016

Development of Xe and Kr empirical potentials for CeO 2 , ThO 2 , UO 2 and PuO 2 , combining DFT with high temperature MD

In this study, the development of embedded atom method (EAM) many-body potentials for actinide oxides and associated mixed oxide (MOX) systems has motivated the development of a complementary parameter set for gas-actinide and gas-oxygen interactions. A comprehensive set of density functional theory (DFT) calculations were used to study Xe and Kr incorporation at a number of sites in CeO 2 , ThO 2 , UO 2 and PuO 2 . These structures were used to fit a potential, which was used to generate molecular dynamics (MD) configurations incorporating Xe and Kr at 300 K, 1500 K, 3000 K and 5000 K. Subsequent matching to the forces predicted by DFT for these MD configurations was used to refine the potential set. This fitting approach ensured weighted fitting to configurations that are thermodynamically significant over a broad temperature range, while avoiding computationally expensive DFT-MD calculations. The resultant gas potentials were validated against DFT trapping energies and are suitable for simulating combinations of Xe and Kr in solid solutions of CeO 2 , ThO 2 , UO 2 and PuO 2 , providing a powerful tool for the atomistic simulation of conventional nuclear reactor fuel UO 2 as well as advanced MOX fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Reflections on one million compounds in the open quantum materials database (OQMD)

Abstract Density functional theory (DFT) has been widely applied in modern materials discovery and many materials databases, including the open quantum materials database (OQMD), contain large collections of calculated DFT properties of experimentally known crystal structures and hypothetical predicted compounds. Since the beginning of the OQMD in late 2010, over one million compounds have now been calculated and stored in the database, which is constantly used by worldwide researchers in advancing materials studies. The growth of the OQMD depends on project-based high-throughput DFT calculations, including structure-based projects, property-based projects, and most recently, machine-learning-based projects. Another major goal of the OQMD is to ensure the openness of its materials data to the public and the OQMD developers are constantly working with other materials databases to reach a universal querying protocol in support of the FAIR data principles.

08 HYDROGEN↗
DOE OSTI2021

HPC for Fission Modeling in Support of Nuclear Data [Slides]

The presentation discusses the background of the High-Performance Computing (HPC) division at the Los Alamos National Laboratory (LANL), with emphasis on its computational needs and computational techniques. It also discusses the highlights of the HPC and the results of calculations, including energy sharing parameterized from DFT calculations, parameterization used in CGMF calculations, angular momentum from microscopic calculations, and correlations between FF angular momenta. The presentation concludes with an outlook including what remains to be done in the HPC division and how it will benefit nuclear data, specifically how it will provide a better understanding of the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Structural and Electrochemical Characteristics of Ca-Doped “Flower-like” Li 4 Ti 5 O 12 Motifs as High-Rate Anode Materials for Lithium-Ion Batteries

Doped motifs offer an intriguing structural pathway toward improving conductivity for battery applications. Specifically, Ca-doped, three-dimensional “flower-like” Li 4–x Ca x Ti 5 O 12 (“x” = 0, 0.1, 0.15, and 0.2) micrometer-scale spheres have been successfully prepared for the first time using a simple and reproducible hydrothermal reaction followed by a short calcination process. The products were experimentally characterized by means of X-ray diffraction (XRD), transmission electron microscopy (TEM), scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS) mapping, inductively coupled plasma optical emission spectrometry (ICP-OES), X-ray photoelectron spectroscopy (XPS), cyclic voltammetry (CV), electrochemical impedance spectroscopy (EIS), and galvanostatic charge–discharge testing. Calcium dopant ions were shown to be uniformly distributed within the LTO structure without altering the underlying “flower-like” morphology. The largest lattice expansion and the highest Ti 3+ ratios were noted with XRD and XPS, respectively, whereas increased charge transfer conductivity and decreased Li + -ion diffusion coefficients were displayed in EIS for the Li 4–x Ca x Ti 5 O 12 (“x” = 0.2) sample. The “x” = 0.2 sample yielded a higher rate capability, an excellent reversibility, and a superior cycling stability, delivering 151 and 143 mAh/g under discharge rates of 20C and 40C at cycles 60 and 70, respectively. In addition, a high cycling stability was demonstrated with a capacity retention of 92% after 300 cycles at a very high discharge rate of 20C. In addition, first-principles calculations based on density functional theory (DFT) were conducted with the goal of further elucidating and understanding the nature of the doping mechanism in this study. The DFT calculations not only determined the structure of the Ca-doped Li 4 Ti 5 O 12 , which was found to be in accordance with the experimentally measured XPD pattern, but also yielded valuable insights into the doping-induced effect on both the atomic and electronic structures of Li 4 Ti 5 O 12 .

25 ENERGY STORAGE↗
DOE OSTI2024

Reactivity of Sulfur Vacancy-Rich MoS 2 to Water Dissociation

Exposed Mo atoms on the surface of MoS 2 can catalyze certain useful chemical reactions, including the dissociation of water to produce hydrogen. However, a comprehensive understanding of water reactivity on defective MoS 2 surfaces remains elusive. Here, we use in situ near ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) to investigate water dissociation reactions on MoS 2 surfaces before and after Ar + ion beam bombardment. To make the surfaces reactive, we treated them with Ar + ion beam sputtering first, which created exposed Mo sites. Using ultraviolet photoelectron spectroscopy (UPS) and density of states calculations conducted using density functional theory (DFT), we verified that stripping the surface of sulfur atoms creates metallic surface states that can catalyze water dissociation. At elevated H 2 O pressures, XPS measurements combined with DFT calculations suggested the presence of four distinct surface species from dissociated water. Specifically, we found that oxides and hydroxides are prominent at the surface, while chemisorbed and physisorbed H 2 O molecules are also present. In conclusion, this study provides new insights that reveal the prospects of surface-engineered MoS 2 as a catalyst for water dissociation.

08 HYDROGEN↗