Dynamics of collisions and uptake of alcohol molecules with hydrated nitric acid clusters.

We investigate the collisions of different alcohol molecules with hydrated nitric acid clusters using a molecular beam experiment and molecular dynamics simulations. The uptake cross sections σp for the molecules evaluated from the experiment are in excellent agreement with the simulations. This suggests that (i) the nontrivial assumptions implemented in the evaluation procedure of the experimental data are valid, and (ii) the simulations describe correctly the major processes in the molecule-cluster collisions. We observe that σp decreases with the increasing alkyl chain length of the alcohol, and also with the branching of the molecules that have the same mass but different structures. These systematic trends can be rationalized based on the accessibility of the hydrophilic OH group, which decreases with the increasing chain length and steric hindrance. The observed trends and their interpretation differ significantly from the simple model of hard-sphere collisions. The obtained data shall be beneficial not only for the fundamental understanding of the molecule-cluster collisions, but also in the modelling of atmospheric new-particle formation and aerosol growth.

Heterogeneous Reactions of Methane with Cl Radicals on Large ArN Clusters.

Heterogeneous chemistry on the surfaces of atmospheric particles has a wide impact on the properties and composition of the Earth's atmosphere. In laboratory studies, clusters can represent proxies to atmospheric aerosols and help to discern the individual steps in reactions on or in aerosols. We investigate the reactivity of Cl and CCl3 radicals with methane on argon clusters using the pickup method. For radical generation, we built a new pyrolysis source partially adapting the design of radical sources that utilize the supersonic expansion into a heated silicon carbide tube. Large ArN, N̅ ≈ 110, clusters were generated in a supersonic expansion, and CH4 molecules were embedded in the clusters via a pickup process followed by the uptake of the radicals produced in the pyrolysis source. The analysis of the mass spectra recorded under different experimental conditions (i.e., with the pyrolysis ON and OFF and with only one or both reactants) allowed us to identify various products of the radical reactions on ArN. We propose a sequence of reactions based on the reaction energetics. It starts with the hydrogen abstraction from CH4 by a Cl radical resulting in HCl and CH3 followed by a halogenation step where CCl4 molecules react with the available CH3 radicals, yielding CH3Cl. By analogy, the CH3Cl enters another hydrogen abstraction by Cl, producing HCl and the CH2Cl radical, which again undergoes a halogenation step with CCl4, generating CH2Cl2. Further reaction of CH2Cl2 with Cl terminates the sequence by the production of HCl and CHCl2.

Pickup and reactions of molecules on clusters relevant for atmospheric and interstellar processes.

In this perspective, we review experiments with molecules picked up on large clusters in molecular beams with the focus on the processes in atmospheric and interstellar chemistry. First, we concentrate on the pickup itself, and we discuss the pickup cross sections. We measure the uptake of different atmospheric molecules on mixed nitric acid-water clusters and determine the accommodation coefficients relevant for aerosol formation in the Earth's atmosphere. Then the coagulation of the adsorbed molecules on the clusters is investigated. In the second part of this perspective, we review examples of different processes triggered by UV-photons or electrons in the clusters with embedded molecules. We start with the photodissociation of hydrogen halides and Freon CF2Cl2 on ice nanoparticles in connection with the polar stratospheric ozone depletion. Next, we mention reactions following the excitation and ionization of the molecules adsorbed on clusters. The first ionization-triggered reaction observed between two different molecules picked up on the cluster was the proton transfer between methanol and formic acid deposited on large argon clusters. Finally, negative ion reactions after slow electron attachment are illustrated by two examples: mixed nitric acid-water clusters, and hydrogen peroxide deposited on large ArN and (H2O)N clusters. The selected examples are discussed from the perspective of the atmospheric and interstellar chemistry, and several future directions are proposed.

Ion and radical chemistry in (H2O2)N clusters.

We investigate the ionization induced chemistry of hydrogen peroxide in (H2O2)N clusters generated after the pickup of individual H2O2 molecules on large free ArM, M[combining macron]≈ 160, nanoparticles in molecular beams. Positive and negative ion mass spectra are recorded after an electron ionization of the clusters at energies 5-70 eV and after a slow electron attachment (below 4 eV), respectively. The spectra demonstrate that (H2O2)N clusters with N≥ 20 are formed on argon nanoparticles. This is the first experimental report on hydrogen peroxide clusters in molecular beams. The major negative cluster ion series (H2O2)nO2- indicates O2- ion formation. The dissociative electron attachment to H2O2 molecules in the gas phase yielded only OH- and O- (Nandi et al., Chem. Phys. Lett., 2003, 373, 454). These ions and the series containing them are much less abundant in the clusters. We propose a sequence of ion-molecule and radical reactions to explain the formation of O2-, HO2- and other ions observed in the negatively charged cluster ion series. Since hydrogen peroxide plays an important role in many areas of chemistry from the Earth's atmosphere to biological tissues, our study opens new horizons for experimental investigations of hydrogen peroxide chemistry in complex environments.

Proton Transfer Reactions between Methanol and Formic Acid Deposited on Free ArN Nanoparticles.

We have sequentially picked up two astrochemically relevant Brønsted acids (methanol and formic acid) on the surface of argon nanoparticles acting as a cold support. Photoionization and electron ionization yield (HCOOH)xH+, (CH3OH)xH+, and mixed protonated clusters. Experiments with perdeuterated methanol CD3OD demonstrate notable proton transfer (PT) to formic acid acting as a proton acceptor in addition to the PT from formic acid which is, perhaps, a more intuitive one. We, therefore, for the first time observed reactions between two different complex molecules adsorbed individually on argon nanoparticles. The experimental results are compared with state-of-the-art quantum chemistry calculations, showing that both CH3OH•+ and HCOOH•+ radical cations resulting from ionization can act as efficient proton donors and neutral CH3OH and HCOOH as proton acceptors. According to the theoretical calculations, the C-H bond cleavage in the radical cation should be more favorable than the O-H bond cleavage. Both channels are observed and distinguished in the experiments with CD3OH and CH3OD. Our detailed mechanism of formation of the CH3O• and CH2OH• radicals contributes to understanding of astrochemistry in the protostellar medium.

Rapid Water Transport through Organic Layers on Ice.

Processes involving atmospheric aerosol and cloud particles are affected by condensation of organic compounds that are omnipresent in the atmosphere. On ice particles, organic compounds with hydrophilic functional groups form hydrogen bonds with the ice and orient their hydrophobic groups away from the surface. The organic layer has been expected to constitute a barrier to gas uptake, but recent experimental studies suggest that the accommodation of water molecules on ice is only weakly affected by condensed short-chain alcohol layers. Here, we employ molecular dynamics simulations to study the water interactions with n-butanol covered ice at 200 K and show that the small effect of the condensed layer is due to efficient diffusion of water molecules along the surface plane while seeking appropriate sites to penetrate, followed by penetration driven by the combined attractive forces from butanol OH groups and water molecules within the ice. The water molecules that penetrate through the n-butanol layer become strongly bonded by approximately three hydrogen bonds at the butanol-ice interface. The obtained accommodation coefficient (0.81 ± 0.03) is in excellent agreement with results from previous environmental molecular beam experiments, leading to a picture where an adsorbed n-butanol layer does not alter the apparent accommodation coefficient but dramatically changes the detailed molecular dynamics and kinetics.

Calcium ions in aqueous solutions: Accurate force field description aided by ab initio molecular dynamics and neutron scattering.

We present a combination of force field and ab initio molecular dynamics simulations together with neutron scattering experiments with isotopic substitution that aim at characterizing ion hydration and pairing in aqueous calcium chloride and formate/acetate solutions. Benchmarking against neutron scattering data on concentrated solutions together with ion pairing free energy profiles from ab initio molecular dynamics allows us to develop an accurate calcium force field which accounts in a mean-field way for electronic polarization effects via charge rescaling. This refined calcium parameterization is directly usable for standard molecular dynamics simulations of processes involving this key biological signaling ion.

Modeling photoionization of aqueous DNA and its components.

Radiation damage to DNA is usually considered in terms of UVA and UVB radiation. These ultraviolet rays, which are part of the solar spectrum, can indeed cause chemical lesions in DNA, triggered by photoexcitation particularly in the UVB range. Damage can, however, be also caused by higher energy radiation, which can ionize directly the DNA or its immediate surroundings, leading to indirect damage. Thanks to absorption in the atmosphere, the intensity of such ionizing radiation is negligible in the solar spectrum at the surface of Earth. Nevertheless, such an ionizing scenario can become dangerously plausible for astronauts or flight personnel, as well as for persons present at nuclear power plant accidents. On the beneficial side, ionizing radiation is employed as means for destroying the DNA of cancer cells during radiation therapy. Quantitative information about ionization of DNA and its components is important not only for DNA radiation damage, but also for understanding redox properties of DNA in redox sensing or labeling, as well as charge migration along the double helix in nanoelectronics applications. Until recently, the vast majority of experimental and computational data on DNA ionization was pertinent to its components in the gas phase, which is far from its native aqueous environment. The situation has, however, changed for the better due to the advent of photoelectron spectroscopy in liquid microjets and its most recent application to photoionization of aqueous nucleosides, nucleotides, and larger DNA fragments. Here, we present a consistent and efficient computational methodology, which allows to accurately evaluate ionization energies and model photoelectron spectra of aqueous DNA and its individual components. After careful benchmarking, the method based on density functional theory and its time-dependent variant with properly chosen hybrid functionals and polarizable continuum solvent model provides ionization energies with accuracy of 0.2-0.3 eV, allowing for faithful modeling and interpretation of DNA photoionization. The key finding is that the aqueous medium is remarkably efficient in screening the interactions within DNA such that, unlike in the gas phase, ionization of a base, nucleoside, or nucleotide depends only very weakly on the particular DNA context. An exception is the electronic interaction between neighboring bases which can lead to sequence-specific effects, such as a partial delocalization of the cationic hole upon ionization enabled by presence of adjacent bases of the same type.

Oxidation half-reaction of aqueous nucleosides and nucleotides via photoelectron spectroscopy augmented by ab initio calculations.

Oxidative damage to DNA and hole transport between nucleobases in oxidized DNA are important processes in lesion formation for which surprisingly poor thermodynamic data exist, the relative ease of oxidizing the four nucleobases being one such example. Theoretical simulations of radiation damage and charge transport in DNA depend on accurate values for vertical ionization energies (VIEs), reorganization energies, and standard reduction potentials. Liquid-jet photoelectron spectroscopy can be used to directly study the oxidation half-reaction. The VIEs of nucleic acid building blocks are measured in their native buffered aqueous environment. The experimental investigation of purine and pyrimidine nucleotides, nucleosides, pentose sugars, and inorganic phosphate demonstrates that photoelectron spectra of nucleotides arise as a spectral sum over their individual chemical components; that is, the electronic interactions between each component are effectively screened from one another by water. Electronic structure theory affords the assignment of the lowest energy photoelectron band in all investigated nucleosides and nucleotides to a single ionizing transition centered solely on the nucleobase. Thus, combining the measured VIEs with theoretically determined reorganization energies allows for the spectroscopic determination of the one-electron redox potentials that have been difficult to establish via electrochemistry.

Ion pairing in aqueous lithium salt solutions with monovalent and divalent counter-anions.

Molecular dynamics simulations of concentrated aqueous solutions of LiCl and Li2SO4 were conducted in order to provide molecular insight into recent neutron scattering data. The structures predicted from the molecular dynamics simulations using standard nonpolarizable force fields provided a very poor fit to the experiment; therefore, refinement was needed. The electronic polarizability of the medium was effectively accounted for by implementing the electronic continuum correction, which practically means rescaling the ionic charges. Consistent with previous studies, we found that this approach in each case provided a significantly improved fit to the experimental data, which was further enhanced by slightly adjusting the radius of the lithium ion. The polarization effect was particularly pronounced in the Li2SO4 solution where the ions in the nonpolarizable simulations tended to cluster unphysically. With the above alterations, the employed force field displayed an excellent fit to the neutron scattering data and provided a useful interpretative framework for the experimental measurements. At the same time, the present study underlines the importance of solvent polarization effects in hydration of ions with high charge density.

Peptide salt bridge stability: from gas phase via microhydration to bulk water simulations.

The salt bridge formation and stability in the terminated lysine-glutamate dipeptide is investigated in water clusters of increasing size up to the limit of bulk water. Proton transfer dynamics between the acidic and basic side chains is described by DFT-based Born-Oppenheimer molecular dynamics simulations. While the desolvated peptide prefers to be in its neutral state, already the addition of a single water molecule can trigger proton transfer from the glutamate side chain to the lysine side chain, leading to a zwitterionic salt bridge state. Upon adding more water molecules we find that stabilization of the zwitterionic state critically depends on the number of hydrogen bonds between side chain termini, the water molecules, and the peptidic backbone. Employing classical molecular dynamics simulations for larger clusters, we observed that the salt bridge is weakened upon additional hydration. Consequently, long-lived solvent shared ion pairs are observed for about 30 water molecules while solvent separated ion pairs are found when at least 40 or more water molecules hydrate the dipeptide. These results have implications for the formation and stability of salt bridges at partially dehydrated surfaces of aqueous proteins.

Uptake of Hydrogen Bonding Molecules by Benzene Nanoparticles.

The uptake of molecules on nanometer-size clusters of polyaromatic hydrocarbons (PAHs) is important for the condensation of water on PAH aerosols in the atmosphere and for ice mantle growth on nanoparticles in the interstellar medium. We generate benzene clusters BzN of mean size N̅ ≈ 300 (radius R̅ ≈ 2.2 Å) as a model system for the PAH nanoparticles. Using molecular beams and mass spectrometry detection, we investigate the uptake of water, methanol, and ethanol by these clusters. All picked up molecules are highly mobile on BzN and generate clusters within <3 ms. The relative uptakes for the different investigated molecules can be directly compared and quantified. Water molecules exhibit the lowest relative pickup probability that is ∼30% lower than those for methanol and ethanol, which are approximately the same.

Chasing charge localization and chemical reactivity following photoionization in liquid water.

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within ~30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H(2)O(+) + H(2)O → OH + H(3)O(+). The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM∕MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with ~40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H(2)O(+) ((aq)) species.

Protein Preferential Solvation in Water:Glycerol Mixtures.

For proteins in solvent mixtures, the relative abundances of each solvent in their solvation shell have a critical impact on their properties. Preferential solvation of a series of proteins in water-glycerol mixtures is studied here over a broad range of solvent compositions via classical molecular dynamics simulations. Our simulation results reveal that the differences between shell and bulk compositions exhibit dramatic changes with solvent composition, temperature, and protein nature. In contrast with the simple and widely used picture where glycerol is completely excluded from the protein interface, we show that for aqueous solutions with less than 50% glycerol in volume, protein solvation shells have approximately the same composition as the bulk solvent and proteins are in direct contact with glycerol. We further demonstrate that at high glycerol concentration, glycerol depletion from the solvation shell is due to an entropic factor arising from the reduced accessibility of bulky glycerol molecules in protein cavities. The resulting molecular picture is important to understand protein activity and cryopreservation in mixed aqueous solvents.

Oxidation Enhances Aerosol Nucleation: Measurement of Kinetic Pickup Probability of Organic Molecules on Hydrated Acid Clusters.

We investigate the uptake of the most prominent biogenic volatile organic compounds (VOCs)-isoprene, α-pinene, and their selected oxidation products-by hydrated acid clusters in a molecular beam experiment and by DFT calculations. Our experiments provide a unique and direct way of determination of the surface accommodation coefficient (αS) on the proxies of ultrafine aerosol particles. Since we are able to determine unambiguously the fraction of the clusters to which the molecules stick upon collisions, our αS is a purely kinetic parameter disentangling the molecule pickup from its evaporation. Oxidation increases the αS of VOCs by more than an order of magnitude, because oxidized compounds form hydrogen bonds with the clusters, whereas the interactions of the parent VOCs are weaker and nonspecific. This work provides molecular-level insight into the condensation of single molecules into atmospheric particles, which has important implications for aerosol nucleation and growth.

Structure and Dynamics of the Hydration Shell: Spatially Decomposed Time Correlation Approach.

Molecular simulations provide insight into solvation structures and dynamics with unparalleled spatial and temporal resolution. Here, we take advantage of this fact and develop a set of generally applicable computational tools for a detailed analysis of the hydration shell around an ionic or molecular solute. These tools allow us to quantify and visualize orientationally resolved radial distribution functions as well as distance-resolved orientational time-correlation functions of water molecules surrounding the solute. Such a detailed view of the hydration shells allows us to unravel important structural and dynamical features, which are not accessible when employing standard analysis techniques.

Size and Origins of Long-Range Orientational Water Correlations in Dilute Aqueous Salt Solutions.

Long-range ordering of water around solutes has been repeatedly invoked as the key to its biological function. Recently, it has been shown that in an 8 mM aqueous NaCl solution the orientational correlation between water molecules extends beyond 8 nm. This was interpreted as arising from ion-induced long-range effects on the water collective hydrogen-bond interactions. Each ion was suggested to affect >10 000 water molecules, leading to a picture involving nanoscopic "ordered domains". Using molecular dynamics simulations, we show that the very small long-range tail in the correlation function is caused primarily by pairs of water molecules belonging to different ions' hydration shells and can be made to practically disappear by rearranging the ionic positions. This means that the ion-induced water orientational ordering in millimolar salt solutions cannot be separated from ion-ion interaction effects, for which the Debye-Hückel screening length shrinks to less than 1 nm at physiological ionic strengths.

Coupled Valence-Bond State Molecular Dynamics Description of an Enzyme-Catalyzed Reaction in a Non-Aqueous Organic Solvent.

Enzymes are widely used in nonaqueous solvents to catalyze non-natural reactions. While experimental measurements showed that the solvent nature has a strong effect on the reaction kinetics, the molecular details of the catalytic mechanism in nonaqueous solvents have remained largely elusive. Here we study the transesterification reaction catalyzed by the paradigm subtilisin Carlsberg serine protease in an organic apolar solvent. The rate-limiting acylation step involves a proton transfer between active-site residues and the nucleophilic attack of the substrate to form a tetrahedral intermediate. We design the first coupled valence-bond state model that simultaneously describes both reactions in the enzymatic active site. We develop a new systematic procedure to parametrize this model on high-level ab initio QM/MM free energy calculations that account for the molecular details of the active site and for both substrate and protein conformational fluctuations. Our calculations show that the reaction energy barrier changes dramatically with the solvent and protein conformational fluctuations. We find that the mechanism of the tetrahedral intermediate formation during the acylation step is similar to that determined under aqueous conditions, and that the proton transfer and nucleophilic attack reactions occur concertedly. We identify the reaction coordinate to be mostly due to the rearrangement of some residual water molecules close to the active site.

Dependence of the Rate of LiF Ion-Pairing on the Description of Molecular Interaction.

We present an analysis of the dynamics of ion-pairing of lithium fluoride (LiF) in aqueous solvent using both detailed molecular simulation as well as reduced models within a generalized Langevin equation (GLE) framework. We explored the sensitivity of the ion-pairing phenomena to the details of descriptions of molecular interaction, comparing two empirical potentials to explicit quantum based density functional theory. We find quantitative differences in the potentials of mean force for ion-pairing as well as time dependent frictions that lead to variations in the rate constant and reactive flux correlation functions. These details reflect differences in solvent response to ion-pairing between different representations of molecular interaction and influence anharmonicity of the dynamic response. We find that the short-time anharmonic response is recovered with a GLE parametrization. Recovery of the details of long time response may require extensions to the reduced model. We show that the utility of using a reduced model leads to a straightforward application of variational transition state theory concepts to the condensed phase system. The significance of this is reflected in the analysis of committor distributions and the variation of planar hypersurfaces, leading to an improved understanding of factors that determine the rate of LiF ion-pairing.

Exploring Ion-Ion Interactions in Aqueous Solutions by a Combination of Molecular Dynamics and Neutron Scattering.

Recent advances in computational and experimental techniques have allowed for accurate description of ion pairing in aqueous solutions. Free energy methods based on ab initio molecular dynamics, as well as on force fields accounting effectively for electronic polarization, can provide quantitative information about the structures and occurrences of individual types of ion pairs. When properly benchmarked against electronic structure calculations for model systems and against structural experiments, in particular neutron scattering, such force field simulations represent a powerful tool for elucidating interactions of salt ions in complex biological aqueous environments.