Classical and quantum gibbs free energies and phase behavior of water using simulation and cell theory

A method to calculate the classical and quantum free energy of a liquid from a computer simulation by using cell theory [J. Chem. Phys. 2007, 126, 064504] is tested for liquid water and ice Ih against experiment as a function of temperature. This fast and efficient method reproduces reasonably well the experimental values of entropy, enthalpy, and free energy of a liquid across the supercooled, stable, and superheated range of temperatures considered. There are small differences between classical and quantum results of water at 298 K, necessitating a small correction term to reproduce water's enthalpy of vaporisation. Only at higher temperatures is entropy underestimated by up to 9 J K(-1) mol(-1) as verified by thermodynamic integration calculations. Satisfactory agreement for ice, however, is only obtained by using the quantum formulation. Even then, at higher temperatures, the entropies exceed experiment by up to 15 J K(-1) mol(-1). Further insight into the quantum nature of water is provided by inspecting the temperature dependence of the frequencies. The harmonic approximation is further supported by the harmonic force and torque distributions and the very similar entropies obtained from the force and torque variances. All these results suggest that the single molecule harmonic oscillator approximation for water, although not exact, provides a rapid, insightful, and useful means to evaluate the thermodynamic properties of water from a computer simulation in a way that can account for the quantization of water's energy levels.     

Thermodynamic anomalies in a lattice model of water

We investigate a lattice-fluid model of water, defined on a three-dimensional body centered cubic lattice. Model molecules possess a tetrahedral symmetry, with four equivalent bonding arms, aiming to mimic the formation of hydrogen bonds. The model is similar to the one proposed by Roberts and Debenedetti [J. Chem. Phys. 105, 658 (1996)], simplified in that no distinction between bond "donors" and "acceptors" is imposed. Bond formation depends both on orientation and local density. In the ground state, we show that two different ordered (ice) phases are allowed. At finite temperature, we analyze homogeneous phases only, working out phase diagram, response functions, the temperature of maximum density locus, and the Kauzmann line. We make use of a generalized first-order approximation on a tetrahedral cluster. In the liquid phase, the model exhibits several anomalous properties observed in real water. In the low temperature region (supercooled liquid), there are evidences of a second critical point and, for some range of parameter values, this scenario is compatible with the existence of a reentrant spinodal.     

Remarks on “A new equation for modelling nonisothermal reactions”

The formula proposed by Agrawal [7] for the approximation of the temperature integral in non-isothermal kinetics is shown to be less accurate than several approximations proposed earlier that are of the same complexity. The domain of applicability of 10 approximate formulae is discussed.     

Decoupling charge transport from the structural dynamics in room temperature ionic liquids

Light scattering and dielectric spectroscopy measurements were performed on the room temperature ionic liquid (RTIL) [C4mim][NTf2] in a broad temperature and frequency range. Ionic conductivity was used to estimate self-diffusion of ions, while light scattering was used to study structural relaxation. We demonstrate that the ionic diffusion decouples from the structural relaxation process as the temperature of the sample decreases toward T(g). The strength of the decoupling appears to be significantly lower than that expected for a supercooled liquid of similar fragility. The structural relaxation process in the RTIL follows well the high-temperature mode coupling theory (MCT) scenario. Using the MCT analysis we estimated the dynamic crossover temperature in [C4mim][NTf2] to be T(c) ~ 225 ± 5 K. However, our analysis reveals no sign of the dynamic crossover in the ionic diffusion process.     

Dynamics of low temperature induced water shedding from AOT reverse micelles

The effects of low temperature and ionic strength on water encapsulated within reverse micelles were investigated by solution NMR. Reverse micelles composed of AOT and pentane and solutions with varying concentrations of NaCl were studied at temperatures ranging from 20 degrees C to -30 degrees C. One-dimensional (1)H solution NMR spectroscopy was used to monitor the quantity and structure of encapsulated water. At low temperatures, e.g., -30 degrees C, reverse micelles lose water at rates that are dependent on the ionic strength of the aqueous nanopool. The final water loading (w0 = [water]/[surfactant]) of the reverse micelles is likewise dependent on the ionic strength of the aqueous phase. Remarkably, water resonance(s) at temperatures between -20 degrees C and -30 degrees C displayed fine structure indicating the presence of multiple transient water populations. Results of this study demonstrate that reverse micelles are an excellent vehicle for studies of confined water across a broad range of conditions, including the temperature range that provides access to the supercooled state.     

Revealing long-range density fluctuations in dialkylimidazolium chloride ionic liquids by dynamic light scattering

In this study, the structures and dynamics of ionic liquids of 1-allyl-3-methylimidazolium chloride ([AMIM][Cl]) and 1-n-butyl-3-methylimidazolium chloride ([BMIM][Cl]) were studied by dynamic light scattering with polarized and depolarized geometries in the temperature range from 300 to 400 K. The temperature range covered supercooled and liquid states for [BMIM][Cl] and covered the liquid state for [AMIM][Cl]. The results show that for these ionic liquids at all chosen temperatures only one ultraslow relaxation is observed in the polarized component of dynamic light scattering, however, the ultraslow relaxation is not observed in the depolarized component. The ultraslow relaxation exhibited several typical features of the "cluster" mode generally found in glass-forming liquids and polymer melts, such as diffusive, strongly scattering-vector-dependent, and nearly exponential characters, which thus corresponded to long-range density fluctuations. The physical origin for long-range density fluctuations was the existence of heterogeneities with large characteristic length scales in these ionic liquids. It was further considered that molecules of these ionic liquids not only tended to aggregate to form dynamic clusters but also possibly formed dynamic networks in the supercooled state and the heterogeneities could exist even at temperatures higher than the melting points.     

Comparison of physicochemical properties of new ionic liquids based on imidazolium, quaternary ammonium, and guanidinium cations

More than 50 ionic liquids were prepared by using imidazolium, quaternary ammonium, and guanidinium cations and various anions. In these series, different cationic structures such as 1-benzyl-3-methylimidazolium [Bzmim]+, 1,3-dibenzylimidazolium [BzmiBz]+, 1-octyl-3-methylimidazolium [C8mim]+, 1-decyl-3-methylimidazolium [C10mim]+, tricapryl-methylammonium [Aliquat]+, benzyltriethylammonium [BzTEA]+, phenyltrimethylammonium [PhTMA]+, and dimethyldihexylguanidinium [DMG]+ were combined with anions, p-toluenesulfonate [TSA](-), dicyanoamide [DCA]-, saccharine (2-sulfobenzoic acid imide sodium salt) [SAC]-, trifluoroacetate [TFA]-, bis(trifluoromethanesulfonyl)imide [Tf2N]-, trifluoromethanesulfonate [TfO]-, and thiocyanate [SCN]-. Important physical data for these ionic liquids are collated, namely solubility in common solvents, viscosity, density, melting point and water content. Apart from the viscosity, the Newtonian and non-Newtonian behavior of these ionic liquids is also disclosed. Stability of these ionic liquids under thermal, basic, acidic, nucleophilic, and oxidative conditions was also studied. The features of the solid-liquid phase transition were analyzed, namely the glass transition temperature and the heat capacity jump associated with the transition from the non-equilibrium glass to the metastable supercooled liquid. A degradation temperature of each ionic liquid was also determined. Comparisons of the properties of various ionic liquids were made.     

A non-polarizable model of water that yields the dielectric constant and the density anomalies of the liquid: TIP4Q

A four-site rigid water model is presented, whose parameters are fitted to reproduce the experimental static dielectric constant at 298 K, the maximum density of liquid water and the equation of state at low pressures. The model has a positive charge on each of the three atomic nuclei and a negative charge located at the bisector of the HOH bending angle. This charge distribution allows increasing the molecular dipole moment relative to four-site models with only three charges and improves the liquid dielectric constant at different temperatures. Several other properties of the liquid and of ice Ih resulting from numerical simulations with the model are in good agreement with experimental values over a wide range of temperatures and pressures. Moreover, the model yields the minimum density of supercooled water at 190 K and the minimum thermal compressibility at 310 K, close to the experimental values. A discussion is presented on the structural changes of liquid water in the supercooled region where the derivative of density with respect to temperature is a maximum.     

The enthalpies and heat capacities of hydration of ammonium and tetraalkylammonium bromides

The heat effects of solution of ammonium, tetraethylammonium, and tetrabutylammonium bromides in water were measured over the temperature range 277.15-328.15 K. The standard enthalpies of solution were calculated in the second approximation of the Debye-Hückel theory. The results were compared with the literature data on the enthalpies of solution of tetrabutylammonium and tetrapentylammonium bromides. The applicability of various equations suggested to describe the temperature dependences of enthalpy characteristics was analyzed. Heat capacity changes caused by the hydration of the salts and the temperature dependences of the entropies of hydration were calculated and discussed.     

Thermodynamic anomalies in a lattice model of water: solvation properties

We investigate a lattice-fluid model of water, defined on a three-dimensional body-centered-cubic lattice. Model molecules possess a tetrahedral symmetry, with four equivalent bonding arms. The model is similar to the one proposed by Roberts and Debenedetti [J. Chem. Phys. 105, 658 (1996)], simplified by removing distinction between "donors" and "acceptors." We focus on the solvation properties, mainly as far as an ideally inert (hydrophobic) solute is concerned. As in our previous analysis, devoted to neat water [J. Chem. Phys. 121, 11856 (2004)], we make use of a generalized first-order approximation on a tetrahedral cluster. We show that the model exhibits quite a coherent picture of water thermodynamics, reproducing qualitatively several anomalous properties observed both in pure water and in solutions of hydrophobic solutes. As far as supercooled liquid water is concerned, the model is consistent with the second critical-point scenario.     

Measurement and analysis of the diffuse reflectance of powdered samples at terahertz frequencies using a quantum cascade laser

We report terahertz (THz) diffuse reflectance measurements of bulk powdered samples at a frequency of 2.83 THz using a narrowband quantum cascade laser. Samples studied comprise polydisperse powders with absorption coefficients extending over two orders of magnitude from ～3 cm(-1) to >200 cm(-1). Diffuse reflectance measurements are used to obtain the effective absorption coefficient of these samples from the backscattering cross-section, predicted under the quasi-crystalline approximation (QCA) in the T-matrix formulation and in conjunction with the Percus-Yevick pair distribution function. Results are compared with effective absorption coefficients obtained from THz time-domain spectroscopy measurements on pressed pellet samples, and show good agreement over the range of effective absorption coefficients studied. We observe that the backscattering cross-section predicted under the QCA is strongly dependent on both the real and imaginary components of the complex permittivity of the sample, and we show that reliable determination of the absorption coefficient from diffuse reflectance measurements therefore requires knowledge of the sample's refractive index. This work demonstrates the applicability of diffuse reflectance measurements, using a THz frequency quantum cascade laser, to the high-resolution spectroscopic analysis of bulk powdered samples at THz frequencies.     

Dielectric and conductivity relaxations in poly(hema) and of water in its hydrogel

The dielectric permittivity ε′ and loss ε″ of anhydrous poly(2-hydroxyethyl methacrylate) and its 38.6 w/w% hydrogel have been measured in the frequency range from 12 Hz to 200 kHz and the temperature range from 77 to 273 K. The former has a sub-T_g relaxation with a half-width of 4.5 decades for the loss spectra, whose strength increases with temperature, and an activation energy of 62.5 kJ/mol. The dielectric relaxation time of the α process of supercooled water in the hydrogel is 53 s at its calorimetric T_g of 135 K. The half-width of the relaxation spectrum is 2.85 decades and, in the narrow temperature range, its apparent activation energy is 60.8 kJ/mol. Heating of the hydrogel causes crystallization of water which begins at about 207 K and becomes readily detectable as a second dielectric loss peak at about 230 K. For each temperature between 207 and 267 K, supercooled water in the hydrogel coexists with its crystallized form, with the amount of the crystallized solid increasing with increasing temperature. These results are discussed in terms of “bound” and “free” states of water in the hydrogel.     

Correlation between quasielastic Raman scattering and configurational entropy in an ionic liquid

Low-frequency (5-200 cm(-1)) Raman spectra are reported for the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, [bmim]PF(6), in glassy, supercooled liquid, and normal liquid phases (77-330 K). Raman spectra of [bmim]PF(6) agree with previous results obtained by optical Kerr effect spectroscopy and molecular dynamics simulation. Both the superposition model and the coupling model give reasonable fit to low-frequency Raman spectra of [bmim]PF(6). The configurational entropy of [bmim]PF(6) has been evaluated as a function of temperature using recently reported data of heat capacity. The calculated configurational entropy is inserted in the Adam-Gibbs theory for supercooled liquids, giving a good fit to non-Arrhenius behavior of viscosity and diffusive process, with the latter revealed by a recent neutron scattering investigation of [bmim]PF(6). There is a remarkable linear dependence between intensity of quasielastic Raman scattering and configurational entropy from 77 K up to the melting point of [bmim]PF(6). This correlation offers insight into the nature of dynamical processes probed by low-frequency Raman spectra of ionic liquids.     

Interaction of water with LiCl, LiBr, and LiI in the deeply supercooled region

The hydration mechanism of lithium halides was studied using time-of-flight secondary ion mass spectrometry as a function of temperature. The lithium halides embedded in thin films of amorphous solid water segregate to the surface at temperatures higher than 135-140 K, with efficiency increasing in the order of LiCl, LiBr, and LiI. A monolayer of LiCl and LiI adsorbed on the surface of amorphous solid water tends to diffuse into the bulk at 160 K. The infrared absorption band revealed that the aqueous lithium-halide solutions and crystals are formed simultaneously at 160 K; these phenomena are explicable as a consequence of the evolution of supercooled liquid water. The strong surfactant effect is inferred to arise from hydration of a contact ion pair having hydrophilic (lithium) and hydrophobic (halide) moieties. Furthermore, bulk diffusion of lithium halides might result from the formation of a solvent-separated ion pair in supercooled liquid water. The presence of two liquid phases of water with different local structures is probably responsible for the formation of these two hydrates, consistent with the calculated result reported by Jungwirth and Tobias[J. Phys. Chem. B 106, 6361 (2002)].     

Breaking of the supercooled state of water by a nanocavity with disordered atomic configuration I: freezing behavior of sorbed water into polymethylmethacrylate film as examined by Fourier transform infrared spectroscopy

A freezing behavior of sorbed water into an annealed atactic polymethylmethacrylate (PMMA) film was examined by Fourier transform infrared (IR) spectroscopy. Temperature dependence of IR spectra of the sorbed water indicated that nonfreezable water in the PMMA film formed associated waters on cooling and were not universally nonfreezable. Its freezing temperature was found to be 249 K, which is much higher than the homogeneous nucleation temperature of water, 236 K. This result suggests that a disordered atomic configuration of a nanocavity wall with hydrogen bond acceptors induces a breaking of the supercooled state of water. The former two findings were visibly clarified by the vibrational spectroscopic method for the first time, and the latest finding was different from the literature on supercooled water in small spaces.     

The influence of temperature on the phase behavior of ionic liquid aqueous two-phase systems

In this article, the equilibrium behavior of solutions of 1-ethyl-3-methylimidazolium dimethyl phosphate ([Emim]DMP) and ethyl acetate or acetone in aqueous two-phase system (ATPS) was discussed to understand the liquid–liquid equilibrium (LLE) behavior of these organic solvents. Thus, we determined phase diagrams and LLE data at 303.15, 308.15, 313.15, and 323.15 K for the investigated biphasic systems. Four empirical equations were used to study the tie lines. The results showed that for the [Emim]DMP?+?acetone?+?water biphasic systems within the investigated temperature range, temperature influences the phase behavior, but for the [Emim]DMP?+?ethyl acetate?+?water biphasic systems within the investigated temperature range, there are no notable changes on the phase behavior with rising temperature. The results may have important applications for the separation of antibiotics and for the recovery of ionic liquids (ILs).     

Application of the transfer-matrix method to equilibria of polymer reactions and adsorption from solution

The transfer-matrix method is used to describe the occupation of a lattice of sites in equilibrium with bulk solution, thus giving a model for adsorption from solution (or binding or reaction). Numerical illustration is given using the zeroth approximation of the quasi-crystalline model for the bulk phase. Application of the method to two-dimensional systems involving orientation effects and phase changes is briefly discussed.     

Liquid-vapor and liquid-liquid phase equilibria of the Brodholt-Sampoli-Vallauri polarizable water model

Liquid-vapor and liquid-liquid phase equilibria of the polarizable Brodholt-Sampoli-Vallauri water model have been investigated by Gibbs ensemble Monte Carlo computer simulations. The coexisting liquid and vapor densities and energy of vaporization of the model is found to be in a reasonable agreement with experimental data in the entire temperature range of liquid-vapor coexistence. The critical temperature and density of the model are found to be 615 K and 0.278 gcm(3), respectively, close to the experimental values of 647.1 K and 0.322 gcm(3). In the supercooled state two distinct liquid-liquid coexistence regions are observed. The existence of liquid-liquid phase separation of a polarizable water model is demonstrated for the first time.     

Heat capacity of tetrahydrofuran clathrate hydrate and of its components, and the clathrate formation from supercooled melt

We report a thermodynamic study of the formation of tetrahydrofuran clathrate hydrate by explosive crystallization of water-deficient, near stoichiometric, and water-rich solutions, as well as of the heat capacity, C(p), of (i) supercooled tetrahydrofuran-H2O solutions and of the clathrate hydrate, (ii) tetrathydrofuran (THF) liquid, and (iii) supercooled water and the ice formed on its explosive crystallization. In explosive freezing of supercooled solutions at a temperature below 257 K, THF clathrate hydrate formed first. The nucleation temperature depends on the cooling rate, and excess water freezes on further cooling. The clathrate hydrate melts reversibly at 277 K and C(p) increases by 770 J/mol K on melting. The enthalpy of melting is 99.5 kJ/mol and entropy is 358 J/mol K. Molar C(p) of the empty host lattice is less than that of the ice, which is inconsistent with the known lower phonon frequency of H2O in the clathrate lattice. Analysis shows that C(p) of THF and ice are not additive in the clathrate. C(p) of the supercooled THF-H2O solutions is the same as that of water at 247 K, but less at lower temperatures and more at higher temperatures. The difference tends to become constant at 283 K. The results are discussed in terms of the hydrogen-bonding changes between THF and H2O.