Fernando Bresme

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.

Equilibrium and nonequilibrium molecular-dynamics simulations of the central force model of water

Equilibrium and nonequilibrium molecular-dynamics simulations of the central force model of water (CFM) [Lemberg and Stillinger, J. Chem. Phys. 62, 1677 (1975)] are presented. We consider a model based on a functional form introduced in theoretical studies of associating systems employing integral equations [F. Bresme, J. Chem. Phys. 108, 4505 (1998)]. Results on thermodynamic, dynamic, dielectric, and coexistence properties are presented. The central force model shows satisfactory agreement with the experimental results in all these cases. In addition, nonequilibrium molecular-dynamics simulations show that the CFM predicts a decrease of the thermal conductivity with temperature, as observed in the experiment, but this dependence is reproduced qualitatively at temperatures characteristic of supercooled states. These results emphasize the need for further studies of the heat conduction and properties of water in these conditions. Overall the present potential should provide a basis for further theoretical...

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.

Molecular dynamics simulations of the surface tension of ionic liquids

We report molecular dynamics computer simulations of the surface tension and interfacial thickness of ionic liquid-vapor interfaces modeled with a soft core primitive model potential. We find that the surface tension shows an anomalous oscillatory behavior with interfacial area. This observation is discussed in terms of finite size effects introduced by the periodic boundary conditions employed in computer simulations. Otherwise we show that the thickness of the liquid-vapor interface increases with surface area as predicted by the capillary wave theory. Data on the surface tension of size-asymmetric ionic liquids are reported and compared with experimental data of molten salts. Our data suggest that the surface tensions of size-asymmetric ionic liquids do not follow a corresponding states law.

Heat transfer in protein-water interfaces.

We investigate using transient non-equilibrum molecular dynamics simulation the temperature relaxation process of three structurally different proteins in water, namely; myoglobin, green fluorescence protein (GFP) and two conformations of the Ca(2+)-ATPase protein. By modeling the temperature relaxation process using the solution of the heat diffusion equation we compute the thermal conductivity and thermal diffusivity of the proteins, as well as the thermal conductance of the protein-water interface. Our results indicate that the temperature relaxation of the protein can be described using a macroscopic approach. The protein-water interface has a thermal conductance of the order of 100-270 MW K(-1) m(-2), characteristic of water-hydrophilic interfaces. The thermal conductivity of the proteins is of the order of 0.1-0.2 W K(-1) m(-1) as compared with approximately 0.6 W K(-1) m(-1) for water, suggesting that these proteins can develop temperature gradients within the biomolecular structures that are larger than those of aqueous solutions. We find that the thermal diffusivity of the transmembrane protein, Ca(2+)-ATPase is about three times larger than that of myoglobin or GFP. Our simulation shows that the Kapitza length of these structurally different proteins is of the order of 1 nm, showing that the protein-water interface should play a major role in defining the thermal relaxation of biomolecules.

Nonequilibrium molecular dynamics simulations of the thermal conductivity of water: A systematic investigation of the SPC/E and TIP4P/2005 models

We report an extensive nonequilibrium molecular dynamics investigation of the thermal conductivity of water using two of the most accurate rigid nonpolarizable empirical models available, SPC/E and TIP4P/2005. Our study covers liquid and supercritical states. Both models predict the anomalous increase of the thermal conductivity with temperature and the thermal conductivity maximum, hence confirming their ability to reproduce the complex anomalous behaviour of water. The performance of the models strongly depends on the thermodynamic state investigated, and best agreement with experiment is obtained for states close to the liquid coexistence line and at high densities and temperatures. Considering the simplicity of these two models the overall agreement with experiments is remarkable. Our results show that explicit polarizability and molecular flexibility are not needed to reproduce the anomalous heat conduction of water.

Stress anisotropy induced by periodic boundary conditions.

Finite size effects due to periodic boundary conditions are investigated using computer simulations in the canonical ensemble. We study liquids with densities corresponding to typical liquid coexistence densities, and temperatures between the triple and critical points. The components of the pressure tensor are computed in order to analyze the finite size effects arising from the size and geometry of the simulation box. Two different box geometries are considered: cubic and parallelepiped. As expected the pressure tensor is isotropic in cubic boxes, but it becomes anisotropic for small noncubic boxes. We argue this is the origin of the anomalous behavior observed recently in the computation of the surface tension of liquid-vapor interfaces. Otherwise, we find that the bulk pressure is sensitive to the box geometry when small simulation boxes are considered. These observations are general and independent of the model liquid considered. We report results for liquids interacting through short range forces, square well and Lennard-Jones, and also long range Coulombic interactions. The effect that small surface areas have on the surface tension is discussed, and some preliminary results at the liquid vapor-interface for the square well potential are given.

The short range anion-H interaction is the driving force for crystal formation of ions in water.

The crystal formation of NaCl in water is studied by extensive molecular dynamics simulations. Ionic solutions at room temperature and various concentrations are studied using the SPC/E and TIP4P/2005 water models and seven force fields of NaCl. Most force fields of pure NaCl fail to reproduce the experimental density of the crystal, and in solution some favor dissociation at saturated conditions, while others favor crystal formation at low concentration. A new force field of NaCl is proposed, which reproduces the experimental phase diagram in the solid, liquid, and vapor regions. This force field overestimates the solubility of NaCl in water at saturation conditions when used with standard Lorentz-Berthelot combining rules for the ion-water pair potentials. It is shown that precipitation of ions is driven by the short range interaction between Cl-H pairs, a term which is generally missing in the simulation of ionic solutions. The effects of intramolecular flexibility of water on the solubility of NaCl ions are analyzed and is found to be small compared to rigid models. A flexible water model, extending the rigid SPC/E, is proposed, which incorporates Lennard-Jones interactions centered on the hydrogen atoms. This force field gives liquid-vapor coexisting densities and surface tensions in better agreement with experimental data than the rigid SPC/E model. The Cl-H, Na-O, and Cl-O pair distribution functions of the rigid and flexible models agree well with experiment. The predicted concentration dependence of the electric conductivity is in fair agreement with available experimental data.

Stability of particles adsorbed at liquid/fluid interfaces: Shape effects induced by line tension

We analyze the stability of particles at liquid–fluid interfaces as a function of their shape and orientation with respect to the interface plane. To do this we apply the generalized Young’s equation, including the line and surface tensions, to the specific shapes of prolate and oblate particles. We show that particles with aspect ratio α larger than a critical value (depending on the surface and line tensions) cannot be adsorbed at the interface. The stability of nonspherical particles at interfaces exhibits a rather complex dependence with line tension and orientation. The interplay between these variables results in a destabilization of prolate particles. In general the stability decreases in the order, oblate>sphere>prolate. It is argued that in experiments aiming to investigate the effect and magnitude of the line tension, the optimum experimental choice corresponds to very elongated particles.

The fluid-solid equilibrium for a charged hard sphere model revisited

The global phase diagram of a system of charged hard spheres, composed of positive and negative ions of the same size, is obtained by means of computer simulations. Thermodynamic integration and Einstein crystal calculations are used to determine the free energies of the different possible solid structures. In this way, the fluid–solid and solid–solid phase transitions are located. Gibbs–Duhem integration is used to trace the full coexistence curves between the different phases involved. Three different solid structures are found to be stable for the model considered; namely, a cesium chloride structure (CsCl), a substitutionally disordered close packed structure which is faced centered cubic (fcc), and a tetragonal ordered structure with a fcc arrangement of atoms if the charge of the ions is not considered. At high temperatures, freezing leads to the substitutionally disordered close packed structure. This solid structure undergoes an order–disorder transition at low temperatures transforming into the t...