Boltzmann bias grand canonical Monte Carlo
arXiv:0801.4622 · doi:10.1063/1.2883683
Abstract
We derive an efficient method for the insertion of structured particles in grand canonical Monte Carlo simulations of adsorption in very confining geometries. We extend this method to path integral simulations and use it to calculate the isotherm of adsorption of hydrogen isotopes in narrow carbon nanotubes (2D confinement) and slit pores (1D confinement) at the temperatures of 20 K and 77 K, discussing its efficiency by comparison to the standard path integral grand canonical Monte Carlo algorithm. We use this algorithm to perform multicomponent simulations in order to calculate the hydrogen isotope selectivity for adsorption in narrow carbon nanotubes and slit pores at finite pressures. The algorithm described here can be applied to the study of adsorption of real oligomers and polymers in narrow pores and channels.
12 pages, 6 figures. Accepted for publication in The Journal of Chemical Physics
References in corpus (1)
Cited by in corpus (5)
- Path-integral calculation of the third virial coefficient of quantum gases at low temperatures
- Ab initio Calculation of Fluid Properties for Precision Metrology
- Path-integral calculation of the third dielectric virial coefficient of noble gases
- Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions
- Donnan equilibrium in charged slit-pores from a hybrid nonequilibrium Molecular Dynamics / Monte Carlo method with ions and solvent exchange