No Arabic abstract
We present the global phase diagram of the extended boson Hubbard model on a simple cubic lattice by quantum Monte Carlo simulation with worm update algorithm. Four kinds of phases are supported by this model, including superfluid, supersolid, Mott, and charge density wave (CDW) states, which are identified in the phase diagram of chemical potential $mu$ versus nearest neighbor interaction V . By changing the chemical potential, a continuous transition is found from the Mott phase to a superfluid phase without breaking the translational symmetry. For an insulating CDW state, adding particles to it gives rise to a continuous transition to a supersolid phase, while removing particles usually leads to a first-order one to either supersolid or superfluid phase. By tuning the nearest neighbor interaction, one can realize the transition between two insulating phases, Mott and CDW with the same particle density, which turns out to be of the first-order. We also demonstrate that a supersolid phase with average particle density less than 1/2 can exist in a small region of $mu$ - V phase diagram.
The `dynamic Hubbard Hamiltonian describes interacting fermions on a lattice whose on-site repulsion is modulated by a coupling to a fluctuating bosonic field. We investigate one such model, introduced by Hirsch, using the determinant Quantum Monte Carlo method. Our key result is that the extended s-wave pairing vertex, repulsive in the usual static Hubbard model, becomes attractive as the coupling to the fluctuating Bose field increases. The sign problem prevents us from exploring a low enough temperature to see if a superconducting transition occurs. We also observe a stabilization of antiferromagnetic correlations and the Mott gap near half-filling, and a near linear behavior of the energy as a function of particle density which indicates a tendency toward phase separation.
We present an unbiased numerical density-matrix renormalization group study of the one-dimensional Bose-Hubbard model supplemented by nearest-neighbor Coulomb interaction and bond dimerization. It places the emphasis on the determination of the ground-state phase diagram and shows that, besides dimerized Mott and density-wave insulating phases, an intermediate symmetry-protected topological Haldane insulator emerges at weak Coulomb interactions for filling factor one, which disappears, however, when the dimerization becomes too large. Analyzing the critical behavior of the model, we prove that the phase boundaries of the Haldane phase to Mott insulator and density-wave states belong to the Gaussian and Ising universality classes with central charges $c=1$ and $c=1/2$, respectively, and merge in a tricritical point. Interestingly we can demonstrate a direct Ising quantum phase transition between the dimerized Mott and density-wave phases above the tricritical point. The corresponding transition line terminates at a critical end point that belongs to the universality class of the dilute Ising model with $c=7/10$. At even stronger Coulomb interactions the transition becomes first order.
We investigate the phase diagram of the half-filled SU(N) Hubbard-Heisenberg model with hopping t, exchange J and Hubbard U, on a square lattice. In the large-N limit, and as a function of decreasing values of t/J, the model shows a transition from a d-density wave state to a spin dimerized insulator. A similar behavior is observed at N=6 whereas at N=2 a spin density wave insulating ground state is stabilized. The N=4 model, has a d-density wave ground state at large values of t/J which as a function of decreasing values of t/J becomes unstable to an insulating state with no apparent lattice and spin broken symmetries. In this state, the staggered spin-spin correlations decay as a power-law,resulting in gapless spin excitations at q = (pi,pi). Furthermore, low lying spin modes with small spectral weight are apparent around the wave vectors q = (0,pi) and q = (pi,0). This gapless spin liquid state is equally found in the SU(4) Heisenberg model in the self-adjoint antisymmetric representation. An interpretation of this state in terms of a pi-flux phase is offered. Our results stem from projective (T=0) quantum Monte-Carlo simulations on lattice sizes ranging up to 24 X 24.
We report large scale determinant Quantum Monte Carlo calculations of the effective bandwidth, momentum distribution, and magnetic correlations of the square lattice fermion Hubbard Hamiltonian at half-filling. The sharp Fermi surface of the non-interacting limit is significantly broadened by the electronic correlations, but retains signatures of the approach to the edges of the first Brillouin zone as the density increases. Finite size scaling of simulations on large lattices allows us to extract the interaction dependence of the antiferromagnetic order parameter, exhibiting its evolution from weak coupling to the strong coupling Heisenberg limit. Our lattices provide improved resolution of the Greens function in momentum space, allowing a more quantitative comparison with time-of-flight optical lattice experiments.
We study thermodynamics of the 3D Hubbard model at half filling on approach to the Neel transition by means of large-scale unbiased Diagrammatic Determinant Monte Carlo simulations. We obtain the transition temperature in the strongly correlated regime, as well as temperature dependence of energy, entropy, double occupancy, and the nearest-neighbor spin correlation function. Our results improve the accuracy of previous unbiased studies and present accurate benchmarks in the ongoing effort to realize the antiferromagnetic state of matter with ultracold atoms in optical lattices.