Reply to Comments on "Universal and Non-Universal Correction Terms of Bose Gases in Dilute Region: A Quantum Monte Carlo Study"
We reply to the comments on our previous paper J. Phys. Soc. Jpn. 91, 024001 (2022)
Discover
Research tools
Network
Opportunities
Account
Source author record
Akiko Masaki-Kato appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We reply to the comments on our previous paper J. Phys. Soc. Jpn. 91, 024001 (2022)
We study dilute gases of interacting bosons at zero-temperature in the region where the system is characterized only by the s-wave scattering length. We carry out quantum Monte Carlo simulation of the Bose Hubbard model and a continuous-space hard-core model. Fitting the extended Lee-Huang-Yang formula to the Monte Carlo results establishes the detailed mapping from the lattice model to the continuous field theory characterized only by the s-wave scattering length. Our estimate of the intrinsic s-wave scattering length $a_s$ of the Bose-Hubbard model is $a_s/a_l = 0.316(2)$ where $a_l$ is the lattice constant. It turned out that inclusion of the universal second correction term of $O(n \log n)$ makes the fitting worse while it is restored by inclusion of the non-universal third correction term, of which the existence was predicted analytically.
Based on the worm algorithm in the path-integral representation, we propose a general quantum Monte Carlo algorithm suitable for parallelizing on a distributed-memory computer by domain decomposition. Of particular importance is its application to large lattice systems of bosons and spins. A large number of worms are introduced and its population is controlled by a fictitious transverse field. For a benchmark, we study the size-dependence of the Bose-condensation order parameter of the hardcore Bose-Hubbard model with $L\times L\times βt = 10240\times 10240\times 16$, using 3200 computing cores, which shows good parallelization efficiency.