Balázs Fábián

Pressure Profile Calculation with Mesh Ewald Methods

The importance of calculating pressure profiles across liquid interfaces is increasingly gaining recognition, and efficient methods for the calculation of long-range contributions are fundamental in addressing systems with a large number of charges. Here, we show how to compute the local pressure contribution for mesh-based Ewald methods, retaining the typical N log N scaling as a function of the lattice nodes N. This is a considerable improvement on existing methods, which include approximating the electrostatic contribution using a large cutoff and the, much slower, Ewald calculation. As an application, we calculate the contribution to the pressure profile across the water/vapor interface, coming from different molecular layers, both including and removing the effect of thermal capillary waves. We compare the total pressure profile with the one obtained using the cutoff approximation for the calculation of the stresses, showing that the stress distributions obtained using the Harasima and Irving-Kirkwood path are quite similar and shifted with respect to each other at most 0.05 nm.

Nonzero Ideal Gas Contribution to the Surface Tension of Water

Surface tension, the tendency of fluid interfaces to behave elastically and minimize their surface, is routinely calculated as the difference between the lateral and normal components of the pressure or, invoking isotropy in momentum space, of the virial tensor. Here we show that the anisotropy of the kinetic energy tensor close to a liquid-vapor interface can be responsible for a large part of its surface tension (about 15% for water, independent from temperature).

Miscibility and Thermodynamics of Mixing of Different Models of Formamide and Water in Computer Simulation

The thermodynamic changes that occur upon mixing five models of formamide and three models of water, including the miscibility of these model combinations itself, is studied by performing Monte Carlo computer simulations using an appropriately chosen thermodynamic cycle and the method of thermodynamic integration. The results show that the mixing of these two components is close to the ideal mixing, as both the energy and entropy of mixing turn out to be rather close to the ideal term in the entire composition range. Concerning the energy of mixing, the OPLS/AA_mod model of formamide behaves in a qualitatively different way than the other models considered. Thus, this model results in negative, while the other ones in positive energy of mixing values in combination with all three water models considered. Experimental data supports this latter behavior. Although the Helmholtz free energy of mixing always turns out to be negative in the entire composition range, the majority of the model combinations tested either show limited miscibility, or, at least, approach the miscibility limit very closely in certain compositions. Concerning both the miscibility and the energy of mixing of these model combinations, we recommend the use of the combination of the CHARMM formamide and TIP4P water models in simulations of water-formamide mixtures.

Local lateral environment of the molecules at the surface of DMSO-water mixtures

Molecular dynamics simulations of the liquid-vapour interface of dimethyl sulphoxide (DMSO)-water mixtures of 11 different compositions, including two neat systems are performed on the canonical (N, V, T) ensemble at 298 K. The molecules constituting the surface layer of these systems are selected by means of the identification of the truly interfacial molecules (ITIM) method, and their local lateral environment at the liquid surface is investigated by performing Voronoi analysis. The obtained results reveal that both molecules prefer to be in a mixed local environment, consisting of both kinds of molecules, at the liquid surface, and this preference is even stronger here than in the bulk liquid phase. Neat-like patches, in which a molecule is surrounded by like neighbours, are not found. However, vacancies that are surrounded solely by water molecules are observed at the liquid surface. Our results show that strongly hydrogen bonded DMSO·H2O complexes, known to exist in the bulk phase of these mixtures, are absent from the liquid surface.

Pytim: A python package for the interfacial analysis of molecular simulations

Pytim is a versatile python framework for the analysis of interfacial properties in molecular simulations. The code implements several algorithms for the identification of instantaneous interfaces of arbitrary shape, and analysis tools written specifically for the study of interfacial properties, such as intrinsic profiles. The code is written in the python language, and makes use of the numpy and scipy packages to deliver high computational performances. Pytim relies on the MDAnalysis library to analyze the trajectory file formats of popular simulation packages such as gromacs, charmm, namd, lammps or Amber, and can be used to steer OpenMM simulations. Pytim can write information about surfaces and surface atomic layers to vtk, cube, and pdb files for easy visualization. The classes of Pytim can be easily customized and extended to include new interfacial algorithms or analysis tools. The code is available as open source and is free of charge.

Single Particle Dynamics at the Intrinsic Surface of Various Apolar, Aprotic Dipolar, and Hydrogen Bonding Liquids As Seen from Computer Simulations

We investigate the single molecule dynamics at the intrinsic liquid/vapor interface of five different molecular liquids (carbon tetrachloride, acetone, acetonitrile, methanol, and water). After assessing that the characteristic residence times in the surface layer are long enough for a meaningful definition of several transport properties within the layer itself, we characterize the dynamics of the individual molecules at the liquid surface by analyzing their normal and lateral mean-square displacements and lateral velocity autocorrelation functions and, in the case of the hydrogen bonding liquids (i.e., water and methanol), also the properties of the hydrogen bonds. Further, dynamical properties as well as the clustering of the molecules residing unusually long in the surface layer are also investigated. The global picture emerging from this analysis is that of a noticeably enhanced dynamics of the molecules at the liquid surface, with diffusion coefficients up to 4 times larger than in the bulk, and the disappearance of the caging effect at the surface of all liquids but water. The dynamics of water is dominated by the strong hydrogen bonding structure also at the liquid surface.