CC BY
18
Bulletin of PFUR Series Mathematics. Information Sciences. Physics. No 2 (2). 2010. Pp. 161-166
UDC 510.676, 519.7
Mean Field Green Function Solution of the Two-Band Hubbard Model in Cuprates
Gh. Adam, S. Adam
Laboratory of Information Technologies Joint Institute for Nuclear Research Joliot-Curie 6, 141980 Dubna, Moscow region, Russia Horia Hulubei National Institute for Physics and Nuclear Engineering (IFIN-HH), 407 Atomistilor, Magurele - Bucharest, 077125, Romania
The paper discusses the generalized mean field solution of the Green function matrix of the effective two-band two-dimensional Hubbard model of the high-Tc superconductivity in cuprates, as recently modified to include appropriate boundary conditions at zero doping. Two main results are found. (i) Hybridization of normal state energy levels preserves the center of gravity of the unhybridized levels. (ii) Hybridization of superconducting state energy levels displaces the center of gravity of the unhybridized normal levels. The whole spectrum is displaced towards lower frequencies.
Key words and phrases: high-Tc superconductivity, cuprates, two-band Hubbard model, Green function, mean field approximation.
1. Basic principles of the theoretical description of
cuprates
The high critical temperature superconductivity in cuprates is still a puzzle of the today solid state physics, in spite of the unprecedented wave of interest and number of publications. The two-band two-dimensional Hubbard model [1-3] provides a description of it based on four basic principles:
(1) Deciding role of the experiment.
The derivation of reliable experimental data on various cuprate properties asks for manufacturing high quality samples, performing high-precision measurements by adequate experimental methods. The next section summarizes the wealth of data provided by experiment.
(2) Hierarchical ordering of the interactions inferred from data.
(3) Derivation of the simplest model Hamiltonian following from the Weiss principle,
1.e., hierarchical implementation of the various interactions into the model.
(4) Mathematical solution by right quantum statistical methods which secure rigorous implementation of the existing physical symmetries and observe the principles of mathematical consistency and simplicity.
2. Experimental input to the theoretical model
There are five kinds of experimental data which are essential for the derivation of a consistent theoretical model.
2.1.Crystal structure characterization points to the occurrence of layered ternary perovskite structures, with an overwhelming contribution to the superconducting pairing coming from the CuO2 planes.
Consequence: An effective two-dimensional (2D) model for the CuO2 plane is requested.
Received 28th November, 2009.
The authors acknowledge partial financial support from project CEX 05-D11-67 (COMCIP), the Hulubei-Meshcheryakov Programme, JINR Orders 324/28.05.09, p. 30, and 328/29.05.09, p. 8, and the LIT-JINR theme 05-6-1060/2005-2010.
2.2. Existence of the Fermi surface was undoubtedly evidenced, first by 2D-ACAR positron spectroscopy at the beginning of the nineties [4-6] and then by ARPES and optical methods [7].
Consequence: The energy bands lying at or nearest to the Fermi level are to be retained in the model.
2.3. The charge-transfer insulator nature of the cuprates occurs with the following relationship between the energy band parameters, U > A > W.
Consequence 1: The hybridization results in the Zhang-Rice singlet subband.
Consequence 2: The Zhang-Rice singlet and the upper Hubbard subbands enter the simplest model.
Consequence 3: Since A ~ 2^, the model is to be developed and solved in the strong correlation limit.
2.4. Tightly bound electrons in the metallic state.
Consequence 1: Occurrence of low density hopping conduction consisting of both fermion and boson (singlet) carriers.
Consequence 2: This asks for the Hubbard operator description [8] of the system states.
2.5. Occurrence of cuprate families which are characterized by specific stoichio-metric reference structures, doped with either holes or electrons.
Consequence 1: The doping parameter Ô is essential in the theoretical description of the cuprates.
Consequence 2: The (Ô,T) phase diagrams arise which have to be accounted for.
3. Input: abstractions, concepts, facts
Besides the straightforward inferences following from the experiment, a number of additional input items need consideration.
3.1. Abstraction of the physical CuO2 plane with doped electron states by a doped effective spin lattice. This is done by a one-to-one mapping from the copper sites inside the CuO2 plane to the spins of the effective spin lattice. There are four possible spin states at each lattice site i in the effective spin lattice: |0) (vacuum), |a) = |f) and |<r) = |4,) (single particle spin states inside the hole subband), and |2) = |f|) (singlet state in the singlet subband).
The spin lattice constants equal ax, ay, the CuO2 lattice constants.
The effective spin lattice is characterized by antiferromagnetic spin ordering at zero doping.
The doping of the electron states inside the CuO2 plane is equivalent to the creation of defects inside the spin lattice, by spin vacancies and/or singlet states.
The occurrence of a hopping conductivity inside the spin lattice is a consequence of the doping. The hopping conductivity may consist either of single spin hopping (fermionic conductivity) or singlet hopping (bosonic conductivity).
3.2. Concept : The global description of the hopping conduction around a spin lattice site can be done by means of the Hubbard 1-forms [9].
3.3. Fact : The hopping induced energy correction effects are finite over the whole range of the doping parameter Ô [10], hence appropriate boundary conditions are to be imposed in the limit of vanishing doping.
4. Model Hamiltonian
The originally derived model Hamiltonian of the 2D Hubbard model [1] was rewritten in terms of Hubbard 1-forms [9] and then put in locally manifest Hermitian form [10]
H = Ho + pHh = ^ (ho,i + phh>i), hlti = hoti, hjl,i = hhii, (1)
i
with the single particle and hopping contributions at the spin lattice site i given respectively by
ho,i = E\ +E2X22 ,
^11 V"^ ( <rO,O<r 0<t,<t0\ , ^22 \ ^ ( 2<t,<T2 a2,2a\ .
T1A - T1,i 1 +
M1 / a0,0a 0a,a0\ . ^22 (
kh,i = "y [T1 , i - T1, i )
') + (
+ ^ E
2 a ,0a _ Oa,2a N + / a Of 2 _ a2,a 0
1,i T1,i J + {T1,i T1,i
The Hubbard 1-forms define the hopping conduction neighbourhood of z-th lattice site,
T1,i = u%m
Xm (TaP,Tn\ ^ = _ % m i i ' m ; = - ' M •
The Hubbard operators at the spin lattice site i are defined as = |ia){i31, where |a) and |3) denote the initial, respectively final spin states.
In (1), p implements boundary conditions at vanishing doping. This Hamiltonian serves to the definition of the quasi-particle spectrum and superconducting properties.
5. Rigorous mean field solution of Green function
matrix
The retarded and advanced Green function (GF) matrices are defined in terms of Nambu operators in the space-time (r, t) representation [11]. Using the equations of motion of the involved Heisenberg operators, differential equations of motion are derived. In this representation, splittings of the higher order correlation functions are also done.
In the dual space-energy (r, w) representation, which is obtained from the (r, t)-representation by appropriate Fourier transforms, the differential equations of motion are transformed in algebraic equations of motion. Analytic extensions in the complex energy plane result in a unique GF in the complex plane. Calculations of statistical averages are done in this representation by use of spectral theorems.
Performing one more Fourier transform from the space variable r to the momentum variable q we get the momentum-energy (q, w)-representation. Within this representation we get: compact functional GF expressions, equations for the energy spectra, statistical average calculations from spectral theorems, spectral distributions inside the Brillouin zone.
The Green function matrices of the model define space-time correlations for the four-component Nambu column operator [2,3] Xi(J = (XJ2 X0cr X2j Xf0)T and its adjoint operator Xja = (X2<7 Xj0 Xj2 X0a) (the superscript T denotes the transposition).
To get a standard eigenvalue problem for the spectrum, the (q, w)-representation of the generalized mean field approximation (GMFA) GF solution is written in terms of the energy matrix [10],
~ , r ~ i -1
G (q,w)=X1/2[Iw -£a (q)] X1/2 , (2)
£(q) = r1/2X(q)x-1/2; x = aX*,xl}), (q) = ^ eiq-r^' Aij J ; r^- = r,- - r ; Ai<T = ({[Xi<T, H],X]ja}) .
The x matrix is diagonal,
* =(x X) ' X=C° x°i) , °=(0 0) '
where *2 = (nia) = (nis) and = (n^) = (n^) = 1 — X2 denote spin and site independent averages of particle number operators. In terms of the doping rate 5, it results that in the hole-doped cuprates, *2 = <5, = 1 — 8, while in the electron-doped cuprates, = X2 = 1 —
The energy matrix can be written in the block .structure form,
^ (q) = I ,t' iqiL) (3)
Ä(q))t -(E*(q))T/
With the choice [10] p = X1X2, the normal and anomalous 2 x 2 energy matrices are obtained as follows
K(q)= f2 _ ^ ,
<1„ (q) = ( -2 ^ ^ T21
(q) ^-T21 -VXÖ2 2aT1 • X2 ) ,
W2 = W2(q) = (E1 + A) + [«22 + d22(q)] • X1, W1 = W1 (q) = E1 + [«22 + dn(q)] • X2, U21 = W21 (q) = [«21 + ^21(q)], T2 = T2(q) = ^2261 + (1 - S)a&2(q) + ^1&s(q), T1 = T1 (q) = £n&1 + (1 - S)^&2(q) + 5&&3(q), T21 = T21(q) = £2161 + (1 - 5)6&2(q) + 56&3(q),
where expressions for the r.h.s. terms have been reported in [10].
6. GMFA energy spectrum
All the 16 matrix elements of the Green function (2) share a same denominator, V = P(q, w) = Iw — (q), with the following monic bi-quadratic dependence in w:
P(q,w) = (w2 — uw + v)(uj2 + uw + u), (4)
where u and v are spin-independent quantities detailed in [1°]. The zeros of P(q, w), provide the GMFA energy spectrum of the system.. In the normal state, the energy spectrum is given by the roots of the second order equation w2 — uw + v = 0 solved previously [1] in the mean field approximation.
There are two main results established from the study of the energy spectrum of the energy matrix (3).
(i) Hybridization of normal state energy levels preserves the center of gravity of the unhybridized levels.
(ii) Hybridization of superconducting state energy levels displaces the center of gravity of the unhybridized normal levels. The whole spectrum is displaced towards lower frequencies.
6.1. Energy spectrum of the normal state
The spectral equation w2 — uOw + vO = 0 has the coefficients given respectively by Uo = W2 + W1 and VO = W2W1 — X1X2IW2112.
Its solutions are n2 = w2 + 3O for the upper subband and = w1 — 3O for the lower subband respectively. The small parameter 3 satisfies the equation /02 + DO/O — X1X2|w2112 = 0 ; Do = W2 — W1, wherefrom
3o = (Do^o/2)/(l + V1 + Vo) ; Vo = (4X1X2|w21|2)/D2 •
6.2. Energy spectrum of the superconducting state
From the factorization (4), it results that Q3 = —Q2 and Q4 = —Q1, with Q2 and Q1 obtained from the secular equation w2 — uw + v = 0, where
V2 = + V ^ V = \! V0 + V ^ <k> = V — Vo = p/(v + Vo ).
u2 = u0 + 2<5v + ^ ^ u = \JV0 + 2<k> + ^ ^ = u — u0 = (2<5i> + ^)/(u + u0).
The hybridization yields Q2 = Q2> + 5u/2 + ^1 and Q1 = Q01 + bu/2 — /31, where the small parameter ft satisfies the equation + D1fj1 + 5v/2 — -0/4 = 0; D1 = — Q0,
wherefrom
ßi = {Dim/2)/( 1 + ^1 + m) ; VI = — 2<to)/D2.
References
1. Plakida N. M., Hayn R., Richard J.-L. Two-Band Singlet-Hole Model for the Copper Oxide Plane // Phys. Rev. B. — 1995. — Vol. 51. — Pp. 16599-16607.
2. Plakida N. M. Superconductivity in the Two-Band Singlet-Hole Model for Copper Oxyde // Physica C. — 1997. — Pp. 1737-1738.
3. Exchange and Spin-Fluctuation Mechanisms of Superconductivity in Cuprates / N. M. Plakida, L. Anton, S. Adam, Gh. Adam // JETP. — 2003. — Vol. 97. — Pp. 331-342.
4. Direct Observation of Fermi Surface in YBa2Cu3O7— / H. Haghighi, J. H. Kaiser, S. Rayner et al. // Phys. Rev. Lett. — 1991. — Vol. 67. — Pp. 382-385.
5. Window Least Squares Method Applied to Statistical Noise Smoothing of Positron Annihilation Data / Gh. Adam, S. Adam, B. Barbiellini et al. // Nucl. Instr. Meth. in Phys. Res. — 1993. — Vol. A337. — Pp. 188-203.
6. Fermi Surface Ridge at Second and Third Umklapp Positron Annihilations in YBa2Cu3O7— / Gh. Adam, S. Adam, B. Barbiellini et al. // Solid State Commun. — 1993. — Vol. 88. — Pp. 739-742.
7. Damascelli A., Hussain Z., Shen Z.-X. Angle-Resolved Photoemission Studies of the Cuprate Superconductors // Rev. Mod. Phys. — 2003. — Vol. 75. — Pp. 473541.
8. Hubbard J. Electron Correlations in Narrow Energy Bands. II. The Degenerate Band Case // Proc. R. Soc. Lond. A. — 1964. — Vol. 277. — Pp. 237-259.
9. Adam Gh., Adam S. Rigorous Derivation of the Mean Field Green Functions of the Two-Band Hubbard Model of Superconductivity // J. Phys. A: Math. Theor. — 2007. — Vol. 40. — Pp. 11205-11219.
10. Adam Gh., Adam S. Finiteness of the Hopping Induced Energy Corrections in Cuprates // Romanian J. Phys. — 2009. — Vol. 54. — Pp. 797-814.
11. Zubarev D. N. Two-time Green Functions in Statistical Physics // Sov. Phys. Uspekhi. — 1960. — Vol. 3. — Pp. 320-366.
УДК 510.676, 519.7
Функции Грина для решения двухзонной модели Хаббарда в приближении среднего поля в купратах
Г. Адам, С. Адам
Лаборатория информационных технологий Объединенный институт ядерных исследований ул. Жолио-Кюри д.6, Дубна, Московская область, 141980, Россия Национальный институт физики и ядерной технологии им. Х. Хулубея Магурель-Бухарест, 077125, Румыния
В работе рассматривается использование техники уравнения движения функций Грина для решения эффективной двухзонной модели Хаббарда высокотемпературной сверхпроводимости в приближении среднего поля в купратах, которая недавно была модифицирована авторами для включения надлежащих граничных условий при нулевом легировании. Получены два новых результата. 1. Гибридизация нормальных уровней состояния энергии сохраняет центр тяжести негибридизированных уровней. 2. Гибридизация энергетических уровней состояния сверхпроводимости смещает центр тяжести негибридизированных уровней. Весь спектр смещается в сторону более низких частот.
Ключевые слова: высокотемпературная сверхпроводимость, купраты, двухзонная модель Хаббарда, функции Грина, приближение среднего поля.
