No abstract

One of the long standing challenges in molecular simulation is the description of interfaces. On the molecular length scale, finite size effects significantly influence the properties of the interface such as its interfacial tension, which can be reliably investigated by molecular dynamics simulation of planar vapor-liquid interfaces. For the Lennard-Jones fluid, finite size effects are examined here by varying the thickness of the liquid slab. It is found that the surface tension and density in the center of the liquid region decreases significantly for thin slabs. The influence of the slab thickness on both the liquid density and the surface tension is found to scale with 1/S 3 in terms of the slab thickness S, and a linear correlation between both effects is obtained. The results corroborate the analysis of Malijevský and Jackson, J. Phys.: Cond. Mat. 24: 464121 (2012), who recently detected an analogous effect for the surface tension of liquid nanodroplets.

The molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular dynamics simulation. With an efficient dynamic load balancing scheme, it delivers high scalability even for challenging heterogeneous configurations. * To whom correspondence should be addressed † High Performance Computing Center Stuttgart (HLRS), Germany ‡ Laboratory of Engineering Thermodynamics (LTD), Univ. of Kaiserslautern, Germany ¶ Scientific Computing in Computer Science (SCCS), TU München, Germany § Thermodynamics and Energy Technology (ThEt), Univ. of Paderborn, Germany 1 Presently, multi-center rigid potential models based on Lennard-Jones sites, point charges and higher-order polarities are supported. Due to its modular design, ls1 mardyn can be extended to new physical models, methods, and algorithms, allowing future users to tailor it to suit their respective needs. Possible applications include scenarios with complex geometries, e.g.for fluids at interfaces, as well as non-equilibrium molecular dynamics simulation of heat and mass transfer.

No abstract

The vapor-liquid equilibrium (VLE) of the Mie potential, where the dispersive exponent is constant (m = 6) while the repulsive exponent n is varied between 9 and 48, is systematically investigated by molecular simulation. For systems with planar vapor-liquid interfaces, long-range correction expressions are derived, so that interfacial and bulk properties can be computed accurately. The present simulation results are found to be consistent with the available body of literature on the Mie fluid which is substantially extended. On the basis of correlations for the considered thermodynamic properties, a multicriteria optimization becomes viable. Thereby, users can adjust the three parameters of the Mie potential to the properties of real fluids, weighting different thermodynamic properties according to their importance for a particular application scenario. In the present work, this is demonstrated for carbon dioxide for which different competing objective functions are studied which describe the accuracy of the model for representing the saturated liquid density, the vapor pressure and the surface tension. It is shown that models can be found which describe simultaneously the saturated liquid density and vapor pressure with good accuracy, and it is discussed to what extent this accuracy can be upheld as the model accuracy for the surface tension is further improved.