Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density Functional calculations”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A combined machine learning and density functional theory study of binary Ti-Nb and Ti-Zr alloys: Stability and Young’s modulus

The multicomponent Ti alloys, specifically the -phase, have experienced a strong growth over the last decades, due to their outstanding properties of ultra-high strength and low Young’s modulus. These properties play a significant role in many aerospace and biomedical applications. Selection and optimization of multicomponent alloys is challenging due to the vast chemical and compositional space. Here we investigate the use of machine learning techniques informed by density functional calculations to guide the selection of Nb- and Zr-based Ti binary alloys. From the cubic structures obtained from high throughput calculations and literature, we identify several structures with Young’s moduli below 40 GPa. The multivariant decision tree methods provide efficient surrogate models to identify structure variables have high influences on the energetic stability and Young’s modulus. Furthermore, we implement a workflow of incorporating DFT provided results and machine learning method to explore the chemical and composition space of other binary and multicomponent alloys, to eventually accelerate the material design via taking advantages of identified key variables.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Isotope effect and critical magnetic fields of superconducting YH 6 : A Migdal-Eliashberg theory approach

The emergence of near-ambient temperature superconductivity under pressure in metal hydride systems has motivated a desire to further understand such remarkable properties, specifically critical magnetic fields. Here, YH 6 is suggested to be a departure from conventional superconductivity, due to apparent anomalous behavior. Using density functional calculations in conjunction with Migdal-Eliashberg theory we show that in YH 6 the critical temperature and the isotope effect under pressure, as well as the high critical fields, are consistent with strong-coupling conventional superconductivity, a property anticipated to extend to other related systems. Furthermore, strong-coupling corrections occur to the expected BCS values for the isotope effect coefficient (α), Ginzburg-Landau parameter [κ 1 ⁡(T)], London penetration depth [λ L (T)], electromagnetic coherence length [ξ⁡(T)], and the energy gap (Δ 0 ).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Magnetic interactions and excitations in SrMnSb 2

The magnetic interactions in the antiferromagnetic (AFM) Dirac semimetal candidate SrMnSb 2 are investigated using ab initio linear response theory and inelastic neutron scattering (INS). Our calculations reveal that the first two nearest in-plane couplings (J 1 and J 2 ) are both AFM in nature, indicating a significant degree of spin frustration, which aligns with experimental observations. The orbital resolution of exchange interactions shows that J 1 and J 2 are dominated by direct and superexchange, respectively. In a broader context, a rigid-band model suggests that electron doping fills the minority spin channel and results in a decrease in the AFM coupling strength for both J 1 and J 2 . To better compare with INS measurements, we calculate the spin-wave spectra within a linear spin-wave theory, utilizing the computed exchange parameters. Although the calculated spin-wave spectra somewhat overestimate the magnon bandwidth, they exhibit overall good agreement with measurements from INS experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Coexisting paramagnetic spins and long-range magnetic order in Ba4⁢(Ru0.92⁢Ir0.08)3⁢O10

We investigate the effect of dilute Ir substitution on the magnetism of the trimer-based ruthenate Ba4⁢Ru3⁢O10 using neutron diffraction, magnetic susceptibility measurements, first-principles calculations, and atomistic simulations. Neutron diffraction shows that Ir doping preserves the zigzag antiferromagnetic structure and the ordered-moment magnitude of the parent compound, in which the moments reside exclusively on the two outer Ru(2) sites of each Ru3⁢O12 trimer, while the central Ru(1) site remains nonmagnetic. The Néel temperature is reduced from ∼105 to 84.0(1) K upon 8% Ir substitution, while magnetic susceptibility reveals a pronounced low-temperature Curie-like upturn, indicating the coexistence of paramagnetic spins with long-range antiferromagnetic order. Density-functional calculations show that Ir preferentially occupies the central Ru(1) site, which disrupts the exchange pathways of the affected Ru3⁢O12 trimer state. Atomistic simulations incorporating this paramagnetic dilution reproduce the suppressed ordering temperature and the coexistence of ordered and paramagnetic components.

Islam, Farhan [Ames Laboratory]↗
DOE OSTI2024

Fine-tuning catalytic selectivity by modulating catalyst-environment interactions: CO 2 hydrogenation over Pd-based catalysts

Capturing catalytic behaviors under operational conditions is pivotal to gaining a mechanistic understanding and promoting the design of robust catalysts. The challenge lies in the difficulty of monitoring real-time surface dynamics driven by catalyst-environment interactions. Here, in this work, we introduce a framework based on density functional calculations and kinetic modeling. This framework significantly improves the accuracy of theoretical models’ descriptions of experimental observations by quantifying environmental impacts on surface phases and active sites. CO 2 hydrogenation over Pd-based catalysts is taken as a showcase. The observed selectivity variations of Pd and Pd-M bimetallic catalysts strongly correlate with hydrogen coverage maintained under typical CO 2 hydrogenation conditions. By reducing the amount of surface hydrogen, the selectivity tuned effectively from formic acid toward CO and methanol. This study not only deepens the comprehension of dynamics of active sites under active chemical conditions but also introduces an alternative opportunity for catalytic tuning by modulating catalyst-environment interactions.

03 NATURAL GAS↗
DOE OSTI2024

Self-compensation of group-V acceptors in CdTe

Cadmium Telluride is at the core of an important thin-film technology for photovoltaics that is already commercially available, yet the CdTe-based solar cell efficiency remains limited at 22%, well below the theoretical limit of ~30%. Increasing the hole concentration is crucial for higher efficiency, and group-V elements such as As, P, and Sb are potential dopants as they are shallow acceptors. Nevertheless, group-V doped p-type CdTe often exhibits low doping activation, and the compensation source remains debated. Here, we performed hybrid density functional calculations with spin-orbit coupling to investigate possible sources of hole compensation in group-V doped CdTe. First, regarding possible self-compensation of the group-V dopants, we find that the formation of AX centers is unlikely since they are found to be unstable relative to the shallow acceptor forms. However, if the group-V dopants come in during growth (such as dimer molecules As2, P2, and Sb2), we find that the impurity atoms, which would occupy nearest neighbor sites, maintain the V-V bonds, limiting the hole density. For the native defects, our study reveals that Cd interstitial is the lowest energy donor defect in p-type CdTe. Still, it has a small migration barrier of 0.5 eV, making it highly mobile and unstable at room temperature. The Te vacancy is the next lowest formation energy donor. The migration barrier of 1.4 eV indicates that the Te vacancy is stable at room temperature. The antisite CdTe is also a donor, with low formation energy and stable at room temperature, potentially limiting the hole concentration. Our results, therefore, shed light on possible compensation centers and some guidance on how to avoid them.

14 SOLAR ENERGY↗
DOE OSTI2024

Natural formation of linear defect structures in MoS 2

Near surface defects can significantly impact the quality of metallic interconnects and other interfaces necessary to create device structures incorporating two-dimensional materials. Furthermore, the impact of such defects can strongly depend on their organization. In this study, we present scanning tunneling microscopy images and tunneling spectroscopy of point and linear defects near the surface of natural MoS 2 . The point defects share similar structural and electronic characteristics and occur with comparable frequency as subsurface sulfur vacancies observed previously on natural MoS 2 . The linear defects observed here occur less frequently than the point defects but share the same depth profile and electronic structure. These data indicate that the linear defects are actually a one-dimensional organization of subsurface sulfur vacancies. Our density functional calculations agree with this assessment in that, for sufficient local defect concentrations, it is energetically more favorable for the defects to be organized in a linear fashion rather than as clusters or even isolated single point defects. Given these measurements were taken from naturally formed MoS 2 , this organization likely occurs during crystal formation. Considering the impact of one-dimensional organization on the local properties of layered materials, and the potential for them to be introduced purposefully during crystal formation, research into the formation mechanism and properties of these defects could enable new paths for defect engineering in MoS 2 -based systems.

2D materials↗
DOE OSTI2023

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of morphology in defect formation and photo-induced carrier instabilities in amorphous indium oxide

Ab initio molecular dynamics liquid-quench simulations and hybrid density functional calculations are performed to model the effects of room-temperature atomic fluctuations and photo-illumination on the structural and electronic properties of amorphous sub-stoichiometric In 2 O 2.96 . A large configurational ensemble is employed to reliably predict the distribution of localized defects as well as their response to the thermal and light activation. The results reveal that the illumination effects on the carrier concentration are greater in amorphous configurations with shorter In–O bond length and reduced polyhedral sharing as compared to the structures with a more uniform morphology. Further, the obtained correlation between the photo-induced carrier density and the reduction in the number of fully coordinated In-atoms implies that metal oxides with a significant fraction of crystalline/amorphous interfaces would show a more pronounced response to illumination. Photo-excitation also produces In–O 2 –In defects that have not been previously found in sub-stoichiometric amorphous oxides; these defects are responsible for carrier instabilities due to overdoping.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Modeling an S N 2 Reaction Mechanism for Hydrolysis of Siloxane Linkages with Density Functional Theory under Basic Conditions and Implications for Dissolution of Quartz

An extended silicate molecular cluster was used to perform density functional theory calculations to determine if an S N 2 mechanism, where OH – attacks a Q 1 Si coupled with Na + charge balancing an adjacent siloxane bridging O atom, could explain the observed energy of activation of quartz dissolution under basic conditions.

chemical calculations↗
DOE OSTI2022

Evidence for electronic signature of a magnetic transition in the topological magnet HoSbTe

Topological insulators with intrinsic magnetic order are emerging as an exciting platform to realize fundamentally new excitations from topological quantum states of matter. To study these systems and their physics, people have proposed a variety of magnetic topological insulator systems, including HoSbTe, an antiferromagnetic weak topological insulator candidate. Here, in this work, we use scanning tunneling microscopy to probe the electronic structure of HoSbTe with antiferromagnetic and ferromagnetic orders that are tuned by applying an external magnetic field. Although around the Fermi energy we find minor differences between the quasiparticle interferences under the ferromagnetic and antiferromagnetic orders, deep inside the valance region, a quasiparticle interference signal emerges with ferromagnetism. This observation is consistent with our first-principles calculations indicating the magnetism-driven transition of the electronic states in this spin-orbit coupled topological magnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(C0)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule 04(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in 04(+) has also been studied using the Brueckner coupled cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗
NASA NTRS1994

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(CO)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule O4(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI,and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in O4(+) has also been studied using the Brueckner coupled-cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗
DOE OSTI2020

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Surface facet dependence of competing alloying mechanisms

Metal alloys are ubiquitous in many branches of heterogeneous catalysis, and it is now fairly well established that the local atomic structure of an alloy can have a profound influence on its chemical reactivity. While these effects can be difficult to probe in nanoparticle catalysts, model studies using well defined single crystal surfaces alloyed with dopants enable these structure–function correlations to be drawn. The first step in this approach involves understanding the alloying mechanism and the type of ensembles formed. In this study, we examined the atomic structure of RhCu single-atom alloys formed on Cu(111), Cu(100), and Cu(110) surfaces. Our results show a striking difference between Rh atoms alloying in Cu(111) vs the more open Cu(100) and Cu(110) surface facets. Unlike Cu(111) on which Rh atoms preferentially place-exchange with Cu atoms in the local regions above step edges leaving the majority of the Cu surface free of Rh, highly dispersed, homogeneous alloys are formed on the Cu(100) and (110) surfaces. These dramatically different alloying mechanisms are understood by quantifying the energetic barriers for atomic hopping, exchange, swapping, and vacancy filling events for Rh atoms on different Cu surfaces through theoretical calculations. Density functional theory results indicate that the observed differences in the alloying mechanism can be attributed to a faster hopping rate, relatively high atomic exchange barriers, and stronger binding of Rh atoms in the vicinity of step edges on Cu(111) compared to Cu(110) and Cu(100). Furthermore, these model systems will serve as useful platforms for examining structure sensitive chemistry on single-atom alloys.

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

Equation of state of boron carbide B 4 ⁢C

We present the results of recent experiments conducted on the Sandia Z machine and a new tabular equation of state for B 4 ⁢C. The equation of state was calibrated to a combination of density functional calculations reported here and fits to preexisting data. It was constructed partly to recover the effects of a shock-driven, polymorphic phase transition of unknown character beginning at particle velocities of just under 3 km/s (shock pressures of 95 GPa). Some of the Z experiments included sound speeds determined by the overtaking rarefaction method, from which we calculate the Grüneisen parameter and compare with previous experiments conducted at the OMEGA laser [Fratanduono et al ., Phys. Rev. B 94 , 184107 (2016)], our own first principles calculations, and another recent tabular equation of state [Zhang et al ., Phys. Rev. E 102 , 053203 (2020)]. We also compare our results with previous static compression, thermophysical, and melt studies, finding mixed consistency. We predict the onset and completion of shock melting at 225 and 265 GPa, respectively, and predict a melt curve that is largely flat to pressures of several hundred GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Li 8 MnO 6 : A Novel Cathode Material with Only Anionic Redox

In Li-excess transition-metal-oxide cathode materials, anionic oxygen redox can offer high capacity and high voltages, although peroxo and superoxo species may cause oxygen loss, poor cycling performance, and capacity fading. Previous work showed that undesirable formation of peroxide and superoxide bonds is controlled to some extent by Mn substitution, and the present work uses density functional calculations to examine the reasons for this by studying the anionic redox mechanism Li 8 MnO 6 . This material is obtained by substituting Mn for Sn in Li 8 SnO 6 or for Zr in Li 8 ZrO 6 , and we also compare to previous work on those materials. The calculations predict that Li 8 MnO 6 is stable at room temperature (with a band gap of 3.19 eV as calculated HSE06 and 1.82 eV as calculated with the less reliable with PBE+U), and they elucidate the chemical and structural effects involved in the inhibition of oxygen release in this cathode. Throughout the whole delithiation process, only O 2- ions are oxidized. The directional Mn-O bonds formed from unfilled 3d orbitals effectively inhibit the formation of O-O bonds, and the layered structure is maintained even after removing 3 Li per Li 8 MnO 6 formula unit. The calculated average voltage for removal of 3 Li is 3.69 V by HSE06, and the corresponding capacity is 389 mAh/g. The high voltage of oxygen anionic redox and the high capacity result in a high energy density of 1436 Wh/kg. The Li-ion diffusion barrier for the dominant interlayer diffusion path along the c-axis is 0.57 eV by PBE+U. Finally, these results help us to understand the oxygen redox mechanism in a new lithium-rich Li 8 MnO 6 cathode material and contribute to the design of high-energy-density lithium-ion-battery cathode materials with favorable electrochemical properties based on anionic oxygen redox.

Li8MnO6↗
DOE OSTI2021

Origin of the magnetic and orbital ordering in α-Sr 2 CrO 4

Motivated by recent experimental progress in transition metal oxides with the K 2 NiF 4 structure, in this study we investigate the magnetic and orbital ordering in α-Sr 2 CrO 4 . Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the ab initio band structure near the Fermi level. The unique reverse splitting of t 2g orbitals in α-Sr 2 CrO 4 , with the 3d 2 electronic configuration for the Cr 4+ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground state phase diagram for the two dimensional compound α-Sr 2 CrO 4 . We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of α-Sr 2 CrO 4 , we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound α-Sr 2 CrO 4 , in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with d xy driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the d xz and d yz orbitals driving the orbital order.

36 MATERIALS SCIENCE↗