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.     

Wetting of mixed OHH(2)O layers on Pt(111)

We describe the effect of growth temperature and OHH(2)O composition on the wetting behavior of Pt(111). Changes to the desorption rate of ice films were measured and correlated to the film morphology using low energy electron diffraction and thermal desorption of chloroform to measure the area of multilayer ice and monolayer OHH(2)O exposed. Thin ice films roughen, forming bare (radical39 x radical39)R16 degrees water monolayer and ice clusters. The size of the clusters depends on growth temperature and determines their kinetic stability, with the desorption rate decreasing when larger clusters are formed by growth at high temperature. Continuous films of more than approximately 50 layers thick stabilize an ordered incommensurate ice film that does not dewet. OH coadsorption pins the first layer into registry with Pt, forming an ordered hexagonal (OH+H(2)O) structure with all the H atoms involved in hydrogen bonding. Although this layer has a similar honeycomb OH(x) skeleton to ice Ih, it is unable to reconstruct to match the bulk ice lattice parameter and does not form a stable wetting layer. Water aggregates to expose bare monolayer (OH+H(2)O), forming bulk ice crystallites whose size depend on preparation temperature. Increasing the proportion of water in the first layer provides free OH groups which stabilize the multilayer. The factors influencing multilayer wetting are discussed using density functional theory calculations to compare water adsorption on top of (OH+H(2)O) and on simple models for commensurate water structures. We show that both the (OH+H(2)O) structure and "H-down" water layers are poor proton acceptors, bonding to the first layer being enhanced by the presence of free OH groups. Formation of an ordered ice multilayer requires a water-metal interaction sufficient to wet the surface, but not so strong as to prevent the first layer relaxing to stabilize the interface between the metal and bulk ice.     

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.     

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.     

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.     

A reexamination of the ice III/IX hydrogen bond ordering phase transition

Ice III is a hydrogen bond disordered crystal which when cooled 1 K / min or faster transforms to an antiferroelectric hydrogen bond ordered structure, ice IX. Throughout its region of stability, experiments indicate that the H bonds in ice III are, in fact, partially ordered, i.e., some proton arrangements are preferred. In addition, there has been evidence that the structure of ice IX retains some residual disorder after the transition. Diffraction experiments and calorimetry apparently conflict with regard to the degree of ordering at the ice III/IX transition. Mean field statistical mechanical theories have been used to link partial occupations from diffraction data with thermodynamics. In this work, we investigate the ice III/IX proton ordering phase transition using electronic density functional theory calculations for small unit cells, extended to simulate the phase transition in a large unit cell using graph invariants. In agreement with experiment, we observe partial ordering over a wide range of temperatures as ice III transforms to partially disordered ice IX, near 126 K, which becomes fully ordered at lower temperatures. We compare our results from full statistical mechanical simulations with mean field models, finding small errors for the low-temperature ice IX phase and much larger errors for the high-temperature ice III phase. The failure of mean field theories may explain the apparent conflict between diffraction experiments and calorimetry.     

Prediction of a phase transition to a hydrogen bond ordered form of ice VI

Ice VI is a hydrogen bond disordered crystal over its known region of stability. In this work, we predict that ice VI will transform into a hydrogen bond ordered phase near 108 K, and have identified the likely low-temperature phase as ferroelectric (space group Cc) with an antiferroelectric structure (space group P2(1)2(1)2(1)) close by in energy. Electronic density functional theory calculations provide input to our calculations, which are extended to cells large enough for statistical simulations by using graph invariants. A significant decrease in the configurational entropy is predicted as hydrogen bonds exhibit partial order above the transition, provided that the hydrogen bonds can equilibrate on an experimental time scale. Conversely, partial disorder is predicted at temperatures below the transition. Although some evidence for ordering of ice VI has been observed in experiments, a low-temperature proton ordered phase has not been identified experimentally.     

Formation of high density amorphous ice by decompression of ice VII and ice VIII at 135 K

Monte Carlo computer simulations of ice VII and ice VIII phases have been undertaken using the four-point transferable intermolecular potential model of water. By following thermodynamic paths similar to those used experimentally, ice is decompressed resulting in an amorphous phase. These phases are compared to the high density amorphous phase formed upon compression of ice Ih and are found to have very similar structures. By cooling liquid water along the water/Ih melting line a high density amorphous phase was also generated.     

Phase changes of CO2 hydrate under high pressure and low temperature

High pressure and low temperature experiments with CO(2) hydrate were performed using diamond anvil cells and a helium-refrigeration cryostat in the pressure and temperature range of 0.2-3.0 GPa and 280-80 K, respectively. In situ x-ray diffractometry revealed that the phase boundary between CO(2) hydrate and water+CO(2) extended below the 280 K reported previously, toward a higher pressure and low temperature region. The results also showed the existence of a new high pressure phase above approximately 0.6 GPa and below 1.0 GPa at which the hydrate decomposed to dry ice and ice VI. In addition, in the lower temperature region of structure I, a small and abrupt lattice expansion was observed at approximately 210 K with decreasing temperature under fixed pressures. The expansion was accompanied by a release of water content from the sI structure as ice Ih, which indicates an increased cage occupancy. A similar lattice expansion was also described in another clathrate, SiO(2) clathrate, under high pressure. Such expansion with increasing cage occupancy might be a common manner to stabilize the clathrate structures under high pressure and low temperature.     

Pressure amorphized ices--an atomistic perspective

We offer our viewpoint on the nature of amorphous ices produced by pressurization of crystalline ice Ih and the inter-relationship between them from an atomistic perspective. We argue that the transformation of high density amorphous (HDA) ice from crystalline ice is due to a mechanical process arising from the instability of the ice Ih structure. The densification of HDA upon thermal annealing under pressure is a relaxation process. The conversion of the densified amorphous ice to a lower density form (LDA) upon the release of pressure can be attributed to a similar process. It is speculated that amorphous ices are metastable frustrated structures due to the large activation barriers associated with proton reorientation in the formation of the underlying stable crystalline ice polymorphs.     

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.     

Phase transitions of AlFeO3 and GaFeO3 from the chiral orthorhombic (Pna2(1)) structure to the rhombohedral (R3c) structure

AlFeO(3) and GaFeO(3), which crystallize in a chiral orthorhombic (Pna2(1)) structure, transform to a rhombohedral (R3c) structure when subjected to ball-milling. There is a distinct difference between the transformations of AlFeO(3) and GaFeO(3). AlFeO(3) first transforms to an orthorhombic P2(1)2(1)2(1) structure followed by its transformation to the R3c structure, while GaFeO(3) goes directly to the R3c structure. The transformations have been characterized by X-ray diffraction and Raman spectroscopy. Magnetic properties of Pna2(1) and the transformed phases show significant differences. It is noteworthy that partial substitution of aluminum by gallium in AlFeO(3) as in Al(0.5)Ga(0.5)FeO(3) eliminates the intermediate P2(1)2(1)2(1) phase, causing direct transformation of the Pna2(1) structure to the R3c structure. All of the transformations are thermodynamically first-order associated with significant changes in volume. We have used first-principles simulations to determine the pressure-dependent properties of AlFeO(3) and GaFeO(3) in orthorhombic and corundum structures and have estimated the critical pressures for the structural phase transition between the two structures. On the basis of this information, we also comment on the differences in the behavior of AlFeO(3) and GaFeO(3) under ball-milling.     

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.     

A plastic phase of water from computer simulation

We report a member of ices called plastic or rotator phase, in which individual water molecules make facile rotations as in liquid state but are held tightly in an ordered structure. Molecular dynamics simulations of three classical models of water show that a plastic ice phase appears at a temperature when ice VII is heated or liquid water is cooled at high pressures above several gigapascals. A large amount of latent heat is absorbed when ice VII is transformed to the rotator phase at 590 K and 10 GPa, which is a typical characteristic of the plastic transitions for nearly spherical molecules. In addition to the spontaneous formation of plastic phase in the simulations, its existence is supported by robustness of plastic phase for hypothetical water with varying degrees of Coulombic interactions.     

Hydration structure of water confined between mica surfaces

We report further molecular dynamics simulations on the structure of bound hydration layers under extreme confinement between mica surfaces. We find that the liquid phase of water is maintained down to 2 monolayer (ML) thick, whereas the structure of the K(+) ion hydration shell is close to the bulk structure even under D = 0.92 nm confinement. Unexpectedly, the density of confined water remains approximately the bulk value or less, whereas the diffusion of water molecules decreases dramatically. Further increase in confinement leads to a transition to a bilayer ice, whose density is much less than that of ice Ih due to the formation of a specific hydrogen-bonding network.     

Structural and dynamical properties of solid ammonia borane under high pressure

The structural and dynamical properties of solid ammonia borane were investigated by means of extensive density functional theory calculation up to 60 GPa. Molecular dynamics simulations suggest that the Cmc2(1) phase found by recent room-temperature x-ray diffraction experiments can be obtained from the Pmn2(1) structure at high pressure and low temperature. Two new high-pressure phases were found on further compression at room temperature. We also found that all three high-pressure phases have proton-ordered structures, and the separation of the NH(3) and BH(3) rotation observed in the simulations can be explained by their distinct rotational energy barriers. The role of dihydrogen bonds in the high-pressure phases is discussed.     

Pressure tuning between NH...N hydrogen-bonded ice analogue and NH...Br polar dabcoHBr complexes

At normal conditions 1,4-diazabicyclo[2.2.2]octane hydrobromide [C(6)H(13)N(2)](+.)Br(-) forms centrosymmetric crystals, space group Pm2, NH(+)...N hydrogen-bonded linear polycationic chains with disordered protons in the structure. As in H(2)O ice Ih, the protons in [C(6)H(13)N(2)](+.)Br(-) crystals remain disordered at low temperatures. Above 0.4 GPa the [C(6)H(13)N(2)](+.)Br(-) crystals transform into a new polar NH(+)...Br(-) hydrogen bonded complex, space group Cmc2. It has been crystallized in-situ in a diamond anvil cell and its structure determined by X-rays. The low-pressure triggering of this transformation indicates that it is a possible source of defects in the real structure at normal conditions, where, along with disproportionation defects, they can be responsible for anomalous dielectric properties, including relaxor-like behavior of NH...N hydrogen-bonded compounds.     

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.     

The effect of hydrogen bonding on the excited-state proton transfer in 2-(2'-hydroxyphenyl)benzothiazole: a TDDFT molecular dynamics study

The dynamics of the excited-state proton transfer (ESPT) in a cluster of 2-(2'-hydroxyphenyl)benzothiazole (HBT) and hydrogen-bonded water molecules was investigated by means of quantum chemical simulations. Two different enol ground-state structures of HBT interacting with the water cluster were chosen as initial structures for the excited-state dynamics: (i) an intramolecular hydrogen-bonded structure of HBT and (ii) a cluster where the intramolecular hydrogen bond in HBT is broken by intermolecular interactions with water molecules. On-the-fly dynamics simulations using time-dependent density functional theory show that after photoexcitation to the S(1) state the ESPT pathway leading to the keto form strongly depends on the initial ground state structure of the HBT-water cluster. In the intramolecular hydrogen-bonded structures direct excited-state proton transfer is observed within 18 fs, which is a factor two faster than proton transfer in HBT computed for the gas phase. Intermolecular bonded HBT complexes show a complex pattern of excited-state proton transfer involving several distinct mechanisms. In the main process the tautomerization proceeds via a triple proton transfer through the water network with an average proton transfer time of approximately 120 fs. Due to the lack of the stabilizing hydrogen bond, intermolecular hydrogen-bonded structures have a significant degree of interring twisting already in the ground state. During the excited state dynamics, the twist tends to quickly increase indicating that internal conversion to the electronic ground state should take place at the sub-picosecond scale.