title: Metropolis-Hastings thermal state sampling for numerical simulations of Bose-Einstein condensates.
abstract: We demonstrate the application of the Metropolis-Hastings algorithm to sampling of classical thermal states of one-dimensional Bose-Einstein quasicondensates in the classical fields approximation, both in untrapped and harmonically trapped case. The presented algorithm can be easily generalized to higher dimensions and arbitrary trap geometry. For truncated Wigner simulations the quantum noise can be added with conventional methods (half a quantum of energy in every mode). The advantage of the presented method over the usual analytical and stochastic ones lies in its ability to sample not only from canonical and grand canonical distributions, but also from the generalized Gibbs ensemble, which can help to shed new light on thermodynamics of integrable systems.
authors.affiliation: Vienna Center for Quantum Science and Technology, Atominstitut TU Wien, 1020 Vienna, Austria.
authors.forename: Pjotrs
authors.initials: P
authors.lastname: Grišins
authors.affiliation: Vienna Center for Quantum Science and Technology, Atominstitut TU Wien, 1020 Vienna, Austria ; Ioffe Physico-Technical Institute, Russian Academy of Sciences, 194021 St. Petersburg, Russia.
authors.forename: Igor E
authors.initials: IE
authors.lastname: Mazets
keywords: 1D; BEC; Gross–Pitaevskii equation; Metropolis–Hastings algorithm; Monte-Carlo simulations
references.citation: Phys Rev Lett. 2007 Oct 12;99(15):150405
references.pmid: 17995147
references.citation: Science. 2012 Sep 14;337(6100):1318-22
references.pmid: 22956685
references.citation: Phys Rev Lett. 2013 Jan 18;110(3):033001
references.pmid: 23373918
references.citation: Nat Commun. 2013;4:2077
references.pmid: 23804159
references.citation: Science. 2012 Mar 2;335(6072):1070-2
references.pmid: 22345397
references.citation: Nature. 2011 Feb 10;470(7333):236-9
references.pmid: 21270797
references.citation: Phys Rev Lett. 2007 Feb 2;98(5):050405
references.pmid: 17358832