Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRUNCATION”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dominant balance-based adaptive mesh refinement for incompressible fluid flows

This work introduces a novel adaptive mesh refinement (AMR) method that utilizes dominant balance analysis (DBA) for efficient and accurate grid adaptation in computational fluid dynamics (CFD) simulations. The proposed method leverages a Gaussian mixture model (GMM) to classify grid cells into active and passive regions based on the dominant physical interactions within the equation space. By modeling truncation error probabilistically from discretized terms, the method identifies regions of high interaction where numerical accuracy is most sensitive to resolution. Unlike traditional AMR strategies, this approach does not rely on heuristic-based sensors or user-defined thresholds, providing a fully automated and problem-independent framework for AMR. Applied to the incompressible Navier-Stokes equations for steady and unsteady flow past a cylinder, the DBA-based AMR method achieves comparable accuracy to high-resolution grids while reducing computational costs by up to 70 %. The validation highlights the method’s effectiveness in capturing complex flow features while minimizing grid cells, directing computational resources toward regions with the most critical dynamics. This modular and scalable strategy is adaptable to a wide range of applications, presenting a promising tool for efficient high-fidelity simulations in CFD and other multiphysics domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗
DOE OSTI2026

Gas-driven short disconnection mitigates thermal runaway in Li-ion batteries under mechanical abuse

Mechanical abuse poses a critical safety risk to Li-ion batteries by inducing internal short circuits that initiate thermal runaway. Remarkably, voltage recovery frequently emerges during thermal runaway initiation, a phenomenon that conventional theories fail to explain. Here, our study develops a new multiphysics mechanism to explain voltage recovery, termed gas-driven short disconnection, whereby internal gas pressure mechanically disengages short-circuit contacts and causes the cell voltage to rebound. This mechanism incorporates gas generation and its structural impact on the short circuit. Real-time optical and thermal imaging and X-ray computed tomography reveal fluid-structure interaction between internal gas flow and adjacent shorting contacts. We establish a mechanistic framework linking gas-driven short disconnection to cell-level voltage and temperature responses, elucidating the extension–truncation pattern of voltage recovery. Furthermore, thermal regime maps show that a voltage recovery duration exceeding 5 s correlates with limited temperature rise below 150 °C, indicating that sustained short-circuit disconnection suppresses Joule heating. Additionally, a dimensionless criterion is deduced from scaling analysis for physical plausibility of gas-driven short disconnection in mechanically abused cells. This finding inspires smart venting control, which regulates gas release to maintain the internal pressure while dissipating gas enthalpy, thereby providing a device-level strategy for thermal runaway mitigation.

Gas-driven short disconnection↗
DOE OSTI2026

Measurements of excitation functions and photoneutron cross sections of 96,98,99,104 Ru and 100 Mo

Photoneutron cross sections were extracted for isotopes of ruthenium and molybdenum on natural abundance targets via a reduced chi-squared analysis on excitation functions, or activation yields, as function of electron beam energy. Bremsstrahlung photons were produced via an electron linear accelerator and a tungsten radiator at the Idaho Accelerator Center. A total of 16 irradiations were performed with electron beam energies ranging from approximately 8 MeV to 23 MeV. Induced radioactivity was measured using a high purity germanium detector. The cross sections were extracted assuming a three parameter Lorentzian fit with a smooth truncated rise at reaction threshold.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2024

How do mirror charge radii constrain density dependence of the symmetry energy?

It has recently been suggested that differences in the charge radii of mirror nuclei (Δ$R^{mirr}_{ch}$) are strongly correlated with the neutron-skin thickness ($R_{skin}$) of neutron-rich nuclei and with the slope of the symmetry energy ($L$). To test this assumption, we present ab initio calculations of $R_{skin}$ in 48 Ca and 208 Pb, Δ$R^{mirr}_{ch}$ in 36 Ca$-$ 36 S, 38 Ca$-$ 38 Ar, 41 Sc$-$ 41 Ca, 48 Ni$-$ 48 Ca, 52 Ni$-$ 52 Cr, and 54 Ni$-$ 54 Fe mirror pairs, and $L$. Employing the recently developed 34 chiral interaction samples, identified by the history matching approach, we conduct rigorous statistical analysis of correlations among Δ$R^{mirr}_{ch}$, $R_{skin}$ and $L$, accounting for quantified uncertainties from low-energy constants of chiral interaction, chiral effective field theory truncation and many-body method approximation. The ab initio results reveal an appreciable Δ$R^{mirr}_{ch}$ $-$ $L$ correlation in $fp$-shell mirror pairs. However, contrary to previous studies, the present calculation finds that the studied $sd$-shell mirror pairs do not exhibit any Δ$R^{mirr}_{ch}$ $-$ $L$ correlation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electrolyte-Induced Restructuring of Acid-Stable Oxygen Evolution Catalysts

Crystalline metal oxide catalysts operating under oxygen evolution reaction (OER) conditions invariably restructure, resulting in active sites with hydroxo/oxo species in an amorphous environment. An increase in the population of terminal hydroxo/oxo species (i.e., edge sites) facilitates proton-coupled electron-transfer (PCET) kinetics for oxygen generation and thus improves catalyst competency. While amorphous films benefit from a greater density of active sites, they suffer from diminished charge transport as compared to that of extended crystalline lattices. Managing this amorphous–crystalline dichotomy is essential when designing OER catalysts, which we highlight with the examination of electrodeposited PbO x materials, which historically are very poor OER catalysts. Along these lines, the presence of phosphate during PbO x electrodeposition truncates the growth of an extended lattice owing to its strong bonding to oxide surfaces to afford an amorphous catalyst film (A-PbO x ) with significant charge-transfer resistance (138 ± 42 Ω) and poor OER kinetics (420 ± 105 mV dec –1 Tafel slope). Conversely, electrodeposition of Pb 2+ in the presence of less coordinating electrolytes such as nitrate affords crystalline β-PbO2 with improved charge-transfer resistance (42.6 ± 1.1 Ω), though still poor OER kinetics (134 ± 36 mV dec –1 Tafel slope). By operating amorphous A-PbOx in less coordinating electrolytes, however, a new partially crystalline material can be generated (μc-PbO x ) with further reduced charge-transfer resistance (33.0 ± 1.4 Ω) and improved OER kinetics (70 ± 15 mV dec –1 Tafel slope). The enhanced OER activity of μc-PbO x is the result of coupling the high edge-site population of an amorphous PbOx phase with crystalline-like charge transport properties. Finally, the ability to use an electrolyte to induce OER activity in an inactive amorphous form of PbO x highlights the benefits of optimizing the amorphous–crystalline phase compositions in the design of active OER catalysts.

catalysts↗
DOE OSTI2024

WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation

We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Finally, other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.

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

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗
DOE OSTI2025

Qubit-Efficient Quantum Chemistry with the ADAPT Variational Quantum Eigensolver and Double Unitary Downfolding

Here, in this work, we combine the recently developed double unitary coupled cluster (DUCC) theory with the adaptive, problem-tailored variational quantum eigensolver (ADAPT-VQE) to explore the accuracy of unitary downfolded Hamiltonians for quantum simulation of chemistry. We benchmark the ability of DUCC effective Hamiltonians to recover dynamical correlation energy outside of an active space. We consider the effects of strong correlation, commutator truncation, higher-body terms, and approximate external amplitudes on the accuracy of these effective Hamiltonians. When combining these DUCC Hamiltonians with ADAPT-VQE, we observe similar convergence of the ground state as compared with bare active space Hamiltonians, demonstrating that DUCC Hamiltonians provide increased accuracy without increasing the load on the quantum processor.

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

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗
DOE OSTI2026

Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.

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

Revisiting the Optical Spectrum of the Plutonyl Ion (PuO 2 ) 2+ in 1 M HClO 4

The analysis of the solution absorption spectrum of the plutonyl ion in an aqueous environment was given by Eisenstein and Pryce (E&P) in 1968. In 2011 a new spectrum was published of the (PuO 2 ) 2+ ion in 1 M HClO 4 . We have been provided with the original data of this spectrum and have found in the data a previously unreported low-lying transition at 7385 cm –1 which we have assigned as a magnetic dipole transition. We have fit most of the near-infrared and optical transitions with Gaussian fits and tabulated a new energy level list up to 22,000 cm –1 which mostly agrees with the data of E&P. We assumed a crystal field of D ∞h (only axial symmetry) and utilized the intensity calculations published for the isoelectronic (NpO 2 ) 1+ ion using a complete basis set for the 5f 2 problem including the Coulombic, spin–orbit as well as the crystal field Hamiltonian. Our results differ substantially from those of E&P. Subsequently, we used a truncated Hamiltonian to try to establish the effects of assuming the σ antibonding orbitals are at such high energies that we can ignore their contributions to the lower lying φ and δ orbitals.

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

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗
DOE OSTI2024

Anatase TiO 2 Surface Hydration: Infrared Vibrational Assignment and Facet-Specific Stability

Anatase TiO 2 nanoparticles support a subtle diversity of surface aqua and hydroxy coordination environments that are characteristic of distinct facets and defect sites. The disparate properties and stability of each molecularly or dissociatively absorbed water environment may provide unique reactivity for facet- and site-selective atomic layer deposition (ALD). Here we probe the temperature-dependent and facet-specific stretching frequency and stability of surface hydroxyls on faceted anatase TiO 2 nanoparticles via in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). Well-defined truncated octahedra, disks, and rods fabricated through hydrothermal synthesis present dominant (101), (001), and (010) facets, respectively. The selective dehydration of discrete local coordination environments as a function of temperature provides a potential pathway to site-selective and facet-selective ALD on metal oxide nanomaterials.

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

When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu 14 and Au 20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps

In view of the tremendous progress in atomically precise metal nanoclusters where electrochemical and optical energetics are routinely supported by computations to establish structure−function correlations, we explore the relationship between different protocols for measuring and computing band gaps of two distinct organic ligand-protected nanoclusters: [Cu 14 H 10 (MBN) 3 (PPh 3 ) 8 ] + and Au 20 (TBBT) 16 . Through UV/visible spectroscopy and differential pulse voltammetry, we measure optical and electrochemical band gaps in those systems. We then compare these experimentally determined gaps to HOMO−LUMO gaps, fundamental gaps, vertical excitation energies, and E o ox − E o red potentials computed using different density functional theory (DFT) or time-dependent DFT (TDDFT) methods. Specifically, in both copper and gold nanoclusters, we test the effect of truncating inert ligands from the model and compare density functionals with varying degrees of Hartree−Fock (HF) exchange from 0 to 50%, range-separated hybrids with a varying long-range tuning parameter, different correlation functionals, basis sets, and (equilibrium and nonequilibrium) continuum solvation models. Despite having different frontier orbital characters (the copper nanocluster has a metal-to-ligand charge transfer character while the gold nanocluster has metal-centered frontier orbitals), both nanoclusters display a similar sensitivity of the HOMO−LUMO gap to the HF exchange that is partially mitigated when computing the fundamental, optical, and electrochemical gaps. Other factors, such as the nature of the correlation functional, basis set, and geometry relaxation, have a considerably smaller effect on computed band gaps in these systems. Overall, this work provides guidelines for factors of varied importance for correlating computed and experimental band gap values.

Chemical calculations↗
DOE OSTI2024

Coherent Phonon Dynamics in Plasmonic Gold Tetrahedral Nanoparticle Ensembles

Coherent phonon modes supported by plasmonic nanoparticles offer prospective applications in chemical and biological sensing. Whereas characterization of these phonon modes often requires single particle measurements, synthetic routes to narrow size distributions of nanoparticles permit ensemble investigations. Recently, the synthesis of highly monodisperse gold tetrahedral nanoparticles with tunable edge length and corner sharpness has been developed. Herein, we characterize a size series of these nanoparticles in colloidal dispersion via transient absorption spectroscopy to examine their mechanical and plasmonic responses upon photoexcitation. Oscillations of transient absorption signals are observed in the plasmon resonance and correspond to the lowest order radial breathing modes of the nanoparticles, the frequencies of which are affected by edge length and truncation of the corners. Furthermore, quality factor values ranging from 10 to 15 are observed for the oscillations that convey potential utility in mass-sensing and plasmon-exciton coupling photonics schemes. Finite-difference time domain and finite element analysis calculations establish the specific optically-relevant phonon modes.

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

Pyrite (001) Interface Chemistry is Controlled by a Sulfoxy Termination

Pyrite (FeS 2 ) is the most common sulfide mineral on Earth, forming through inorganic reactions in the crust and oceanic hydrothermal systems and via microbially driven processes in anaerobic sediments. The pyrite–water interface is the site of a wide range of adsorption and reaction processes in Earth systems including oxidation that dramatically affects the geochemistry of surface waters and influences global carbon and oxygen cycles. Mechanistic geochemical models of pyrite interfacial reactivity, however, are limited by the lack of experimentally derived atomistic structures of the reduced and reacting surfaces. Here, in this work, we reveal the atomic-scale structure of the pyrite (001)-water interface that forms at very low oxygen partial pressures, relevant to suboxic environments in Earth. The interface structure and surface speciation were obtained using the crystal truncation rod method supported by ambient-pressure photoelectron spectroscopy and density functional theory calculations. The surface is dominantly composed of disulfide groups bound to a single oxygen atom, forming a sulfoxy group that has no known molecular or bulk mineral analog. This surface is interpreted as the first step in the oxidative dissolution of pyrite. The sulfoxy group is readily protonated through surface acid–base reactions that alter the structure of interfacial water and the free energy of interfacial reactions. Surface iron sites are not oxidized. Surprisingly, this interface can likely develop in equilibrium with bulk pyrite in some reducing and acidic solutions. This termination is therefore likely representative of pyrite surfaces under a vast range of experimental, industrial and Earth conditions.

oxidation↗