Patrice Malfreyt

Calculation of the surface tension from Monte Carlo simulations: Does the model impact on the finite-size effects?

We report two-phase Monte Carlo simulations of the liquid-vapor interface of the Lennard-Jones (LJ) fluids in order to study the impact of the methodology used for the energy calculation on the oscillatory behavior of the surface tension with the system sizes. The surface tension values are illustrated through the LJ parameters of methane. The first methodology uses a standard truncated LJ potential, the second one adds a long range correction (LRC) contribution to the energy into the Metropolis scheme, and the third one uses a LJ potential modified by a polynomial function in order to remove the discontinuities at the cutoff distance. The surface tension is calculated from the mechanical and thermodynamic routes and the LRCs to the surface tension are systematically calculated from appropriate expressions within these definitions. The oscillatory behavior has been studied as a function of the size of the interfacial area and of the length of the dimension perpendicular to the surface. We show that the methodology has an important effect on the oscillatory variation in the surface tension with the system size. This oscillatory variation in the surface tension with the system size is investigated through its intrinsic and LRC contributions. We complete this work by studying the dependence of the surface tension with respect to the cutoff distance when the LRC part to the energy is considered into the Metropolis scheme.

Surface tension of water and acid gases from Monte Carlo simulations

We report direct Monte Carlo (MC) simulations on the liquid-vapor interfaces of pure water, carbon dioxide, and hydrogen sulfide. In the case of water, the recent TIP4P/2005 potential model used with the MC method is shown to reproduce the experimental surface tension and to accurately describe the coexistence curves. The agreement with experiments is also excellent for CO(2) and H(2)S with standard nonpolarizable models. The surface tensions are calculated by using the mechanical and the thermodynamic definitions via profiles along the direction normal to the surface. We also discuss the different contributions to the surface tension due to the repulsion-dispersion and electrostatic interactions. The different profiles of these contributions are proposed in the case of water.

Direct Monte Carlo simulations of the equilibrium properties of n-pentane liquid–vapor interface

Direct MC calculations have been carried out to study the liquid–vapor equilibrium of n-pentane. We have used the local long range correction to the configurational energy within the Metropolis scheme and an algorithm allowing us to select randomly with equal probability two different maximum displacements. The thermal equilibrium conditions are checked by calculating the profiles of the configurational temperature in the vapor and liquid phases. We also check the mechanical equilibrium by calculating the profiles of the normal and tangential pressure components. The normal pressure profile is constant through the interface and in both phases on the conditions that the two parts of the long range corrections to the normal pressure area included in the calculations. The critical densities and temperatures are well predicted and the vapor pressures agree satisfactory with the experimental values within larger statistical fluctuations. The long range corrections to the surface tension are calculated using va...

Monte Carlo versus molecular dynamics simulations in heterogeneous systems: an application to the n-pentane liquid-vapor interface.

The Monte Carlo (MC) and molecular dynamics (MD) methodologies are now well established for computing equilibrium properties in homogeneous fluids. This is not yet the case for the direct simulation of two-phase systems, which exhibit nonuniformity of the density distribution across the interface. We have performed direct MC and MD simulations of the liquid-gas interface of n-pentane using a standard force-field model. We obtained density and pressure components profiles along the direction normal to the interface that can be very different, depending on the truncation and long range correction strategies. We discuss the influence on predicted properties of different potential truncation schemes implemented in both MC and MD simulations. We show that the MD and MC profiles can be made in agreement by using a Lennard-Jones potential truncated via a polynomial function that makes the first and second derivatives of the potential continuous at the cutoff distance. In this case however, the predicted thermodynamic properties (phase envelope, surface tension) deviate from experiments, because of the changes made in the potential. A further readjustment of the potential parameters is needed if one wants to use this method. We conclude that a straightforward use of bulk phase force fields in MD simulations may lead to some physical inconsistencies when computing interfacial properties.

Surface tension of water–alcohol mixtures from Monte Carlo simulations

Monte Carlo simulations are reported to predict the dependence of the surface tension of water-alcohol mixtures on the alcohol concentration. Alcohols are modeled using the anisotropic united atom model recently extended to alcohol molecules. The molecular simulations show a good agreement between the experimental and calculated surface tensions for the water-methanol and water-propanol mixtures. This good agreement with experiments is also established through the comparison of the excess surface tensions. A molecular description of the mixture in terms of density profiles and hydrogen bond profiles is used to interpret the decrease of the surface tension with the alcohol concentration and alcohol chain length.

Electrostatic Interactions in Dissipative Particle Dynamics: Toward a Mesoscale Modeling of the Polyelectrolyte Brushes

We report mesoscopic simulations of bulk electrolytes and polyelectrolyte brushes using the dissipative particle dynamics (DPD) method. The calculation of the electrostatic interactions is carried out using both the Ewald summation method and the particle-particle particle-mesh technique with charges distributed over the particles. The local components of the pressure tensor are calculated using the Irving and Kirkwood, and the method of planes and mechanical equilibrium is demonstrated. The profiles of the normal component of the pressure tensor are shown to be similar for both the Ewald and particle-particle particle-mesh methods for a single polyelectrolyte brush. We show that the PPPM method with the MOP technique is the appropriate choice for simulations of this type. The mesoscale modeling of a strongly stretched polylectrolyte brush formed by strong charged polymer chains at a high grafting density shows that the polyelectrolyte follows the nonlinear osmotic regime, as expected from the calculation of the Gouy-Chapman length and the dimensionless Manning ratio.

Mesoscopic simulation of entanglements using dissipative particle dynamics : Application to polymer brushes

We use a simple spring-spring repulsion to model entanglements between polymers in dissipative particle dynamics (DPD) simulations. The model is applied to a polymer brushes system to study lubrication. We demonstrate that this method leads to mechanical equilibrium in polymer brushes using the normal DPD time step. The number of bond crossings is calculated to provide a quantitative description of the entanglement. We demonstrate that it is possible to avoid 99% of the bond crossings with the values of spring-spring repulsion that can be used without significantly decreasing the time step. A shear force is applied to the system to study the effect of the decrease in the bond crossings on the structure and rheological properties of the brushes. In particular, we show how the friction coefficient increases with the decrease in the bond crossings of the polymers.

Molecular Dynamics Simulations of the Liquid Surface of the Ionic Liquid 1-Hexyl-3-methylimidazolium Bis(trifluoromethanesulfonyl)amide: Structure and Surface Tension

Molecular dynamics simulations of the liquid-vacuum interface of the ionic liquid 1-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)amide were performed with an all-atom force field. Structural properties of the interface, such as orientational ordering and density profiles, were calculated. The hexyl side chain of the cation is likely to protrude outward from the surface layer. There is a region with enhanced density from that of the bulk where the cation preferably slants with the imidazolium ring tending to be perpendicular to the interface. The surface tensions are calculated using mechanical and thermodynamic definitions via profiles along the direction normal to the interface. We also discuss the different contributions to the surface tension due to the repulsion-dispersion and electrostatic interactions. The use of local pressure profiles provides an explanation to the systematic problems encountered by several researchers to obtain accurate values of the surface tension at low temperature. Even when macroscopically the system looks in equilibrium, locally this is not accomplished.

Monte Carlo calculation of the methane-water interfacial tension at high pressures

Monte Carlo simulations have been performed in the Np(N)AT statistical ensemble to study the methane-water mixture as a function of pressure. The interfacial tensions are calculated with different definitions and are reported for pressures from 1 to 50 MPa. The interfacial tensions, coexisting densities, and composition of the methane and water phases are shown to be in good agreement with the corresponding experimental properties. The interfacial region has been described through the profiles of the number of hydrogen bonds, the coordination number of each species, and the different energy contributions. We complete this study by a theoretical investigation of the thermal and mechanical equilibria in the binary mixture. We have also examined the profile of the intrinsic and long range correction parts of the interfacial tension along the normal to the water surface.

Mesoscale modeling of polyelectrolyte brushes with salt.

We report dissipative particle dynamics (DPD) simulations of a polyelectrolyte brush under athermal solvent conditions. The electrostatic interactions are calculated using the particle-particle particle-mesh (PPPM) method with charges distributed over the particles. The polymer beads, counterions, co-ions, and solvent particles are modeled explicitly. The DPD simulations show a dependence of the brush height on the grafting density and the charge fraction that is typical of the nonlinear osmotic brush regime. We report the effect of the addition of salt on the structural properties of the brush. In the case of a polyelectrolyte brush with a high surface coverage, the simulations reproduce the transition between the nonlinear osmotic brush regime where the thickness of the brush is independent of the salt concentration and the salted regime where the brush height decreases weakly with the salt concentration.