Lifetimes of cagelike water clusters immersed in bulk liquid water: a molecular dynamics study on gas hydrate nucleation mechanisms

Molecular dynamics simulations were performed to observe the evolution of cagelike water clusters immersed in bulk liquid water at 250 and 230 K. Totally, we considered four types of clusters--dodecahedral (5(12)) and tetrakaidecahedral (5(12)6(2)) cagelike water clusters filled with or without a methane molecule, respectively. The lifetimes of these clusters were calculated according to their Lindemann index (delta) using the criterion of delta> or =0.07. The lifetimes of the clusters at 230 K are longer than that at 250 K, and their ratios are the same as the ratio of structure relaxation times of bulk water at these temperatures. For both the filled and empty clusters, the lifetimes of 5(12)6(2) cagelike clusters are similar to that of 5(12) cagelike clusters. Although the methane molecules indeed make the filled cagelike water clusters live longer than the empty ones, the empty cagelike water clusters still have the chance of being long lived. These observations support the cluster nucleation hypothesis for the formation mechanisms of gas hydrates.     

Melting, freezing, sublimation, and phase coexistence in sodium chloride nanocrystals

Calorimetry measurements, performed by multicollision induced dissociation, have been used to probe the melting of a number of (NaCl)nNa+ clusters with n=22-37. The clusters anneal at 225-325 K and melt at 750-850 K. (NaCl)22Na+ and (NaCl)37Na+, which can adopt geometries that are perfect fragments of the bulk lattice melt at around 850 K. The other clusters, which (except for n=31) must have defects, melt at temperatures which are up to 100 K lower than the perfect nanocrystals. The internal energy distributions become bimodal near the melting temperature. This is the signature of slow dynamic phase coexistence where clusters spontaneously jump back and forth between the solid and liquid states with an average period that is longer than required for thermal equilibration. The jump frequency must be between 10(4) and 10(7) s(-1) for the bimodal distribution to be observable in our experiments. The (NaCl)nNa+ clusters can dissociate by an unusual thermally activated process where melting and freezing raise the internal energy to generate hot solid clusters that can sublime before they cool to the ambient temperature.     

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.     

Mimicking solvent shells in the gas phase. II. Solvation of K+

The observed gas-phase coordination number of K+ in K+(H2O)m clusters is smaller than that observed in bulk solution, where the coordination number has been reported to be between 6 and 8. Both theoretical and gas-phase studies of K+(H2O)m cluster ions point to a coordination number closer to 4. In the gas phase, the coordination number is determined by a variety of factors-the most critical being the magnitude of the K+...ligand pairwise interaction. Decreasing the magnitude of the ion...ligand interaction allows more ligands to directly interact with the cation. One method for decreasing the ion...ligand interaction in K+(H2O)m clusters is to systematically substitute weakly bound ligands for the more strongly bound water molecules. The systematic introduction of para-difluorobenzene (DFB) to K+(H2O)m clusters was monitored using infrared photodissociation spectroscopy in the OH stretching region. By varying the ratio of DFB molecules to water molecules present in K+(H2O)m(DFB)n clusters, the observed coordination number of gas-phase K+ was increased to 8, similar to that reported for bulk solution.     

Temperature and structural changes of water clusters in vacuum due to evaporation

This paper presents a study on evaporation of pure water clusters. Molecular dynamics simulations between 20 ns and 3 micros of clusters ranging from 125 to 4096 molecules in vacuum were performed. Three different models (SPC, TIP4P, and TIP5P) were used to simulate water, starting at temperatures of 250, 275, and 300 K. We monitored the temperature, the number of hydrogen bonds, the tetrahedral order, the evaporation, the radial distribution functions, and the diffusion coefficients. The three models behave very similarly as far as temperature and evaporation are concerned. Clusters starting at a higher temperature show a higher initial evaporation rate and therefore reach the point where evaporation stop (around 240 K) sooner. The radius of the clusters is decreased by 0.16-0.22 nm after 0.5 micros (larger clusters tend to decrease their radius slightly more), which corresponds to around one evaporated molecule per nm(2). The cluster temperature seems to converge towards 215 K independent of cluster size, when starting at 275 K. We observe only small structural changes, but the clusters modeled by TIP5P show a larger percentage of molecules with low diffusion coefficient as t-->infinity, than those using the two other water models. TIP4P seems to be more structured and more hydrogen bonds are formed than in the other models as the temperature falls. The cooling rates are in good agreement with experimental results, and evaporation rates agree well with a phenomenological expression based on experimental observations.     

Water adsorption on Rh(111) at 20 K: from monomer to bulk amorphous ice

The adsorption of water (D(2)O) molecules on Rh(111) at 20 K was investigated using infrared reflection absorption spectroscopy (IRAS). At the initial stage of adsorption, water molecules exist as monomers on Rh(111). With increasing water coverage, monomers aggregate into dimers, larger clusters (n = 3-6), and two-dimensional (2D) islands. Further exposure of water molecules leads to the formation of three-dimensional (3D) water islands and finally to a bulk amorphous ice layer. Upon heating, the monomer and dimer species thermally migrate on the surface and aggregate to form larger clusters and 2D islands. Based on the temperature dependence of OD stretching peaks, we succeeded in distinguishing water molecules inside 2D islands from those at the edge of 2D islands. From the comparison with the previous vibrational spectra of water clusters on other metal surfaces, we conclude that the number of water molecules at the edge of 2D islands is comparable with that of water molecules inside 2D islands on the Rh(111) surface at 20 K. This indicates that the surface migration of water molecules on Rh(111) is hindered as compared with the cases on Pt(111) and Ni(111) and thus the size of 2D islands on Rh(111) is relatively small.     

Protonated adenine: Tautomers,solvated clusters,and dissociation mechanisms

Ab initio and density functional theory calculations at the B3-MP2 and CCSD(T)/6-311 + G(3df,2p) levels of theory are reported that address the protonation of adenine in the gas phase, water clusters, and bulk aqueous solution. The calculations point to N-1-protonated adenine (1+) as the thermodynamically most stable cationic tautomer in the gas phase, water clusters, and bulk solution. This strongly indicates that electrospray ionization of adenine solutions produces tautomer 1+ with a specificity calculated as 97-90% in the 298-473 K temperature range. The mechanisms for elimination of hydrogen atoms and ammonia from 1+ have also been studied computationally. Ion 1+ is calculated to undergo fast migrations of protons among positions N-1, C-2, N-3, N-10, N-7, and C-8 that result in an exchange of five hydrogens before loss of a hydrogen atom forming adenine cation radical at 415 kJ mol(-1) dissociation threshold energy. The elimination of ammonia is found to be substantially endothermic requiring 376-380 kJ mol(-1) at the dissociation threshold and depending on the dissociation pathway. The overall dissociation is slowed by the involvement of ion-molecule complexes along the dissociation pathways. The competing isomerization of 1+ proceeds by a sequence of ring opening, internal rotations, imine flipping, ring closures, and proton migrations to effectively exchange the N-1 and N-10 atoms in 1+, so that either can be eliminated as ammonia. This mechanism explains the previous N-1/N-10 exchange upon collision-induced dissociation of protonated adenine.     

Unusual properties of water at hydrophilic/hydrophobic interfaces

The behaviour of water at mosaic hydrophilic/hydrophobic surfaces of different silicas and in biosystems (biomacromolecules, yeast cells, wheat seeds, bone and muscular tissues) was studied in different dispersion media over wide temperature range using 1H NMR spectroscopy with layer-by-layer freezing-out of bulk water (close to 273 K) and interfacial water (180 < T < 273 K), thermally stimulated depolarization current (TSDC) (90 < T < 270 K), infrared (IR) spectroscopy, and quantum chemical methods. Bulk water and water bound to hydrophilic/hydrophobic interfaces can be assigned to different structural types. There are (i) weakly associated interfacial water (1H NMR chemical shift delta(H) = 1.1-1.7 ppm) that can be assigned to high-density water (HDW) with collapsed structure (CS), representing individual molecules in hydrophobic pockets, small clusters and interstitial water with strongly distorted hydrogen bonds or without them, and (ii) strongly associated interfacial water (delta(H) = 4-5 ppm) with larger clusters, nano- and microdomains, and continuous interfacial layer with both HDW and low-density water (LDW). The molecular mobility of weakly associated bound water is higher (because hydrogen bonds are distorted and weakened and their number is smaller than that for strongly associated water) than that of strongly associated bound water (with strong hydrogen bonds but nevertheless weaker than that in ice Ih) that results in the difference in the temperature dependences of the 1H NMR spectra at T < 273 K. These different waters are also appear in changes in the IR and TSDC spectra.     

Structured water in partially dehydrated yeast cells and at partially hydrophobized fumed silica surface

Nonfreezable water structured due to interaction with endocellular objects in yeast cells (endocellular water) or with partially hydrophobic fumed silica (interfacial water) was studied by means of (1)H NMR spectroscopy with layer-by-layer freezing-out of bulk water and quantum chemical methods applied to water clusters in the gas and liquid (chloroform and cyclohexene) phases and adsorbed on silylated silica. Variation in cell hydration as well as in amount of water adsorbed on modified fumed silica leads to changes in the ratio between contributions of two water states characterized by NMR chemical shifts at delta(H)=1.1-1.7 and 4-5 ppm. Lowering of hydration and temperature results in an enhancement of the first signal, and the opposite result is observed for the second signal. These effects may be explained by structured water distribution in the form of relatively large nanodroplets (delta(H)=4-5 ppm is close to that for bulk water) and small clusters of the 2D structure, in which the fraction of hydrogen atoms out of the hydrogen bonds (delta(H)=1.1-1.7 ppm) is larger than that in nanodroplets.     

Heat capacity of water in nanopores

Heat capacity of controlled amounts of water in Vycor's 2 nm radius pores has been determined in real time during the course of water's isothermal nanoconfinement from bulk state at 358 K, by using temperature-modulated calorimetry. As water transfers from bulk to nanopores via the vapor phase, its heat capacity per molecule increases asymptotically toward a limiting value of 1.4 times the heat capacity of bulk water for 1.8 wt % water in Vycor and 1.04 times for 10.0 wt %. The observations indicate that vibrational and configurational contributions to the heat capacity are highest when the amount of water is insufficient to completely cover the pore wall, and they decrease as more water is present in the nanopores and water clusters form. The heat capacity of water in completely filled nanopores approaches the value for bulk water, thus indicating that the heat capacity varies with the water molecules' position in the nanopores.     

Far-infrared absorption of water clusters by first-principles molecular dynamics

Based on first-principle molecular dynamic simulations, we calculate the far-infrared spectra of small water clusters (H(2)O)(n) (n = 2, 4, 6) at frequencies below 1000 cm(-1) and at 80 K and at atmospheric temperature (T>200 K). We find that cluster size and temperature affect the spectra significantly. The effect of the cluster size is similar to the one reported for confined water. Temperature changes not only the shape of the spectra but also the total strength of the absorption, a consequence of the complete anharmonic nature of the classical dynamics at high temperature. In particular, we find that in the frequency region up to 320 cm(-1), the absorption strength per molecule of the water dimer at 220 K is significantly larger than that of bulk liquid water, while tetramer and hexamer show bulklike strengths. However, the absorption strength of the dimer throughout the far-infrared region is too small to explain the measured vapor absorption continuum, which must therefore be dominated by other mechanisms.     

Free energy of solvation of simple ions: molecular-dynamics study of solvation of Cl- and Na+ in the ice/water interface

Molecular-dynamics simulations of Cl(-) and Na(+) ions are performed to calculate ionic solvation free energies in both bulk simple point-charge/extended water and ice 1 h at several different temperatures, and at the basal ice 1 h/water interface. For the interface we calculate the free energy of "transfer" of the ions across the ice/water interface. For the ions in bulk water in the NPT ensemble at 298 K and 1 atm, results are found to be in good agreement with experiments, and with other simulation results. Simulations performed in the NVT ensemble are shown to give equivalent solvation free energies, and this ensemble is used for the interfacial simulations. Solvation free energies of Cl(-) and Na(+) ions in ice at 150 K are found to be approximately 30 and approximately 20 kcal mol(-1), respectively, less favorable than for water at room temperature. Near the melting point of the model the solvation of the ions in water is the same (within statistical error) as that measured at room temperature, and in the ice is equivalent and approximately 10 kcal mol(-1) less favorable than the liquid. The free energy of transfer for each ion across ice/water interface is calculated and is in good agreement with the bulk observations for the Cl(-) ion. However, for the model of Na(+) the long-range electrostatic contribution to the free energy was more negative in the ice than the liquid, in contrast with the results observed in the bulk calculations.     

X-ray absorption spectroscopy of liquid methanol microjets: bulk electronic structure and hydrogen bonding network

We have measured the X-ray absorption (XA) spectrum of liquid (298 K) methanol at the oxygen and carbon K edges. The 4a(1) orbital at the O K edge exhibits a pronounced sensitivity to the formation of intermolecular hydrogen bonds, with significant differences observed between the vapor and bulk spectra, whereas the C K edge reveals only subtle corresponding spectral changes. Comparison with DFT computed spectra of model methanol clusters indicates that the bulk liquid comprises long chains (n > 6) and rings of hydrogen-bonded monomers.     

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.     

Effect of solvation on the vertical ionization energy of thymine: from microhydration to bulk

The effect of hydration on the vertical ionization energy (VIE) of thymine was characterized using equation-of-motion ionization potential coupled-cluster (EOM-IP-CCSD) and effective fragment potential (EFP) methods. We considered several microsolvated clusters as well as thymine solvated in bulk water. The VIE in bulk water was computed by averaging over solvent-solute configurations obtained from equilibrium molecular dynamics trajectories at 300 K. The effect of microsolvation was analyzed and contrasted against the combined effect of the first solvation shell in bulk water. Microsolvation reduces the ionization energy (IE) by about 0.1 eV per water molecule, while the first solvation shell increases the IE by 0.1 eV. The subsequent solvation lowers the IE, and the bulk value of the solvent-induced shift of thymine's VIE is approximately -0.9 eV. The combined effect of the first solvation shell was explained in terms of specific solute-solvent interactions, which were investigated using model structures. The convergence of IE to the bulk value requires the hydration sphere of approximately 13.5 ? radius. The performance of the EOM-IP-CCSD/EFP scheme was benchmarked against full EOM-IP-CCSD using microhydrated structures. The errors were found to be less than 0.01-0.02 eV. The relative importance of the polarization and higher multipole moments in EFP model was also investigated.     

Interaction of water with LiCl, LiBr, and LiI in the deeply supercooled region

The hydration mechanism of lithium halides was studied using time-of-flight secondary ion mass spectrometry as a function of temperature. The lithium halides embedded in thin films of amorphous solid water segregate to the surface at temperatures higher than 135-140 K, with efficiency increasing in the order of LiCl, LiBr, and LiI. A monolayer of LiCl and LiI adsorbed on the surface of amorphous solid water tends to diffuse into the bulk at 160 K. The infrared absorption band revealed that the aqueous lithium-halide solutions and crystals are formed simultaneously at 160 K; these phenomena are explicable as a consequence of the evolution of supercooled liquid water. The strong surfactant effect is inferred to arise from hydration of a contact ion pair having hydrophilic (lithium) and hydrophobic (halide) moieties. Furthermore, bulk diffusion of lithium halides might result from the formation of a solvent-separated ion pair in supercooled liquid water. The presence of two liquid phases of water with different local structures is probably responsible for the formation of these two hydrates, consistent with the calculated result reported by Jungwirth and Tobias[J. Phys. Chem. B 106, 6361 (2002)].     

Single-walled carbon nanotubes filled with M OH (M = K, Cs) and then washed and refilled with clusters and molecules

Heating single-walled carbon nanotubes (SWNTs) with molten hydroxides MOH (M = K, Cs) gave MOH@SWNT in good yield; high resolution transmission electron microscopy (HRTEM) indicated that CsOH in CsOH@SWNT often adopts twisted 1D crystal structures inside SWNTs; treating MOH@SWNT with water at room temperature removes the soluble hydroxide filling and the resulting SWNTs may then be filled using aqueous solutions of uranyl acetate or uranyl nitrate at rt giving SWNTs filled with UO(2) clusters and uranyl acetate molecules.     

Molecular dynamics studies of ion distributions for DNA duplexes and DNA clusters: salt effects and connection to DNA melting

We present extensive molecular dynamics simulations of the ion distributions for DNA duplexes and DNA clusters using the Amber force field with implicit water. The distribution of ions and the electrostatic energy of ions around an isolated DNA duplex and clusters of DNA duplexes in different salt (NaCl) concentrations over the range 0.2-1.0 mol/L are determined on the basis of the simulation results. Using the electrostatic energy profile, we determine a local net charge fraction phi, which is found to increase with increasing of salt concentration. For DNA clusters containing two DNA duplexes (DNA pair) or four DNA duplexes, phi increases as the distance between the duplexes decreases. Combining this result with experimental results for the dependence of the DNA melting temperature on bulk salt concentration, we conclude that for a pair of DNA duplexes the melting temperature increases by 5-10 K for interaxis separations of 25-40 A. For a cluster of four DNA duplexes, an even larger melting temperature increase should occur. We argue that this melting temperature increase in dense DNA clusters is responsible for the cooperative melting mechanism in DNA-linked nanoparticle aggregates and DNA-linked polymer aggregates.     

Structural and electric field effects of ions in aqueous nanodrops

Ensemble infrared photodissociation (IRPD) spectra in the hydrogen stretch region (～2950-3800 cm(-1)) are reported for M(H(2)O)(35-37), with M = I(-), Cl(-), HCO(3)(-), OH(-), tetrabutyl-, tetrapropyl-, and tetramethylammonium, Cs(+), Na(+), Li(+), H(+), Ba(2+), Ca(2+), Co(2+), Mg(2+), La(3+), and Tm(3+), at 133 K. A single, broad feature is observed in the bonded-OH region of the spectra that indicates that the water network in these clusters is bulk-like and likely resembles liquid water more strongly than ice. The free-OH region for all of these clusters is dominated by peaks corresponding to water molecules that accept two and donate one hydrogen bond (AAD water molecules), indicating that AAD water molecules are more abundant at the surface of these ions than AD water molecules. A-only water molecules are present in significant abundance only for the trivalent metal cations. The frequency of the AAD free-OH stretch band shifts nearly linearly with the charge state of the ion, consistent with a Stark shift attributable to the ion's electric field. From these data, a frequency range of 3704.9-3709.7 cm(-1) is extrapolated for the free-OH of AAD water molecules at the (uncharged) bulk liquid water surface, consistent with sum-frequency generation spectroscopy experiments. Differences in both the bonded- and the free-OH regions of the spectra for these ions are attributable to ion-induced patterning of the water network that extends to the surface of the clusters, which includes water molecules in the third and fourth solvation shells; that is, these ions pattern water molecules at long distance to various extents. These spectra are simulated using two different electrostatic models previously used to calculate OH-stretch spectra of bulk water and aqueous solutions and parametrized for bonded-OH frequencies. These models qualitatively reproduce a number of features in the experimental spectra, although it is evident that more sophisticated treatment of water molecule and ion polarizability and vibrational coupling is necessary for more quantitative comparisons.