Ab Initio Molecular Dynamics Simulations of the Liquid/Vapor Interface of Sulfuric Acid Solutions

Chemistrees: Data-Driven Identification of Reaction Pathways via Machine Learning

We propose to analyze molecular dynamics (MD) output via a supervised machine learning (ML) algorithm, the decision tree. The approach aims to identify the predominant geometric features which correlate with trajectories that transition between two arbitrarily defined states. The data-driven algorithm aims to identify these features without the bias of human “chemical intuition”. We demonstrate the method by analyzing the proton exchange reactions in formic acid solvated in small water clusters. The simulations were performed with ab initio MD combined with a method to efficiently sample the rare event, path sampling. Our ML analysis identified relevant geometric variables involved in the proton transfer reaction and how they may change as the number of solvating water molecules changes.

Atmospheric Spectroscopy and Photochemistry at Environmental Water Interfaces

The air-water interface is ubiquitous in nature, as manifested in the form of the surfaces of oceans, lakes, and atmospheric aerosols. The aerosol interface, in particular, can play a crucial role in atmospheric chemistry. The adsorption of atmospheric species onto and into aerosols modifies their concentrations and chemistries. Moreover, the aerosol phase allows otherwise unlikely solution-phase chemistry to occur in the atmosphere. The effect of the air-water interface on these processes is not entirely known. This review summarizes recent theoretical investigations of the interactions of atmosphere species with the air-water interface, including reactant adsorption, photochemistry, and the spectroscopy of reactants at the water surface, with an emphasis on understanding differences between interfacial chemistries and the chemistries in both bulk solution and the gas phase. The results discussed here enable an understanding of fundamental concepts that lead to potential air-water interface effects, providing a framework to understand the effects of water surfaces on our atmosphere.

Ab initio molecular dynamics studies of formic acid dimer colliding with liquid water

Ab initio molecular dynamics simulations of formic acid (FA) dimer colliding with liquid water at 300 K have been performed using density functional theory. The two energetically lowest FA dimer isomers were collided with a water slab at thermal and high kinetic energies up to 68kBT. Our simulations agree with recent experimental observations of nearly a complete uptake of gas-phase FA dimer: the calculated average kinetic energy of the dimers immediately after collision is 5 ± 4% of the incoming kinetic energy, which compares well with the experimental value of 10%. Simulations support the experimental observation of no delayed desorption of FA dimers following initial adsorption. Our analysis shows that the FA dimer forms hydrogen bonds with surface water molecules, where the hydrogen bond order depends on the dimer structure, such that the most stable isomer possesses fewer FA-water hydrogen bonds than the higher energy isomer. Nevertheless, even the most stable isomer can attach to the surface through one hydrogen bond despite its reduced hydrophilicity. Our simulations further show that the probability of FA dimer dissociation is increased by high collision energies, the dimer undergoes isomerization from the higher energy to the lowest energy isomer, and concerted double-proton transfer occurs between the FA monomers. Interestingly, proton transfer appears to be driven by the release of energy arising from such isomerization, which stimulates those internal vibrational degrees of freedom that overcome the barrier of a proton transfer.

Microstructures of the Sulfonic Acid-Functionalized Ionic Liquid/Sulfuric Acid and Their Interactions: A Perspective from the Isobutane Alkylation

The all-atom force field for concentrated sulfuric acid (98.30 wt %) was developed in this work based on ab initio calculations. The structural and dynamical properties of sulfuric acid and the mixing behaviors of sulfuric acid with ionic liquids (ILs), i.e., SFIL (1-methyl-3-(propyl-3-sulfonate) imidazolium bisulfate ([PSMim][HSO₄])) and non-SFIL (1-methyl-3-propyl imidazolium bisulfate ([PMim][HSO₄])), were investigated using a molecular dynamics simulation. For sulfuric acid, most H₃O⁺ ions were found beside HSO₄^- ions, forming a contact ion pair with the HSO₄^- ions, and three-dimensional hydrogen-bonding networks existed in the sulfuric acid. Analyses indicate that both ILs could be miscible with sulfuric acid with a strong exothermic character. The new strong interaction site between the sulfonic acid group of SFIL and an H₂SO₄ molecule through a strong hydrogen-bonding interaction was observed, which was beneficial to the catalytic activity and stability of the sulfuric acid. This observation is in good agreement with the experimental results that indicate SFILs could enhance the reusability of sulfuric acid for the isobutane alkylation about 4-fold compared to that of non-SFILs. Hopefully this work will provide insights into the screening and designing of new isobutane alkylation catalysts based on sulfuric acid and SFILs.

Sulphur Kβ emission spectra reveal protonation states of aqueous sulfuric acid

In this paper we report an X-ray emission study of bulk aqueous sulfuric acid. Throughout the range of molarities from 1 M to 18 M the sulfur Kβ emission spectra from H2SO4 (aq) depend on the molar fractions and related deprotonation of H2SO4. We compare the experimental results with results from emission spectrum calculations based on atomic structures of single molecules and structures from ab initio molecular dynamics simulations. We show that the S Kβ emission spectrum is a sensitive probe of the protonation state of the acid molecules. Using non-negative matrix factorization we are able to extract the fractions of different protonation states in the spectra, and the results are in good agreement with the simulation for the higher part of the concentration range.

Deprotonation of formic acid in collisions with a liquid water surface studied by molecular dynamics and metadynamics simulations

Formic acid has a lower barrier to deprotonation at the air–water interface than in bulk liquid water.

Computational studies of atmospherically-relevant chemical reactions in water clusters and on liquid water and ice surfaces

CONSPECTUS: Reactions on water and ice surfaces and in other aqueous media are ubiquitous in the atmosphere, but the microscopic mechanisms of most of these processes are as yet unknown. This Account examines recent progress in atomistic simulations of such reactions and the insights provided into mechanisms and interpretation of experiments. Illustrative examples are discussed. The main computational approaches employed are classical trajectory simulations using interaction potentials derived from quantum chemical methods. This comprises both ab initio molecular dynamics (AIMD) and semiempirical molecular dynamics (SEMD), the latter referring to semiempirical quantum chemical methods. Presented examples are as follows: (i) Reaction of the (NO(+))(NO3(-)) ion pair with a water cluster to produce the atmospherically important HONO and HNO3. The simulations show that a cluster with four water molecules describes the reaction. This provides a hydrogen-bonding network supporting the transition state. The reaction is triggered by thermal structural fluctuations, and ultrafast changes in atomic partial charges play a key role. This is an example where a reaction in a small cluster can provide a model for a corresponding bulk process. The results support the proposed mechanism for production of HONO by hydrolysis of NO2 (N2O4). (ii) The reactions of gaseous HCl with N2O4 and N2O5 on liquid water surfaces. Ionization of HCl at the water/air interface is followed by nucleophilic attack of Cl(-) on N2O4 or N2O5. Both reactions proceed by an SN2 mechanism. The products are ClNO and ClNO2, precursors of atmospheric atomic chlorine. Because this mechanism cannot result from a cluster too small for HCl ionization, an extended water film model was simulated. The results explain ClNO formation experiments. Predicted ClNO2 formation is less efficient. (iii) Ionization of acids at ice surfaces. No ionization is found on ideal crystalline surfaces, but the process is efficient on isolated defects where it involves formation of H3O(+)-acid anion contact ion pairs. This behavior is found in simulations of a model of the ice quasi-liquid layer corresponding to large defect concentrations in crystalline ice. The results are in accord with experiments. (iv) Ionization of acids on wet quartz. A monolayer of water on hydroxylated silica is ordered even at room temperature, but the surface lattice constant differs significantly from that of crystalline ice. The ionization processes of HCl and H2SO4 are of high yield and occur in a few picoseconds. The results are in accord with experimental spectroscopy. (v) Photochemical reactions on water and ice. These simulations require excited state quantum chemical methods. The electronic absorption spectrum of methyl hydroperoxide adsorbed on a large ice cluster is strongly blue-shifted relative to the isolated molecule. The measured and calculated adsorption band low-frequency tails are in agreement. A simple model of photodynamics assumes prompt electronic relaxation of the excited peroxide due to the ice surface. SEMD simulations support this, with the important finding that the photochemistry takes place mainly on the ground state. In conclusion, dynamics simulations using quantum chemical potentials are a useful tool in atmospheric chemistry of water media, capable of comparison with experiment.

First and second deprotonation of H2SO4on wet hydroxylated (0001) α-quartz

We present anab initiomolecular dynamics study of deprotonation of sulfuric acid on wet quartz, a topic of atmospheric interest.

NOx Reactions on Aqueous Surfaces with Gaseous HCl: Formation of a Potential Precursor to Atmospheric Cl Atoms

Chlorine atoms are highly reactive free radicals known to catalyze ozone depletion in the stratosphere and organic oxidation in the troposphere. They are readily produced photolytically upon irradiation of some stable Cl containing species, for instance, nitrosyl chloride, ClNO. We predict the formation of ClNO using ab initio molecular dynamics (AIMD) simulations of an NO2 dimer on the surface of a thin film of water upon which gaseous HCl impinges. The reactant is chloride ion formed when HCl ionizes on the water film. The same mechanism for ClNO production may occur in humid environments when ONONO2 (the asymmetric NO2 dimer examined here) comes in contact with either HCl or sea salt. The film of water serves to (1) stabilize ONONO2 on the film surface so that it is localized and physically accessible for reaction, (2) provide the medium to ionize HCl, and (3) activate ONONO2 making it more susceptible to nucleophilic attack by chloride. This substitution/elimination mechanism is new for NOx chemistry on thin water films and could not be derived from studies on small clusters.