6533b820fe1ef96bd1279157

RESEARCH PRODUCT

Finite-size effects of Kirkwood–Buff integrals from molecular simulations

Pieter KrugerPieter KrugerDick BedeauxThijs J. H. VlugtSondre K. SchnellNoura DawassJean-marc SimonSigne KjelstrupSigne Kjelstrup

subject

Thermodynamic stateGeneral Chemical EngineeringMonte Carlo methodInverse02 engineering and technology01 natural sciencesMolecular dynamicsthermodynamicsKirkwood–Buff integrals0103 physical sciencesfinite-size effectsGeneral Materials ScienceStatistical physicsScalingPhysics010304 chemical physicsIsotropyGeneral Chemistry021001 nanoscience & nanotechnologyCondensed Matter Physicsmolecular dynamicsModeling and SimulationThermodynamic limitPartial derivative0210 nano-technologyInformation Systems

description

The modelling of thermodynamic properties of liquids from local density fluctuations is relevant to many chemical and biological processes. The Kirkwood–Buff (KB) theory connects the microscopic structure of isotropic liquids with macroscopic properties such as partial derivatives of activity coefficients, partial molar volumes and compressibilities. Originally, KB integrals were formulated for open and infinite systems which are difficult to access with standard Molecular Dynamics (MD) simulations. Recently, KB integrals for finite and open systems were formulated (J Phys Chem Lett. 2013;4:235). From the scaling of KB integrals for finite subvolumes, embedded in larger reservoirs, with the inverse of the size of these subvolumes, estimates for KB integrals in the thermodynamic limit are obtained. Two system size effects are observed in MD simulations: (1) effects due to the size of the simulation box and the size of the finite subvolume embedded in the simulation box, and (2) effects due to computing radial distribution functions (RDF) from a closed and finite system. In this study, we investigate the two effects in detail by computing KB integrals using the following methods: (1) Monte Carlo simulations of finite subvolumes of a liquid with an analytic RDF and (2) MD simulations of a WCA mixture for various simulation box sizes, but at the same thermodynamic state. We investigate the effect of the size of the simulation box and quantify the differences compared to KB integrals computed in the thermodynamic limit. We demonstrate that calculations of KB integrals should not be extended beyond half the size of the simulation box. For finite-size effects related to the RDF, we find that the Van der Vegt correction (J Chem Theory Comput. 2013;9:1347) yields the most accurate results. © 2017 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License (http://creativecommons.org/licenses/by-nc-nd/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited, and is not altered, transformed, or built upon in any way.

http://hdl.handle.net/11250/2588274