Understanding population annealing Monte Carlo simulations
Cornell University , 2021. No. 2102.06611.
Weigel M., Barash L., Shchur L., Janke W.
Population annealing is a recent addition to the arsenal of the practitioner in computer simulations in statistical physics and beyond that is found to deal well with systems with complex free-energy landscapes. Above all else, it promises to deliver unrivaled parallel scaling qualities, being suitable for parallel machines of the biggest caliber. Here we study population annealing using as the main example the two-dimensional Ising model which allows for particularly clean comparisons due to the available exact results and the wealth of published simulational studies employing other approaches. We analyze in depth the accuracy and precision of the method, highlighting its relation to older techniques such as simulated annealing and thermodynamic integration. We introduce intrinsic approaches for the analysis of statistical and systematic errors and provide a detailed picture of the dependence of such errors on the simulation parameters. The results are benchmarked against canonical and parallel tempering simulations.
Algorithm for replica redistribution in an implementation of the population annealing method on a hybrid supercomputer architecture
, , , Computer Physics Communications 2021 Vol. 261 P. 107786
The population annealing method is a promising approach for large-scale simulations because it is potentially scalable on any parallel architecture. We present an implementation of the algorithm on a hybrid program architecture combining CUDA and MPI. The problem is to keep all general-purpose graphics processing unit devices as busy as possible by efficiently redistributing replicas. ...
Added: December 28, 2020
, , et al., European Physical Journal: Special Topics 2017 Vol. 226 No. 4 P. 595-604
Population annealing is a hybrid of sequential and Markov chain Monte Carlo methods geared towards the efficient parallel simulation of systems with complex free-energy landscapes. Systems with first-order phase transitions are among the problems in computational physics that are difficult to tackle with standard methods such as local-update simulations in the canonical ensemble, for example ...
Added: October 20, 2017
, , Physical Review E - Statistical, Nonlinear, and Soft Matter Physics 2018 Vol. 98 No. 022218 P. 1-15
We examine the question of the influence of sparse long-range communications on the synchronization in parallel discrete event simulations (PDES). We build a model of the evolution of local virtual times (LVT) in a conservative algorithm including several choices of local links. All network realizations belong to the small-world network class. We find that synchronization ...
Added: July 13, 2018
GPU-Accelerated Population Annealing Algorithm: Frustrated Ising Antiferromagnet on the Stacked Triangular Lattice
, , et al., EPJ Web of Conferences 2016 Vol. 108 P. 02016-p.1-02016-p.6
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 ...
Added: January 31, 2018
, , et al., JETP Letters 2016 Vol. 103 No. 2 P. 131-136
This article presents new features of the MCSANC v.1.20 program, a Monte Carlo tool for calculation of the next-to-leading order electroweak and QCD corrections to various Standard Model processes. The extensions concern implementation of Drell--Yan-like processes and include a systematic treatment of the photon-induced contribution in proton--proton collisions and electroweak corrections beyond NLO approximation. There ...
Added: June 28, 2018
, , Computer Physics Communications 2013 No. 184 P. 2343-2350
mcsanc is a Monte-Carlo tool based on the SANC (Support for Analytic and Numeric Calculations for experiments at colliders) modules for higher order calculations in hadron collider physics. It allows to evaluate NLO QCD and EW cross sections for Drell-Yan processes (inclusive), associated Higgs and gauge boson production and single-top quark production in s- and ...
Added: June 27, 2018
Algorithm for the replica redistribution in the implementation of parallel annealing method on the hybrid supercomputer architecture
, , , / Cornell University. Series arXiv "math". 2020. No. 2006.00561.
The parallel annealing method is one of the promising approaches for large scale simulations as potentially scalable on any parallel architecture. We present an implementation of the algorithm on the hybrid program architecture combining CUDA and MPI. The problem is to keep all general-purpose graphics processing unit devices as busy as possible redistributing replicas and ...
Added: June 2, 2020
, , , Physical Review E - Statistical, Nonlinear, and Soft Matter Physics 2017 Vol. 96 No. 4 P. 043307-1-043307-7
The Wang-Landau (WL) algorithm has been widely used for simulations in many areas of physics. Our analysis of the WL algorithm explains its properties and shows that the difference of the largest eigenvalue of the transition matrix in the energy space from unity can be used to control the accuracy of estimating the density of ...
Added: October 20, 2017
Повышение производительности программного комплекса для моделирования внутренних гравитационных волн IGW Research с помощью Intel® Parallel Studio XE 2013
, , et al., Фундаментальная и прикладная гидрофизика 2012 Т. 5 № 3 С. 89-95
The new version of the program complex intended for numerical modeling of propagation and transformation of internal gravity waves in the ocean, with a finalized unit calculation of a ray of internal waves and with a paralleling of the program, which can significantly speed up the ongoing computation is presented. As a practical example of ...
Added: October 8, 2012
, , et al., Computer Physics Communications 2017 Vol. 220 P. 341-350
Population annealing is a promising recent approach for Monte Carlo simulations in statistical physics, in particular for the simulation of systems with complex free-energy landscapes. It is a hybrid method, combining importance sampling through Markov chains with elements of sequential Monte Carlo in the form of population control. While it appears to provide algorithmic capabilities ...
Added: October 20, 2017
, , et al., / Cornell University. Series cond-mat "arxiv.org". 2018. No. 1808.04340.
Estimating the density of states of systems with rugged free energy landscapes is a notoriously difficult task of the utmost importance in many areas of physics ranging from spin glasses to biopolymers to quantum computing. Some of the standard approaches suffer from a spurious convergence of the estimates to metastable minima, and these cases are ...
Added: October 17, 2018
The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip
, , , / Cornell University. Series math "arxiv.org". 2013. No. arxiv: 1303.3471.
We consider an initial-boundary value problem for a generalized 2D time-dependent Schrödinger equation on a semi-infinite strip. For the Crank-Nicolson finite-difference scheme with approximate or discrete transparent boundary conditions (TBCs), the Strang-type splitting with respect to the potential is applied. For the resulting method, the uniform in time L2-stability is proved. Due to the ...
Added: March 16, 2013
, , et al., Odessa Astronomical Publications 2013 No. 26/2 P. 311-315
Now a day the task of comparative analysis of sample sources from different astronomical catalogs is becoming very topical for astronomers. It is useful both for the same spectral range data and for their cross-analysis for different spectral ranges. It is also important, that the real data from astronomical instruments is placed on-line. For the ...
Added: December 11, 2013
, , et al., Physical Review E - Statistical, Nonlinear, and Soft Matter Physics 2021 Vol. 103 P. 053301-1-053301-24
Population annealing is a recent addition to the arsenal of the practitioner in computer simulations in statistical physics and it proves to deal well with systems with complex free-energy landscapes. Above all else, it promises to deliver unrivaled parallel scaling qualities, being suitable for parallel machines of the biggest caliber. Here we study population annealing using as the ...
Added: May 5, 2021
, , Odessa Astronomical Publications 2014 Vol. 27 No. 2 P. 72-73
Added: November 24, 2014
, , Чебышевский сборник 2012 Т. 13 № 2 (42) С. 91-105
We consider different parallel algortihms for operations in prime fields and their applications for operations on points of elliptic curves. The work provides results for implementations of these algorithms on NVIDIA graphical processors. ...
Added: February 25, 2013
Acceleration of plasma in current sheet during substorm dipolarizations in the Earth’s magnetotail: Comparison of different mechanisms
, , et al., Physics of Plasmas 2019 Vol. 26 No. 4 P. 1-9
This work is devoted to the investigation of particle acceleration during magnetospheric dipolarizations. A numerical model is presented taking into account the four scenarios of plasma acceleration that can be realized: (A) total dipolarization with characteristic time scales of 3 min; (B) single peak value of the normal magnetic component Bz occurring on the time scale of ...
Added: April 16, 2019
Пенза : ПГУ, 2016
В сборник трудов включены доклады ХХI-го Международного симпозиума «Надежность и качество», проходившего с 23 по 29 мая 2016 г. в городе Пензе. Рассмотрены актуальные проблемы теории и практики повышения надежности и качества; эффективности внедрения инновационных и информационных технологий в фундаментальных научных и прикладных исследованиях, образовательных и коммуникативных системах и средах, экономике и юриспруденции; методов и средств ...
Added: May 27, 2016
, , Computational Mathematics and Modeling 2016 Vol. 27 No. 2 P. 247-253
Added: December 22, 2016
«Суперкомпьютеры» vs «параллельное программирование». «Параллельное программирование» vs «совместная деятельность». Как изучать тему «параллельные вычисления» в средней школе?
, Современные информационные технологии и ИТ-образование 2015 Т. 1 № 11 С. 302-309
Статья посвящена вопросу о включении в школьный курс информатики темы «параллельные вычисления». Упоминается ряд возникающих при этом проблем, рассматривается цель изучения темы, отбор материала, некоторые предложения по методике обучения, механизмы апробации предложенной методики и накопленный опыт. Не затрагивается вопрос о месте этого материала в учебной программе. ...
Added: November 29, 2015
, Научно-технический вестник Поволжья 2016 № 5 С. 21-24
This paper presents an analysis of the results of numerical calculation of Bayes criterion of efficiency for some financial strategies used in the series of games with nature. A large number of simulations of the stochastic process of the game are formed with probability characteristics coincide with the similar values of the task. The results ...
Added: October 19, 2016
, , et al., Journal of High Energy Physics 2019 Vol. 1904 No. 77 P. 1-16
Heavy Neutral Leptons (HNLs) are hypothetical particles predicted by many extensions of the Standard Model. These particles can, among other things, explain the origin of neutrino masses, generate the observed matter-antimatter asymmetry in the Universe and provide a dark matter candidate. The SHiP experiment will be able to search for HNLs produced in decays of heavy ...
Added: April 30, 2019
ЗАЧЕМ и КАКИЕ нужны суперкомпьютеры эксафлопсного класса? Предсказательное моделирование свойств и многомасштабных процессов в материаловедении
, , et al., Программные системы: теория и приложения 2014 Т. 5 № 1 С. 191-244
Рассматривается подход, позволяющий выявить, для каких задач нужны суперкомпьютеры эксафлопсного класса. Возможности подхода рассмотрены на примерах актуальных задач материаловедения, физики конденсированного вещества и плотной плазмы, для решения которых необходимо атомистическое моделирование на современных и создаваемых в настоящее время суперкомпьютерах. Для каждой задачи проведено соответствие между набором изучаемых явлений и требуемым уровнем быстродействия (числа ядер) вычислительной ...
Added: August 11, 2014
, , et al., International Journal of Polymer Analysis and Characterization 2015 Vol. 20 No. 8 P. 743-753
This study produced for the first time composites based on polystyrene microspheres with nano-sized titania synthesized by sulfate and sol-gel methods, as well as commercial Hombifine N, Hombikat UV100, and Degussa P25. The starting samples and the composites were characterized by X-ray diffraction, scanning electron microscopy, nitrogen adsorption capacity, and IR spectroscopy. The highest content ...
Added: December 2, 2015