We analyse the one-dimensional (1D) and the two-dimensional (2D) repulsive Hubbard models (HM) for densities slightly away from half-filling through the behaviour of two central quantities of a system: the uniform charge and spin susceptibilities. We point out that a consistent renormalization group treatment of them can only be achieved within a two-loop approach or beyond. In the 1D HM, we show that this scheme reproduces correctly the metallic behaviour given by the well-known Luttinger liquid fixed-point result. Then, we use the same approach to deal with the more complicated 2D HM. In this case, we are able to show that both uniform susceptibilities become suppressed for moderate interaction parameters as one takes the system towards the Fermi surface. Therefore, this result adds further support to the interpretation that those systems are in fact insulating spin liquids. Later, we perform the same calculations in 2D using the conventional random phase approximation, and establish clearly a comparison between the two schemes.
We study the phase diagram of a one-dimensional extended Hubbard model with antiferromagnetic exchange interaction analytically and numerically. The bosonization and transfer-matrix renormalization group methods are used in the corresponding coupling regimes. At half-filling, the system is a Mott insulator with a finite spin excitation gap if the on-site Coulomb repulsion is fairly smaller than the antiferromagnetic exchange
J. This Mott insulator is characterized by the bond-charge-density-wave order or spontaneously dimerization. In the weak-coupling regime where the spin-charge separation holds approximately, the critical point separating the gapless and gapped spin liquid phases is ${U}_{c}\ensuremath{\sim}J∕2$. However, as
J increases, the spin-charge couplings become important and the critical point
Uc is significantly suppressed and eventually tends to zero as $J\ensuremath{\rightarrow}\ensuremath{\infty}$. Away from half-filling, the charge gap completely collapses but the spin gap persists.
We study the dimensional crossover from one to two dimensions in half-filled and lightly doped, weakly interacting N-leg Hubbard ladders. In this case, the Hubbard ladders are equivalent to an N-band model. Using renormalization-group techniques, we find that the half-filled ladders exhibit (in the spin sector) an odd-even effect only below a crossover energy ${E}_{c}\ensuremath{\propto}\mathrm{exp}[\ensuremath{-}\ensuremath{\alpha}\mathrm{exp}(\ensuremath{\gamma}N)]$ $(\ensuremath{\alpha}\ensuremath{\ll}1$ and $\ensuremath{\gamma}\ensuremath{\sim}1$ depend on the interaction strength and on the hopping matrix elements): Below
Ec, the dominant interactions take place within band pairs $(j,N+1\ensuremath{-}j)$ [and within the band
(N+1)/2 for N odd], such that even-leg ladders are an insulating spin liquid, while odd-leg ladders have one gapless spinon mode. In contrast, above the energy scale
Ec, all bands are interacting with each other and the system is a two-dimensional- (2D)-like (insulating) antiferromagnet; we obtain an analytical expression for the Hamiltonian which is similar to the 2D Heisenberg antiferromagnet. Bosonization techniques show that in the charge sector the Mott insulator is as well below and above
Ec of the same type as in
N/2 half-filled two-leg ladders. Doping away from half-filling, we find that the effect of an increasing doping $\ensuremath{\delta}$ is very similar to decreasing the number of legs N: In both cases interactions between unpaired bands are suppressed and thereby the antiferromagnetic correlations reduced. The resulting band pairs form then insulating spin liquids and when doped, there is a spin gap, but phase coherence exists only within the band pairs. At higher doping levels ${\ensuremath{\delta}}_{c}={\ensuremath{\delta}}_{c}(N),$ phase coherence between all band pairs sets in and the system becomes a 2D-like d-wave superconductor $[{\ensuremath{\delta}}_{c}(N)\ensuremath{\rightarrow}0$ for $\stackrel{\ensuremath{\rightarrow}}{N}\ensuremath{\infty}].$
We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near
1/8 doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about
1/5 doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.
The Hubbard model on a square lattice in the ground state is investigated. Various aspects of Mott transition at half-filling are clarified from the quantum Monte Carlo study. Critical exponents of the transition for charge and spin correlations are estimated. For the doping concentration δ, the charge susceptibility is proportional to δ -1 indicating the divergence of charge mass for δ→0, while the spin susceptibility is always finite and not singular near δ=0. Several theoretical consequences of the above results are discussed. Incommensurate peak amplitude of spin correlation scales as (δ-δ c ) -1 with δ c <0.01 in disagreement with RPA results ln (δ-δ c ). The antiferromagnetic order at the half-filling also shows strong correlation character and is not sensitive to the shape of the fermi surface, while the incommensurate peak observed away from the half-filling is sensitively suppressed by the loss of partial nesting consistently with the weak coupling picture. Comparisons with experimental indications of high- T c oxides are made.
Abstract We present a T = 0 K renormalization group (RG) phase diagram for the electronic Hubbard model in two dimensions on the square lattice at half filling. The RG procedure treats quantum fluctuations in the single particle occupation number nonperturbatively via the unitarily decoupling of one electronic state at every RG step. The resulting phase diagram thus possesses the quantum fluctuation energy scale (ω) as one of its axes. A relation is derived between ω and the effective temperature scale upto which gapless, as well as emergent gapped phases can be obtained. We find that the normal and insulating phases of the half-filled Hubbard model correspond, for any on-site repulsion, to a marginal Fermi liquid normal phase and a topologically-ordered gapped Mott insulating liquid respectively. The marginal Fermi liquid is found to arise from singular forward scattering in directions normal to the nested Fermi surface, while singular backscattering events lead to Mott liquid state. The transition between these two phases involves passage through a pseudogapped phase bookended by Fermi surface topology-changing Lifshitz transitions. The pseudogap phase is observed to arise from the electronic differentiation encoded within the nested Fermi surface, and involves the gradual gapping of the Fermi surface from antinodes to nodes via charge and spin excitations that are mutually entangled. We obtain effective Hamiltonians for various phases, as well as wavefunctions for the low-energy many-body eigenstates of the Mott liquid. Benchmarking of the ground-state energy per particle and the double-occupancy fraction for the Mott liquid against existing numerical results yields excellent agreement. Presence of a Néel ordering symmetry-breaking perturbation in the RG leads to an antiferromagnetic spin-ordered charge insulating Mott state. Our results thus offer novel insights on a variety of aspects of the Mott–Hubbard problem, and can be extended to the doped system.