Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Agreements”

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

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium–Ligand Bonds

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy D e and bond length r e are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15–0.19. For the single-reference methods, Kohn–Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the D e and r e of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core–valence terms; this results in highly accurate dissociation energies and good bond lengths. In conclusion, averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

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

On the Rate Constant for NH 2 +HO 2 and Third-Body Collision Efficiencies for NH 2 +H(+M) and NH 2 +NH 2 (+M)

In low-temperature flash photolysis of NH 3 /O 2 /N 2 mixtures, the NH 2 consumption rate and the product distribution is controlled by the reactions NH 2 + HO 2 → products (R1), NH 2 + H (+M) → NH 3 (+M) (R2), and NH 2 + NH 2 (+M) → N 2 H 4 (+M) (R3). In the present work, published flash photolysis experiments by, among others, Cheskis and co-workers, are re-interpreted using recent direct measurements of NH 2 + H (+N 2 ) and NH 2 + NH 2 (+N 2 ) from Altinay and Macdonald. To facilitate analysis of the FP data, relative third-body collision efficiencies compared to N 2 for R2 and R3 were calculated for O 2 and NH 3 as well as for other selected molecules. We report results were in good agreement with the limited experimental data. Based on reported NH 2 decay rates in flash photolysis of NH 3 /O 2 /N 2 , a rate constant for NH 2 + HO 2 → NH 3 + O 2 (R1a) of $k_{1\text{a}}$ = 1.5(±0.5) × 10 14 cm 3 mol –1 s –1 at 295 K was derived. This value is higher than earlier determinations based on the FP results but in good agreement with recent theoretical work. Kinetic modeling of reported N 2 O yields indicates that NH 2 + HO 2 → H 2 NO + O (R1c) is competing with R1a, but perturbation experiments with addition of CH4 indicate that it is not a dominating channel. Measured HNO profiles indicate that this component is formed directly by NH 2 + HO 2 → HNO + H 2 O (R1b), but theoretical work indicates that R1b is only a minor channel. Based on this analysis, we estimate $k_{1\text{c}}$ = 2.5 × 10 13 cm 3 mol –1 s –1 and $k_{1\text{b}}$ = 2.5 × 10 12 cm 3 mol –1 s –1 at 295 K, with significant uncertainty margins.

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

Cryogenic Vibrationally Resolved Photoelectron Spectroscopy of OH–(H2O): Confirmation of Multidimensional Franck-Condon Simulation Results for the Transition State of the OH + H2O Reaction

We present a transition state spectroscopic study of the OH + H2O reaction using the experimental technique of cryogenic Negative Ion Photoelectron Spectroscopy (NIPES). The recorded NIPE spectrum at 193 nm exhibits multiple vibrational progressions that include excitations to the shared H atom antisymmetric stretching mode with an interval of 0.32 eV as well as other progressions, mainly involving the H bending and O…O symmetric stretching modes. The Vertical Detachment Energy (VDE) was measured at 3.53 eV, whereas an upper limit for the Adiabatic Detachment Energy (ADE) was estimated at 2.90 eV. These values are in excellent agreement with the theoretically computed values of 3.51 eV and 2.83 eV respectively, obtained at the CCSD(T)/aug-cc-pV5Z level of theory. The recorded NIPE spectrum is in very good agreement when compared to the one recently reported from four-dimensional Franck-Condon simulations, in which a similar spectral profile was predicted. Besides observing the ground state, a charge-transfer excited state in the form of [OH–(H2O)+] is identified with a relative energy of 1.39 eV, well matching the previous prediction of 1.36 eV. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer. This research also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

Cao, Wenjin↗
DOE OSTI2022

Thermodynamics of the Metal Carbonates and Bicarbonates of Mn, Co, Ni, Cu, and Zn Relevant to Mineral Energetics

The gas phase heats of formation of ground-state MCO 3 , M(HCO 3 ) 2 , and M(HCO 3 )(OH), where M = Mn, Co, Ni, Cu, and Zn, have been predicted using the correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit using the Feller–Peterson–Dixon (FPD) approach. Cohesive energies of the carbonates were predicted based on the calculated gas phase and experimental solid heats of formation. Coulombic dissociation energies (CDEs) between metal cations and anions show a near-linear correlation with Shannon metal cation atomic radii, yet no correlation is found with the hardness of these cations. The total reaction dissociation energies (TRDEs) of transition metals are higher than their CDEs for the di-bicarbonates, in contrast to those for Mg and Ca based on our prior work. In addition to differences in the energies needed to prepare the transition metal dications, electron donation from the ligands to the 3d orbitals of open-shell transition metal dications from lone pairs of adjacent O atoms also plays a role. No electron donation from the ligands to the fully occupied 3d orbitals of Zn and Cd was found. Decomposition energies for generating MO, CO 2 , and/or H 2 O were calculated. Gas phase metal exchange energies only partially correlate with the electrochemical series for M(s) → M 2+( aq). The FPD heats of formation were used to benchmark a range of density functional theory exchange–correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. So, the best agreement with the FPD results is predicted for the $\tau$-HCTH functional with an average unsigned error of 8.3 kcal/mol.

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

Theory-Based Mechanism for Fluoromethane Combustion I: Thermochemistry and Abstraction Reactions

A new detailed chemical kinetic mechanism is presented for small fluorinated hydrocarbons. Ab initio electronic structure theory is used to provide heats of formation with subchemical accuracy. The ANL0 method is extended to include fluorine. The resulting heats of formation at 0 K are in excellent agreement with 36 benchmark species in the Active Thermochemical Tables, with a mean error of μ = –0.02 kJ/mol and a standard deviation of σ = 0.91 kJ/mol. The thermophysical properties for 92 small-molecule H/C/O/F species are computed. The rate coefficients for 40+ H-abstraction reactions involving H, O, F, OH, OF, HO 2 , and various methyl radicals with CH 4 , CH 3 F, CH 2 F 2 , CHF 3 , CH 2 O, and CHFO are discussed. Here, the computed rate constants are in excellent agreement with the available literature. Additionally, 30+ rate constants are provided for F abstraction, which are several orders of magnitude smaller than H abstraction. The thermophysical properties and rate constants are provided in a mechanism. This mechanism is the first in a series of theory-based investigations into the thermal destruction of per- and polyfluorinated species.

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

Importance of Dispersion in the Molecular Geometries of Mn(III) Spin-Crossover Complexes

The computational investigation of the molecular geometries of a pair of manganese(III) spin-crossover complexes is reported. For the geometry of the quintet high-spin state, density functionals significantly overestimate Mn–Namine bond distances, although the geometry for the triplet intermediate-spin state is well described. Here, comparisons with several wave function-based methods demonstrate that this error is due to the limited ability of commonly used density functionals to recover dispersion beyond a certain extent. Among the methods employed for geometry optimization, restricted open-shell Møller–Plesset perturbation theory (MP2) appropriately describes the high-spin geometry but results in a slightly shorter Mn–O distance in both spin states. On the other hand, extended multistate complete active space second-order perturbation theory (XMS-CASPT2) provides a good description of the geometry for the intermediate-spin state but also sufficiently recovers dispersion, performing well for the high-spin state. Despite the fact that the electronic structure of both spin states is dominated by one-electron configuration, XMS-CASPT2 offers a balanced approach, leading to molecular geometries with much better agreement with experiment than MP2 and DFT. A scan along the Mn–N amine bond demonstrates that for these complexes coupled cluster methods (i.e., DLPNO-CCSD(T)) also yield bond distances in agreement with experiment while multiconfiguration pair density functional theory (MC-PDFT) is unable to recover dispersion well enough, analogous to single-reference DFT.

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

G 0 W 0 Ionization Potentials of First-Row Transition Metal Aqua Ions

Here, we report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V–Cu) aqua ions in their most common oxidation states, i.e., V 3+ , Cr 3+ , Cr 2+ , Mn 2+ , Fe 3+ , Fe 2+ , Co 2+ , Ni 2+ , and Cu 2+ . The d-orbital occupancy of these systems spans a broad range from d 2 to d 9 . All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G 0 W 0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G 0 W 0 @PBE, G 0 W 0 @PBE0, and G 0 W 0 @r 2 SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G 0 W 0 @PBE0 and G 0 W 0 @r 2 SCAN approaches, with the r 2 SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.

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

R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

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

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

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

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗
DOE OSTI2020

Prediction of Structures and Atomization Energies of Coinage Metals, (M)n, n < 20: Extrapolation of Normalized Clustering Energies to Predict the Cohesive Energy

The geometries of the group 11 coinage metals (n = 2-20) were optimized to determine the lowest energy isomers for each cluster size, singlets for even numbers and doublets for odd numbers. For copper and silver, 2-D (planar) geometries were favored up to n = 6. For gold, 2D (planar) geometries were favored up to n = 13. Normalized clustering energies were plotted as a function of cluster size (n-1/3, for n = 4-20) with various DFT functionals and the CCSD(T)-F12b method and were extrapolated to predict the bulk cohesive energy. In the case of copper and silver, there is excellent agreement between the cohesive energies predicted at the CCSD(T)-F12b level of theory and the experimental values. For gold, the CCSD(T)-F12b values needed to be corrected for spin-orbit relativistic effects to obtain good agreement with experiment. Electronic properties including the HOMO-LUMO gaps for the even clusters and the spin densities for the odd clusters were calculated. The lowest gap is predicted to occur for n = 16 where the HOMO and LUMO are very similar in shape.

Persuad, Rudradatt R.↗
DOE OSTI2021

First-Principles Simulations of CuCl in High-Temperature Water Vapor

Experimental data suggest that the solubility of copper in high-temperature water vapor is controlled by the formation of hydrated clusters of the form CuCl(H 2 O) n , where the average number of water molecules in the cluster generally increases with increasing density [Migdisov, A. A.; et al. Geochim. Cosmochim. Acta2014, 129, 33–53]. However, the precise nature of these clusters is difficult to probe experimentally. Moreover, there are some discrepancies between experimental estimates of average cluster size and prior simulation work [Mei, Y. Geofluids2018, 2018, 4279124]. Here, we have performed first-principles Monte Carlo (MC) and molecular dynamics (MD) simulations to explore these clusters in finer detail. We find that molecular dynamics is not the most appropriate technique for studying aggregation in vapor phases, even at relatively high temperatures. Specifically, our MD simulations exhibit substantial problems in adequately sampling the equilibrium cluster size distribution. In contrast, MC simulations with specialized cluster moves are able to accurately sample the phase space of hydrogen-bonding vapors. At all densities, we find a stable, slightly distorted linear H 2 O–Cu–Cl structure, which is in agreement with the earlier simulations, surrounded by a variable number of water molecules. The surrounding water molecules do not form a well-defined second solvation shell but rather a loose network of hydrogen-bonded water with molecular CuCl on the outside edge of the water cluster. We also find a broad distribution of hydration numbers, especially at higher densities. In contrast to previous simulation work but in agreement with experimental data, we find that the average hydration number substantially increases with increasing density. Moreover, the value of the hydration number depends on the choice of cluster definition.

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

Equipartition and the Temperature of Maximum Density of TIP4P/2005 Water

Here, we simulate TIP4P/2005 water in the temperature range of 257 to 318 K with time-steps δ = 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density–temperature behavior obtained using 0.25 or 0.50 fs is in excellent agreement with each other but differs from those obtained using time steps that have been shown earlier to lead to a breakdown of equipartition. For δt = 0.25 or 0.50 fs, the temperature of maximum density (TMD) is 277.15 K and the density value is in close agreement with experiments. For δt = 1.00 fs, the TMD is 277.15 K, but the density value is shifted higher. For the other time steps considered here, the TMD is shifted to progressively lower values for longer time steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water–water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid–vapor phase envelope. We present a simple physically transparent explanation that highlights the separation of the time scales between translational and rotational motion. We also develop a metric, χ, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating χ is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time steps (≤0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct ensemble sampling is fundamental to ensure reproducibility across codes and simulation algorithms.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
DOE OSTI2026

Reparameterizing a Lipid Force Field Using Small-Angle X-ray Scattering to Improve Predictions of Multicomponent Membranes under Organic Solvent Stress

Understanding and predicting properties of lipid-bilayer membranes are essential to elucidating organismal physiology and pathophysiology. Therefore, substantial efforts have been undertaken to derive accurate molecular-mechanics force fields (FF) to allow simulation of their properties; however, much of these past efforts focused on tuning force fields to accurately reproduce the properties of model single-component membranes and not multicomponent or stressed membranes. Here, we tuned the CHARMM36 FF for a simple 2-component model of a Gram-positive bacterium. This updated force field is parametrized against a single condensed-phase property, the small-angle X-ray scattering (SAXS) intensities of the 2-component lipid system, using a version of ForceBalance implemented previously (named ForceBalance-SAS) with SAXS and small-angle neutron scattering (SANS) intensities as optimization targets. After tuning, we observe improved agreement with the experimental SAXS of the membrane under 1-butanol and tetrahydrofuran solvent stresses, with a factor of 10.0 and 7.5 reduction, respectively, in χ 2 , the measure of the discrepancy between experimental and computed SAXS intensities. Furthermore, this reparametrized force field yielded improved agreement between experimental and simulated membrane thicknesses for the pure lipid systems. However, we note limitations in transferability and diagnose the source of such limitations, particularly the need for SANS in addition to SAXS.

Demerdash, Omar N. A. [Oak Ridge National Laborato↗
DOE OSTI2020

Interpreting the Operando XANES of Surface-Supported Subnanometer Clusters: When Fluxionality, Oxidation State, and Size Effect Fight.

X-ray absorption near edge structure (XANES) spectroscopy is widely used for operando catalyst characterization. We show that, for highly fluxional supported nanoclusters, the customary extraction of the oxidation state of the metal from the XANES data by fitting to the bulk standards is highly questionable. The XANES signatures as well as the apparent oxidation state for such clusters arise from a complex combination of many factors, and not only from the chemical composition in reaction conditions (e.g., oxygen content in oxidizing atmosphere). The thermally accessible isomerization and population of several structurally distinct cluster forms, cluster-support interaction, and intrinsic size effects all impact the metal oxidation state and XANES signal. We demonstrate this on copper oxide clusters with different compositions, Cu4Ox(x = 2-5) and Cu5Oy (y = 3, 5), deposited on amorphous alumina and ultrananocrystalline diamond, for which we computed the XANES spectra and compare the results to the experiment. We show in addition that fitting the experimental spectrum to calculated spectra of supported clusters can, in contrast, provide good agreement and insight into the spectrum-composition-structure relation. Experimental XANES interpreted using the proposed fitting scheme shows the partial reduction of Cu oxide clusters at rising temperatures, and pinpoints the specific stoichiometries that dominate in the ensemble of cluster states as the temperature changes.

Zandkarimi, Borna↗
DOE OSTI2020

Closing the Gap Between Experiment and Theory: Reactive Scattering of HCl from Au(111)

Accurate simulation of molecules reacting on metal surfaces, which can help in improving heterogeneous catalysts, remains out of reach for several reactions. For example, a large disagreement between theory and experiment for HCl reacting on Au(111) still remains, despite many efforts. In this work, the dissociative chemisorption of HCl on Au(111) is investigated with a recently developed MGGA density functional (MS-RPBEl) and a high-dimensional neural network potential. Additionally, previous experimental sticking probabilities are re-examined. A considerably improved agreement between experiment and theory is obtained, although theory still overestimates experimental sticking probabilities by a factor of 2-7 at the highest incidence energy. Computed and measured vibrational transition probabilities are also in improved agreement. Several dynamical effects such as angular steering and energy transfer from the molecule to the surface are found to play an important role.

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

A Comparative Study of Electrical Double Layer Effects for CO Reduction Reaction Kinetics

Solvation models describe how the interactions between the solutes and solvents affect the reactivity and selectivity in electrochemical processes. In this study, we developed a framework for evaluating the effects of applied potential and electrical double layer on CO reduction (COR), comparing fully explicit, implicit, and hybrid solvation models at the standard hydrogen electrode (SHE) scale. We analyzed all crucial intermediates leading to the production of C 1 and C 2 products and found good agreement across these models. Some notable differences were observed in the implicit description of *C and *CO adsorption at higher overpotentials and overall trends in the adsorption energies within the hybrid model. Using this unified SHE framework, we built comparable microkinetic models for COR kinetics and rates. Despite small differences in thermodynamics, solvation-model-based microkinetic simulations showed good agreement for onset potentials against the benchmark experiment. Only qualitative difference was observed for C 1+ versus hydrogen evolution at high overpotentials for the implicit model. Finally, we constructed a generalized C 2 selectivity map in descriptor space (U SHE , ΔG electrolyte CO ), which highlights the limitations of a copper-based COR catalyst and guides the search for optimal descriptor parameters to maximize C 2 selectivity. In conclusion, these findings demonstrate the importance of considering electrical double layer effects in reduction reactions and offer a useful framework for comparing solvation models and predicting optimal electrochemical conditions for specific applications.

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

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗