An efficient algorithm for numerical computations of continuous densities of states

Langfeld, K.; Lucini, Biagio; Pellegrini, R.; Rago, A.

doi:10.1140/epjc/s10052-016-4142-5

Journal article 945 views 145 downloads

An efficient algorithm for numerical computations of continuous densities of states

K. Langfeld, Biagio Lucini

, R. Pellegrini, A. Rago

The European Physical Journal C, Volume: 76, Issue: 6

Swansea University Author: Biagio Lucini

PDF | Version of Record

This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY).
Download (953.65KB)

Check full text

DOI (Published version): 10.1140/epjc/s10052-016-4142-5

Abstract

In Wang-Landau type algorithms, Monte-Carlo updates are performed with respect to the density of states, which is iteratively refined during simulations. The partition function and thermodynamic observables are then obtained by standard integration. In this work, our recently introduced method in th...

Full description

Published in:	The European Physical Journal C
ISSN:	1434-6044 1434-6052
Published:	Springer Science and Business Media LLC 2016
Online Access:	Check full text
URI:	https://cronfa.swan.ac.uk/Record/cronfa28023
Tags:	Add Tag No Tags, Be the first to tag this record!

first_indexed	2016-05-18T01:18:58Z
last_indexed	2020-08-04T02:44:43Z
id	cronfa28023
recordtype	SURis
fullrecord	<?xml version="1.0"?><rfc1807><datestamp>2020-08-03T14:12:27.8403269</datestamp><bib-version>v2</bib-version><id>28023</id><entry>2016-05-17</entry><title>An efficient algorithm for numerical computations of continuous densities of states</title><swanseaauthors><author><sid>7e6fcfe060e07a351090e2a8aba363cf</sid><ORCID>0000-0001-8974-8266</ORCID><firstname>Biagio</firstname><surname>Lucini</surname><name>Biagio Lucini</name><active>true</active><ethesisStudent>false</ethesisStudent></author></swanseaauthors><date>2016-05-17</date><deptcode>SMA</deptcode><abstract>In Wang-Landau type algorithms, Monte-Carlo updates are performed with respect to the density of states, which is iteratively refined during simulations. The partition function and thermodynamic observables are then obtained by standard integration. In this work, our recently introduced method in this class (the LLR approach) is analysed and further developed. Our approach is a histogram free method particularly suited for systems with continuous degrees of freedom giving rise to a continuum density of states, as it is commonly found in Lattice Gauge Theories and in some Statistical Mechanics systems. We show that the method possesses an exponential error suppression that allows us to estimate the density of states over several orders of magnitude with nearly-constant relative precision. We explain how ergodicity issues can be avoided and how expectation values of arbitrary observables can be obtained within this framework. We then demonstrate the method using Compact U(1) Lattice Gauge Theory as a show case. A thorough study of the algorithm parameter dependence of the results is performed and compared with the analytically expected behaviour. We obtain high precision values for the critical coupling for the phase transition and for the peak value of the specific heat for lattice sizes ranging from 84 to 204. Our results perfectly agree with the reference values reported in the literature, which covers lattice sizes up to 184. Robust results for the 204 volume are obtained for the first time. This latter investigation, which, due to strong metastabilities developed at the pseudo-critical coupling of the system, so far has been out of reach even on supercomputers with importance sampling approaches, has been performed to high accuracy with modest computational resources. This shows the potential of the method for studies of first order phase transitions. Other situations where the method is expected to be superior to importance sampling techniques are pointed out.</abstract><type>Journal Article</type><journal>The European Physical Journal C</journal><volume>76</volume><journalNumber>6</journalNumber><publisher>Springer Science and Business Media LLC</publisher><issnPrint>1434-6044</issnPrint><issnElectronic>1434-6052</issnElectronic><keywords>Wilson Loop; Order Phase Transition; Importance Sampling; Lattice Gauge Theory; Tunnelling Time</keywords><publishedDay>2</publishedDay><publishedMonth>6</publishedMonth><publishedYear>2016</publishedYear><publishedDate>2016-06-02</publishedDate><doi>10.1140/epjc/s10052-016-4142-5</doi><url/><notes>EPJC is an open-access journal funded by SCOAP3 and licensed under CC BY 4.0</notes><college>COLLEGE NANME</college><department>Mathematics</department><CollegeCode>COLLEGE CODE</CollegeCode><DepartmentCode>SMA</DepartmentCode><institution>Swansea University</institution><degreesponsorsfunders>SCOAP3</degreesponsorsfunders><apcterm/><lastEdited>2020-08-03T14:12:27.8403269</lastEdited><Created>2016-05-17T22:08:24.6576121</Created><path><level id="1">Faculty of Science and Engineering</level><level id="2">School of Biosciences, Geography and Physics - Physics</level></path><authors><author><firstname>K.</firstname><surname>Langfeld</surname><order>1</order></author><author><firstname>Biagio</firstname><surname>Lucini</surname><orcid>0000-0001-8974-8266</orcid><order>2</order></author><author><firstname>R.</firstname><surname>Pellegrini</surname><order>3</order></author><author><firstname>A.</firstname><surname>Rago</surname><order>4</order></author></authors><documents><document><filename>0028023-14102016123454.pdf</filename><originalFilename>Lucini2016.pdf</originalFilename><uploaded>2016-10-14T12:34:54.3930000</uploaded><type>Output</type><contentLength>1094309</contentLength><contentType>application/pdf</contentType><version>Version of Record</version><cronfaStatus>true</cronfaStatus><documentNotes>This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY).</documentNotes><copyrightCorrect>true</copyrightCorrect><language>eng</language><licence>http://creativecommons.org/licenses/by/4.0</licence></document></documents><OutputDurs/></rfc1807>
spelling	2020-08-03T14:12:27.8403269 v2 28023 2016-05-17 An efficient algorithm for numerical computations of continuous densities of states 7e6fcfe060e07a351090e2a8aba363cf 0000-0001-8974-8266 Biagio Lucini Biagio Lucini true false 2016-05-17 SMA In Wang-Landau type algorithms, Monte-Carlo updates are performed with respect to the density of states, which is iteratively refined during simulations. The partition function and thermodynamic observables are then obtained by standard integration. In this work, our recently introduced method in this class (the LLR approach) is analysed and further developed. Our approach is a histogram free method particularly suited for systems with continuous degrees of freedom giving rise to a continuum density of states, as it is commonly found in Lattice Gauge Theories and in some Statistical Mechanics systems. We show that the method possesses an exponential error suppression that allows us to estimate the density of states over several orders of magnitude with nearly-constant relative precision. We explain how ergodicity issues can be avoided and how expectation values of arbitrary observables can be obtained within this framework. We then demonstrate the method using Compact U(1) Lattice Gauge Theory as a show case. A thorough study of the algorithm parameter dependence of the results is performed and compared with the analytically expected behaviour. We obtain high precision values for the critical coupling for the phase transition and for the peak value of the specific heat for lattice sizes ranging from 84 to 204. Our results perfectly agree with the reference values reported in the literature, which covers lattice sizes up to 184. Robust results for the 204 volume are obtained for the first time. This latter investigation, which, due to strong metastabilities developed at the pseudo-critical coupling of the system, so far has been out of reach even on supercomputers with importance sampling approaches, has been performed to high accuracy with modest computational resources. This shows the potential of the method for studies of first order phase transitions. Other situations where the method is expected to be superior to importance sampling techniques are pointed out. Journal Article The European Physical Journal C 76 6 Springer Science and Business Media LLC 1434-6044 1434-6052 Wilson Loop; Order Phase Transition; Importance Sampling; Lattice Gauge Theory; Tunnelling Time 2 6 2016 2016-06-02 10.1140/epjc/s10052-016-4142-5 EPJC is an open-access journal funded by SCOAP3 and licensed under CC BY 4.0 COLLEGE NANME Mathematics COLLEGE CODE SMA Swansea University SCOAP3 2020-08-03T14:12:27.8403269 2016-05-17T22:08:24.6576121 Faculty of Science and Engineering School of Biosciences, Geography and Physics - Physics K. Langfeld 1 Biagio Lucini 0000-0001-8974-8266 2 R. Pellegrini 3 A. Rago 4 0028023-14102016123454.pdf Lucini2016.pdf 2016-10-14T12:34:54.3930000 Output 1094309 application/pdf Version of Record true This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY). true eng http://creativecommons.org/licenses/by/4.0
title	An efficient algorithm for numerical computations of continuous densities of states
spellingShingle	An efficient algorithm for numerical computations of continuous densities of states Biagio Lucini
title_short	An efficient algorithm for numerical computations of continuous densities of states
title_full	An efficient algorithm for numerical computations of continuous densities of states
title_fullStr	An efficient algorithm for numerical computations of continuous densities of states
title_full_unstemmed	An efficient algorithm for numerical computations of continuous densities of states
title_sort	An efficient algorithm for numerical computations of continuous densities of states
author_id_str_mv	7e6fcfe060e07a351090e2a8aba363cf
author_id_fullname_str_mv	7e6fcfe060e07a351090e2a8aba363cf_***_Biagio Lucini
author	Biagio Lucini
author2	K. Langfeld Biagio Lucini R. Pellegrini A. Rago
format	Journal article
container_title	The European Physical Journal C
container_volume	76
container_issue	6
publishDate	2016
institution	Swansea University
issn	1434-6044 1434-6052
doi_str_mv	10.1140/epjc/s10052-016-4142-5
publisher	Springer Science and Business Media LLC
college_str	Faculty of Science and Engineering
hierarchytype
hierarchy_top_id	facultyofscienceandengineering
hierarchy_top_title	Faculty of Science and Engineering
hierarchy_parent_id	facultyofscienceandengineering
hierarchy_parent_title	Faculty of Science and Engineering
department_str	School of Biosciences, Geography and Physics - Physics{{{_:::_}}}Faculty of Science and Engineering{{{_:::_}}}School of Biosciences, Geography and Physics - Physics
document_store_str	1
active_str	0
description	In Wang-Landau type algorithms, Monte-Carlo updates are performed with respect to the density of states, which is iteratively refined during simulations. The partition function and thermodynamic observables are then obtained by standard integration. In this work, our recently introduced method in this class (the LLR approach) is analysed and further developed. Our approach is a histogram free method particularly suited for systems with continuous degrees of freedom giving rise to a continuum density of states, as it is commonly found in Lattice Gauge Theories and in some Statistical Mechanics systems. We show that the method possesses an exponential error suppression that allows us to estimate the density of states over several orders of magnitude with nearly-constant relative precision. We explain how ergodicity issues can be avoided and how expectation values of arbitrary observables can be obtained within this framework. We then demonstrate the method using Compact U(1) Lattice Gauge Theory as a show case. A thorough study of the algorithm parameter dependence of the results is performed and compared with the analytically expected behaviour. We obtain high precision values for the critical coupling for the phase transition and for the peak value of the specific heat for lattice sizes ranging from 84 to 204. Our results perfectly agree with the reference values reported in the literature, which covers lattice sizes up to 184. Robust results for the 204 volume are obtained for the first time. This latter investigation, which, due to strong metastabilities developed at the pseudo-critical coupling of the system, so far has been out of reach even on supercomputers with importance sampling approaches, has been performed to high accuracy with modest computational resources. This shows the potential of the method for studies of first order phase transitions. Other situations where the method is expected to be superior to importance sampling techniques are pointed out.
published_date	2016-06-02T03:34:04Z
_version_	1763751435880628224
score	11.035634

An efficient algorithm for numerical computations of continuous densities of states

Similar Items