Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review A

Monte Carlo simulation of lattice bosons in three dimensions

Allan S. Blaer, Jining Han

Published Sep 1, 1992• 5 citations
Publisher

Abstract

We present an algorithm for calculating the thermodynamic properties of a system of nonrelativistic bosons on a three-dimensional spatial lattice. The method, which maps the three-dimensional quantum system onto a four-dimensional classical system, uses Monte Carlo sampling of configurations in either the canonical or the grand canonical ensemble. Our procedure is applicable to any system of lattice bosons with arbitrary short-range interactions. We test the algorithm by computing the temperature dependence of the energy, the heat capacity, and the condensate fraction of the free Bose gas.

DOI
10.1103/physreva.46.3225

Related papers

Quantum Monte Carlo algorithm for softcore boson systems

Jurij Šmakov, Kenji Harada, Naoki Kawashima

An efficient quantum Monte Carlo algorithm for the simulation of bosonic systems on a lattice in a grand canonical ensemble is proposed. It is based on the mapping of bosonic models to the spin models in the limit of the infinite total spin quantum number. It is demonstrated how this limit may be taken explicitly in the algorithm, eliminating the systematic errors. The efficiency of the algorithm is examined for the noninteracting lattice boson model and compared with the stochastic series expansion method with the heat-bath-type scattering probability of the random walker.

arXiv:cond-mat/0301416•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 21, 2003•3 citations

Monte Carlo simulation of hard-core bosons on a three-dimensional lattice

Morten Holm Pedersen, T. Schneider

A model of hard-core bosons on a three-dimensional (3D) lattice has been studied using both quantum Monte Carlo and classical Monte Carlo techniques. We have calculated the transition temperature for Bose-Einstein condensation, the zero-temperature penetration depth and the specific-heat critical amplitude, and have explored the hyperuniversal scaling relationship between these quantities pertaining to the expected 3D XY critical point. As a conclusion we draw a phenomenological analogy between the model studied and extreme type-II superconductors. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Mar 1, 1996•5 citations

Finite-size effects in canonical and grand-canonical quantum Monte Carlo simulations for fermions

Zhenjiu Wang, Fakher F. Assaad, Francesco Parisen Toldin

We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.

arXiv:1706.01874•Physical review. E•Oct 16, 2017•24 citations

Path integral Monte Carlo on a lattice. II. Bound states

Mark O'Callaghan, Bruce N. Miller

The equilibrium properties of a single quantum particle (qp) interacting with a classical gas for a wide range of temperatures that explore the system's behavior in the classical as well as in the quantum regime is investigated. Both the qp and the atoms are restricted to sites on a one-dimensional lattice. A path integral formalism developed within the context of the canonical ensemble is utilized, where the qp is represented by a closed, variable-step random walk on the lattice. Monte Carlo methods are employed to determine the system's properties. To test the usefulness of the path integral formalism, the Metropolis algorithm is employed to determine the equilibrium properties of the qp in the context of a square well potential, forcing the qp to occupy bound states. We consider a one-dimensional square well potential where all atoms on the lattice are occupied with one atom with an on-site potential except for a contiguous set of sites of various lengths centered at the middle of the lattice. Comparison of the potential energy, the energy fluctuations, and the correlation function are made between the results of the Monte Carlo simulations and the numerical calculations.

Physical review. E•Jul 18, 2016•0 citations

Thermodynamics of the noninteracting Bose gas in a two-dimensional box

Heqiu Li, Qiujiang Guo, Ji Jiang +1 more

Bose-Einstein condensation (BEC) of a noninteracting Bose gas of N particles in a two-dimensional box with Dirichlet boundary conditions is studied. Confirming previous work, we find that BEC occurs at finite N at low temperatures T without the occurrence of a phase transition. The conventionally-defined transition temperature T(E) for an infinite three-dimensional (3D) system is shown to correspond in a 2D system with finite N to a crossover temperature between a slow and rapid increase in the fractional boson occupation N(0)/N of the ground state with decreasing T. We further show that T(E)∼1/logN at fixed area per boson, so in the thermodynamic limit there is no significant BEC in 2D at finite T. Thus, paradoxically, BEC only occurs in 2D at finite N with no phase transition associated with it. Calculations of thermodynamic properties versus T and area A are presented, including Helmholtz free energy, entropy S, pressure p, ratio of p to the energy density U/A, heat capacity at constant volume (area) C(V) and at constant pressure C(p), isothermal compressibility κ(T) and thermal expansion coefficient α(p), obtained using both the grand-canonical ensemble (GCE) and canonical ensemble (CE) formalisms. The GCE formalism gives acceptable predictions for S, p, p/(U/A), κ(T) and α(p) at large N, T and A but fails for smaller values of these three parameters for which BEC becomes significant, whereas the CE formalism gives accurate results for all thermodynamic properties of finite systems even at low T and/or A where BEC occurs.

arXiv:1508.07656•Physical Review E•Dec 4, 2015•6 citations

Finite-Temperature Monte Carlo Calculations for Systems with Fermions

Shiwei Zhang

We present a quantum Monte Carlo method which allows calculations on many-fermion systems at finite temperatures without any sign decay. This enables simulations of the grand-canonical ensemble at large system sizes and low temperatures. Both diagonal and off-diagonal expectations can be computed straightforwardly. The sign decay is eliminated by a constraint on the fermion determinant. The algorithm is approximate. Tests on the Hubbard model show that accurate results on the energy and correlation functions can be obtained.

arXiv:cond-mat/9909091•Physical Review Letters•Oct 4, 1999•71 citations
We present an algorithm for calculating the thermodynamic properties of a system of nonrelativistic bosons on a three-dimensional spatial lattice. The method, which maps the three-dimensional quantum system onto a four-dimensional classical system, uses Monte Carlo sampling of configurations in either the canonical or the grand canonical ensemble. Our procedure is applicable to any system of lattice bosons with arbitrary short-range interactions. We test the algorithm by computing the temperature dependence of the energy, the heat capacity, and the condensate fraction of the free Bose gas.
An efficient quantum Monte Carlo algorithm for the simulation of bosonic systems on a lattice in a grand canonical ensemble is proposed. It is based on the mapping of bosonic models to the spin models in the limit of the infinite total spin quantum number. It is demonstrated how this limit may be taken explicitly in the algorithm, eliminating the systematic errors. The efficiency of the algorithm is examined for the noninteracting lattice boson model and compared with the stochastic series expansion method with the heat-bath-type scattering probability of the random walker.
A model of hard-core bosons on a three-dimensional (3D) lattice has been studied using both quantum Monte Carlo and classical Monte Carlo techniques. We have calculated the transition temperature for Bose-Einstein condensation, the zero-temperature penetration depth and the specific-heat critical amplitude, and have explored the hyperuniversal scaling relationship between these quantities pertaining to the expected 3D XY critical point. As a conclusion we draw a phenomenological analogy between the model studied and extreme type-II superconductors. \textcopyright{} 1996 The American Physical Society.
We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.
The equilibrium properties of a single quantum particle (qp) interacting with a classical gas for a wide range of temperatures that explore the system's behavior in the classical as well as in the quantum regime is investigated. Both the qp and the atoms are restricted to sites on a one-dimensional lattice. A path integral formalism developed within the context of the canonical ensemble is utilized, where the qp is represented by a closed, variable-step random walk on the lattice. Monte Carlo methods are employed to determine the system's properties. To test the usefulness of the path integral formalism, the Metropolis algorithm is employed to determine the equilibrium properties of the qp in the context of a square well potential, forcing the qp to occupy bound states. We consider a one-dimensional square well potential where all atoms on the lattice are occupied with one atom with an on-site potential except for a contiguous set of sites of various lengths centered at the middle of the lattice. Comparison of the potential energy, the energy fluctuations, and the correlation function are made between the results of the Monte Carlo simulations and the numerical calculations.
Bose-Einstein condensation (BEC) of a noninteracting Bose gas of N particles in a two-dimensional box with Dirichlet boundary conditions is studied. Confirming previous work, we find that BEC occurs at finite N at low temperatures T without the occurrence of a phase transition. The conventionally-defined transition temperature T(E) for an infinite three-dimensional (3D) system is shown to correspond in a 2D system with finite N to a crossover temperature between a slow and rapid increase in the fractional boson occupation N(0)/N of the ground state with decreasing T. We further show that T(E)∼1/logN at fixed area per boson, so in the thermodynamic limit there is no significant BEC in 2D at finite T. Thus, paradoxically, BEC only occurs in 2D at finite N with no phase transition associated with it. Calculations of thermodynamic properties versus T and area A are presented, including Helmholtz free energy, entropy S, pressure p, ratio of p to the energy density U/A, heat capacity at constant volume (area) C(V) and at constant pressure C(p), isothermal compressibility κ(T) and thermal expansion coefficient α(p), obtained using both the grand-canonical ensemble (GCE) and canonical ensemble (CE) formalisms. The GCE formalism gives acceptable predictions for S, p, p/(U/A), κ(T) and α(p) at large N, T and A but fails for smaller values of these three parameters for which BEC becomes significant, whereas the CE formalism gives accurate results for all thermodynamic properties of finite systems even at low T and/or A where BEC occurs.
We present a quantum Monte Carlo method which allows calculations on many-fermion systems at finite temperatures without any sign decay. This enables simulations of the grand-canonical ensemble at large system sizes and low temperatures. Both diagonal and off-diagonal expectations can be computed straightforwardly. The sign decay is eliminated by a constraint on the fermion determinant. The algorithm is approximate. Tests on the Hubbard model show that accurate results on the energy and correlation functions can be obtained.