Doping induced metal-insulator transition in two-dimensional Hubbard, t−Ut-U, and extended Hubbard, t−U−Wt-U-W, models

We show numerically that the nature of the doping induced metal-insulator transition in the two-dimensional Hubbard model is radically altered by the inclusion of a term, $W$, which depends upon a square of a single-particle nearest-neighbor hopping. This result is reached by computing the localization length, $\xi_l$, in the insulating state. At finite values of $W$ we find results consistent with $\xi_l \sim | \mu - \mu_c|^{- 1/2}$ where $\mu_c$ is the critical chemical potential. In contrast, $\xi_l \sim | \mu - \mu_c|^{-1/4}$ for the Hubbard model. At finite values of $W$, the presented numerical results imply that doping the antiferromagnetic Mott insulator leads to a $d_{x^2 - y ^2}$ superconductor.Comment: 19 pages (latex) including 7 figures in encapsulated postscript format. Submitted for publication in Phys. Rev.

Charge and Spin Structures of a dx2−y2d_{x^2 - y^2} Superconductor in the Proximity of an Antiferromagnetic Mott Insulator

To the Hubbard model on a square lattice we add an interaction, $W$, which depends upon the square of a near-neighbor hopping. We use zero temperature quantum Monte Carlo simulations on lattice sizes up to $16 \times 16$, to show that at half-filling and constant value of the Hubbard repulsion, the interaction $W$ triggers a quantum transition between an antiferromagnetic Mott insulator and a $d_{x^2 -y^2}$ superconductor. With a combination of finite temperature quantum Monte Carlo simulations and the Maximum Entropy method, we study spin and charge degrees of freedom in the superconducting state. We give numerical evidence for the occurrence of a finite temperature Kosterlitz-Thouless transition to the $d_{x^2 -y^2}$ superconducting state. Above and below the Kosterlitz-Thouless transition temperature, $T_{KT}$, we compute the one-electron density of states, $N(\omega)$, the spin relaxation rate $1/T_1$, as well as the imaginary and real part of the spin susceptibility $\chi(\vec{q},\omega)$. The spin dynamics are characterized by the vanishing of $1/T_1$ and divergence of $Re \chi(\vec{q} = (\pi,\pi), \omega = 0)$ in the low temperature limit. As $T_{KT}$ is approached $N(\omega)$ develops a pseudo-gap feature and below $T_{KT}$ $Im \chi(\vec{q} = (\pi,\pi), \omega)$ shows a peak at finite frequency.Comment: 46 pages (latex) including 14 figures in encapsulated postscript format. Submitted for publication in Phys. Rev.

Critical Exponents of the Metal-Insulator Transition in the Two-Dimensional Hubbard Model

We study the filling-controlled metal-insulator transition in the two-dimensional Hubbard model near half-filling with the use of zero temperature quantum Monte Carlo methods. In the metallic phase, the compressibility behaves as $\kappa \propto |\mu - \mu_c|^{-0.58\pm0.08}$ where $\mu_c$ is the critical chemical potential. In the insulating phase, the localization length follows $\xi_l \propto |\mu - \mu_c|^{-\nu_l}$ with $\nu_l = 0.26 \pm 0.05$. Under the assumption of hyperscaling, the compressibility data leads to a correlation length exponent $\nu_\kappa = 0.21 \pm 0.04$. Our results show that the exponents $\nu_\kappa$ and $\nu_l$ agree within statistical uncertainty. This confirms the assumption of hyperscaling with correlation length exponent $\nu = 1/4$ and dynamical exponent $z = 4$. In contrast the metal-insulator transition in the generic band insulators in all dimensions as well as in the one-dimensional Hubbard model satisfy the hyperscaling assumption with exponents $\nu = 1/2$ and $z = 2$.Comment: Two references added. The DVI file and PS figure files are also available at http://www.issp.u-tokyo.ac.jp/labs/riron/imada/furukawa/; to appear in J. Phys. Soc. Jpn 65 (1996) No.

Dynamic Exponent of t-J and t-J-W Model

Drude weight of optical conductivity is calculated at zero temperature by exact diagonalization for the two-dimensional t-J model with the two-particle term, $W$. For the ordinary t-J model with $W$=0, the scaling of the Drude weight $D \propto \delta^2$ for small doping concentration $\delta$ is obtained, which indicates anomalous dynamic exponent $z$=4 of the Mott transition. When $W$ is switched on, the dynamic exponent recovers its conventional value $z$=2. This corresponds to an incoherent-to-coherent transition associated with the switching of the two-particle transfer.Comment: LaTeX, JPSJ-style, 4 pages, 5 eps files, to appear in J. Phys. Soc. Jpn. vol.67, No.6 (1998

Quantum Transition between an Antiferromagnetic Mott Insulator and dx2−y2d_{x^2 - y^2} Superconductor in Two Dimensions

We consider a Hubbard model on a square lattice with an additional interaction, $W$, which depends upon the square of a near-neighbor hopping. At half-filling and a constant value of the Hubbard repulsion, increasing the strength of the interaction $W$ drives the system from an antiferromagnetic Mott insulator to a $d_{x^2 -y^2}$ superconductor. This conclusion is reached on the basis of zero temperature quantum Monte Carlo simulations on lattice sizes up to $16 \times 16$.Comment: 4 pages (latex) and 4 postscript figure

Effect of the W-term for a t-U-W Hubbard ladder

Antiferromagnetic and d_{x2-y2}-pairing correlations appear delicately balanced in the 2D Hubbard model. Whether doping can tip the balance to pairing is unclear and models with additional interaction terms have been studied. In one of these, the square of a local hopping kinetic energy H_W was found to favor pairing. However, such a term can be separated into a number of simpler processes and one would like to know which of these terms are responsible for enhancing the pairing. Here we analyze these processes for a 2-leg Hubbard ladder

Insulator-Metal Transition in the One and Two-Dimensional Hubbard Models

We use Quantum Monte Carlo methods to determine $T=0$ Green functions, $G(\vec{r}, \omega)$, on lattices up to $16 \times 16$ for the 2D Hubbard model at $U/t =4$. For chemical potentials, $\mu$, within the Hubbard gap, $|\mu | < \mu_c$, and at {\it long} distances, $\vec{r}$, $G(\vec{r}, \omega = \mu) \sim e^{ -|\vec{r}|/\xi_l}$ with critical behavior: $\xi_l \sim | \mu - \mu_c |^{-\nu}$, $\nu = 0.26 \pm 0.05$. This result stands in agreement with the assumption of hyperscaling with correlation exponent $\nu = 1/4$ and dynamical exponent $z = 4$. In contrast, the generic band insulator as well as the metal-insulator transition in the 1D Hubbard model are characterized by $\nu = 1/2$ and $z = 2$.Comment: 9 pages (latex) and 5 postscript figures. Submitted for publication in Phys. Rev. Let

Spin and charge dynamics of the ferromagnetic and antiferromagnetic two-dimensional half-filled Kondo lattice model

We present a detailed numerical study of spin and charge dynamics of the two-dimensional Kondo lattice model with hopping t and exchange J. At T=0 and J > 0, the competition between the RKKY interaction and Kondo effect triggers a quantum phase transition between magnetically ordered and disordered insulators: J_c/t = 1.45(5). The quasiparticle gap scales as |J|. S(q,\omega), evolves smoothly from its strong coupling form with spin gap at q = (\pi,\pi) to a spin wave form. At J>0, A(\vec{k},\omega) shows a dispersion relation following that of hybridized bands. For J < J_c this feature is supplemented by shadows thus pointing to a coexistence of Kondo screening and magnetism. For J < 0 A(\vec{k},\omega) is similar to that of non-interacting electrons in a staggered magnetic field. Spin, T_S, and charge, T_C, scales are defined. For weak to intermediate couplings, T_S marks the onset of antiferromagnetic fluctuations and follows a J^2 law. At strong couplings T_S scales as J. T_C scales as J both at weak and strong couplings. At and slightly below T_C we observe i) a rise in the resistivity as a function of decreasing temperature, ii) a dip in the integrated density of states at the Fermi energy and iii) the occurrence of hybridized bands in A(k,\omega). It is shown that in the weak coupling limit, the charge gap of order J is of magnetic origin. The specific heat shows a two peak structure, the low temperature peak being of magnetic origin. Our results are compared to various mean-field theories.Comment: 30 pages, 24 figure

Enhancement of Pairing Correlation and Spin Gap through Suppression of Single-Particle Dispersion in One-Dimensional Models

We investigate the effects of suppression of single-particle dispersion near the Fermi level on the spin gap and the singlet-pairing correlation by using the exact diagonalization method for finite-size systems. We consider strongly correlated one-dimensional models, which have flat band dispersions near wave number k=\pi/2, if the interactions are switched off. Our results for strongly correlated models show that the spin gap region expands as the single-particle dispersion becomes flatter. The region where the singlet-pairing correlation is the most dominant also expands in models with flatter band dispersions. Based on our numerical results, we propose a pairing mechanism induced by the flat-band dispersion.Comment: 5 pages, including 5 eps figures, to appear in J.Phys.Soc.Jpn Vol.69 No.