Abstract:The Grand Canonical Monte Carlo (GCMC) ensemble defined by the excess chemical potential, μex , volume, and temperature, in the context of molecular simulations allows for variations in the number of particles in the system. In practice, GCMC simulations have been widely applied for the sampling of rare gasses and water, but limited in the context of larger molecules. To overcome this limitation, the oscillating μex GCMC method was introduced and shown to be of utility for sampling small solutes, such as formamide, propane, and benzene, as well as for ionic species such as monocations, acetate, and methylammonium. However, the acceptance of GCMC insertions is low, and the method is computationally demanding. In the present study, we improved the sampling efficiency of the GCMC method using known cavity-bias and configurational-bias algorithms in the context of GPU architecture. Specifically, for GCMC simulations of aqueous solution systems, the configurational-bias algorithm was extended by applying system partitioning in conjunction with a random interval extraction algorithm, thereby improving the efficiency in a highly parallel computing environment. The method is parallelized on the GPU using CUDA and OpenCL, allowing for the code to run on both Nvidia and AMD GPUs, respectively. Notably, the method is particularly well suited for GPU computing as the large number of threads allows for simultaneous sampling of a large number of configurations during insertion attempts without additional computational overhead. In addition, the partitioning scheme allows for simultaneous insertion attempts for large systems, offering considerable efficiency. Calculations on the BK Channel, a transporter, including a lipid bilayer with over 760,000 atoms, show a speed up of ~53-fold through the use of system partitioning. The improved algorithm is then combined with an enhanced μex oscillation protocol and shown to be of utility in the context of the site-identification by ligand competitive saturation (SILCS) co-solvent sampling approach as illustrated through application to the protein CDK2.

Rapid screening of gas solubility in ionic liquids using biased particle insertions with pre-sampled liquid trajectories

Kinetic Study Of Rapid Transfer Of Tetraethylammonium At The 1,2-Dichloroethane/Water Interface By Nanopipet Voltammetry Of Common Ions

Rapid Measurement of Gas Solubility in Ionic Liquids with a Simple Tube-in-Tube Reactor

Rapid Measurement of Gas Solubility in Liquids Using a Membrane Dispersion Microcontactor

High-Throughput Screening and Accurate Prediction of Ionic Liquid Viscosities Using Interpretable Machine Learning

Non-destructive molecular-dynamics simulation of the chemical potential of a fluid

Probing gas adsorption in MOFs using an efficient ab initio widom insertion Monte Carlo method

Spatially Resolved Thermodynamic Integration: An Efficient Method to Compute Chemical Potentials of Dense Fluids

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

Chemical potential of quadrupolar two-centre Lennard-Jones fluids by gradual insertion

A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals

Accelerated convergence via adiabatic sampling for adsorption and desorption processes

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman-Rubin Convergence Diagnostics

Improving the Efficiency of Non-equilibrium Sampling in the Aqueous Environment via Implicit-Solvent Simulations.

Molecular simulation of concurrent gas-liquid interfacial adsorption and partitioning in gas-liquid chromatography

Comparative assessment of physics-based in silico methods to calculate relative solubilities

GPU-specific algorithms for improved solute sampling in grand canonical Monte Carlo simulations

A Thin Film Approach to Rapid, Quantitative Measurements of Mixed-Gas Adsorption Equilibrium in Nanoporous Materials

Electronic screening using a virtual Thomas-Fermi fluid for predicting wetting and phase transitions of ionic liquids at metal surfaces

On the use of QM/MM Surface Hopping simulations to understand thermally-activated rare event nonadiabatic transitions in the condensed phase

Lattice Boltzmann based discrete simulation for gas-solid fluidization