A Multi-scale Monte Carlo Method for Electrolytes

TL;DR

This paper introduces a multi-scale Monte Carlo method for electrolyte simulations that avoids artifacts of traditional periodic boundary conditions by explicitly modeling ions inside a cavity and implicitly outside, improving accuracy and efficiency.

Contribution

The paper presents a novel multi-scale Monte Carlo approach combining explicit cavity simulation with continuum outside treatment, derived effective interactions, and a fast algorithm for improved electrolyte modeling.

Findings

Accurately captures electrolyte physics with smaller simulation scales.

Reduces artifacts caused by periodic boundary conditions.

Demonstrates improved efficiency over traditional methods.

Abstract

Artifacts arise in the simulations of electrolytes using periodic boundary conditions (PBC). We show the origin of these artifacts are the periodic image charges and the constraint of charge neutrality inside the simulation box, both of which are unphysical from the view point of real systems. To cure these problems, we introduce a multi-scale Monte Carlo method, where ions inside a spherical cavity are simulated explicitly, whilst ions outside are treated implicitly using continuum theory. Using the method of Debye charging, we explicitly derive the effective interactions between ions inside the cavity, arising due to the fluctuations of ions outside. We find that these effective interactions consist of two types: 1) a constant cavity potential due to the asymmetry of the electrolyte, and 2) a reaction potential that depends on the positions of all ions inside. Combining the Grand Canonical Monte Carlo (GCMC) with a recently developed fast algorithm based of image charge method, we perform a multi-scale Monte Carlo simulation of symmetric electrolytes, and compare it with other simulation methods, including PBC+GCMC method, as well as large scale Monte Carlo simulation. We demonstrate that our multi-scale MC method is capable of capturing the correct physics of a large system using a small scale simulation.