## Mott transition in cuprate high-temperature superconductors

Takashi Yanagisawa a and Mitake Miyazaki b a Electronics and Photonics Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba Central 2, 1-1-1 Umezono, Tsukuba 305-8568, Japan b

Hakodate National College of Technology, 14-1 Tokura, Hakodate, Hokkaido 042-8501, Japan

In this study, we investigate the metal-insulator transition of charge transfer type in hightemperature cuprates. We first show that we must introduce a new band parameter in the three-band d-p model to reproduce the Fermi surface of high temperature cuprates such as BSCCO, YBCO and Hg1201. We present a new wave function of a Mott insulator based on the improved Gutzwiller function, and show that there is a transition from a metal to a charge-transfer insulator for such parameters by using the variational Monte Carlo method. This transition occurs when the level difference ∆ dp ≡ /epsilon1 p -/epsilon1 d between d and p orbitals reaches a critical value (∆ dp ) c . The energy gain ∆ E , measured from the limit of large ∆ dp , is proportional to 1 / ∆ dp for ∆ dp &gt; (∆ dp ) c : ∆ E ∝ -t 2 dp / ∆ dp . We obtain (∆ dp ) c /similarequal 2 t dp using the realistic band parameters.

PACS numbers: 71.10.-w, 71.27.+a

Introduction The study of high-temperature superconductors has been intensively addressed since the discovery of cuprate high-temperature superconductors. The research of mechanism of superconductivity (SC) in hightemperature superconductors has attracted much attention and has been extensively studied using various models. It has been established that the Cooper pairs of high-temperature cuprates have the d -wave symmetry in the hole-doped materials[1]. Therefore the electron correlation plays an important and the CuO 2 plane in cuprates plays a key role for the appearance of superconductivity[2-5]. The three-band d-p model is the most fundamental model for high-temperature cuprates[3-12].

The purpose of this paper is to investigate the effect of electron correlation in the half-filled case, that is, to discuss the metal-insulator transition due to the on-site Coulomb repulsion. As was discussed in Ref.[2], insulators are classified in terms of charge-transfer insulator or Mott insulator. The cuprates belong to the class of charge-transfer insulators. When the level difference ∆ dp ≡ /epsilon1 p -/epsilon1 d between d and p orbitals is large, the ground state will be insulating when the Coulomb repulsion U d on copper sites is large. (In this paper, we use the hole picture.) When U d is large, there will be a transition form a metal to an insulator as ∆ dp is increased. This is the Mott transition of charge-transfer type.

We will investigate this transition by using a variational Monte Carlo (VMC) method. In correlated electron systems, we must take into account the electron correlation correctly. Using the VMC method we can treat the electron systems properly from weakly to strongly correlated regions. We propose a wave function for an insulator on the basis of the Gutzwiller ansatz and examine the ground state within the space of variational functions. The expectation values are evaluated by using the variational Monte Carlo algorithm[13-18].

We first discuss the Mott state of the single-band Hubbard model by proposing a Mott-state wave function. We show that there is a metal-insulator transition as the on- site Coulomb repulsion U is increased. The energy gain, compared to the limit of large U , is proportional to the exchange interaction J ≡ 4 t 2 /U in the insulating state. The wave function is generalized straightforwardly to the d-p model. In the localized region, where ∆ dp is greater than the critical value (∆ dp ) c , the energy gain ∆ E is proportional to -1 / ∆ dp , that is, ∆ E ∝ -t 2 dp / ∆ dp . In this region we have the insulating ground state. (∆ dp ) c is of the order of the transfer integral t dp between holes in adjacent copper and oxygen atoms. This value is consistent with the result obtained by the dynamical mean-field theory[19].

Hamiltonian The three-band model that explicitly includes oxygen p and copper d orbitals contains the parameters U d , U p , t dp , t pp , /epsilon1 d and /epsilon1 p . The Hamiltonian is written as

$$\text{ effect } \text{ parameters } U _ { d }, U _ { p }, t _ { d p }, t _ { p p }, \epsilon _ { d } \text{ and } \epsilon _ { p }. \text{ The Hamiltonian is} \\ \text{ is, to } \text{ written as} \\ \text{n-site} \\ \text{nsu-} \\ \text{me class} \\ \text{t-road} \\ \text{t-ice} \\ \text{n-rot} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t- ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-Ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ICE} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-nce} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-rice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ie} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ies} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ise} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-i-} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-race} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-iere} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-iser} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-is-} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t+ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t+ise} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t+iser} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t+ie} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t+ies} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice} \\ \text{t-ice}$$

d iσ and d † iσ are the operators for the d holes. p i ± ˆ x/ σ 2 and p † i ± ˆ x/ σ 2 denote the operators for the p holes at the site R i ± ˆ x/ 2 , and in a similar way p i ± ˆ 2 y/ σ and p † i ± ˆ 2 y/ σ are defined. t dp is the transfer integral between adjacent Cu and O orbitals and t pp is that between nearest p orbitals. 〈〈 ij 〉〉 denotes a next nearest-neighbor pair of copper sites. U d is the strength of the on-site Coulomb energy between d holes. In this paper we neglect U p

among p holes because U p is small compared to U d [2022]. In the low-doping region, U p will be of minor importance because the p-hole concentration is small[23]. The parameter values were estimated as, for example, U d = 10 5, . U p = 4 0 and . U dp = 1 2 in eV[20] where . U dp is the nearest-neighbor Coulomb interaction between holes on adjacent Cu and O orbitals. In this paper we neglect U dp because U dp is small compared to U d . We use the notation ∆ dp = /epsilon1 p -/epsilon1 d . The number of sites is denoted as N , and the total number of atoms is N a = 3 N . Our study is done within the hole picture where the lowest band is occupied up to the Fermi energy µ .

The single-band Hubbard model is also important in the study of strongly correlated electron systems[24]. This model is regarded as an approximation to the threeband model. The Hamiltonian is given by

$$H \ & = \ - t \sum _ { \langle i j \rangle } ( c _ { i \sigma } ^ { \dagger } c _ { j \sigma } + \text{h.c.} ) - t ^ { \prime } \sum _ { \langle \langle j \ell \rangle \rangle } ( c _ { j \sigma } ^ { \dagger } c _ { \ell \sigma } + h. c. ) \\ & + \ U \sum _ { i } n _ { i \uparrow } n _ { i \downarrow },$$

where 〈 ij 〉 and 〈〈 j/lscript 〉〉 indicate the nearest neighbor and next-nearest neighbor pairs of sites, respectively. c iσ and c † iσ indicate the operators of d electrons. U is the on-site Coulomb repulsion.

Band parameters and the Fermi surface We need an additional band parameter because we cannot reproduce the deformed Fermi surface for cuprates by means of only t dp and t pp [25]. Thus we have introduced the parameter t ′ d in the Hamiltonian in the previous section. We show that the inclusion of t ′ d will enable us to reproduce the curvature of the Fermi surface. The non-zero t ′ d may be attributed to the integral between d x 2 -y 2 and d 3 z 2 -r 2 or d s orbitals. As we will show, the large t ′ d is not required to discribe the curvature of the Fermi surface. We must mention that there is another method to explain the curvature of the Fermi surface. For example, the inclusion of the O p x -O p x transfer integrals leads to the deformed Fermi surface[26].

Typical Fermi surfaces of cuprate superconductors have been reported for, for example, (La,Sr) CuO 2 4 (LSCO) and Bi Sr CaCu O 2 2 2 8+ δ . The Fermi surface for Bi Sr CaCu O 2 2 2 8+ δ (Bi2212)[27] and Tl Ba CuO 2 2 6+ δ [28] is deformed considerably, while that for LSCO is likely the straight line.

For LSCO, the band parameter is estimated as t ′ ∼ -0 12 when we fit by using the single-band model. . On the other hand, Tl2201 ( T c = 93K) and Hg1201 ( T c = 98K) band calculations by Singh and Pickett[29] give very much deformed Fermi surfaces that can be fitted by large | t ′ | such as t ′ ∼ -0 4. . For Tl2201, an Angular Magnetoresistance Oscillations (AMRO) work[28] gives information of the Fermi surface, which allows to get t ′ ∼ -0 2 . and t ′′ ∼ 0 165. There is also an Angle-Resolved Photoe-. mission Study (ARPES)[30], which provides similar values. In the case of Hg1201, there is an ARPES work[31], from which we obtain by fitting t ′ ∼ -0 2 and . t ′′ ∼ 0 175. .

We show the Fermi surface for the d-p model in Fig.1, where we set t pp = 0 and t ′ d = 0. The Fermi surface shown in Fig.1 is consistent with the Fermi surface for (La,Sr) CuO . 2 4 However, the deformed Fermi surfaces cannot be well fitted by using only t dp and t pp . We show the Fermi surface with t ′ d in Fig.2; the figure indicates that the inclusion of t ′ d gives a deformed Fermi surface. This Fermi surface agrees with that for Bi2212, Tl2201 and Hg1201.

FIG. 1: Fermi surface of the 2D d-p model for t pp = 0, t ′ d = 0 and ∆ dp = 2 0 in units of . t dp . The doped-hole density is n h ∼ 0 13. . The dashed line is for n h ∼ 0 0 (half-filled case). . This Fermi surface is for LSCO.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-04ee0f85-370c-4609-95cc-309c4ef03ef0.jpg)

FIG. 2: Fermi surface of the 2D d-p model for t pp = 0 4 and . ∆ dp = 2 0 where . t ′ d = -0 1, . t ′ d = -0 2 and . t ′ d = -0 3 in units . of t dp . The carrier density is n h ∼ 0 1. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-de3539a0-a9df-43df-85e2-93c5097b3c92.jpg)

Wave function of Mott state Here we propose a wave function to represent a Mott insulator. We first discuss it for the single-band Hubbard model[24]. The energy is measured in units of t in this section. The charge-transfer Mott state of the three-band model will be given by a generalization of the one-band Mott state. The Gutzwiller

function ψ G itself does not describe the insulating state because this function has no kinetic energy gain in the limit g → 0. Wave functions for the Mott transition have been proposed for the single-band Hubbard model by considering the doublon-holon correlation[32, 33] or backflow correlations[34]. In the latter, a variational wave function with a Jastrow factor is considered. It seems, however, not straightforward to generalize these wave functions to the three-band case.

The Gutzwiller wave function is given as

$$\psi _ { G } = P _ { G } \psi _ { 0 }, \text{ \quad \ \ } ( 3 ) \text{ \ \ } \text{$\quad$\nolimits$}$$

where P G is the Gutzwiller projection operator given by P G = ∏ i [1 -(1 -g n ) i ↑ n i ↓ ] with the variational parameter g in the range from 0 to unity. P G controls the on-site electron correlation and ψ 0 is the Fermi sea in this paper. We consider the Gutzwiller function with an optimization operator[35]:

$$\psi _ { \lambda } = \exp ( \lambda K ) \psi _ { G }, \text{ \quad \ \ } ( 4 )$$

where K is the kinetic part of the Hamiltonian and λ is a variational parameter. This type of wave function is an approximation to the wave function in quantum Monte Carlo method[12, 36, 37]. The operator e λK lowers the energy considerably. We have finite energy gain with this function even in the limit g → 0 due to the kinetic operator K . We show that ψ λ with vanishing g describes a Mott insulator.

In Fig.3 we show the energy per site as a function of U on a 10 × 10 lattice with the periodic boundary conditions. The upper curve shows the energy for the Gutzwiller function and the lower one is for the optimized function ψ λ . It is seen from Fig.3 that the ground-state energy changes the curvature near U ∼ U c . This suggests that there is a transition from a metal to an insulator. We show the parameter g on 10 × 10 lattice in Fig.4. The parameter g vanishes at a critical value U c ∼ 8. The energy for ψ λ is well approximated by a function C/U with a constant C when U is large:

$$\frac { E } { N } \sim - C \frac { t ^ { 2 } } { U } \in - C J. \quad \ \ ( 5 ) \quad \ \dot { \psi } _ { G } \\ \quad \ \ e r \varepsilon$$

This means that the energy gain mainly comes from the exchange interaction which is of the order of 1 /U , showing that the ground state is insulating. The effective interaction in the limit of large U is given by the effective interaction, given by J ∑ 〈 ij 〉 S i · S j with J = 4 t 2 /U , and the three-site interactions[39] if we consider only the nearest-neighbor transfer t . In our calculation, we have C ∼ 3. This means that the ground-state energy per site is approximately given by E/N ∼ -0 75 . J . This value, -0 75 . J , will become better by improving the wave function ψ λ [35]. The inset in Fig.4 exhibits that there is a singularity in g at U /similarequal U c as a function of U . There is a small jump in g . This indicates that the transition is first order.

The Fig.5 shows the momentum distribution function n k :

$$n _ { \mathbf k } = \frac { 1 } { 2 } \sum _ { \sigma } \langle c _ { \mathbf k \sigma } ^ { \dagger } c _ { \mathbf k \sigma } \rangle.$$

There is clearly the gap in n k at the Fermi surface for small U , namely, U ≤ 6. In contrast, for large U being greater than 7, the gap at the Fermi surface disappears. This indicates that the ground state is an insulator for large U . This is consistent with the VMC study in Ref.[33] where the different trial wave function is adopted. Other quantities are also consistent. The critical value of U is consistent; both have given U c ∼ 7 . t The ground-state energy is also well approximated by a curve t 2 /U when U is large beyond the critical value.

FIG. 3: Ground state energy of the 2D Hubbard model per site as a function U at half-filling on 10 × 10 lattice. We set t ′ = -0 2. . The upper curve is for the Gutzwiller function and the lower curve is for ψ λ . The dotted line shows a curve C/U where C &lt; ( 0) is a constant.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-f5305677-bbd7-41a2-b963-274c883b1896.jpg)

Charge-transfer Mott state In this section, we consider the ground state of the three-band d-p model in the halffilled case. The energy unit is given by t dp in this section. The Gutzwiller wave function for the d-p model is ψ G = P G ψ 0 , where P G is the Gutzwiller projection operator for d electrons. We neglect the on-site Coulomb repulsion on the oxygen site because it is not important when the number of p holes is small. ψ 0 is a one-particle wave function given by the Fermi sea. ψ 0 contains the variational parameters ˜ t dp , ˜ t pp , ˜ , ˜ t ′ d /epsilon1 d and ˜ : /epsilon1 p

$$\psi _ { 0 } = \psi _ { 0 } ( \tilde { t } _ { d p }, \tilde { t } _ { p p }, \tilde { t } _ { d } ^ { \prime }, \tilde { \epsilon } _ { d }, \tilde { \epsilon } _ { p } ). \quad \quad ( 7 )$$

In the non-interacting case, ˜ t dp , ˜ t pp and ˜ t ′ d coincide with t dp , t pp and t ′ d in the Hamiltonian, respectively. As we have shown, t ′ d plays an important role to determine the Fermi surface within the d-p model. The Fermi surface is determined by ˜ t dp , ˜ t pp and ˜ t ′ d in the correlated wave function.

The optimized wave function is

$$\psi _ { \lambda } = \exp ( \lambda K ) \psi _ { G },$$

FIG. 4: Gutzwiller parameter g as a function of U at halffilling on 10 × 10 lattice with t ′ = -0 2. . The parameter g almost vanishes at U ∼ 8 as for 6 × 6 lattice. The upper line is for the Gutzwiller function and the lower is for ψ λ . The inset shows g using a logarithmic scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-2ac65c8d-bfdb-4660-b506-6405e7d6a902.jpg)

FIG. 5: Momentum distribution function n k for ψ λ at halffilling on 10 × 10 lattice. We show n k for U = 3 5 6 7 9 and , , , , 10. n k shows the insulating behavior for U ≥ 7.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-d43a1424-470c-42f4-a7e1-f116da09fb9f.jpg)

where K is the kinetic part of the total Hamiltonian H dp and λ is a variational parameter. In general, the band parameters t pp , t ′ d , /epsilon1 d and /epsilon1 p in K are regarded as variational parameters:

$$K = K ( \hat { t } _ { p p }, \hat { t } _ { d } ^ { \prime }, \hat { \epsilon } _ { d }, \hat { \epsilon } _ { p } ). \quad \quad ( 9 ) \quad \text{Th} _ { \text{too} }$$

For simplicity, we take ˆ t pp = ˜ t pp , ˆ t ′ d = ˜ t ′ d , ˆ /epsilon1 d = ˜ /epsilon1 d and ˆ /epsilon1 p = ˜ . /epsilon1 p The energy expectation value is minimized for variational parameters g , ˜ t dp , ˜ t pp , ˜ , ˜ t ′ d /epsilon1 p -˜ /epsilon1 d and λ .

We show the parameter g as a function of ∆ dp for U d = 8, t pp = 0 4 and . t ′ d = 0 on 6 × 6 lattice in Fig.6. g of ψ λ is decreasing rapidly for positive ∆ dp and vanishes at ∆ dp ∼ 0 5 . t dp while that in ψ G decreases gradually as a function of ∆ dp . The figure 8 exhibits the ground-state energy per site ∆ E ≡ E/N -/epsilon1 d as a function of ∆ dp for t pp = 0, t ′ d = 0 and U d = 8. This set of parameters correspond

FIG. 6: Parameters g and λ as a function of ∆ dp at halffilling on 6 × 6 lattice. We used t pp = 0 4, . t ′ d = -0 0 . and U d = 8 (in units of t dp ). The upper line is for the Gutzwiller function and the lower is for ψ λ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-2ce3328b-7cf8-4ecb-bca1-186b74731fd0.jpg)

FIG. 7: Ground state energy of the 2D d-p model as a function of ∆ dp for t pp = 0 0, . t ′ d = -0 0 and . U d = 8 (in units of t dp ) in the half-filled case. The calculations were performed on 6 × 6 lattice. The dotted curve is for the Gutzwiller function, namely λ = 0. The dashed curve indicates a curve given by a constant times 1 / /epsilon1 ( p -/epsilon1 d ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-cdf081a9-b7e9-495e-9008-6cf7aa8a39f3.jpg)

to that for LSCO. We can find that the curvature of the energy, as a function of ∆ dp , is changed near ∆ dp ∼ 2 t dp . This means that the region ∆ dp &gt; 2 t dp is a large-∆ dp region. The energy is well fitted by 1 ∆ / dp shown by the dashed curve in Fig.7 when ∆ dp is greater than the critical value (∆ dp ) c ∼ 2 t dp . A similar behavior is also observed for the other set of parameters such as t pp = 0 4, . t ′ d = -0 2 and . U d = 8 for Bi2212, Tl2201 and Hg1201, as shown in Fig.8. We show the momentum distribution for d electrons defined by n k = 〈 d † kσ d kσ 〉 in Fig.9. This was calculated on 8 × 8 lattice and exhibits the effect of correlation in the localized region.

The results indicate that the energy gain ∆ E mostly comes from the second-order perturbation with the exci-

FIG. 8: Ground state energy of the 2D d-p model as a function of ∆ dp for t pp = 0 4, . t ′ d = -0 2 and . U d = 8 (in units of t dp ) in the half-filled case. The calculations were performed on 6 × 6 lattice. The dotted curve is for the Gutzwiller function with λ = 0. The dashed curve indicates a curve given by a constant times 1 / /epsilon1 ( p -/epsilon1 d ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-ec82d61f-38be-4a4b-9264-05a33a420eb9.jpg)

FIG. 9: Momentum distribution function of d holes in the 2D d-p model for t pp = 0 4, . t ′ d = -0 0 and . U d = 8 (in units of t dp ) in the half-filled case on 8 × 8 lattice. From the top, /epsilon1 p -/epsilon1 d = 6, 4, 2 for ψ λ and the bottom is for the Gutzwiller function with /epsilon1 p -/epsilon1 d = 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0617v1-figure-92a4d391-3451-495f-94d1-0866a6985176.jpg)

tation energy ∆ dp . Thus, in general, the energy gain ∆ E is expanded in terms of ∆ -1 dp when ∆ dp is large:

$$\Delta E \simeq \frac { A _ { 1 } } { \Delta _ { d p } } + \frac { A _ { 2 } } { \Delta _ { d p } ^ { 2 } } + \cdots. \quad \quad ( 1 0 ) \quad \text{$\det$} \\ \text{tra}$$

As seen from Figs.7 and 8, A 1 is negative and A 2 is positive. It is known that the antiferromagnetic exchange interaction works between d electrons on neighboring copper atoms. The coupling J Cu -Cu is given as[40, 41]

$$J _ { C u - C u } = \frac { 4 t _ { d p } ^ { 2 } } { ( \Delta _ { d p } + U _ { d p } ) ^ { 2 } } \left ( \frac { 1 } { U _ { d } } + \frac { 2 } { 2 \Delta _ { d p } + U _ { p } } \right ). \quad ( 1 1 ) \quad \stackrel { \dots } { \text{by a} } },$$

The exchange coupling J Cu -Cu will give the energy gain being proportional to 1 / ∆ 2 dp when the d electrons are antiferromagnetically aligned on copper atoms. Our results show that this contribution is small keeping A 2 positive. For large ∆ dp the energy gain is proportional to 1 / ∆ dp :

$$\Delta E = \frac { E } { N } - \epsilon _ { d } \simeq - C \frac { t _ { d p } ^ { 2 } } { \Delta _ { d p } }, \text{ \quad \ \ } ( 1 2 )$$

for a constant C . This indicates that the ground state is an insulator of charge-transfer type. In our calculation we obtain C ∼ 3.

Summary We have proposed a wave function of Mott insulator based on an optimized Gutzwiller function in strongly correlated electron systems. We have investigated Mott transition at half-filling in the single-band Hubbard model first and generalized it to the three-band d-p model by employing the variational Monte Carlo method. The metal-insulator transition occurs as a result of strong correlation. The wave function in this paper describes a first-order transition from a metal to a Mott insulator.

Our wave function has the form

$$\psi _ { \lambda } = e ^ { \lambda K } \psi _ { G } ( g ), \text{ \quad \ \ } ( 1 3 )$$

where g is the Gutzwiller parameter. The limit g → 0 indicates no double occupancy of d holes in ψ G . In this limit the energy of the single-band Hubbard model is given by that of the strong-coupling limit U /greatermuch t , namely, E ∝ -t 2 /U . This means that ψ λ is an insulator in the limit g → 0. This state indeed becomes stable when U is as large as 7 . t This shows that there is a metal-insulator transition at U = U c ∼ 7 . t There is a singularity in g at U /similarequal U c as a function of U , indicating that the transition is first order. The same discussion also holds for the three-band d-p model except that the cuprates exhibit charge-transfer transition. In the localized region with large ∆ dp , the energy of ψ λ with vanishing g is given by E/N -/epsilon1 d ∼ -t 2 dp / ∆ dp , indicating that ψ λ is a chargetransfer insulator. The stabilization of ψ λ for large U d and ∆ dp shows the existence of a metal-insulator transition in the d-p model. The transition occurs for the band parameters that are suitable for high temperature cuprates. Our result shows that (∆ dp ) c ∼ 2 t dp . If we use t dp /similarequal 1 5eV[20-22], . the charge-transfer insulator has a gap of 3eV.

Finally, we give a discussion on magnetism. The competition between magnetic state and paramagnetic state would depend on band parameters. There would be a transition from a magnetic insulator to a paramagnetic insulator as the band parameters are varied. We expect that t ′ d plays an important role in this transition because t ′ d would play a similar role to the next-nearest neighbor transfer integral t ′ in the single-band Hubbard model.

We express sincere thanks to J. Kondo, K. Yamaji and I. Hase for useful discussions. This work was supported by a Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan. A part of the numerical calculations was

- [1] The Physics of Superconductor (Vol.I and Vol.II) edited by K. H. Bennemann and J. B. Ketterson (SpringerVerlag, Berlin, 2003).
- [2] J. Zaanen, G. A. Sawatzky and J. W. Allen, Phys. Rev. Lett. 55, 418 (1985).
- [3] J. E. Hirsch, E. Y. Loh, D. J. Scalapino and S. Tang: Phys. Rev. B39, 243 (1989).
- [4] R. T. Scalettar, D. J. Scalapino, R. L. Sugar, and S. R. White: Phys. Rev. B44, 770 (1991).
- [5] M. Guerrero, J. E. Gubernatis and S. Zhang: Phys. Rev. B57, 11980 (1998).
- [6] S. Koikegami and K. Yamada: J. Phys. Soc. Jpn. 69, 768 (2000).
- [7] T. Yanagisawa, S. Koike and K. Yamaji: Phys. Rev. B64, 184509 (2001); ibid. B67, 132408 (2003).
- [8] T. Yanagisawa: New J. Physics 10, 023014 (2008).
- [9] T. Yanagisawa, M. Miyazaki and K. Yamaji: J. Phys. Soc. Jpn. 78, 013706 (2009).
- [10] C. Weber, A. Lauchi, F. Mila and T. Giamarchi: Phys. Rev. Lett. 102, 017005 (2009).
- [11] B. Lau, M. Berciu and G. A. Sawatzky: Phys. Rev. Lett. 106, 036401 (2011).
- [12] T. Yanagisawa, M. Miyazaki and K. Yamaji: J. Mod. Phys. 4, 33 (2013).
- [13] C. Gros, R. Joynt, and T. M. Rice: Phys. Rev. B36, 381 (1987).
- [14] H. Yokoyama and H. Shiba: J. Phys. Soc. Jpn. 56, 1490 (1987).
- [15] T. Nakanishi, K. Yamaji and T. Yanagisawa: J. Phys. Soc. Jpn. 66, 294 (1997).
- [16] K. Yamaji, T. Yanagisawa, T. Nakanishi and S. Koike: Physica C 304, 225 (1998).
- [17] K. Yamaji, T. Yanagisawa and S. Koike: Physica B284288, 415 (2000).
- [18] M. Miyazaki, T. Yanagisawa and K. Yamaji: J. Phys, Soc. Jpn. 73, 1643 (2004); ibid. 78, 043706 (2009).
- [19] C. Weber, K. Haule and G. Kotliar: Phys. Rev. B78, 134519 (2008).
- [20] M. S. Hybertsen, M. Schluter and N. E. Christensen: Phys. Rev. B39, 9028 (1989).
- [21] H. Eskes, G. A. Sawatzky and L. F. Feiner: Physica C160, 424 (1989).
- [22] A. K. McMahan, J. F. Annett and R. M. Martin: Phys. Rev. B42, 6268 (1990).
- [23] H. Eskes and G. Sawatzky: Phys. Rev. B43, 119 (1991).
- [24] J. Hubbard: Proc. Roy. Soc. A276, 237 (1963).
- [25] T. Yanagisawa, M. Miyazaki and K. Yamaji: JPS Conf. Proc. 3, 015046 (2014).
- [26] O. K. Andersen, A. I. Liechtenstein, O. Jepsen and F. Paulsen: J. Phys. Chem. Solids 56, 1573 (1995).
- [27] K. McElroy, R. W. Simmonds, J. E. Hoffman, D.-H. Lee, J. Orenstein, H. Eisaki, S. Uchida and J. C. Davis: Nature 422, 592 (2003).
- [28] N. E. Hussey, M. Abdel-Jawad, A. Carrington, A. P. Mackenzie and L. Balicas: Nature 425, 814 (2003).
- [29] D. J. Singh and W. E. Pickett: Physica C203, 193 (1992).
- [30] M. Plate, D. F. Mottershead, I. B. Elfimov, D. C. Peets, R. Liang, D. A. Bonn, W. H. Hardy, S. Chluzbalan, M. Falub, M. Shi, L. Patthey and A. Damascelli: Phys. Rev. Lett. 95, 07001 (2005).
- [31] W. S. Lee et al.: arXiv cond-mat/0606347 (2006).
- [32] H. Yokoyama, Y. Tanaka, M. Ogata and H. Tsuchiura: J. Phys. Soc. Jpn. 73, 1119 (2004).
- [33] H. Yokoyama, M. Ogata and Y. Tanaka: J. Phys. Soc. Jpn. 75, 114706 (2006).
- [34] L. F. Tocchio, F. Becca and C. Gross: Phys. Rev. B83, 195138 (2011).
- [35] T. Yanagisawa, S. Koike and K. Yamaji: J. Phys. Soc. Jpn. 67, 3867 (1998); ibid. 68, 3608 (1999).
- [36] J. E. Hirsch, D. J. Scalapino and R. L. Sugar: Phys. Rev. Lett. 47, 1628 (1981).
- [37] T. Yanagisawa: Phys. Rev. B75, 224503 (2007).
- [38] T. Yanagisawa: New J. Phys. 15, 033012 (2013).
- [39] A. B. Harris and R. V. Range: Phys. Rev. 157, 295 (1967).
- [40] H. Eskes and J. H. Jefferson: Phys. Rev. B48, 9788 (1993).
- [41] A. P. Kampf: Phys. Rep. 249, 219 (1994).
Takashi Yanagisawa Mitake Miyazaki
Mott transition in cuprate high-temperature superconductors
In this study, we investigate the metal-insulator transition of charge
transfer type in high-temperature cuprates. We first show that we must
introduce a new band parameter in the three-band d-p model to reproduce the
Fermi surface of high temperature cuprates such as BSCCO, YBCO and Hg1201. We
present a new wave function of a Mott insulator based on the improved
Gutzwiller function, and show that there is a transition from a metal to a
charge-transfer insulator for such parameters by using the variational Monte
Carlo method. This transition occurs when the level difference
$\Delta_{dp}\equiv \epsilon_p-\epsilon_d$ between d and p orbitals reaches a
critical value $(\Delta_{dp})_c$. The energy gain $\Delta E$, measured from the
limit of large $\Delta_{dp}$, is proportional to $1/\Delta_{dp}$ for
$\Delta_{dp}>(\Delta_{dp})_c$: $\Delta E\propto -t_{dp}^2/\Delta_{dp}$. We
obtain $(\Delta_{dp})_c\simeq 2t_{dp}$ using the realistic band parameters.