Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 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 OSTI2025

Monolayer TiS 2 Nanosheets on Au(111)–Structural Characterization and Effect of Edge Stability for Shape Control

Transition metal dichalcogenides are promising alternatives to noble metal catalysts, e.g., for (photo-)activation of greenhouse gases or hydrogenations. Herein, a direct synthetic route for 2D TiS 2 nanosheets on Au(111) by titanium deposition in the presence of a mild, organic, non-oxidizing sulfur source is presented. High-resolution scanning tunneling microscopy (STM) is used to gain atomic-level insights into the TiS 2 nanosheet morphology. In contrast to the literature, this protocol gains mostly hexagonal and truncated triangular nanosheets with an increased edge contrast in STM, analog to metallic edge states in MoS 2 . Synchrotron-based photoelectron spectroscopy allows insights into compositional details, specifically to distinguish different S sites on the TiS 2 sheets and other S species on the sample. Further, a minimum size is identified (9 S atoms side length), which underlines the importance of moiré reconstructions for stress relief. The TiS 2 sheets coexist with [Au]Ti 1 S 3 clusters, in which a single gold atom is alloyed into the surface and capped by three S atoms. Together with the finding of a critical sheet size, this points toward on-surface Ostwald ripening as a relevant process in the sheet formation. Ab-initio calculations (density functional theory) underscore that the chemical potential of S is an essential descriptor to maintain shape control.

2D nanomaterials↗
DOE OSTI2024

Stabilizing lithium superoxide formation in lithium-air batteries by Janus chalcogenide catalysts

Solid lithium peroxide (Li 2 O 2 ) is the major discharge product in Li-air batteries. However, the electronically insulating nature of Li 2 O 2 tends to affect the battery’s performance such as the polarization gap and cyclability. On the other hand, lithium superoxide (LiO 2 ), generated through a one-electron transfer process, offers greater electronic conductivity, lower charge transfer resistance, and thus reduced charge potential. Nevertheless, LiO 2 long-term stabilization as a final product remains a significant challenge. Here, in this study, we present the molybdenum (Mo)-based Janus chalcogenide family featuring asymmetric structures as a new generation of cathode catalysts for Li-air batteries. These catalysts demonstrate remarkable efficacy in stabilizing LiO 2 discharge products, even under high current densities of 5000 mA/g (corresponding to 0.5 mA/cm 2 ). Our density functional calculations provide an understanding of why the asymmetric Mo-Janus chalcogenides result in LiO 2 formation whereas the symmetric Mo-dichalcogenides produce Li 2 O 2 as the discharge product. These results pave the way to explore a new generation of advanced catalysts for superoxide-based Li-air batteries.

Chalcogenide↗
DOE OSTI2024

High-Throughput Screening for Boride Superconductors

A high-throughput screening using density functional calculations is performed to search for stable boride superconductors from the existing materials database. The workflow employs the fast frozen-phonon method as the descriptor to evaluate the superconducting properties quickly. Twenty-three stable candidates were identified during the screening. The superconductivity was obtained earlier experimentally or computationally for almost all found binary compounds. Previous studies on ternary borides are very limited. Here our extensive search among ternary systems confirmed superconductivity in known systems and found several new compounds. Among these discovered superconducting ternary borides, TaMo 2 B 2 shows the highest superconducting temperature of ∼12 K. Most predicted compounds were synthesized previously; therefore, our predictions can be examined experimentally. Our work also demonstrates that the boride systems can have diverse structural motifs that lead to superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Defect Stability in CdTe Based on Formation Energies and Migration Barriers

Native point defects are thought to play a key role in CdTe, either as compensation centers in intentionally doped material, as a source of conductivity in nominally undoped material, or as electron–hole recombination centers. Here, the discussion of their concentration and impact has often centered only on formation energies and transition levels. Using hybrid density functional calculations, including the effects of spin–orbit coupling (SOC), we discuss the stability of native point defects in CdTe based on their formation energies and migration barriers. We show that although Cd interstitials are the lowest energy donor defects, they are unstable at room temperature due to a low migration barrier. They are important for maintaining charge neutrality during growth or annealing at high temperatures, but once the material is brought to room temperature, they are not frozen in as often assumed and are expected to anneal out, leaving the other more stable defects to determine the conductivity. Taking this into account in the solution of the charge neutrality equation, we are able to predict the conductivity type and carrier concentrations that are in good agreement with experimental observations.

14 SOLAR ENERGY↗
DOE OSTI2025

Selective Adsorption of Thiol-Containing Molecules on Copper Sulfide Surfaces via Molecule–Surface Disulfide Bridges

Recent results in the fields of nanoenhanced agriculture and expanding interest in prebiotic chemistry have placed increased emphasis on understanding the chemically selective interaction of small molecules with the surfaces of metal sulfides. Here, we present an integrated experimental and computational study of the interaction of thiol-containing molecules with copper sulfide (covellite) surfaces in aqueous media. In situ Fourier-transform infrared (FTIR) measurements and ex situ X-ray photoelectron spectroscopy (XPS) measurements show that molecules bearing free thiol groups, including glutathione and cysteine, bind strongly to CuS (covellite) nanoparticles and to CuS (001) single crystals, while control studies show that similar molecules lacking the free thiol group exhibit much less binding. Additional experiments show that these thiol-containing molecules interact transiently with CuO nanoparticle surfaces but are readily removed by rinsing. The FTIR and XPS experiments demonstrate that adsorption of molecular thiols to CuS surfaces occurs in a chemically selective manner. Further experimental studies and density functional calculations show that the preferred mode of binding is through the surface S atoms, forming a Solid–S–S–Molecule disulfide linkage. While the role of disulfide linkages in controlling structure and function of proteins and other biomolecules is widely known, the formation of surface disulfide linkages as a motif for covalent molecular binding at surfaces has not been established previously.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Untangling Sources of Error in the Density-Functional Many-Body Expansion

The many-body expansion provides a framework for data-driven applications of electronic structure theory, including parametrization of classical force fields and machine learning. In this article, we demonstrate that its use significantly amplifies quadrature grid errors when modern density-functional approximations are employed. Standard grids that work well in conventional density-functional calculations result in runaway error accumulation when used with the many-body expansion. At the same time, delocalization error is also exacerbated, leading to exaggerated estimates of nonadditive n-body interactions. This is illustrated for anion–water clusters using the SCAN, r2SCAN, ωB97X-V and ωB97M-V functionals. By employing dense quadrature grids, the inherent self-interaction error is exposed, which can then be mitigated using a variety of other strategies.

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

Controlling Noncollinear Ferromagnetism in van der Waals Metal–Organic Magnets

Van der Waals (vdW) magnets both allow exploration of fundamental 2D physics and offer a route toward exploiting magnetism in next generation information technology, but vdW magnets with complex, noncollinear spin textures are currently rare. We report here the syntheses, crystal structures, magnetic properties and magnetic ground states of four bulk vdW metal–organic magnets (MOMs): FeCl 2 (pym), FeCl 2 (btd), NiCl 2 (pym), and NiCl 2 (btd), pym = pyrimidine and btd = 2,1,3-benzothiadiazole. Using a combination of neutron diffraction and bulk magnetometry we show that these materials are noncollinear magnets. Although only NiCl 2 (btd) has a ferromagnetic ground state, we demonstrate that low-field hysteretic metamagnetic transitions produce states with net magnetization in zero-field and high coercivities for FeCl 2 (pym) and NiCl 2 (pym). By combining our bulk magnetic data with diffuse scattering analysis and broken-symmetry density-functional calculations, we probe the magnetic superexchange interactions, which when combined with symmetry analysis allow us to suggest design principles for future noncollinear vdW MOMs. These materials, if delaminated, would prove an interesting new family of 2D magnets.

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

Hidden magnetism and split off flat bands in the insulator metal transition in VO 2

Transition metal d -electron oxides with an odd number of electrons per unit cell are expected to form metals with partially occupied energy bands, but exhibit in fact a range of behaviors, being either insulators, or metals, or having insulator-metal transitions. Traditional explanations involved predominantly electron-electron interactions in fixed structural symmetry. The present work focuses instead on the role of symmetry breaking local structural motifs. Viewing the previously observed V-V dimerization in VO 2 as a continuous knob, reveals in density functional calculations the splitting of an isolated flat band from the broad conduction band. This leads past a critical percent dimerization to the formation of the insulating phase while lowering the total energy. In VO 2 this transition is found to have a rather low energy barrier approaching the thermal energy at room temperature, suggesting energy-efficient switching in neuromorphic computing. Interestingly, sufficient V-V dimerization suppresses magnetism, leading to the nonmagnetic insulating state, whereas magnetism appears when dimerization is reduced, forming a metallic state. This study opens the way to design novel functional quantum materials with symmetry breaking-induced flat bands.

Chemistry↗
DOE OSTI2026

Magnetization orientation-dependent Shubnikov-de Haas oscillations in ferromagnetic Weyl semimetal Co 3 Sn 2 S 2

We report a study of Shubnikov–de Haas oscillations in high-quality single crystals of ferromagnetic Weyl semimetal Co 3 Sn 2 S 2 . The Fermi surfaces resolved in our experiments are three-dimensional and reflect an underlying trigonal crystallographic symmetry. Combined with density functional calculations, we identify that multiple Fermi surfaces in the system—of both electron and hole nature—arise from the energy dispersion of the (spin-orbit gapped) mirror-protected nodal rings. We observe an evolution of the Fermi surfaces with in-plane magnetic fields, in contrast to field perpendicular to the kagome lattice planes, which has little effect. Viewed alongside the easy-axis anisotropy of the system, our observation reveals an evolution of the electronic structure of Co 3 Sn 2 S 2 —including the Weyl points—with the ferromagnetic moment orientation. Through the case study of Co 3 Sn 2 S 2 , our results provide concrete experimental evidence of an anisotropic interplay via spin-orbit coupling between the magnetic degrees of freedom and electronic band singularities, which has long been expected in semimetallic and metallic magnetic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pressure induced modification of electronic and magnetic properties of MnCrNbAl and MnCrTaAl

Spin-gapless semiconductor (SGS) is a new class of material that has been studied recently for potential applications in spintronics. This material behaves as an insulator for one spin channel, and as a gapless semiconductor for the opposite spin. In this work, we present results of a computational study of two quaternary Heusler alloys, MnCrNbAl and MnCrTaAl that have been recently reported to exhibit spin-gapless semiconducting electronic structure. In particular, using density functional calculations we analyze the effect of external pressure on electronic and magnetic properties of these compounds. It is shown that while these two alloys exhibit nearly SGS behavior at optimal lattice constants and at negative pressure (expansion), they are half-metals at equilibrium, and magnetic semiconductors at larger lattice constant. At the same time, reduction of the unit cell volume has a detrimental effect on electronic properties of these materials, by modifying the exchange splitting of their electronic structure and ultimately destroying their half-metallic/semiconducting behavior. Thus, our results indicate that both MnCrNbAl and MnCrTaAl may be attractive for practical device applications in spin-based electronics, but a potential compression of the unit cell volume (e.g. in thin-film applications) should be avoided.

Materials Science↗
DOE OSTI2024

Utilizing integrated neutron diffraction and elastoplastic self-consistent crystal plasticity model to quantitatively assess the strengthening mechanism in Al–12.5Ce and Al–12.5Ce–0.4Mg alloys

An integrated in-situ neutron diffraction and elastic plastic self-consistent crystal plasticity (EPSC-CP) modeling scheme is performed on a binary Al–12Ce alloy and a ternary Al–12Ce–0.4Mg alloys. Using this scheme, the constitutive parameters, i.e. elastic constants and slip system parameters of individual phases can be calibrated which can be used in microstructure-based CP models to predict materials performance. From this study, it is shown that the elastic constants of Al 11 Ce 3 intermetallics calculated from density function theory calculation in the literature are rather accurate. When applied to the EPSC-CP model, the lattice strains of both the binary and ternary alloys are correctly predicted as compared with experiments, and large lattice strain differences between Al (100) plane and Al 11 Ce 3 (010) plane are demonstrated. The slip system parameters calibrated by the scheme shows that the addition of 0.4 wt% Mg in the alloy has little influence on the critical resolved shear stress of initial dislocation glide in the Al matrix which caused plastic yield in the material. This can be explained by the very dilute Mg solute content in the Al solid solution, causing large spacing of Al–Mg lattice misfit sites and little impact on resistance of dislocation glide at initial yield. The 0.4 wt% Mg addition, on the other hand, has a large influence on the hardening term in the slip system parameters, indicating those Al–Mg misfit sites do help dislocation accumulation during the deformation. The impact of dilute Mg addition on the Al slip system parameters is also reflected in the flow behavior of the ternary alloy: little impact on the yield stress, but a large impact on working hardening and tensile strength of the materials which is consistent with the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dipolar and quadrupolar correlations in the 5⁢𝑑 2 Re-based double perovskites Ba 2 ⁢ YReO 6 and Ba 2 ⁢ ScReO 6

Double perovskites containing heavy transition-metal ions are an important family of compounds for the study of the interplay between electron correlation and spin-orbit coupling. Here, by combining magnetic susceptibility, heat capacity, and neutron scattering measurements, we investigate the dipolar and quadrupolar correlations in two prototypal rhenium-based double perovskite compounds Ba 2 ⁢YReO 6 and Ba 2 ⁢ ScReO 6 . A type-I dipolar antiferromagnetic ground state with a propagation vector 𝐪 = (0,0,1) is observed in both compounds. At temperatures above the magnetic transitions, a quadrupolar ordered phase is identified. Weak spin excitations, which are gapped at low temperatures and softened in the correlated paramagnetic phase, are explained using a minimal model that considers both the dipolar and quadrupolar interactions. In conclusion, at larger wave vectors, we observe dominant phonon excitations that are well described by density functional calculations.

Omar, Otkur [University of Science and Technology ↗
DOE OSTI2024

Influence of temperature, doping, and amorphization on the electronic structure and magnetic damping of iron

Hybrid magnonic quantum systems have drawn increased attention in recent years for coherent quantum information processing, but too large magnetic damping is a persistent concern when metallic magnets are used. Their intrinsic damping is largely determined by electron-magnon scattering induced by spin-orbit interactions. In the low scattering limit, damping is dominated by intra-band electronic transitions, which has been theoretically shown to be proportional to the electronic density of states at the Fermi level. In this work, we focus on body-centered-cubic iron as a paradigmatic ferromagnetic material. We comprehensively study its electronic structure using first-principles density functional theory simulations and account for finite lattice temperature, boron (B) doping, and structure amorphization. Our results indicate that temperature induced atomic disorder and amorphous atomic geometries only have a minor influence. Instead, boron doping noticeably decreases the density of states near the Fermi level with an optimal doping level of 6.25%. In addition, we show that this reduction varies significantly for different atomic geometries and report that the highest reduction correlates with a large magnetization of the material. Furthermore, this may suggest materials growth under external magnetic fields as a route to explore in experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Application of machine learning to discover new intermetallic catalysts for the hydrogen evolution and the oxygen reduction reactions

The adsorption energies for hydrogen, oxygen, and hydroxyl were calculated by means of density functional theory on the lowest energy surface of 24 pure metals and 332 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B taking into account the effect of biaxial elastic strains. This information was used to train two random forest regression models, one for the hydrogen adsorption and another for the oxygen and hydroxyl adsorption, based on 9 descriptors that characterized the geometrical and chemical features of the adsorption site as well as the applied strain. All the descriptors for each compound in the models could be obtained from physico-chemical databases. The random forest models were used to predict the adsorption energy for hydrogen, oxygen, and hydroxyl of ≈2700 binary intermetallic compounds with stoichiometries AB, A 2 B, and A 3 B made of metallic elements, excluding those that were environmentally hazardous, radioactive, or toxic. This information was used to search for potential good catalysts for the HER and ORR from the criteria that their adsorption energy for H and O/OH, respectively, should be close to that of Pt. Further, this investigation shows that the suitably trained machine learning models can predict adsorption energies with an accuracy not far away from density functional theory calculations with minimum computational cost from descriptors that are readily available in physico-chemical databases for any compound. Moreover, the strategy presented in this paper can be easily extended to other compounds and catalytic reactions, and is expected to foster the use of ML methods in catalysis.

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

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗
DOE OSTI2024

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗