Chemistry. Building solutions--one molecule at a time

Molecular vibrational spectral simulation connects theoretical cluster structure identification and vibrational spectral evidence

Structure identification of molecular clusters has long been a fundamental and challenging issue for cluster science. The traditional theoretical optimization on the potential energy surface heavily depends on the levels of theory and sometimes diverse identifications were reported. A solution to these disputations is to reinspect the theoretical results with the experimental data such as vibrational predissociation spectra with high sensitivity to the molecular cluster structures. Herein, the combination of global low-lying structure search and vibrational predissociation spectral simulation is proposed as an accurate and reliable approach for cluster structure identification, by which the assignments can be validated using experimental measurements. The qualitative agreement between simulated and measured vibrational spectra lends solid experimental evidence to the assignment of the cluster structures. Taking NH₄⁺(H₂O)_n (n = 2-4) as an example, we have unambiguously identified their structures and directly demonstrated the coexistence of two NH₄⁺(H₂O)₄ isomers (with 3 and 4 water molecules directly linked to NH₄⁺, respectively), which were debatable in previous studies. The developed methods would pave the way to the structure determination of the molecular clusters.

Grotthuss Molecular Dynamics Simulations for Modeling Proton Hopping in Electrosprayed Water Droplets

Excess protons in water exhibit unique transport properties because they can rapidly hop along H-bonded water wires. Considerable progress has been made in unraveling this Grotthuss diffusion mechanism using quantum mechanical-based computational techniques. Unfortunately, high computational cost tends to restrict those techniques to small systems and short times. Molecular dynamics (MD) simulations can be applied to much larger systems and longer time windows. However, standard MD methods do not permit the dissociation/formation of covalent bonds, such that Grotthuss diffusion cannot be captured. Here, we bridge this gap by combining atomistic MD simulations (using Gromacs and TIP4P/2005 water) with proton hopping. Excess protons are modeled as hydronium ions that undergo H₃O⁺ + H₂O → H₂O + H₃O⁺ transitions. In accordance with ab initio MD data, these Grotthuss hopping events are executed in "bursts" with quasi-instantaneous hopping across one or more waters. The bursts are separated by regular MD periods during which H₃O⁺ ions undergo Brownian diffusion. The resulting proton diffusion coefficient agrees with the literature value. We apply this Grotthuss MD technique to highly charged water droplets that are in a size regime encountered during electrospray ionization (5 nm radius, ∼17,000 H₂O). The droplets undergo rapid solvent evaporation and occasional H₃O⁺ ejection, keeping them at ca. 81% of the Rayleigh limit. The simulated behavior is consistent with phase Doppler anemometry data. The Grotthuss MD technique developed here should be useful for modeling the behavior of various proton-containing systems that are too large for high-level computational approaches. In particular, we envision future applications related to electrospray processes, where earlier simulations used metal cations while in reality excess protons dominate.

Accessing different binding sites of a multifunctional molecule: IR spectroscopy of propargyl alcohol⋯water complexes in helium droplets

Helium droplets facilitate the formation of two distinct local minimum structures for 1 : 1 propargyl alcohol⋯water complex.

A THz/FTIR fingerprint of the solvated proton: evidence for Eigen structure and Zundel dynamics

Zundel (orange), Eigen (red) and hydration water (light blue) contributions to the THz/FIR extinction of the solvated proton.

Rovibrational energies and spectroscopic constants for H2O-Ng complexes

In this work, rovibrational energies and spectroscopic constants for the water -Ng complexes (Ng = He, Ne, Ar, Kr and Xe) were calculated through two different approaches (by solving the Nuclear Schrödinger equation and by applying the Dunham's method) and using two different potential energy curves (PEC). These PEC were determined using potential parameters obtained through molecular beam scattering experiments and accurate theoretical calculation, respectively. It was found that the theoretical rovibrational energies are in a good agreement (only for the lowest numbers of vibrational states) with those obtained through experimental PEC. Another important conclusions was regarding the calculated first two rovibrational energies for the H 2 O-Ar system, that are in a good agreement with the experimental data.

Accurate predictions of water cluster formation, (H₂O)(n=2-10)

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2-10 water molecules. This approach reveals new low energy conformers for (H(2)O)(n=7,9,10). The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Water-chloride and water-bromide hydrogen-bonded networks: influence of the nature of the halide ions on the stability of the supramolecular assemblies

Two compounds, namely, [TTPH(2)](Cl)(2) x 4 H(2)O (1) and [TTPH(2)](Br)(2) x 4 H(2)O (2), (TTP = 4'-p-tolyl-2,2':6',2''-terpyridine) were synthesized from purely aqueous media and characterized by physical techniques. In the solid-state structures of these compounds, interesting supramolecular assemblies are observed. In 1, an unusual staircase-like architecture of the tape of edge-shared planar water hexamer is of importance, where the chloride ions are at the two edges of the tape. In 2, the polymeric nature of the water-bromide assembly is of interest, where discrete open-cube water octamers are doubly bridged by bromide ions. Semiempirical and DFT calculations confirm that the nature of the anion indeed affects the topology of the water-halide assemblies. We conclude that the protonated [TTPH(2)](2+) species can act as appropriate receptors for halide ions, which in turn act as a matrix for the formation of polymeric 1D water-halide assemblies.

Hydration of the Calcium Dication: Direct Evidence for Second Shell Formation from Infrared Spectroscopy

Infrared laser action spectroscopy in a Fourier-transform ion cyclotron resonance mass spectrometer is used in conjunction with ab initio calculations to investigate doubly charged, hydrated clusters of calcium formed by electrospray ionization. Six water molecules coordinate directly to the calcium dication, whereas the seventh water molecule is incorporated into a second solvation shell. Spectral features indicate the presence of multiple structures of Ca(H2O)(7)2+ in which outer-shell water molecules accept either one (single acceptor) or two (double acceptor) hydrogen bonds from inner-shell water molecules. Double-acceptor water molecules are predominantly observed in the second solvent shells of clusters containing eight or nine water molecules. Increased hydration results in spectroscopic signatures consistent with additional second-shell water molecules, particularly the appearance of inner-shell water molecules that donate two hydrogen bonds (double donor) to the second solvent shell. This is the first reported use of infrared spectroscopy to investigate shell structure of a hydrated multiply charged cation in the gas phase and illustrates the effectiveness of this method to probe the structures of hydrated ions.

Proton transfer 200 years after von Grotthuss: insights from ab initio simulations

In the last decade, ab initio simulations and especially Car-Parrinello molecular dynamics have significantly contributed to the improvement of our understanding of both the physical and chemical properties of water, ice, and hydrogen-bonded systems in general. At the heart of this family of in silico techniques lies the crucial idea of computing the many-body interactions by solving the electronic structure problem "on the fly" as the simulation proceeds, which circumvents the need for pre-parameterized potential models. In particular, the field of proton transfer in hydrogen-bonded networks greatly benefits from these technical advances. Here, several systems of seemingly quite different nature and of increasing complexity, such as Grotthuss diffusion in water, excited-state proton-transfer in solution, phase transitions in ice, and protonated water networks in the membrane protein bacteriorhodopsin, are discussed in the realms of a unifying viewpoint.

In search of CS2(H2O)n=1–4 clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS(2)(H2O)n, where n=1-4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS(2)(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS(2)(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS(2) in the atmosphere, the concentrations of these clusters are predicted to be on the order of 10(5) CS(2)(H2O) clusters.cm(-3) and 10(2) CS(2)(H2O)(2) clusters.cm(-3) at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS(2)(H(2)O)(n) structures predicted to be present in a low-temperature molecular beam experiment.

Ab Initio Molecular Dynamics Studies of Ionic Dissolution and Precipitation of Sodium Chloride and Silver Chloride in Water Clusters, NaCl(H2O)n and AgCl(H2O)n,n = 6, 10, and 14

An ab initio molecular dynamics method was used to compare the ionic dissolution of soluble sodium chloride (NaCl) in water clusters with the highly insoluble silver chloride (AgCl). The investigations focused on the solvation structures, dynamics, and energetics of the contact ion pair (CIP) and of the solvent-separated ion pair (SSIP) in NaCl(H(2)O)(n) and AgCl(H(2)O)(n) with cluster sizes of n = 6, 10 and 14. We found that the minimum cluster size required to stabilize the SSIP configuration in NaCl(H(2)O)(n) is temperature-dependent. For n = 6, both configurations are present as two distinct local minima on the free-energy profile at 100 K, whereas SSIP is unstable at 300 K. Both configurations, separated by a low barrier (<10 kJ mol(-1)), are identifiable on the free energy profiles of NaCl(H(2)O)(n) for n = 10 and 14 at 300 K, with the Na(+)/Cl(-) pairs being internally solvated in the water cluster and the SSIP configuration being slightly higher in energy (<5 kJ mol(-1)). In agreement with the low bulk solubility of AgCl, no SSIP minimum is observed on the free-energy profiles of finite AgCl(H(2)O)(n) clusters. The AgCl interaction is more covalent in nature, and is less affected by the water solvent. Unlike NaCl, AgCl is mainly solvated on the surface in finite water clusters, and ionic dissolution requires a significant reorganization of the solvent structure.

Molecular-beam study of the water-helium system: Features of the isotropic component of the intermolecular interaction and a critical test for the available potential-energy surfaces

We report molecular-beam measurements of the total integral cross sections for the scattering of water molecules by helium atoms. A combined analysis of the new experimental data together with available differential cross section results has allowed an accurate determination of the isotropic component of the interaction potential for this prototypical system. The potential well shows a depth of 0.265 +/- 0.010 kJ/mol at a distance between He and the center of mass of the water molecule of 0.345 +/- 0.02 nm. An effective isotropic long-range attraction constant C(LR) = (6.3+/-0.3) x 10(-4) kJ mol(-1) nm(-6), including both dispersion and induction contributions, has also been determined. The most recent and accurate ab initio potential-energy surfaces have been tested against these new experimental results.

Hydrogen-Bonded Networks in Ethanol Proton Wires: IR Spectra of (EtOH)qH+−Ln Clusters (L = Ar/N2, q ≤ 4, n ≤ 5)

Isolated and microsolvated protonated ethanol clusters, (EtOH)qH+-Ln with L = Ar and N2, are characterized by infrared photodissociation (IRPD) spectroscopy in the 3 microm range and quantum chemical calculations. For comparison, also the spectrum of the protonated methanol dimer, (MeOH)2H+, is presented. The IRPD spectra carry the signature of H-bonded (EtOH)qH+ chain structures, in which the excess proton is either strongly localized on one or (nearly) equally shared between two EtOH molecules, corresponding to Eigen-type ion cores (EtOH2+ for q = 1, 3) or Zundel-type ion cores (EtOH-H+-HOEt for q = 2, 4), respectively. In contrast to neutral (EtOH)q clusters, no cyclic (EtOH)qH+ isomers are detected in the size range investigated (q < or = 4), indicative of the substantial impact of the excess proton on the properties of the H-bonded ethanol network. The acidity of the two terminal OH groups in the (EtOH)qH+ chains decreases with the length of the chain (q). Comparison between (ROH)qH+ with R = CH3 and C2H5 shows that the acidity of the terminal O-H groups increases with the length of the aliphatic rest (R). The most stable (EtOH)qH+-Ln clusters with n < or = 2 feature intermolecular H-bonds between the inert ligands and the two available terminal OH groups of the (EtOH)qH+ chain. Asymmetric microsolvation of (EtOH)qH+ with q = 2 and 4 promotes a switch from Zundel-type to Eigen-type cores, demonstrating that the fundamental structural motif of the (EtOH)qH+ proton wire sensitively depends on the environment. The strength of the H-bonds between L and (EtOH)qH+ is shown to provide a rather sensitive probe of the acidity of the terminal OH groups.

Structure of protonated water clusters: Low-energy structures and finite temperature behavior

The structure of protonated water clusters H+(H2O)n (n=5-22) are examined by two Monte Carlo methods in conjunction with the OSS2 potential [L. Ojamae, I. Shavitt, and S. J. Singer J. Chem. Phys. 109, 5547 (1998)]. The basin-hopping method is employed to explore the OSS2 potential energy surface and to locate low-energy structures. The topology of the "global minimum," the most stable low-energy structure, changes from single ring to multiple ring to polyhedral cage as the cluster size grows. The temperature dependence of the cluster geometry is examined by carrying out parallel tempering Monte Carlo simulations. Over the temperature range we studied (25-330 K), all water clusters undergo significant structural changes. The trends are treelike structures dominating at high temperature and single-ring structures appearing in slightly lower temperatures. For n> or =7, an additional transition from single ring to multiple rings appears as the temperature decreases. Only for n> or =16 do polyhedral structures dominate the lowest temperature range. Our results indicate very dynamic structural changes at temperature range relevant to atmospheric chemistry and current experiments. The structures and properties of medium-sized protonated clusters in this temperature range are far from their global minimum cousins. The relevance of these findings to recent experiments and theoretical simulations is also discussed.

Electroosmotic enhancement of the binding of a neutral molecule to a transmembrane pore

The flux of solvent water coupled to the transit of ions through protein pores is considerable. The effect of this electroosmotic solvent flow on the binding of a neutral molecule [beta-cyclodextrin (betaCD)] to sites within the staphylococcal alpha-hemolysin pore was investigated. Mutant alpha-hemolysin pores were used to which betaCD can bind from either entrance and through which the direction of water flow can be controlled by choosing the charge selectivity of the pore and the polarity of the applied potential. The Kd values for betaCD for individual mutant pores varied by >100-fold with the applied potential over a range of -120 to +120 mV. In all cases, the signs of the changes in binding free energy and the influence of potential on the association and dissociation rate constants for betaCD were consistent with an electroosmotic effect.

Thermodynamic treatment of oligonucleotide duplex-simplex equilibria

Thermodynamic formulations have been devised to obtain Delta G degrees values directly from spectroscopic data at a fixed common temperature in nucleic acid duplex-simplex melting curves. In addition, the dependence of melting on salt concentration has been expressed in terms of a stepwise stoichiometric representation, which leads to a specific equation for the partition of the added sodium ions between the different oligonucleotide forms.

Atmospheric detection of water dimers via near-infrared absorption

Weakly bonded pairs of water molecules (H2O)2, or water dimers (WDs), may play an important role in photochemistry and climate, but the overlap of most of its spectral features with the water monomer (WM) has made detection difficult. We report on WD absorption measurements by means of atmospheric long-path (18.34 kilometers) differential optical absorption spectroscopy of the near-infrared OH stretching mode 0>f 4>b overtone transition predicted to be located near 746 nanometers. Our observation is in reasonable agreement with the known thermochemistry, calculated and measured structure, and spectroscopy (band strength, shape, and width) of the WD. The observation implies that the WD 0>f 4>b band is located at 749.5 nanometers, with a full width at half maximum of approximately 19.4 wave numbers, and that its band strength ranges between 1.23 x 10(-22) and 5.25 x 10(-22) centimeters per molecule.