We show that backflow correlations in the variational wave function for the Hubbard model greatly improve the previous results given by the Slater-Jastrow state, usually considered in this context. We provide evidence that, within this approach, it is possible to have a satisfactory connection with the strong-coupling regime. Moreover, we show that, for the Hubbard model on the lattice, backflow correlations are essentially short range, inducing an effective attraction between empty (holons) and doubly occupied sites (doublons). In the presence of frustration, we report the evidence that the metal to Mott-insulator transition is marked by a discontinuity of the double occupancy, together with a similar discontinuity of the kinetic term that does not change the number of holons and doublons, while the other kinetic terms are continuous across the transition. Finally, we show the estimation of the charge gap, obtained by particle-hole excitations \`a la Feynman over the ground-state wave function.
We show that backflow correlations in the variational wave function for the Hubbard model greatly improve the previous results given by the Slater-Jastrow state, usually considered in this context. We provide evidence that, within this approach, it is possible to have a satisfactory connection with the strong-coupling regime. Moreover, we show that, for the Hubbard model on the lattice, backflow correlations are essentially short range, inducing an effective attraction between empty (holons) and doubly occupied sites (doublons). In the presence of frustration, we report the evidence that the metal to Mott-insulator transition is marked by a discontinuity of the double occupancy, together with a similar discontinuity of the kinetic term that does not change the number of holons and doublons, while the other kinetic terms are continuous across the transition. Finally, we show the estimation of the charge gap, obtained by particle-hole excitations `a la Feynman over the ground-state wave function.
We study the phase diagram of the frustrated t - t ' Hubbard model on the square lattice by using a novel variational wave function. Taking the clue from the backflow correlations that have been introduced long-time ago by Feynman and Cohen and have been used for describing various interacting systems on the continuum (like liquid 3 He, the electron jellium, and metallic Hydrogen), we consider many-body correlations to construct a suitable approximation for the ground state of this correlated model on the lattice. In this way, a very accurate ansatz can be achieved both at weak and strong coupling. We present the evidence that an insulating and non-magnetic phase can be stabilized at strong coupling and sufficiently large frustrating ratio t ' = t .
In this thesis we consider a variational wave function approach as a possible route to describe the competition between disorder and strong electron-electron interaction in two dimensions. In particular we aim to obtain a transparent and physically intuitive understanding of the competition between these two localizing forces within the simplest model where they both are active, namely the disordered Hubbard model at half filling and in a square lattice. Our approach is based on an approximate form of the ground-state wave function, which we believe contains the physically relevant ingredients for a correct description of both the Mott and the Anderson insulators, where electrons are localized by the Coulomb repulsion and by disorder, respectively. For strongly interacting fermionic systems, a standard variational wave function is constructed by a correlation term acting on a Slater determinant, the latter being an uncorrelated metallic state. Previous variational calculations showed that a long-range density-density correlation factor, so called Jastrow factor, is needed to correctly describe the Mott insulator [9]. This term, which is collective by definition, correlates spatially charge uctuations, thus preventing their free motion that would otherwise imply metallic conductance. For this reason, our variational wave function does include such a term. Anderson localization is instead mostly a matter of single-particle wave functions, hence it pertains to the uncorrelated Slater determinant which the Jastrow factor acts onto. We consider both the case in which we enforce paramagnetism in the wave function and the case in which we allow for magnetic ordering. Summarizing briefly our results, we find that, when the variational wave function is forced to be paramagnetic, the Anderson insulator to Mott insulator transition is continuous. This transition can be captured by studying several quantities. In particular, a novel one that we have identified and that is easily accessible variationally is the disconnected density-density fluctuation at long wavelength, defined by lim where ^nq is the Fourier transform of the charge density at momentum q, (...) denotes quantum average at fixed disorder and the overbar represents the average over disorder configurations. We find that Ndisc q!0 is everywhere finite in the Anderson insulator and vanishes critically at the Mott transition, staying zero in the Mott insulator. When magnetism is allowed and the hopping only connects nearest neighbor sites, upon increasing interaction the paramagnetic Anderson insulator first turns antiferromagnetic and finally the magnetic and compressible Anderson insulator gives way to an incompressible antiferromagnetic Mott insulator. The optimized uncorrelated Slater determinant is always found to be the eigenstate of a disordered non-interacting effective Hamiltonian, which suggests that the model is never metallic. Finally, when magnetism is frustrated by a next to nearest neighbor hopping, the overall sequence of phases does not change. However, the paramagnetic to magnetic transition within the Anderson insulator basin of stability turns first order. Indeed, within the magnetically ordered phase, we find many almost degenerate paramagnetic states with well defined local moments. This is suggestive of an emerging glassy behavior when the competition between disorder and strong correlation is maximum.
The Mott metal-insulator transition in the Hubbard model is studied by constructing a dynamical slave-boson mean-field theory in the limit of large lattice coordination number z that incorporates the binding between doubly occupied (doublon) and empty (holon) sites. On the Mott insulating side where all doublons and holons bond in real space into excitonic pairs leading to the charge gap, the theory simplifies considerably to leading order in 1/z and becomes exact on the infinite-z Bethe lattice. An asymptotic solution is obtained for a continuous Mott transition associated with the closing of the charge gap at a critical value of the Hubbard Uc and the corresponding doublon density ndc, hopping ${\ensuremath{\chi}}_{d}^{c}$ and doublon-holon pairing ${\mathrm{\ensuremath{\Delta}}}_{d}^{c}$ amplitudes. We find ${U}_{c}={U}_{\mathrm{BR}}[1\ensuremath{-}2{n}_{d}^{c}\ensuremath{-}\sqrt{z}({\ensuremath{\chi}}_{d}^{c}+{\mathrm{\ensuremath{\Delta}}}_{d}^{c})]\ensuremath{\simeq}0.8{U}_{\mathrm{BR}}$, where UBR is the critical value for the Brinkman-Rice transition in the Gutzwiller approximation captured in the static mean-field solution of the slave-boson formulation of Kotliar and Ruckenstein. Thus, the Mott transition can be viewed as the quantum correction to the Brinkman-Rice transition due to doublon-holon binding. Quantitative comparisons are made to the results of the dynamical mean-field theory, showing good agreement. In the absence of magnetic order, the Mott insulator is a U(1) quantum spin liquid with nonzero intersite spinon hopping that survives the large-z limit and lifts the 2N-fold degeneracy of the local moments. We show that the spinons are coupled to the doublons/holons by a dissipative compact U(1) gauge field in the deconfined phase, realizing the spin-charge separated gapless spin liquid Mott insulator.
We consider the one-band Hubbard model on the square lattice by using variational and Green's function Monte Carlo methods, where the variational states contain Jastrow and backflow correlations on top of an uncorrelated wave function that includes BCS pairing and magnetic order. At half-filling, where the ground state is antiferromagnetically ordered for any value of the on-site interaction U, we can identify a hidden critical point UMott, above which a finite BCS pairing is stabilized in the wave function. The existence of this point is reminiscent of the Mott transition in the paramagnetic sector and determines a separation between a Slater insulator (at small values of U), where magnetism induces a potential energy gain, and a Mott insulator (at large values of U), where magnetic correlations drive a kinetic energy gain. Most importantly, the existence of UMott has crucial consequences when doping the system: We observe a tendency for phase separation into hole-rich and hole-poor regions only when doping the Slater insulator, while the system is uniform by doping the Mott insulator. Superconducting correlations are clearly observed above UMott, leading to the characteristic dome structure in doping. Furthermore, we show that the energy gain due to the presence of a finite BCS pairing above UMott shifts from the potential to the kinetic sector by increasing the value of the Coulomb repulsion.
The mechanism of Mott transition in the Hubbard model on the square lattice is studied without explicit introduction of magnetic and superconducting correlations, using a variational Monte Carlo method. In the trial wave functions, we consider various types of binding factors between a doubly-occupied site (doublon, D) and an empty site (holon, H), like a long-range type as well as a conventional nearest-neighbor type, and add independent long-range D-D (H-H) factors. It is found that a wide choice of D-H binding factor leads to Mott transitions at critical values near the band width. We renew the D-H binding picture of Mott transitions by introducing two characteristic length scales, the D-H binding length l_{DH} and the minimum D-D distance l_{DD}, which we appropriately estimate. A Mott transition takes place at l_{DH}=l_{DD}. In the metallic regime (l_{DH}>l_{DD}), the domains of D-H pairs overlap with one another, thereby doublons and holons can move independently by exchanging the partners one after another. In contrast, the D-D factors give only a minor contribution to 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.