Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory modeling”

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

Fast Proton Insertion in Layered H 2 W 2 O 7 via Selective Etching of an Aurivillius Phase

H 2 W 2 O 7 , a metastable material synthesized via selective etching of the Aurivillius-related Bi 2 W 2 O 9 , is shown as an electrode for high power proton-based energy storage. Comprehensive structural characterization is performed to obtain a high-fidelity crystal structure of H 2 W 2 O 7 using an iterative approach that combines X-ray diffraction, neutron pair distribution function, scanning transmission electron microscopy, Raman spectroscopy, and density functional theory modeling. Electrochemical characterization shows a capacity retention of ≈80% at 1000 mV s –1 (1.5-s charge/discharge time) as compared to 1 mV s –1 (≈16-min charge/discharge time) with cyclability for over 100 000 cycles. Energetics from density functional theory calculations indicate that proton storage occurs at the terminal oxygen sites within the hydrated interlayer. Last, optical micrographs collected during in situ Raman spectroscopy show reversible, multicolor electrochromism, with color changes from pale yellow to blue, purple, and last, orange as a function of proton content. These results highlight the use of selective etching of layered perovskites for the synthesis of metastable transition metal oxide materials and the use of H 2 W 2 O 7 as an anode material for proton-based energy storage or electrochromic applications.

25 ENERGY STORAGE↗

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

Modeling neutral defects in III-V ternary alloys with a special quasirandom structure: Analysis of As- and III-site point defects in InGaAs

While first-principles density functional theory modeling has become a vital tool to investigate defect properties in semiconductors, the lack of crystalline periodicity in pseudobinary random composition alloys, such as In 1−𝑥 ⁢Ga 𝑥 ⁢As, complicates such analyses. We present a simulation strategy to systematically take into account the variability in the local defect environment in order to predict statistical properties of neutral intrinsic defects in In 1−𝑥⁢ Ga 𝑥 ⁢As. We use a comprehensive sampling from a modest-sized 64-atom special quasirandom structure (SQS) to define a statistically representative set of defects, and use a 512-atom hypercell, a 2 × 2 × 2 supercell of SQS supercells, to achieve cell-size convergence. We articulate an equivalent site principle and describe how it constrains atomic chemical reference energies in computation of defect formation energies in pseudobinary alloys. A simple protocol for estimating reference energies for the Ga and In atoms sharing the III site succeeds in obtaining the equivalence of defects at Ga-sites and In sites in the SQS supercell, (<30 meV differences in average formation energies). For III-site defects, such as the As antisite As III , the statistical variability in formation energies is modest, ≈ 0.1–0.2 eV. The variability in formation energy at As-site defects, such as the As vacancy 𝑣 As , can be much larger, >1 eV. The As antisite is shown to be a low-energy defect and the most likely to be present in as-grown materials, just as in GaAs. All other defects are higher-energy defects unlikely to be important in native material, but potentially important in radiation-damaged material. With a strong variability in defect energies, especially on the As-site, explicit consideration of statistical variability due to compositional randomness will be imperative for meaningful and quantitative comparisons to experiment.

Density functional theory↗

Understanding the Role of CdTe in Polycrystalline CdSe x Te 1– x /CdTe-Graded Bilayer Photovoltaic Devices

Grading of bandgap by alloying CdTe with selenium to form a CdSe x Te 1–x /CdTe-graded bilayer device has led to a device efficiency over 19%. A CdSe x Te 1–x absorber would increase the short-circuit current due to its lower bandgap but at the expense of open-circuit voltage. It has been demonstrated that adding a CdTe layer at the back of such a CdSe x Te 1–x film reduces the voltage deficit caused by the lower bandgap of absorber from selenium alloying while maintaining the higher short-circuit current. This leads to a photovoltaic device that draws advantage from both materials with an efficiency greater than either of them. In this report there is a detailed account using device data, ultraviolet photoelectron spectroscopy, electron microscopy, and first-principles density functional theory modeling is provided, which shows that CdTe acts as an electron reflector for CdSe x Te 1–x .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic Investigation of Molybdenum Disulfide Defect Photoluminescence Quenching by Adsorbed Metallophthalocyanines

Lattice defects play an important role in determining the optical and electrical properties of monolayer semiconductors such as MoS 2 . Although the structures of various defects in monolayer MoS 2 are well studied, little is known about the nature of the fluorescent defect species and their interaction with molecular adsorbates. In this study, the quenching of the low-temperature defect photoluminescence (PL) in MoS 2 is investigated following the deposition of metallophthalocyanines (MPcs). The quenching is found to significantly depend on the identity of the phthalocyanine metal, with the quenching efficiency decreasing in the order CoPc > CuPc > ZnPc, and almost no quenching by metal-free H 2 Pc is observed. Futher, time-correlated single photon counting (TCSPC) measurements corroborate the observed trend, indicating a decrease in the defect PL lifetime upon MPc adsorption, and the gate voltage-dependent PL reveals the suppression of the defect emission even at large Fermi level shifts. Density functional theory modeling argues that the MPc complexes stabilize dark negatively charged defects over luminescent neutral defects through an electrostatic local gating effect. These results demonstrate the control of defectbased excited-state decay pathways via molecular electronic structure tuning, which has broad implications for the design of mixeddimensional optoelectronic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding surfaces and interfaces in nanocomposites of silicone and barium titanate through experiments and modeling

Barium titanate (BTO) is a ferroelectric perovskite used in electronics and energy storage systems because of its high dielectric constant. Decreasing the BTO particle size was shown to increase the dielectric constant of the perovskite, which is an intriguing but contested result. We investigated this result by fabricating silicone-matrix nanocomposite specimens containing BTO particles of decreasing diameter. Furthermore, density functional theory modeling was used to understand the interactions at the BTO particle surface. Combining results from experiments and modeling indicated that polymer type, particle surface interactions, and particle surface structure can influence the dielectric properties of polymer-matrix nanocomposites containing BTO.

composite↗

An Analytical Tool to Evaluate Defect Thermodynamics of (La,Ba)Fe1-xMxO3-δ Perovskites for Solid-Oxide Cell Applications

A modeling tool of the defect thermodynamics of (La,Ba)Fe1-xMxO3-δ perovskites which includes energetic information about oxygen vacancy formation, hydration, hydride formation, and charge disproportionation reactions has been developed1. This tool incorporates defect energies and entropies expressed as sixth order polynomial functions to allow refinements of the defect reaction equilibrium constants in the thermodynamic analysis. Calculation of (La,Ba)Fe1-xMxO3-δ Brouwer diagrams as a function of pO2/pH2O and pH2/H2O in a range of temperatures of interest is facilitated by this modeling tool. The results obtained can provide direct guidance how the electronic and ionic defect concentrations of the triple conducting perovskite materials can be used to optimize performance of solid oxide cells for energy applications. The impact of magnetic and electronic structures of the perovskites on the defect reaction energies and entropies as obtained from density function theory modeling and the role played by hydride defect species will also be discussed.

Lee, Yueh-Lin↗

Explaining the structure sensitivity of Pt and Rh for aqueous-phase hydrogenation of phenol

Phenol is an important model compound to understand the thermocatalytic (TCH) and electrocatalytic hydrogenation (ECH) of biomass to biofuels. Although Pt and Rh are among the most studied catalysts for aqueous-phase phenol hydrogenation, the reason why certain facets are active for ECH and TCH is not fully understood. Herein, we identify the active facet of Pt and Rh catalysts for aqueous-phase hydrogenation of phenol and explain the origin of the size-dependent activity trends of Pt and Rh nanoparticles. Phenol adsorption energies extracted on the active sites of Pt and Rh nanoparticles on carbon by fitting kinetic data show that the active sites adsorb phenol weakly. We predict that the turnover frequencies (TOFs) for the hydrogenation of phenol to cyclohexanone on Pt(111) and Rh(111) terraces are higher than those on (221) stepped facets based on density functional theory modeling and mean-field microkinetic simulations. The higher activities of the (111) terraces are due to lower activation energies and weaker phenol adsorption, preventing high coverages of phenol from inhibiting hydrogen adsorption. We measure that the TOF for ECH of phenol increases as the Rh nanoparticle diameter increases from 2 to 10 nm at 298 K and -0.1 V vs the reversible hydrogen electrode, qualitatively matching prior reports for Pt nanoparticles. The increase in experimental TOFs as Pt and Rh nanoparticle diameters increase is due to a larger fraction of terraces on larger particles. These findings clarify the structure sensitivity and active site of Pt and Rh for the hydrogenation of phenol and will inform the catalyst design for the hydrogenation of bio-oils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic insights into nitrogen activation on atomic Ru clusters in self-pillared pentasil using operando atomistic models and experimental kinetics

Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.

36 MATERIALS SCIENCE↗

Constraints on triple oxygen isotope kinetics

The general assumption that systems are in equilibrium is convenient, but imperfect. Instead, the partial or full expression of kinetic reaction steps may need to be determined. Using an isotope fractionation Monte-Carlo model under the harmonic approximation (Figure 1 envelopes), we find the range of kinetic isotope effects (i.e., resulting from unidirectional chemical reactions) contains and exceeds the range of accessible equilibrium isotope effects. For the 16O-17O-18O system, accessible kinetic isotope effects can yield Δ′17O variations as large 0.8‰ from a single kinetic fractionation step. Reactions that yield the largest deviations from equilibrium are those with the largest molecular mass difference between decomposition fragments, typically the two products, such as a gas and solid generated from the thermal decomposition of an initial solid. For example, previously published experiments on the thermal decomposition of calcite (CaCO3) and dehydroxylation of brucite (Mg(OH)2) are matched by density functional theory models of each kinetic isotope effect, thus revealing their origin in mass-dependent reaction. Beyond oxygen, our approach applies to any system with more than two isotopes.

Oxygen isotopes↗

Aqueous-phase effects on ethanol decomposition over Ru-based catalysts

The effects of an aqueous phase on ethanol decomposition for hydrogen production over a Ru(0001) catalyst surface model have been investigated from first principles. Solvent effects on the reaction mechanism and kinetic parameters have been quantified with the help of a microkinetic reactor model, density functional theory, and an implicit solvation scheme (iSMS). Our calculations indicate that in both vapor- and aqueous-phase reaction environments, the ethanol decomposition starts with acetaldehyde formation on the surface, some of which further dehydrogenates to ketenyl species (CHCO), where the C–C bond cleaves to form methylidyne (CH) and CO. In the vapor phase, adsorbed CH gets hydrogenated to methane, and CO desorbs or undergoes methanation reducing the amount of hydrogen produced. In contrast, under aqueous phase reaction conditions, the methanation is inhibited, and the water–gas shift (WGS) reaction is accelerated, leading to complete conversion of CO to CO 2 and H 2 . Furthermore, calculations indicate that the observed reaction behavior under aqueous phase reforming conditions originates primarily from the higher water chemical potential, and implicit solvent models predict only a small solvation effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atomically Revealing Bulk Point Defect Dynamics in Hydrogen‐Driven γ‐Fe 2 O 3 → Fe 3 O 4 → FeO Transformation

Understanding how point defects in the bulk govern redox transformations is essential for advancing hydrogen-based metal production and designing high-performance oxide materials. This study reveals the atomic-scale mechanisms driving hydrogen-induced reduction of γ-Fe 2 O 3 to Fe 3 O 4 , focusing on how bulk vacancy dynamics dictate structural evolution and reaction kinetics. A key finding is the pronounced contrast in defect behavior between the two oxides: in γ-Fe 2 O 3 , intrinsic Fe vacancies promote oxygen vacancy clustering, destabilizing the local lattice and driving nanopore formation. In contrast, Fe 3 O 4 exhibits a higher oxygen vacancy formation energy and lacks intrinsic Fe vacancies, suppressing vacancy aggregation and maintaining a dense, pore-free structure. This divergence governs distinct reduction pathways—γ-Fe 2 O 3 undergoes an interface-reaction-limited transformation confined to the γ-Fe 2 O 3 /Fe 3 O 4 boundary, while Fe 3 O 4 supports a uniform increase in oxygen vacancy concentration, enabling bulk-phase reduction to lower-oxide FeO. Integrated in situ electron microscopy and density functional theory modeling uncover a vacancy-mediated mechanism, where synergistic cation-anion vacancy dynamics steer microstructure evolution and phase progression. These insights highlight the critical role of vacancy dynamics in controlling oxide reactivity and offer a pathway toward vacancy engineering to enhance reduction kinetics in hydrogen metallurgy and to tailor porosity, reactivity, and structural resilience in oxide-based catalysts and energy materials.

36 MATERIALS SCIENCE↗

Vapor Phase Installation of CpCo(CO) 2 in MOF‐808

Metal-organic frameworks (MOFs), including MOF-808 present an opportunity for vapor phase installation of reactive centers that may serve as uniform sites for precision catalysis. The mechanism of cobalt installation in MOF-808 through atomic layer deposition (ALD) upon CpCo(CO )2 exposure is investigated through in situ FTIR spectroscopy complemented by density functional theory modeling. The role of subsequent H 2 O exposure in the hydrolysis of carbonyl and cyclopentadienyl ligands to complete the installation is also investigated. In situ FTIR study reveals that upon exposure of CpCo(CO) 2 to MOF-808 at 115 °C, a long-lived, stable intermediate is formed that presents carbonyl stretching vibrations similar to the undissociated precursor that are attributed to multiple chemisorbed carbonyl complexes. In MOF-808, cobalt species are initially absorbed near the BTC linker rather than directly on the nodes, revealing the importance of non-covalent interactions in the installation process. DFT suggests that H 2 O exposure promotes carbonyl elimination but leaves a cyclopentadienyl-capped cobalt that is more stable than the complete hydrolysis product. The simulations are consistent with cobalt installation in MOF-808 using CpCo(CO) 2 and H 2 O at 115 °C and further consistent with the lack of ALD thin film deposition of cobalt oxide thin films on planar supports.

36 MATERIALS SCIENCE↗

Mechanisms and kinetics of citrate-promoted dissolution of a uranyl phosphate mineral

The formation of uranyl phosphate precipitates is a remediation strategy to limit the mobility of uranium in contaminated soils. However, exposure to organic ligands, like the plant exudate citrate, can remobilize the uranium. The purpose of this study is to provide a more thorough comprehension of citrate-promoted dissolution of a uranyl phosphate mineral, chernikovite [(H 3 O)(UO 2 )(PO 4 )·3H 2 O], by determining the extent of uranium release from chernikovite at a wide range of citrate concentrations. Here we have quantified the kinetics of dissolution and proposed potential mechanisms of chernikovite dissolution to gain a better understanding of the fate of uranyl phosphate precipitates in the environment. Batch dissolution and continuously stirred tank reactor (CSTR) experiments indicate that increasing citrate concentrations from 0.1 mM to 50 mM increases the concentration of dissolved uranium. However, at citrate concentrations of 10 mM and greater, the effectiveness of the ligand to enhance the dissolution decreases. Density functional theory modeling and Raman spectroscopy indicate a strong interaction between citrate and the uranyl ions at the surface of the uranyl phosphate. The interaction between citrate and uranyl suggests the formation of an alteration layer of uranyl-citrate surface complexes, which could impede uranium release at citrate concentrations of 10 mM or higher. Together, these results indicate that the citrate-promoted dissolution of chernikovite is a fast reaction that is hindered by a combination of surface saturation and a secondary-phase precipitation reaction at higher concentrations of citrate and reaction times greater than or equal to 72 min based on the hydraulic residence time of CSTR experiments. These findings improve our understanding of the stability of uranyl phosphate for use in sequestration of uranium groundwater plumes at contaminated sites.

58 GEOSCIENCES↗

The u 3 si 2 -H system

U 3 Si 2 is of interest to the nuclear industry as a candidate fuel material due to its high uranium density and high thermal conductivity. However, it has been observed to react with hydrogen, resulting in material decrepitation. As a result, it is important to understand the thermodynamics of the U 3 Si 2 -H system. In this study, the thermodynamics of the hydrogen absorption reaction of USi were determined experimentally using Sievert’s gas absorption and related to crystallographic evolution with hydrogen content using X-ray diffraction. Experimentally-determined thermodynamic parameters were compared with results from density functional theory modeling. Results from this study were also compared with those determined in previous work. Sievert’s gas absorption results were used to develop the pressure-composition-temperature (PCT) curves of the U 3 Si 2 -H system. It was found that the hydride phase exhibited a maximum stoichiometry between U 3 Si 2 H 1.8 and U 3 Si 2 H 2 . The two-phase region for hydride formation from U 3 Si 2 exhibited a miscibility gap with a critical temperature between 623 and 673 K, as calculated from the PCT curves. Analysis of the PCT curves also showed that both the enthalpy and entropy of the hydrogen absorption reaction increased with hydrogen content but were lower than the values for uranium trihydride formation from uranium metal. The enthalpy of reaction for hydrogen absorption was calculated to range between -86.9 and -94.8 kJ mol -1 , while the entropy of reaction was calculated to range between 101.9 and 138.8 J mol -1 K -1 . Furthermore, DFT modeling of the thermoydnamic stability of the USi hydride phases yielded a decomposition temperature of U 3 Si 2 H 2 of approximately 600 K, which was consistent with the experimental results. Similarly, the DFT-calculated enthalpy and entropy of reaction to form USiH were determined to be -106.5kJ mol and 121.8J mol -1 K -1 , respectively, which were both in close agreement with the experimentally-determined values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unlocking enhanced gas capture via core scrambling of porous-organic cages

The demand for low-cost, low-energy, and highly selective gas capture and separations is an ongoing driver of porous material development. Porous liquids have been identified as a promising gas separation material by creating permanent porosity in inorganic solvents through inclusion of nanoporous materials that sterically exclude solvent from their internal porosity. Among the nanoporous materials that can be used to form porous liquids, porous-organic cages (POCs) have been one of the most popular due to the inherent tunability of POCs. “Scrambled” POCs with varying functionalities on the POC vertices have been developed and incorporated into porous liquid compositions, increasing their gas adsorption capacity. An unexplored avenue to tailor the properties of porous liquids is through scrambling the functionality of the core of the POC. Here, therefore, we have synthesized a new POC, a CC3-OH derivative with scrambled hydroxides on the core and evaluated the impact on the CO 2 uptake capacity in silicon oil-based porous liquids. Core scrambling of the POC resulted in a twofold increase CO 2 adsorption capacity in the porous liquid, an emergent property that is a dramatic increase beyond a linear combination of the gas adsorption capacity of the neat solvent and the POC. Density functional theory modeling of the CC3 POC and its hydroxide-based derivatives identified that free rotation of the linker hydroxide allowed for forced interaction between the CO 2 molecule and the hydroxide in the pore window. Solvation of the POC may release scrambled core hydroxides from intramolecular bonding with a neighboring imine, allowing for increased gas uptake in the porous liquid over the neat POC. These results identify a key structural relationship of POCs that enables emergent properties in porous liquids and can guide future development of liquid phase gas capture and separation materials for environmental and industrial applications.

Gas capture↗

Atomic-Scale View of Redox Induced Changes for Monolayer MoO x on α-TiO 2 (110) with Chemical-State Sensitivity

Supported molybdenum oxide (MoO x ) plays an important role in catalytic transformations from alcohol dehydrogenation to transesterification. During these reactions, Mo and oxygen surface species undergo structural and chemical changes. A detailed, chemical-state specific, atomic-scale structural analysis of the catalyst under redox conditions is important for improving catalytic properties. In this study, monolayer of MoO x grown on α-TiO 2 (110) by atomic-layer deposition is analyzed by X-ray standing wave (XSW) excited X-ray photoelectron spectroscopy (XPS). The chemical shifts for Mo 2p 3/2 and O 1s peaks are used to distinguish Mo 6+ from Mo 4+ and surface O from bulk O. Excitation of XPS by XSW allows pin-pointing the location of these surface species relative to the underlying substrate lattice. Finally, measured 3D composite atomic density maps for the oxidized and reduced interfaces compare well with our density functional theory models, and collectively create a unique view of the redoxdriven dynamics for this complex catalytic structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗