УДК 53.01, 538.955, 538.915
Роль четырехспинового обмена в магнитном механизме сверхпроводимости в купратах
Александр В. Шнуренко*
Институт инженерной физики и радиоэлектроники, Сибирский федеральный университет, Свободный, 79, Красноярск, 660041,
Россия
Сергей Г. Овчинников^
Институт физики им. Л.В.Киренского СО РАН, Академгородок, 50/38, Красноярск, 660036,
Россия
Елена И. Шнейдер*
Сибирский государственный аэрокосмический университет,
им. ак. М.Ф.Решетнева, Красноярский рабочий, 31, Красноярск, 660014,
Россия
Получена 18.09.2010, окончательный вариант 25.11.2010, принята к печати 10.12.2010 Поправки к обменному взаимодействию могут изменять температуру сверхпроводящего перехода, поэтому на основе численного самосогласованного решения системы двух уравнений в рамках Ь-1-модели мы исследуем влияние четырехспинового обмена на концентрационную зависимость температуры перехода в сверхпроводящую фазу с йх2_у2 — типом симметрии параметра порядка.
Ключевые слова: сильные электронные корреляции, модель Хаббарда, Ь-1-модель.
Принято считать [1—3], что сильные электронные корреляции (СЭК) играют существенную роль в механизме высокотемпературной сверхпроводимости (ВТСП). Одной из центральных моделей, отражающих СЭК, является модель Хаббарда [4], содержащая в простейшем случае два энергетических параметра: интеграл электронного перескока 1 между ближайшими узлами и энергию кулоновского отталкивания и двух электронов, находящихся на одном узле с противоположными проекциями спиновых моментов (двойки). В режиме СЭК выполняется неравенство И^|1|. При электронных концентрациях п<1 переходят к описанию модели Хаббарда в усеченном гильбертовом пространстве, не содержащем двоек.
При построении эффективного гамильтониана из модели Хаббарда, помимо обычного гейзенберговского обменного слагаемого с J = /ц, который получается при выводе модели, делаем унитарное преобразование, где учитываем поправки в четвертом порядке теории возмущений, порождающие четырехспиновый кольцевой обмен. Необходимость учета четырехспинового обмена видна в экспериментах по нейтронному рассеиванию [5], т.к. это слагаемое дает правильное описание закона дисперсии магнонов в недопированном антиферромагнетике.
Гамильтониан в представлении Гейзенберга для четырехспинового обменного взаимодействия выглядит согласно рис.1.
* [email protected] [email protected] [email protected] © Siberian Federal University. All rights reserved
Рис. 1. Медно-кислородная плоскость атомных орбиталей (с симметрией для Си-йх2_у2 и для 0-2Рх,у), вовлеченных в магнитные взаимодействия. Ле — 4-спиновый кольцевой обмен. Стрелочки показывают спины валентных электронов, вовлеченных в обмен [5]
Н = .7С ^ { (ББ,) (Бк$) + (Б$) (БкБ,) - (БЛ) (бз $) }, (1)
где Б — это перенормированный вектор Я на узлах 1, _].
Эффективный гамильтониан 1-Л-модели [6,7], проецируя на пространство состояний, описывающих нижнюю хаббардовскую подзону, записываем в виде суммы гамильтонианов, описывающих собственное значение энергии. Кинетический член описывает перескоки электронов по узлам решетки с матричным элементом 1 при условии, что на одном узле не может быть больше одного электрона. Гейзенберговский обменный член описывает эффективное обменное взаимодействие электронов на соседних узлах с обменным интегралом
2 ■ г2
J = и . И добавляем гамильтониан с четырехспиновым обменным взаимодействием. Эффективный гамильтониан в представлении Х-операторов:
н = £ (£ -у)ху + £ ^хс0х0с +1 £ JijхГ -х-1
гсгО™Осг _ V"4 J (^аа _ ^ааг^а<т\
_ Ь1 xj +2 ^ "
+ 2 { £ Jc[(xfxJс - хсСхсс) < БкБ > + (хсСхСс - хссхСс) < Бк^ > - (2) 2 <^,к,1> с,а
-(хСсхСс - хСсхССсс) < БSj > ]} ,
где < Sj > и < SkSj > — спиновые перенормированные корреляторы.
Для исследования одновременного влияния магнитных корреляций и четырехспиновых взаимодействий на условия реализации ¿х2-у2-сверхпроводимости воспользуемся методом неприводимых функций Грина, построенных на операторах Хаббарда. Схема расчета обычна для метода неприводимых функций Грина [8, 9], где проводится линеаризация уравнения с введением аномальной и нормальной функций. В функции Грина ((Ху?с0 )) выделяем нормальные средние в приближении Хаббард-1 и аномальные средние с ¿х2_у2 симметрией, составляем уравнение движения в узельном представлении. После перехода в квазиимпульсное представление получается замкнутая система уравнений, аналогичная уравнениям Горькова. Разрешаем эти уравнения и применяем спектральную теорему, обычным образом приходим к перенормированному уравнению самосогласования для сверхпроводящего параметра порядка Ак с учетом четырехспинового обмена:
1 __А Е
Ак = Е (к + р) + 71(к - Р)2Ер ^ (3)
р 1
где ренормированное обменное взаимодействие
1 — х 2
. + (1 — х).с5о1 ) •
(4)
Здесь зависимость от волнового вектора описывается функцией 71 (р) = ^2 егН1Р, вектор — Н1 Ъ,1 соединяет ближайших соседей.
Для сравнения приведем уравнение на энергетическую щель для и модели:
Дк = N Е (.(к + Р) + .(к — Р)2Е th ЕР
2Ер 2кТ
В итоге получено уравнение на температуру сверхпроводящего перехода:
г!
N ■
1 = Е ((соврх — совРу)2 ^ Ер
Еп
2кТс
Запишем также уравнение на химический потенциал:
1
N
= N Е РоаfF
2
(5)
(6)
(7)
Численное решений этих уравнений позволяет вычислить значения зависимости температуры сверхпроводящего перехода от допирования, используя конкретные значения корреляторов [10]. Тем самым определяются ренормировки энергетического спектра за счет магнитных флуктуаций и их влияние на температуру сверхпроводящего перехода.
0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28
X
Рис. 2. Зависимость критической температуры от концентрации для 1-^модели (график пунктир) и для 1-^модели с учетом четырехспинового слагаемого (график сплошная)
Как видно из графиков зависимости (рис. 2), четырехспиновое обменное слагаемое влияет на перенормировку значения критической температуры примерно на 15%. Без учета четырехцентрового слагаемого максимум температуры имел значение, равное 182 К. А с учетом четырехцентрового слагаемого критическая температура равняется 156 К.
Итоги
Получена система самосогласованных уравнений, определяющих температуру сверхпроводящего перехода и химический потенциал.
200
т
150
100
50
Показано, что учет четырехспинового кольцевого слагаемого приводит к перенормировке константы связи, уменьшая ее, что вызывает существенное уменьшение области реализации сверхпроводящей фазы с типом симметрии параметра порядка.
Критическая температура при учете кольцевого обмена понижается по сравнению с температурой, которая определяется обычным гейзенберговским обменом.
Список литературы
[1] P.W.Anderson, Science, 235(1987), 1196.
[2] N.M.Plakida, High-Temperature Superconductivity, Springer, Berlin, 1995.
[3] С.Г.Овчинников, УФН, 167(1997), 1043.
[4] J.Hubbard, Proc. Roy. Soc. A, 276(1963), 238.
[5] R.Coldea et al., Physical Review Letters, 86(2001), №23.
[6] N.M.Plakida, V.S.Oudovenko, Phys. Rev. B, 59(1999).
[7] Yu.A.Izuymov, B.M.Letfullov, J. Phys: condens. Matter., 3(1991), 5373.
[8] N.M.Plakida, V.Yu.Yushankhay, I.V.Stasyuk, Physica C, 787(1989), 162-164.
[9] V.Yu.Yushankhay, N.M.Plakida, P.Kalinay, Physica C, 174(1991), 401.
[10] M.M.Korshunov, S.G.Ovchinnikov, The European Physical Journal B, 57(2007).
Role of Four-spin Exchange in a Magnetic Mechanism of Superconductivity in Cuprates
Alexander V. Shnurenko Sergey G. Ovchinnikov Elena I. Shneyder
Amendments to the exchange interaction can alter the superconducting transition temperature, so on the basis of numerical self-consistent solution of two equations in the framework of the tJ model, we investigate the influence of the four spin exchange on the concentration dependence of the transition temperature to the superconducting phase with d(X2 — y2) — the type of symmetry of the order.
Keywords: strong electron correlation, Hubbard model, t-J model.