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.

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

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.

Asymmetric spin-1/2 two-leg ladders

We consider asymmetric spin-1/2 two-leg ladders with non-equal antiferromagnetic (AF) couplings J_|| and \kappa J_|| along legs (\kappa <= 1) and ferromagnetic rung coupling, J_\perp. This model is characterized by a gap \Delta in the spectrum of spin excitations. We show that in the large J_\perp limit this gap is equivalent to the Haldane gap for the AF spin-1 chain, irrespective of the asymmetry of the ladder. The behavior of the gap at small rung coupling falls in two different universality classes. The first class, which is best understood from the case of the conventional symmetric ladder at \kappa=1, admits a linear scaling for the spin gap \Delta ~ J_\perp. The second class appears for a strong asymmetry of the coupling along legs, \kappa J_|| << J_\perp << J_|| and is characterized by two energy scales: the exponentially small spin gap \Delta ~ J_\perp \exp(-J_|| / J_\perp), and the bandwidth of the low-lying excitations induced by a Suhl-Nakamura indirect exchange ~ J_\perp^2 /J_|| . We report numerical results obtained by exact diagonalization, density matrix renormalization group and quantum Monte Carlo simulations for the spin gap and various spin correlation functions. Our data indicate that the behavior of the string order parameter, characterizing the hidden AF order in Haldane phase, is different in the limiting cases of weak and strong asymmetry. On the basis of the numerical data, we propose a low-energy theory of effective spin-1 variables, pertaining to large blocks on a decimated lattice.Comment: 18 pages, 11 figure

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

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.

Interaction driven metal-insulator transition in strained graphene

The question of whether electron-electron interactions can drive a metal to insulator transition in graphene under realistic experimental conditions is addressed. Using three representative methods to calculate the effective long-range Coulomb interaction between $\pi$-electrons in graphene and solving for the ground state using quantum Monte Carlo methods, we argue that without strain, graphene remains metallic and changing the substrate from SiO$_2$ to suspended samples hardly makes any difference. In contrast, applying a rather large -- but experimentally realistic -- uniform and isotropic strain of about $15\%$ seems to be a promising route to making graphene an antiferromagnetic Mott insulator.Comment: Updated version: 6 pages, 3 figure

The role of electron-electron interactions in two-dimensional Dirac fermions

The role of electron-electron interactions on two-dimensional Dirac fermions remains enigmatic. Using a combination of nonperturbative numerical and analytical techniques that incorporate both the contact and long-range parts of the Coulomb interaction, we identify the two previously discussed regimes: a Gross-Neveu transition to a strongly correlated Mott insulator, and a semi-metallic state with a logarithmically diverging Fermi velocity accurately described by the random phase approximation. Most interestingly, experimental realizations of Dirac fermions span the crossover between these two regimes providing the physical mechanism that masks this velocity divergence. We explain several long-standing mysteries including why the observed Fermi velocity in graphene is consistently about 20 percent larger than the best values calculated using ab initio and why graphene on different substrates show different behavior.Comment: 11 pages, 4 figure

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 gap and string order parameter in the ferromagnetic Spiral Staircase Heisenberg Ladder: a quantum Monte Carlo study

We consider a spin-1/2 ladder with a ferromagnetic rung coupling J_\perp and inequivalent chains. This model is obtained by a twist (\theta) deformation of the ladder and interpolates between the isotropic ladder (\theta=0) and the SU(2) ferromagnetic Kondo necklace model (\theta=\pi). We show that the ground state in the (\theta,J_\perp) plane has a finite string order parameter characterising the Haldane phase. Twisting the chain introduces a new energy scale, which we interpret in terms of a Suhl-Nakamura interaction. As a consequence we observe a crossover in the scaling of the spin gap at weak coupling from \Delta/J_\| \propto J_\perp/J_\| for \theta < \theta_c \simeq 8\pi/9 to \Delta/J_\| \propto (J_\perp/J_\|)^2 for \theta > \theta_c. Those results are obtained on the basis of large scale Quantum Monte Carlo calculations.Comment: 4 page