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

arXiv:cond-mat/9609138 · Physical review. B, Condensed matter

Critical behavior of the supersolid transition in Bose-Hubbard models

Erwin Frey, Leon Balents

Published Jan 1, 1997• 51 citations
arXivHTMLPDF

Abstract

We study the phase transitions of interacting bosons at zero temperature between superfluid (SF) and supersolid (SS) states. The latter are characterized by simultaneous off-diagonal long-range order and broken translational symmetry. The critical phenomena is described by a long-wavelength effective action, derived on symmetry grounds and verified by explicit calculation. We consider two types of supersolid ordering: checkerboard (X) and collinear (C), which are the simplest cases arising in two dimensions on a square lattice. We find that the SF-CSS transition is in the three-dimensional XY universality class. The SF-XSS transition exhibits nontrivial critical behavior, and appears, within a d=3-\ensuremath{\epsilon} expansion, to be driven generically first order by fluctuations. However, within a one-loop calculation directly in d=2 a strong-coupling fixed point with striking ``non-Bose-liquid'' behavior is found. At special isolated multicritical points of particle-hole symmetry, the system falls into the three-dimensional Ising universality class.

DOI
10.1103/physrevb.55.1050

Related papers

Thermal phase transitions of supersolids in the extended Bose-Hubbard model

Kwai-Kong Ng

We investigate numerically the finite-temperature phase diagrams of the extended Bose-Hubbard model in a two-dimensional square lattice. In particular, we focus on the melting of supersolid phases of two different crystal orderings, stripe and star orders, arising from the competition of the nearest- and next-nearest-neighbor interactions in the vicinity of quarter filling. The two crystal orders are the result of broken translational symmetry in either one or in both $x$ and $y$ directions. The broken gauge symmetry of the supersolids are found to be restored via a Kosterlitz-Thouless transition while the broken translational symmetries are restored via a single second-order phase transition, instead of two second-order transitions in the Ising universality class. On the other hand, the phase transitions between the star and stripe orders are first order in nature.

arXiv:1011.1068•Physical Review B•Nov 3, 2010•14 citations

Generic first-order versus continuous quantum nucleation of supersolidity

Lars Bonnes, Stefan Weßel

We analyze the nucleation of supersolid order out of the superfluid ground state of bosons on the triangular lattice. While the stability of supersolidity against phase separation in this system is by now well established for nearest-neighbor and long-range dipolar interactions, relevant for two-dimensional arrays of ultracold polar molecules, here we address directly the nature of the superfluid-to-supersolid transition. Based on symmetry arguments and quantum Monte Carlo simulations, we conclude that this quantum phase transition is driven first order beyond the line of particle-hole symmetry. Along this line, the transition is continuous and its scaling behavior is consistent with the three-dimensional (3D) XY universality class. We relate this finding to a 3D ${\mathbb{Z}}_{6}$ clock model description of the enlarged symmetry of the solid order parameter field. In the generic case, however, the symmetry reduces to that of a 3D ${\mathbb{Z}}_{3}$ clock model, which reflects the first-order nature of the generic superfluid-to-supersolid quantum phase transition on the triangular lattice.

arXiv:1105.1682•Physical Review B•Aug 5, 2011•30 citations

Phase diagrams, critical, and multicritical behavior of hard-core Bose-Hubbard models

Christian Pich, Erwin Frey

We determine the zero-temperature phase diagram of the hard-core Bose-Hubbard model on a square lattice by mean-field theory supplemented by a linear spin-wave analysis. Due to the interplay between nearest- and next-nearest-neighbor interaction and cubic anisotropy several supersolid phases with checkerboard, stripe-domain, or intermediate symmetry are stabilized. The phase diagrams show three different topologies depending on the relative strength of nearest- and next-nearest-neighbor interaction. We also find a rich variety of new quantum critical behavior and multicritical points and discuss the corresponding effective actions and universality classes.

arXiv:cond-mat/9711209•Physical review. B, Condensed matter•Jun 1, 1998•24 citations

Extended Bose-Hubbard model on a honeycomb lattice

Jing Yu Gan, Yu Chuan Wen, Jinwu Ye +3 more

We study the extended Bose-Hubbard model on a two-dimensional honeycomb lattice by using large-scale quantum Monte Carlo simulations. We present the ground-state phase diagrams for both the hard-core and the soft-core bosons. For the hard-core case, the transition between the $\ensuremath{\rho}=1∕2$ solid and the superfluid is first order, and the supersolid state is unstable toward phase separation. For the soft-core case, due to the presence of multiple occupation, a stable particle-induced supersolid (SS-p) phase emerges when $1∕2<\ensuremath{\rho}<1$. The transition from the solid at $\ensuremath{\rho}=1∕2$ to the SS-p phase is second order with the superfluid density scaling as ${\ensuremath{\rho}}_{s}\ensuremath{\sim}\ensuremath{\rho}\ensuremath{-}1∕2$. The SS-p phase has the same diagonal order as the solid at $\ensuremath{\rho}=1∕2$. As the chemical potential increases further, the SS-p phase turns into a solid where two bosons occupy each site of one sublattice through a first-order transition. We also calculate the critical exponents of the transition between the $\ensuremath{\rho}=1∕2$ solid and superfluid at the Heisenberg point for the hard-core case. We find the dynamical critical exponent $z=0.15$, which is smaller than results obtained on smaller lattices. This indicates that $z$ approaches zero in the thermodynamic limit, and thus the transition is also first order even at the Heisenberg point.

arXiv:cond-mat/0701120•Physical Review B•Jun 15, 2007•37 citations

Quantum phase transitions in disordered Bose systems

Ranjan Mukhopadhyay

NOTE: Text or symbols not renderable in plain ASCII are indicated by [...]. Abstract is included in .pdf document. We study the nature of various quantum phase transitions corresponding to the onset of superfluidity, at zero temperature, of bosons in a quenched medium. Particle-hole symmetry plays an essential role in determining the universality class of the transitions. To obtain a model with an exact particle-hole symmetry it is necessary to use the Josephson junction array Hamiltonian, which may include disorder in the Josephson couplings between phases at different sates. The functional integral formulation of this problem in d spatial dimensions yields a (d + 1)-dimensional classical XY-model with extended disorder, constant along the extra imaginary time dimension -- the so-called random rod problem. Particle-hole symmetry may then be broken by adding nonzero chemical potentials or site energies, which may also be site dependent and random. We may then distinguish three cases: (i) exact particle-hole symmetry, in which the site energies all vanish, (ii) statistical particle-hole symmetry in which the site energy distribution is symmetric about zero and hence vanishes on average, and (iii) complete absence of particle-hole symmetry in which the distribution is generic. We explore in each case the nature of the excitations in the nonsuperfluid Pose glass phase. We find, for example, that the compressibility, which has the interpretation of a temporal spin stiffness or superfluid density, is positive in cases (ii) and (iii), but that it vanishes with an essential singularity as full particle-hole symmetry is restored. We then focus on the critical point and discuss the validity of various scaling arguments. In particular, we argue that the dynamical exponent z could be different from d, and the arguments leading to their equality are incorrect. We then discuss the relevance of a type (ii) particle-hole symmetry breaking perturbation to the random rod critical behavior, identifying a nontrivial crossover exponent. This exponent cannot be calculated exactly but is argued to be positive and the symmetry breaking perturbation therefore relevant. We argue next that a perturbation of type (iii) is irrelevant to the resulting type (ii) critical behavior: the statistical symmetry is restored on large scales close to the critical point, and case (ii) therefore describes the dirty boson fixed point. Using various duality transformations we verify all of these ideas in one dimension. To study higher dimensions we attempt, with partial success, to generalize the Dorogovtsev-Cardy-Boyonovsky double epsilon expansion technique to this problem. We find that when the dimension of time [...] is sufficiently small the symmetry breaking perturbation of type (ii) is irrelevant, but that for sufficiently large [...] this is a relevant perturbation and a new stable commensurate fixed point appears. We speculate that this new fixed point becomes the dirty boson fixed point when [...] = 1. We also show that for [...], there exists a particle-hole asymmetric fixed point of type (iii), but we provide evidence that it merges with the commensurate fixed point for some finite [...]. This tends to confirm symmetry restoration at the physical [...] = 1.

•Feb 15, 2007•0 citations

Quantum critical behavior of ultracold bosons in the two-dimensional Bose-Hubbard lattice

Ulli Pohl, Sayak Ray, Johann Kroha

We investigate the temperature-dependent behavior emerging in the vicinity of the superfluid (SF) to Mott insulator (MI) transition of interacting bosons in a two-dimensional optical lattice, described by the Bose-Hubbard model. The equilibrium phase diagram at finite temperatures is computed by means of the cluster mean-field theory (CMF) where the effect of non-local correlations is analyzed systematically by finite-size scaling of the cluster size. The phase diagram exhibits a rich structure including a transition and a crossover of the SF and MI phases respectively to a normal fluid (NF) state at finite temperature. In order to characterize these phases, and the NF transition and crossover scales, we calculate, in addition to the condensate amplitude, the superfluid fraction, sound velocity and compressibility. The phase boundaries obtained by CMF with finite-size scaling agree quantitatively with quantum Monte Carlo (QMC) results as well as with experiments. The von Neumann entanglement entropy of a cluster exhibits critical enhancement near the SF-MI quantum critical point (QCP). We also discuss the behavior of the transition lines near this QCP at the particle-hole symmetric point located at the tip of a Mott lobe as well as away from particle-hole symmetry.

•Jun 28, 2021•0 citations
We study the phase transitions of interacting bosons at zero temperature between superfluid (SF) and supersolid (SS) states. The latter are characterized by simultaneous off-diagonal long-range order and broken translational symmetry. The critical phenomena is described by a long-wavelength effective action, derived on symmetry grounds and verified by explicit calculation. We consider two types of supersolid ordering: checkerboard (X) and collinear (C), which are the simplest cases arising in two dimensions on a square lattice. We find that the SF-CSS transition is in the three-dimensional XY universality class. The SF-XSS transition exhibits nontrivial critical behavior, and appears, within a d=3-\ensuremath{\epsilon} expansion, to be driven generically first order by fluctuations. However, within a one-loop calculation directly in d=2 a strong-coupling fixed point with striking ``non-Bose-liquid'' behavior is found. At special isolated multicritical points of particle-hole symmetry, the system falls into the three-dimensional Ising universality class.
We investigate numerically the finite-temperature phase diagrams of the extended Bose-Hubbard model in a two-dimensional square lattice. In particular, we focus on the melting of supersolid phases of two different crystal orderings, stripe and star orders, arising from the competition of the nearest- and next-nearest-neighbor interactions in the vicinity of quarter filling. The two crystal orders are the result of broken translational symmetry in either one or in both xxx and yyy directions. The broken gauge symmetry of the supersolids are found to be restored via a Kosterlitz-Thouless transition while the broken translational symmetries are restored via a single second-order phase transition, instead of two second-order transitions in the Ising universality class. On the other hand, the phase transitions between the star and stripe orders are first order in nature.
We analyze the nucleation of supersolid order out of the superfluid ground state of bosons on the triangular lattice. While the stability of supersolidity against phase separation in this system is by now well established for nearest-neighbor and long-range dipolar interactions, relevant for two-dimensional arrays of ultracold polar molecules, here we address directly the nature of the superfluid-to-supersolid transition. Based on symmetry arguments and quantum Monte Carlo simulations, we conclude that this quantum phase transition is driven first order beyond the line of particle-hole symmetry. Along this line, the transition is continuous and its scaling behavior is consistent with the three-dimensional (3D) XY universality class. We relate this finding to a 3D Z6{\mathbb{Z}}_{6}Z6​ clock model description of the enlarged symmetry of the solid order parameter field. In the generic case, however, the symmetry reduces to that of a 3D Z3{\mathbb{Z}}_{3}Z3​ clock model, which reflects the first-order nature of the generic superfluid-to-supersolid quantum phase transition on the triangular lattice.
We determine the zero-temperature phase diagram of the hard-core Bose-Hubbard model on a square lattice by mean-field theory supplemented by a linear spin-wave analysis. Due to the interplay between nearest- and next-nearest-neighbor interaction and cubic anisotropy several supersolid phases with checkerboard, stripe-domain, or intermediate symmetry are stabilized. The phase diagrams show three different topologies depending on the relative strength of nearest- and next-nearest-neighbor interaction. We also find a rich variety of new quantum critical behavior and multicritical points and discuss the corresponding effective actions and universality classes.
We study the extended Bose-Hubbard model on a two-dimensional honeycomb lattice by using large-scale quantum Monte Carlo simulations. We present the ground-state phase diagrams for both the hard-core and the soft-core bosons. For the hard-core case, the transition between the $\ensuremath{\rho}=1∕2$ solid and the superfluid is first order, and the supersolid state is unstable toward phase separation. For the soft-core case, due to the presence of multiple occupation, a stable particle-induced supersolid (SS-p) phase emerges when $1∕2<\ensuremath{\rho}<1$. The transition from the solid at $\ensuremath{\rho}=1∕2$ to the SS-p phase is second order with the superfluid density scaling as ${\ensuremath{\rho}}_{s}\ensuremath{\sim}\ensuremath{\rho}\ensuremath{-}1∕2$. The SS-p phase has the same diagonal order as the solid at $\ensuremath{\rho}=1∕2$. As the chemical potential increases further, the SS-p phase turns into a solid where two bosons occupy each site of one sublattice through a first-order transition. We also calculate the critical exponents of the transition between the $\ensuremath{\rho}=1∕2$ solid and superfluid at the Heisenberg point for the hard-core case. We find the dynamical critical exponent z=0.15z=0.15z=0.15, which is smaller than results obtained on smaller lattices. This indicates that zzz approaches zero in the thermodynamic limit, and thus the transition is also first order even at the Heisenberg point.
NOTE: Text or symbols not renderable in plain ASCII are indicated by [...]. Abstract is included in .pdf document. We study the nature of various quantum phase transitions corresponding to the onset of superfluidity, at zero temperature, of bosons in a quenched medium. Particle-hole symmetry plays an essential role in determining the universality class of the transitions. To obtain a model with an exact particle-hole symmetry it is necessary to use the Josephson junction array Hamiltonian, which may include disorder in the Josephson couplings between phases at different sates. The functional integral formulation of this problem in d spatial dimensions yields a (d + 1)-dimensional classical XY-model with extended disorder, constant along the extra imaginary time dimension -- the so-called random rod problem. Particle-hole symmetry may then be broken by adding nonzero chemical potentials or site energies, which may also be site dependent and random. We may then distinguish three cases: (i) exact particle-hole symmetry, in which the site energies all vanish, (ii) statistical particle-hole symmetry in which the site energy distribution is symmetric about zero and hence vanishes on average, and (iii) complete absence of particle-hole symmetry in which the distribution is generic. We explore in each case the nature of the excitations in the nonsuperfluid Pose glass phase. We find, for example, that the compressibility, which has the interpretation of a temporal spin stiffness or superfluid density, is positive in cases (ii) and (iii), but that it vanishes with an essential singularity as full particle-hole symmetry is restored. We then focus on the critical point and discuss the validity of various scaling arguments. In particular, we argue that the dynamical exponent z could be different from d, and the arguments leading to their equality are incorrect. We then discuss the relevance of a type (ii) particle-hole symmetry breaking perturbation to the random rod critical behavior, identifying a nontrivial crossover exponent. This exponent cannot be calculated exactly but is argued to be positive and the symmetry breaking perturbation therefore relevant. We argue next that a perturbation of type (iii) is irrelevant to the resulting type (ii) critical behavior: the statistical symmetry is restored on large scales close to the critical point, and case (ii) therefore describes the dirty boson fixed point. Using various duality transformations we verify all of these ideas in one dimension. To study higher dimensions we attempt, with partial success, to generalize the Dorogovtsev-Cardy-Boyonovsky double epsilon expansion technique to this problem. We find that when the dimension of time [...] is sufficiently small the symmetry breaking perturbation of type (ii) is irrelevant, but that for sufficiently large [...] this is a relevant perturbation and a new stable commensurate fixed point appears. We speculate that this new fixed point becomes the dirty boson fixed point when [...] = 1. We also show that for [...], there exists a particle-hole asymmetric fixed point of type (iii), but we provide evidence that it merges with the commensurate fixed point for some finite [...]. This tends to confirm symmetry restoration at the physical [...] = 1.
We investigate the temperature-dependent behavior emerging in the vicinity of the superfluid (SF) to Mott insulator (MI) transition of interacting bosons in a two-dimensional optical lattice, described by the Bose-Hubbard model. The equilibrium phase diagram at finite temperatures is computed by means of the cluster mean-field theory (CMF) where the effect of non-local correlations is analyzed systematically by finite-size scaling of the cluster size. The phase diagram exhibits a rich structure including a transition and a crossover of the SF and MI phases respectively to a normal fluid (NF) state at finite temperature. In order to characterize these phases, and the NF transition and crossover scales, we calculate, in addition to the condensate amplitude, the superfluid fraction, sound velocity and compressibility. The phase boundaries obtained by CMF with finite-size scaling agree quantitatively with quantum Monte Carlo (QMC) results as well as with experiments. The von Neumann entanglement entropy of a cluster exhibits critical enhancement near the SF-MI quantum critical point (QCP). We also discuss the behavior of the transition lines near this QCP at the particle-hole symmetric point located at the tip of a Mott lobe as well as away from particle-hole symmetry.