To understand the mechanism of Mott transitions in case of no magnetic influence, superfluid-insulator (Mott) transitions in the S=0 Bose Hubbard model at unit filling are studied on the square and triangular lattices, using a variational Monte Carlo method. In trial many-body wave functions, we introduce various types of attractive correlation factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), which play a central role for Mott transitions, in addition to the onsite repulsive (Gutzwiller) factor. By optimizing distance-dependent parameters, we study various properties of this type of wave functions. With a hint from the Mott transition arising in a completely D-H bound state, we propose an improved picture of Mott transitions, by introducing two characteristic length scales, the D-H binding length
ξdh and the minimum D-D exclusion length
ξdd. Generally, a Mott transition occurs when
ξdh becomes comparable to
ξdd. In the conductive (superfluid) state, domains of D-H pairs overlap with each other (
ξdh>ξdd); thereby D and H can propagate independently as density carriers by successively exchanging the partners. In contrast, intersite repulsive Jastrow (D-D and H-H) factors have little importance for the Mott transition.
We argue that the binding between doubly occupied (doublon) and empty (holon) sites governs incoherent excitations and plays a key role in the Mott transition in strongly correlated Mott-Hubbard systems. We construct a new saddle-point solution with doublon-holon binding in the Kotliar-Ruckenstein slave-boson functional integral formulation of the Hubbard model. On a half-filled honeycomb lattice and square lattice, the ground state is found to exhibit a continuous transition from a paramagnetic semimetal/metal to an antiferromagnetic-ordered Slater insulator with coherent quasiparticles at
Uc1, followed by a Mott transition into an electron-fractionalized AF
∗ phase without coherent excitations at
Uc2. Such a phase structure appears to be generic for bipartite lattices without frustration. We show that doublon-holon binding unites the three important ideas of strong correlation: coherent quasiparticles, incoherent Hubbard bands, and deconfined Mott insulator.
An analytical solution of the Mott transition is obtained for the Hubbard model on the Bethe lattice in the large coordination number (
z) limit. The excitonic binding of doublons (doubly occupied sites) and holons (empty sites) is shown to be the origin of a continuous Mott transition between a metal and an emergent quantum spin liquid insulator. The doublon-holon binding theory enables a different large-
z limit and a different phase structure than the dynamical meanfield theory by allowing intersite spinon correlations to lift the
2N-fold degeneracy of the local moments in the insulating phase. We show that the spinons are coupled to doublons/holons by a dissipative compact U(1) gauge field that is in the deconfined phase, stabilizing the spin-charge separated gapless spin liquid Mott insulator.