A GPU-based Large-scale Monte Carlo Simulation Method for Systems with Long-range Interactions
arXiv:1602.05716 · doi:10.1016/j.jcp.2017.02.069
Abstract
In this work we present an efficient implementation of Canonical Monte Carlo simulation for Coulomb many body systems on graphics processing units (GPU). Our method takes advantage of the GPU Single Instruction, Multiple Data (SIMD) architectures. It adopts the sequential updating scheme of Metropolis algorithm, and makes no approximation in the computation of energy. It reaches a remarkable 440-fold speedup, compared with the serial implementation on CPU. We use this method to simulate primitive model electrolytes. We measure very precisely all ion-ion pair correlation functions at high concentrations, and extract renormalized Debye length, renormalized valences of constituent ions, and renormalized dielectric constants. These results demonstrate unequivocally physics beyond the classical Poisson-Boltzmann theory.
References in corpus (8)
- Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model
- Massively parallel Monte Carlo for many-particle simulations on GPUs
- Gravitational tree-code on graphics processing units: implementation in CUDA
- The Chamomile Scheme: An Optimized Algorithm for N-body simulations on Programmable Graphics Processing Units
- Simulating Lattice Spin Models on Graphics Processing Units
- A Multi-scale Monte Carlo Method for Electrolytes
- Efficient Implementation of the Barnes-Hut Octree Algorithm for Monte Carlo Simulations of Charged Systems
- Charge Renormalization and Charge Oscillation in Asymmetric Primitive Model
Cited by in corpus (4)
- A fast GPU Monte Carlo Radiative Heat Transfer Implementation for Coupling with Direct Numerical Simulation
- A High-Performance Implementation of Atomistic Spin Dynamics Simulations on x86 CPUs
- Fractional Brownian motion meets topology: statistical and topological properties of globular macromolecules with volume interactions
- Harmonic surface mapping algorithm for electrostatic potentials in an atomistic/continuum hybrid model for electrolyte solutions