Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spectroscopic Evidence of Hyponitrite Radical Intermediate in NO Disproportionation at a MOF-Supported Mononuclear Copper Site

Dianionic hyponitrite (N 2 O 2 2- ) is often proposed, based on model complexes, as the key intermediate in reductive coupling of nitric oxide to nitrous oxide at the bimetallic active sites of heme-copper oxidases and nitric oxide reductases. In this work, we examine the gas-solid reaction of nitric oxide with the metal–organic framework Cu I -ZrTpmC* with a suite of in situ spectroscopies and density functional theory simulations, and identify an unusual chelating N 2 O 2 .- intermediate. These results highlight the advantage provided by site-isolation in metal–organic frameworks (MOFs) for studying important reaction intermediates, and provide a mechanistic scenario compatible with the proposed one-electron couple in these enzymes.

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

Spectroscopic Evidence of Hyponitrite Radical Intermediate in NO Disproportionation at a MOF‐Supported Mononuclear Copper Site

Abstract Dianionic hyponitrite (N 2 O 2 2− ) is often proposed, based on model complexes, as the key intermediate in reductive coupling of nitric oxide to nitrous oxide at the bimetallic active sites of heme‐copper oxidases and nitric oxide reductases. In this work, we examine the gas‐solid reaction of nitric oxide with the metal–organic framework Cu I ‐ZrTpmC* with a suite of in situ spectroscopies and density functional theory simulations, and identify an unusual chelating N 2 O 2 .− intermediate. These results highlight the advantage provided by site‐isolation in metal–organic frameworks (MOFs) for studying important reaction intermediates, and provide a mechanistic scenario compatible with the proposed one‐electron couple in these enzymes.

Sun, Chenyue↗
DOE OSTI2024

3D Lead-Organoselenide-Halide Perovskites and their Mixed-Chalcogenide and Mixed-Halide Alloys

Here, we incorporate Se into the 3D halide perovskite framework using the zwitterionic ligand: SeCYS ( + NH 3 (CH 2 ) 2 Se - ), which occupies both the X - and A + sites in the prototypical ABX 3 perovskite. The new organoselenide-halide perovskites: (SeCYS)PbX 2 (X=Cl, Br) expand upon the recently discovered organosulfide-halide perovskites. Single-crystal X-ray diffraction and pair distribution function analysis reveal the average structures of the organoselenide-halide perovskites, whereas the local lead coordination environments and their distributions were probed through solid-state 77 Se and 207 Pb NMR, complemented by theoretical simulations. Density functional theory calculations illustrate that the band structures of (SeCYS)PbX 2 largely resemble those of their S analogs, with similar band dispersion patterns, yet with a considerable band gap decrease. Optical absorbance measurements indeed show band gaps of 2.07 and 1.86 eV for (SeCYS)PbX 2 with X=Cl and Br, respectively. We further demonstrate routes to alloying the halides (Cl, Br) and chalcogenides (S, Se) continuously tuning the band gap from 1.86 to 2.31 eV–straddling the ideal range for tandem solar cells or visible-light photocatalysis. The comprehensive description of the average and local structures, and how they can fine-tune the band gap and potential trap states, respectively, establishes the foundation for understanding this new perovskite family, which combines solid-state and organo-main-group chemistry.

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

Cu Evolution over Bimetallic Cu‐Y/Beta Zeolite Under H 2 and Ethanol Atmospheres: Unveiling the Role of Diatomic Metal–Metal Interactions

Understanding the dynamic evolution of Cu species under varying environmental conditions is critical for addressing challenges related to the activity and the stability of copper‐based catalysts in thermo‐, photo‐, and electrocatalysis. However, metal–metal interactions between dual single atoms and their effects on Cu evolution after exposure to different environmental molecules remain underexplored. Herein, we synthesized bimetallic Cu‐Y/Beta catalysts with dual single‐atom Cu and Y sites and monometallic Cu‐Beta catalysts with isolated Cu sites in dealuminated Beta zeolites. By varying Cu and Y compositions, diatomic interactions were studied under H 2 and ethanol atmospheres. With 6 wt% Y loading, approximately 0.4 wt% of Cu species in Cu‐Y/Beta remained partially oxidized as Cu(I) after reduction in pure H 2 at 350 °C, in contrast to the full transition to metallic Cu observed in Cu‐Beta. Combining X‐ray absorption spectroscopy with kinetic studies revealed that metallic Cu became the predominant species after reduction with H 2 as Cu loading increased from 0.4 to 1.7 wt%, quadrupling the initial ethanol dehydrogenation rate and demonstrating the dominant role of Cu(0) sites. In conclusion, scanning transmission electron microscopy and density functional theory simulations indicated spatial proximity between dual single‐atom Cu and Y sites and elucidated Cu speciation controlled by diatomic interactions.

Copper↗
DOE OSTI2024

Interpretable Machine Learning Models for Practical Antimonate Electrocatalyst Performance

Computationally predicting the performance of catalysts under reaction conditions is a challenging task due to the complexity of catalytic surfaces and their evolution in situ, different reaction paths, and the presence of solid-liquid interfaces in the case of electrochemistry. We demonstrate here how relatively simple machine learning models can be found that enable prediction of experimentally observed onset potentials. Inputs to our model are comprised of data from the oxygen reduction reaction on non-precious transition-metal antimony oxide nanoparticulate catalysts with a combination of experimental conditions and computationally affordable bulk atomic and electronic structural descriptors from density functional theory simulations. From human-interpretable genetic programming models, we identify key experimental descriptors and key supplemental bulk electronic and atomic structural descriptors that govern trends in onset potentials for these oxides and deduce how these descriptors should be tuned to increase onset potentials. Here, we finally validate these machine learning predictions by experimentally confirming that scandium as a dopant in nickel antimony oxide leads to a desired onset potential increase. Macroscopic experimental factors are found to be crucially important descriptors to be considered for models of catalytic performance, highlighting the important role machine learning can play here even in the presence of small datasets.

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

First-principles analysis of the Al-rich corner of Al-Li-Cu phase diagram

The phase diagram of Al-Li-Cu system in the Al-rich region was determined by means of first-principles calculations and statistical mechanics. The mixing enthalpies of many configurations for different lattices in the whole Al-Li-Cu system were determined by density functional theory simulations to find the stable phases in the convex hull. They were fitted with a cluster expansion to calculate the free energy of the configurations with different compositions as a function of temperature in the Al-rich region (Al content > 40 at.%) by means of Monte Carlo simulations. It was found that the ground state phases in the Al-rich part of the Al-Li-Cu phase diagram were α-Al, θ' (Al 2 Cu), $δ$' (Al 3 Li), $δ$ (AlLi) and T 1 (Al 6 Cu 4 Li 3 ), while θ'' (Al 3 Cu), T 1' (Al 2 CuLi) and Al 3 Cu 2 Li were found on the lowest mixing enthalpy surfaces of their lattices and were metastable. α-Al, $δ$ and T 1 are stable phases in the whole temperature range while $δ$' becomes metastable at very low temperature and θ (Al 2 Cu) replaces θ' as the stable phase at approximately 550 K due to the vibrational entropic contribution. In addition, the phase diagram in the Al-rich region was built and it was shown in isothermal sections from 100 to 900 K. They were in good agreement with the limited experimental data in the literature and provided new information regarding the stability, solubility and stoichiometry of the different phases. This information is important to understand the precipitation mechanisms during high temperature aging.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formation of pyrophosphates across grain boundaries induces the formation of mismatched but oriented interfaces in silver phosphate polypods

The interfaces and their misfit defects determine the materials properties for a wide range of applications, such as electronic devices, photocatalysis, and mechanical engineering, etc. However, current understanding of atomic interfacial structures is limited. Here we discover a special interfacial structure, mismatched but oriented interface via two distinct facets. Transmission and scanning electron microscopy results suggest that, in Ag 3 PO 4 polypods structures, interfaces of {1 0 0} and {1 1 0}, ({1 0 0}/(1 1 0}), {1 0 0}/{1 1 1}, {1 1 0}/{1 1 1}, and {1 0 0}/{1 0 0}, etc., have a certain orientation relationship, corresponding to the energy minima and coincident site lattice of interfacial atoms as demonstrated by molecular dynamics simulations. Density functional theory demonstrates that the formation of pyrophosphate and/or phosphates rotation to bond across the interface compensate the lattice mismatch at the interfaces, as well as deformations of Ag-O bonds. Furthermore, our work opens up a new avenue for a much wider range of interfacial structures, allow for a higher diversity of structures, and shine light on tailoring crystal structures, morphologies, and the resulting properties.

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

Hydrogen Localization and Cluster Formation in α-Zr from First-Principles Investigations

We report alpha-zirconium (α-Zr) is a hydrogen atom host material used in nuclear energy and semiconducting applications due to Zr's excellent properties such as corrosion resistance, low thermal neutrons cross-section, and thermodynamically favorable interstitial sequestering of H atoms. Herein we use density functional theory simulations to systematically analyze the initial stages of H insertion in the Zr matrix, and determine optimal hydrogen localization, distribution, and clustering tendencies from a thermodynamics and kinetics perspective. Our results suggest that H occupation of tetrahedral sites is strongly favored, consistent with previous experimental and theoretical studies. Further, a cooperative effect in H localization is shown where H atoms favor clustering. The H pair is found to be the most stable by occupying tetrahedral sites with pairs at second nearest-neighbor positions (TTn2). Among the four H atom clusters, the most favorable configuration has two TTn2 pairs forming "box " or "zig-zag" shaped-cluster with the pairs separated by 3.03 or 3.25 Å, respectively. Nudged elastic band simulations suggest that the formation of TTn2 is facile and is limited by the small diffusion barrier of the monomer of less than 0.1 eV. Further, the formation routes for the 4H atom clusters are also presented. The current study can serve as a foundation to systematically model H atom clustering in metals and alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The role of moisture in MgCl 2 salt: A multiscale approach to TES performance

This study presents a comprehensive multiscale analysis to evaluate the influence of moisture on the thermal performance of thermal energy storage systems using magnesium chloride (MgCl 2 ) as the phase-changing material. The system uses graphite foam with 90% relative density to enhance thermal conductivity. The analysis includes thermal conductivity calculations and specific heat capacity for different hydrate phases of MgCl 2 : Anhydrous, Mono, Di, Tetra, and Hexa. These properties were evaluated for the first time using the phonon density of states from Density Functional Theory simulations. Results showed that thermal conductivity decreased, while specific heat capacity increased by a factor of two as the phase changed from anhydrous to hexahydrate. Meso-scale models were created to account for the anisotropy of graphite foam and property variations of the MgCl 2 hydrate phases. Asymptotic Expansion homogenization simulations determined the anisotropic thermal conductivity for all phases. This unique methodology improved simulation accuracy, which matched experimental data for anhydrous MgCl 2 . A parametric study examined various operating conditions and their effect on TES performance. It revealed that higher charging temperatures did not enhance exergy efficiency, but increased discharge mass flow rates improved it due to better heat transfer. Thermal performance evaluated by the exergy efficiency remained consistent across hydrate systems under the tested conditions. Furthermore, the study suggests that this uniformity is linked to missing key data, particularly latent heat and melting point for different hydrates. Overall, it highlights the importance of multiscale effects and accurate material properties in designing and optimizing TES systems.

Molten salt degradation↗
DOE OSTI2026

A generative machine learning model for designing metal hydrides applied to hydrogen storage

Developing new metal hydrides is a critical step toward efficient hydrogen storage in carbon-neutral energy systems. However, existing materials databases, such as the Materials Project, contain a limited number of well-characterized hydrides, which constrains the discovery of optimal candidates. This work presents a framework that integrates causal discovery with a lightweight generative machine learning model to generate novel metal hydride candidates that may not exist in current databases. Using a dataset of 450 samples (270 training, 90 validation, and 90 testing), the model generates 1000 candidates. After ranking and filtering, six previously unreported chemical formulas and crystal structures are identified, four of which are validated by density functional theory simulations and show strong potential for future experimental investigation. Overall, the proposed framework provides a scalable and time-efficient approach for expanding hydrogen storage datasets and accelerating materials discovery.

generative model↗
DOE OSTI2025

Pd-promoted reduction and restructuring of an In 2 O 3 -based catalyst for CO 2 hydrogenation at room temperature

An unconventional reaction mechanism in an In 2 O 3 /Pd(1 1 1) inverse model catalyst for the CO 2 hydrogenation reaction has been uncovered: In 2 O 3 is partially reduced at room temperature in a reaction atmosphere as a result of its direct contact with Pd(1 1 1), which is an efficient H 2 splitter. The reduction induces changes in surface free energy, leading to a dynamical restructuring at the In 2 O 3 /Pd(1 1 1) interface via formation of InO x and outward diffusion of Pd, as revealed by ambient pressure X-ray photoelectron spectroscopy, X-ray absorption spectroscopy and density functional theory simulations. This dynamical restructuring eventually promotes the growth of 2D InPd y O x nanodomains as the catalytically active phase and the exclusive formation of methanol upon hydrogenation of CO 2 at room temperature. A comparable high selectivity toward CH 3 OH was found in more realistic bulk catalytic systems (2 wt% Pd/In 2 O 3 catalyst and commercial CZA catalyst). Scanning tunneling microscopy under ultrahigh vacuum and ambient pressure reaction atmospheres further reveals the structural dynamics at the InO x /Pd(1 1 1) interface, where we follow in situ the evolution of the InO x particles on Pd(1 1 1) and the mobility of the InPd y O x nanodomains in a CO 2 + H 2 environment. The present findings of the formation of a mixed oxide phase in a dynamically restructuring metal/reducible-oxide interface indicate further implications for other heterogeneous catalytic systems beyond the present CO 2 hydrogenation example and highlight the importance of in situ investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Magnetism at the interface of MoSe 2 /V 2 O 5 heterostructures

In this study, we report the existence of magnetism in a heterostructure made up of an MoSe 2 transition metal dichalcogenide monolayer and a V 2 O 5 substrate. Using density functional theory simulations, we find that ferromagnetic ordering can be found in the MoSe 2 /V 2 O 5 heterostructure even though the individual components are nonmagnetic. The role of van der Waals interactions is examined.

2D magnet↗
DOE OSTI2020

Engineering the electronic and strained interface for high activity of PdM core @Pt monolayer electrocatalysts for oxygen reduction reaction

Alloyed nanoparticles with core-shell structures provide a favorable model to modulate interfacial interaction and surface structures at the atomic level, which is important for designing electrocatalysts with high activity and durability. Herein, core-shell structured Pd 3 M@Pt/C nanoparticles with binary PdM alloy cores (M = Fe, Ni, and Co) and a monolayer Pt shell were successfully synthesized with diverse interfaces. Among these, Pd 3 Fe@Pt/C exhibited the best oxygen reduction reaction catalytic performance, roughly 5.4 times more than that of the commercial Pt/C catalyst used as reference. The significantly enhanced activity is attributed to the combined effects of strain engineering, interfacial electron transfer, and improved Pt utilization. Density functional theory simulations and extended X-ray absorption fine structure analysis revealed that engineering the alloy core with moderate lattice mismatch and alloy composition (Pd 3 Fe) optimizes the surface oxygen adsorption energy, thereby rendering excellent electrocatalytic activity. Finally, future researches may use this study as a guide on the construction of highly effective core-shell electrocatalysts for various energy conversions and other applications.

25 ENERGY STORAGE↗
DOE OSTI2022

Effect of single atom Platinum (Pt) doping and facet dependent on the electronic structure and light absorption of Lanthanum Titanium Oxide (La 2 Ti 2 O 7 ): A Density Functional Theory study

Charge generation and separation are two key processes for semiconductor photocatalysis. Here, we use Pt as a single atom catalyst to systematically examine the facet-dependent electronic band structure and light absorption of the layered perovskite-type wide-gap semiconductor, lanthanum titanium oxide (La 2 Ti 2 O 7 , LTO) by means of density functional theory simulations. In this study, it is found that single Pt atom doping of different LTO surfaces (here, (100), (101) and (001)) can not only create states in the bandgap that would promote the formation of recombination centers, but also shift the optical absorption edge to the visible region. Interestingly, the Pt doping forms a heterojunction, with the valence band maximum consisting of Pt states and the conduction band minimum consisting of LTO states, respectively. To study the facet dependent surface activity, acetic acid (CH 3 COOH) was used as a model molecule to investigate the adsorption and charge transfer on the (101), (100) and (001) Pt-LTO surface facets. The results show that the (101) facet could enable stronger adsorption of CH 3 COOH by promoting more electron transfer during the interfacial interaction. Our theoretical findings aim to promote the design and optimization of the single atom catalysts (SACs) for photocatalytic applications and other broad catalysis systems.

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

First-Principles Modeling of Sodium Ion and Water Intercalation into Titanium Disulfide Interlayers for Water Desalination

Recent experiments revealed the possibility of using titanium disulfide (TiS 2 ) as the cathode material in capacitive deionization (CDI) devices for water desalination. Although it performed stably up to 70 cycles with a salt removal capacity of 14 mg/g (corresponding to a sodium ion removal capacity of 35.8 mg/g) at high molar concentration (600 mM NaCl), the maximum capacity of TiS 2 as a CDI electrode was much lower (~70 mAh/g) than as a supercapacitor (239 mAh/g). Understanding why the ion capacities of these two configurations of the same material differ will entail elucidating detailed charge/discharge mechanisms at the atomic scale. Here, we present density functional theory simulations of sodium intercalation into TiS 2 interlayers in order to gain such understanding. We systematically investigated TiS 2 stacking patterns and ion intercalation sites, energetics of the intercalated compounds, and phase transformation during sodium intercalation. The calculated structural evolution and capacitance–voltage curve agree quite well with previous measurements. We conclude that the different maximum capacities of TiS 2 measured in aqueous and dry environments originate from weaker interlayer interactions with respect to shear strain after ~33% intercalation of Na + –H 2 O pairs, which is detrimental to the mechanical stability of TiS 2 . In conclusion, this study sheds light on the underlying mechanisms of ion intercalation into layered materials and contributes to understanding requirements for future design and optimization of CDI electrode materials for water desalination.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Facile Synthesis of Oxyhydrides by Reaction with NaBH 4 in an Open System

Oxyhydrides are an intriguing class of materials in which there is partial replacement of the oxide ion with hydride ions and oxygen vacancies. Conventional synthesis relies on vacuum sealed ampules, using long reaction times at high temperatures, limiting accessibility. Here, we demonstrate a rapid, ambient-pressure route to oxyhydride formation using NaBH 4 under flowing argon in just 1 h. This approach significantly lowers experimental barriers, enabling broader exploration of these materials. The maximum hydride incorporation, obtained using a reaction temperature of 400 °C, is given by the formula SrTiO 2.945 H 0.049 , where the hydride, oxygen vacancy, and unpaired electron concentrations are determined through thermogravimetric analysis, quantitative solid-state nuclear magnetic resonance (NMR) spectroscopy, and electron paramagnetic resonance (EPR) spectroscopy. Density functional theory simulations of the 1 H NMR shifts for candidate point defects support the assignment of the observed hydride peak, validating the efficacy of the synthetic approach. A combination of in situ and ex situ studies of the reaction pathway reveal that hydride incorporation into the perovskite occurs via direct solid-solid reaction, with higher reaction temperatures favoring NaBH 4 decomposition and H 2 (g) release over oxyhydride formation. The electronic defect structure established from the EPR and NMR studies indicate that, at ambient temperature, anion vacant sites are occupied by single electrons, whereas hydride sites do not trap electrons. As a result, this work establishes a scalable synthesis strategy and provides a computational-experimental framework for understanding defect chemistry in oxyhydrides, opening pathways for their integration into energy and electronic applications.

Anions↗
DOE OSTI2025

Structurally Driven Selective Adsorption of Hydrocarbons by Metal Substitution in Isostructural Rare-Earth Metal–Organic Frameworks

The design and realization of highly selective nanoporous materials are necessary to target critical separations across industries. By leveraging pore size, pore shape, and linker functionalization, the design of nanoporous solid adsorbents will enable the rapid production of energy efficient separation materials for high-value gas mixtures. This study uses a combination of modeling, synthesis, and gas adsorption testing to investigate a new class of small-pore isostructural rare-earth (RE) 2,5-dihydroxyterephthalic acid (DOBDC) metal–organic frameworks (MOFs) (RE: Pr-, Gd-, Er-, Yb; DOBDC = 2,5-dihydroxyterephthalic acid) and their adsorption selectivity for acetylene/ethylene mixtures. Density functional theory simulations identified that selective binding of acetylene over ethylene in the Gd-, Er-, and Yb-DOBDC MOFs was due to hydrogen-bonding between acetylene and the linker hydroxyl. Adsorption experiments validated the computational results by identifying mechanisms that control the acetylene/ethylene adsorption selectivity and high acetylene adsorption. Furthermore, dynamic column breakthrough experiments with the Gd-DOBDC MOF validated the simulations and indicated that ethylene can be separated from acetylene in a mixture containing 1 vol % acetylene and 39 vol % ethylene (balance argon). In conclusion, the results highlight the complexity of gas binding in functional porous materials and how combining modeling and experiment enables a fundamental understanding of gas–framework interactions that can be leveraged for the design of future separation materials.

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

First-Principles Study of Molecular Hydrogen Activation by Defects in Boron Nitride

Here, we used density functional theory simulations in combination with ab initio thermodynamics to determine the H 2 partial pressure (p H2 )-dependent energetics associated with H 2 activation and recovery at various defect sites in hexagonal boron nitride (h-BN). We found that some defects are very reactive with hydrogen, thereby definitely trapping hydrogen in defective h-BN. However, depending on hydrogen partial pressure, less reactive defect sites can be populated. Because of the lower binding capability of these sites, they would allow hydrogen to be recycled and recovered. For small defect sizes, we found that hydrogen preferentially binds to nitrogen sites by forming N–H bonds, and if no N sites are available then boron sites would be the next to bind hydrogen. Hydrogen dissociation via frustrated Lewis pair is found to be more favorable than forming only N–H bonds but only if the defect size is large enough to accommodate steric effects. For specific conditions such as T = 400 K, p H2 = 1 bar, and only considering one molecular H 2 per defect, three defects, namely, the N monovacancy, 3V(1B2N), and hexagonal 6V(3B3N) could play a role in both the activation and recycling of H 2 as they would be reacting enough to allow a favorable splitting of H 2 while not binding too strongly to allow its recovery. More broadly, a range of p H2 and hydrogen loading conditions were investigated for different types of defects and the finding suggests that p H2 could be used to fine-tune the Gibbs free energy of hydrogenation, thereby allowing several types of defects at different hydrogen loading contents to play a role in the activation/recovery process of H 2 in defective h-BN.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗