КРАТКИЕ СООБЩЕНИЯ ТЕОРЕТИЧЕСКАЯ И МАТЕМАТИЧЕСКАЯ ФИЗИКА
Расчет потенциала обменного взаимодействия при наличии флуктуирующего поля неупорядоченной системы электронных спинов
М. К. Бахнян1, A.M. Савченко1,0, Б. И. Садовников1,2,6
1 Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра
квантовой статистики и теории поля. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2. 2Институт нанотехнологий микроэлектроники РАН. Россия, 119991, Москва, Ленинский пр-т, д. 32А. E-mail: а a.m.savchenko@gmail.com, ь sadovnikov@phys.msu.ru
Статья поступила 15.06.2011, подписана в печать 27.07.2011
В работе для описания взаимодействий в приближении сильной связи используется контурное представление операторов. На основе такого представления найден потенциал обменного взаимодействия, показано, что эффективное обменное взаимодействие уже не определяется только парными корреляциями, а формируется более сложным образом с учетом четырехчастичного взаимодействия. Ключевые слова: потенциал обменного взаимодействия, контурное представление, функция Грина. УДК: 538.91. PACS: 71.10.-w.
Для описания взаимодействий в приближении сильной связи мы воспользуемся контурным представлением операторов, т. е. будем рассматривать нашу систему как систему «взаимодействующих» контуров, на которых определены операторы А„, фа{х), 0а1/к, 07К [1-4]. Перейдем теперь в этом формализме к рассмотрению обменного взаимодействия между электронами. Рассмотрим взаимодействие пары электронов в данной системе. Они взаимодействуют с неупорядоченной спиновой подсистемой и самосогласованно взаимодействуют с кристаллической решеткой и между собой через кулоновское отталкивание.
Ориентация их спинов может быть параллельной или антипараллельной, чему соответствуют симметричная и антисимметричная волновая функция. Поскольку в данном случае система описывается совокупностью взаимодействующих контуров, то соответствующие волновые функции реализуются в контурном представлении:
Ф,(хь ЩХ1, о = о ЪВ(Т)Ф1(Х1, 0) +
+ №/(х1,0ТгВ(ГШх/,0))сг3, (1) Фа(хь цщ, 0 = 0 ТгВ(Г№/(*/, 0) -
^(ф/(х1,0ТгВ(Г)ф1(х/,0))Ста, (2)
где (та — спиновые волновые функции, описывающие состояние пары электронов с полным спином
й + Я/.
Функции «свободных» электронов Фг(хг,0 удовлетворяют системе самосогласованных уравнений в контурном представлении
1
'А
5с t
д(ш) dt
2тгСк
DK dx" gph(x)gph(x") х
х (ф{х)ф{х")фЦх", Г')фк{х", ПЪВ(ТШх, о) -
х (Ф1(х", Пфк{х", П7гВ(ТШх, 0). (3)
Надо отметить, что интеграл здесь берется по еДГ) — малой петле, «приклеенной» к контуру Г:
Ф,(*,*|Г) = /тг^Гехр
()dxv Аъ
■ф-Хх, t)
ФгМГ|) + Тг[£>2,Фг(лс,*|Г)] =
Уравнения (3) можно легко замкнуть, если предположить, что 1-й электрон взаимодействует с конечным числом электронов. На языке контуров это означает, что внутри петли еДГ) сосредоточено малое число электронов, взаимодействующих с электроном под номером /. Если же внутрь петли вообще попадает лишь один электрон, то мы получаем систему двух уравнений.
Используя выражения (1), (2), можно, как обычно, построить эффективное кулоновское взаимодействие пары электронов, которое распадается на прямое кулоновское взаимодействие (2(хг,х;|Г) и так называемое обменное взаимодействие /(х;,х;|Г).
Средние значения потенциалов (2(х;,х;|Г) и /(х;,х;|Г) тогда будут соответственно равны
(Я(х1, х,|Г)) = N йх1 йх;(ф1{х{)ф1{х1)У{х1 - X/) X
хТгВ(ГШ*/№/(*/)>, (4)
йХ1 йх!(ф*{х1)ф*{х^У{х1 -X,) X хТгВСГШх/ЩхО). (5)
Средние выражения (4), (5) следует понимать как средние хронологические произведения операторных
{/ (х;,х;-|Г)) = N
волновых функций. Поскольку система предполагается пространственно неоднородной, то средние значения потенциалов (<Э(л:/,л:,|Г)), {/(ж;,*;-|Г)) приходятся на выделенную элементарную ячейку с индексом к. Размер ячейки 4 при этом должен быть много меньше характерного размера неоднородности Ьк: 1К/Ьк <с 1.
Тогда
(QK(xl,xi\T)) = Nf2(N)
dxi d.Xj ip*K(xi)ip*K(xj) x
х У(Х1 - х))ф1к{хдфы,(хд{1гВ(Г„Л) Тг б(Гд- Л). (6)
Контуры ТХьХ1, непрерывной деформацией могут быть стянуты в точки, тогда выражение (6) с точностью до коэффициента переходит в выражение, впервые полученное Гайтлером и Лондоном [5].
Проводя усреднение по конфигурациям ячеек, получим
При этом потенциал обменного взаимодействия можно выразить через функции Грина, а затем преобразовать к следующему виду:
(41К2(л;г,*/|Г)) = —МММ) х
(N2
1)2
DKl(x>>)DK2(x>>)
-кь-кгбГ*..*,.;*
х G*iKi (xi\xj)v(xi
(2тгск)(2тгск) dXidXiX
Xj)GjiKJxi\xj) х
х exp
Ar — 1
dx'l pl0u{Ki] +
dx"p20vlK2]
"Kl (*,")
«2 (■*,")
х (Тгтс11В(Г)тс11В(к1))(Тгтс12В(Г)тс12В(к2)), (7)
,€¡¿2 = 1
0!1К1 №) = (хдФич (*/)),
Г/ .,. ДА';1А',:) = (ХгЦ-кг (*/))■
Отметим, что экспоненты в формуле (7) возникают из-за контурного Фурье-представления («-представления) функций Грина.
Подчеркнем физический смысл этих функций. Неупорядоченная система спинов не является «замерзшей». В ней даже при нуле температур всегда существуют флуктуации. Флуктуации электронных спинов связаны с их локальным «переворотом», который осуществляется за характерное время гЧт, т' ~ с, где т — время свободного пробега электронов.
Распределение электронных спинов в пространстве определяется пространственным распределением зарядов. Действительно, за время «переворота» спина выделенного электрона, остальные электроны успевают перераспределиться в пространстве, т. е. происходит изменение распределения зарядовой плотности системы. Вследствие локальных «переворотов» происходит и изменение спиновой плотности в данной точке пространства х. На этом основании можно сделать вывод, что флуктуации спиновой и зарядовой плотностей связаны друг с другом.
Неупорядоченная система электронных спинов описывается потенциалом калибровочного поля А„, кото-
рый преобразуется как в координатном, так и в спиновом пространстве. Это поле ответственно также и за флуктуации и поэтому его можно представить в виде Av = A°v + dv. Тогда функции ршЛк<] описывают поле импульса, который электрон приобретает за счет взаимодействия с флуктуирующим полем av. При этом корреляционная функция полей ршЛк<]> PjOv'[Kj\ ПР°~ порциональна (Тга^к^аЦк,]), причем коэффициент пропорциональности выражается через магнитную восприимчивость системы, так как флуктуации спиновой плотности связаны с появлением в выделенном микроскопическом объеме неравновесной намагниченности пг(к).
Из выражения (7) следует, что эффективное обменное взаимодействие уже не определяется только парными корреляциями, а формируется более сложным образом. Действительно, можно представить себе, что два связанных в ячейке к\ контура Г электрона обмениваются импульсом с парой связанных электронов в ячейке к2 того же контура. Таким образом возникает так называемое четырехчастичное взаимодействие, при котором электроны взаимодействуют как внутри ячеек, так и между ячейками. Необходимо подчеркнуть, что на важность многочастичных, в частности четырехчастич-ных, взаимодействий в магнитных системах указывал в своих основополагающих работах и в последующих обобщающих исследованиях H.H. Боголюбов [6-8]. Им было показано, что представление о коллективных электронных возбуждениях как куперовских парах является концептуально важным первым приближением, но оно не может быть абсолютизировано, и, вообще говоря, мы имеем дело с коллективными возбуждениями всего электронного конденсата как с целостной системой.
В работе [9] впервые обсуждалась возможность спаривания электронов с неравными по модулю импульсами. Такое размытие импульсов объяснялось многочастичными корреляциями. В данном же случае мы имеем аналогичную ситуацию, однако за соответствующее размытие импульсов электронов ответственны флуктуации спиновой и зарядовой плотностей.
Перейдем к вычислению интегралов в (7). Для этого перепишем (7) в более компактной форме:
(4к* (*ь*/|Г)) =
N2
(N2 - I)2
öNf2(N) х
-к* er.
DK(x>>)DK(x>>) ^
(2жск)(2жск) dXidX>X
х JKK. exp
dx'l Pqv[k] +
dx'lp0l/[K*
Ф'/)
лг-i
x (Тгта1В(Г)та1В(к))(Тгта2В(Г)та2В(к*)),
<2| ,<22 = 1
dxi dxj G*iK(xi\xj)V{xi - Xj)GjjK*
Здесь мы заменили к\, к2 на к, к*, так как только в этом случае импульс получаемой пары электронов,
ТЕОРЕТИЧЕСКАЯ И МАТЕМАТИЧЕСКАЯ ФИЗИКА
133
находящихся в данной ячейке «2, совпадает с импульсом, испускаемым электронами в ячейке «1.
Ячейки к, к* отличаются только направлением обхода. В случае, если размер ячейки 1К много меньше размера контура Г, интегралы в (7) удается вычислить. Тогда имеем
= {4кОр(«#/2(^)(ТгВ(Г))2{сЬ(7г4Л))Р(к,Ро), (8) где р(к) — функция распределения конфигурации ячеек; р0 — средний импульс, передаваемый между ячейками к, к*; р0^2ж/Ьк; (1кк*)р{к) = /0; (ТгВ(Г)) следует понимать как среднее по всем конфигурациям и ориентациям контура Г в пространстве. Параметр (сЬ(тг1кро)} является сложным функционалом от (ТгВ(Г)). При ХфО, С(ТгВ(Г)) имеет численное значение такое, что при 1к/ьк <С 1
(гс) = ( Тг
dxx27(x\Tx) / (Тг
dxJ{x\Tx)
((гс) — обменный радиус корреляции).
Функция (ТгВ(Г)) может быть получена как решение полевого уравнения и зависит от среднего размера контура Г.
Таким образом, выражение (8) можно также понимать как потенциал взаимодействия между электронными парами различных ячеек, т. е. ({/кк.играет роль обменного взаимодействия между ячейками.
Если функция (ТгВ(Г)) является осциллирующей в пространстве, это означает, что в системе могут возникнуть кластеры с дальним магнитным порядком: ферромагнитный, если /< О, и антиферромагнитный, если / > 0. л
Если (ТгВ(Г)) — медленно меняющаяся функция координат, это означает, что на фоне неупорядоченных спинов может возникнуть длинноволновая магнитная структура с модулированным по величине магнитным моментом.
Так как (ТгВ(Г)) ^ 1, то при (ТгВ(Г)) -¥ 1 обменное взаимодействие может оказаться существенно выше, чем в обычных магнитных системах, т. е. обычное обменное взаимодействие может быть усилено дополнительно обменным взаимодействием между ячейками.
При (ТгВ(Г)) -¥ 0 обменное взаимодействие между электронами отсутствует, т. е. они сильно экранированы разупорядоченной системой спинов и между ними существует только кулоновское отталкивание. В этом случае дальний магнитный порядок будет отсутствовать.
Итак, физический смысл параметра ((JKK* состоит в том, что он описывает обменное взаимодействие пары электронов, которые в свою очередь взаимодействуют с эффективным полем остальных электронов и с флуктуирующим полем неупорядоченной системы электронных спинов. Как уже отмечалось, из-за флуктуаций такое парное взаимодействие эквивалентно четырехчастичному. Если {гс) -¥ оо, то это означает, что на каждый электронный спин действует постоянное самосогласованное поле. Но тогда этот случай соответствует отсутствию всяких флуктуаций.
Список литературы
1. Sadovnikov B.I., Savchenko A.M. 11 Phys. A. 1999. 271. P. 411.
2. Sadovnikova M.B., Savchenko A.M., Scarpetta G. 11 Phys. Lett. A. 2000. 274. P. 236.
3. Savchenko M.A., Stefanovich A.V. Fluctuational Superconductivity of Magnetic Systems. Springer-Verlag, 1990.
4. Савченко A.M., Сорокина E.M. I ! Вестн. Моск. ун-та. Физ. Астрон. 2009. № 3. С. 12.
5. Гайтлер В. Элементарная квантовая механика. М., 1948.
6. Боголюбов H.H. Избр. тр.: В 3 т. Т. 3. Киев, 1971.
7. Боголюбов H.H. // Успехи физ. наук. 1959. 67, № 4. С. 549.
8. Боголюбов H.H., Толмачев В.В., Ширков Д.В. Новый метод в теории сверхпроводимости. М., 1958.
9. Heine V., Pippard A.B. I ! Philos. Mag. 1958. 3. P. 1046.
The calculation of the exchange interaction potential in a fluctuating field of disordered system of electron spins
M.K. Bakhnyan1 , A.M. Savchenko1B.I. Sadovnikov1
1 Department of Quantum Statistics and Field Theory, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
2 Institute of Nanotechnologies in Microelectronics, Russian Academy of Sciences, Leninsky ave. 32A, Moscow 119991, Russia.
E-mail: a a.m.savchenko@gmail.com, bsadovnikov@phys.msu.ru.
To describe the interaction in the strong coupling approximation the contour representation of the operators is under consideration. On the basis of such a representation the exchange interaction potential is found. It is shown, that effective exchange interaction is not determined by pairing correlations only, but forms in a more complex way on a basis of four-particles interaction.
Keywords: exchange interaction potential, contour representation, Green's function. PACS: 73.43.Nq. Received 16 June 2011.
English version: Moscow University Physics Bulletin 1(2012).
Сведения об авторах
1. Бахнян Михаил Константинович — аспирант; тел. (495) 939-12-90, e-mail: bakhnyan@grnail.com.
2. Савченко Александр Максимович — докт. физ.-мат. наук, доцент, ст. преподаватель; e-mail: a.m.savchenko@gmail.com.
3. Садовников Борис Иосифович — докт. физ.-мат. наук, профессор; (495) 932-80-10, e-mail: sadovnikov@phys.msu.ru.