Thermal conductivity of ice polymorphs: a computational study.

Thermal transport in ice features an unusual response. In addition to its intrinsic scientific interest, an understanding of the mechanisms determining the thermal conductivity of ice might be relevant in climate modelling and planetary science. Accurate microscopic models can provide important molecular insight into these mechanisms. In this work, we quantify using molecular simulations and state of the art forcefields, the thermal conductivity of ice Ih, VI, VII and a plastic phase that has been proposed very recently at pressures in the GPa range. The TIP4P models used in this study underestimate significantly the thermal conductivity of ice Ih and ice VII, while they show good agreement with experimental measurements of ice VI. The discrepancies observed are examined by investigating the temperature dependence of the thermal conductivity. The simulations indicate that the models are too anharmonic and they potentially feature a higher structural disorder than the experimental systems. We suggest that at high pressures the simulated thermal conductivities can be rationalized in terms of the performance of the models in predicting the equation of state of ice. The thermal conductivity of the plastic phase is very similar to that of the coexisting ice VII. Since the water molecules in the plastic phase feature orientational disorder, these results indicate that the hydrogen bond network does not play a significant role in defining the thermal transport mechanisms of ice at high pressures.

The rich phase behavior of the thermopolarization of water: from a reversal in the polarization, to enhancement near criticality conditions.

We investigate using non-equilibrium molecular dynamics simulations the polarization of water induced by thermal gradients using the accurate TIP4P/2005 water model. The full dependence of the polarization covering a wide range of thermodynamic states, from near supercritical to ambient conditions, is reported. Our results show a strong dependence of the thermo-polarization field with the thermodynamic state. The field features a strong enhancement near the critical point, which can be rationalized in terms of the large increase and ultimately the divergence of the thermal expansion of the fluid at the critical temperature. We also show that the TIP4P/2005 model features a reversal in the sign of the thermal polarization at densities ∼1 g cm(-3). The latter result is consistent with the recent observation of this reversal phenomenon in SPC/E water and points the existence of this general physical phenomenon in water.

Mechanical behavior and interphase structure in a silica-polystyrene nanocomposite under uniaxial deformation.

The mechanical behavior of polystyrene and a silica-polystyrene nanocomposite under uniaxial elongation has been studied using a coarse-grained molecular dynamics technique. The Young's modulus, the Poisson ratio and the stress-strain curve of polystyrene have been computed for a range of temperatures, below and above the glass transition temperature. The predicted temperature dependence of the Young's modulus of polystyrene is compared to experimental data and predictions from atomistic simulations. The observed mechanical behavior of the nanocomposite is related to the local structure of the polymer matrix around the nanoparticles. Local segmental orientational and structural parameters of the deforming matrix have been calculated as a function of distance from nanoparticle's surface. A thorough analysis of these parameters reveals that the segments close to the silica nanoparticle's surface are stiffer than those in the bulk. The thickness of the nanoparticle-matrix interphase layer is estimated. The Young's modulus of the nanocomposite has been obtained for several nanoparticle volume fractions. The addition of nanoparticles results in an enhanced Young's modulus. A linear relation describes adequately the dependence of Young's modulus on the nanoparticle volume fraction.