Connecting density fluctuations and Kirkwood-Buff integrals for finite-size systems

J Chem Phys. 2022 Jan 28;156(4):044502. doi: 10.1063/5.0076744.

Abstract

Kirkwood-Buff integrals (KBIs) connect the microscopic structure and thermodynamic properties of liquid solutions. KBIs are defined in the grand canonical ensemble and evaluated by assuming the thermodynamic limit (TL). In order to reconcile analytical and numerical approaches, finite-size KBIs have been proposed in the literature, resulting in two strategies to obtain their TL values from computer simulations. (i) The spatial block analysis method in which the simulation box is divided into subdomains of volume V to compute density fluctuations. (ii) A direct integration method where a corrected radial distribution function and a kernel that accounts for the geometry of the integration subvolumes are combined to obtain KBI as a function of V. In this work, we propose a method that connects both strategies into a single framework. We start from the definition of finite-size KBI, including the integration subdomain and an asymptotic correction to the radial distribution function, and solve them in Fourier space where periodic boundary conditions are trivially introduced. The limit q → 0, equivalent to the value of the KBI in the TL, is obtained via the spatial block-analysis method. When compared to the latter, our approach gives nearly identical results for all values of V. Moreover, all finite-size effect contributions (ensemble, finite-integration domains, and periodic boundary conditions) are easily identifiable in the calculation. This feature allows us to analyze finite-size effects independently and extrapolates the results of a single simulation to different box sizes. To validate our approach, we investigate prototypical systems, including SPC/E water and aqueous urea mixtures.