The melting point of ice Ih for common water models calculated from direct coexistence of the solid-liquid interface

In this work we present an implementation for the calculation of the melting point of ice I(h) from direct coexistence of the solid-liquid interface. We use molecular dynamics simulations of boxes containing liquid water and ice in contact. The implementation is based on the analysis of the evolution of the total energy along NpT simulations at different temperatures. We report the calculation of the melting point of ice I(h) at 1 bar for seven water models: SPC/E, TIP4P, TIP4P-Ew, TIP4P/ice, TIP4P/2005, TIP5P, and TIP5P-E. The results for the melting temperature from the direct coexistence simulations of this work are in agreement (within the statistical uncertainty) with those obtained previously by us from free energy calculations. By taking into account the results of this work and those of our free energy calculations, recommended values of the melting point of ice I(h) at 1 bar for the above mentioned water models are provided.     

Structure of coexisting liquid phases of supercooled water: analogy with ice polymorphs

The structural changes occurring in supercooled liquid water upon moving from one coexisting liquid phase to the other have been investigated by computer simulation using a polarizable interaction potential model. The obtained results favorably compare with recent neutron scattering data of high and low density water. In order to assess the physical origin of the observed structural changes, computer simulation of several ice polymorphs has also been carried out. Our results show that there is a strict analogy between the structure of various disordered (supercooled) and ordered (ice) phases of water, suggesting that the occurrence of several different phases of supercooled water is rooted in the same physical origin that is responsible for ice polymorphism.     

Molecular-dynamics studies of surface of ice Ih

We performed molecular dynamics calculations of surface of ice Ih in order to investigate formation mechanism of melting layer on the surface. The results showed that the vibrational amplitude of the atoms in the surface layer greatly depends on the crystal orientation, whereas that in the ice bulk is isotropic. The anisotropy of the vibration is due to a dangling motion of the free O-H bonds exist at the surface layer. The dangling motion enhances the rotational motion of the water molecules. The vibrational density of state showed a coupling between the rotational vibration and the lattice vibration of the water molecules in the surface layer. The coupling of the vibrations causes a distortion of ice lattice. Through the hydrogen-bonding network, the distortion transmits to the interior of the crystal. We conclude that the dangling motion of the free O-H bonds exist at the surface layer is one of the dominant factors governing the surface melting of ice crystal.     

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.     

Quantum path integral simulation of isotope effects in the melting temperature of ice Ih

The isotope effect in the melting temperature of ice Ih has been studied by free energy calculations within the path integral formulation of statistical mechanics. Free energy differences between isotopes are related to the dependence of their kinetic energy on the isotope mass. The water simulations were performed by using the q-TIP4P/F model, a point charge empirical potential that includes molecular flexibility and anharmonicity in the OH stretch of the water molecule. The reported melting temperature at ambient pressure of this model (T=251?K) increases by 6.5±0.5 and 8.2±0.5?K upon isotopic substitution of hydrogen by deuterium and tritium, respectively. These temperature shifts are larger than the experimental ones (3.8 and 4.5 K, respectively). In the classical limit, the melting temperature is nearly the same as that for tritiated ice. This unexpected behavior is rationalized by the coupling between intermolecular interactions and molecular flexibility. This coupling makes the kinetic energy of the OH stretching modes larger in the liquid than in the solid phase. However, the opposite behavior is found for intramolecular modes, which display larger kinetic energy in ice than in liquid water.     

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.     

Transient absorption of vibrationally excited ice Ih

The ultrafast dynamics of HDO:D2O ice Ih at 180 K is studied by midinfrared ultrafast pump-probe spectroscopy. The vibrational relaxation of HDO:D2O ice is observed to proceed via an intermediate state, which has a blueshifted absorption spectrum. Polarization resolved measurements reveal that the intermediate state is part of the intramolecular relaxation pathway of the HDO molecule. In addition, slow dynamics on a time scale of the order of 10-100 ps is observed, related to thermally induced collective reorganizations of the ice lattice. The transient absorption line shape is analyzed within a Lippincott-Schroeder model for the OH-stretch potential. This analysis identifies the main mechanism behind the strong spectral broadening of the v(OH)=1-->2 transition.     

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.     

Observability of coexisting phases of clusters

That two or more phases of small clusters can coexist in thermodynamic equilibrium over ranges of temperature and pressure has become well established. Moreover the explanation for this apparent violation of the Gibbs phase rule is also now well known. The origin of the phenomenon lies entirely with the difference between systems of small numbers of component atoms or molecules and those made of large numbers, e.g., tens, vs. 10²⁰. However little has been said about the maximum sizes of clusters for which such coexistence may be expected to be observable. Here we show how one can estimate that maximum size for observable coexisting phases, in which the unfavored minority phase constitutes a detectable fraction of the total sample. In addition, the role of atom thermal motion in the phase transition is analyzed.     

A van der Waals density functional study of ice Ih

Density functional theory with the van der Waals density functional (vdW-DF) is used to calculate equilibrium crystal structure, binding energy, and bulk modulus of ice Ih. It is found that although it overestimates the equilibrium volume, vdW-DF predicts accurate binding energy of ice Ih, as compared with high level quantum chemistry calculations and experiment. Inclusion of the nonlocal correlation, i.e., van der Waals interaction, leads to an overall improvement over the standard generalized gradient approximation in describing water ice.     

Direct dynamics study of ultrafast vibrational energy relaxation in ice Ih

Car-Parrinello molecular dynamics (CPMD) and a previously developed wave packet model are used to study ultrafast relaxation in water clusters. Water clusters of 15 water molecules are used to represent ice Ih. The relaxation is studied by exciting a symmetric or an asymmetric stretch mode of the central water molecule. The CPMD results suggest that relaxation occurs within 100 fs. This is in agreement with experimental work by Woutersen and Bakker and the earlier wave packet calculations. The CPMD results further indicate that the excitation energy is transferred both intramolecularly and intermolecularly on roughly the same time scale. The intramolecular energy transfer occurs predominantly between the symmetric and asymmetric modes while the bend mode is largely left unexcited on the short time scale studied here.     

Melting temperature: from nanocrystalline to amorphous phase

By extrapolating the mean grain size of nanocrystal to an infinitesimal value, an amorphous phase has been obtained from the Voronoi construction. The molecular dynamics simulations indicated that for nanocrystal, the grain size variation of melting temperature exhibits two characteristic regions. As mean grain size above about 4 nm for Ag, the melting temperatures decrease with decreasing grain size. However, with grain size further shrinking, the melting temperatures almost keep a constant. This is because the dominant factor on the melting temperature of nanocrystal shifts from grain phase to grain boundary. As a result of fundamental difference in structure, the amorphous phase has a much lower solid-to-liquid transformation temperature than that of nanocrystal.     

Melting temperature of trans-polyoctenamer

The effect of microstructure on crystallizability of polyoctenamers prepared by R₃Al-WCl₆ catalyst was studied. The results indicate that polyoctenamers with a broad range of trans-vinylene content do crystallize. The measured melting points are dependent on the trans-vinylene content. From the dependence of melting temperature on copolymer composition, a value of 73 ± 2°C. for the melting point and a molar heat of fusion ΔH_u of 3520 cal./mole are calculated for 100% trans-polyoctenamer. From the melting point depression in the presence of diluent, a value for ΔH_u of 4800 cal./mole is obtained.     

Quasi-harmonic approximation of thermodynamic properties of ice Ih, II, and III

Several thermodynamic properties of ice Ih, II, and III are studied by a quasi-harmonic approximation and compared to results of quantum path integral and classical simulations. This approximation allows to obtain thermodynamic information at a fraction of the computational cost of standard simulation methods, and at the same time permits studying quantum effects related to zero-point vibrations of the atoms. Specifically, we have studied the crystal volume, bulk modulus, kinetic energy, enthalpy, and heat capacity of the three ice phases as a function of temperature and pressure. The flexible q-TIP4P/F model of water was employed for this study, although the results concerning the capability of the quasi-harmonic approximation are expected to be valid independently of the employed water model. The quasi-harmonic approximation reproduces with reasonable accuracy the results of quantum and classical simulations showing an improved agreement at low temperatures (T< 100 K). This agreement does not deteriorate as a function of pressure as long as it is not too close to the limit of mechanical stability of the ice phases.     

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.     

Distribution of coexisting solid and fluid phases alters the kinetics of collapse from phospholipid monolayers

To determine how coexistence of liquid-expanded (LE) and tilted-condensed (TC) phases in phospholipid monolayers affects collapse from the air/water interface, we studied binary films containing dioleoyl phosphatidylcholine-dipalmitoyl phosphatidylcholine (DPPC) mixtures between 10 and 100% DPPC. Previously published results established that this range of compositions represents the LE-TC coexistence region at the equilibrium spreading pressure of 47 mN/m. When held at 49.5 mN/m on a captive bubble, the extent of total collapse fit with the LE area predicted by the phase diagram. The kinetics of collapse, however, when normalized for changes in the LE area, slowed with increasing mole fraction of DPPC. Surface area expressed as stretched exponential functions of time yielded an Avrami exponent that decreased from 1 for the homogeneously LE film to 0.3 for DPPC > or = 70%. Microscopic studies showed that the largest changes in kinetics occurred when either alterations of the initial composition or the process of collapse induced the films to cross the percolation threshold, so that the LE phase became divided into isolated domains. Our results show that although coexisting solid and fluid phases collapse to extents that are independent, the kinetics of collapse, corrected for differences in LE area, depend on the distribution of the two phases.     

Polymer modifies the critical region of the coexisting liquid phases

Several recent conceptual advances, which take advantage of the polymer conformation in the near critical point of coexisting liquid phases and practical techniques of some unique molecular interactions between polymer chain and the solvent molecules, have been made to allow the investigation of the effect of the well-defined polymer in phase separation of binary mixtures. The behavior of a flexible linear or branched chain polymer (polyethylene oxide, PEO, MW = 9 x 10(5), as an impurity) in the critical binary mixture of isobutyric acid (I) + water (W) was studied by the refractive index (n) measurements using a very accurate and sensitive refractometer. The refractive index in each phase of IW as well as three different PEO concentrations (C = 0.395, 0.796, and 1.605 mg/cm(3)) in the near critical composition of IW have been measured at temperatures below the system's upper critical point. We observed that the polymer was significantly affected in the critical region of IW and these various concentrations of PEO show an important behavior on the critical exponents (beta), the critical temperatures (T(c)), and critical composition (phi(c)), which are depicting the shape of the coexistence curve. The phase-transition region of coexisting phases of IW shifts down with the addition of PEO and T(c) decreases linearly with increasing PEO concentrations. This may indicate that the polymer chain entangles with each phase, thereby the polymer monomers strongly interact with neighbor solvent particles and also intrachain interaction between the polymer segments. At such conditions, the collapse of polymer chain is possible in the vicinity of the critical point. At temperatures T close enough to T(c), the critical exponent beta (defined by the relation (n(1) - n(2)) proportional, variant (T(c) - T)(beta), with n(1) and n(2) being the refractive indices of the coexisting phases) was found to decrease from 0.382 to 0.360 when the PEO concentration changes from 0.395 to 1.605 mg/cm(3). These values are higher than that of 0.326 +/- 0.005 of pure IW, which is compatible with the three-dimensional Ising value beta = 0.325. The observed critical exponents for the PEO in IW are fully renormalized Ising critical exponents. Besides, the phi(c) values decrease with increasing the C values in the mixture of IW. It appears that the shape of the PEO in IW coexistence curves is similar from that of pure IW.     

Melting temperature of Pb nanostructural materials from free energy calculation

The thermodynamic properties of lead, including the entropy, heat capacity, Gibbs free energy, and surface free energy have been studied. Based on bulk thermodynamic properties of lead, Gibbs free energy for nanostructural materials is obtained and used to calculate the size-dependent melting point depression for lead nanostructural materials. The studies indicate that the surface free energy difference between solid phase and liquid phase is a decisive factor for the size-dependent melting of nanostructural materials. The calculated results are in agreement with recent experimental values and the available molecular dynamics simulation data.