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.

The uptake of trace gases from the air into aerosol particles impacts a wide range of environmental systems. Among other things, such multiphase processes help to determine the oxidative power of the atmosphere by acting as sinks for nitrogen oxides. Of particular long-standing interest is the reactive uptake of N2O5 in aqueous aerosol, which is estimated to account for 15–50% of the loss of NOx in the troposphere. Despite the significant study, basic questions remain concerning the mechanism of N2O5 uptake. Molecular dynamics simulations can be used to obtain a molecular perspective on gaseous uptake, free of underlying rate limitation assumptions. However, studying such processes theoretically imposes challenges, since uptake coefficients are exponentially sensitive to free energy differences and the simulations involve large systems and long times to model the complex dynamics. While qualitative predictions of mechanisms can be typically studied with conventional empirical force fields or density functional theory-based models, quantitative predictions require higher levels of accuracy. To address this challenge, a many-body potential, MB-nrg, has recently been parameterized from coupled-cluster calculations, providing the capability of making quantitative predictions of the thermodynamics and kinetics leading to the N2O5 uptake.

The validity of the traditional resistor model for the reactive uptake of N2O5 has been recently called into question due to the difficulty of reconciling the kinetics with field measurements, combined with theoretical work providing indications of interfacial stability and reactivity. The mechanism of uptake has been recently explored directly using a neural network-based reactive model, and it was found that interfacial rather than bulk phase processes dictate the observed uptake coefficient. Using training data obtained from density functional theory, this study found that the hydrolysis rate was sufficiently fast at the interface that bulk phase partitioning cannot kinetically compete, and the uptake was determined by a competition between interfacial hydrolysis and evaporation. These calculations found modest agreement with experimental uptake coefficient values, consistent with the expected qualitative accuracy of the model employed. As direct experimental confirmation of the importance of the interface is difficult, an alternative means of validating it is to employ models with higher chemical accuracy. This is the aim of the current work, to apply a quantitatively accurate potential to extract the thermodynamic and kinetic properties underpinning the uptake of N2O5 into water.