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

Theoretical Assessment of the Ligand/Metal/Quadruplex Recognition in the Non-Canonical Nucleic Acids Structures

Quadruplexes (GQs), peculiar DNA/RNA motifs concentrated in specific genomic regions, play a vital role in biological processes including telomere stability and, hence, represent promising targets for anticancer therapy. GQs are formed by folding guanine-rich sequences into square planar G-tetrads which stack onto one another. Metal cations, most often potassium, further stabilize the architecture by coordinating the lone electron pairs of the O atoms. The presence of additional nucleic acid bases, however, has been recently observed experimentally and contributes substantially to the structural heterogeneity of quadruplexes. Therefore, it is of paramount significance to understand the factors governing the underlying complex processes in these structures. The current study employs DFT calculations to model the interactions between metal cations (K+, Na+, Sr2+) and diverse tetrads composed of a guanine layer in combination with a guanine (G)-, adenine (A)-, cytosine (C)-, thymine (T)-, or uracil (U)-based tetrad layer. Moreover, the addition of 4-(3,4-dihydroisoquinolin-2-yl)-2-(quinolin-2-yl)quinazoline to the modeled quadruplexes as a possible mechanism of its well-exerted antitumor effect is assessed. The calculations imply that the metal cation competition and ligand complexation are influenced by the balance between electronic and implicit/explicit solvation effects, the composition of the tetrad layers, as well as by the solvent exposure to the surrounding environment expressed in terms of different dielectric constant values. The provided results significantly enhance our understanding of quadruplex diversity, ligand recognition, and the underlying mechanisms of stabilization at an atomic level.

Momentum-dependent sum-frequency vibrational spectroscopy of bonded interface layer at charged water interfaces

Interface-specific hydrogen (H)-bonding network of water directly controls the energy transfer and chemical reaction pathway at many charged aqueous interfaces, yet to characterize these bonded water layer structures remains a challenge. We now develop a sum-frequency spectroscopic scheme with varying photon momenta as an all-optic solution for retrieving the vibrational spectra of the bonded water layer and the ion diffuse layer and, hence, microscopic structural and charging information about an interface. Application of the method to a model surfactant-water interface reveals a hidden weakly donor H-bonded water species, suggesting an asymmetric hydration-shell structure of fully solvated surfactant headgroups. In another application to a zwitterionic phosphatidylcholine lipid monolayer-water interface, we find a highly polarized bonded water layer structure associating to the phosphatidylcholine headgroup, while the diffuse layer contribution is experimentally proven to be negligible. Our all-optic method offers an in situ microscopic probe of electrochemical and biological interfaces and the route toward future imaging and ultrafast dynamics studies.

Solvation structures of calcium and magnesium ions in water with the presence of hydroxide: a study by deep potential molecular dynamics

The solvation structures of calcium (Ca²⁺) and magnesium (Mg²⁺) ions with the presence of hydroxide (OH^-) ion in water are essential for understanding their roles in biological and chemical processes but have not been fully explored. Ab initio molecular dynamics (AIMD) is an important tool to address this issue, but two challenges exist. First, an accurate description of OH^- from AIMD needs an appropriate exchange-correlation functional. Second, a long trajectory is needed to reach an equilibrium state for the Ca²⁺-OH^- and Mg²⁺-OH^- ion pairs in aqueous solutions. Herein, we adopt a deep potential molecular dynamics (DPMD) method to simulate 1 ns trajectories for the Ca²⁺-OH^- and Mg²⁺-OH^- ion pairs in water; the DPMD method provides efficient machine-learning-based models that have the accuracy of the SCAN exchange-correlation functional within the framework of density functional theory. The solvation structures of the cations and the OH^- in terms of three different species have been systematically investigated. On the one hand, we find that OH^- have more significant effects on the solvation structure of Ca²⁺ than that of Mg²⁺. We observe that the OH^- substantially affects the orientation angles of water molecules surrounding the cation. Through the time correlation functions, we conclude that the water molecules in the first solvation shell of Ca²⁺ change their preferred orientation faster than those of Mg²⁺. On the other hand, with the presence of the cation in the first solvation shell of OH^-, we find that the hydrogen bonds of OH^- are severely altered, and the adjacent water molecules of OH^- are squeezed. The two cations have substantially different effects on the solvation structure of OH^-. Our work provides new insight into the solvation structures of Ca²⁺ and Mg²⁺ in water with the presence of OH^-.

Homopolymer nanobowls with controlled size and denting degree

Homopolymer nanobowls hold promising potential applications in many fields because of their designability, large specific surface area and high packing density. However, it is still challenging to prepare nanobowls with...

Ground-State Structures of Hydrated Calcium Ion Clusters From Comprehensive Genetic Algorithm Search

We searched the lowest-energy structures of hydrated calcium ion clusters Ca2+(H2O)n (n = 10-18) in the whole potential energy surface by the comprehensive genetic algorithm (CGA). The lowest-energy structures of Ca2+(H2O)10-12 clusters show that Ca2+ is always surrounded by six H2O molecules in the first shell. The number of first-shell water molecules changes from six to eight at n = 12. In the range of n = 12-18, the number of first-shell water molecules fluctuates between seven and eight, meaning that the cluster could pack the water molecules in the outer shell even though the inner shell is not full. Meanwhile, the number of water molecules in the second shell and the total hydrogen bonds increase with an increase in the cluster size. The distance between Ca2+ and the adjacent water molecules increases, while the average adjacent O-O distance decreases as the cluster size increases, indicating that the interaction between Ca2+ and the adjacent water molecules becomes weaker and the interaction between water molecules becomes stronger. The interaction energy and natural bond orbital results show that the interaction between Ca2+ and the water molecules is mainly derived from the interaction between Ca2+ and the adjacent water molecules. The charge transfer from the lone pair electron orbital of adjacent oxygen atoms to the empty orbital of Ca2+ plays a leading role in the interaction between Ca2+ and water molecules.

Infrared photodissociation spectroscopic investigation on VO+ and NbO+ hydrolysis catalyzed by water molecules

We investigate the hydrolysis of vanadium/niobium monoxide cation (VO⁺/NbO⁺) with water molecules in the gas phase.

Local vibrational mode analysis of ion-solvent and solvent-solvent interactions for hydrated Ca2+ clusters

Hydrated calcium ion clusters have received considerable attention due to their essential role in biological processes such as bone development, hormone regulation, blood coagulation, and neuronal signaling. To better understand the biological role of the cation, the interactions between the Ca²⁺ ions and water molecules have been frequently investigated. However, a quantitative measure for the intrinsic Ca-O (ion-solvent) and intermolecular hydrogen bond (solvent-solvent) interactions has been missing so far. Here, we report a topological electron density analysis and a natural population analysis to analyze the nature of these interactions for a set of 14 hydrated calcium clusters via local mode stretching force constants obtained at the ωB97X-D/6-311++G(d,p) level of theory. The results revealed that the strength of inner Ca-O interactions for Ca(H₂O)_n ²⁺ (n = 1-8) clusters correlates with the electron density. The application of a second hydration shell to Ca(H₂O)_n ²⁺ (n = 6-8) clusters resulted in stronger Ca-O interactions where a larger electron charge transfer between lp(O) of the first hydration shell and the lower valence of Ca prevailed. The strength of the intermolecular hydrogen bonds, formed between the first and second hydration shells, became stronger when the charge transfers between hydrogen bond (HB) donors and HB acceptors were enhanced. From the local mode stretching force constants of implicitly and explicitly solvated Ca²⁺, we found the six-coordinated cluster to possess the strongest stabilizations, and these results prove that the intrinsic bond strength measures for Ca-O and hydrogen bond interactions form new effective tools to predict the coordination number for the hydrated calcium ion clusters.

Molecular-level origin of the carboxylate head group response to divalent metal ion complexation at the air-water interface

We exploit gas-phase cluster ion techniques to provide insight into the local interactions underlying divalent metal ion-driven changes in the spectra of carboxylic acids at the air-water interface. This information clarifies the experimental findings that the CO stretching bands of long-chain acids appear at very similar energies when the head group is deprotonated by high subphase pH or exposed to relatively high concentrations of Ca2+ metal ions. To this end, we report the evolution of the vibrational spectra of size-selected [Ca2+·RCO2-]+·(H2O) n=0to12 and RCO2-·(H2O) n=0to14 cluster ions toward the features observed at the air-water interface. Surprisingly, not only does stepwise hydration of the RCO2- anion and the [Ca2+·RCO2-]+ contact ion pair yield solvatochromic responses in opposite directions, but in both cases, the responses of the 2 (symmetric and asymmetric stretching) CO bands to hydration are opposite to each other. The result is that both CO bands evolve toward their interfacial asymptotes from opposite directions. Simulations of the [Ca2+·RCO2-]+·(H2O) n clusters indicate that the metal ion remains directly bound to the head group in a contact ion pair motif as the asymmetric CO stretch converges at the interfacial value by n = 12. This establishes that direct metal complexation or deprotonation can account for the interfacial behavior. We discuss these effects in the context of a model that invokes the water network-dependent local electric field along the C-C bond that connects the head group to the hydrocarbon tail as the key microscopic parameter that is correlated with the observed trends.

Aqueous solvation of Mg(ii) and Ca(ii): A Born-Oppenheimer molecular dynamics study of microhydrated gas phase clusters

The hydration features of [Mg(H₂O)_n]²⁺ and [Ca(H₂O)_n]²⁺ clusters with n = 3-6, 8, 18, and 27 were studied by means of Born-Oppenheimer molecular dynamics simulations at the B3LYP/6-31+G** level of theory. For both ions, it is energetically more favorable to have all water molecules in the first hydration shell when n ≤ 6, but stable lower coordination average structures with one water molecule not directly interacting with the ion were found for Mg²⁺ at room temperature, showing signatures of proton transfer events for the smaller cation but not for the larger one. A more rigid octahedral-type structure for Mg²⁺ than for Ca²⁺ was observed in all simulations, with no exchange of water molecules to the second hydration shell. Significant thermal effects on the average structure of clusters were found: while static optimizations lead to compact, spherically symmetric hydration geometries, the effects introduced by finite-temperature dynamics yield more prolate configurations. The calculated vibrational spectra are in agreement with infrared spectroscopy results. Previous studies proposed an increase in the coordination number (CN) from six to eight water molecules for [Ca(H₂O)_n]²⁺ clusters when n ≥ 12; however, in agreement with recent measurements of binding energies, no transition to a larger CN was found when n > 8. Moreover, the excellent agreement found between the calculated extended X-ray absorption fine structure spectroscopy spectra for the larger cluster and the experimental data of the aqueous solution supports a CN of six for Ca²⁺.

On the aqueous solvation of AsO(OH)3vs. As(OH)3. Born–Oppenheimer molecular dynamics density functional theory cluster studies

BOMD simulations were used to reveal the hydration features of As(OH)₃ and (for the first time) AsO(OH)₃ in aqueous solution.

Sequential water molecule binding enthalpies for aqueous nanodrops containing a mono-, di- or trivalent ion and between 20 and 500 water molecules

Sequential water molecule binding enthalpies, ΔH_n,n-1, are important for a detailed understanding of competitive interactions between ions, water and solute molecules, and how these interactions affect physical properties of ion-containing nanodrops that are important in aerosol chemistry. Water molecule binding enthalpies have been measured for small clusters of many different ions, but these values for ion-containing nanodrops containing more than 20 water molecules are scarce. Here, ΔH_n,n-1 values are deduced from high-precision ultraviolet photodissociation (UVPD) measurements as a function of ion identity, charge state and cluster size between 20-500 water molecules and for ions with +1, +2 and +3 charges. The ΔH_n,n-1 values are obtained from the number of water molecules lost upon photoexcitation at a known wavelength, and modeling of the release of energy into the translational, rotational and vibrational motions of the products. The ΔH_n,n-1 values range from 36.82 to 50.21 kJ mol^-1. For clusters containing more than ∼250 water molecules, the binding enthalpies are between the bulk heat of vaporization (44.8 kJ mol^-1) and the sublimation enthalpy of bulk ice (51.0 kJ mol^-1). These values depend on ion charge state for clusters with fewer than 150 water molecules, but there is a negligible dependence at larger size. There is a minimum in the ΔH_n,n-1 values that depends on the cluster size and ion charge state, which can be attributed to the competing effects of ion solvation and surface energy. The experimental ΔH_n,n-1 values can be fit to the Thomson liquid drop model (TLDM) using bulk ice parameters. By optimizing the surface tension and temperature change of the logarithmic partial pressure for the TLDM, the experimental sequential water molecule binding enthalpies can be fit with an accuracy of ±3.3 kJ mol^-1 over the entire range of cluster sizes.

Experimental measurements of water molecule binding energies for the second and third solvation shells of [Ca(H2O) n ]2+ complexes

Further understanding of the biological role of the Ca²⁺ ion in an aqueous environment requires quantitative measurements of both the short- and long-range interactions experienced by the ion in an aqueous medium. Here, we present experimental measurements of binding energies for water molecules occupying the second and, quite possibly, the third solvation shell surrounding a central Ca²⁺ ion in [Ca(H₂O) _n ]²⁺ complexes. Results for these large, previously inaccessible, complexes have come from the application of finite heat bath theory to kinetic energy measurements following unimolecular decay. Even at n = 20, the results show water molecules to be more strongly bound to Ca²⁺ than would be expected just from the presence of an extended network of hydrogen bonds. For n > 10, there is very good agreement between the experimental binding energies and recently published density functional theory calculations. Comparisons are made with similar data recorded for [Ca(NH₃) _n ]²⁺ and [Ca(CH₃OH) _n ]²⁺ complexes.

Experimental measurements of water molecule binding energies for the second and third solvation shells of [Ca(H2O) n ]2+ complexes

Further understanding of the biological role of the Ca²⁺ ion in an aqueous environment requires quantitative measurements of both the short- and long-range interactions experienced by the ion in an aqueous medium. Here, we present experimental measurements of binding energies for water molecules occupying the second and, quite possibly, the third solvation shell surrounding a central Ca²⁺ ion in [Ca(H₂O) _n ]²⁺ complexes. Results for these large, previously inaccessible, complexes have come from the application of finite heat bath theory to kinetic energy measurements following unimolecular decay. Even at n = 20, the results show water molecules to be more strongly bound to Ca²⁺ than would be expected just from the presence of an extended network of hydrogen bonds. For n > 10, there is very good agreement between the experimental binding energies and recently published density functional theory calculations. Comparisons are made with similar data recorded for [Ca(NH₃) _n ]²⁺ and [Ca(CH₃OH) _n ]²⁺ complexes.

Effects of electronic structure on the hydration of PbNO3(+) and SrNO3(+) ion pairs

Hydration of PbNO3(+) and SrNO3(+) with up to 30 water molecules was investigated with infrared photodissociation (IRPD) spectroscopy and with theory. These ions are the same size, yet the IRPD spectra of these ion pairs for n = 2-8 are significantly different. Bands in the bonded O-H region (∼3000-3550 cm(-1)) indicate that the onset of a second hydration shell begins at n = 5 for PbNO3(+) and n = 6 for SrNO3(+). Spectra for [PbNO3](+)(H2O)2-5 and [SrNO3](+)(H2O)3-6 indicate that the structures of clusters with Pb(ii) are hemidirected with a void in the coordinate sphere. A natural bond orbital analysis of [PbNO3](+)(H2O)5 indicates that the anisotropic solvation of the ion is due to a region of asymmetric electron density on Pb(ii) that can be explained by charge transfer from the nitrate and water ligands into unoccupied p-orbitals on Pb(ii). There are differences in the IRPD spectra of PbNO3(+) and SrNO3(+) with up to 25 water molecules attached. IR intensity in the bonded O-H region is blue-shifted by ∼50 cm(-1) in nanodrops containing SrNO3(+) compared to those containing PbNO3(+), indicative of a greater perturbation of the water H-bond network by strontium. The free O-H stretches of surface water molecules in nanodrops containing 10, 15, 20, and 25 water molecules are red-shifted by ∼3-8 cm(-1) for PbNO3(+) compared to those for SrNO3(+), consistent with more charge transfer between water molecules and Pb(ii). These results demonstrate that the different electronic structure of these ions significantly influences how they are solvated.

Hydrated copper ion chemistry: guided ion beam and computational investigation of Cu2+(H2O)n (n = 7-10) complexes

Cross sections for the threshold collision-induced dissociation of Cu(2+)(H(2)O)(n), where n = 8 - 10, are measured using a guided ion beam tandem mass spectrometer. The primary dissociation pathway is found to be loss of a single water molecule followed by the sequential loss of additional water molecules until n = 8, at which point charge separation to form CuOH(+)(H(2)O)(4) (+) H(+)(H(2)O)(3) is observed to occur at a slightly lower energy than loss of a water molecule. Competition from charge separation prohibits the formation of appreciable amounts of the n = 7 or smaller complexes as reactants in the source. These findings indicate that Cu(2+) has a critical size of 8. Analysis of the data using statistical modeling techniques that account for energy distributions and lifetime effects yields primary and sequential bond dissociation energies (BDEs) for loss of one and two water molecules from n = 8 - 10 complexes as well as the barrier for charge separation from n = 8. More speculative analysis extends the thermochemistry obtained down to n = 5 and 6. Theoretical BDEs are determined from quantum chemical calculations using structures optimized at the B3LYP/6 311(+)G(d,p) level along with the lowest-energy isomers suggested by single point energies at the MP2(full), M06, B3LYP, and B3P86 levels of theory using a 6- 311(+)G(2d,2p) basis set. BDEs at 0K are converted to 298 K thermodynamic values using a rigid rotor/harmonic oscillator approximation. Experimental and theoretical entropies of activation suggest that a third solvent shell forms at n = 9, in accord with previous findings. The present work represents the first experimentally determined hydration enthalpies for the Cu(2+)(H(2)O)n system.

Gas phase structure and reactivity of doubly charged microhydrated calcium(II)-catechol complexes probed by infrared spectroscopy

Doubly charged microhydrated adducts formed from catechol and calcium(II) were produced in the gas phase using electrospray ionization (ESI) appearing as the most important ions in the mass spectra recorded. The gas phase structures of [Ca(catechol)2(H2O)](2+) and [Ca(catechol)2(H2O)2](2+) have been assayed by IR multiphoton dissociation (IRMPD) spectroscopy, recording their vibrational spectra in the 3450-3750 cm(-1) range (OH stretching region) and in the 900-1700 cm(-1) fingerprint spectral region. The agreement between experimental and calculated IR spectra of the selected cluster ions confirmed the suitability of the proposed geometries. In addition, quantum chemical calculations at the B3LYP/6-311+G(d,p) level of theory were performed for [Ca(catechol)2(H2O)](2+) to gain insight into the major routes of dissociation. The results suggest that loss of the water molecule is the lowest energy fragmentation channel followed by charge separation products and neutral loss of one catechol molecule, in agreement with the product ions observed upon collision-induced dissociation (CID).

Microsolvation and hydration enthalpies of CaC₂O₄(H₂O) n (n=0-16) and C₂O₄²⁻(H₂O) n (n=0-14): an ab initio study

We studied hydrated calcium oxalate and its ions at the restricted Hartree-Fock RHF/6-31G* level of theory. Performing a configurational search seems to improve the fit of the HF/6-31G* level to experimental data. The first solvation shell of calcium oxalate contains 13 water molecules, while the first solvation shell of oxalate ion is formed by 14 water molecules. The first solvation shell of Ca(II) is formed by six water molecules, while the second shell contains five. At 298.15 K, we estimate the asymptotic limits (infinite dilution) of the total standard enthalpies of hydration for Ca(II), oxalate ion and calcium oxalate as -480.78, -302.78 and -312.73 kcal mol(-1), resp. The dissociation of hydrated calcium oxalate is an endothermic process with an asymptotic limit of +470.84 kcal mol(-1).

Invited review article: laser vaporization cluster sources

The laser vaporization cluster source has been used for the production of gas phase atomic clusters and metal-molecular complexes for 30 years. Numerous experiments in the chemistry and physics of clusters have employed this source. Its operation is simple in principle, but there are many subtle design features that influence the number and size of clusters produced, as well as their composition, charge state, and temperature. This article examines all aspects of the design of these cluster sources, discussing the relevant chemistry, physics, and mechanical aspects of experimental configurations employed by different labs. The principles detailed here provide a framework for the design and implementation of this source for new applications.

Hydration of calcium oxide surface predicted by reactive force field molecular dynamics

In this work, we present the parametrization of Ca-O/H interactions within the reactive force field ReaxFF, and its application to study the hydration of calcium oxide surface. The force field has been fitted using density functional theory calculations on gas phase calcium-water clusters, calcium oxide bulk and surface properties, calcium hydroxide, bcc and fcc Ca, and proton transfer reactions in the presence of calcium. Then, the reactive force field has been used to study the hydration of the calcium oxide {001} surface with different water contents. Calcium oxide is used as a catalyzer in many applications such as CO(2) sequestration and biodiesel production, and the degree of surface hydroxylation is a key factor in its catalytic performance. The results show that the water dissociates very fast on CaO {001} bare surfaces without any defect or vacancy. The surface structure is maintained up to a certain amount of water, after which the surface undergoes a structural rearrangement, becoming a disordered calcium hydroxyl layer. This transformation is the most probable reason for the CaO catalytic activity decrease.

Coordination numbers of hydrated divalent transition metal ions investigated with IRPD spectroscopy

Hydration of the divalent transition metal ions, Mn, Fe, Co, Ni, Cu, and Zn, with 5-8 water molecules attached was investigated using infrared photodissociation spectroscopy and photodissociation kinetics. At 215 K, spectral intensities in both the bonded-OH and free-OH stretch regions indicate that the average coordination number (CN) of Mn(2+), Fe(2+), Co(2+), and Ni(2+) is ~6, and these CN values are greater than those of Cu(2+) and Zn(2+). Ni has the highest CN, with no evidence for any population of structures with a water molecule in a second solvation shell for the hexa-hydrate at temperatures up to 331 K. Mn(2+), Fe(2+), and Co(2+) have similar CN at low temperature, but spectra of Mn(2+)(H(2)O)(6) indicate a second population of structures with a water molecule in a second solvent shell, i.e., a CN < 6, that increases in abundance at higher temperature (305 K). The propensity for these ions to undergo charge separation reactions at small cluster size roughly correlates with the ordering of the hydrolysis constants of these ions in aqueous solution and is consistent with the ordering of average CN values established from the infrared spectra of these ions.

Ions in size-selected aqueous nanodrops: sequential water molecule binding energies and effects of water on ion fluorescence

The effects of water on ion fluorescence were investigated, and average sequential water molecule binding energies to hydrated ions, M(z)(H(2)O)(n), at large cluster size were measured using ion nanocalorimetry. Upon 248-nm excitation, nanodrops with ~25 or more water molecules that contain either rhodamine 590(+), rhodamine 640(+), or Ce(3+) emit a photon with average energies of approximately 548, 590, and 348 nm, respectively. These values are very close to the emission maxima of the corresponding ions in solution, indicating that the photophysical properties of these ions in the nanodrops approach those of the fully hydrated ions at relatively small cluster size. As occurs in solution, these ions in nanodrops with 8 or more water molecules fluoresce with a quantum yield of ~1. Ce(3+) containing nanodrops that also contain OH(-) fluoresce, whereas those with NO(3)(-) do not. This indirect fluorescence detection method has the advantages of high sensitivity, and both the size of the nanodrops as well as their constituents can be carefully controlled. For ions that do not fluoresce in solution, such as protonated tryptophan, full internal conversion of the absorbed 248-nm photon occurs, and the average sequential water molecule binding energies to the hydrated ions can be accurately obtained at large cluster sizes. The average sequential water molecule binding energies for TrpH(+)(H(2)O)(n) and a doubly protonated tripeptide, [KYK + 2H](2+)(H(2)O)(n), approach asymptotic values of ~9.3 (n ≥ 11) and ~10.0 kcal/mol (n ≥ 25), respectively, consistent with a liquidlike structure of water in these nanodrops.

Infrared spectroscopy of Mn+ (H2O) and Mn2+ (H2O) via argon complex predissociation

Singly and doubly charged manganese-water cations, and their mixed complexes with attached argon atoms, are produced by laser vaporization in a pulsed nozzle source. Complexes of the form Mn(+)(H(2)O)Ar(n) (n = 1-4) and Mn(2+)(H(2)O)Ar(4) are studied via mass-selected infrared photodissociation spectroscopy, detected in the mass channels corresponding to the elimination of argon. Sharp resonances are detected for all complexes in the region of the symmetric and asymmetric stretch vibrations of water. With the guidance of density functional theory computations, specific vibrational band resonances are assigned to complexes having different argon attachment configurations. In the small singly charged complexes, argon adds first to the metal ion site and later in larger clusters to the hydrogens of water. The doubly charged complex has argon only on the metal ion. Vibrations in all of these complexes are shifted to lower frequencies than those of the free water molecule. These shifts are greater when argon is attached to hydrogen and also greater for the dication compared to the singly charged species. Cation binding also causes the IR intensities for water vibrations to be much greater than those of the free water molecule, and the relative intensities are greater for the symmetric stretch than the asymmetric stretch. This latter effect is also enhanced for the dication complex.

Average sequential water molecule binding enthalpies of M(H2O)(19-124)2+ (M = Co, Fe, Mn, and Cu) measured with ultraviolet photodissociation at 193 and 248 nm

The average sequential water molecule binding enthalpies to large water clusters (between 19 and 124 water molecules) containing divalent ions were obtained by measuring the average number of water molecules lost upon absorption of an UV photon (193 or 248 nm) and using a statistical model to account for the energy released into translations, rotations, and vibrations of the products. These values agree well with the trend established by more conventional methods for obtaining sequential binding enthalpies to much smaller hydrated divalent ions. The average binding enthalpies decrease to a value of ~10.4 kcal/mol for n > ~40 and are insensitive to the ion identity at large cluster size. This value is close to that of the bulk heat of vaporization of water (10.6 kcal/mol) and indicates that the structure of water in these clusters may more closely resemble that of bulk liquid water than ice, owing either to a freezing point depression or rapid evaporative cooling and kinetic trapping of the initial liquid droplet. A discrete implementation of the Thomson equation using parameters for liquid water at 0 °C generally fits the trend in these data but provides values that are ~0.5 kcal/mol too low.

Structures, stability, vibration entropy and IR spectra of hydrated calcium ion clusters [Ca(H(2)O)(n)](2+) (n = 1-20, 27): a systematic investigation by density functional theory

The low-lying candidates of hydrated calcium ion clusters, [Ca(H(2)O)(n)](2+) with n = 1-20 and 27, have been extensively sought by using density functional theory (DFT) at BLYP/6-311+G(d,p) level. The results showed that the first hydration shell around the calcium ion was fully occupied by six water molecules, whereas the second hydration shell might be fully occupied with different numbers of water molecules. This just corresponds to different growth patterns of the hydrated calcium ion clusters. Furthermore, we revealed that the vibration entropy contributed to the free energy of an isomer significantly. As a result, the stability of some low-lying candidates at zero-temperature was not maintained at finite temperatures. Therefore, we suggested that, at finite temperatures, the realistic products of [Ca(H(2)O)(n)](2+) should be a mixture of the best candidate and some of metastable isomers for a given cluster size. For a cluster having second and/or third shell of water molecules, we found structural transitions between a low-lying structure and the lowest-energy structure undergoing much lower energy barriers. In addition, the IR spectra of the best candidates were predicted, in which the evolution of hydrogen-bond configurations with the cluster size was revealed.

Effect of magnesium cation on the interfacial properties of aqueous salt solutions

Sodium chloride solutions have been used extensively as a model of seawater in both theoretical and experimental studies of the chemistry of sea salt aerosol. Many groups have found that chloride anions are present at the air-solution interface. This observation has been important for the development of a mechanism for the heterogeneous production of molecular chlorine from chloride in sea salt aerosol. However, while sodium chloride is a major constituent of seawater, it is by no means the only salt present. Seawater contains one Mg(2+) for every eight Na(+). Mg(2+) is naturally occurring in ocean waters from mineral deposits in the Earth's crust and biological sources. Mg(2+) forms a hexahydrate structure, rather than contact ion pairs with chloride anion, and this impacts the ordering of water in solution. In this study, we use molecular dynamics simulations, ab initio calculations, and vibrational sum frequency generation (SFG) spectroscopy to explore the effect of the Mg(2+) cation and its tightly bound solvation shell on the surface propensity of chloride, ion-ion interactions, and water structure of the air-solution interface of concentrated chloride salt solutions. In addition, we provide molecular level details that may be relevant to the heterogeneous reactions of chloride in deliquesced sea salt aerosols. In particular, we show that the presence of the divalent Mg(2+) cation does not modify the surface propensity of chloride compared to Na(+) and hence, its availability to interfacial reaction, although some differences in the behavior of chloride may occur due to specific ion interactions. In this work, we also discuss the SFG free OH band at the surface of salt solutions and conclude that it is often not straightforward to interpret.

Zn2+ has a primary hydration sphere of five: IR action spectroscopy and theoretical studies of hydrated Zn2+ complexes in the gas phase

Complexes of Zn(2+)(H(2)O)(n), where n = 6-12, are examined using infrared photodissociation (IRPD) spectroscopy, blackbody infrared radiative dissociation (BIRD), and theory. Geometry optimizations and frequency calculations are performed at the B3LYP/6-311+G(d,p) level along with single point energy calculations for relative energetics at the B3LYP, B3P86, and MP2(full) levels with a 6-311+G(2d,2p) basis set. The IRPD spectrum of Zn(2+)(H(2)O)(8) is most consistent with the calculated spectrum of the five-coordinate MP2(full) ground-state (GS) species. Results from larger complexes also point toward a coordination number of five, although contributions from six-coordinate species cannot be ruled out. For n = 6 and 7, comparisons of the individual IRPD spectra with calculated spectra are less conclusive. However, in combination with the BIRD and laser photodissociation kinetics as well as a comparison to hydrated Cu(2+) and Ca(2+), the presence of five-coordinate species with some contribution from six-coordinate species seems likely. Additionally, the BIRD rate constants show that Zn(2+)(H(2)O)(6) and Zn(2+)(H(2)O)(7) complexes are less stable than Zn(2+)(H(2)O)(8). This trend is consistent with previous work that demonstrates the enthalpic favorability of the charge separation process forming singly charged hydrated metal hydroxide and protonated water complexes versus loss of a water molecule for complexes of n ≤ 7. Overall, these results are most consistent with the lowest-energy structures calculated at the MP2(full) level of theory and disagree with those calculated at B3LYP and B3P86 levels.