Search NASA⌕ Search

SEARCH · Search NASA

Results for “X-ray dynamical theory”

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

Real-Space Constrained Density Functional Theory Investigation of Site-Specific, Interfacial Charge Recombination Dynamics Across the Au Nanoparticle/TiO 2 Heterojunction

Au nanoparticle (NP)/TiO 2 heterojunction is a representative system to study interfacial charge transfer in photocatalysis and photovoltaics, where suppressing recombination from TiO 2 to Au can enhance hot carrier extraction. We apply real-space constrained density functional theory (CDFT) with Marcus theory to quantify charge recombination time scales across Au/TiO 2 . This approach enables direct control and visualization of charge-separated states, aligning with site-specific probes like time-resolved X-ray photoelectron spectroscopy (trXPS). We find that the charge-separated state features a bipolaron, with recombination dominated by TiO 2 LUMO to Au HOMO transitions, primarily at interfacial Au sites. Marcus rate predictions are benchmarked with surface hopping methods, quantifying differences in time scales and computational efficiency. Lastly, we examine how the Au cluster size affects the free energy change (ΔG) and reorganization energy (λ), explaining trends in closed-shell systems and highlighting challenges for open-shell extrapolations. Overall, CDFT + Marcus theory provides efficient, mechanistically transparent interfacial charge transfer modeling, and we clearly defined its applicability and limitation.

Glenna, Drew M. [Univ. of Idaho, Idaho Falls, ID (↗
DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

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

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2026

Pursuing Heteroleptic Ligand Design Principles for Photoactive Fe Complexes with Ultrafast X-ray Emission and Variable-Temperature Optical Spectroscopies

Understanding the key parameters that govern the photophysical and photochemical properties of transition metal complexes is essential for the development of efficient photosensitizers for photocatalytic applications. Achieving this objective necessitates clear and detailed investigations of their electronic excited states, for which time-resolved metal Kβ X-ray emission spectroscopy (XES) has proven highly effective. Here, we present a time-resolved Fe Kβ XES study of a heteroleptic Fe(II) polypyridyl carbene complex, [Fe(phen) 2 (C 4 H 10 N 4 )] 2+ (1; phen = 1,10-phenanthroline), utilizing both the valence-to-core and Kβ mainline spectral regions, complemented by variable-temperature transient optical absorption (VT-TA) spectroscopy. Detailed analysis of the time-resolved Kβ XES data, supported by density functional theory (DFT) calculations and an Eyring analysis of the VT-TA data, reveals parallel excited state relaxation dynamics that support an assignment of the long-lived excited state to a triplet metal-centered state. Placing these results in the context of prior studies of heteroleptic Fe(II) polypyridyl cyanide complexes motivated a series of DFT calculations to investigate the effects of ligand structural flexibility and arrangement. These calculations reinforce the experimentally derived conclusion that constraining structural flexibility with multidentate ligands significantly impacts the excited state relaxation dynamics. Furthermore, our study emphasizes that the arrangement of strong field ligands in heteroleptic complexes substantially affects the energy of Jahn–Teller active triplet metal-centered states in low-spin d 6 metal complexes. Together, these findings provide synthetic design principles for extending metal-to-ligand charge transfer excited state lifetimes of heteroleptic Fe complexes.

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

Signatures of Thermoreversible Associations in X-ray and Neutron Scattering from Dilute Polyzwitterion Solutions

In aqueous solutions of polyzwitterions (PZs), an interplay between dipole–dipole interactions and hydration of zwitterionic groups can lead to thermoreversible associations, which have been difficult to detect in experiments. Here, in this study, we investigated dilute aqueous solutions of poly(1-(3-sulphonatopropyl)-2-vinylpyridinium) (P2VPPS) using small-angle X-ray and neutron scattering (SAXS and SANS) to probe the structure and neutron spin-echo (NSE) spectroscopy to probe dynamics. The SAXS and SANS data show that the correlation length increases with an increase in the concentration of P2VPPS for three different molecular weights. The addition of 0.1 M NaCl to one of the solutions led to almost no dependence of the correlation length on the concentration. Such a concentration dependence of the correlation length suggests the formation of clusters driven by thermoreversible associations in the solutions. The NSE measurements show that the solutions with a larger correlation length display slower relaxation, reflecting the reduced mobility of larger clusters. These results should be considered as signatures of thermoreversible associations in the solutions of P2VPPS. To establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS, we combined a thermoreversible gelation theory for the structure of PZ solutions (Li, S.-F.; Muthukumar, M. Theory of Thermoreversible Gelation and Anomalous Concentration Fluctuations in Polyzwitterion Solutions. J. Chem. Phys. 2024, 161, 024903) with a model for the dynamics of the clusters (generalized Zimm model), developed in this work. Using such a theoretical framework, we have predicted the distribution of clusters in the solutions probed with SANS and SAXS. With the distributions, the generalized Zimm model has been used to extract the diffusion constant of the clusters and their characteristic size from the NSE data, where the latter agrees with the values estimated from the SAXS/SANS data. These findings confirm the presence of thermoreversible associations and establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS. Furthermore, with growing interest in technological applications, this work can provide useful insights into the structural and dynamical properties of other PZ solutions.

field theory↗
DOE OSTI2025

Digital Twin for Chemical Science: a case study on water interactions on the Ag(111) surface

Directly visualizing chemical trajectories offers insights into catalysis, gas-phase reactions and photoinduced dynamics. Tracking the transformation of chemical species is best achieved by coupling theory and experiment. Here we developed Digital Twin for Chemical Science (DTCS) v.01, which integrates theory, experiment and their bidirectional feedback loops into a unified platform for chemical characterization. DTCS addresses a core question: given a set of experimental conditions, what is the expected outcome and why? It consists of a forward solver that takes a chemical reaction network and predicts spectra under experimental conditions, and an inverse solver that infers kinetics from measured spectra. We applied DTCS to ambient-pressure X-ray photoelectron spectroscopy measurements of the Ag–H2O interface as an example. This approach enables real-time knowledge extraction and guides experiments until a stopping condition is met based on accuracy and degeneracy. As a step toward autonomous chemical characterization, DTCS provides mechanistic knowledge in a verified, standardized manner.

Chemistry↗
DOE OSTI2025

Exploring the Coexistence of Spin States in [Fe(tpy-Ph) 2 ] 2+ Complexes on Au(111) Using DFT Calculations

In this work, we systematically study the electronic structure and stability of spin states of the [Fe-(tpy-ph) 2 ] 2+ molecule in both the gas phase and on a Au(111) substrate using density functional theory + U (DFT+ U ) calculations. We find that the stability of the Fe 2+ ion’s spin states predicted by the computations is significantly influenced by the Hubbard U parameter. In the gas phase, the low-spin (LS, S = 0) state is found to be energetically favorable for U (Fe) ≤ 3 eV, whereas the high-spin (HS, S = 2) state is stabilized for U (Fe) > 3 eV. Interaction with the Au(111) substrate is found to elevate the critical U for the spin-state transition to 3.5 eV. Additionally, we perform L-edge X-ray absorption spectroscopy (XAS) calculations for both HS and LS states. The calculated XAS suggests that the HS state more closely aligns with the experimental observations, indicating the potential coexistence of the HS state as the initial state during the X-ray excitation process. These findings enrich our understanding of spin-state dynamics in [Fe(tpy-Ph) 2 ] 2+ .

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

Pressure‐Induced Structural and Magnetic Evolution in Layered Antiferromagnet YbMn 2 Sb 2

Electronic states under pressure exhibit unconventional spin and charge dynamics that provide a powerful route to uncover exotic phases in quantum materials. Here, we present the structural, magnetic, and electronic evolution of YbMn 2 Sb 2 under pressure. Single-crystal X-ray diffraction reveals a pressure-induced structural transition from the space group trigonal P$\overline{3}$m1 to the monoclinic P2 1 /m phase near 3.5 GPa, which remains stable up to 10 GPa. Magnetization measurements display an anomalously weak net magnetic moment and the absence of Curie–Weiss behavior up to 400 K, suggesting the formation of short-range Mn moment pairs that cancel macroscopically and subsequently evolve into long-range order upon cooling. Temperature-dependent resistivity shows semiconducting behavior with a transition at ∼119 K at ambient pressure, while pressure induces a dramatic suppression of resistance and the emergence of metallic-like temperature dependence, stabilized beyond 5 GPa. This pressure-driven semiconductor-metal transition is consistent with our density functional theory calculations, confirming the closing of the band gap under compression. Neutron diffraction under pressure identifies an incommensurate magnetic structure with antiparallel correlations between paired spins. Together, these results demonstrate how pressure-driven structural tuning and competing exchange interactions stabilize unconventional magnetic states in this low-dimensional magnetic semiconductor.

antiferromagnetic semiconductor↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Structural Propensities in Cs2MBiX6 (M=Na, Ag; X=Cl, Br) Bismuth Halide Double Perovskites

A previously unreported low-temperature phase transition in the bismuth halide double perovskite Cs2AgBiCl6 is reported, thereby establishing trends in the structural ground state across Cs2NaBiCl6, Cs2AgBiCl6, and Cs2AgBiBr6. Using the combined toolkit of variable-temperature synchrotron X-ray and neutron powder diffraction, Raman spectroscopy, and density-functional theory–based electronic structure modeling, we demonstrate a cubic Fm¯3m → tetragonal I4/m transition upon cooling with distinct onset temperatures. Neutron powder diffraction refinements permit the unambiguously assignment of the low-temperature phase of Cs2NaBiCl6 to I4/m, correcting prior reports of an I4/mmm ground state. Cs2AgBiCl6 is also found to transforms to a structure crystallizing in the I4/m space group at low temperatures. Temperaturedependent Raman data and density-functional theory-based modeling capture the softening and freezing of out-of-phase octahedral-tilt modes and quantify relative instabilities. Solid-state nuclear magnetic resonance spectroscopy at room temperature completes the characterization and helps underpin the subtle differences in covalency across the compounds. Trends in the phase transition temperature Ts and tilt magnitudes emerge from coupled effects of halide identity, M(I)–site bonding character, and a mismatch between interatomic distances. These results establish the structure– dynamics–bonding framework for tuning tilt-driven instabilities in halide double perovskites.

Tian, Haowen↗
DOE OSTI2025

Heavy boron doping effects on biaxially tensile strained germanium (>1.5%) investigated via structural characterization, effective lifetime assessment and atomistic modeling

Highly tensile strained germanium (ε-Ge) represents an essential material system for emerging electronic and photonics applications. Moreover, adjusting the doping levels to moderate or high concentrations can effectively tailor the properties of ε-Ge for specific applications. This article combines experimental characterization with a theoretical framework to examine the effects of heavy elemental boron (B) doping on pseudomorphic sub-50 nm ε-Ge. High resolution X-ray diffractometry is used to validate tensile strain levels of 1.53% and 1.68% in Ge epilayers, surpassing the indirect-to-direct band gap crossover point at ∼1.5% biaxial tensile strain. Cross-sectional transmission electron microscopy revealed visual evidence of stacking faults and surface roughening in 1.68% ε-Ge, although a coherent and abrupt Ge/III–V heterointerface is observed, devoid of interfacial misfit dislocations. Effective lifetime measurements demonstrated approximately twofold enhancement in 1.53% B-doped ε-Ge (N B ∼7 × 10 19 cm −3 ) compared to its unstrained B-doped counterpart, while no such improvement was observed in 1.68% B-doped ε-Ge. This lack of enhancement is attributed to the presence of stacking faults and surface roughness within the ε-Ge epilayer. Through density functional theory calculations, we independently demonstrate that substitutional B atoms induce local deformation of Ge–Ge bonds in both unstrained Ge and ε-Ge epilayers, resulting in an additive tensile strain. This phenomenon could potentially lead to dynamic reduction and overcoming of the critical layer thickness for the system, facilitating the nucleation and subsequent glide of 90° leading Shockley partial dislocations, thereby generating stacking faults. In essence, these findings establish an upper limit on the B-doping concentration that can be achieved in highly ε-Ge epilayers, and collectively, offer valuable insights into the significance of heavy doping in Ge-based heterostructures. As such, this study delineates a fundamental constraint for integrating heavily doped ε-Ge in high-performance optoelectronic systems, necessitating precise strain-doping co-optimization to avoid performance degradation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Solid solutions limited by grain-boundary solute clustering in ultrafine-grained alloys

Immiscible Ag-Cu alloys exhibit complex behavior due to varying Cu solid solubilities reported under equilibrium and metastable conditions. In ultrafine-grained alloys, these limits are further complicated by a high fraction of grain boundaries, where solute atoms tend to segregate and, in some cases, form clusters. Here, this study investigates the influence of Cu solute segregation and clustering on solid-solution limits in ultrafine-grained Ag-Cu alloys synthesized by magnetron sputtering with varying Cu content. X-ray diffraction peak shifts reveal a solid-solution concentration plateau for Cu contents from 4.9 to 11.7 at %, in contrast to the peak shifts predicted by density-functional theory for Ag-Cu alloys. Scanning transmission electron microscopy further reveals limited solid solubility and the formation of numerous Cu-rich clusters at grain boundaries. Atomistic simulations demonstrate that such limited solubility does not arise from grain boundary segregation alone, but only when strong solute-solute interactions promote the formation of grain-boundary Cu solute clusters.

Density-functional theory↗
DOE OSTI2025

LDRD conclusion poster - Synthesizing Heterometallic Uranium Single Crystals to Understand the Influence of the Secondary Metals on Uranyl Axial Bond Strength

Understanding how transition metals influenced the chemistry of lanthanide and actinide (f-element) materials is critical for advancing separation technologies, materials design, and coordination chemistry. This project examined how incorporating first-row transition metals affected the structural and spectroscopic properties of f-element coordination polymers. In uranium(VI)-based systems synthesized with 2,6-pyridinedicarboxylic acid (PDC) ligands, single-crystal X-ray diffraction and Raman spectroscopy revealed that the presence of transition metals shortened the uranyl axial bond and induced a blue shift in its symmetric stretching vibration—evidence of increased bond strength. Electronic structure analysis, including Density of States (DOS) calculations using density functional theory (DFT), revealed altered orbital overlaps and highlighted the role of transition metal d-orbitals in modulating bonding. Raman modes were modeled using truncated structural fragments in collaboration with the University of Notre Dame, and although the predicted frequencies were lower than experimental values, they remained within expected ranges. In parallel, similar experiments with cerium (Ce) in the presence of cobalt (Co) and PDC demonstrated multi-step single-crystal-to-single-crystal transformations—behavior not observed in the uranium systems. Initial products included light yellow, orange, and polycrystalline materials. Single-crystal X-ray diffraction studies, conducted in collaboration with the Colorado School of Mines, identified the yellow phase as monometallic Ce(PDC)2(H2O)2·4H2O and the orange phase as heterometallic Ce2Co(PDC)4(H2O)6. After standing in solution for one week, both phases fully transformed into a dark yellow crystalline phase, [Ce3(PDC)5(H2O)8].6(H2O). Remarkably, this transformation was reversible—disturbing the equilibrium by removing some crystals caused reversion to the initial Ce(PDC)2(H2O)2·4H2O phase, highlighting dynamic behavior. All three structures were previously unreported. Solid-state UV-visible and Raman spectroscopy further distinguished these phases, revealing ligand-to-metal charge transfer involving Ce and characteristic d–d transitions from Co(II). The precise mechanism driving these transformations remained unclear; however, pH-dependent experiments confirmed that the transformation did not occur when the pH decreased. Overall, the project demonstrated that transition metals could be employed to tune bonding interactions, structural dimensionality, and optical properties in f-element materials, establishing new pathways for designing functional heterometallic systems. The work resulted in several novel structural discoveries and fostered productive collaborations with the University of Notre Dame and the Colorado School of Mines.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2026

Structure of Complex Liquid–Liquid Extraction Organic Phases for Rare Earth Separations

Complex, multicomponent liquids with hierarchical structure and phase transitions are encountered in many natural and industrial processes, including in chemical separations. One notable example is aggregation and organic phase splitting in liquid–liquid extraction (LLE) of metal ions. While these two phenomena that have long been closely associated, a mechanistic link between mesoscale structure and the capacity-limiting organic phase splitting remains elusive due to complexity of these systems. Here, in this study, we combine small-angle X-ray scattering (SAXS), X-ray photon correlation spectroscopy (XPCS), and molecular dynamics simulation to reveal a comprehensive picture of structure at the nano- and mesoscale in these complex solutions. For the representative case of rare earth extraction from an acidic aqueous phase by a malonamide extractant in dodecane, we investigate a wide range of process-relevant extractant and acid concentrations to provide a complete picture of how aggregation depends on composition. We decompose organic phase structure from SAXS into two contributions, which together can capture the scattering at all compositions: composition fluctuations described by the Ornstein–Zernike equation at low wavenumber Q, and nanostructure modeled by a “pre-peak” at intermediate Q. The former contains information about the thermodynamics of demixing, while the latter reflects nanoscopic self-assembly of the extractant and extracted solutes. While fluctuations have typically not been considered in the literature, we find they in fact dominate the total structure for nearly all practical conditions. As only the fluctuations have a strong temperature response, we confirm this attribution with temperature-dependent SAXS measurements, including for extracted europium nitrate complexes. SAXS and XPCS measurements near the critical point find static and dynamic scaling consistent with theory. Overall, this new paradigm for understanding LLE organic phases connects composition, nanoscale, and mesoscale structuring to phase behavior, providing both a comprehensive picture of solution structure and a quantitative link between aggregation and third phase formation.

Peroutka, Allison A. [Argonne National Laboratory ↗
DOE OSTI2025

Pt Nanoparticle Disintegration at Oxide Interfaces Enhances CO Oxidation Catalysis

Understanding how supported metal nanoparticles dynamically evolve under reaction conditions is critical for controlling their catalytic function. Here, the mechanism behind the dynamic disintegration of Pt nanoparticles (NPs) supported on CeOx-TiO2 (CT) during CO oxidation is elucidated, leading to the formation of single atoms (SAs) and/or sub-nanometer clusters. Density functional theory (DFT) calculations reveal that strong Pt-CO interactions weaken Pt─Pt cohesion, while electronic coupling between Pt and Ce ions stabilizes Pt-CO* intermediates at the oxide interface. Surface oxygen vacancies kinetically trap Pt-CO*, but the vacancies are replenished under oxygen-rich conditions, enabling Pt-CO* surface diffusion and subsequent structural reorganization. In situ spectroscopic analyses confirm the oxygen-driven transformation of Pt NPs, correlating with a threefold increase in mass-specific activity at 150 °C. These findings demonstrate that interfacial oxygen dynamics and metal-support interactions can be leveraged to induce nanoparticle disintegration and optimize catalytic performance, highlighting the catalytic potential of interface-engineered Pt nanostructures.

CO oxidation↗
DOE OSTI2026

Machine learning-based interatomic potential development and phase transition analysis of ferroelectric hafnium dioxide

The ferroelectric phase (𝑃⁢𝑐⁢𝑎⁢2 1 , which is in orthorhombic symmetry) of hafnium dioxide (HfO 2 ) has gained much attention due to its potential applications in nanoelectronics and advanced memory devices. However, its complex phase behavior under external stimuli, such as pressure and temperature, remains a subject of intense investigation. This study focuses on developing a machine learning-based interatomic potential (MLIP) that is trained with data from density-functional theory (DFT) calculations to simulate phase transitions and mechanical properties of HfO 2 . The developed MLIP predicts lattice parameters, equations of state, bulk and shear moduli, and elastic constants that closely align with DFT predictions for several phases and at various pressures. Once validated, the MLIP is used to investigate the phase transitions of ferroelectric HfO 2 (𝑃⁢𝑐⁢𝑎⁢2 1 ) under both isobaric and constant stress conditions at elevated temperatures ranging from 200 to 2500 K. We used several complementary methods, including local symmetry identification, radial distribution function, and x-ray diffraction characterization, to identify interesting phase transitions among several competitive hafnia phases predicted from our simulations. The suggested methods uniformly reveal that under pure deviatoric condition, the system favors a transition from the orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase to a tetragonal (𝑃⁢4 2 /𝑛⁢𝑚⁢𝑐) phase, whereas a zero stress condition drives the system from the 𝑃⁢𝑐⁢𝑎⁢2 1 phase to another orthorhombic (𝑃⁢𝑏⁢𝑐⁢𝑛) phase. These findings provide crucial insights into stress and temperature-induced phase behavior of hafnia, guiding future experimental and theoretical studies for optimizing hafnia-based ferroelectric devices.

Ferroelectric HfO2↗
DOE OSTI2025

Barrierless nucleation in glassy precursors drives zeolite formation

Zeolites are crystalline, microporous silicates widely used in catalysis and separations, yet the molecular mechanisms of their formation remain unresolved. Experiments indicate that hydrothermal synthesis of silica zeolites from clear solution proceeds through amorphous nanoaggregates that gradually develop zeolite order in an apparently continuous amorphous to crystal transformation. Here, we combine molecular simulations with advanced algorithms that identify zeolite order and computer vision to elucidate the pathway from clear solution to zeolite nanocrystal. We show that at conditions of hydrothermal synthesis of silica zeolites, the transformation of precursor aggregate into zeolite is not limited by nucleation barriers but by the slow dynamics of reorganization in the glassy precursor matrix. The negligible nucleation barriers result in spinodal-like crystallization that leads to a gradual formation of a mosaic of small crystallites that explain the seemingly continuous character of zeolite crystallization and the catalytic activity of X-ray amorphous, protozeolites and embryonic zeolites. We find that zeolite-like porosity and short-range order emerge early within glassy precursors, well before crystallinity is detected in transmission electron microscopy (TEM) images or X-ray diffraction. The nanoaggregate’s temperature-size phase diagram reveals a convergence of the zeolite–amorphous equilibrium and maximum crystallization rate at ~3 nm diameter nanoparticle diameters and ~200 °C. This convergence signals the termination of the first-order amorphous-to-zeolite transition. Our results provide a unifying framework for understanding nucleation of silica zeolites from solution and suggest that barrierless nucleation may govern the formation of other nanoparticle systems, including minerals and oxides synthesized far below their bulk melting points.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic Features of Cu-Ceria Interface under CO 2 Hydrogenation to Methanol

It is generally accepted that metal–support interaction is very important for the hydrogenation of CO 2 to methanol, but little has been revealed about the feature of interfacial active sites under real reaction conditions since there are only limited techniques that can be applied under high-pressure conditions. Here, in this work, by combining multiple in situ and operando techniques on a model Cu/ceria catalyst, we have tracked Cu and ceria sites for methanol formation. Under the reaction condition, it is found that upon reaching the reaction temperature, oxidized Cu species in the as-synthesized catalyst immediately change into metallic Cu species. Following this, it is the gradual formation of methanol, the changing rate of which coincides with the formation of a unique Ce 3+ species. The combined experimental results and density functional theory (DFT) calculations have determined that the formed Ce 3+ sites driven by the reaction conditions are bound to hydrides, adsorbed carbonate species, and interfacial active Cu sites. The Cu-ceria interaction in this complex moiety is weak and can be easily disturbed with reaction environment variations, leading to dynamic changes at the interface upon the hydrogenation of active carbonate intermediates, which are precursors for the formation of methanol. The formation of this unique Cu–Ce 3+ interface and its dynamicity lead to an increase of methanol selectivity from less than 20% to 60%. These results suggest that reactant-derived species (H – and carbonate in this work) can be essential components of the active center with the functions of manipulating the metal−oxide interaction and directing reaction pathways.

CO2 hydrogenation to methanol↗