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

Deprotonation of organic acids at aqueous surfaces has important implications in atmospheric chemistry and other disciplines, yet it is not well-characterized or understood. This article explores the interactions of formic acid (FA), including ionization, in collisions at the air-water interface. Ab initio molecular dynamics simulations with dispersion-corrected density functional theory were used. The 8-50 picosecond duration trajectories all resulted in the adsorption of FA within the interfacial region, with no scattering, absorption into the bulk or desorption into the vapor. Despite the known weak acidity of FA, spontaneous deprotonation of the acid was observed at the interface on a broad picosecond timescale, ranging from a few picoseconds typical for stronger acids to tens of picoseconds. Deprotonation occurred in 4% of the trajectories, and was followed by Grotthuss proton transfer through adjacent water molecules. Both sequential and ultrafast concerted proton transfer were observed. The formation of contact ion pairs and solvent-separated ion pairs, and finally the reformation of neutral FA, both trans and cis conformers, occurred in different stages of the dynamics. To better understand the deprotonation mechanisms at the interface compared with the process in bulk water, we used well-tempered metadynamics to obtain deprotonation free energy profiles. While in bulk water FA deprotonation has a free energy barrier of 14.8 kJ mol(-1), in fair agreement with the earlier work, the barrier at the interface is only 7.5 kJ mol(-1). Thus, at the air-water interface, FA may dissociate more rapidly than in the bulk. This finding can be rationalized with reference to the dissimilar aqueous solvation and hydrogen-bonding environments in the interface compared to those in bulk liquid water.Peer reviewe

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 68k(B)T. 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.Peer reviewe

Catalytic Role for Water in the Atmospheric Production of ClNO

High level ab initio calculations of clusters comprised of water, HCl, and ON-ONO2 are used to study nitrosyl chloride (ClNO) formation in gas phase water clusters, which are also mimics for thin water films present at environmental interfaces. Two pathways are considered, direct formation from the reaction of gaseous HCl with ON-ONO2 and an indirect pathway involving the hydrolysis of ON-ONO2 to form HONO, followed by the reaction of HONO with HCl to form ClNO. Surprisingly, direct formation of ClNO is found to be the dominant channel in the presence of water despite the possibility of a competing hydrolysis of ON-ONO2 to form HONO. A single water molecule effectively catalyzes the ON-ONO2 + HCl reaction, and in the presence of two or more water molecules the reaction to form ClNO becomes spontaneous. Direct formation of ClNO is fast at room and ice temperatures, indicating the possible significance of this pathway for chlorine activation chemistry in both the polar and midlatitude troposphere, in volcanic plumes and indoors. The reaction enthalpies, activation energies, and rate constants for all studied reactions are reported. The results are discussed in light of recent experiments

Computed Pre-reactive Complex Association Lifetimes Explain Trends in Experimental Reaction Rates for Peroxy Radical Recombinations

The lifetimes of pre-reactive complexes, although implicitly part of the equations used to model many gas-phase bimolecular reactions, have seldom been included in quantitative calculations of rate coefficients. Here, we demonstrate the application of empirical molecular dynamics simulations of collisions between peroxy radicals to model association lifetimes. With the exception of the methyl peroxy−acetyl peroxy system, measurements of the lifetimes based on a phenomenological model are shown to correlate well with available experimental data for recombination reactions of peroxy radicals in cases where the rate-limiting transition state lies below the reactants in energy. Further, we predict reaction rates for larger α-pinene-derived peroxy radicals, and we interpret our results in tandem with available experimental data on these systems, which are of great relevance to improve our understanding of atmospheric aerosol formation.Peer reviewe

Structure of Large Nitrate−Water Clusters at Ambient Temperatures: Simulations with Effective Fragment Potentials and Force Fields with Implications for Atmospheric Chemistry

Structural properties of large NO3−·(H2O)n (n = 15−500) clusters are studied by Monte Carlo simulations using effective fragment potentials (EFPs) and by classical molecular dynamics simulations using a polarizable empirical force field. The simulation results are analyzed with a focus on the description of hydrogen bonding and solvation in the clusters. In addition, a comparison between the electronic structure based EFP and the classical force field description of the 32 water cluster system is presented. The EFP simulations, which focused on the cases of n = 15 and 32, show an internal, fully solvated structure and a “surface adsorbed” structure for the 32 water cluster at 300 K, with the latter configuration being more probable. The internal solvated structure and the “surface adsorbed” structure differ considerably in their hydrogen bonding coordination numbers. The force field based simulations agree qualitatively with these results, and the local geometry of NO3− and solvation at the surface-adsorbed site in the force field simulations are similar to those predicted using EFPs. Differences and similarities between the description of hydrogen bonding of the anion in the two approaches are discussed. Extensive classical force field based simulations at 250 K predict that long time scale stability of “internal” NO3−, which is characteristic of extended bulk aqueous interfaces, emerges only for n \u3e 300. Ab initio Møller−Plesset perturbation theory is used to test the geometries of selected surface and interior anions for n = 32, and the results are compared to the EFP and MD simulations. Qualitatively, all approaches agree that surface structures are preferred over the interior structures for clusters of this size. The relatively large aqueous clusters of NO3− studied here are of comparable size to clusters that lead to new particle formation in air. Nitrate ions on the surface of such clusters may have significantly different photochemistry than the internal species. The possible implications of surface-adsorbed nitrate ions for atmospheric chemistry are discussed

Hygroscopic Growth and Deliquescence of NaCl Nanoparticles Mixed with Surfactant SDS

Several complementary experimental and theoretical methodologies were used to explore water uptake on sodium chloride (NaCl) particles containing varying amounts of sodium dodecyl sulfate (SDS) to elucidate the interaction of water with well-defined, environmentally relevant surfaces. Experiments probed the hygroscopic growth of mixed SDS/NaCl nanoparticles that were generated by electrospraying aqueous 2 g/L solutions containing SDS and NaCl with relative NaCl/SDS weight fractions of 0, 5, 11, 23, or 50 wt/wt %. Particles with mobility-equivalent diameters of 14.0(±0.2) nm were size selected and their hygroscopic growth was monitored by a tandem nano-differential mobility analyzer as a function of relative humidity (RH). Nanoparticles generated from 0 and 5 wt/wt % solutions deliquesced abruptly at 79.1(±1.0)% RH. Both of these nanoparticle compositions had 3.1(±0.5) monolayers of adsorbed surface water prior to deliquescing and showed good agreement with the Brunauer−Emmett−Teller and the Frenkel−Halsey−Hill isotherms. Above the deliquescence point, the growth curves could be qualitatively described by Köhler theory after appropriately accounting for the effect of the particle shape on mobility. The SDS/NaCl nanoparticles with larger SDS fractions displayed gradual deliquescence at a RH that was significantly lower than 79.1%. All compositions of SDS/NaCl nanoparticles had monotonically suppressed mobility growth factors (GFm) with increasing fractions of SDS in the electrosprayed solutions. The Zdanovskii−Stokes−Robinson model was used to estimate the actual fractions of SDS and NaCl in the nanoparticles; it suggested the nanoparticles were enhanced in SDS relative to their electrospray solution concentrations. X-ray photoelectron spectroscopy (XPS), FTIR, and AFM were consistent with SDS forming first a monolayer and then a crystalline phase around the NaCl core. Molecular dynamics simulations of water vapor interacting with SDS/NaCl slabs showed that SDS kinetically hinders the initial water uptake. Large binding energies of sodium methyl sulfate (SMS)−(NaCl)4, H2O−(NaCl)4, and SMS−H2O−(NaCl)4 calculated at the MP2/cc-pVDZ level suggested that placing H2O in between NaCl and surfactant headgroup is energetically favorable. These results provide a comprehensive description of SDS/NaCl nanoparticles and their properties

Chemically-bound xenon in fibrous silica

High-level quantum chemical calculations reported here predict the existence and remarkable stability, of chemically-bound xenon atoms in fibrous silica. The results may support the suggestion of Sanloup and coworkers that chemically-bound xenon and silica account for the problem of "missing xenon" (by a factor of 20!) from the atmospheres of Earth and Mars. So far, the host silica was assumed to be quartz, which is in contradiction with theory. The xenon-fibrous silica molecule is computed to be stable well beyond room temperature. The calculated Raman spectra of the species agree well with the main features of the experiments by Sanloup et al. The results predict computationally the existence of a new family of noble-gas containing materials. The fibrous silica species are finite molecules, their laboratory preparation should be feasible, and potential applications are possible

When a proton attacks cellobiose in the gas phase: ab initio molecular dynamics simulations

Investigations of reaction pathways between a proton and cellobiose (CB), a glucose disaccharide of importance, were carried out in cis and trans CB using Ab Initio Molecular Dynamics (AIMD) simulations starting from optimized configurations where the proton is initially placed near groups with affinity for it. Near and above 300 K, protonated CB (H(+)CB) undergoes several transient reactions including charge transfer to the sugar backbone, water formation and dehydration, ring breaking and glycosidic bond breaking events as well as mutarotation and ring puckering events, all on a 10 ps timescale. cis H(+)CB is energetically favoured over trans H(+)CB in vacuo, with an energy gap larger than for the neutral CB