Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic & molecular clusters”

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

Domain Separated Density Functional Theory for Reaction Energy Barriers and Optical Excitations

We recently proposed Domain Separated Density Functional Theory (DS-DFT), a framework that allows for the combination of different levels of theory for the computation of the electronic structure of molecules. Here, this work discusses the application of DS-DFT to the computation of transition-state energy barriers and optical absorption spectra. We considered several hydrogen abstraction reactions and optical spectra of molecule/metal cluster systems, including the absorption of individual species such as carbon monoxide, methane, and molecular hydrogen to a Li6 cluster. We present and discuss two domain-separated methods: i), the screened-density approximation (SDA), and, ii), linearly weighted exchange (LWE). We find that SDA, which is applied as a hybridization based on atomic domains, could be useful to computing energy barriers, whereas LWE is suited for the analysis of electronic properties such as ground-state gaps, excitation energies, and oscillator strengths.

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

Immobilization of “Capping Arene” Cobalt(II) Complexes on Ordered Mesoporous Carbon for Electrocatalytic Water Oxidation

We report the synthesis, characterization, and electrocatalytic water oxidation activity of two cobalt complexes, (6-FP)Co(NO 3 ) 2 (1) (6-FP = 8,8'-(1,2-phenylene)diquinoline) and (5-FP)Co(NO 3 ) 2 (2) (5-FP = 1,2-bis(N-7-azaindolyl)benzene), containing "capping arene" bidentate ligands with nitrogen atom donors. The cobalt complexes 1 and 2 were supported on ordered mesoporous carbon (OMC) by π-π stacking, resulting in heterogenized cobalt materials 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, respectively, and studied for electrocatalytic water oxidation. We find that 6-FP-Co-OMC-1 exhibits an overpotential of 355 mV for a current density of 10 mA cm -2 and a turnover frequency (TOF) of ~0.53 s -1 at an overpotential of 400 mV at pH 14. 6-FP-Co-OMC-1 exhibits activity that is ~1.6 times that of 5-FP-Co-OMC-2, which gives a TOF of 0.32 s -1 at 400 mV overpotential. The structural stability of the single-atom Co site was demonstrated for 6-FP-Co-OMC-1 using X-ray absorption spectroscopy for the molecular complex supported on OMC, but slow degradation in catalyst activity can be attributed to eventual formation of Co oxide clusters. DFT computations of electrocatalytic water oxidation using the molecular complexes as models provide a description of the catalytic mechanism. These studies reveal that the mechanism for O-O bond formation involves an intermediate Co IV oxo complex that undergoes an intramolecular reductive O-O coupling to form a Co II -OOH species. Further, the calculations predict that the molecular 6-FP-Co structure is more active for electrocatalytic water oxidation than 5-FP-Co, which is consistent with experimental studies of 6-FP-Co-OMC-1 and 5-FP-Co-OMC-2, highlighting the possibility that the ligand structure influences the catalytic activity of the supported molecular catalysts.

25 ENERGY STORAGE↗
DOE OSTI2021

Exploring electron beam induced atomic assembly via reinforcement learning in a molecular dynamics environment

We report atom-by-atom assembly of functional materials and devices is perceived as one of the ultimate targets of nanotechnology. Recently it has been shown that the beam of a scanning transmission electron microscope can be used for targeted manipulation of individual atoms. However, the process is highly dynamic in nature rendering control difficult. One possible solution is to instead train artificial agents to perform the atomic manipulation in an automated manner without need for human intervention. As a first step to realizing this goal, we explore how artificial agents can be trained for atomic manipulation in a simplified molecular dynamics environment of graphene with Si dopants, using reinforcement learning. We find that it is possible to engineer the reward function of the agent in such a way as to encourage formation of local clusters of dopants under different constraints. This study shows the potential for reinforcement learning in nanoscale fabrication, and crucially, that the dynamics learned by agents encode specific elements of important physics that can be learned.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Molecular dynamics study of interstitial He clusters in nickel

This study presents a molecular dynamics analysis focusing on the behavior of interstitial helium (He) clusters in nickel (Ni), examining their formation, stability, and migration energetics. Consistent with previous research, we found that the binding energies of interstitial helium within a helium cluster are positive and increase with the cluster size, indicating a preference for helium atoms to cluster together. However, our findings also reveal that while the formation energy increases monotonically with cluster size, the increase in binding energy is non-monotonic. Importantly, small He clusters were observed to be thermally unstable at reactor operational temperatures (approximately 600 K), with the He 2 cluster exhibiting instability even at room temperature. With a binding energy of 0.49 eV for a He 4 cluster, we hypothesize that for helium bubbles to form via homogeneous nucleation (i.e., through trap mutation) at reactor operating temperatures, the helium concentration must be high enough to facilitate the formation of helium clusters of at least size 4 or larger. As expected, interstitial helium and small helium clusters are highly mobile. This mobility was observed not only at room temperature but also at temperatures as low as approximately 200 K. Furthermore, the mean squared displacement method has been utilized to determine the migration barriers and the corresponding prefactors for clusters ranging from He 1 to He 6

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The Reactivity and Structural Dynamics of Supported Metal Nanoclusters Using Electron Microscopy, in-situ X-Ray Spectroscopy, Electronic Structure Theories, and Molecular Dynamics Simulations

This research will provide new methods to determine the structures and compositions of heterogeneous catalysts as they exist in reaction conditions. The new tools will make it possible to understand how extremely small material structures—ones in a size range spanning from single atoms to larger polyatomic clusters—make it possible to carry out highly selective chemical transformations required to synthesize useful materials as well as for efficient energy production and utilization. The tools will make it possible to see how the atomic structures of catalytic materials evolve over time under the conditions of their use, and in this way provide understandings to guide the design of new energy-efficient catalysts. The work will make it possible to answer long standing questions that include the following. How does catalyst structure help guide molecules along specific reaction pathways? What are the important intermediates involved in these transformations? How do the structures and reactivity of catalytic materials, as well as the chemistry of their interfaces, vary with reaction conditions? Is it possible to develop new catalytic materials using new knowledge of the active structures and reaction mechanisms? The work will provide heretofore unavailable experimental tools that can be applied in advanced programs of research addressing such questions.

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

Structural adaptation of oxygen tolerance in 4-hydroxybutyrl-CoA dehydratase, a key enzyme of archaeal carbon fixation

Abstract Autotrophic microorganisms that convert inorganic carbon into organic matter were key players in the evolution of life on Earth. As the early atmosphere became oxygenated, microorganisms needed to develop mechanisms for oxygen protection, especially those relying on enzymes containing oxygen-sensitive metal clusters (e.g., Fe-S). Here we investigated how 4-hydroxybutyryl-CoA dehydratase (4HBD) - the key enzyme of the 3-hydroxypropionate/4-hydroxybutyrate (HP/HB) cycle for CO 2 -fixation - adapted as conditions shifted from anoxic to oxic. 4HBD is found in both anaerobic bacteria and aerobic ammonia-oxidizing archaea (AOA). The oxygen-sensitive bacterial 4HBD and oxygen-tolerant archaeal 4HBD share 59 % amino acid identity. To examine the structural basis of oxygen tolerance in archaeal 4HBD, we determined the atomic resolution structure of the enzyme. Two tunnels providing access to the canonical [4Fe-4S] cluster in oxygen-sensitive bacterial 4HBD were closed with four conserved mutations found in all aerobic AOA and other archaea. Further biochemical experiments and molecular dynamics simulations support our findings that restricting access to the active site is the key to oxygen tolerance, explaining how active site evolution drove a major evolutionary transition. Significance statement Autotrophy (primary production) was the first life strategy on Earth. Before photosynthesis (using solar energy to fix carbon dioxide), life relied on chemical reactions for energy. These chemosynthetic reactions are present in all domains of life, including archaea possessing the most energy-efficient carbon fixation pathway - the 3-hydroxypropionate/4-hydroxybutyrate cycle. This efficiency results from enzyme modifications, including enhanced enzyme stability and catalysis of multiple reactions. We reveal the first structure of aerobic 4-hydroxybutyryl-CoA dehydratase (4HBD) from ammonia-oxidizing archaea. These archaea are among the most abundant organisms on the planet, and their 4HBD active site evolved oxygen tolerance to support aerobic metabolism. This modification can provide further insight into enzyme evolution on early earth, as photosynthesis developed and began oxygenating the atmosphere.

DeMirci, Hasan (ORCID:0000000291355397)↗
DOE OSTI2020

Materials Data on H5ClO6 by Materials Project

H5O2(H2O)4(HClO4)2H5O6ClClO4 crystallizes in the triclinic P1 space group. The structure is zero-dimensional and consists of one molecular hydrogen;dihydrate molecule, four water molecules, one ClO4 cluster, one H5O6Cl cluster, and two HClO4 clusters. In the ClO4 cluster, there are four inequivalent O sites. In the first O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.47 Å. In the second O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.47 Å. In the third O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.47 Å. In the fourth O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. Cl is bonded in a tetrahedral geometry to four O atoms. In the H5O6Cl cluster, there are five inequivalent H sites. In the first H site, H is bonded in a single-bond geometry to two O atoms. There is one shorter (1.00 Å) and one longer (1.71 Å) H–O bond length. In the second H site, H is bonded in a linear geometry to two O atoms. There is one shorter (1.20 Å) and one longer (1.24 Å) H–O bond length. In the third H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 0.99 Å. In the fourth H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 0.99 Å. In the fifth H site, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 0.98 Å. There are six inequivalent O sites. In the first O site, O is bonded in a distorted single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. In the second O site, O is bonded in a distorted bent 120 degrees geometry to one H and one Cl atom. The O–Cl bond length is 1.47 Å. In the third O site, O is bonded in a trigonal non-coplanar geometry to three H atoms. In the fourth O site, O is bonded in a trigonal non-coplanar geometry to three H atoms. In the fifth O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.49 Å. In the sixth O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.45 Å. Cl is bonded in a tetrahedral geometry to four O atoms. In one of the HClO4 clusters, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 1.00 Å. There are four inequivalent O sites. In the first O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.43 Å. In the second O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.44 Å. In the third O site, O is bonded in a distorted water-like geometry to one H and one Cl atom. The O–Cl bond length is 1.63 Å. In the fourth O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.43 Å. Cl is bonded in a tetrahedral geometry to four O atoms. In one of the HClO4 clusters, H is bonded in a single-bond geometry to one O atom. The H–O bond length is 1.00 Å. There are four inequivalent O sites. In the first O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.43 Å. In the second O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.43 Å. In the third O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.43 Å. In the fourth O site, O is bonded in a distorted water-like geometry to one H and one Cl atom. The O–Cl bond length is 1.64 Å. Cl is bonded in a tetrahedral geometry to four O atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗
DOE OSTI2024

CO 2 hydrogenation over rhodium cluster catalyst nucleated within a manganese oxide framework

Rhodium-based manganese oxide frameworks were explored as a prototype for carbon dioxide reactive capture and conversion. Three-dimensional frameworks of MnOx were utilized as support structures to isolate Rh metal centers. V, Na, and Zn were introduced as counterions to stabilize the structure and for their beneficial effect as promoters. Here, with this multicomponent catalyst, Rh active centers with MnOxs and varied counterions, we were able to selectively tune the catalytic performance of the material via the choice of counterion and structure of the host material. With cryptomelane-type tunnel manganese oxides octahedral molecular sieve (OMS2), we found that Rh-V-OMS2 was highly stable even after 48 hours on stream with a reaction rate of around 1.5x10 -4 mol CO 2 /g Rh /s, surpassing the net reactivity of other initially more active combinations. Furthermore, during CO 2 hydrogenation, in situ XAFS showed that single Rh atoms nucleated into nanoparticles/ sub-nanometer clusters with a coordination number of 5.5 or less. Our finding of the correlation between the reaction rate and particle size offers the potential for enhanced control over the reaction rate by tuning particle size. Our activity study with control experiments demonstrates that the activities of the catalysts are proved due to the unique metal support interaction offered by the Rh-X-MnO.

10 SYNTHETIC FUELS↗
DOE OSTI2024

Spin Decoherence in VOPc@graphene Nanoribbon Complexes

Carbon nanoribbons or nanographene qubit arrays can facilitate quantum-to-quantum transduction between light, charge, and spin, making them an excellent testbed for fundamental science in quantum coherent systems and for the construction of higher-level qubit circuits. Here, in this work, we study spin decoherence due to coupling with a surrounding nuclear spin bath of an electronic molecular spin of a vanadyl phthalocyanine (VOPc) molecule integrated onto an armchair-edged graphene nanoribbon (GNR). Density functional theory (DFT) is used to obtain ground-state atomic configurations. Decay of spin coherence in Hahn echo experiments is then simulated using the cluster correlation expansion method with a spin Hamiltonian involving hyperfine and electric field gradient tensors calculated from DFT. We find that the decoherence time T 2 is anisotropic with respect to magnetic field orientation and determined only by the hydrogen nuclear spins on both VOPc and GNR. Large electron spin echo envelope modulation (ESEEM) due to nitrogen and vanadium nuclear spins is present at specific field ranges and can be completely suppressed by tuning the magnetic field. The relation between these field ranges and the hyperfine interactions is analyzed. The effects of interactions with the nuclear quadrupole moments are also studied, validating the applicability and limitations of the spin Hamiltonian when they are disregarded.

Hamiltonians↗
DOE OSTI2021

Modifying Li + and Anion Diffusivities in Polyacetal Electrolytes: A Pulsed-Field-Gradient NMR Study of Ion Self-Diffusion

Polyacetal electrolytes have been demonstrated as promising alternatives to liquid electrolytes and poly(ethylene oxide) (PEO) for rechargeable lithium-ion batteries; however, the relationship between polymer structure and ion motion is difficult to characterize. Here, we study structure-property trends in ion diffusion with respect to polymer composition for a systematic series of five polyacetals with varying ratios of ethylene oxide (EO) to methylene oxide (MO) units, denoted as P(xEO-yMO), and PEO. We first use 7 Li and 19 F pulsed-field-gradient NMR spectroscopy to measure cation and anion self-diffusion, respectively, in polymer/lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt mixtures. At 90 °C, we observe modest changes in Li + diffusivity across all polymer compositions, while anion (TFSI - ) self-diffusion coefficients decrease significantly with increasing MO content. At a given reduced temperature (T - T g ), all polyacetal electrolytes exhibit faster Li + self-diffusion than PEO. Intriguingly, P(EO-MO) and P(EO-2MO) also show slower TFSI - anion self-diffusion than PEO at a given reduced temperature. Molecular dynamics simulations reveal that shorter distances between acetal oxygen atoms (O-CH 2 -O) compared to ether oxygens (O-CH 2 -CH 2 -O) promote more diverse, often asymmetric, Li + coordination environments. Finally, Raman spectra reveal that anion-rich ion clusters in P(EO-MO) and P(EO-2MO) lead to decreased anion diffusivity, which along with increased cation diffusivity, support the viability of polyacetals as high-performance polymer electrolytes.

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

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

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

Structural, Electronic, and Photophysical Insights into a Few Atom Copper-Sulfur Cluster in the Solid and Solution States

Coinage-metal chalcogenide clusters are widely studied for their attractive photoluminescence properties. Copper chalcogenides are especially promising, but are often confined to solid-state investigations due to their limited solution stability and the difficulty of synthesizing stable, well-defined clusters. Here, we investigate copper–sulfur clusters incorporating a small number of Cu atoms to elucidate fundamental atomic interactions, ground- and excited-state characteristics, and photophysical behavior in both solid and solution. We have synthesized the Cu6(4,6-dimethyl-2-mercaptopyrimidine)6 cluster in both neutral and charged states, Cu6 and Cu6-2H2+, respectively, by selective ligand protonation. The molecular structures are determined using single-crystal X-ray diffraction, while Cu K-edge X-ray absorption spectroscopy is used to probe Cu electronic structure differences arising from the ligand modification. Steady-state and pump-probe optical spectroscopy is used to investigate photophysical properties, interpreted using density functional theory methods. Both clusters exhibit good stability in the solid state and in solution and show characteristic near-infrared emission with microsecond lifetimes. Overall, the Cu6S6 clusters display favorable charge–transfer characteristics and show potential for further use in driving photochemical transformations.

Copper-sulfur clusters↗
DOE OSTI2025

Exact closed-form unitary transformations of fermionic operators

Unitary transformations play a fundamental role in many-body physics, and except for special cases, they are not expressible in closed form. We present closed-form expressions for unitary transformations generated by a single fermionic operator for Hermitian and anti-Hermitian generators. We demonstrate the usefulness of these expressions in formal analyses of unitary transformations and numerical applications to Hamiltonian downfolding in quantum computing and Heisenberg dynamics. Furthermore, this work paves the way for new analytical treatments of unitary transformations and numerical many-body methods for fermions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Unsupervised learning of representative local atomic arrangements in molecular dynamics data

Molecular dynamics (MD) simulations present a data-mining challenge, given that they can generate a considerable amount of data but often rely on limited or biased human interpretation to examine their information content. By not asking the right questions of MD data we may miss critical information hidden within it. Here we combine dimensionality reduction (UMAP) and unsupervised hierarchical clustering (HDBSCAN) to quantitatively characterize prevalent coordination environments of chemical species within MD data. By focusing on local coordination, we significantly reduce the amount of data to be analyzed by extracting all distinct molecular formulas within a given coordination sphere. We then efficiently combine UMAP and HDBSCAN with alignment or shape-matching algorithms to partition these formulas into structural isomer families indicating their relative populations. The method was employed to reveal details of cation coordination in electrolytes based on molecular liquids.

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

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

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

Unveiling the Role of Surface Ir-Oxo Species in O 2 Evolution at IrO 2 Electrocatalysts via Embedded Cluster Multireference Calculations

Understanding the mechanisms driving the oxygen evolution reaction (OER) on iridium oxide (IrO 2 )-based catalysts is essential to improving their performance and enabling an actual scale-up of water-splitting photoelectrochemical cells. The mechanistic pathways at IrO 2 interfaces have been extensively investigated computationally using density functional theory (DFT), which predicts a high-energy barrier for the last step of the OER of molecular oxygen detachment and release from the catalyst surface. Nevertheless, surface O 2 over- and under-binding results by standard generalized gradient approximation and hybrid density functionals, respectively, call for further analysis of this crucial step via multireference methods. Aiming at unveiling the nature of such a barrier, we hereby address the formation of O 2 from the most-stable IrO 2 (110) surface with both periodic DFT and an electrostatic embedded cluster approach at the n-electron valence-state perturbation theory. With this multireference approach, we find a value for the aforementioned energy barrier that is much closer to experimental indications than DFT ones. Here, an in-depth analysis of the involved molecular orbitals suggests that the origin of this barrier is related to the breaking of a π interaction between O 2 and Ir surface atom and to a significant additional O 2 interaction with adjacent electrophilic Ir-oxo species, which is present under experimental operating conditions. Besides shedding light on the mechanism of the OER on IrO 2 , these findings point out the importance of multireference methods for dissecting complex reactions at electrocatalytic interfaces and pave the route for further investigations with effective embedding approaches.

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

Structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquids close to Earth's liquid core conditions

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid close to the Earth's liquid core conditions. Here, the simulation results reveal that the short-range structural order (SRO) in the $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid is very strong. About 80% of the atoms are arranged in crystallinelike SRO motifs. In particular, ~ $70$% of Fe-centered clusters can be classified as either hexagonal-close-packed-like or icosahedrallike SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O-diffusion coefficient is about two to three times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗