Issue |
EPJ Web of Conferences
Volume 108, 2016
Mathematical Modeling and Computational Physics (MMCP 2015)
|
|
---|---|---|
Article Number | 02016 | |
Number of page(s) | 6 | |
Section | Conference Contributions | |
DOI | https://doi.org/10.1051/epjconf/201610802016 | |
Published online | 09 February 2016 |
https://doi.org/10.1051/epjconf/201610802016
GPU-Accelerated Population Annealing Algorithm: Frustrated Ising Antiferromagnet on the Stacked Triangular Lattice
1 Department of Theoretical Physics and Astrophysics, Faculty of Science, P.J. Šafárik University, Park Angelinum 9, 040 01 Košice, Slovak Republic
2 Applied Mathematics Research Centre, Coventry University, Coventry, CV1 5FB, United Kingdom
3 Landau Institute for Theoretical Physics, 142432 Chernogolovka, Russia
4 Science Center in Chernogolovka, 142432 Chernogolovka, Russia
a e-mail: borovsky.michal@gmail.com
b e-mail: Martin.Weigel@coventry.ac.uk
c e-mail: barash@itp.ac.ru
d e-mail: milan.zukovic@upjs.sk
Published online: 9 February 2016
The population annealing algorithm is a novel approach to study systems with rough free-energy landscapes, such as spin glasses. It combines the power of simulated annealing, Boltzmann weighted differential reproduction and sequential Monte Carlo process to bring the population of replicas to the equilibrium even in the low-temperature region. Moreover, it provides a very good estimate of the free energy. The fact that population annealing algorithm is performed over a large number of replicas with many spin updates, makes it a good candidate for massive parallelism. We chose the GPU programming using a CUDA implementation to create a highly optimized simulation. It has been previously shown for the frustrated Ising antiferromagnet on the stacked triangular lattice with a ferromagnetic interlayer coupling, that standard Markov Chain Monte Carlo simulations fail to equilibrate at low temperatures due to the effect of kinetic freezing of the ferromagnetically ordered chains. We applied the population annealing to study the case with the isotropic intra- and interlayer antiferromagnetic coupling (J2/|J1| = −1). The reached ground states correspond to non-magnetic degenerate states, where chains are antiferromagnetically ordered, but there is no long-range ordering between them, which is analogical with Wannier phase of the 2D triangular Ising antiferromagnet.
© Owned by the authors, published by EDP Sciences, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.