Investigation of the hydrogen bonding in ice Ih by first-principles density function methods

It is a well recognized difficult task to simulate the vibrational dynamics of ices using the density functional theory (DFT), and there has thus been rather limited success in modelling the inelastic neutron scattering (INS) spectra for even the simplest structure of ice, ice Ih, particularly in the translational region below 400 cm(-1). The reason is partly due to the complex nature of hydrogen bonding (H-bond) among water-water molecules which require considerable improvement of the quantum mechanical simulation methods, and partly owing to the randomness of protons in ice structures which often requires simulation of large super-lattices. In this report, we present the first series of successful simulation results for ice Ih using DFT methods. On the basis of the recent advancement in the DFT programs, we have achieved for the first time theoretical outcomes that not only reproduce the rotational frequencies between 500 to 1200 cm(-1) for ice Ih, but also the two optic peaks at ～240 and 320 cm(-1) in the translational region of the INS spectra [J. C. Li, J. Chem. Phys 105, 6733 (1996)]. Besides, we have also investigated the impact of pairwise configurations of H(2)O molecules on the H-bond and found that different proton arrangements of pairwise H(2)O in the ice Ih crystal lattice could not alter the nature of H-bond as significantly as suggested in an early paper [J. C. Li and D. K. Ross, Nature (London) 365, 327 (1993)], i.e., reproducing the two experimental optic peaks do not need to invoke the two H-bonds as proposed in the previous model which led to considerable debates. The results of this work suggest that the observed optic peaks may be attributed to the coupling between the two bands of H-O stretching modes in H(2)O. The current computational work is expected to shed new light on the nature of the H-bonds in water, and in addition to offer a new approach towards probing the interaction between water and biomaterials for which H-bond is essential.     

Raman spectra of proton ordered phase XI of ICE I. Translational vibrations below 350 cm(-1)

Polarized Raman spectra from single crystals of ice XI (proton ordered phase of ice Ih) were measured and assigned for the modes below 350 cm(-1) in the translational vibration region. In contrast to the proton disordered ice Ih, the spectra in ice XI show clear polarization dependence and several new peaks are observed. Most of the vibrational modes were successfully assigned by the simplified point mass model with the symmetry C(2v) (12)(Cmc2(1)) and by the depolarization effect. In particular, LO-TO splitting of the mode near 240 cm(-1) was experimentally confirmed for the first time, which indicates that the long range force effect appears distinctly in ice XI.     

Vibrational modes of hydrogens in the proton ordered phase XI of ice: Raman spectra above 400 cm(-1)

Polarized Raman spectra of the proton ordered phase of ice Ih, i.e., ice XI, were measured above 400 cm(-1) in the range of librational, bending, and stretching vibrations. Vibrational modes in ice XI, of which symmetry is C(2v) (12)(Cmc2(1)), were discussed from the group theoretical point of view. In the librational mode spectra below 1200 cm(-1), several new peaks and clear polarization dependencies were observed. Assignments of the librational modes agree reasonably well with the recent MD calculations by Iwano et al. (J. Phys. Soc. Jpn. 79, 063601 (2010)). In contrast, the spectra for bands above 1200 cm(-1) show no distinct polarization dependencies and the spectra resemble those in ice Ih. In ice XI, however, fine structure composed of several weak peaks appear on the broad bending and the combination band. No direct evidence of the LO-TO splitting of the ν(3) anti-symmetric stretching mode was obtained. It is contrary to the case of the translational modes Abe and Shigenari (J. Chem. Phys. 134, 104506 (2011)). Present results suggest that the influence of the proton ordering in ice XI is weaker than the effect of inter- and intra-molecular couplings in the stretching vibrations of ice Ih.     

Raman H-bond pair volume for water

The dispersion of the H-bond pair volume Delta V over the decoupled OD and coupled OH-stretching contours from HDO in H(2)O was determined from Raman intensities at pressures to 9700 bar at 301 K. The dispersion of Delta V was determined from -RT[partial differential ln(I(i)/I(REF))/ partial differential P](T) versus omega (in cm(-1)), where i refers to omega's over the stretching contours and I(REF) refers to the reference intensity at the isosbestic frequency. The maximum H-bond pair volume (defined for breakage) is 1.4+/-0.1 cm(3)/mol H-bond, which corresponds to the volume difference between a large dispersion maximum at 2,675 cm(-1) (near the OD stretch omega of HDO in dense supercritical water) and a large, broad minimum centered near 2,375 cm(-1) (just below the OD stretch omega of HDO in lda ice). The average DeltaV is 0.71+/-0.10 cm(3)/mol H-bond. Other minima near 2,625 cm(-1) (OD) and 3550 cm(-1) (OH) refers to bent H-bonds whose angles are approximately 150 deg. Isothermal pressurization of water lowers the molal volume by decreasing the concentration of long, weak H-bonds, and increasing the concentrations of bent H-bonds and short, strong, linear H-bonds. Such bending, shortening, and strengthening produces freezing to ice VI near 10 kbar at 301 K. The isobaric temperature derivative of the maximum H-bond volume is (partial differential Delta V/partial differential T)(P)< or =(2-5) x 10(-3) cm(3)/deg mol H-bond. The OH enthalpy dispersion curve for saturated NaBF(4) in water, yields a large maximum at 3,530-3,540 cm(-1) indicating that BF(4) (-) interacts preferentially with the dangling or "free" OH groups of water thus producing weak, strongly bent H-bonds having angles similar to those of the 3,550 cm(-1) high-pressure H-bond bending feature.     

Sticking of CO to crystalline and amorphous ice surfaces

We present results of classical trajectory calculations on the sticking of hyperthermal CO to the basal plane (0001) face of crystalline ice Ih and to the surface of amorphous ice Ia. The calculations were performed for normal incidence at a surface temperature Ts = 90 K for ice Ia, and at Ts = 90 and 150 K for ice Ih. For both surfaces, the sticking probability can be fitted to a simple exponentially decaying function of the incidence energy, Ei: Ps = 1.0e(-Ei(kJ/mol)/90(kJ/mol)) at Ts = 90 K. The energy transfer from the impinging molecule to the crystalline and the amorphous surface is found to be quite efficient, in agreement with the results of molecular beam experiments on the scattering of the similar molecule, N2, from crystalline and amorphous ice. However, the energy transfer is less efficient for amorphous than for crystalline ice. Our calculations predict that the sticking probability decreases with Ts for CO scattering from crystalline ice, as the energy transfer from the impinging molecule to the warmer surfaces becomes less efficient. At high Ei (up to 193 kJ/mol), no surface penetration occurs in the case of crystalline ice. However, for CO colliding with the amorphous surface, a penetrating trajectory was observed to occur into a large water pore. The molecular dynamics calculations predict that the average potential energy of CO adsorbed to ice Ih is -10.1 +/- 0.2 and -8.4 +/- 0.2 kJ/mol for CO adsorbed to ice Ia. These values are in agreement with previous experimental and theoretical data. The distribution of the potential energy of CO adsorbed to ice Ia was found to be wider (with a standard deviation sigma of 2.4 kJ/mol) than that of CO interacting with ice Ih (sigma = 2.0 kJ/mol). In collisions with ice Ia, the CO molecules scatter at larger angles and over a wider distribution of angles than in collisions with ice Ih.     

Simple computer experiments with ordinary ice

Simple computer experiments in which various fractions of hydrogen bonds (H-bonds) in ice are allowed to break are presented in this paper. First, up to six million water molecules were used to build an artificial piece of ordinary hexagonal ice in the form of a cube, a monolayer, a bilayer, a trilayer, and thicker layers. Then, certain percentages of H-bonds were broken, and the obtained structures were examined. It was found that a large percentage of H-bonds must be broken in order to completely fragment the network of ice into clusters. For a cubic piece of ice, which can be considered bulk ordinary ice, this percentage is equal to 61% H-bonds, a figure also predicted as the threshold of the percolation theory for ice. If, as usually assumed, 13-20% of H-bonds are broken during melting (estimates based on the comparison between the heats of melting and sublimation of ice), the H-bond network of ice is not fragmented and the overwhelming majority of water molecules (>99%) belong to a new, distorted but unbroken network. The percentage of broken H-bonds required for full fragmentation of layers increases with the number of layers and reaches the bulk value of ice for 5-8 layers. This value is consistent with the literature observation that films of water thicker than 20-30 A have properties close to those of the bulk structure.     

Proton Ordering in Cubic Ice. II. Structure of Optimal Configurations

The configurations of the ground energy state of cubic ice obtained previously using the model of strong and weak H-bonds are analyzed. It is established that the stepwise ordered structure is a defectless axial twin. A peculiarity of the twin lies in the fact that strict crystal ordering in its regions is due to the presence of the boundary and its structure. The difference between the structures of the hexagonal layers of ices Ih and Ic is revealed. It is shown that the one-dimensionally disordered layered state in the base planes has antiferroelectric ordering.     

Two-dimensional infrared spectroscopy of isotope-diluted ice Ih

We present experimental 2D IR spectra of isotope diluted ice Ih (i.e., the OH stretch mode of HOD in D(2)O and the OD stretch mode of HOD in H(2)O) at T = 80 K. The main spectral features are the extremely broad 1-2 excited state transition, much broader than the corresponding 0-1 groundstate transition, as well as the presence of quantum beats. We do not observe any inhomogeneous broadening that might be expected due to proton disorder in ice Ih. Complementary, we perform simulations in the framework of the Lippincott-Schroeder model, which qualitatively reproduce the experimental observations. We conclude that the origin of the observed line shape features is the coupling of the OH-vibrational coordinate with crystal phonons and explain the beatings as a coherent oscillation of the O···O hydrogen bond degree of freedom.     

A dispersive model of radical accumulation in irradiated solids

- A dispersive model of radical accumulation in irradiated solids is discussed by the example of hexagonal ice (Ih) and cubic ice (Ic) obtained by annealing of hyperquenched glassy water at 160 K. The model assumes that radical production upon -irradiation is accompanied by their second-order decay, which proceeds according to dispersive kinetics with the time-dependent specific reaction rate k(t) = Bt ^-1, where B = constant and is the dispersion parameter equal to 0.40 for both the hexagonal and cubic ices. The radicals, OH in Ih, and OH plus HO₂ in Ic are produced at the same rate upon -irradiation. The observed enhancement of radical accumulation in ice Ic in comparison with Ih is due to the marked decrease of radical recombination in ice Ic. This is rationalized by the hindrance of radical transport in the lattice of ice Ic containing substantial amounts of Bjerrum-type L-defects.     

Raman spectra of gas hydrates--differences and analogies to ice 1h and (gas saturated) water

It is generally accepted that Raman spectroscopic investigations of gas hydrates provide vital information regarding the structure of the hydrate, hydrate composition and cage occupancies, but most research is focused on the vibrational spectra of the guest molecules. We show that the shape and position of the Raman signals of the host molecules (H(2)O) also contain useful additional information. In this study, Raman spectra (200-4000 cm(-1)) of (mixed) gas hydrates with variable compositions and different structures are presented. The bands in the OH stretching region (3000-3800 cm(-1)), the O-H bending region (1600-1700 cm(-1)) and the O-O hydrogen bonded stretching region (100-400 cm(-1)) are compared with the corresponding bands in Raman spectra of ice Ih and liquid water. The interpretation of the differences and similarities with respect to the crystal structure and possible interactions between guest and host molecules are presented.     

Unusual Gruneisen and Bridgman parameters of low-density amorphous ice and their implications on pressure induced amorphization

The low-temperature limiting value of the Grüneisen parameter for low-frequency phonons and the density dependence of the thermal conductivity (Bridgman parameter) of low-density amorphous (LDA) ice, high-density amorphous (HDA) ice, hexagonal ice Ih, and cubic ice Ic were calculated from high-pressure sound velocity and thermal conductivity measurements, yielding negative values for all states except HDA ice. LDA ice is the first amorphous state to exhibit a negative Bridgman parameter, and negative Grüneisen parameters are relatively unusual. Since Ih, Ic, and LDA ice all transform to HDA upon pressurization at low temperatures and share the unusual feature of negative Grüneisen parameters, this seems to be a prerequisite for pressure induced amorphization. We estimate that the Grüneisen parameter increases at the ice Ih to XI transition, and may become positive in ice XI, which indicates that proton-ordered ice XI does not amorphize like ice Ih on pressurization.     

Low temperature FTIR spectroscopy and hydrogen bonding in cytosine polycrystals

The FTIR spectra of both the pure NH and isotopically substituted ND (<10% and >90% D) polycrystalline cytosine were recorded in the range 400-4000 cm(-1) as a function of temperature (10-300 K). For the first time, uncoupled NH(D) stretching mode bands of amine and imine groups were observed in the spectra of isotopically diluted cytosine at low temperatures. These bands correspond to the three distinct H-bonds that are present in the crystal, in agreement with the available data obtained by structural methods. At least nine bands were observed below 1000 cm(-1) and, in consonance with their temperature and isotopic exchange behavior, were assigned to the NH proton out-of-the-plane bending modes. Six of these bands were found to correspond to additional "disordered" H-bonds, which could not be observed by structural methods. Empirical correlations of spectral and thermodynamic parameters enabled to estimate the contribution of the H-bonds to the sublimation enthalpy of the crystal, in agreement with independent experimental data.     

Simulations of proton order and disorder in ice Ih

Computer simulations of ice Ih with different proton orientations are presented. Simulations of proton disordered ice are carried out using a Monte Carlo method which samples over proton degree of freedom, allowing for the calculation of the dielectric constant and for the examination of the degree of proton disorder. Simulations are also presented for two proton ordered structures of ice Ih, the ferroelectric Cmc2(1) structure or ice XI and the antiferroelectric Pna2(1) structure. These simulations indicate that a transition to a proton ordered phase occurs at low temperatures (below 80 K). The symmetry of the ordered phase is found to be dependent on the water potential. The stability of the two proton ordered structures is due to a balance of short-ranged interactions which tend to stabilize the Pna2(1) structure and longer-range interactions which stabilize the Cmc2(1) structure.     

Low temperature FTIR spectra and hydrogen bonds in polycrystalline cytidine

FTIR spectra of polycrystalline samples of cytidine, pure and containing a small quantity of N(O)H or N(O)D groups (<20%), were measured in KBr pellets from 4000 to 400 cm(-1) at temperatures from 300 to 20K. For the first time the bands of the narrow isotopically decoupled proton stretching vibration mode (nu(1)) of OH- and NH- groups were found; their number corresponds to the number of H-bonds in crystal according to structural data. The FTIR spectra at low temperature in the out-of-plane bending nu(4) proton mode range (lower than 1000 cm(-1)) of N(O)H groups revealed narrow bands, which correspond to nu(1) bands together with several "extra" bands, which are influenced by the isotopic exchange and (or) cooling. All of them have their counterparts in the N(O)D-substance spectrum with an isotopic frequency ratio of 1.30-1.40. The "extra" bands are assigned to the H-bound OH and NH protons, which are disordered and cannot be seen with X-ray crystal structure analysis. The peak positions of both mode bands (expressed as the red shift of nu(1) or blue shift of nu(4) modes relatively free molecules) were used for the estimation of the energy of different H-bonds using previously established empirical correlations between spectral and thermodynamic parameters of hydrogen bonds. The correlation of the red shift and H-bond length is also confirmed for all five H-bonds of cytidine.     

Ultrastructure in frozen/etched saline solutions: on the internal cleansing of ice

Seawater, with its 3.5% salt content, freezes into hexagonal ice (Ih) that encloses concentrated brine within its matrix. When unsubmerged sea ice reaches a certain height and temperature, the brine drains downward through narrow channels. This mechanism was now modeled by frozen 2-3.5% saline as investigated by cryo-etch high-resolution secondary electron microscopy. Thus, saline was either plunge-frozen in liquid ethane at -183 degrees C or else high-pressure frozen to -105 degrees C in 5-6 ms. Ice from a freshly exposed surface was then subjected to a high-vacuum sublimation ("etching"), a procedure that removes pure bulk ice in preference to ice from frozen hydrated salt. After chromium-coating the etched surface with a 2-nm film, the sample was examined by cryo-HRSEM. Granular icy "fences" were seen surrounding empty areas where amorphous ice had originally resided. Since the fences, about 1-2 mum high, survived the etching, it is likely that they consist of frozen brine. The presence of such fences suggests that, during freezing, saline can purge itself of salt with remarkable speed (5-6 ms). Alternatively, channels (perhaps routed around submicroscopic crystallites of cubic ice (Ic) embedded in the amorphous ice at -105 degrees C) can guide the migration of salt to the periphery of ice patches. Macromolecules fail to form fences because they diffuse too slowly or because they are too large to pass through the channels.     

An application of flexible constraints in Monte Carlo simulations of the isobaric--isothermal ensemble of liquid water and ice Ih with the polarizable and flexible mobile charge densities in harmonic oscillators model

The method of flexible constraints was implemented in a Monte Carlo code to perform numerical simulations of liquid water and ice Ih in the constant number of molecules, volume, and temperature and constant pressure, instead of volume ensembles, using the polarizable and flexible mobile charge densities in harmonic oscillators (MCDHO) model. The structural and energetic results for the liquid at T=298 K and rho=997 kg m(-3) were in good agreement with those obtained from molecular dynamics. The density obtained at P=1 atm with flexible constraints, rho=1008 kg m(-3), was slightly lower than with the classical sampling of the intramolecular vibrations, rho=1010 kg m(-3). The comparison of the structures and energies found for water hexamers and for ice Ih with six standard empirical models to those obtained with MCDHO, show this latter to perform better in describing water far from ambient conditions: the MCDHO minimum lattice energy, density, and lattice constants were in good agreement with experiment. The average angle HOH of the water molecule in ice was predicted to be slightly larger than in the liquid, yet 1.2% smaller than the experimental value.     

Proton ordering in cubic ice and hexagonal ice; a potential new ice phase--XIc

Ordinary water ice forms under ambient conditions and has two polytypes, hexagonal ice (Ih) and cubic ice (Ic). From a careful comparison of proton ordering arrangements in Ih and Ic using periodic density functional theory (DFT) and diffusion Monte Carlo (DMC) approaches, we find that the most stable arrangement of water molecules in cubic ice is isoenergetic with that of the proton ordered form of hexagonal ice (known as ice XI). We denote this potential new polytype of ice XI as XIc and discuss a possible route for preparing ice XIc.     

Electrostatic calculation of linear and non-linear optical properties of ice Ih,II, IX and VIII

Macroscopic first- and third-order susceptibilities of ice Ih, ice II, ice IX and ice VIII are calculated using static and frequency-dependent electronic and static vibrational molecular (hyper)polarizabilities at the MP2 level. The molecular properties are in good agreement with experiment and with high-level ab initio calculations. Intermolecular electrostatic and polarization effects due to induced dipoles are taken into account using a rigorous local-field theory. The electric field due to permanent dipoles is used to calculate effective in-crystal (hyper)polarizabilities. The polarizability depends only weakly on the permanent field, but the dipole moment and the hyperpolarizabilities are strongly affected. The calculated linear susceptibility is in good agreement with available experimental data for ice Ih, and the third-order susceptibility for a third harmonic generation experiment is in reasonable agreement with experimental values for liquid water. The molecular vibrational contributions have a small effect on the susceptibilities. The electric properties of a water tetramer are calculated and used to estimate the effect of non-dipolar interactions on the susceptibilities of ice Ih, which are found to be small.     

The bend angle of water in ice Ih and liquid water: The significance of implementing the nonlinear monomer dipole moment surface in classical interaction potentials

The implementation of the physically accurate nonlinear dipole moment surface of the water monomer in the context of the Thole-type, polarizable, flexible interaction potential results in the only classical potential, which, starting from the gas phase value for the bend angle (104.52 degrees), reproduces its experimentally observed increase in the ice Ih lattice and in liquid water. This is in contrast to all other classical potentials to date, which predict a decrease of the monomer bend angle in ice Ih and in liquid water with respect to the gas phase monomer value. Simulations under periodic boundary conditions of several supercells consisting of up to 288 molecules of water used to sample the proton disorder in the ice Ih lattice yield an average value of vartheta(HOH)(I(h))=108.4 degrees +/-0.2 degrees for the minimized structures (T=0 K) and 108.1 degrees +/-2.8 degrees at T=100 K. Analogous simulations for liquid water predict an average value of vartheta(HOH)(liquid)=106.3 degrees +/-4.9 degrees at T=300 K. The increase of the monomer bend angle of water in condensed environments is attributed to the use of geometry-dependent charges that are used to describe the nonlinear character of the monomer's dipole moment surface. Our results suggest a new paradigm in the development of classical interaction potential models of water that can be used to describe condensed aqueous environments.