Parallel kinetic Monte Carlo simulation of Coulomb glasses
CONICET, Centro Atomico Bariloche, 8400 Bariloche, Argentina
arXiv:1407.5026 [cond-mat.dis-nn], (18 Jul 2014)
@article{2014arXiv1407.5026F,
author={Ferrero}, E.~E. and {Kolton}, A.~B. and {Palassini}, M.},
title={"{Parallel kinetic Monte Carlo simulation of Coulomb glasses}"},
journal={ArXiv e-prints},
archivePrefix={"arXiv"},
eprint={1407.5026},
primaryClass={"cond-mat.dis-nn"},
keywords={Condensed Matter – Disordered Systems and Neural Networks, Condensed Matter – Statistical Mechanics, Physics – Computational Physics},
year={2014},
month={jul},
adsurl={http://adsabs.harvard.edu/abs/2014arXiv1407.5026F},
adsnote={Provided by the SAO/NASA Astrophysics Data System}
}
We develop a parallel rejection algorithm to tackle the problem of low acceptance in Monte Carlo methods, and apply it to the simulation of the hopping conduction in Coulomb glasses using Graphics Processing Units, for which we also parallelize the update of local energies. In two dimensions, our parallel code achieves speedups of up to two orders of magnitude in computing time over an equivalent serial code. We find numerical evidence of a scaling relation for the relaxation of the conductivity at different temperatures.
July 24, 2014 by hgpu
