Nonlocal Dielectric Response of Water in Nanoconfinement.

Recent experiments reporting a very low dielectric permittivity for nanoconfined water have renewed the interest in the structure and dielectric properties of water in narrow gaps. Here, we describe such systems with a minimal Landau-Ginzburg field theory composed of a nonlocal bulk-determined term and a local water-surface interaction term. We show how the interplay between the boundary conditions and intrinsic bulk correlations encodes the dielectric properties of confined water. Our theoretical analysis is supported by molecular dynamics simulations and comparison with the experimental data.

Electrotunable lubricity with ionic liquid nanoscale films.

One of the main challenges in tribology is finding the way for an in situ control of friction without changing the lubricant. One of the ways for such control is via the application of electric fields. In this respect a promising new class of lubricants is ionic liquids, which are solvent-free electrolytes, and their properties should be most strongly affected by applied voltage. Based on a minimal physical model, our study elucidates the connection between the voltage effect on the structure of the ionic liquid layers and their lubricating properties. It reveals two mechanisms of variation of the friction force with the surface charge density, consistent with recent AFM measurements, namely via the (i) charge effect on normal and in-plane ordering in the film and (ii) swapping between anion and cation layers at the surfaces. We formulate conditions that would warrant low friction coefficients and prevent wear by resisting "squeezing-out" of the liquid under compression. These results give a background for controllable variation of friction.

Communication: minimum in the thermal conductivity of supercooled water: a computer simulation study.

We report the results of a computer simulation study of the thermodynamic properties and the thermal conductivity of supercooled water as a function of pressure and temperature using the TIP4P-2005 water model. The thermodynamic properties can be represented by a two-structure equation of state consistent with the presence of a liquid-liquid critical point in the supercooled region. Our simulations confirm the presence of a minimum in the thermal conductivity, not only at atmospheric pressure, as previously found for the TIP5P water model, but also at elevated pressures. This anomalous behavior of the thermal conductivity of supercooled water appears to be related to the maximum of the isothermal compressibility or the minimum of the speed of sound. However, the magnitudes of the simulated thermal conductivities are sensitive to the water model adopted and appear to be significantly larger than the experimental thermal conductivities of real water at low temperatures.

Note: local thermal conductivities from boundary driven non-equilibrium molecular dynamics simulations.

We report non-equilibrium molecular dynamics simulations of heat transport in models of molecular fluids. We show that the "local" thermal conductivities obtained from non-equilibrium molecular dynamics simulations agree within numerical accuracy with equilibrium Green-Kubo computations. Our results support the local equilibrium hypothesis for transport properties. We show how to use the local dependence of the thermal gradients to quantify the thermal conductivity of molecular fluids for a wide range of thermodynamic states using a single simulation.

Water polarization induced by thermal gradients: the extended simple point charge model (SPC/E).

We investigate the non-equilibrium response of extended simple point charge (SPC/E) water to thermal gradients. Using non-equilibrium molecular dynamics simulations, we show that SPC/E water features the thermo-polarization orientation effect, namely, water becomes polarized as a response to a thermal gradient. The polarization field increases linearly with the thermal gradient, in agreement with predictions of non-equilibrium thermodynamics theory. This observation confirms the generality of the thermo-polarization effect, first reported using the Modified Central Force Model (MCFM), and shows this physical effect is present irrespective of the water model details, in particular, dipole moment magnitude and model flexibility. The magnitude of the effect is the same for both models, although the sign of the electrostatic field is reversed in going from the MCFM to the SPC/E model. We further analyze the impact that the molecular geometry and mass distribution has on the magnitude of the polarization. Our results indicate that the thermo-polarization effect should be observed in a wide range of polar fluids, including fluids where hydrogen bonding is not present. Using various molecular models, we show that the polarization of these fluids under appropriate thermodynamic conditions can be of the same order or stronger than in water.

Molecular dynamics computations of brine-CO2 interfacial tensions and brine-CO2-quartz contact angles and their effects on structural and residual trapping mechanisms in carbon geo-sequestration.

In the context of carbon geo-sequestration projects, brine-CO(2) interfacial tension γ and brine-CO(2)-rock surface water contact angles θ directly impact structural and residual trapping capacities. While γ is fairly well understood there is still large uncertainty associated with θ. We present here an investigation of γ and θ using a molecular approach based on molecular dynamics computer simulations. We consider a system consisting of CO(2)/water/NaCl and an α-quartz surface, covering a brine salinity range between 0 and 4 molal. The simulation models accurately reproduce the dependence of γ on pressure below the CO(2) saturation pressure at 300 K, and over predict γ by ~20% at higher pressures. In addition, in agreement with experimental observations, the simulations predict that γ increases slightly with temperature or salinity. We also demonstrate that for non-hydroxylated quartz surfaces, θ strongly increases with pressure at subcritical and supercritical conditions. An increase in temperature significantly reduces the contact angle, especially at low-intermediate pressures (1-10 MPa), this effect is mitigated at higher pressures, 20 MPa. We also found that θ only weakly depends on salinity for the systems investigated in this work.

Intrinsic profiles and the structure of liquid surfaces.

We present a brief review of the advances made in the characterization of liquid surfaces over the last decade. We focus particularly on the links between the capillary wave theory, the density functional formalism and the direct evaluation of the intrinsic density profiles from computer simulations. A new perspective of the liquid surfaces is appearing, with a sharper view of their molecular structure, which opens new challenges for theoretical and experimental studies. Novel results on the intrinsic interfacial structure of molten salt liquid-vapor interfaces are presented.

Nanoparticles at fluid interfaces.

Nanoparticles at fluid interfaces are becoming a central topic in colloid science studies. Unlike in the case of colloids in suspensions, the description of the forces determining the physical behavior of colloids at interfaces still represents an outstanding problem in the modern theory of colloidal interactions. These forces regulate the formation of complex two-dimensional structures, which can be exploited in a number of applications of technological interest; optical devices, catalysis, molecular electronics or emulsions stabilization. From a fundamental viewpoint and typical for colloidal systems, nanoparticles and microparticles at interfaces are ideal experimental and theoretical models for investigating questions of relevance in condensed matter physics, such as the phase behavior of two-dimensional fluids. This review is a topical survey of the stability, self-assembly behavior and mutual interactions of nanoparticles at fluid interfaces. Thermodynamic models offer an intuitive approach to explaining the interfacial stability of nanoparticles in terms of a few material properties, such as the surface and line tensions. A critical discussion of the theoretical basis, accuracy, limitations, and recent predictions of the thermodynamic models is provided. We also review recent work concerned with nanoparticle self-assembly at fluid interfaces. Complex two-dimensional structures varying considerably with the particle nature have been observed in a number of experiments. We discuss the self-assembly behavior in terms of nanoparticle composition, focusing on sterically stabilized, charged and magnetic nanoparticles. The structure of the two-dimensional assemblies is a reflection of complex intercolloidal forces. Unlike the case for bulk colloidal suspensions, which often can be described reasonably well using DLVO (Derjaguin-Landau-Verwey-Overbeek) theory, the description of particles at interfaces requires the consideration of interfacial deformations as well as interfacial thermal fluctuations. We analyze the importance of both deformation and fluctuations, as well as the modification of electrostatic and van der Waals interactions. Finally, we discuss possible future directions in the field of nanoparticles at interfaces.

Algorithms for the computation of solutions of the Ornstein-Zernike equation.

We introduce a robust and efficient methodology to solve the Ornstein-Zernike integral equation using the pseudoarc length (PAL) continuation method that reformulates the integral equation in an equivalent but nonstandard form. This enables the computation of solutions in regions where the compressibility experiences large changes or where the existence of multiple solutions and so-called branch points prevents Newton's method from converging. We illustrate the use of the algorithm with a difficult problem that arises in the numerical solution of integral equations, namely the evaluation of the so-called no-solution line of the Ornstein-Zernike hypernetted chain (HNC) integral equation for the Lennard-Jones potential. We are able to use the PAL algorithm to solve the integral equation along this line and to connect physical and nonphysical solution branches (both isotherms and isochores) where appropriate. We also show that PAL continuation can compute solutions within the no-solution region that cannot be computed when Newton and Picard methods are applied directly to the integral equation. While many solutions that we find are new, some correspond to states with negative compressibility and consequently are not physical.

Interactions of polarizable media in water: a molecular dynamics approach.

We investigate the interactions of polarizable solutes in water as a function of the solute permittivity. A generic and computationally efficient simulation methodology for the investigation of systems involving dielectric discontinuities is introduced. We report results for interactions between two polarizable cylindrical solutes of nanometer dimensions, which demonstrate that the interactions between the solutes strongly depend on the solute permittivity epsilon. For low permittivity, epsilon approximately 1-2, the interactions are dominated by surface tension forces whose origin lies in the formation of a vapor cavity between the two hydrophobic solutes. This effect leads to a drying transition, where the intersolute force changes discontinuously at a specific solute-solute separation. We find that a moderate permittivity, epsilon approximately 20, enhances the solvation of the polarizable objects inhibiting this drying transition. In the limit of moderately high permittivity, the interactions are dominated by solvation forces. These forces are much larger than those calculated using macroscopic models of dielectrics, which consider water as a continuum dielectric medium. Our results emphasize the importance of including the solvent explicitly to investigate dielectric discontinuities and interactions between polarizable media in water.

Order-disorder transition in the solid phase of a charged hard sphere model.

We investigate the solid phases of the restricted primitive model (RPM). Monte Carlo simulations show the existence of an order-disorder transition from a substitutionally disordered face centered cubic lattice (fcc) to a new ordered fcc structure which is proposed as the ground state of the RPM at the close packing density. Our results suggest that the new phase might turn out in a new triple point in the RPM phase diagram involving three solid phases: CsCl, fcc ordered and fcc disordered structures. The order-disorder transition is also studied using the cell theory. The theory shows good agreement with the simulation results and suggests that the transition is weakly first order.