Billion-atom synchronous parallel kinetic Monte Carlo simulations of critical 3D Ising systems
- IMDEA-Materiales, Madrid 28040 (Spain)
- Massachusetts Institute of Technology, Cambridge, MA 02139 (United States)
- Lawrence Livermore National Laboratory, Livermore, CA 94551 (United States)
An extension of the synchronous parallel kinetic Monte Carlo (spkMC) algorithm developed by Martinez et al. [J. Comp. Phys. 227 (2008) 3804] to discrete lattices is presented. The method solves the master equation synchronously by recourse to null events that keep all processors' time clocks current in a global sense. Boundary conflicts are resolved by adopting a chessboard decomposition into non-interacting sublattices. We find that the bias introduced by the spatial correlations attendant to the sublattice decomposition is within the standard deviation of serial calculations, which confirms the statistical validity of our algorithm. We have analyzed the parallel efficiency of spkMC and find that it scales consistently with problem size and sublattice partition. We apply the method to the calculation of scale-dependent critical exponents in billion-atom 3D Ising systems, with very good agreement with state-of-the-art multispin simulations.
- OSTI ID:
- 21499769
- Journal Information:
- Journal of Computational Physics, Vol. 230, Issue 4; Other Information: DOI: 10.1016/j.jcp.2010.11.006; PII: S0021-9991(10)00608-X; Copyright (c) 2010 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Synchronous parallel kinetic Monte Carlo for continuum diffusion-reaction systems
New multispin coding algorithm for Monte Carlo simulation of the Ising model