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:1203.5629 · Physical Review B

Quantum and thermal transitions out of the pair-supersolid phase of two-species bosons on a square lattice

Chia-Min Chung, Shiang Fang, Pochung Chen

Published Jun 13, 2012• 14 citations
arXivHTMLPDF

Abstract

We investigate two-species bosons in a two-dimensional square lattice by the quantum Monte Carlo method. We show that the interspecies attraction and nearest-neighbor intraspecies repulsion result in the pair-supersolid phase, where a diagonal solid order coexists with an off-diagonal pair-superfluid order. The quantum and thermal transitions out of the pair-supersolid phase are characterized. It is found that there is a direct first-order transition from the pair-supersolid phase to the double-superfluid phase without an intermediate region. Furthermore, the melting of the pair supersolid occurs in two steps. Upon heating, first the pair superfluid is destroyed via a Kosterlitz-Thouless transition, then the solid order melts via an Ising transition.

DOI
10.1103/physrevb.85.214513

Related papers

Quantum phase transitions of two-species bosons

邱俊嘉

We investigate the phase of two-species bosons in square lattice by the quantum Monte Carlo simulations. The Mott insulator phase becomes solid by increasing the nearest-neighbor repulsion. The phase changing can be driven by increasing the strength of hopping. It is found that there is a direct first-order transition. There is also the transition form the superfluid phase to the supersolid phase by increasing the nearest-neighbor repulsion. And the transition can be found when we drive the strength of hopping. We found that there is a second-order transition.

•Jan 1, 2014•0 citations

Quantum Monte Carlo Study of Two Component Bosons and Entanglement

鍾佳民

This thesis is divided into two parts. In the first part we study quantum phases in bosonic systems in an optical lattice by using numerical unbiased Quantum Monte Carlo (QMC) method. Ultracold atoms experiments in optical lattices provide a possibility to simulate the materials of condensed matter in a clean and well-controled way. We focus on two components bosonic systems where exotic phases such as pair-superfluid and pair-supersolid appear. In the second part we provide numerical methods based on QMC simulations to study the entanglement properties of strong correlated systems. By employing the replica trick, we reconstruct the entanglement spectrum from the Renyi entropies. Furthermore, we study the trace of the power of the partial transposed reduced density matrix, which is related to the entanglement measurement negativity for mixed states. In the following we briefly summarize each chapter. In Chapter 1 we very briefly introduce the cold atom systems and the connection with QMC simulations. In Chapter 2 we review the QMC method based on the path integral (world line) representation. We focus on the directed worm algorithm which is an efficient algorithm with global updtae for bosonic and spin systems. All the following methods in this thesis – the method of two component bosonic systems, the reconstruction of entanglement spectrum and the measuring of partial transposed quantities, are based on directed worm algorithm introduced in this chapter. In Chapter 3 we discuss two component bosons in a square lattice. We show that the interspecies attraction and nearest-neighbor intraspecies repulsion result in the pair-supersolid phase, where a diagonal solid order coexists with an off-diagonal pair-superfluid order. The quantum and thermal transitions out of the pair-supersolid phase are characterized. It is found that there is a direct first-order transition from the pair-supersolid phase to the double-superfluid phase without an intermediate region. Furthermore, the melting of the pair-supersolid occurs in two steps. Upon heating, first the pair-superfluid is destroyed via a Kosterlitz-Thouless transition, then the solid order melts via an Ising transition. In Chapter 4 we represent a new method to reconstruct a subset of the entanglement spectrum of quantum many body systems by QMC, where the method can in principle be applied to two or higher dimension. The approach builds on the replica trick to evaluate particle number resolved traces of the first n of powers of a reduced density matrix. From this information we reconstruct first n entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector. In Chapter 5 we devise a Quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. In particular, we study several scale invariant combinations of the moments for the 1D hard-core boson model. For two adjacent intervals unusual finite size corrections are present, showing parity effects that oscillate with a filling dependent period. For large chains we find perfect agreement with conformal field theory (CFT) calculations. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. The finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderate large sizes.

•Jan 1, 2014•0 citations

Pair supersolid with atom-pair hopping on the state-dependent triangular lattice

Wanzhou Zhang, Ruoxi Yin, Yan-Cheng Wang

We systematically study an extended Bose-Hubbard model with atom hopping and atom-pair hopping in the presence of a three-body constraint on the triangular lattice. By means of large-scale quantum Monte Carlo simulations, the ground-state phase diagram is studied. We find a first-order transition between the atomic superfluid phase and the pair superfluid phase when the ratio of the atomic hopping and the atom-pair hopping is adapted. The first-order transition remains unchanged under various conditions. We then focus on the interplay among the atom-pair hopping, the on-site repulsion, and the nearest-neighbor repulsion. With on-site repulsion present, we observe first-order transitions between the Mott insulators and pair superfluid driven by the pair hopping. With the nearest-neighbor repulsion turning on, three typical solid phases with $2/3$, 1, and $4/3$ filling emerge at small atom-pair hopping region. A stable pair supersolid phase is found at small on-site repulsion. This is due to the three-body constraint and the pair hopping, which essentially make the model a quasihardcore boson system. Thus the pair supersolid state emerges basing on the order-by-disorder mechanism, by which hardcore bosons avoid classical frustration on the triangular lattice. Without on-site repulsion, the transitions between the pair supersolid and the atom superfluid or pair superfluid are first order, except for the particle-hole symmetric point. With weak on-site repulsion and atom hopping turning on, the transition between the pair supersolid and pair superfluid phase becomes continuous. The transition between solid and pair supersolid is three-dimensional $XY$ university, with dynamical exponent $z=1$ and correlation exponent $\ensuremath{\nu}=0.671\phantom{\rule{0.16em}{0ex}}55$. The thermal melting of pair supersolid belongs to the two-dimensional Ising university. We check both energetic and mechanical balance of pair supersolid phase. Lowering the three-body constraint, no pair supersolid is found due to the absence of degeneracy of pair solids in classical limits. We describe the experimental realization of pair tunneling on state dependent lattice.

arXiv:1208.6506•Physical Review B•Nov 19, 2013•17 citations

Ground-state phase diagram of spin- 12 bosons in a two-dimensional optical lattice

Laurent de Forges de Parny, F. Hébert, V. G. Rousseau +2 more

We study a two-species bosonic Hubbard model on a two-dimensional square lattice by means of quantum Monte Carlo simulations. In addition to the usual contact repulsive interactions between the particles, the Hamiltonian has an interconversion term which allows the transformation of two particles from one species to the other. The phases are characterized by their solid or superfluid properties and by their polarization, i.e., the difference in the populations. When interspecies interactions are smaller than the intraspecies ones, the system is unpolarized, whereas in the opposite case the system is unpolarized in even Mott insulator lobes and polarized in odd Mott lobes and also in the superfluid phase. We show that in the latter case the transition between the Mott insulator of total density 2 and the superfluid can be of either second or first order depending on the relative values of the interactions, whereas the transitions are continuous in all other cases.

arXiv:1105.6141•Physical Review B•Aug 26, 2011•20 citations

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

Two-species hard-core bosons on the triangular lattice: A quantum Monte Carlo study

Jian-Ping Lv, Qing‐Hu Chen, Youjin Deng

Using worm-type quantum Monte Carlo simulations, we investigate bosonic mixtures on a triangular lattice of two species of bosons, which interact via nearest-neighbor intraspecies ($V$) and on-site interspecies ($U$) repulsions. For the case of symmetric hopping amplitude (${t}_{A}/V={t}_{B}/V$) and $U/V=1$, we determine a rich ground-state phase diagram that contains double-solid, double-superfluid, supersolid (SS), solid-superfluid (solid-SF), and counterflow supersolid (CSS) states. The SS, solid-SF, and CSS states exhibit spontaneous symmetry breaking among the three sublattices of the triangular lattice and between the two species, which leads to a nonzero crystalline density wave order in each species. We, furthermore, show that the CSS and the SS states are present for ${t}_{A}/V\ensuremath{\ne}{t}_{B}/V$, and the latter even survives up to the ${t}_{A}/V\ensuremath{\rightarrow}\ensuremath{\infty}$ or ${t}_{B}/V\ensuremath{\rightarrow}\ensuremath{\infty}$ limit. The effects induced by the variation of $U/V$ and by the imbalance of particle numbers of the two species are also explored.

arXiv:1311.1604•Physical Review A•Jan 30, 2014•18 citations
We investigate two-species bosons in a two-dimensional square lattice by the quantum Monte Carlo method. We show that the interspecies attraction and nearest-neighbor intraspecies repulsion result in the pair-supersolid phase, where a diagonal solid order coexists with an off-diagonal pair-superfluid order. The quantum and thermal transitions out of the pair-supersolid phase are characterized. It is found that there is a direct first-order transition from the pair-supersolid phase to the double-superfluid phase without an intermediate region. Furthermore, the melting of the pair supersolid occurs in two steps. Upon heating, first the pair superfluid is destroyed via a Kosterlitz-Thouless transition, then the solid order melts via an Ising transition.
We investigate the phase of two-species bosons in square lattice by the quantum Monte Carlo simulations. The Mott insulator phase becomes solid by increasing the nearest-neighbor repulsion. The phase changing can be driven by increasing the strength of hopping. It is found that there is a direct first-order transition. There is also the transition form the superfluid phase to the supersolid phase by increasing the nearest-neighbor repulsion. And the transition can be found when we drive the strength of hopping. We found that there is a second-order transition.
This thesis is divided into two parts. In the first part we study quantum phases in bosonic systems in an optical lattice by using numerical unbiased Quantum Monte Carlo (QMC) method. Ultracold atoms experiments in optical lattices provide a possibility to simulate the materials of condensed matter in a clean and well-controled way. We focus on two components bosonic systems where exotic phases such as pair-superfluid and pair-supersolid appear. In the second part we provide numerical methods based on QMC simulations to study the entanglement properties of strong correlated systems. By employing the replica trick, we reconstruct the entanglement spectrum from the Renyi entropies. Furthermore, we study the trace of the power of the partial transposed reduced density matrix, which is related to the entanglement measurement negativity for mixed states. In the following we briefly summarize each chapter. In Chapter 1 we very briefly introduce the cold atom systems and the connection with QMC simulations. In Chapter 2 we review the QMC method based on the path integral (world line) representation. We focus on the directed worm algorithm which is an efficient algorithm with global updtae for bosonic and spin systems. All the following methods in this thesis – the method of two component bosonic systems, the reconstruction of entanglement spectrum and the measuring of partial transposed quantities, are based on directed worm algorithm introduced in this chapter. In Chapter 3 we discuss two component bosons in a square lattice. We show that the interspecies attraction and nearest-neighbor intraspecies repulsion result in the pair-supersolid phase, where a diagonal solid order coexists with an off-diagonal pair-superfluid order. The quantum and thermal transitions out of the pair-supersolid phase are characterized. It is found that there is a direct first-order transition from the pair-supersolid phase to the double-superfluid phase without an intermediate region. Furthermore, the melting of the pair-supersolid occurs in two steps. Upon heating, first the pair-superfluid is destroyed via a Kosterlitz-Thouless transition, then the solid order melts via an Ising transition. In Chapter 4 we represent a new method to reconstruct a subset of the entanglement spectrum of quantum many body systems by QMC, where the method can in principle be applied to two or higher dimension. The approach builds on the replica trick to evaluate particle number resolved traces of the first n of powers of a reduced density matrix. From this information we reconstruct first n entanglement spectrum levels using a polynomial root solver. We illustrate the power and limitations of the method by an application to the extended Bose-Hubbard model in one dimension where we are able to resolve the quasidegeneracy of the entanglement spectrum in the Haldane-insulator phase. In general, the method is able to reconstruct the largest few eigenvalues in each symmetry sector. In Chapter 5 we devise a Quantum Monte Carlo (QMC) method to calculate the moments of the partially transposed reduced density matrix at finite temperature. These are used to construct scale invariant combinations that are related to the negativity, a true measure of entanglement for two intervals embedded in a chain. In particular, we study several scale invariant combinations of the moments for the 1D hard-core boson model. For two adjacent intervals unusual finite size corrections are present, showing parity effects that oscillate with a filling dependent period. For large chains we find perfect agreement with conformal field theory (CFT) calculations. Oppositely, for disjoint intervals corrections are more severe and CFT is recovered only asymptotically. Furthermore, we provide evidence that their exponent is the same as that governing the corrections of the mutual information. Additionally we study the 1D Bose-Hubbard model in the superfluid phase. The finite-size effects are smaller and QMC data are already in impressive agreement with CFT at moderate large sizes.
We systematically study an extended Bose-Hubbard model with atom hopping and atom-pair hopping in the presence of a three-body constraint on the triangular lattice. By means of large-scale quantum Monte Carlo simulations, the ground-state phase diagram is studied. We find a first-order transition between the atomic superfluid phase and the pair superfluid phase when the ratio of the atomic hopping and the atom-pair hopping is adapted. The first-order transition remains unchanged under various conditions. We then focus on the interplay among the atom-pair hopping, the on-site repulsion, and the nearest-neighbor repulsion. With on-site repulsion present, we observe first-order transitions between the Mott insulators and pair superfluid driven by the pair hopping. With the nearest-neighbor repulsion turning on, three typical solid phases with 2/32/32/3, 1, and 4/34/34/3 filling emerge at small atom-pair hopping region. A stable pair supersolid phase is found at small on-site repulsion. This is due to the three-body constraint and the pair hopping, which essentially make the model a quasihardcore boson system. Thus the pair supersolid state emerges basing on the order-by-disorder mechanism, by which hardcore bosons avoid classical frustration on the triangular lattice. Without on-site repulsion, the transitions between the pair supersolid and the atom superfluid or pair superfluid are first order, except for the particle-hole symmetric point. With weak on-site repulsion and atom hopping turning on, the transition between the pair supersolid and pair superfluid phase becomes continuous. The transition between solid and pair supersolid is three-dimensional XYXYXY university, with dynamical exponent z=1z=1z=1 and correlation exponent $\ensuremath{\nu}=0.671\phantom{\rule{0.16em}{0ex}}55$. The thermal melting of pair supersolid belongs to the two-dimensional Ising university. We check both energetic and mechanical balance of pair supersolid phase. Lowering the three-body constraint, no pair supersolid is found due to the absence of degeneracy of pair solids in classical limits. We describe the experimental realization of pair tunneling on state dependent lattice.
We study a two-species bosonic Hubbard model on a two-dimensional square lattice by means of quantum Monte Carlo simulations. In addition to the usual contact repulsive interactions between the particles, the Hamiltonian has an interconversion term which allows the transformation of two particles from one species to the other. The phases are characterized by their solid or superfluid properties and by their polarization, i.e., the difference in the populations. When interspecies interactions are smaller than the intraspecies ones, the system is unpolarized, whereas in the opposite case the system is unpolarized in even Mott insulator lobes and polarized in odd Mott lobes and also in the superfluid phase. We show that in the latter case the transition between the Mott insulator of total density 2 and the superfluid can be of either second or first order depending on the relative values of the interactions, whereas the transitions are continuous in all other cases.
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.
Using worm-type quantum Monte Carlo simulations, we investigate bosonic mixtures on a triangular lattice of two species of bosons, which interact via nearest-neighbor intraspecies (VVV) and on-site interspecies (UUU) repulsions. For the case of symmetric hopping amplitude (tA/V=tB/V{t}_{A}/V={t}_{B}/VtA​/V=tB​/V) and U/V=1U/V=1U/V=1, we determine a rich ground-state phase diagram that contains double-solid, double-superfluid, supersolid (SS), solid-superfluid (solid-SF), and counterflow supersolid (CSS) states. The SS, solid-SF, and CSS states exhibit spontaneous symmetry breaking among the three sublattices of the triangular lattice and between the two species, which leads to a nonzero crystalline density wave order in each species. We, furthermore, show that the CSS and the SS states are present for ${t}_{A}/V\ensuremath{\ne}{t}_{B}/V$, and the latter even survives up to the ${t}_{A}/V\ensuremath{\rightarrow}\ensuremath{\infty}$ or ${t}_{B}/V\ensuremath{\rightarrow}\ensuremath{\infty}$ limit. The effects induced by the variation of U/VU/VU/V and by the imbalance of particle numbers of the two species are also explored.