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

Abstract

We use quantum Monte Carlo methods to determine $T\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}0$ Green functions, $G(\stackrel{\ensuremath{\rightarrow}}{r},\ensuremath{\omega})$, on lattices up to $16\ifmmode\times\else\texttimes\fi{}16$ for the 2D Hubbard model at $U/t\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}4$. For chemical potentials $\ensuremath{\mu}$ within the Hubbard gap $|\ensuremath{\mu}|<{\ensuremath{\mu}}_{c}$ and at long distances $\stackrel{\ensuremath{\rightarrow}}{r}$, $G(\stackrel{\ensuremath{\rightarrow}}{r},\ensuremath{\omega}\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}\ensuremath{\mu})\ensuremath{\sim}{e}^{\ensuremath{-}|\stackrel{\ensuremath{\rightarrow}}{r}|/{\ensuremath{\xi}}_{l}}$ with critical behavior ${\ensuremath{\xi}}_{l}\ensuremath{\sim}|\ensuremath{\mu}\ensuremath{-}{\ensuremath{\mu}}_{c}{|}^{\ensuremath{-}\ensuremath{\nu}}$, $\ensuremath{\nu}\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}0.26\ifmmode\pm\else\textpm\fi{}0.05$. This result stands in agreement with the assumption of hyperscaling with correlation exponent $\ensuremath{\nu}\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}1/4$ and dynamical exponent $z\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}4$. In contrast, the generic band insulator as well as the metal-insulator transition in the 1D Hubbard model are characterized by $\ensuremath{\nu}\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}1/2$ and $z\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}2$.