Alchemical osmostat for MC simulation: Sampling aqueous electrolyte solution in open systems

10 November 2022, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

Molecular simulations involving electrolytes are usually performed at a fixed amount of salt ions in the simulation box, reproducing macroscopic concentration. Although this statement is valid in the bulk, the concentration of an electrolyte confined in a nanoporous materials such as MOFs or zeolites is greatly affected and remains a priori unknown. The nanoporous material in equilibrium with the bulk electrolyte exchange water and ions at a given chemical potential in the semigrand canonical ensemble that must be calibrated in order to determine the concentration in the nanoporous material. In this work, we propose an algorithm based on non-equilibrium candidate Monte Carlo (NCMC) moves to ultimately perform MC simulations in contact with a saline reservoir. First, we adapt the Widom insertion technique to calibrate the chemical potential by alchemically transmuting water molecules into ions by using NCMC moves. The chemical potential defines a Monte carlo osmostat in the semigrand constant volume and temperature ensemble to be added in a Monte Carlo simulation where the number of ions fluctuates. In order to validate the method, we adapted the NCMC move to determine the free energy of water solvation and subsequently explore thermodynamics of electrolyte solvation at infinite dilution in water. Finally, we implemented the osmostat in MC simulations initialized with bulk water that are driven towards electrolytes of similar concencentration as the saline reservoir. Our results demonstrate that alchemical osmostat MC simulation is a promising tool for use to sample electrolyte insertion in nanoporous materials.

Supplementary materials

Title
Description
Actions
Title
Supporting information
Description
Supporting information
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.