Free energy study of H2O, N2O5, SO2, and O3 gas sorption by water droplets/slabs

Reparameterization of Polarizable Force Fields for Studying Ion Transfer across Liquid-Liquid Interfaces

We have developed a general scheme for refining classical polarizable molecular dynamics (MD) force fields that can accurately describe the molecular interactions in systems with liquid-liquid interfaces. While ab initio MD (AIMD) simulations can naturally describe molecular interactions, they are often so computationally expensive that simulating large system sizes and/or long time scales is usually infeasible. To resolve this, we parameterized efficient and accurate classical polarizable force fields that use AIMD reference data by minimizing both the relative entropy and the root mean squared deviation in atomic forces. We utilized our new multiscale models to study chloride ion transfer across the water-dichloromethane (DCM) interface with and without the tetraethylammonium cation as the phase-transfer catalyst. Our calculated free-energy barrier for the water-DCM interface is consistent with the other reported simulation results. We further analyzed the ion-transfer process by studying the hydration shell structures around the chloride ion and the ion-pair formation to better understand the mechanism. We observed that electronic polarizability is an important factor for the studied phase-transfer catalyst to lower the free-energy barrier of the ion transfer. Using the water-benzene interface system as an additional example, we show that our parameterization scheme provides a general route for modeling different liquid-liquid interface systems even when the experimental data or force field parameters are not readily available.

Flexible Polarizable Water Model Parameterized via Gaussian Process Regression

Water is one of the most common components in molecular dynamics (MD) simulations. Using Gaussian process regression for predicting the properties of a water model without the need of running a simulation whenever the parameters are changed, we obtained a flexible polarizable water model, named SWM4/Fw, that is able to reproduce many reference water properties. The added flexibility is critical for modeling chemical reactions in which chemical bonds can be stretched or even broken and for directly calculating vibrational spectra. In addition to being one of the few water models that are both flexible and polarizable, SWM4/Fw is also efficient thanks to the extended Lagrangian scheme with Drude oscillators. The overall accuracy is on par with or better than the related SWM4-NDP model.

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O₃. However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. To examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients (trans-aconitic acid (C₆H₆O₆), maleic acid (C₄H₄O₄), and sodium nitrite (NaNO₂)). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O₃ partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O₃ and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O₃, compared to gas and bulk phases is rapidly maintained and remains constant even as O₃ is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O₃ loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O₃ arrival in the bulk by diffusion. As a result, bulk-phase O₃ becomes depleted from its Henry's law solubility. This bulk-phase O₃ depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k_rxn ≈ 10³-10⁴ M^-1 s^-1), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k_rxn ≈ 10⁵ M^-1 s^-1) where surface reactions strongly impact the observed decay kinetics.

Uptake of N2O5 by aqueous aerosol unveiled using chemically accurate many-body potentials

The reactive uptake of N2O5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N2O5 in water. The free energy profile highlights that N2O5 is selectively adsorbed to the liquid-vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction-diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction-diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N2O5 in more complex solutions.

Highly Accurate Many-Body Potentials for Simulations of N2O5 in Water: Benchmarks, Development, and Validation

Dinitrogen pentoxide (N₂O₅) is an important intermediate in the atmospheric chemistry of nitrogen oxides. Although there has been much research, the processes that govern the physical interactions between N₂O₅ and water are still not fully understood at a molecular level. Gaining a quantitative insight from computer simulations requires going beyond the accuracy of classical force fields while accessing length scales and time scales that are out of reach for high-level quantum-chemical approaches. To this end, we present the development of MB-nrg many-body potential energy functions for nonreactive simulations of N₂O₅ in water. This MB-nrg model is based on electronic structure calculations at the coupled cluster level of theory and is compatible with the successful MB-pol model for water. It provides a physically correct description of long-range many-body interactions in combination with an explicit representation of up to three-body short-range interactions in terms of multidimensional permutationally invariant polynomials. In order to further investigate the importance of the underlying interactions in the model, a TTM-nrg model was also devised. TTM-nrg is a more simplistic representation that contains only two-body short-range interactions represented through Born-Mayer functions. In this work, an active learning approach was employed to efficiently build representative training sets of monomer, dimer, and trimer structures, and benchmarks are presented to determine the accuracy of our new models in comparison to a range of density functional theory methods. By assessing the binding curves, distortion energies of N₂O₅, and interaction energies in clusters of N₂O₅ and water, we evaluate the importance of two-body and three-body short-range potentials. The results demonstrate that our MB-nrg model has high accuracy with respect to the coupled cluster reference, outperforms current density functional theory models, and thus enables highly accurate simulations of N₂O₅ in aqueous environments.

Understanding the properties of methanesulfinic acid at the air-water interface

Methanesulfinic acid (MSIA), an organic sulfur compound, is mainly produced in the oxidation process of dimethyl sulfide in the atmosphere. The properties of MSIA at the air-water interface were studied using molecular dynamics (MD) simulations. The result shows that the lowest system free energy is located at the interface. Because the free energy difference between the interface and water phase is 3.2 kJ mol^-1, the MSIA molecule can easily get out of the free energy well and travel to water phase by the thermal motion, leading to only a 21% probability of its occurrence at the interface. The MSIA molecule tends to tilt at the interface with the sulfino group (-S(O)-OH) pointing toward the water phase. The feature of hydration status at the air-water interface may be favorable to the heterogeneous oxidation of MSIA.

Mass Transfer Thermodynamics through a Gas-Liquid Interface

Molecular level information about thermodynamic variations (enthalpy, entropy, and free energy) of a gas molecule as it crosses a gas-liquid interface is strongly lacking from an experimental perspective under equilibrium conditions. Herein, we perform in situ measurements of water interacting with the ionic liquid (IL) 1-butyl-3-methylimidazolium acetate, [C₄mim][Ace], using ambient pressure X-ray photoelectron spectroscopy in order to assess the interfacial uptake of water quantitatively as a function of temperature, pressure, and water mole fraction ( x_w). The surface spectroscopy results are compared to existing bulk water absorption experiments, showing that the amount of water in the interfacial region is consistently greater than that in the bulk. The enthalpy and entropy of water sorption vary significantly between the gas-liquid interface and the bulk as a function of x_w, with a crossover that occurs near x_w = 0.6 where the water-IL mixture converts from being homogeneous ( x_w < 0.6) to nanostructured ( x_w > 0.6). Free energy results reveal that water at the gas-IL interface is thermodynamically more favorable than that in the bulk, consistent with the enhanced water concentration in the interfacial region. The results herein show that the efficacy for an ionic liquid to absorb a gas phase molecule is not merely a function of bulk solvation parameters but also is significantly influenced by the thermodynamics occurring across the gas-IL interface during the mass transfer process.