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:1805.05042 · Physical Review A

Out-of-equilibrium dynamics of multiple second-order quantum phase transitions in an extended Bose-Hubbard model: Superfluid, supersolid, and density wave

Keita Shimizu, Takahiro Hirano, Jonghoon Park, Yoshihito Kuno, Ikuo Ichinose

Published Dec 4, 2018• 11 citations
arXivHTMLPDF

Abstract

In this paper, we study the dynamics of the Bose-Hubbard model with the nearest-neighbor repulsion by using time-dependent Gutzwiller methods. Near the unit filling, the phase diagram of the model contains density wave (DW), supersolid (SS), and superfluid (SF). The three phases are separated by two second-order phase transitions. We study ``slow-quench'' dynamics by varying the hopping parameter in the Hamiltonian as a function of time. In the phase transitions from the DW to SS and from the DW to SF, we focus on how the SF order forms and study scaling laws of the SF correlation length, vortex density, etc. The results are compared with the Kibble-Zurek scaling. On the other hand from the SF to DW, we study how the DW order evolves with generation of the domain walls and vortices. Measurement of first-order SF coherence reveals interesting behavior in the DW regime.

DOI
10.1103/physreva.98.063603

Related papers

Dynamics of first-order quantum phase transitions in extended Bose–Hubbard model: from density wave to superfluid and vice versa

Keita Shimizu, Takahiro Hirano, Jonghoon Park +2 more

In this paper, we study the nonequilibrium dynamics of the Bose–Hubbard model with the nearest-neighbor repulsion by using time-dependent Gutzwiller (GW) methods. In particular, we vary the hopping parameters in the Hamiltonian as a function of time, and investigate the dynamics of the system from the density wave (DW) to the superfluid (SF) crossing a first-order phase transition and vice versa. From the DW to SF, we find scaling laws for the correlation length and vortex density with respect to the quench time. This is a reminiscence of the Kibble–Zurek scaling for continuous phase transitions and contradicts the common expectation. We give a possible explanation for this observation. On the other hand from SF to DW, the system evolution depends on the initial SF state. When the initial state is the ground state obtained by the static GW methods, a coexisting state of the SF and DW domains forms after passing through the critical point. Coherence of the SF order parameter is lost as the system evolves. This is a phenomenon similar to the glass transition in classical systems. When the state starts from the SF with small local phase fluctuations, the system obtains a large size DW domain structure with thin domain walls.

arXiv:1803.02548•New Journal of Physics•Jul 26, 2018•20 citations

Dynamics of a quantum phase transition in the Bose-Hubbard model: Kibble-Zurek mechanism and beyond

Keita Shimizu, Yoshihito Kuno, Takahiro Hirano +1 more

In this paper, we study the dynamics of the Bose-Hubbard model by using time-dependent Gutzwiller methods. In particular, we vary the parameters in the Hamiltonian as a function of time, and investigate the temporal behavior of the system from the Mott insulator to the superfluid (SF) crossing a second-order phase transition. We first solve a time-dependent Schr\"odinger equation for the experimental setup recently done by Braun et al. [Proc. Natl. Acad. Sci. USA 112, 3641 (2015)] and show that the numerical and experimental results are in fairly good agreement. However, these results disagree with the Kibble-Zurek scaling. From our numerical study, we reveal a possible source of the discrepancy. Next, we calculate the critical exponents of the correlation length and vortex density in addition to the SF order parameter for a Kibble-Zurek protocol. We show that beside the ``freeze'' time $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{t}$, there exists another important time, ${t}_{\mathrm{eq}}$, at which an oscillating behavior of the SF amplitude starts. From calculations of the exponents of the correlation length and vortex density with respect to a quench time ${\ensuremath{\tau}}_{\mathrm{Q}}$, we obtain a physical picture of a coarsening process. Finally, we study how the system evolves after the quench. We give a global picture of dynamics of the Bose-Hubbard model.

arXiv:1711.08882•Physical Review A•Mar 28, 2018•29 citations

Route to supersolidity for the extended Bose-Hubbard model

M. Iskin

We use the Gutzwiller ansatz and analyze the phase diagram of the extended Bose-Hubbard Hamiltonian with on-site ($U$) and nearest-neighbor ($V$) repulsions. For $d$-dimensional hypercubic lattices, when $2\mathit{dV}<U$, it is well known that the ground state alternates between the charge-density-wave (CDW) and Mott insulators, and the supersolid (SS) phase occupies small regions around the CDW insulators. However, when $2\mathit{dV}>U$, in this Rapid Communication, we show that the ground state has only CDW insulators, and more importantly, the SS phase occupies a much larger region in the phase diagram, existing up to very large hopping values which could be orders of magnitude higher than that of the well-known case. In particular, the SS-superfluid phase boundary increases linearly as a function of hopping when $2\mathit{dV}\ensuremath{\gtrsim}1.5U$, for which the prospects of observing the SS phase with dipolar Bose gases loaded into optical lattices is much higher.

arXiv:1101.6067•Physical Review A•May 31, 2011•72 citations

Superfluid, Mott-insulator, and mass-density-wave phases in the one-dimensional extended Bose-Hubbard model

Ramesh V. Pai, Rahul Anil Pandit

We use the finite-size, density-matrix-renormalization-group (FSDMRG) method to obtain the phase diagram of the one-dimensional $(d=1)$ extended Bose-Hubbard model for density $\ensuremath{\rho}=1$ in the $U\ensuremath{-}V$ plane, where $U$ and $V$ are, respectively, onsite and nearest-neighbor interactions. The phase diagram comprises three phases: superfluid (SF), Mott insulator (MI), and mass-density-wave (MDW). For small values of $U$ and $V$, we get a reentrant SF-MI-SF phase transition. For intermediate values of interactions the SF phase is sandwiched between MI and MDW phases with continuous SF-MI and SF-MDW transitions. We show, by a detailed, finite-size scaling analysis, that the MI-SF transition is of Kosterlitz-Thouless (KT) type whereas the MDW-SF transition has both KT and two-dimensional Ising characters. For large values of $U$ and $V$ we get a direct, first-order, MI-MDW transition. The MI-SF, MDW-SF, and MI-MDW phase boundaries join at a bicritical point at $(U,V)=(8.5\ifmmode\pm\else\textpm\fi{}0.05,4.75\ifmmode\pm\else\textpm\fi{}0.05)$.

arXiv:cond-mat/0407011•Physical Review B•Mar 15, 2005•40 citations

Superfluid, phase-separation, supersolid and density wave phases in extended Bose-Hubbard model

Pallavi P. Gaude, Ananya Das, Ramesh V. Pai

We study the one-dimensional extended Bose-Hubbard model in incommensurate densities using the density matrix renormalization group with the mean-filed (CMFT+DMRG)and density matrix renormalization group (DMRG) methods. We focus on the region of the phase diagram where onsite interaction (U) is smaller than the nearest-neighbor interaction (V). We determine the superfluid and the density wave order parameters to identify superfluid, phase-separated, supersolid, density wave phases, and obtain the phase diagram in the density-onsite interaction(ρ, U) plane for fixed nearest-neighbor interaction. The supersolid and phase-separation phases exist for a wide range of densities. The density ranges for which the phase-separation is seen shrinking with theincrease in the onsite interaction.

AIP conference proceedings•Jan 1, 2024•0 citations

Successive phase transitions at finite temperatures toward the supersolid state in a three-dimensional extended Bose-Hubbard model

Keisuke Yamamoto, Synge Todo, Seiji Miyashita

We study the finite temperature properties of the extended Bose-Hubbard model on a cubic lattice. This model exhibits the so-called supersolid state. To start with, we investigate the ground-state phase diagram and find the supersolid state in the region of underdoped (less than the half-filled) density in contrast to the cases of one- and two-dimensional systems where the superfluid state and the solid state show the phase separation. Next, we investigate ordering processes at finite temperatures by quantum Monte Carlo simulations and find successive superfluid and solid phase transitions. There, we find that the two order parameters compete with each other. We establish a finite temperature phase diagram, which contains the superfluid, the solid, the supersolid, and the disordered phases. We develop a mean-field theory to analyze the ordering processes and compare the result with that obtained by simulations and discuss the mechanism of the competition of these two orders. We also study how the supersolid region shrinks as the on-site repulsion becomes strong.

arXiv:0807.0948•Physical Review B•Mar 3, 2009•35 citations
In this paper, we study the dynamics of the Bose-Hubbard model with the nearest-neighbor repulsion by using time-dependent Gutzwiller methods. Near the unit filling, the phase diagram of the model contains density wave (DW), supersolid (SS), and superfluid (SF). The three phases are separated by two second-order phase transitions. We study ``slow-quench'' dynamics by varying the hopping parameter in the Hamiltonian as a function of time. In the phase transitions from the DW to SS and from the DW to SF, we focus on how the SF order forms and study scaling laws of the SF correlation length, vortex density, etc. The results are compared with the Kibble-Zurek scaling. On the other hand from the SF to DW, we study how the DW order evolves with generation of the domain walls and vortices. Measurement of first-order SF coherence reveals interesting behavior in the DW regime.
In this paper, we study the nonequilibrium dynamics of the Bose–Hubbard model with the nearest-neighbor repulsion by using time-dependent Gutzwiller (GW) methods. In particular, we vary the hopping parameters in the Hamiltonian as a function of time, and investigate the dynamics of the system from the density wave (DW) to the superfluid (SF) crossing a first-order phase transition and vice versa. From the DW to SF, we find scaling laws for the correlation length and vortex density with respect to the quench time. This is a reminiscence of the Kibble–Zurek scaling for continuous phase transitions and contradicts the common expectation. We give a possible explanation for this observation. On the other hand from SF to DW, the system evolution depends on the initial SF state. When the initial state is the ground state obtained by the static GW methods, a coexisting state of the SF and DW domains forms after passing through the critical point. Coherence of the SF order parameter is lost as the system evolves. This is a phenomenon similar to the glass transition in classical systems. When the state starts from the SF with small local phase fluctuations, the system obtains a large size DW domain structure with thin domain walls.
In this paper, we study the dynamics of the Bose-Hubbard model by using time-dependent Gutzwiller methods. In particular, we vary the parameters in the Hamiltonian as a function of time, and investigate the temporal behavior of the system from the Mott insulator to the superfluid (SF) crossing a second-order phase transition. We first solve a time-dependent Schr\"odinger equation for the experimental setup recently done by Braun et al. [Proc. Natl. Acad. Sci. USA 112, 3641 (2015)] and show that the numerical and experimental results are in fairly good agreement. However, these results disagree with the Kibble-Zurek scaling. From our numerical study, we reveal a possible source of the discrepancy. Next, we calculate the critical exponents of the correlation length and vortex density in addition to the SF order parameter for a Kibble-Zurek protocol. We show that beside the ``freeze'' time $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{t}$, there exists another important time, teq{t}_{\mathrm{eq}}teq​, at which an oscillating behavior of the SF amplitude starts. From calculations of the exponents of the correlation length and vortex density with respect to a quench time ${\ensuremath{\tau}}_{\mathrm{Q}}$, we obtain a physical picture of a coarsening process. Finally, we study how the system evolves after the quench. We give a global picture of dynamics of the Bose-Hubbard model.
We use the Gutzwiller ansatz and analyze the phase diagram of the extended Bose-Hubbard Hamiltonian with on-site (UUU) and nearest-neighbor (VVV) repulsions. For ddd-dimensional hypercubic lattices, when 2dV<U2\mathit{dV}<U2dV<U, it is well known that the ground state alternates between the charge-density-wave (CDW) and Mott insulators, and the supersolid (SS) phase occupies small regions around the CDW insulators. However, when 2dV>U2\mathit{dV}>U2dV>U, in this Rapid Communication, we show that the ground state has only CDW insulators, and more importantly, the SS phase occupies a much larger region in the phase diagram, existing up to very large hopping values which could be orders of magnitude higher than that of the well-known case. In particular, the SS-superfluid phase boundary increases linearly as a function of hopping when $2\mathit{dV}\ensuremath{\gtrsim}1.5U$, for which the prospects of observing the SS phase with dipolar Bose gases loaded into optical lattices is much higher.
We use the finite-size, density-matrix-renormalization-group (FSDMRG) method to obtain the phase diagram of the one-dimensional (d=1)(d=1)(d=1) extended Bose-Hubbard model for density $\ensuremath{\rho}=1$ in the $U\ensuremath{-}V$ plane, where UUU and VVV are, respectively, onsite and nearest-neighbor interactions. The phase diagram comprises three phases: superfluid (SF), Mott insulator (MI), and mass-density-wave (MDW). For small values of UUU and VVV, we get a reentrant SF-MI-SF phase transition. For intermediate values of interactions the SF phase is sandwiched between MI and MDW phases with continuous SF-MI and SF-MDW transitions. We show, by a detailed, finite-size scaling analysis, that the MI-SF transition is of Kosterlitz-Thouless (KT) type whereas the MDW-SF transition has both KT and two-dimensional Ising characters. For large values of UUU and VVV we get a direct, first-order, MI-MDW transition. The MI-SF, MDW-SF, and MI-MDW phase boundaries join at a bicritical point at $(U,V)=(8.5\ifmmode\pm\else\textpm\fi{}0.05,4.75\ifmmode\pm\else\textpm\fi{}0.05)$.
We study the one-dimensional extended Bose-Hubbard model in incommensurate densities using the density matrix renormalization group with the mean-filed (CMFT+DMRG)and density matrix renormalization group (DMRG) methods. We focus on the region of the phase diagram where onsite interaction (U) is smaller than the nearest-neighbor interaction (V). We determine the superfluid and the density wave order parameters to identify superfluid, phase-separated, supersolid, density wave phases, and obtain the phase diagram in the density-onsite interaction(ρ, U) plane for fixed nearest-neighbor interaction. The supersolid and phase-separation phases exist for a wide range of densities. The density ranges for which the phase-separation is seen shrinking with theincrease in the onsite interaction.
We study the finite temperature properties of the extended Bose-Hubbard model on a cubic lattice. This model exhibits the so-called supersolid state. To start with, we investigate the ground-state phase diagram and find the supersolid state in the region of underdoped (less than the half-filled) density in contrast to the cases of one- and two-dimensional systems where the superfluid state and the solid state show the phase separation. Next, we investigate ordering processes at finite temperatures by quantum Monte Carlo simulations and find successive superfluid and solid phase transitions. There, we find that the two order parameters compete with each other. We establish a finite temperature phase diagram, which contains the superfluid, the solid, the supersolid, and the disordered phases. We develop a mean-field theory to analyze the ordering processes and compare the result with that obtained by simulations and discuss the mechanism of the competition of these two orders. We also study how the supersolid region shrinks as the on-site repulsion becomes strong.