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.     

The vaporization rate of ice at temperatures near its melting point

The first study of free vaporization kinetics of ice at temperatures near its melting point is reported. The experimental approach employed is based on a unique combination of thermal desorption spectroscopy, microcalorimetry, and time-of-flight mass spectrometry, making it possible to overcome challenges associated with the introduction of volatile solids into a high vacuum environment. Measurements of the vaporization rate of polycrystalline ice demonstrate that the vaporization kinetics deviate dramatically from those predicted by a simple mobile precursor mechanism. The vaporization rate follows Arrhenius behavior from -40 to 0 degrees C with an effective activation energy of 50+/-4 kJ/mol, which is significantly higher than the value predicted by the simple mobile precursor mechanism. Extrapolation of earlier measurements conducted below -40 degrees C yields a value of approximately 0.02 at 0 degrees C for the vaporization coefficient alphav. In contrast, experimentally determined vaporization coefficient is found to be 0.7+/-0.3 and shows a weak dependence on temperature up to the bulk melting point. The role of possible surface phase transitions in the mechanisms of release and uptake of H2O and other chemical species by ice surfaces is discussed.     

H/D isotopic exchange between water molecules at ice surfaces

H/D isotopic exchange between H(2)O and D(2)O molecules was studied at the surface of ice films at 90-140 K by the technique of Cs(+) reactive ion scattering. Ice films were deposited on a Ru(0001) substrate in different compositions of H(2)O and D(2)O and in various structures to study the kinetics of isotopic exchange. H/D exchange was very slow on an ice film at 95-100 K, even when H(2)O and D(2)O were uniformly mixed in the film. At 140 K, H/D exchange occurred in a time scale of several minutes on the uniform mixture film. Kinetic measurement gave the rate coefficient for the exchange reaction, k(140 K)=1.6(+/-0.3) x 10(-19) cm(2) molecule(-1) s(-1) and k(100 K)< or =5.7(+/-0.5) x 10(-21) cm(2) molecule(-1) s(-1) and the Arrhenius activation energy, E(a)> or =9.8 kJ mol(-1). Addition of HCl on the film to provide excess protons greatly accelerated the isotopic exchange reaction such that it went to completion very quickly at the surface. The rapid reaction, however, was confined within the first bilayer (BL) of the surface and did not readily propagate to the underlying sublayer. The isotopic exchange in the vertical direction was almost completely blocked at 95 K, and it slowly occurred only to a depth of 3 BLs from the surface at 140 K. Thus, the proton transfer was highly directional. The lateral proton transfer at the surface was attributed to the increased mobility of protonic defects at the molecularly disordered and activated surface. The slow, vertical proton transfer was probably assisted by self-diffusion of water molecules.     

Diffusion kinetics for methanol in polycrystalline ice

Quantitative analyses of the isothermal desorption kinetics from methanol-doped H2O films on Pt(111) reveal that transport kinetics for CH3OH in polycrystalline ice are much slower than previously reported. They also indicate that MeOH displays first-order desorption kinetics with respect to its instantaneous surface concentration below 0.1 mole fraction in ice. These observations allow isothermal desorption rate measurements to be interpreted in terms of a depth profiling analysis providing one-dimensional concentration depth profiles from methanol-doped polycrystalline ice films. Using a straightforward approach to inhibit ice sublimation, transport properties are extracted from the evolution of concentration depth profiles obtained after thermal annealing of binary ice films at high temperature. Heterodiffusion coefficients for methanol in polycrystalline (cubic) ice Ic films are reported for temperatures between 145 and 195 K and for concentrations below 10(-3) mole fraction. Finally, diffusion kinetics for methanol in ice are shown to display a very strong concentration dependence that may contribute, in addition to variations in laboratory samples microstructure, to the disagreements reported in the literature regarding the transport properties of ice.     

Interactions of oxalic acid and ice on Cu surface

The interactions between oxalic acid (C 2H 2O 4) and H 2O on a polycrystalline Cu surface have been investigated by reflection-absorption infrared spectroscopy (RAIRS) and temperature-programmed desorption (TPD) methods. The desorption of H 2O and C 2H 2O 4 was studied; we found that the ice desorption temperature increases with the ice-film thickness. Desorption of the C 2H 2O 4 layer involves a structural modification and sublimation. The H 2O/C 2H 2O 4 and C 2H 2O 4/H 2O interfaces and the codeposited C 2H 2O 4+H 2O were prepared on the Cu surface by varying deposition sequences of gaseous C 2H 2O 4 and H 2O at 155 K. We found that the interaction between ice and C 2H 2O 4 does not lead to the H 2O-induced deprotonation of C 2H 2O 4 in a temperature range 155-283 K. However, H-bonding interactions between H 2O and C 2H 2O 4 can lead to the formation of a metastable oxalic acid-ice complex in the C 2H 2O 4/H 2O and C 2H 2O 4+H 2O systems during the TPD process. Desorption of H 2O from the C 2H 2O 4/H 2O/Cu system is suggested to involve the diffusion of H 2O through the top C 2H 2O 4 layer. H 2O desorption is followed by a rearrangement of C 2H 2O 4 to form a C 2H 2O 4 adlayer on Cu in the C 2H 2O 4+H 2O system. These experimental findings suggest that C 2H 2O 4 is not ionized on snow and ice in the polar boundary layer and at upper tropospheric temperatures ( approximately 240 K).     

Molecular dynamics simulations of D2O ice photodesorption

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) photodissociation of D(2)O in an amorphous D(2)O ice surface at 10, 20, 60, and 90 K, in order to investigate the influence of isotope effects on the photodesorption processes. As for H(2)O, the main processes after UV photodissociation are trapping and desorption of either fragments or D(2)O molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption occurs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D(2)O molecule photodesorption. D(2)O desorption takes places either by direct desorption of a recombined D(2)O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding D(2)O molecule out of the surface by transferring part of its momentum. Desorption probabilities are calculated for photoexcitation of D(2)O in the top four monolayers and are compared quantitatively with those for H(2)O obtained from previous MD simulations of UV photodissociation of amorphous water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The main conclusions are the same, but the average D atom photodesorption probability is smaller than that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas the average OD radical photodesorption probability is larger than that of OH (by about a factor of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for every ice temperature studied. The average D(2)O photodesorption probability is larger than that of H(2)O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger contribution of the D(2)O kick-out mechanism. This is an isotope effect: the kick-out mechanism is more efficient for D(2)O ice, because the D atom formed after D(2)O photodissociation has a larger momentum than photogenerated H atoms from H(2)O, and D transfers momentum more easily to D(2)O than H to H(2)O. The total (OD + D(2)O) yield has been compared with experiments and the total (OH + H(2)O) yield from previous simulations. We find better agreement when we compare experimental yields with calculated yields for D(2)O ice than when we compare with calculated yields for H(2)O ice.     

Diffusion of HDO in pure and acid-doped ice films

In these experiments, a few bilayers of D(2)O were vapor-deposited on a pure crystalline H(2)O ice film or an ice film doped with a small amount of HCl. Upon deposition, H/D isotopic exchange quickly converted the D(2)O layer into an HDO-rich mixture layer. Infrared absorption spectroscopy followed the changes of the HDO from the initial HDO mixture layer to HDO isolated in the H(2)O ice film. This was possible because isolated HDO in H(2)O ice has a unique, sharp peak in the O-D stretch region that can be distinguished from the broad peak due to the initial HDO mixture layer. The absorbance of isolated HDO displayed first-order kinetics and was attributed to diffusion of HDO from the HDO-rich mixture layer into the underlying H(2)O ice film. While negligible diffusion was observed for pure ice films and for ice films with HCl concentrations up to 1 x 10(-4) mole fraction, diffusion of HDO occurred for higher concentrations of (2-20) x 10(-4) mole fraction HCl with a concentration-independent rate constant. The diffusion under these conditions followed Arrhenius behavior for T = 135-145 K yielding E(a) = 25 +/- 5 kJ/mol. The mechanism for the HDO diffusion involves either (i) molecular self-diffusion or (ii) long-range H/D diffusion by a series of multiple proton hop and orientational turn steps. While these spectroscopic results compare favorably with recent studies of molecular self-diffusion in low-temperature ice films, the diffusion results from all the ice film studies at low temperatures (ca. T < 170 K) differ from earlier bulk ice studies at higher temperatures (ca. T > 220 K). A comparison and discussion of the various diffusion studies are included in this report.     

Hydrogen atom formation from the photodissociation of water ice at 193 nm

The TOF spectra of photofragment hydrogen atoms from the 193 nm photodissociation of amorphous ice at 90-140 K have been measured. The spectra consist of both a fast and a slow components that are characterized by average translational energies of 2k(B)T(trans)=0.39+/-0.04 eV (2300+/-200 K) and 0.02 eV (120+/-20 K), respectively. The incident laser power dependency of the hydrogen atom production suggests one-photon process. The electronic excitation energy of a branched cluster, (H(2)O)(6+1), has been theoretically calculated, where (H(2)O)(6+1) is a (H(2)O)(6) cyclic cluster attached by a water molecule with the hydrogen bond. The photoabsorption of this branched cluster is expected to appear at around 200 nm. The source of the hydrogen atoms is attributed to the photodissociation of the ice surface that is attached by water molecules with the hydrogen bond. Atmospheric implications are estimated for the photodissociation of the ice particles (Noctilucent clouds) at 190-230 nm in the region between 80 and 85 km altitude.     

Low-energy electron-stimulated luminescence of thin H2O and D2O layers on Pt(111)

The electron-stimulated luminescence (ESL) from amorphous solid water and crystalline ice films deposited on Pt(111) at 100 K is investigated as a function of the film thickness, incident electron energy (5-1000 eV), isotopic composition, and film structure. The ESL emission spectrum has a characteristic double-peaked shape that has been attributed to a transition between a superexcited state (C) and the dissociative, first excited state (A) in water: C --> A. Comparing the electron-stimulated luminescence and O2 electron-stimulated desorption (ESD) yields versus incident electron energy, we find the ESL threshold is approximately 3 eV higher than the O2 ESD threshold, which is close to the center of the emission spectrum near 400 nm and supports the C --> A assignment for the ESL. For thin films, radiative and nonradiative interactions with the substrate tend to quench the luminescence. The luminescence yield increases with coverage since the interactions with the substrate become less important. The ESL yield from D2O is approximately 4-5 times higher than that from H2O. With use of layered films of H2O and D2O, this sizable isotopic effect on the ESL is exploited to spatially profile the luminescence emission within the ASW films. These experiments show that most of the luminescence is emitted from within the penetration depth of the incident electron. However, the results depend on the order of the isotopes in the film and can be modeled by assuming some migration of the electronically excited states within the film. The ESL is very sensitive to defects and structural changes in solid water, and the emission yield is significantly higher from amorphous films than from crystalline ice.     

An online method combining a thermal conversion elemental analyzer with isotope ratio mass spectrometry for the determination of hydrogen isotope composition and water concentration in geological samples

An online continuous-flow method, combining a thermal conversion elemental analyzer (TC/EA) with isotope ratio mass spectrometry (MS), is evaluated for the determination of both the hydrogen isotope composition and the water concentration of hydrous and nominally anhydrous minerals. The technique involves reduction of hydrous minerals or nominally anhydrous minerals by reaction with glassy carbon at 1450 degrees C in a helium stream. The product gases, H2 and CO, are separated on a gas chromatographic column prior to analysis in the mass spectrometer. Calibration curves for the H concentration analysis were generated from a standard of benzoic acid (C7H6O2) that has an H concentration of 5.0 wt%; the analytical uncertainties were better than +/-0.05% in our runs. Two standards of material with given D values, polyethylene IAEA-CH-7 and biotite NBS-30, were tested for the purpose of calibrating a natural garnet 04BXL02 representing nominally anhydrous minerals. Preheating at 90 degrees C for 12 h was found to be suitable for removing adsorption water on the sample surface. This results in constant D values and total H2O contents for the garnet, with weighted means of -94 +/- 1 and 522 +/- 11 ppm (wt), respectively. The TC/EA-MS technique allows routine analysis of sample sizes as small as 0.01 microL H2O. For natural minerals, absolute reproducibilities for D values are +/-0.5 to +/-2 (1) and relative uncertainties for total H2O concentrations are at levels of +/-1% to +/-3% (1). Therefore, this online method can be used for the quantitative determination of H isotope composition and H2O concentration of either hydrous or anhydrous minerals.     

Interstellar ice surface site modification induced by dicyanoacetylene adsorption

Dicyanoacetylene adsorbed on amorphous ice water at 10 K presents an interaction with the dangling H site and induces a s(4) adsorption site formation due to the restructuring of the ice bulk. Warming up the sample provokes the dicyanoacetylene desorption from the H(2)O ice film, which could be due to the beginning of the ice crystallization process. The desorption activation energy measured by temperature-programmed desorption (E(d) = 42 +/- 5 kJ x mol(-1)) is in good agreement with that calculated (E(d) = 46 kJ x mol(-1)) and gives evidence of a hydrogen-bonded adsorbed state on amorphous ice films.     

Unusual decrease in conductivity upon hydration in acceptor doped, microcrystalline ceria

The impact of hydration on the transport properties of microcrystalline Sm(0.15)Ce(0.85)O(1.925) has been examined. Dense, polycrystalline samples were obtained by conventional ceramic processing and the grain boundary regions were found, by high resolution transmission electron microscopy, to be free of impurity phases. Impedance spectroscopy measurements were performed over the temperature range 250 to 650 °C under dry, H(2)O-saturated, and D(2)O-saturated synthetic air; and over the temperature range 575 to 650 °C under H(2)-H(2)O atmospheres. Under oxidizing conditions humidification by either H(2)O or D(2)O caused a substantial increase in the grain boundary resistivity, while leaving the bulk (or grain interior) properties unchanged. This unusual behavior, which was found to be both reversible and reproducible, is interpreted in terms of the space-charge model, which adequately explains all the features of the measured data. It is found that the space-charge potential increases by 5-7 mV under humidification, in turn, exacerbating oxygen vacancy depletion in the space-charge regions and leading to the observed reduction in grain boundary conductivity. It is proposed that the heightened space-charge potential reflects a change in the relative energetics of vacancy creation in the bulk and at the grain boundary interfaces as a result of water uptake into the grain boundary core. Negligible bulk water uptake is detected under both oxidizing and reducing conditions.     

Hydrogen peroxide formation following the vacuum ultraviolet photodissociation of water ice films at 90 K

The production of gaseous OH radicals from the 300-350 nm photodissociation of H(2)O(2) that was photolytically produced on a water ice surface following the 157 nm photolysis of water ice at 90 K was directly monitored using resonance-enhanced multiphoton ionization. The translational energy distribution estimated by the time-of-flight spectrum of the OH products is represented by a Maxwell-Boltzmann energy distribution with a translational temperature of 3750+/-250 K. The rotational temperature was estimated by a spectral simulation to be 225+/-25 K. Surface defects produced by HCl deposition on the water ice contributed to the higher production rate of H(2)O(2) in the 157 nm photoirradiation of water ice while surface coverage caused by CD(3)OH deposition decreased the H(2)O(2) production rate.     

The interaction of water with Ni(111) and H/Ni(111) studied by TPD and HREELS

We have used temperature-programmed desorption in combination with specular and off-specular high resolution electron energy loss spectroscopy to study the interaction of H(2)O and D(2)O with the bare and hydrogen-covered Ni(111) surface. Our results for the bare metal surface agree with previous reports and we are able to relate two prominent features in vibrational spectra to nuclear motions at the surface. Pre-covering Ni(111) with hydrogen alters both adsorption and desorption of water significantly. The strong H-Ni bond does not allow for isotopic exchange with co-adsorbed D(2)O. Strong resemblance of desorption traces and vibrational spectra of submonolayer coverages on H-covered Ni(111) and multilayers on bare Ni(111) suggests that adsorption of hydrogen makes this nickel surface hydrophobic.     

A new high-temperature multinuclear-magnetic-resonance probe and the self-diffusion of light and heavy water in sub- and supercritical conditions

A high-resolution nuclear-magnetic-resonance probe (500 MHz for 1H) has been developed for multinuclear pulsed-field-gradient spin-echo diffusion measurements at high temperatures up to 400 degrees C. The convection effect on the self-diffusion measurement is minimized by achieving the homogeneous temperature distributions of +/-1 and +/-2 degrees C, respectively, at 250 and 400 degrees C. The high temperature homogeneity is attained by using the solid-state heating system composed of a ceramic (AlN) with high thermal conductivity comparable with that of metal aluminium. The self-diffusion coefficients D for light (1H2O) and heavy (2H2O) water are distinguishably measured at subcritical temperatures of 30-350 degrees C with intervals of 10-25 degrees C on the liquid-vapor coexisting curve and at a supercritical temperature of 400 degrees C as a function of water density between 0.071 and 0.251 gcm3. The D value obtained for 1H2O is 10%-20% smaller than those previously reported because of the absence of the convection effect. At 400 degrees C, the D value for 1H2O is increased by a factor of 3.7 as the water density is reduced from 0.251 to 0.071 gcm3. The isotope ratio D(1H2O)D(2H2O) decreases from 1.23 to approximately 1.0 as the temperature increases from 30 to 400 degrees C. The linear hydrodynamic relationship between the self-diffusion coefficient divided by the temperature and the inverse viscosity does not hold. The effective hydrodynamic radius of water is not constant but increases with the temperature elevation in subcritical water.     

Segregation of hydroxide ions to an ice surface

Hydroxide ions that are initially buried within an ice film segregate to the ice film surface at elevated temperatures. This process was observed by conducting experiments with an ice film constructed with a bottom H(2)O layer and an upper D(2)O layer, with an excess of hydroxide ions trapped at the H(2)O/D(2)O interface as they were generated by Na hydrolysis. The transport of hydroxide ions from the interfacial layer to the surface was examined as a function of time using a low energy sputtering method. The progress of the H/D exchange reaction in surface water molecules was also monitored with the Cs(+) reactive ion scattering technique. At 90 K, only a small portion of buried hydroxide ions moved to the surface in the form of OD(-) species. This was due to hydroxide transport via proton hopping through a D(2)O layer, 3 BL thick, in the surface region. At 135 K, at which point water self-diffusion is active in the ice film, the majority of the buried hydroxide ions segregated to the surface after ～1 h. Both OH(-) and OD(-) species were produced at the surface, at an OH(-)/OD(-) population ratio ≥1. Based on kinetic measurements for the transport of OH(-) and OD(-) species and the H/D exchange of surface water molecules, we concluded that the major transport channel for hydroxide ions in this regime is the migration of molecular hydroxide species. H/D exchange reactions also occur between surface hydroxide ions and water molecules. No evidence was observed for the occurrence of the hop-and-turn process at 135 K, although it is known as an important mechanism of proton transport in ice.     

Glycine-ice nanolayers: morphology and surface energetics

Ultrathin glycine-ice films (nanolayers) have been prepared in ultrahigh vacuum by condensation of H(2)O and glycine at 110 K and 150 K on single crystalline Al(2)O(3) surfaces and have been investigated by temperature programed thermal desorption, x-ray photoelectron spectroscopy, and work function measurements. Various layer architectures have been considered, including glycine-on-ice, ice-on-glycine, and mixed glycine-ice nanolayers. Low coverages of adsorbed glycine molecules on amorphous ice surfaces suppress the amorphous-to-crystalline phase transition in the temperature range 140-160 K in near-surface regions and consequently lead to a lower desorption temperature of H(2)O molecules than from pure ice layers. Thicker glycine overlayers on ice provide a kinetic restriction to H(2)O desorption from the underlying ice layers until the glycine molecules become mobile and develop pathways for water desorption at higher temperature (>170 K). Ice overlayers do not wet glycine film surfaces, but the glycine molecules on ice are sufficiently immobile at 110 K, so that continuous glycine overlayers form. In mixed glycine-ice nanolayers the glycine phase displays hydrophobic behavior and a phase separation takes place, with the accumulation of glycine near the surfaces of the films.     

Mechanism for Solvent Exchange on trans-[Os(en)(2)(eta(2)-H(2))S](2+)

Solvent exchange on trans-[Os(en)(2)(eta(2)-H(2))S](2+) (S = H(2)O, CH(3)CN) has been studied in neat solvent as a function of temperature and pressure by (17)O NMR line-broadening and isotopic labeling experiments (S = H(2)O) and by (1)H NMR isotopic labeling experiments (S = CH(3)CN). Rate constants and activation parameters are as follows for S = H(2)O and CH(3)CN, respectively: k(ex)(298) = 1.59 +/- 0.04 and (2.74 +/- 0.03) x 10(-)(4) s(-)(1); DeltaH() = 72.4 +/- 0.5 and 98.0 +/- 1.4 kJ mol(-)(1); DeltaS() = +1.7 +/- 1.8 and +15.6 +/- 4.9 J mol(-)(1) K(-)(1); DeltaV() = -1.5 +/- 1.0 and -0.5 +/- 1.0 cm(3) mol(-)(1). The present investigation of solvent exchange when compared with a previous study on substitution reactions on the same complexes leads to the conclusion that substitution reactions on these compounds undergo an interchange dissociative, I(d), or dissociative, D, reaction mechanism, where solvent dissociation is the rate-limiting step.     

Pyridine- and phosphonate-containing ligands for stable Ln complexation. Extremely fast water exchange on the Gd(III) chelates

Two novel ligands containing pyridine units and phosphonate pendant arms, with ethane-1,2-diamine (L2) or cyclohexane-1,2-diamine (L3) backbones, have been synthesized for Ln complexation. The hydration numbers obtained from luminescence lifetime measurements in aqueous solutions of the Eu(III) and Tb(III) complexes are q = 0.6 (EuL2), 0.7 (TbL2), 0.8 (EuL3), and 0.4 (TbL3). To further assess the hydration equilibrium, we have performed a variable-temperature and -pressure UV-vis spectrophotometric study on the Eu(III) complexes. The reaction enthalpy, entropy, and volume for the hydration equilibrium EuL <--> EuL(H2O) were calculated to be DeltaH degrees = -(11.6 +/- 2) kJ mol(-1), DeltaS degrees = -(34.2 +/- 5) J mol(-1) K(-1), and = 1.8 +/- 0.3 for EuL2 and DeltaH degrees = -(13.5 +/- 1) kJ mol(-1), DeltaS degrees = -(41 +/- 4) J mol(-1) K(-1), and = 1.7 +/- 0.3 for EuL3, respectively. We have carried out variable-temperature 17O NMR and nuclear magnetic relaxation dispersion (NMRD) measurements on the GdL2(H2O)q and GdL3(H2O)q systems. Given the presence of phosphonate groups in the ligand backbone, a second-sphere relaxation mechanism has been included for the analysis of the longitudinal (17)O and (1)H NMR relaxation rates. The water exchange rate on GdL2(H2O)q, = (7.0 +/- 0.8) x 10(8) s(-1), is extremely high and comparable to that on the Gd(III) aqua ion, while it is slightly reduced for GdL3(H2O)q, = (1.5 +/- 0.1) x 10(8) s(-1). This fast exchange can be rationalized in terms of a very flexible inner coordination sphere, which is slightly rigidified for L3 by the introduction of the cyclohexyl group on the amine backbone. The water exchange proceeds via a dissociative interchange mechanism, evidenced by the positive activation volumes obtained from variable-pressure 17O NMR for both GdL2(H2O)q and GdL3(H2O)q (DeltaV = +8.3 +/- 1.0 and 8.7 +/- 1.0 cm(3) mol(-1), respectively).     

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.