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.

Generalization of Biasing Particle Insertion and Deletion for Grand Canonical Monte Carlo Simulation

Boltzmann bias grand canonical Monte Carlo

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

Monte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm

Unbiased Markov chain Monte Carlo methods with couplings

Generalized HMC using Nambu mechanics for lattice QCD

Enabling Ab-Initio Molecular Dynamics under Bias: The CP2K+SMEAGOL Interface for Integrating Density Functional Theory and Non-Equilibrium Green Functions

Unbiasing fermionic quantum Monte Carlo with a quantum computer

Local bias potential in hyper molecular dynamics method

Hyper Molecular Dynamics with A Local Bias Potential

Multilevel Asymptotic-Preserving Monte Carlo for Particle Simulations

Monte Carlo simulation of granular fluids

Green's function technique for a two-electrode mesoscopic system under bias

Reflection coupling for unadjusted generalized Hamiltonian Monte Carlo in the nonconvex stochastic gradient case

Theory and Monte Carlo simulation of the ideal gas with shell particles in the canonical, isothermal-isobaric, grand canonical, and Gibbs ensembles

Fluctuation without dissipation: Microcanonical Langevin Monte Carlo

A quantum Monte Carlo study of systems with effective core potentials and node nonlinearities

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

A Configurational Bias Monte Carlo Method for Linear and Cyclic Peptides

Introduction to Configuration Path Integral Monte Carlo

Multigrid Monte Carlo Revisited: Theory and Bayesian Inference