Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive molecular dynamics”

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

Analysis of an all-solid state nanobattery using molecular dynamics simulations under an external electric field

Present Li-ion battery (LIB) technology requires strong improvements in performance, energy capacity, charging-time, and cost to expand their application to e-mobility and grid storage. Li-metal is one of the most promising materials to replace commercial anodes such as graphite because of its 10 times higher specific capacity. However, Li-metal has high reactivity with commercial liquid electrolytes; thus, new solid materials are proposed to replace liquid electrolytes when Li-metal anodes are used. We present a theoretical analysis of the charging process in a full nanobattery, containing a LiCoO 2 cathode, a Li 7 P2S 8 I solid-state electrolyte (SSE), a Li-metal anode as well as Al and Cu collectors for the cathode and anode, respectively. In addition, we added a Li 3 P/Li 2 S film as a solid electrolyte interphase (SEI) layer between the Li-anode and SSE. Thus, we focus this study on the SEI and SSE. We simulated the charging of the nanobattery with an external voltage by applying an electric field. We estimated temperature profiles within the nanobattery and analyzed Li-ion transport through the SSE and SEI. Here, we observed a slight temperature rise at the SEI due to reactions forming $PS_{3}^{–}$ and $P_{2}S_{7}^{4}$$^{–}$ fragments at the interfaces; however, this temperature profile changes due to the charging current under the presence of the external electric field ε = 0.75 V Å –1 . Without the external field, the calculated open-circuit voltage (OCV) was 3.86 V for the battery, which is within the range of values of commercial cobalt-based LIBs. This voltage implies a spontaneous fall of available Li-ions from the anode to the cathode (during discharge). The charge of this nanobattery requires overcoming the OCV plus an additional voltage that determines the charging current. Thus, we applied an external potential able to neutralize the OCV, plus an additional 1.6 V to induce the transport of Li + from the cathode up to the anode. Several interesting details about Li + transport paths through the SSE and SEI are discussed.

25 ENERGY STORAGE↗
DOE OSTI2020

Mechanistic understanding of catalytic conversion of ethanol to 1-butene over 2D-pillared MFI zeolite

Ethanol is an important C2 platform molecule for producing value-added chemicals and distillate hydrocarbon fuels (e.g., jet and diesel). Among these, catalytic upgrading of ethanol to butenes can generate valuable commodity chemicals (e.g., 1-butene) and provide C4 olefin intermediates that can be further upgraded to jet/diesel fuels. 2D zeolites offer hierarchical mesoporous structures, leading to improved mass transport and reduced diffusion length, which can help to address the coking challenges faced by ethanol conversion to hydrocarbons over 3D zeolites. In this study, we investigate the acid-catalyzed conversion of ethanol to 1-butene over the Brønsted acid sites (BAS) in 2D-pillared MFI zeolite (2D-PMFI) using ab initio molecular dynamics (AIMD) simulations, in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and calorimetric measurements. A detailed thermodynamic analysis, using quasi-harmonic approximation (QHA), on the Gibbs free energy pathway of ethanol conversion shows that the consideration of entropy is critical to accurately capture the detailed thermodynamic profiles. The formation of bulky ethoxy-ethene complex is found to be a potential rate-determining step (RDS), proceeding via stepwise mechanism. The reactivity of 2D-PMFI can be further tuned by manipulating RDS through a careful control of number of BAS and operating temperatures. The accuracy of computational model is confirmed through the comparison of calculated vibrational density of states (VDOS) against the experimental DRIFTS measurements. Overall, our study provides mechanistic insights into ethanol upgrading over the 2D-PMFI and shows the importance of evaluating entropic effects in such confined system.

Ethanol conversion, 2D-pillared zeolite, MFI, Anha↗
DOE OSTI2025

De Novo Design of Proteins That Bind Naphthalenediimides, Powerful Photooxidants with Tunable Photophysical Properties

De novo protein design provides a framework to test our understanding of protein function and build proteins with cofactors and functions not found in nature. Here, we report the design of proteins designed to bind powerful photooxidants and the evaluation of the use of these proteins to generate diffusible small-molecule reactive species. Because excited-state dynamics are influenced by the dynamics and hydration of a photooxidant’s environment, it was important to not only design a binding site but also to evaluate its dynamic properties. Thus, we used computational design in conjunction with molecular dynamics (MD) simulations to design a protein, designated NBP (NDI Binding Protein), that held a naphthalenediimide (NDI), a powerful photooxidant, in a programmable molecular environment. Solution NMR confirmed the structure of the complex. We evaluated two NDI cofactors in this de novo protein using ultrafast pump–probe spectroscopy to evaluate light-triggered intra- and intermolecular electron transfer function. Moreover, we demonstrated the utility of this platform to activate multiple molecular probes for protein labeling.

carbonyls↗
DOE OSTI2023

Unexpected Kinetic Solvent Effects Enhance Activity and Selectivity in Biphasic Systems

Biphasic dehydration of fructose to 5-hydroxymethylfurfural (HMF) has shown unprecedented increases in productivity, but a mechanistic understanding is lacking. Herein, we couple fast experimental reaction kinetics, multiscale modeling (phase behavior, classical molecular dynamics(MD), and quantum mechanics/molecular mechanics MD), in situ sampling, and IR and 13 C-NMR spectroscopy to elucidate the complex effects of nonpolar extracting organic solvents on the kinetics of fructose dehydration. We show that these organic solvents can reach significant mutual solubility with water at reaction temperatures, enabling the partition of the sugar and catalyst into the extracting phase. In the organic-rich environment, the dehydration of fructose proceeds faster and more selectively than in water due to increased relative abundance of the reactive furanose isomer, enhanced water–catalyst–substrate interactions driven by nanophase separation, and higher product stability stemming from preferential solvation. Furthermore, we demonstrate that these solvent effects impact other critical biphasic reactions in biomass upgrading and provide qualitative principles for solvent selection.

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

Distinct Gas-Particle Partitioning and Viscosity Characteristics of Secondary Organic Aerosols Derived from α-Pinene versus Ocimene

Secondary organic aerosols (SOA) have complex, multicomponent composition that controls particle viscosity and gas-particle partitioning, key factors to their atmospheric evolution. This study investigates the chemical composition, volatility and viscosity of SOA formed by ozonolysis of cyclic α-pinene (PSOA) and acyclic ocimene (OSOA) monoterpenes. Using Temperature-Programmed Desorption combined with Direct Analysis in Real-Time ionization and High-Resolution Mass Spectrometry, we determined the molecular composition and saturation mass concentration of individual SOA constituents. These data enabled gas-particle partitioning and viscosity estimates under varied atmospheric conditions. PSOA, composed of higher molecular weight and less oxidized species, shows higher condensability and viscosity under high total organic mass (tOM) loadings. Here, in contrast, OSOA, consisting of more oxidized, lower molecular weight species, exhibits greater sensitivity to tOM, with viscosity increasing significantly upon dilution. Poke-flow experiments support this trend, indicating that OSOA undergoes more dynamic compositional and phase changes during atmospheric aging. These observations reveal distinct dynamic trends in the atmospheric transformations and reactivity of SOA from cyclic and acyclic monoterpenes, with the latter showing greater compositional changes during aging that alter viscosity and diffusion. This highlights the importance of incorporating such dynamic transformations into atmospheric models to improve predictions of SOA atmospheric loadings, lifetimes, and impacts.

cyclic and acyclic monoterpenes↗
DOE OSTI2023

Atomic-scale modeling of the dissolution of oxidized platinum nanoparticles in an explicit water environment

Pt nanoparticles (NPs) are currently being investigated for use in fuel cells; however, Pt NP oxidation as a function of size, morphology, and temperature is not well understood or currently quantified. In this study, the stability and dissolution of oxidized platinum NPs is examined via classical molecular dynamics in an explicit water environment. The NPs considered range in size from 1.35 to 2.92 nm in diameter and included five different monolayer (ML) coverages of O*. The simulations were performed at 300, 450, and 600 K with the many-body, reactive third-generation charge-optimized many-body or COMB3 potentials and examine the kinetics of NP dissolution in water. The Pt–O layer, which reduces the kinetic activity for Pt atom dissolution, is projected to make dissolution more favorable for O* MLs smaller than 0.5. The simulations further indicate that the Pt NPs' kinetic rates of dissolution are slowed by an increase in the number of adsorbed species caused by the dissociation of water molecules at the reconstructed facets of the Pt NPs. Furthermore, these findings quantify the effect of oxygen and temperature on the stability and dissolution of oxidized platinum NP in an explicit water environment similar to the conditions in fuel cells and electrocatalysis.

08 HYDROGEN↗
DOE OSTI2022

Identification of efficient and stable solid-solid electrochemical interfaces for Li-S batteries. Final Report

The objective of this work is to develop a better understanding of the interactions between battery components that drive the observed battery performance. Using computational tools at the atomistic level, we elucidate molecular, crystal, and amorphous structures, interfacial structures and reactions, and electrochemical behavior. We focus on Li metal anodes, liquid electrolytes, solid electrolytes, S cathodes, and metal-oxide cathodes, and their corresponding interfaces. Within liquid electrolytes, we study their bulk structure and dynamics, as a function of salt concentration. The analysis allows to understand the formation of localized high concentrated regions as the solvent is substituted by a weakly interacting diluent species. Analysis of the electron distribution yields the reactive behavior at interfaces. Moreover, we employ tools that allow the characterization of ion and electron transport at interfaces, during the processes of solvation and desolvation occurring while cycling. For Sulfur cathodes we elucidate synthesis processes and the role of carbon-sulfur and nitrogen-sulfur interactions not only to provide electronic conductivity but also to define alternative reaction pathways during discharge of Li-S batteries. Full-cell simulations offer new insights regarding the integration of all the battery components. Our studies also include solid electrolytes and their interfaces with Li metal, S-C, and metal-oxide electrodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Designing a Dynamic Material via Interlocking Kraft Lignin with Ultrahigh Molecular Weight Poly(ethylene oxide)

The environmental issues stemming from plastic waste and the excessive use of petroleum-based chemicals are both concerning and urgent. To effectively address this problem, we need to adopt a comprehensive approach involving developing more sustainable and renewable materials. These materials should be capable of demonstrating similar or even better performance than functional polymers synthesized from petroleum-based chemicals. Our research aims to develop more sustainable materials by leveraging the intrinsic characteristics of a natural polymer, lignin, a byproduct of the biorefinery industries. We utilized the three-dimensional branching structure of lignin as a material framework. We employed an approach to co-reactive melt processing of kraft lignin with aliphatic flexible chains of soft cross-linkers through the reaction of the cross-linker epoxy chain ends with lignin functional groups. Our study demonstrates that the interlocked structure formed from the co-blending of kraft lignin with ultrahigh molecular weight poly(ethylene oxide) and the cross-linking of lignin chains through a solvent-free process can manipulate the macromolecular interactions and relaxation. This manipulation results in materials that exhibit a wide range of thermomechanical properties and self-healing and shape memory effects, with drastically improved stiffness in a single material. We have explored the fundamental understanding of the macromolecular chain relaxation dynamics originating from the interlocking structure formation. In conclusion, this investigation utilized various techniques, including thermal, mechanical, rheology, and quasi-elastic neutron scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗
DOE OSTI2026

pH-Driven Restructuring of Hydration Layers and Cation Adsorption at the Alumina–Water Interface

Oxide−water interfaces underpin ion separations, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. However, the molecular-scale links between pHdependent surface protonation, hydration-layer structure, and counterion adsorption remain poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)−water interface in 10 mM Rb + over pH 3−12. CTR measurements reveal two distinct adsorbed water layers at ∼2.2 and ∼3.5 Å above the surface. Each water layer shifts toward the substrate at transition pHs near 6.5 and 10.6, respectively, reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb + coverage and a decrease in mean adsorption height from ∼3.5 to ∼2.7 Å with increasing pH, indicating enhanced counterion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., the potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥ 3.5. The negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb + coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid−base chemistry to ion binding distances at an oxide−water interface.

Adsorption↗
DOE OSTI2023

Applying the core-satellite species concept: Characteristics of rare and common riverine dissolved organic matter

Introduction: Dissolved organic matter (DOM) composition varies over space and time, with a multitude of factors driving the presence or absence of each compound found in the complex DOM mixture. Compounds ubiquitously present across a wide range of river systems (hereafter termed core compounds) may differ in chemical composition and reactivity from compounds present in only a few settings (hereafter termed satellite compounds). Here, we investigated the spatial patterns in DOM molecular formulae presence (occupancy) in surface water and sediments across 97 river corridors at a continental scale using the “Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems—WHONDRS” research consortium. Methods: We used a novel data-driven approach to identify core and satellite compounds and compared their molecular properties identified with Fourier-transform ion cyclotron resonance mass spectrometry (FT-ICR MS). Results: In this work, we found that core compounds clustered around intermediate hydrogen/carbon and oxygen/carbon ratios across both sediment and surface water samples, whereas the satellite compounds varied widely in their elemental composition. Within surface water samples, core compounds were dominated by lignin-like formulae, whereas protein-like formulae dominated the core pool in sediment samples. In contrast, satellite molecular formulae were more evenly distributed between compound classes in both sediment and water molecules. Core compounds found in both sediment and water exhibited lower molecular mass, lower oxidation state, and a higher degree of aromaticity, and were inferred to be more persistent than global satellite compounds. Higher putative biochemical transformations were found in core than satellite compounds, suggesting that the core pool was more processed. Discussion: The observed differences in chemical properties of core and satellite compounds point to potential differences in their sources and contribution to DOM processing in river corridors. Overall, our work points to the potential of data-driven approaches separating rare and common compounds to reduce some of the complexity inherent in studying riverine DOM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Isomer-specific spectroscopy of aromatics and their resonance-stabilized radical intermediates

There are many circumstances in which complex gas-phase reaction mixtures must be interrogated and understood, including combustion, plasmas, atmospheres, and molecular clouds. Spectroscopic tools are essential to selective detection and characterization of the widening array of precursors and the reactive intermediates they generate. This research program has focused on the development and implementation of sensitive spectroscopic detection schemes that are isomer-specific and even conformation-specific. In particular, we used an array of laser-based and broadband microwave methods to carry out isomer-specific and conformation-specific spectroscopy on large flexible molecules and reactive intermediates formed by photolysis, pyrolysis, combustion, or electron-driven chemistry.

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

The Role of Nonequilibrium Solvent Effects in Enhancing Direct CO 2 Capture at the Air–Aqueous Amino Acid Interface

Direct air capture (DAC) technologies are limited by the poor understanding of the dynamic role of interfaces in modulating the chemisorption of CO 2 from air into solutions. While the reactivity of aqueous amine-based solvents in the bulk environment is strongly inhibited by nonequilibrium solvent effects, promoting DAC at interfaces posits a possibility to reduce the coupling with the solvent and significantly accelerate DAC. Building on an experimentally proven concept to bring an anionic glycine absorbent to the interface through ion-pairing interactions with a positively charged surfactant, we establish the fundamental time scales for key elementary steps involved in DAC with rate theory and enhanced-sampling ab initio molecular dynamics simulations. Here, we elucidate the mechanism by which water influences the free energy barriers and dynamical crossing-recrossing of those barriers, affecting the reaction rates. We find that water reorganizes to partially dehydrate [-NH 2 ], facilitating S N 2-based CO 2 conversion to a zwitterion, which then releases a proton via overhydration of [-NH 2 ]. The low-density interfacial water favors dehydration over overhydration, leading to a comparatively higher barrier (slower kinetics) for proton release. The barrier-recrossing events neutralize this effect, letting both steps occur at the same time scale (sub-microseconds) and making the overall DAC process faster at the interface than in the bulk water. Such an understanding of environment-sensitive solvent effects on the reaction kinetics will help design tailored interfaces for enhanced CO 2 capture kinetics via control of solvation and ion paring.

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

Nanoscale wetting controls reactive Pd ensembles in synthesis of dilute PdAu alloy catalysts

The performance of bimetallic dilute alloy catalysts is largely determined by the size of minority metal ensembles on the nanoparticle surface. By analyzing the synthesis of catalysts comprising Pd 8 Au 92 nanoparticles supported on silica using surface-sensitive techniques, we report that whether Pd overgrowth occurs before or after Au nanoparticle deposition onto the support controls the surface Pd ensemble size and abundance. These differences in Pd ensembles influence catalytic reactivity in H 2 –D 2 isotope exchange and benzaldehyde hydrogenation, which, in correlation with theoretical calculations, is used to elucidate the active site(s) in each reaction. To clarify how the synthetic sequence controls the formation of Pd ensembles, we combine numerical wetting calculations and molecular dynamics simulations (with a machine-learned force field) to visualize Pd deposition and migration on the nanoparticle surface, respectively. Our results suggest that the nanoparticle–support interface restricts nanoparticle accessibility to Pd deposition, which consequently controls the Pd ensemble size, illustrating the critical role of nanoscale wetting phenomena during bimetallic catalyst preparation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗
DOE OSTI2021

Persistent Homology Metrics Reveal Quantum Fluctuations and Reactive Atoms in Path Integral Dynamics

Nuclear quantum effects (NQEs) are known to impact a number of features associated with chemical reactivity and physicochemical properties, particularly for light atoms and at low temperatures. In the imaginary time path integral formalism, each atom is mapped onto a “ring polymer” whose spread is related to the quantum mechanical uncertainty in the particle’s position, i.e., its thermal wavelength. A number of metrics have previously been used to investigate and characterize this spread and explain effects arising from quantum delocalization, zero-point energy, and tunneling. Many of these shape metrics consider just the instantaneous structure of the ring polymers. However, given the significant interest in methods such as centroid molecular dynamics and ring polymer molecular dynamics that link the molecular dynamics of these ring polymers to real time properties, there exists significant opportunity to exploit metrics that also allow for the study of the fluctuations of the atom delocalization in time. Here we consider the ring polymer delocalization from the perspective of computational topology, specifically persistent homology, which describes the 3-dimensional arrangement of point cloud data, (i.e. atomic positions). We employ the Betti sequence probability distribution to define the ensemble of shapes adopted by the ring polymer. The Wasserstein distances of Betti sequences adjacent in time are used to characterize fluctuations in shape, where the Fourier transform and associated principal components provides added information differentiating atoms with different NQEs based on their dynamic properties. We demonstrate this methodology on two representative systems, a glassy system consisting of two atom types with dramatically different de Broglie thermal wavelengths, and ab initio molecular dynamics simulation of an aqueous 4 M HCl solution where the H-atoms are differentiated based on their participation in proton transfer reactions.

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

Exploring the Unusual Reactivity of the Hydrated Electron with CO 2

Many questions remain about the reactions of the hydrated electron despite decades of study. Of particular note is the fact that they do not appear to follow the Marcus theory of electron transfer reactions, a feature that has yet to be explained. To investigate these issues, we use ab initio molecular dynamics (AIMD) simulations to investigate a one of the better studied reactions, the hydrated electron reduction of CO 2 . The rate constant for the hydrated electron-CO 2 reaction complex to react to form CO 2 - is, for the first time, estimated from AIMD simulations. Results at 298 and 373 K show the rate constant is insensitive to temperature, consistent with the low measured activation energy for the reaction, and the implications of this behavior are examined. The sampling provided by the simulations yields insight into the reaction mechanism. The reaction is found to involve both solvent reorganization and changes in the carbon dioxide structure. The latter lead to significant vibrational excitation of the bending and symmetric stretch vibrations in the CO 2 - product, indicating the reaction is vibrationally nonadiabatic. The former is estimated from calculation of an approximate collective solvent coordinate and the free energy in this coordinate is determined. Furthermore, these results indicate that AIMD simulations can reasonably estimate hydrated electron reaction activation energies and provide new insight into the mechanism that can help illuminate the features of this unusual chemistry.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Modeling prebiotic chemistries with quantum accuracy at classical costs

Molecular Dynamics (MD) simulations using classical force-fields are commonly employed in numerous scientific investigations. However, many natural processes involve bond breaking and quantum forces. This complexity is compounded by the presence of multiple competing length and timescales. For example, accurately modeling the thermodynamics and dynamics of a chemical reaction requires accounting for the concerted movements of numerous solvent molecules and ions with their own fast or slow timescales. While widely used static Density Functional Theory (DFT) calculations at 0 temperature can be beneficial for such investigations, they do not account for dynamics, and lack precision in describing the molecular environments. They particularly fail at correct, rigorous treatments of finite-temperature fluctuations, and thus generalization to experimentally relevant conditions. In PNAS Benayad et al develop a scalable, generalizable approach for designing Neural Network Potentials (NNPs) that can handle chemical reactivity in solvated systems with quantum accuracy at classical costs. Specifically, they study phosphoester bond formation and rupture, which is fundamentally relevant to the Phosphorus-Oxygen bond formation central to life, and especially for the RNA world hypothesis. The framework developed here has the potential to generalize to different chemical reactions of energy and biological relevance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗