José Alejandre

Molecular dynamics simulation of the orthobaric densities and surface tension of water

Molecular dynamics simulations have been performed to study the liquid–vapor equilibrium of water as a function of temperature. The orthobaric densities and the surface tension of water are reported for temperatures from 316 K until 573 K. The extended simple point charge (SPC/E) interaction potential for water molecules is used with full Ewald summation. The normal and tangential components of the pressure tensor were calculated and are presented at 328 K. The nature of the long‐range contribution to the surface tension has been studied in detail. At 328 K the calculated surface tension is 66.0±3.0 mN m−1 in comparison with the experimental value of 67 mN m−1. The simulated surface tensions between 316 K and 573 K are in good agreement with experiment. The orthobaric densities are in better agreement with experimental values than those obtained from the Gibbs ensemble calculation for the SPC model of water.

Computer simulations of liquid/vapor interface in Lennard-Jones fluids: Some questions and answers

Canonical molecular dynamics (MD) and Monte Carlo (MC) simulations for liquid/vapor equilibrium in truncated Lennard-Jones fluid have been carried out. Different results for coexistence properties (orthobaric densities, normal and tangential pressure profiles, and surface tension) have been reported in each method. These differences are attributed in literature to different set up conditions, e.g., size of simulation cell, number of particles, cut-off radius, time of simulations, etc., applied by different authors. In the present study we show that observed disagreement between simulation results is due to the fact that different authors inadvertently simulated different model fluids. The origin of the problem lies in details of truncation procedure used in simulation studies. Care has to be exercised in doing the comparison between both methods because in MC calculations one deals with the truncated potential, while in MD calculations one uses the truncated forces, i.e., derivative of the potential. The ...

A Liouville-operator derived measure-preserving integrator for molecular dynamics simulations in the isothermal-isobaric ensemble

The constant-pressure, constant-temperature (NPT) molecular dynamics approach is re-examined from the viewpoint of deriving a new measure-preserving reversible geometric integrator for the equations of motion. The underlying concepts of non-Hamiltonian phase-space analysis, measure-preserving integrators and the symplectic property for Hamiltonian systems are briefly reviewed. In addition, current measure-preserving schemes for the constant-volume, constant-temperature ensemble are also reviewed. A new geometric integrator for the NPT method is presented, is shown to preserve the correct phase-space volume element and is demonstrated to perform well in realistic examples. Finally, a multiple time-step version of the integrator is presented for treating systems with motion on several time scales.

Electrostatic interactions in dissipative particle dynamics using the Ewald sums

The electrostatic interactions in dissipative particle dynamics (DPD) simulations are calculated using the standard Ewald [Ann. Phys. 64, 253 (1921)] sum method. Charge distributions on DPD particles are included to prevent artificial ionic pair formation. This proposal is an alternative method to that introduced recently by Groot [J. Chem. Phys. 118, 11265 (2003)] where the electrostatic field was solved locally on a lattice. The Ewald method is applied to study a bulk electrolyte and polyelectrolyte-surfactant solutions. The structure of the fluid is analyzed through the radial distribution function between charged particles. The results are in good agreement with those reported by Groot for the same systems. We also calculated the radius of gyration of a polyelectrolyte in salt solution as a function of solution pH and degree of ionization of the chain. The radius of gyration increases with the net charge of the polymer in agreement with the trend found in static light scattering experiments of polystyrene sulfonate solutions.

Oscillatory surface tension due to finite-size effects

The simulation results of surface tension at the liquid-vapor interface are presented for fluids interacting with Lennard Jones and square-well potentials. From the simulation of liquids we have reported [M. González-Melchor et al., J. Chem. Phys. 122, 4503 (2005)] that the components of pressure tensor in parallelepiped boxes are not the same when periodic boundary conditions and small transversal areas are used. This fact creates an artificial oscillatory stress anisotropy in the system with even negative values. By doing direct simulations of interfaces we show in this work that surface tension has also an oscillatory decay at small surface areas; this behavior is opposite to the monotonic decay reported previously for the Lennard Jones fluid. It is shown that for small surface areas, the surface tension of the square-well potential artificially takes negative values and even increases with temperature. The calculated surface tension using a direct simulation of interfaces might have two contributions: one from finite-size effects of interfacial areas due to box geometry and another from the interface. Thus, it is difficult to evaluate the true surface tension of an interface when small surface areas are used. Care has to be taken to use the direct simulation method of interfaces to evaluate the predicted surface tension as a function of interfacial area from capillary-wave theory. The oscillations of surface tension decay faster at temperatures close to the critical point. It is also discussed that a surface area does not show any important effect on coexisting densities, making this method reliable to calculate bulk coexisting properties using small systems.

Molecular dynamics for discontinuous potentials

A method to perform molecular dynamics simulations for systems of particles interacting with discontinuous potentials is presented. It is a generalization of Alder and Wainwrights algorithm for potentials represented by an arbitrary number of discontinuous horizontal line segments. The method is applied to several hard molecular fluids of various shapes, discrete Lennard-Jonesium, surfaces and mixtures in order to test its generality and flexibility. Results are compared, when possible, with previously published material, which includes numerical, theoretical and experimental results. These comparisons show that the method is a powerful tool in the study of molecular fluids and looks like a promising alternative approach to perform molecular dynamics for continuous potentials represented by a series of discontinuous line segments. The flexibility of the method will allow the testing of theories based on integral equations, like RISM and RAM, and to develop equations of state for complex hard molecules to...

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.

The surface tension of TIP4P/2005 water model using the Ewald sums for the dispersion interactions

The liquid-vapor phase equilibria and surface tension of the TIP4P/2005 water model is obtained by using the Ewald summation method to determine the long range Lennard-Jones and electrostatic interactions. The method is implemented in a straightforward manner into standard simulation programs. The computational cost of using Ewald sums in dispersion interactions of water is estimated in direct simulation of interfaces. The results of this work at 300 K show a dramatic change in surface tension with an oscillatory behavior for surface areas smaller than 5x5sigma(2), where sigma is the Lennard-Jones oxygen diameter. The amplitude of such oscillations substantially decreases with temperature. Finite size effects are less important on coexisting densities. Phase equilibria and interfacial properties can be determined using a small number of water molecules; their fluctuations are around the same size of simulation error at all temperatures, even in systems where the interfaces are separated a few molecular diameters only. The difference in surface tension of this work compared to the results of other authors is not significant (on the contrary, there is a good agreement). What should be stressed is the different and more consistent approach to obtain the surface tension using the Ewald sums for dispersion interactions. There are two relevant aspects at the interface: An adsorption of water molecules is observed at small surface areas and its thickness systematically increases with system size.

The intrinsic structure of the water surface

An operational procedure to obtain the intrinsic structure of liquid surfaces is applied here to a molecular dynamics simulation of water, with a model of point charges for the molecular interactions. The method, which had been recently proposed and used for simple fluids, is successfully extended to a molecular liquid with the complex bond structure of water. The elimination of the capillary wave fluctuations, in the intrinsic density and orientation profiles, gives a new overall view of the water surface, at the sharpest molecular level, and without the size-dependent broadening observed in the mean profiles. The molecules belonging to the outer liquid layer are clearly identified, and we find that only these molecules exhibit a clear preferential orientation to lie flat on the surface. Moreover, there is a strong correlation between the dipolar structure and the local curvatures of the intrinsic surface, so that at the extrusions of the intrinsic surface the molecular dipoles point preferentially toward the vapor side of the interface. Finally, we have found an intrinsic density layering structure, although the inner structure is strongly damped beyond the second layer.

Surface tension of a square well fluid

We performed Monte Carlo simulations in the canonical ensemble on the liquid–vapor interface of a square well fluid with interaction range of λ=1.5σ. The system contains a liquid slab surrounded by vapor. The surface tension is calculated during simulations by using an original procedure that allows the calculation of the pressure tensor components. The surface tension decreases monotonically with temperature. Coexisting densities and pressure along the liquid–vapor coexistence line have also been obtained and good agreement is found with results calculated from bulk simulations.