Discrete stages in the solvation and ionization of hydrogen chloride adsorbed on ice particles

Electric nuclear quadrupole coupling reveals dissociation of HCl with a few water molecules

Investigating the dissociation of acids in the presence of a limited number of water molecules is crucial for understanding various elementary chemical processes. In our study, focusing on HCl(H2O)n clusters (where HCl is hydrogen chloride and H2O is water) formed in a cold and isolated jet expansion, we used the nuclear quadrupole coupling tensor obtained through rotational spectroscopy to decipher the nature of the hydrogen-chlorine (H-Cl) chemical bond in a microaqueous environment. For n = 1 to 4, the H-Cl bond is covalent. At n = 5 and 7, the contact ion pair of H3O+Cl- is spontaneously formed within the hydrogen bond networks of book and cube acid-water clusters, respectively.

Imaging surface structure and premelting of ice Ih with atomic resolution

Ice surfaces are closely relevant to many physical and chemical properties, such as melting, freezing, friction, gas uptake and atmospheric reaction1-8. Despite extensive experimental and theoretical investigations9-17, the exact atomic structures of ice interfaces remain elusive owing to the vulnerable hydrogen-bonding network and the complicated premelting process. Here we realize atomic-resolution imaging of the basal (0001) surface structure of hexagonal water ice (ice Ih) by using qPlus-based cryogenic atomic force microscopy with a carbon monoxide-functionalized tip. We find that the crystalline ice-Ih surface consists of mixed Ih- and cubic (Ic)-stacking nanodomains, forming 19 × 19 periodic superstructures. Density functional theory reveals that this reconstructed surface is stabilized over the ideal ice surface mainly by minimizing the electrostatic repulsion between dangling OH bonds. Moreover, we observe that the ice surface gradually becomes disordered with increasing temperature (above 120 Kelvin), indicating the onset of the premelting process. The surface premelting occurs from the defective boundaries between the Ih and Ic domains and can be promoted by the formation of a planar local structure. These results put an end to the longstanding debate on ice surface structures and shed light on the molecular origin of ice premelting, which may lead to a paradigm shift in the understanding of ice physics and chemistry.

Structure and dynamics of dissociated and undissociated forms of nitric acid and their implications in interfacial mass transfer: insights from molecular dynamics simulations

Nitric acid (HNO3) is widely used in various chemical and nuclear industries. Therefore, it is important to develop an understanding of the different forms of nitric acid for its practical applications. Molecular dynamics (MD) simulation is one of the best tools to investigate the behavior of concentrated nitric acid in aqueous solution with various forms together with pure nitric acid to identify a suitable model of nitric acid for use in simulations of biphasic systems for interfacial mass transfer. The Mulliken partial charge embedded OPLS-AA force field was used to model the neutral nitric acid, hydronium ion and nitrate ion, and it was found that the Mulliken partial charge embedded force field works quite well. The computed density of the dissociated and mixed-form acid was in good agreement with the experimental values. In water, the HNO3 molecule was seen to be coordinated with three water molecules in the first sphere of coordination. The distribution of water surrounding the HNO3 molecule and nitrate ion was corroborated by the DFT-optimized hydrated cluster. The calculated diffusivity values of the neutral acid and ions were significantly higher in the mixed form of nitric acid, which is an important dynamic quantity controlling the kinetics of the liquid-liquid interfacial extraction. The structural analysis revealed that the local aggregation is minimized when both forms of acid are present together in the solution. The water-ion and water-neutral acid interactions were predicted to be enhanced, as confirmed by H-bond studies. The shear viscosity of the mixed acid exhibited excellent agreement with the experimental values, which again confirms the consideration of the mixed form of nitric acid. The simulated value of surface tension for the mixed form of acid also appeared to be quite accurate based on the surface tension of water. The mixed form of nitric acid comprising both forms of acid is the best representation for nitric acid to be considered for MD simulations of biphasic systems. The mixed form of nitric acid established that the concentrated nitric acid may not be present either in the fully dissociated form or fully undissociated form in the solution.

Acid Dissociation in (HX)n (H2 O)n Clusters (X=F, Cl, Br, I; n=2, 3)

In this work, we analyze the interactions between two or three hydrogen halide molecules and the same number of water moieties through a systematic exploration of their potential energy surfaces. Our results indicate that the most stable HF and HCl aggregates do not experience dissociation of any of the acid fragments, even with three water molecules. In contrast, in the HBr and HI clusters, one of the acid fragments does dissociate. While the global minimum of (HBr)₃ (H₂ O)₃ is a hydrogen-bridged bihalide anion (BrHBr^- ), which is persistent at temperatures up to 203 K, the lowest energy structure of (HI)₃ (H₂ O)₃ has a separated ion pair, but the motif with a bihalide anion (IHI^- ) is only 0.2 kcal mol^-1 above the global minimum. Among the more stable structures is a broad spectrum of contacts, including water⋯water, HX⋯water, and HX⋯HX hydrogen bonds, halogen bonds, ionic and long-range X⋯H contacts.

Evidence for Oxonium Ions in Ethereal "Hydrogen Chloride"

Although solutions of hydrogen chloride in ethereal solvents like diethyl ether or dioxane are commonly employed in the laboratory, the solution structure of these reagents has yet to be firmly established. Here, we analyze solutions of ethereal hydrogen chloride or deuterium chloride in toluene, in dichloromethane, or in the absence of a co-solvent by in situ infrared spectroscopy. The resulting spectra are inconsistent with free HCl or often-proposed 1:1 HCl-ether complexes but closely match the predicted spectra of oxonium ions generated via protonation of diethyl ether. Molecular dynamics simulation of the oxonium chloride complexes provides evidence for an outer-sphere contact ion pair. These results suggest new approaches for tuning the acidity of strong Brønsted acids in organic solvents and demonstrate the importance of conducting spectroscopic measurements under reaction-relevant conditions.

Effect of Hydration on Electron Attachment to Methanesulfonic Acid Clusters

We report an experimental and computational study of the electron-induced chemistry of methanesulfonic acid (MSA, MeSO₃H) in clusters. We combine the mass spectra after the 70 eV electron ionization with the negative ion spectra after electron attachment (EA) at low electron energies of 0-15 eV of the MSA molecule, small MSA clusters, and microhydrated MSA clusters to reveal the solvation effects. The MSA/He coexpansion only generates small MSA clusters with up to four molecules, but adding water substantially hydrates the MSA clusters, resulting in clusters composed of 1-2 MSA molecules accompanied by quite a few water molecules. The clustering strongly suppresses the fragmentation of the MSA molecules upon both the positive ionization and EA. The electron-energy-dependent ion yield for different negative ions is measured. For the MSA molecule and pure MSA clusters, EA leads to an H-abstraction yielding MeSO₃^-. It proceeds efficiently at low electron energies below 2 eV with a shoulder at 3-4 eV and a broad, almost 2 orders of magnitude weaker, peak around 8 eV. The hydrated (H₂O)_nMeSO₃^- ions with n ≤ 3 exhibit only a broad peak around 7 eV similar to EA of pure water clusters. Thus, for the small clusters, the electron attachment and hydrogen abstraction from water occur. On the other hand, the larger clusters with n > 4 display a peak below 2 eV, which quickly dominates the spectrum with increasing n. This peak is related to the formation of the H₃O⁺·MeSO₃^- ion pair upon hydration and subsequent dipole-supported electron attachment followed by the hydronium neutralization and H₃O^• radical dissociation. The size-resolved experimental data indicate that the ionic dissociation of MSA starts to occur in the neutral MeSO₃H(H₂O)_N clusters with about four water molecules.

Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective

As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.

Generating Excess Protons in Microsolvated Acid Clusters under Ambient Conditions: An Issue of Configurational Entropy versus Internal Energy

Acid dissociation, and thus liberation of excess protons in small water droplets, impacts on diverse fields such as interstellar, atmospheric or environmental chemistry. At cryogenic temperatures below 1 K, it is now well established that as few as four water molecules suffice to dissociate the generic strong acid HCl, yet temperature-driven recombination sets in simply upon heating that cluster. Here, the fundamental question is posed of how many more water molecules are required to stabilize a hydrated excess proton at room temperature. Ab initio path integral simulations disclose that not five, but six water molecules are needed at 300 K to allow for HCl dissociation independently from nuclear quantum effects. In order to provide the molecular underpinnings of these observations, the classical and quantum free energy profiles were decomposed along the dissociation coordinate in terms of the corresponding internal energy and entropy profiles. What decides in the end about acid dissociation, and thus ion pair formation, in a specific microsolvated water cluster at room temperature is found to be a fierce competition between classical configurational entropy and internal energy, where the former stabilizes the undissociated state whereas the latter favors dissociation. It is expected that these are generic findings with broad implications on acid-base chemistry depending on temperature in small water assemblies.

Interfacial Vibrational Dynamics of Ice Ih and Liquid Water

Insights into energy flow dynamics at ice surfaces are essential for understanding chemical dynamics relevant to atmospheric and geographical sciences. Here, employing ultrafast surface-specific spectroscopy, we report the interfacial vibrational dynamics of ice I_h. A comparison to liquid water surfaces reveals accelerated vibrational energy relaxation and dissipation at the ice surface for hydrogen-bonded OH groups. In contrast, free-OH groups sticking into the vapor phase exhibit substantially slower vibrational dynamics on ice. The acceleration and deceleration of vibrational dynamics of these different OH groups at the ice surface are attributed to enhanced intermolecular coupling and reduced rotational mobility, respectively. Our results highlight the unique properties of free-OH groups on ice, putatively linked to the high catalytic activities of ice surfaces.

Acid solvation versus dissociation at "stardust conditions": Reaction sequence matters

Chemical reactions at ultralow temperatures are of fundamental importance to primordial molecular evolution as it occurs on icy mantles of dust nanoparticles or on ultracold water clusters in dense interstellar clouds. As we show, studying reactions in a stepwise manner in ultracold helium nanodroplets by mass-selective infrared (IR) spectroscopy provides an avenue to mimic these "stardust conditions" in the laboratory. In our joint experimental/theoretical study, in which we successively add H₂O molecules to HCl, we disclose a unique IR fingerprint at 1337 cm^-1 that heralds hydronium (H₃O⁺) formation and, thus, acid dissociation generating solvated protons. In stark contrast, no reaction is observed when reversing the sequence by allowing HCl to interact with preformed small embryonic ice-like clusters. Our ab initio simulations demonstrate that not only reaction stoichiometry but also the reaction sequence needs to be explicitly considered to rationalize ultracold chemistry.

Size-Dependent Formation of an Ion Pair in HSO4-(H2O) n: A Molecular Model for Probing the Microsolvation of Acid Dissociation

With a p K_a of 2.0, HSO₄^- is not a strong acid as its dissociation percentage is only ∼10% in a solution of 1 M. However, our ab initio molecular dynamics and density functional theory calculations show that acid dissociation in the hydrate clusters, HSO₄^-(H₂O) _n, is quite facile at moderate sizes. It starts at around n = 12 and is completed by n = 16 when the energetics becomes very favorable. The dissociation explains the significant broadening at n = 16 in the S═O stretching region of the previously reported infrared photodissociation spectra for HSO₄^-(H₂O) _n as the solvation shell is tightened around the sulfate dianion and the proton. More importantly, HSO₄^-(H₂O) _n should provide an ideal model to probe the molecular details involved in an acid dissociation process.

Far-Infrared Synchrotron Spectroscopy and Quantum Chemical Calculations of the Potentially Important Interstellar Molecule, 2-Chloroethanol

The high brightness of the Australian synchrotron allowed for detailed spectra to be collected at high resolution (0.00096 cm^-1) in the vicinity of the a/ b/ c-type ν₁₉ band of 2-chloroethanol, which involves O-H torsional motion about the C-O bond. A rovibrational analysis was performed for both chlorine isotopologues in the ν₁₉ fundamental (centered at ∼344 cm^-1) which involved the assignment of 7153 lines ( J ≤ 90, K _a ≤ 41). A global fit to these lines in addition to 119 microwave lines ( J ≤ 29, K _a ≤ 11) led to the determination of spectroscopic constants up to the sextic level in both the ground and excited states using Watson's A-reduction Hamiltonian. The constants agree well with those calculated at the anharmonic MP2/cc-pVTZ level and allow for spectroscopically accurate predictions of rotational transitions in the ground vibrational state to be made over a broad range of rotational energies ( T_R < 1000 K). We explored the role that 2-chloroethanol might play in interstellar molecular clouds by performing calculations on the substitution reaction between HCl and ethylene glycol, and the addition reaction between HCl and oxirane, all of which have been observed in Sagittarius B2(N) and are expected to play important roles in the chemistry that occurs on the icy mantles of interstellar dust grains. While both reactions have relatively high activation barriers, the HCl + oxirane reaction was found be much more exothermic; further calculations on it indicate that a water-like environment significantly reduces the barrier while slightly increasing its exothermicity. These results suggest that 2-chloroethanol could be efficiently produced from the cosmic ray bombardment of common interstellar ices.

A systematic ab initio optimization of monohydrates of HCl•HNO3•H2SO4 aggregates

Hydrates of HCl, HNO₃ and H₂SO₄ involved in polar stratospheric clouds capture the attention of researchers due to the mixtures composed with them. The molecular aggregates generated with these strong acids show different behaviors, geometries and nucleation reactions at atmospheric temperatures. Here is presented a systematic ab initio optimization study of monohydrates of HCl•HNO₃•H₂SO₄ using the Density Functional Theory, by means of geometry optimizations carried out with B3LYP hybrid method and aug-cc-pVTZ basis set, a high level of theory, within Gaussian 09 program. This systematic optimization procedure consists to situate systematically the H₂O molecule around the cluster in study, on the favorable positions to develop higher quantity of hydrogen bonds as possible, in order to obtain major quantity of different electronic structures of these monohydrates. Applying this systematic optimization methodology over previously optimized complexes of HCl, HNO₃ and H₂SO₄, the present theoretical approach provides thirty-two different optimized electronic structures of monohydrates that were yielded from seven initial groups of (HCl•HNO₃•H₂SO₄)-complex, placing the H₂O in eight positions around them. Moreover, their Infrared spectra have been predicted for all (HCl•HNO₃•H₂SO₄)-monohydrates achieved. Likewise, It is shown the outcomes of the electronic energies, relative Gibbs free energies, Infrared spectra, the wavenumbers of hydrogen bonds, inter-monomeric parameters, electronic structures of (HCl•HNO₃•H₂SO₄)-monohydrates. These monohydrates could be considered precursors of the atmospheric heterogeneous nucleation reactions. These results can be useful to experimentalists of Catalysis, Astrophysics, Corrosion of metals and ceramics, aromatic compounds reactions, even environmental pollution and industrial smog.

Water Structure Controls Carbonic Acid Formation in Adsorbed Water Films

Reaction pathways and kinetics in highly structured H₂O adsorbed as Ångstrom to nanometer thick layers on mineral surfaces are distinct from those facilitated by bulk liquid water. We investigate the role of the interfacial H₂O structure in the reaction of H₂O and CO₂ to form carbonic acid (H₂CO₃) in thin H₂O films condensed onto silica nanoparticles from humidified supercritical CO₂. Rates of carbonic acid formation are correlated with spectroscopic signatures of H₂O structure using oxygen isotopic tracers and infrared spectroscopy. While carbonic acid virtually does not form in the supercritical phase, the silica surface catalyzes this reaction by concentrating H₂O through adsorption at hydrophilic silanol groups. Within measurement uncertainty, we found no evidence that carbonic acid forms when exclusively ice-like structured H₂O is detected at the silica surface. Instead, formation of H₂C¹⁸O¹⁶O₂ from H₂¹⁸O and C¹⁶O₂ was found to be linearly correlated with liquid-like structured H₂O that formed on the ice-like layer.

Mass spectrometry of aerosol particle analogues in molecular beam experiments

Nanometer-size particles such as ultrafine aerosol particles, ice nanoparticles, water nanodroplets, etc, play an important, however, not yet fully understood role in the atmospheric chemistry and physics. These species are often composed of water with admixture of other atmospherically relevant molecules. To mimic and investigate such particles in laboratory experiments, mixed water clusters with atmospherically relevant molecules can be generated in molecular beams and studied by various mass spectrometric methods. The present review demonstrates that such experiments can provide unprecedented details of reaction mechanisms, and detailed insight into the photon-, electron-, and ion-induced processes relevant to the atmospheric chemistry. After a brief outline of the molecular beam preparation, cluster properties, and ionization methods, we focus on the mixed clusters with various atmospheric molecules, such as hydrated sulfuric acid and nitric acid clusters, N_x O_y and halogen-containing molecules with water. A special attention is paid to their reactivity and solvent effects of water molecules on the observed processes.

Hydrogen-bonding behavior of various conformations of the HNO3…(CH3OH)2 ternary system

Nine minima were found on the intermolecular potential energy surface for the ternary system HNO₃(CH₃OH)₂ at the MP2/aug-cc-pVDZ level of theory. The cooperative effect, which is a measure of the hydrogen-bonding strength, was probed in these nine conformations of HNO₃…(CH₃OH)₂. The results are discussed here in terms of structures, energetics, infrared vibrational frequencies, and topological parameters. The cooperative effect was observed to be an important contributor to the total interaction energies of the cyclic conformers of HNO₃…(CH₃OH)₂, meaning that it cannot be neglected in simulations in which the pair-additive potential is applied. Graphical abstract The H-bonding behavior of various conformations of the HNO₃(CH₃OH)₂ trimer was investigated.

Acid Dissociation in HCl-Water Clusters is Temperature Dependent and Cannot be Detected Based on Dipole Moments

The dissociation of acids in aqueous environments at low temperatures in the presence of a limited amount of water is underlying a wealth of processes from atmospheric to interstellar science. For the paradigmatic case of HCl(H_{2}O)_{n} clusters, our extensive ab initio path integral simulations quantify in terms of free energy differences and barriers that n=4 water molecules are indeed required to dissociate HCl at low temperatures. Increasing the temperature, however, reverses the process and thus counteracts dissociation by fluctuation-driven recombination. The size of the electric dipole moment is shown to not correlate with the acid being in its dissociated or molecular state, thus rendering its measurement as a function of n unable to detect the dissociation transition.

Coexistence of Physisorbed and Solvated HCl at Warm Ice Surfaces

The interfacial ionization of strong acids is an essential factor of multiphase and heterogeneous chemistry in environmental science, cryospheric science, catalysis research and material science. Using near ambient pressure core level X-ray photoelectron spectroscopy, we directly detected a low surface coverage of adsorbed HCl at 253 K in both molecular and dissociated states. Depth profiles derived from XPS data indicate the results as physisorbed molecular HCl at the outermost ice surface and dissociation occurring upon solvation deeper in the interfacial region. Complementary X-ray absorption measurements confirm that the presence of Cl^- ions induces significant changes to the hydrogen bonding network in the interfacial region. This study gives clear evidence for nonuniformity across the air-ice interface and questions the use of acid-base concepts in interfacial processes.

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

Atmospheric chemical reactions of alternatives of polybrominated diphenyl ethers initiated by OH: A case study on triphenyl phosphate

Many studies have been performed to evaluate the environmental risk caused by alternative flame retardants (AFRs) of polybrominated diphenyl ethers due to their ubiquitous occurrence in the environment. However, as an indispensable component of the environmental risk assessment, the information on atmospheric fate of AFRs is limited although some AFRs have been frequently and highly detected in the atmosphere. Here, a combined quantum chemical method and kinetics modeling were used to investigate atmospheric transformation mechanism and kinetics of AFRs initiated by OH in the presence of O2, taking triphenyl phosphate (TPhP) as a case. Results show that the pathway involving initial OH addition to phenyl of TPhP to form TPhP-OH adduct, and subsequent reaction of the TPhP-OH adduct with O2 to finally form phenol phosphate, is the most favorable for the titled reaction. The calculated overall reaction rate constant is 1.6×10(-12)cm(3) molecule(-1)s(-1), translating 7.6days atmospheric lifetime of TPhP. This clarifies that gaseous TPhP has atmospheric persistence. In addition, it was found that ice surface, as a case of ubiquitous water in the atmosphere, has little effect on the kinetics of the rate-determining step in the OH-initiated TPhP reaction.

Adsorption isotherms for hydrogen chloride (HCl) on ice surfaces between 190 and 220 K

The interaction of hydrogen chloride (HCl) with ice surfaces at temperatures between 190 and 220 K was investigated using a coated-wall flow-tube connected to a chemical ionization mass spectrometer. Equilibrium surface coverages of HCl were determined at gas phase concentrations as low as 2 × 10(9) molecules cm(-3) (∼4 × 10(-8) Torr at 200 K) to derive Langmuir adsorption isotherms. The data are described by a temperature independent partition coefficient: KLang = (3.7 ± 0.2) × 10(-11) cm(3) molecule(-1) with a saturation surface coverage Nmax = (2.0 ± 0.2) × 10(14) molecules cm(-2). The lack of a systematic dependence of KLang on temperature contrasts the behaviour of numerous trace gases which adsorb onto ice via hydrogen bonding and is most likely related to the ionization of HCl at the surface. The results are compared to previous laboratory studies, and the equilibrium partitioning of HCl to ice surfaces under conditions relevant to the atmosphere is evaluated.

How Many Water Molecules Does it Take to Dissociate HCl?

The potential energy surfaces of the HCl(H2O)n (n is the number of water molecules) clusters are systematically explored using density functional theory and high-level ab initio computations. On the basis of electronic energies, the number of water molecules needed for HCl dissociation is four as reported by some experimental groups. However, this number is five owing to the inclusion of entropic factors. Wiberg bond indices are calculated and analyzed, and the results provide a quadratic correlation and classification of clusters according to the nondissociated, partially dissociated, and fully dissociated character of the H-Cl bond. Our computations show that if temperature is not controlled during the experiment, the values obtained for the dipole moment (or for any measurable property) are susceptible to change, providing a different picture of the number of water molecules needed for HCl dissociation in a nanoscopic droplet.

Active Site Dynamical Effects in the Hydrogen Transfer Rate-limiting Step in the Catalysis of Linoleic Acid by Soybean Lipoxygenase-1 (SLO-1): Primary and Secondary Isotope Contributions

Using ab initio molecular dynamics (AIMD) simulations that facilitate the treatment of rare events, we probe the active site participation in the rate-determining hydrogen transfer step in the catalytic oxidation of linoleic acid by soybean lipoxygenase-1 (SLO-1). The role of two different active site components is probed. (a) On the hydrogen atom acceptor side of the active site, the hydrogen bonding propensity between the acceptor side hydroxyl group, which is bound to the iron cofactor, and the backbone carboxyl group of isoleucine (residue number 839) is studied toward its role in promoting the hydrogen transfer event. Primary and secondary (H/D) isotope effects are also probed and a definite correlation with subtle secondary H/D isotope effects is found. With increasing average nuclear kinetic energy, the increase in transfer probability is enhanced due to the presence of the hydrogen bond between the backbone carbonyl of I839 and the acceptor oxygen. Further increase in average nuclear kinetic energy reduces the strength of this secondary hydrogen bond which leads to a deterioration in hydrogen transfer rates and finally embrances an Arrhenius-like behavior. (b) On the hydrogen atom donor side, the coupling between vibrational modes predominantly localized on the donor-side linoleic acid group and the reactive mode is probed. There appears to be a qualitative difference in the coupling between modes that belong to linoleic acid and the hydrogen transfer mode, for hydrogen and deuterium transfer. For example, the donor side secondary hydrogen atom is much more labile (by nearly a factor of 5) during deuterium transfer as compared to the case for hydrogen transfer. This appears to indicate a greater coupling between the modes belonging to the linoleic acid scaffold and the deuterium transfer mode and also provides a new rationalization for the abnormal (nonclassical) secondary isotope effect results obtained by Knapp, Rickert, and Klinman in J. Am. Chem. Soc. , 2002 , 124 , 3865 . To substantiate our findings noted in point a above, we have suggested an I839 → A839 or I839 → V839 mutation. This will modify the bulkiness of hydrogen the bonding residue, allowing greater flexibility in the secondary hydrogen bond formation highlighted above and adversely affecting the reaction rate.

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.

Microhydration of contact ion pairs in M(2+)OH(-)(H2O)(n=1-5) (M = Mg, Ca) clusters: spectral manifestations of a mobile proton defect in the first hydration shell

Vibrational predissociation spectra of D2-"tagged" Mg(2+)OH(-)(H2O)n=1-6 and Ca(2+)OH(-)(H2O)n=1-5 clusters are reported to explore how the M(2+)OH(-) contact ion pairs respond to stepwise formation of the first hydration shell. In both cases, the hydroxide stretching frequency is found to red-shift strongly starting with addition of the third water molecule, quickly becoming indistinguishable from nonbonded OH groups associated with solvent water molecules by n = 5. A remarkably broad feature centered around 3200 cm(-1) and spanning up to ∼1000 cm(-1) appears for the n ≥ 4 clusters that we assign to a single-donor ionic hydrogen bond between a proximal first solvent shell water molecule and the embedded hydroxide ion. The extreme broadening is rationalized with a theoretical model that evaluates the range of local OH stretching frequencies predicted for the heavy particle configurations available in the zero-point vibrational wave function describing the low-frequency modes. The implication of this treatment is that extreme broadening in the vibrational spectrum need not arise from thermal fluctuations in the ion ensemble, but can rather reflect combination bands based on the OH stretching fundamental that involve many quanta of low-frequency modes whose displacements strongly modulate the OH stretching frequency.

Quantum chemical study of atmospheric aggregates: HCl•HNO3•H2SO4

HCl, HNO3 and H2SO4 are implicated in atmospheric processes in areas such as polar stratospheric clouds in the stratosphere. Ternary complexes of HCl, HNO3 and H2SO4 were investigated by ab initio calculations at B3LYP level of theory with aug-cc-pVTZ and aug-cc-pVQZ basis sets, taking into account basis set superposition error (BSSE). The results were assessed in terms of structures (five hexagonal cyclic structures and two quasi-pentagonal cyclic structures), inter-monomeric parameters (all ternary complexes built three hydrogen bonds), energetics (seven minima obtained), infrared harmonic vibrational frequencies (red shifting of complexes from monomers), and relative stability of complexes, which were favorable when the temperature decreases under stratospheric conditions, from 298 K to 188 K, and in concrete, at 210 K, 195 K and 188 K.

Dynamical aspects of intermolecular proton transfer in liquid water and low-density amorphous ices

The microscopic dynamics of an excess proton in water and in low-density amorphous ices has been studied by means of a series of molecular dynamics simulations. Interaction of water with the proton species was modelled using a multistate empirical valence bond Hamiltonian model. The analysis of the effects of low temperatures on proton diffusion and transfer rates has been considered for a temperature range between 100 and 298 K at the constant density of 1 g cm(-3). We observed a marked slowdown of proton transfer rates at low temperatures, but some episodes are still seen at 100 K. In a similar fashion, mobility of the lone proton gets significantly reduced when temperature decreases below 273 K. The proton transfer in low-density amorphous ice is an activated process with energy barriers between 1-10 kJ/mol depending of the temperature range considered and eventually showing Arrhenius-like behavior. Spectroscopic data indicated the survival of both Zundel and Eigen structures along the whole temperature range, revealed by significant spectral frequency shifts.

Photoreaction of rac-leucine in ice by circularly polarized synchrotron radiation: temperature-induced mechanism switching from Norrish Type II to deamination

The delivery of extraterrestrial organics to primitive Earth is considered to have triggered the origin and subsequent evolution of life. Indeed, enantiomerically enriched amino acids of nonterrestrial origin have been found in carbonaceous meteorites, and enantioselective photodecomposition by circularly polarized light (CPL) in outer space has been proposed to have played some role in the initial enantiomeric bias. To experimentally examine this possibility and elucidate the photoreaction mechanisms, we have studied the photolysis of racemic leucine (rac-Leu) in acidic and neutral ice/water media at 21-298 K with left- and right-CPL in an attempt to detect enantiomerically enriched D- and L-Leu, respectively. Comprehensive product analyses revealed that the CPL-induced deracemization of Leu proceeds in both acidic and neutral ice matrices even at 21 K, and that the main mechanism switches from Norrish-type II γ-hydrogen abstraction to SN i deamination on lowering the temperature. The potential role of the CPL-induced photodecomposition of amino acids as a source of the enantiomer imbalance in meteorites is discussed.

Ionization of Nitric Acid on Crystalline Ice: The Role of Defects and Collective Proton Movement

Ionization of nitric acid (HNO3) on a model ice surface is studied using ab initio molecular dynamics at temperatures of 200 and 40 K with a surface slab model that consists of the ideal ice basal plane with locally optimized and annealed defects. Pico- and subpicosecond ionization of nitric acid can be achieved in the defect sites. Key features of the rapid ionization are (a) the efficient solvation of the polyatomic nitrate anion, by stealing hydrogen bonds from the weakened hydrogen bonds at defect sites, (b) formation of contact ion pairs to stable "presolvated" molecular species that are present at the defects,

Probing Small Protonated Water Clusters in Acetonitrile Solutions by 1H NMR

In a previous publication by Kalish et al. (J. Phys. Chem. A 115 (2011) 4063) the existence of well defined small protonated water clusters in acetonitrile has been established by IR spectroscopy. Here we report on a 1H NMR study of triflic acid, CF3SO3H, in acetonitrile-water solutions. Using NMR we are able to corroborate the general solvation scheme we have proposed for the hydrated proton in acetonitrile as a function of the molar ratio between the strong mineral acid and water, n = [H2O]/[acid]. According to this scheme, backed now by both IR absorption spectroscopy and NMR measurements, the very strong triflic acid completely dissociates in acetonitrile/water solutions to yield the aqueous proton and the triflate anion when n > 1. Furthermore, increasing n results in the proton solvated in increasingly larger water clusters formed within the acetonitrile solution.

Clearly distinguishable by NMR are the smallest protonated water clusters, the protonated water monomer, H3 +O, and the protonated water dimer, H5 +O2, which dominate the solution for n = 1,2,3. For larger n the NMR study indicates the gradual increase of the average protonated water cluster size as a function of n while the proton inner solvation core more closely retaining the characteristics of a deformed protonated water dimer, (H2O-H+⋯OH2) s than that of the protonated water monomer (H3 +O) s .

Probing molecular solids with low-energy ions

Ion/surface collisions in the ultralow- to low-energy (1-100-eV) window represent an excellent technique for investigation of the properties of condensed molecular solids at low temperatures. For example, this technique has revealed the unique physical and chemical processes that occur on the surface of ice, versus the liquid and vapor phases of water. Such instrument-dependent research, which is usually performed with spectroscopy and mass spectrometry, has led to new directions in studies of molecular materials. In this review, we discuss some interesting results and highlight recent developments in the area. We hope that access to the study of molecular solids with extreme surface specificity, as described here, will encourage investigators to explore new areas of research, some of which are outlined in this review.

Gas-phase calorimetry of protonated water clusters

Protonated water clusters with 60 to 79 molecules have been studied by nanocalorimetry. The technique is based on multi-collision excitations of the accelerated clusters with helium. The caloric curves indicate transitions that resemble those of water clusters charged by an excess electron, but the transition temperatures of the protonated clusters are higher.

Structure, energetics, and infrared spectra of weakly bound HC(2n+1)N···HCl complexes. A theoretical study

Three model systems, HCN···HCl, HC(3)N···HCl, and HC(5)N···HCl, have been investigated computationally with the use of the second-order Möller-Plesset (MP2) and the coupled cluster (with single and double excitations and noniterative inclusion of triples) methods. The global minima are linear hydrogen-bonded structures with HCl as a proton donor. Bent structures are proton-side complexes with HCl as an electron donor, while the bifurcated hydrogen bonds are predicted for t-shape complexes. One of the most important findings in this paper is that, according to symmetry-adapted perturbation analysis, the induction-to-dispersion ratios are the biggest for linear complexes, and it is the most noticeable difference between linear, bent, and t-shape structures.

Comment on "HCl adsorption on ice at low temperature: a combined X-ray absorption, photoemission and infrared study" by P. Parent, J. Lasne, G. Marcotte and C. Laffon, Phys. Chem. Chem. Phys., 2011, 13, 7142

On the basis of NEXAFS, photoemission and FTIR spectra of ice films with low doses of adsorbed HCl, the authors of the PCCP paper "HCl adsorption on ice at low temperature: a combined X-ray absorption, photoemission and infrared study", Phys. Chem. Chem. Phys., 2011, 13, 7142, have come to conclusions regarding the behavior of submonolayer amounts of HCl at 50 K that contradict published results of the authors of this Comment. Our purpose is to argue that the conclusion, attributed going forward to PLML (authors' initials), that nearly 100% of HCl ionizes for dosage levels near to 0.16 monolayer (ML) or 0.3 Langmuir (L) at 50 K is questionable. Rather, we reaffirm our conclusions of much lower levels of ionization for similar temperatures and HCl dosages based on reactive ion scattering (RIS) and low energy sputtering (LES) data for ice films and FTIR spectra of ice nanocrystals. A second current paper by Ayotte et al., J. Phys. Chem. A, 2011, 115, 6002, that largely parallels in method and results the RAIR spectroscopy of PLML, is also given special notice.