We study spin and charge dynamics of stripes in doped Mott insulators by considering a two-dimensional Hubbard model with N fermion flavors. For N = 2 we recover the normal one-band model while for N → ∞ a spin density wave mean-field solution. For all band fillings, lattice topologies and N = 4 n , the model may be solved by means of quantum Monte Carlo methods without encountering the sign problem. At N = 4 and in the vicinity of the Mott insulator we find a stripe phase. This phase has a quasiparticle gap but conducts due to long-wavelength low-lying collective charge modes.
We study spin and charge dynamics of stripes in doped Mott insulators by considering a two-dimensional Hubbard model with N fermion flavors. For N = 2 we recover the normal one-band model while for N → ∞ a spin density wave mean-field solution. For all band fillings, lattice topologies and N = 4 n , the model may be solved by means of quantum Monte Carlo methods without encountering the sign problem. At N = 4 and in the vicinity of the Mott insulator we find a stripe phase. This phase has a quasiparticle gap but conducts due to long-wavelength low-lying collective charge modes.
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.
We introduce a slave-fermion formulation in which to study the charge dynamics of the half-filled Hubbard model on the square lattice. In this description, the charge degrees of freedom are represented by fermionic holons and doublons and the Mott-insulating characteristics of the ground state are the consequence of holon–doublon bound-state formation. The bosonic spin degrees of freedom are described by the antiferromagnetic Heisenberg model, yielding long-ranged (Néel) magnetic order at zero temperature. Within this framework and in the self-consistent Born approximation, we perform systematic calculations of the average double occupancy, the electronic density of states, the spectral function and the optical conductivity. Qualitatively, our method reproduces the lower and upper Hubbard bands, the spectral-weight transfer into a coherent quasiparticle band at their lower edges and the renormalisation of the Mott gap, which is associated with holon–doublon binding, due to the interactions of both quasiparticle species with the magnons. The zeros of the Green function at the chemical potential give the Luttinger volume, the poles of the self-energy reflect the underlying quasiparticle dispersion with a spin-renormalised hopping parameter and the optical gap is directly related to the Mott gap. Quantitatively, the square-lattice Hubbard model is one of the best-characterised problems in correlated condensed matter and many numerical calculations, all with different strengths and weaknesses, exist with which to benchmark our approach. From the semi-quantitative accuracy of our results for all but the weakest interaction strengths, we conclude that a self-consistent treatment of the spin-fluctuation effects on the charge degrees of freedom captures all the essential physics of the antiferromagnetic Mott–Hubbard insulator. We remark in addition that an analytical approximation with these properties serves a vital function in developing a full understanding of the fundamental physics of the Mott state, both in the antiferromagnetic insulator and at finite temperatures and dopings.
The high-temperature superconducting cuprates are governed by intertwined spin, charge, and superconducting orders. While various state-of-the-art numerical methods have demonstrated that these phases also manifest themselves in doped Hubbard models, they differ on which is the actual ground state. Finite cluster methods typically indicate that stripe order dominates while embedded quantum cluster methods, which access the thermodynamic limit by treating long-range correlations with a dynamical mean field, conclude that superconductivity does. Here, we report the observation of fluctuating spin and charge stripes in the doped single-band Hubbard model using a quantum Monte Carlo dynamical cluster approximation (DCA) method. By resolving both the fluctuating spin and charge orders using DCA, we demonstrate for the first time that they survive in the doped Hubbard model in the thermodynamic limit. This discovery also provides a new opportunity to study the influence of fluctuating stripe correlations on the model's pairing correlations within a unified numerical framework.
Though most fermionic Mott insulators order at low temperatures, ordering is ancillary to their insulating behavior. Our emphasis here is on disentangling ordering from the intrinsic strongly correlated physics of a doped half-filled band. To this end, we focus on the two-dimensional Hubbard model. Because the charge gap arises from on-site correlations, we have been refining the nonperturbative approach of Matsumoto and Mancini [Phys. Rev. B 55, 2095 (1997)] which incorporates local physics. Crucial to this method is a self-consistent two-site dynamical cluster expansion which builds in the nearest-neighbor energy scale J. At half-filling, we find that the spectral function possesses a gap of order U and is devoid of any coherent quasiparticle peaks although ordering or charge fractionalization are absent. At low temperatures, local antiferromagnetic correlations emerge. In the doped case, we find that the Fermi surface exceeds the Luttinger volume. The breakdown of Luttinger's theorem in the underdoped regime is traced both to the dynamically generated Mott gap as well to a nonvanishing of the imaginary part of the self-energy at the Fermi level. Spectral weight transfer across the Mott gap also emerges as a ubiquitous feature of a doped Mott insulator and suggests that high- and low-energy scales are inseparable. Additionally in the underdoped regime, we find that a pseudogap exists in the single-particle density of states as well as in the heat capacity. The pseudogap (which is set by the energy scale t2/U) is argued to be a ubiquitous feature of a lightly doped Mott state and simply represents the fact that hole transport involves double occupancy. In analogy with the Mott gap and antiferromagnetism, we propose that ordering may also accompany the formation of a pseudogap. We suggest a current pattern within a one-band model that preserves translational but breaks time-reversal symmetry along the canonical x and y axes but not along $x=\ifmmode\pm\else\textpm\fi{}y$ that is consistent with the experimental observations. Finally, we show that the Hall coefficient in a doped Mott insulator must change sign at a doping level x<1/3. The sign change is tied to a termination of strong correlation physics in the doped Mott state.
With the aid of both a semianalytical and a numerically exact method, we investigate the charge dynamics in the vicinity of half-filling in the one- and two-dimensional $t\text{\ensuremath{-}}J$ model derived from a Fermi-Hubbard model in the limit of large interaction U and hence small exchange coupling J. The spin degrees of freedom are taken to be disordered. So we consider the limit $0<J\ensuremath{\ll}T\ensuremath{\ll}W$, where W is the bandwidth. We focus on evaluating the local spectral density of a single hole excitation and the charge gap that separates the upper and the lower Hubbard band. We find indications that no band edges exist if the magnetic exchange is taken into account; instead of band edges, Gaussian tails seem to appear. A discussion of the underlying physics is provided.
Experiments on layered materials call for a study of the influence of short-range spin correlations on the Mott transition. To this end, we solve the cellular dynamical mean-field equations for the Hubbard model on a plaquette with continuous-time quantum Monte Carlo calculations. The normal-state phase diagram as a function of temperature T, interaction strength U, and filling n reveals that upon increasing n towards the insulator, there is a surface of first-order transition between two metals at nonzero doping. For T above the critical end line there is a maximum in scattering rate.