Thermal conductivity of supercooled water

Thermal Transport and Thermal Polarization of Water in the Supercooled Regime

The behavior of water in the deep supercooled regime has attracted significant interest, motivated by the hypothesis of the second critical point of water. Previous studies indicated the existence of a water anomaly, characterized by a minimum in the thermal conductivity of water. Here, we employ nonequilibrium molecular dynamics computer simulation and the TIP4P/2005 water force field to investigate the thermal conductivity of supercooled water targeting four different isobars, 1, 200, 700, and 1200 bar. We demonstrate using NEMD simulations the existence of minima in thermal conductivity associated with the maximum isothermal compressibility and the minimum speed of sound in water, hence establishing a firm connection with the second critical point of liquid water. Moreover, we demonstrate that thermal gradients polarize supercooled water with a thermal polarization coefficient of several mV/K. We explain the thermal polarization effect using a theoretical formulation introduced recently that connects the thermal polarization effect to the isobaric thermal expansion coefficient.

Thermal Conductivity of Ionic Liquids and IoNanofluids. Can Molecular Theory Help?

Ionic liquids have been suggested as new engineering fluids, specifically in the area of heat transfer, and as alternatives to current biphenyl and diphenyl oxide, alkylated aromatics and dimethyl polysiloxane oils, which degrade above 200 °C, posing some environmental problems. Addition of nanoparticles to produce stable dispersions/gels of ionic liquids has proved to increase the thermal conductivity of the base ionic liquid, potentially contributing to better efficiency of heat transfer fluids. It is the purpose of this paper to analyze the prediction and estimation of the thermal conductivity of ionic liquids and IoNanofluids as a function of temperature, using the molecular theory of Bridgman and estimation methods previously developed for the base fluid. In addition, we consider methods that emphasize the importance of the interfacial area IL-NM in modelling the thermal conductivity enhancement. Results obtained show that it is not currently possible to predict or estimate the thermal conductivity of ionic liquids with an uncertainty commensurate with the best experimental values. The models of Maxwell and Hamilton are not capable of estimating the thermal conductivity enhancement of IoNanofluids, and it is clear that the Murshed, Leong and Yang model is not practical, if no additional information, either using imaging techniques at nanoscale or molecular dynamics simulations, is available.

Thermal Conductivity of Metastable Ionic Liquid [C2mim][CH3SO3]

Ionic liquids have been suggested as new engineering fluids, namely in the area of heat transfer, as alternatives to current biphenyl and diphenyl oxide, alkylated aromatics and dimethyl polysiloxane oils, which degrade above 200 °C and pose some environmental problems. Recently, we have proposed 1-ethyl-3-methylimidazolium methanesulfonate, [C₂mim][CH₃SO₃], as a new heat transfer fluid, because of its thermophysical and toxicological properties. However, there are some interesting points raised in this work, namely the possibility of the existence of liquid metastability below the melting point (303 K) or second order-disorder transitions (λ-type) before reaching the calorimetric freezing point. This paper analyses in more detail this zone of the phase diagram of the pure fluid, by reporting accurate thermal-conductivity measurements between 278 and 355 K with an estimated uncertainty of 2% at a 95% confidence level. A new value of the melting temperature is also reported, T_melt = 307.8 ± 1 K. Results obtained support liquid metastability behaviour in the solid-phase region and permit the use of this ionic liquid at a heat transfer fluid at temperatures below its melting point. Thermal conductivity models based on Bridgman theory and estimation formulas were also used in this work, failing to predict the experimental data within its uncertainty.

Liquid-liquid critical point in a simple analytical model of water

A statistical model for a simple three-dimensional Mercedes-Benz model of water was used to study phase diagrams. This model on a simple level describes the thermal and volumetric properties of waterlike molecules. A molecule is presented as a soft sphere with four directions in which hydrogen bonds can be formed. Two neighboring waters can interact through a van der Waals interaction or an orientation-dependent hydrogen-bonding interaction. For pure water, we explored properties such as molar volume, density, heat capacity, thermal expansion coefficient, and isothermal compressibility and found that the volumetric and thermal properties follow the same trends with temperature as in real water and are in good general agreement with Monte Carlo simulations. The model exhibits also two critical points for liquid-gas transition and transition between low-density and high-density fluid. Coexistence curves and a Widom line for the maximum and minimum in thermal expansion coefficient divides the phase space of the model into three parts: in one part we have gas region, in the second a high-density liquid, and the third region contains low-density liquid.

Thermodynamic Analysis of the Two-Liquid Model for Anomalies of Water, HDL-LDL Fluctuations, and Liquid-Liquid Transition

After reviewing the protocol-dependent properties of HDA, which thermally anneals to LDA, and the data gap over an unusually large T-range between HDA, LDA, and water, we investigate whether or not, despite HDA's ill-defined state and distinction from a glass, the HDL-LDL fluctuations view of the two-liquid model can explain water's anomalous behavior. An analysis of the density, ρ, compressibility, β, heat capacity, Cp, and thermal conductivity, κ, of water over a monotonic (continuous) path bridging this data gap shows the following: (i) Such a path between ρwater at 320 K and ρHDA yields an untenable thermal expansion coefficient of water. (ii) There is neither a continuous path between βwater at 353 K and βHDA, nor between Cp,water at 363 K and Cp,HDA. (iii) The same value of ρwater, of βwater, or of Cp,water at two temperatures separated by a maxima or a minima is incompatible with the HDL-LDL fluctuations view. (iv) κLDA at ∼140 K is about twice that of κ water at 253 K. (v) κHDA at 120 K is incompatible with κwater at T > 320 K. Thus, there is an internal inconsistency between the thermodynamics of HDA seen as a glass and that of water seen as an HDL-LDL mixture, which is incompatible with both the HDL-LDL fluctuations view and the liquid-liquid transition.

Existence of a liquid-liquid phase transition in methanol

A simple model is constructed to study the phase diagram and thermodynamic properties of methanol, which is described as a dimer of an apolar sphere mimicking the methyl group and a sphere with core-softened potential as the hydroxyl group. Performing classical Monte Carlo simulations, we obtained the phase diagram, showing a second critical point between two different liquid phases. Evaluating systems with a different number of particles, we extrapolate to infinite size in accordance with Ising universality class to obtain bulk values for critical temperature, pressure, and density. Strong evidence that the structure of the liquid changes upon transition from high- to low-density phase was provided. From the experimentally determined hydrogen bond strength and length in methanol and water, we propose where the second critical point of methanol should be.

Two-state thermodynamics of the ST2 model for supercooled water

Thermodynamic properties of the ST2 model for supercooled liquid water exhibit anomalies similar to those observed in real water. A possible explanation of these anomalies is the existence of a metastable, liquid-liquid transition terminated by a critical point. This phenomenon, whose possible existence in real water is the subject of much current experimental work, has been unambiguously demonstrated for this particular model by most recent simulations. In this work, we reproduce the anomalies of two versions of the ST2 model with an equation of state describing water as a non-ideal "mixture" of two different types of local molecular order. We show that the liquid-liquid transition in the ST2 water is energy-driven. This is in contrast to another popular model, mW, in which non-ideality in mixing of two alternative local molecular orders is entropy-driven, and is not sufficiently strong to induce a liquid-liquid transition.

Thermal Conductivity of Supercooled Water: An Equilibrium Molecular Dynamics Exploration

The thermal conductivity of both supercooled and ambient-temperature water at atmospheric pressure has been computed over the 140-270 K temperature range for three popular water models via equilibrium molecular dynamics in the Green-Kubo setting. No strong temperature dependence of thermal conductivity was observed. The underlying phonon modes contributing to thermal conduction processes have been examined in the present work, and it has been established that (translational) acoustic modes dominate in supercooled water.

Thermodynamics and the hydrophobic effect in a core-softened model and comparison with experiments

A simple and computationally inexpensive core-softened model, originally proposed by Franzese [G. Franzese, J. Mol. Liq. 136, 267 (2007)], was adopted to show that it exhibits properties of waterlike fluid and hydrophobic effect. The potential used between particles is spherically symmetric with two characteristic lengths. Thermodynamics of nonpolar solvation were modeled as an insertion of a modified Lennard-Jones particle. It was investigated how the anomalous predictions of the model as well as the nonpolar solvation compare with the experimental data for water anomalies and the temperature dependence of noble gases hydration. It was shown that the model qualitatively follows the same trends as water. The model is able to reproduce waterlike anomalous properties (density maximum, heat capacity minimum, isothermal compressibility, etc.) and hydrophobic effect (minimum solubility for nonpolar solutes near ambient conditions, increased solubility of larger noble gases, etc.). It is argued that the model yields similar results as more complex and computationally expensive models.

The hydrophobic effect in a simple isotropic water-like model: Monte Carlo study

Using Monte Carlo computer simulations, we show that a simple isotropic water-like model with two characteristic lengths can reproduce the hydrophobic effect and the solvation properties of small and large non-polar solutes. Influence of temperature, pressure, and solute size on the thermodynamic properties of apolar solute solvation in a water model was systematically studied, showing two different solvation regimes. Small particles can fit into the cavities around the solvent particles, inducing additional order in the system and lowering the overall entropy. Large particles force the solvent to disrupt their network, increasing the entropy of the system. At low temperatures, the ordering effect of small solutes is very pronounced. Above the cross-over temperature, which strongly depends on the solute size, the entropy change becomes strictly positive. Pressure dependence was also investigated, showing a "cross-over pressure" where the entropy and enthalpy of solvation are the lowest. These results suggest two fundamentally different solvation mechanisms, as observed experimentally in water and computationally in various water-like models.

Core-softened fluids as a model for water and the hydrophobic effect

An interaction model with core-softened potential in three dimensions was studied by Monte Carlo computer simulations and integral equation theory. We investigated the possibility that a fluid with a core-softened potential can reproduce anomalies found experimentally in liquid water, such as the density anomaly, the minimum in the isothermal compressibility as a function of temperature, and others. Critical points of the fluid were also determined. We provided additional arguments that the old notion, postulating that only angular-dependent interactions result in density anomaly, is incorrect. We showed that potential with two characteristic distances is sufficient for the system to exhibit water-like behavior and anomalies, including the famous density maximum. We also found that this model can properly describe the hydrophobic effect.

Possible relation of water structural relaxation to water anomalies

The anomalous behavior of thermodynamic response functions is an unsolved problem in the physics of water. The mechanism that gives rise to the dramatic indefinite increase at low temperature in the heat capacity, the compressibility, and the coefficient of thermal expansion, is unknown. We explore this problem by analyzing both new and existing experimental data on the power spectrum S(Q, ω) of bulk and confined water at ambient pressure. When decreasing the temperature, we find that the liquid undergoes a structural transformation coinciding with the onset of an extended hydrogen bond network. This network onset seems to give rise to the marked viscoelastic behavior, consistent with the interesting possibility that the sound velocity and response functions of water depend upon both the frequency and wave vector.