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 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.     

Two-dimensional infrared spectroscopy of neat ice Ih

The OH stretch line shape of ice Ih exhibits distinct peaks, the assignment of which remains controversial. We address this longstanding question using two dimensional infrared (2D IR) spectroscopy of the OH stretch of H(2)O and the OD stretch of D(2)O of ice Ih at T = 80 K. The isotropic response is dominated by a 2D line shape component which does not depend on the pump pulse frequency. The decay time of the component that does depend on the pump frequency is calculated using singular value decomposition (bi-exponential decay H(2)O: 30 fs, 490 fs; D(2)O: 40 fs, 690 fs). The anisotropic contribution exhibits on-diagonal peaks, which decay on a very fast timescale (H(2)O: 85 fs; D(2)O: 65 fs), with no corresponding anisotropic cross-peaks. Both isotropic and anisotropic results indicate that randomization of excited dipoles occurs with a very rapid rate, just like in neat liquid water. We conclude that the underlying mechanism relates to the complex interplay between exciton migration and exciton-phonon coupling.     

Inelastic neutron scattering study of water in hydrated LTA-type zeolites

The vibrational dynamics of water molecules encapsulated in synthetic Na-A and Mg-exchanged A zeolites were studied versus temperature by inelastic neutron scattering (INS) measurements (30-1200 cm(-1)) as a function of the induced ion-exchange percentage by using the indirect geometry tof spectrometer TOSCA at the ISIS pulse neutron facility (RAL, UK). The experimental INS spectra were compared with those of ice Ih to characterize the structural changes induced by confinement on the H2O hydrogen-bonded network. We observed, after increasing the Mg2+ content, a tendency of water molecules to restore the bulklike arrangements together with more hindered dynamics. These results are confirmed by the analysis of the evaluated one-phonon amplitude-weighted proton vibrational density of states aimed, in particular, to follow the evolution of the water molecules librational mode region.     

Effects of equilibrium H-bond distance and angle changes on Raman intensities from water

The H-bond energy dispersion over the inhomogeneously broadened OD stretching contour from dilute HDO in H(2)O was determined from absolute Raman intensities; it displays a large minimum near omega=2440 cm(-1) from short, strong H bonds (in agreement with the peak omega from lda ice) and a large maximum near 2650-2675 cm(-1) due to extremely weak or broken H bonds (in agreement with the peak omega from dense, supercritical HDO in H(2)O, 0.9 g/cm(3), 673 K). The difference between extrema is the maximum H-bond DeltaE, 5100+/-500 cal/mol, in excellent agreement with Pauling's limiting value. A pressure of 1500 bars yields an additional maximum and shoulder between the two dispersion extrema from pure water; saturated NaCl in water shows the additional maximum. The maxima near 3350 cm(-1) (1500 bar) or near 3360 cm(-1) (NaCl-H(2)O) arise from bent H bonds; 3350 cm(-1) (1500 bar) corresponding to an angle of approximately 170 degrees in the joint frequency/bend, probability of Lawrence and Skinner. Rising omega refers to a higher probability of larger O-O distances, bent H bonds, and H-bond weakening and breakage. A approximately 50-80 cm(-1) difference between the 2727 cm(-1) OD peak from HDO in steam, and the 2650-2675 cm(-1) dispersion maximum is explained via the very broad approximately 60 cm(-1) liquid peak observed at 342 degrees C and 2000 bar.     

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.     

Isotope effects in ice Ih: a path-integral simulation

Ice Ih has been studied by path-integral molecular dynamics simulations, using the effective q-TIP4P/F potential model for flexible water. This has allowed us to analyze finite-temperature quantum effects in this solid phase from 25 to 300 K at ambient pressure. Among these effects we find a negative thermal expansion of ice at low temperatures, which does not appear in classical molecular dynamics simulations. The compressibility derived from volume fluctuations gives results in line with experimental data. We have analyzed isotope effects in ice Ih by considering normal, heavy, and tritiated water. In particular, we studied the effect of changing the isotopic mass of hydrogen on the kinetic energy and atomic delocalization in the crystal as well as on structural properties such as interatomic distances and molar volume. For D(2)O ice Ih at 100 K we obtained a decrease in molar volume and intramolecular O-H distance of 0.6% and 0.4%, respectively, as compared to H(2)O ice.     

H-bond network in amino acid cocrystals with H2O or H2O2. The DFT study of serine-H2O and serine-H2O2

The structure, IR spectrum, and H-bond network in the serine-H(2)O and serine-H(2)O(2) crystals were studied using DFT computations with periodic boundary conditions. Two different basis sets were used: the all-electron Gaussian-type orbital basis set and the plane wave basis set. Computed frequencies of the IR-active vibrations of the titled crystals are quite different in the range of 10-100 cm(-1). Harmonic approximation fails to reproduce IR active bands in the 2500-2800 frequency region of serine-H(2)O and serine-H(2)O(2). The bands around 2500 and 2700 cm(-1) do exist in the anharmonic IR spectra and are caused by the first overtone of the OH bending vibrations of H(2)O and a combination vibration of the symmetric and asymmetric bendings of H(2)O(2). The quantum-topological analysis of the crystalline electron density enables us to describe quantitatively the H-bond network. It is much more complex in the title crystals than in a serine crystal. Appearance of water leads to an increase of the energy of the amino acid-amino acid interactions, up to ~50 kJ/mol. The energy of the amino acid-water H-bonds is ~30 kJ/mol. The H(2)O/H(2)O(2) substitution does not change the H-bond network; however, the energy of the amino acid-H(2)O(2) contacts increases up to 60 kJ/mol. This is caused by the fact that H(2)O(2) is a much better proton donor than H(2)O in the title crystals.     

Cis-trans isomeric and polymorphic effects on the magnetic properties of water-bridged copper coordination chains

The synthesis and magnetic characterization of vanillin-based Cu(II) mononuclear complexes of formula [Cu(van)(2)(H(2)O)(2)](H(2)O)(x) (van = vanillinate; x = 0, compound 1; x = 2, compounds 2 and 3) were performed. Despite the presence of very similar [Cu(van)(2)(H(2)O)(2)] moieties, the crystal structures exhibit distinct Cu···Cu contacts and display three different through-H-bond exchange-coupling pathways. As a result of the relative positions of the water molecules, the experimental (MAGSUS) exchange-coupling constants are dissimilar, i.e., J(1) = -3.0 cm(-1) (the data have been fitted to the Bleaney-Bowers equation considering a dimer; 2J = -6.0 cm(-1)), J(2) = -4.0 cm(-1) (the data have been fitted to the Bonner-Fischer equation for a chain of monomeric copper(II) units), whereas compound 3 is paramagnetic. Subsequently, the theoretical density functional theory (DFT) and wave function theory-based (DDCI) calculations were carried out to better understand the role of the water molecule as a mediator of the magnetic coupling. The use of localized orbitals allows one to elucidate the role of the H-bonds in generating exchange interactions. Since the exchange-coupling constants are strongly dependent on the mechanisms selectively introduced, the role of the H-bond is demonstrated.     

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.     

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.     

Correlation between the hydrogen-bond structures and the C=O stretching frequencies of carboxylic acids as studied by density functional theory calculations: theoretical basis for interpretation of infrared bands of carboxylic groups in proteins

Carboxylic groups (COOH) of Asp and Glu side chains often function as key components in enzymatic reactions, and identifying their H-bond structures in the active sites is essential for understanding the reaction mechanisms. In this study, the correlation between the H-bond structures and the C=O stretching (nuC=O) frequencies of COOH groups was studied using density functional theory calculations. The nuC=O frequencies and their shifts upon OH deuteration were calculated for model complexes of acetic acid and propionic acid H bonded at different sites with various compounds. Calculation results together with some experimental data showed that, upon direct H bonding at the C=O group, the nuC=O frequencies downshift from the free value (1770-1780 cm(-1) in an Ar matrix) to 1745-1760 cm(-1), while H bonding at the OH hydrogen induce even larger downshifts to provide the frequencies at 1720-1745 cm(-1). In contrast, when the COH oxygen is H-bonded, the nuC=O frequencies upshift to 1785-1800 cm(-1). In double and multiple H-bond forms, H-bonding effects at individual sites are basically additive, and complexes in which the C=O and the OH hydrogen are simultaneously H bonded exhibit significantly low nuC=O frequencies at 1725-1700 cm(-1), while complexes H bonded at the oxygen of the COH in addition to either at the C=O or the OH hydrogen exhibit medium frequencies of 1740-1765 cm(-1). The nuC=O frequencies linearly correlate with the C=O lengths, which are changed by H bonding at different sites. Upon OH deuteration, all the complexes showed nuC=O downshifts mostly by approximately 10 cm(-1) and in some cases as large as approximately 20 cm(-1), and hence deuteration-induced downshifts can be a good indicator, irrespective of H-bond forms, for assignments of the nuC=O bands of carboxylic groups. The results in this study provide the criteria for determining the H-bond structures of Asp and Glu side chains in proteins using their nuC=O bands in Fourier transform infrared spectra.     

Competing interactions in the tetranuclear spin cluster [Ni[(OH)(2)Cr(bispictn)](3)]I(5).5H(2)O. An inelastic neutron scattering and magnetic study

The properties of the spin state manifold of the tetranuclear cluster Ni[(OH)(2)Cr(bispictn)](3)]I(5).5H(2)O (bispictn = N,N'-bis(2-pyridylmethyl)-1,3-propanediamine) are investigated by combining magnetic susceptibility and magnetization measurements with an inelastic neutron scattering (INS) study on an undeuterated sample of Ni[(OH)(2)Cr(bispictn)](3)]I(5).5H(2)O. The temperature dependence of the magnetic susceptibility indicates an S = (1)/(2) ground state, which requires antiferromagnetic interactions both between Cr(3+) and Ni(2+) ions and among the Cr(3+) ions. INS reveals potential single-ion anisotropies to be negligibly small and enables an accurate determination of the exchange parameters. The best fit to the experimental energy level diagram is obtained by an isotropic spin Hamiltonian H = J(CrNi)(S(1)().S(4)() + S(2)().S(4)() + S(3)().S(4)()) + J(CrCr)(S(1)().S(2)() + S(1)().S(3)() + S(2)().S(3)()) with J(CrNi) = 1.47 cm(-)(1) and J(CrCr) = 1.25 cm(-)(1). With this model, the experimental intensities of the observed INS transitions as well as the temperature dependence of the magnetic data are reproduced. The resulting overall antiferromagnetic exchange is rationalized in terms of orbital exchange pathways and compared to the situation in oxalato-bridged clusters.     

High pressure-temperature Raman measurements of H2O melting to 22 GPa and 900 K

The melting curve of H(2)O has been measured by in situ Raman spectroscopy in an externally heated diamond anvil cell up to 22 GPa and 900 K. The Raman-active OH-stretching bands and the translational modes of H(2)O as well as optical observations are used to directly and reliably detect melting in ice VII. The observed melting temperatures are higher than previously reported x-ray measurements and significantly lower than recent laser-heating determinations. However, our results are in accord with earlier optical determinations. The frequencies and intensities of the OH-stretching peaks change significantly across the melting line while the translational mode disappears altogether in the liquid phase. The observed OH-stretching bands of liquid water at high pressure are very similar to those obtained in shock-wave Raman measurements.     

Synthesis and physical properties of K4[Fe(C5O5)2(H2O)2](HC5O5)2·4H2O (C5O5(2-) = croconate): a rare example of ferromagnetic coupling via H-bonds

The reaction of the croconate dianion (C(5)O(5))(2-) with a Fe(III) salt has led, unexpectedly, to the formation of the first example of a discrete Fe(II)-croconate complex without additional coligands, K(4)[Fe(C(5)O(5))(2)(H(2)O)(2)](HC(5)O(5))(2)·4H(2)O (1). 1 crystallizes in the monoclinic P2(1)/c space group and presents discrete octahedral Fe(II) complexes coordinated by two chelating C(5)O(5)(2-) anions in the equatorial plane and two trans axial water molecules. The structure can be viewed as formed by alternating layers of trans-diaquabis(croconato)ferrate(II) complexes and layers containing the monoprotonated croconate anions, HC(5)O(5)(-), and noncoordinated water molecules. Both kinds of layers are directly connected through a hydrogen bond between an oxygen atom of the coordinated dianion and the protonated oxygen atom of the noncoordinated croconate monoanion. A H-bond network is also formed between the coordinated water molecule and one oxygen atom of the coordinated croconate. This H-bond can be classified as strong-moderate being the O···O bond distance (2.771(2) ?) typical of moderate H-bonds and the O-H···O bond angle (174(3)°) typical of strong ones. This H-bond interaction leads to a quadratic regular layer where each [Fe(C(5)O(5))(2)(H(2)O)(2)](2-) anion is connected to its four neighbors in the plane through four equivalent H-bonds. From the magnetic point of view, these connections lead to an S = 2 quadratic layer. The magnetic properties of 1 have been reproduced with a 2D square lattice model for S = 2 ions with g = 2.027(2) and J = 4.59(3) cm(-1). This model reproduces quite satisfactorily its magnetic properties but only above the maximum. A better fit is obtained by considering an additional antiferromagnetic weak interlayer coupling constant (j) through a molecular field approximation with g = 2.071(7), J = 2.94(7) cm(-1), and j = -0.045(2) cm(-1) (the Hamiltonian is written as H = -JS(i)S(j)). Although this second model might still be improved since there is also an extra contribution due to the presence of ZFS in the Fe(II) ions, it confirms the presence of weak ferromagnetic Fe-Fe interactions through H-bonds in compound 1 which represents one of the rare examples of ferromagnetic coupling via H-bonds.     

Magnetic polyoxometalates: anisotropic exchange interactions in the moiety of [(NaOH2)Co3(H2O)(P2W15O56)2]17-

The magnetic exchange interactions in a C0(3)(11) moiety encapsulated in Na(17) [(NaOH(2))Co(3)(H(2)O)(P(2)W(15)O(56))(2)] (NaCo(3)) were studied by a combination of magnetic measurements (magnetic susceptibility and low-temperature magnetization), with a detailed Inelastic Neutron Scattering (INS) investigation. The novel structure of the salt was determined by X-ray crystallography. The ferromagnetic Co(3)O(14) triangular cluster core consists of three octahedrally oxo-coordinated Co(II) ions sharing edges. According to the single-ion anisotropy and spin-orbit coupling usually assumed for octahedral Co(II) ions, the appropiate exchange Hamiltonian to describe the ground-state properties of the isosceles triangular Co(3) spin cluster is anisotropic and is expressed as H = - 2sigma(alpha)(=)(x,y,z)(J(alpha)(12)S(1alpha)S(2alpha) + J(alpha)(23)S(2alpha)S(3alpha) + J(alpha)(13)S(1alpha)S(3alpha)), where J(alpha) are the components of the exchange interactions between the Co(II) ions. To reproduce the INS data, nonparallel anisotropic exchange tensors needed to be introduced, which were directly connected to the molecular symmetry of the complex. The following range of parameters (value +/- 0.5 cm(-1)) was found to reproduce all experimental information while taking magnetostructural relations into account: J(x)(12) = J(y)(13) = 8.6 cm(-1); J(y)(12) = J(x)(13) = 1.4 cm(-1); J(z)(12) = J(z)(13) = 10.0 cm(-1); J(x)(23) = J(y)(23) = 6.5 cm(-1) and = 3.4 cm(-1).     

Inelastic neutron scattering and DFT study of 2-amino-3-hydroxymethyl-1,3-propane diol (TRIS)

We report an inelastic neutron scattering (INS) study of 2-amino-3-hydroxymethyl-1,3-propane diol (TRIS). The assignment of the experimental vibrational spectra measured using several incident neutrons’ energies on HRMECS spectrometer has been made by means of DFT calculations. To simulate crystal environment both molecular cluster and solid state models were used. The study has been completed by an alternative approach, molecular dynamics (MD) calculations, done at the same level of the DFT theory. The INS spectra calculated with the solid state models (normal mode analysis, and MD) gave a better fit of the experiment than the cluster model. On the other hand, the peaks between 650 and 850 cm⁻¹ in the experimental INS spectra assigned to OH torsional modes were reproduced better by the cluster calculations. The nature of the stretching frequency of unusually long O–H bond (1.012 Å) was interpreted by means of MD calculations. The interpretation of the spectrum below 100 cm⁻¹ was based on Fourier transform of the velocity autocorrelation function of centre of mass of a molecule of TRIS.     

Theoretical study of the magnetic behavior of [Fe8] and [Fe16] wheels

The reaction of the trimetallic species [Fe(3)O(PhCOO)(6)(H(2)O)(3)]NO(3) with 1,1,1-tris(hydroxymethyl)ethane (H(3)thme) affords either the octametallic species [Fe(8)(PhCOO)(12)(thme)(4)] 1 or the hexadecametallic species [Fe(16)(EtO)(4)(PhCOO)(16)(Hthme)(12)](NO(3))(4) 2, depending on the nature of the solvent used for crystallization. The structure of 1 can be described as a nonplanar wheel of eight Fe(III) ions bridged by a combination of PhCOO(-) and thme(3)(-) ligands, and 2 as a nonplanar wheel of sixteen Fe(III) ions bridged by PhCOO(-), Hthme(2)(-), and EtO(-) ligands. Both compounds can be broken down into simple units of two metal ions and the bridging ligands that connect them. The best fits of the chi vs T curves in the 300-10 K temperature range were obtained with the parameters g = 2.0, J(1) = -24.0 cm(-1), and J(2) = -8.59 cm(-1) for [Fe(8)] and g = 2.0, J(1) = -25.0 cm(-1), J(2) = -11.73 cm(-1), and J(3) = -69.3 cm(-1) for [Fe(16)]. Density functional theory (DFT) calculations show that the antiferromagnetic interactions between the metals in the dinuclear units decrease when two types of bridging ligands are present, as expected for an orbital counter-complementarity effect.     

Influence of monolayer amounts of HNO3 on the evaporation rate of H2O over ice in the range 179 to 208 K: a quartz crystal microbalance study

The evaporation flux J(ev) of H2O from thin H2O ice films containing between 0.5 and 7 monolayers of HNO3 has been measured in the range 179 to 208 K under both molecular and stirred flow conditions in isothermal experiments. FTIR absorption of the HNO3/H2O condensate revealed the formation of metastable alpha-NAT (HNO(3).3H2O) converting to stable beta-NAT at 205 K. After deposition of HNO3 for 16-80 s on a 1 mum thick pure ice film at a deposition rate in the range (6-60) x 10(12) molecules s(-1) the initial evaporative flux J(ev)(H2O) was always that of pure ice. J(ev)(H2O) gradually decreased with the evaporation of H2O and the concomitant increase of the average mole fraction of HNO3, chi(HNO3), indicating the presence of an amorphous mixture of H2O/HNO3 that is called complexed or (c)-ice whose vapor pressure is that of pure ice. The final value of J(ev) was smaller by factors varying from 2.7 to 65 relative to pure ice. Depending on the doping conditions and temperature of the ice film the pure ice thickness d(D) of the ice film for which J(ev) < 0.85J(ev)(pure ice) varied between 130 and 700 nm compared to the 1000 nm thick original ice film at 208 and 191 K, respectively, in what seems to be an inverse temperature dependence. There exist three different types of H2O molecules under the present experimental conditions, namely (a) free H2O corresponding to pure ice, (b) complexed H2O or c-ice, and (c) H2O molecules originating from the breakup of NAT or amorphous H2O/HNO3 mixtures. The significant decrease of J(ev)(H2O) with increasing chi(HNO3) leads to an increase of the evaporative lifetime of atmospheric ice particles in the presence of HNO3 and may help explain the occurrence of persistent and/or large contaminated ice particles at certain atmospheric conditions.     

Fast thermal desorption spectroscopy study of morphology and vaporization kinetics of polycrystalline ice films

Fast thermal desorption spectroscopy was used to investigate the vaporization kinetics of thin (50-100 nm) H(2)O(18) and HDO tracer layers from 2-5 microm thick polycrystalline H(2)O(16) ice films at temperatures ranging from -15 to -2 degrees C. The isothermal desorption spectra of tracer species demonstrate two distinct peaks, alpha and beta, which we attribute to the vaporization of H(2)O(18) initially trapped at or near the grain boundaries and in the crystallites of the polycrystalline ice, respectively. We show that the diffusive transport of the H(2)O(18) and HDO tracer molecules in the bulk of the H(2)O(16) film is slow as compared to the film vaporization. Thus, the two peaks in the isothermal spectra are due to unequal vaporization rates of H(2)O(18) from grain boundary grooves and from the crystallites and, therefore, can be used to determine independently the vaporization rate of the single crystal part of the film and rate of thermal etching of the film. Our analysis of the tracer vaporization kinetics demonstrates that the vaporization coefficient of single crystal ice is significantly greater than those predicted by the classical vaporization mechanism at temperatures near ice melting point. We discuss surface morphological dynamics and the bulk transport phenomena in single crystal and polycrystalline ice near 0 degrees C.