Equilibrium molecular dynamics simulations have been carried out in the microcanonical ensemble at 300 and 255 K on the extended simple point charge (SPC/E) model of water [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)]. In addition to a number of static and dynamic properties, thermal conductivity lambda has been calculated via Green-Kubo integration of the heat current time correlation functions (CF's) in the atomic and molecular formalism, at wave number k=0. The calculated values (0.67 +/- 0.04 W/mK at 300 K and 0.52 +/- 0.03 W/mK at 255 K) are in good agreement with the experimental data (0.61 W/mK at 300 K and 0.49 W/mK at 255 K). A negative long-time tail of the heat current CF, more apparent at 255 K, is responsible for the anomalous decrease of lambda with temperature. An analysis of the dynamical modes contributing to lambda has shown that its value is due to two low-frequency exponential-like modes, a faster collisional mode, with positive contribution, and a slower one, which determines the negative long-time tail. A comparison of the molecular and atomic spectra of the heat current CF has suggested that higher-frequency modes should not contribute to lambda in this temperature range. Generalized thermal diffusivity D-T(k) decreases as a function of k, after an initial minor increase at k = k(min). The k dependence of the generalized thermodynamic properties has been calculated in the atomic and molecular formalisms. The observed differences have been traced back to intramolecular or intermolecular rotational effects and related to the partial structure functions. Finally, from the results we calculated it appears that the SPC/E model gives results in better agreement with experimental data than the transferable intermolecular potential with four points TIP4P water model [Jorgensen et al., J. Chem. Phys. 79, 926 (1983)], with a larger improvement for, e.g., diffusion, viscosities, and dielectric properties and a smaller one for thermal conductivity. The SPC/E model shares, to a smaller extent, the insufficient slowing down of dynamics at low temperature already found for the TIP4P water model.
Thermal conductivity of water: Molecular dynamics and generalized hydrodynamics results
TANI, ALESSANDRO
1997-01-01
Abstract
Equilibrium molecular dynamics simulations have been carried out in the microcanonical ensemble at 300 and 255 K on the extended simple point charge (SPC/E) model of water [Berendsen et al., J. Phys. Chem. 91, 6269 (1987)]. In addition to a number of static and dynamic properties, thermal conductivity lambda has been calculated via Green-Kubo integration of the heat current time correlation functions (CF's) in the atomic and molecular formalism, at wave number k=0. The calculated values (0.67 +/- 0.04 W/mK at 300 K and 0.52 +/- 0.03 W/mK at 255 K) are in good agreement with the experimental data (0.61 W/mK at 300 K and 0.49 W/mK at 255 K). A negative long-time tail of the heat current CF, more apparent at 255 K, is responsible for the anomalous decrease of lambda with temperature. An analysis of the dynamical modes contributing to lambda has shown that its value is due to two low-frequency exponential-like modes, a faster collisional mode, with positive contribution, and a slower one, which determines the negative long-time tail. A comparison of the molecular and atomic spectra of the heat current CF has suggested that higher-frequency modes should not contribute to lambda in this temperature range. Generalized thermal diffusivity D-T(k) decreases as a function of k, after an initial minor increase at k = k(min). The k dependence of the generalized thermodynamic properties has been calculated in the atomic and molecular formalisms. The observed differences have been traced back to intramolecular or intermolecular rotational effects and related to the partial structure functions. Finally, from the results we calculated it appears that the SPC/E model gives results in better agreement with experimental data than the transferable intermolecular potential with four points TIP4P water model [Jorgensen et al., J. Chem. Phys. 79, 926 (1983)], with a larger improvement for, e.g., diffusion, viscosities, and dielectric properties and a smaller one for thermal conductivity. The SPC/E model shares, to a smaller extent, the insufficient slowing down of dynamics at low temperature already found for the TIP4P water model.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
