А. В. Марченко, В. П. Волков, П. П. Серегин
ИЗМЕНЕНИЕ ЭЛЕКТРОННОЙ ПЛОТНОСТИ В УЗЛАХ КРИСТАЛЛИЧЕСКОЙ РЕШЕТКИ ПРИ СВЕРХПРОВОДЯЩЕМ ФАЗОВОМ ПЕРЕХОДЕ
Представлены результаты расчета изменения электронной плотности в узлах одномерной кристаллической решетки Кронига— Пенни при сверхпроводящем фазовом переходе: переход от нормальной к сверхпроводящей фазе сопровождается возрастанием электронной плотности в центре элементарной ячейки, что согласуется с данными мессбауэровской спектроскопии для М1,85Сво: 15Си04, Ьа1,8з8го,15Сп04, УБа2Спз07-х, УВа2Сп408,
Т12Ва2Сап-1Сип02п+4 (п = 1, 2), Вг28г2Сап-1Сип02п+4 (п = 1, 2) и ЩВа2Сап-1Сип02п+2 (п = 1, 2, 3).
Согласно микроскопической теории Бардина—Купера—Шриффера (БКШ) в основе явления сверхпроводимости лежат эффекты возникновения при температуре ниже критической температуры Тс связанных состояний электронов (так называемых куперовских пар) и последующего образования бозе-конденсата куперовских пар [1]. Иными словами, переход от нормального к сверхпроводящему состоянию есть переход от системы электронов, описываемой блоховскими волновыми функциями, к системе электронов, описываемой единой когерентной волновой функцией, так что распределение электронной плотности в кристаллической решетке сверхпроводника должно различаться при температурах выше и ниже температуры Тс.
Мессбауэровская спектроскопия позволяет обнаружить изменение электронной плотности в узле кристалла, содержащего мессбауэровский зонд, при переводе кристалла из нормального в сверхпроводящее состояние путем измерения центра тяжести мессбауэровского спектра зонда при температуре выше и ниже температуры перехода вещества в сверхпроводящее состояние. Является очевидным, что сравнение экспериментальных (данные мессбауэровской спектроскопии) и теоретических величин изменения электронной плотности при сверхпроводящем фазовом переходе может служить критерием выбора моделей, описывающих явление сверхпроводимости.
В настоящей работе представлены результаты расчета изменения электронной плотности в узлах одномерной кристаллической решетки Кронига— Пенни при сверхпроводящем фазовом переходе, а также результаты экспериментального определения изменения электронной плотности в узлах кристаллических решеток металлоксидов меди Кё1.85Се0.15Си04, Ьа1 85 8г015Си04, УВа2Сиэ07-х, УВа2Сщ08, ТЬВа2Сап-1Сип02п+4 (п = 1, 2), В12Бг2Сап-1Сип02п+4 (п = 1, 2) и Н§Ва2Сап-1Сип02п+2 (п = 1, 2, 3) с использованием данных эмиссионной мессбауэровской спектроскопии на изотопах 67Си(672п) и 67Оа(672п).
Расчет изменения электронной плотности в узлах кристаллической решетки Кронига—Пенни
Электронная плотность фермионного газа в нормальном состоянии представлялась в виде
Pf(x) =Z I V(x)|2. (!)
k=-kf
где k — волновое число; оно находится из периодических граничных условий V[ x + N (a + b)] = v( x) и равно
i n 2n 1
k =--------+---------1, (2)
a + b N (a + b)
здесь l = 1, 2, ..., N), kf — волновое число, отвечающее энергии Ферми 8/,
V(x) = uk(x)exp(ikx) (3)
— одноэлектронные блоховские волновые функции, причем для области 0< x < а справедлива зависимость
uk = A exp[i(a - k) x] + B exp[-i(a + k) x], (4)
где a — есть функция от k, и она определяет энергию электрона е, а коэффициенты А и В связаны условием uk(0) = uk(a) и условием нормировки
а
N j1 uk(x)|2 =!-
0
В модели Кронига—Пенни задаются ширина ямы а, ширина потенциального барьера между ямами b, высота потенциального барьера Vo, число атомов в
дг _ \2m(V0 -E) J32ab п п
решетке N, в = Л--------0-------------------------------------------------------> , и-> P, где m — масса электрона, а Р
V n2 2
характеризует энергию связи электрона с решеткой. Функция a(k) определяется из уравнения
sin aa ...
P-+ cosaa = cos ka, (5)
aa
закон дисперсии имеет вид
«(к) = , (6)
2т
а для u*(x) справедливо выражение
ик (x) = A{ cos(a - к)x + ух cos(a + к)x + y2 sin(a + k)x + +z'[sin(a- к)x + y2 cos(a + к)x -yx sin(a + к)x]},
(7)
где
- sin(a + к)a - sin(a - к)a + sin 2aa
y =-------------------------------------; (8)
2[1 - cos(a + к) a]
1 - cos(a + к )a - cos(a- к )a + cos2a
Y2 =------- ----- ------- ---- -----------------------------. (9)
2[1 - cos(a + к) a]
В сверхпроводящем состоянии электронная плотность считалась равной сумме электронных плотностей бозе-конденсата и оставшегося ферми-газа
к=к/-Зк
Рв = ncPc (x) + Z k(x)f , (10)
к=-кf+Зк
где Nc — число куперовских пар, 2Зк = 2П- —---------область в пространстве вол-
a N
новых векторов около уровня Ферми, в которой происходит спаривание,
Na
pC (x) = 11 Y(q, x) |2 dq — зарядовая плотность куперовской пары, состоящей из
0
двух электронов с координатами q и x, ¥(q,x) — ее волновая функция, которая находится из решения уравнения Шредингера
[ H0 (q) + H0 (x) + U (q, x)]Y (q, x) = EY (q, x), (11)
где U(q,x) — энергия взаимодействия двух электронов куперовской пары, H0(x)
— гамильтониан электрона, который не взаимодействует со вторым электроном куперовской пары
Ho(x)k (x) =| £(к) 1 кк (x), (12)
где |£(к)| — закон дисперсии, <рк (x) — одноэлектронная волновая функция в
модели Кронига—Пенни, причем согласно теории ферми-жидкости Ландау отсчет энергии ведется от уровня Ферми 8/ (т. е. от значения s при к = к/). Волновая функция куперовской пары представлялась в виде
Y fe x) = Z скк (q)k- к (x), (13)
к
и, подставив (13) в (11), получали уравнение для коэффициентов разложения ск
волновой функции Y(q,x) по собственным функциям гамильтониана (12). В
этом уравнении имеются матричные элементы межэлектронного парного взаимодействия
икк. = -g, при (к/- ¿к) < к,к < (к/ + ¿к)
и
икк. = 0, если к и к лежат вне указанного интервала.
Иными словами, предполагалось, что эффективное взаимодействие между электронами носит характер притяжения только в узкой области 5к вблизи уровня Ферми.
Решая (13) относительно ск , имеем
&
ск =■
2|£(к )|-Е
к = к ,+Зк
■> = I с, ,
к—ку-Зк
и, подставив (14) в (13), получаем уравнение для энергии Е:
к =ку +Зк
1 = I --------------1-------•
8 к'=кг-Зк 2|т(k )|- Е
(14)
(15)
(16)
В случае положительных 8 существует решение с Е < 0, т. е. существует
т 1 Д - - Ш связанное состояние с т = 1. Для плотности состояний, равной
2п
де(к )
дк
энергия связанного состояния
Е =-----
2Зк
да(к )
1 - ехр
2п да (к ) дк
8Ыа дк к=к/
к-ке ’
к-к,
(17)
где
д^( к)
дк
к-к,
есть значение производной одночастичной энергии по к при
к=ку, и оно связано с законом дисперсии соотношением
да(к )
дк
а коэффициенты ск определялись как
ск =■
к= ке
к - кг
(18)
(19)
2 | т(к) | -Е
Таким образом, для расчета волновых функций ¥(^,х) куперовских пар:
• е(к) определялось по соотношению (6), значения а2(к) брались из (5), а значения к — из равенства (2);
да (к )
учитывалось, что
дк
к=к.
Ь да(к)
= —а
т
дк
к = ке
(здесь а, есть значение
а при к = к,);
• |4(к)| определялись из (18);
• ск определялись из выражения (14);
• коэффициенты ск считались для к в интервале (к, - Зк) < к < (к, + Зк) ;
Е =
учитывалось, что 25к
да(к)
где в =
1 - ехр 4п
2п да(к )
дк к=к/
дк
2Ь2 агЗк
да(к)
к=ке
т <¡1 - ехр
Р-±
8
дк
да( к )
Ыааг дк
к=ке
Задавались следующие параметры: число атомов в периодической структуре N, постоянная решетки а, безразмерный параметр взаимодействия электрона и решетки Р, отношение уровня Ферми к матричному элементу межэлек-
тронного парного взаимодействия , диапазон энергетических уровней, в ко-
8
тором следует учитывать притяжение электронов друг к другу 5к. Нормировка ¥(д,х) проводилась из условия
МаМа
| 11 ¥(д, х) |2 dqdx = 1 с2к = 2.
(20)
о о
На рис.1 приведены зависимости рАх) и Рв(х) для случая N =500, а =1 А,
Р = 4, — = 3,2 и Зк = — . Видно, что переход от нормальной к сверхпрово-
8 а N
дящей фазе сопровождается возрастанием зарядовой плотности в центре элементарной ячейки, где методом мессбауэровской спектроскопии измеряется зарядовая плотность |¥(0)|2. Естественно, на периферии элементарной ячейки зарядовая плотность уменьшается при переходе от нормальной к сверхпроводящей фазе.
<=><
6
к
Л
Расстояние от центра элементарной ячейки, А
Рис. 1. Зависимости электронной плотности от расстояния от центра элементарной ячейки для обычной (штриховая линия)
и сверхпроводящей (сплошная линия) фаз (М =500, а =1 А, Р = 4, — = 3,2 и Зк = )
а N
На рис. 2 приведена зависимость относительной разности электронной плотности в центре элементарной ячейки для сверхпроводящей и нормальной
р (0) — р (0) £г £г
фаз —--------- --- от — . Видно, что уменьшение — (т. е. переход к фазам с
Р-(0) 8 8
более высокими значениями Тс) сопровождается возрастанием
Рв (0) -Р- (0) Р- (0)
Экспериментальное определение изменения электронной плотности
Как отмечалось, изменение электронной плотности в узле кристалла, содержащего мессбауэровский зонд, при переводе кристалла из нормального в сверхпроводящее состояние возможно обнаружить путем измерения центра тяжести мессбауэровского спектра зонда при температуре выше и ниже температуры перехода вещества в сверхпроводящее состояние. Объясняется это тем, что температурная зависимость центра тяжести мессбауэровского спектра Б (при условии постоянства давления Р) определяется выражением [2]:
(21)
где первый член описывает температурную зависимость изомерного сдвига при постоянном объеме, и третий член описывает зависимость допплеровского
сдвига второго порядка Б от температуры. Появление члена вызвано
изменением электронной плотности кристалла, и именно этот эффект ожидается при переходе кристалла из нормального в сверхпроводящее состояние.
I
£
ё
Рис. 2. Зависимость относительной разности электронной плотности в центре элементарной ячейки для сверхпроводящей и нормальной фаз
Рв (0) - Р- (0)
р-(0)
от отношения уровня Ферми к матричному элементу
межэлектронного парного взаимодеиствия
2п 80
— N =500, а =1 А, Р = 4, и 8к =--)
8
а N
8
Обнаружение куперовских пар методом мессбауэровской спектроскопии должно быть наиболее благоприятным для случая высокотемпературных сверхпроводников (имеющих минимальный масштаб куперовской корреляции), если
используется зонд, для которого > 10 (здесь Д — максимально достижимая
20
разность изомерных сдвигов мессбауэровских спектров в обычной и сверхпроводящей фазах, 0= —--------естественная ширина ядерного уровня, то — среднее вре-
Тс
мя жизни ядерного уровня). Выбор объектов для исследования должен также учитывать необходимость введения в узлы решетки мессбауэровского зонда.
Все эти условия выполняются для случая мессбауэровского зонда 672и в решетках металлоксидов меди при использовании эмиссионного варианта мес-сбауэровской спектроскопии на изотопах 67Си(672и) (возможно введение мате-
67
ринского изотопа Си в процессе синтеза в узлы меди, так что дочерний изотоп 672и также оказывается в медном узле решетки [3]) и 670а(672и) (для металло-ксидов меди, содержащих в своем составе редкоземельные металлы (РЗМ), материнский изотоп 670а оказывается в узлах РЗМ, так что дочерний изотоп 672и также оказывается в РЗМ узле решетки [4]), при этом возникает возможность исследования пространственной неоднородности электронной плотности, создаваемой бозе-конденсатом куперовских пар. Описание приготовления безно-сительных радиоактивных изотопов 67Си и 670а можно найти в работе [5], а детали синтеза мессбауэровских источников металлоксидов меди, легированных 67Си, описаны в работе [6].
Алгоритм обработки экспериментальных данных был следующим: определялось положение центра тяжести спектра Б для данного соединения (сверхпроводящего и контрольного, причем в качестве контрольных объектов, для которых не наблюдалось перехода в сверхпроводящее состояние, были выбраны соединения, полученные путем термообработки соответствующих сверхпроводящих образцов) при различных температурах; строились температурные зависимости Б для контрольного и сверхпроводящего соединения; для контрольных соединений согласно формуле (21) определялись температуры Дебая и температурная зависимость дебаевского сдвига Б(Т) аппроксимировалась полиномом; для сверхпроводящего соединения для каждой температуры, при которой измерялся спектр, определялась разность I = [Б - Б] (здесь Б — расчетное значение допплеровского сдвига), строилась температурная зависимость I, которая аппроксимировалась полиномом, и путем экстраполяции этой зависимости на Т ^ 0 К определялось предельное значение изомерного сдвига 1с.
В качестве типичного примера на рис. 3 приведены температурные зависимости Б для узлов Си(1), Си(2) и У в решетках сверхпроводящего и контрольного соединения Т12Ва2Са2Си3010. Зависимости Б(Т) существенно разли-чаются для контрольного и сверхпроводящих материалов: для контрольного образца зависимости для узлов Си(1) и Си(2) описываются теоретической зависимостью дебаевского сдвига Б(Т) с 9 = 410 К во всей области температур, тогда как аналогичные зависимости для сверхпроводящего образца в области температур Т > Тс описываются такой же функцией Б(Т), а при Т < Тс зависимость Б(Т) становится более резкой, хотя при переходе через Тс для всех соединений резких скачков в вели-
чине Б не наблюдается. Значение 1о тем больше, чем выше Тс, а для соединений, включающих позицию иттрия и две позиции меди, величина 1о — наибольшая для узлов Си(2), существенно меньше для узлов Си(1) и минимальна для узлов У. На рис. 4 приведены зависимости 1с = Г(Тс-1) — для узлов меди Си(1), Си(2) и У — для случая экспоненциальной аппроксимирующей функции.
г,
й
А
Й
I
¡/1
т,к
Рис. 3. Температурные зависимости центра тяжести Б мессбауэровских спектров 672п в узлах Си(1) и Си(2) контрольного и сверхпроводящего соединения Т12Ва2Са2Си30]0: 1 — узлы Си(1) и Си(2) контрольного соединения; 2 — узлы Си(1) сверхпроводящего соединения; 3 — узлы Си(2) сверхпроводящего соединения.
Сплошная линия соответствует теоретической температурной зависимости Б для случая допплеровского сдвига второго порядка при 9 = 400 К
Рис. 4. Зависимость 1с от Тс 1 для узлов Си(1), Си(2) и Я в соединениях К^.85Се0.15Си04, Ьа1.85Бг0.15Си04, УВа2Сиз07-х, УВа2Си408, Т12Ва2Сап-1Сип0211+4 (п = 1, 2),
В12 Бг2Сап-1 Сип02п+4 (п = 1, 2) и Ы§Ва2Сап-1Сип02п+2 (п = 1, 2, 3)
Величина изомерного сдвига мессбауэровских спектров 67/и непосредст-
67
венно связана с изменением электронной плотности на ядрах /и
I = аДр ,
(22)
(здесь Др — разность релятивистских электронных плотностей на ядрах 67/и в нормальном и сверхпроводящем образцах, а — калибровочная постоянная изотопа 67/и), причем величины 10 характеризуют электронную плотность, создаваемую бозе-конденсатом в условиях, когда все электроны проводимости находятся в составе бозе-конденсата куперовских пар.
Следовательно, зависимость на рис. 4 свидетельствует, что с ростом Тс возрастает величина Др, и это отражает факт возрастания электронной плотности на ядрах 67/и при переходе от несверхпроводящей к сверхпроводящей фазе. Зависимость Др от Тс-1 возникает вследствие того, что изменение электронной плотности при сверхпроводящем переходе должно зависеть от размера куперовских пар. Действительно, стандартная корреляционная длина («размер» куперовской пары при Т ^ 0 К) в теории БКШ определяется как ~ Тс-1, и поэтому на рис. 4, по сути, приведена зависимость 10 от ^о. Видно, что эта зависимость носит экспоненциальный характер для всех узлов:
10 (мкм/с) = [6,5 ± 0,4] ехр
10 (мкм/с) = [9,8 ± 0,4] ехр
10 (мкм/с) = [3,2 ± 0,4] ехр
56( К ) Тс( К )
33( К )'
Тс( К)
35( К )
Тс( К)
— для узлов Си(1),
— для узлов Си(2),
— для узлов Я.
Оказалось, что максимально возможное изменение электронной плотности на ядрах 67/и при сверхпроводящем переходе составляет для узлов Си(1) Др = 0,10(1) ат. ед. и для узлов Си(2) Др = 0,18(1) ат. ед., и эти величины соответствуют минимально возможным «размерам» куперовской пары <^0шт в под-
решетках Си(1) и Си(2). Существование минимального «размера» куперовских пар связано с невозможностью существования куперовских пар с расстоянии-ем между компонентами, меньшим некоторой критической длины. Если считать, что в области малых имеет место соотношение Др2 ~ ^о-3, то получим £0шт ~ 5Д_ — для подрешетки Си(1), и ^0Ш1П ~ 2,5А — для подрешетки Си(2). Эти
величины согласуются с общепринятыми значениями ^0Ш1П ~ (0,5-30)А.
Таким образом, установлено, что для сверхпроводников в области Т>Тс температурная зависимость Б определяется допплеровским сдвигом второго порядка, тогда как в области Т < Тс на величину Б влияет также процесс появления куперовских пар и их бозе-конденсация. Понижение температуры увеличивает влияние указанного процесса на величину Б, поскольку с понижением температуры увеличивается доля бозе-конденсата.
В общем случае температурная зависимость эффективной плотности сверхтекучих электронов р(Т) может быть записана выражением [1]
Р(Т) = 1 - в (23)
^ Ч Г (ехрв +1)2 ' '
где Ер =-^-----энергия Ферми, т — эффективная масса электрона, к — волно-
кр_
2т
вой вектор, кр — значение волнового вектора на поверхности Ферми, Ек — энергия к-состояния, в = , к0 — постоянная Больцмана.
КТ
С другой стороны, следовало ожидать, что р ~ — (здесь 1Т и 10 — изо-
^0
мерные сдвиги при температуре Т и Т = О К). Поэтому на рис. 5 приведена тео-
к Т I--------
ретическая зависимость р от параметра х = 1,76-^ [здесь А = 3, О6ко^/Тс(Тс - Т)
— энергетическая щель в спектре элементарных возбуждений сверхпроводника], взятая из работы [1], вместе с нашими данными по зависимости — от парато
метра х. Видно, что имеется удовлетворительное согласие расчетных и экспериментальных величин.
I к Т
Рис. 5. Зависимость у- от параметра х =1,76—.
Сплошной кривой показана теоретическая зависимость эффективной плотности сверхтекучих электронов от параметра х [1]. Точками представлены данные для узлов Си(1), Си(2) и Я в соединениях Ка1>85Се0,15СиО4, Ьа1858г015СиО4, УБа2СизО7-х, УВа2СщО8, ТЬБа2Сап-1СипО2п+4 (п = 1, 2), Б12Бг2Сап-1СипО2п+4 (п = 1, 2)
и Ы§Ва2Сап-1СипО 2п+2 (п = 1, 2, 3)
Является неожиданным удовлетворительное согласие расчетных и экспериментальных температурных зависимостей эффективной плотности сверхтекучих электронов высокотемпературных сверхпроводников. В литературе об-
суждалась проблема применимости или неприменимости модели БКШ для описания свойств высокотемпературных сверхпроводников, и поэтому следует с осторожностью относиться к обнаруженному нами согласию между теоретической и экспериментальной зависимостями эффективной плотности сверхтекучих электронов от параметра х и к возможности использовать теорию БКШ в ее немодифицированном виде для описания высокотемпературной сверхпроводимости. Это согласие следует рассматривать как доказательство того, что процессы образования куперовских пар и их бозе-конденсация должны быть обязательными в любой теории высокотемпературной сверхпроводимости.
* * *
Для металлоксидов меди установлены соотношения между микроскопическими параметрами эксперимента (изменением электронной плотности) и теорией (стандартной корреляционной длиной): изменение электронной плотности в узлах кристаллической решетки данной подрешетки при сверхпроводящем фазовом переходе тем больше, чем меньше стандартная корреляционная длина для данной подрешетки, и в сверхпроводящей фазе существует минимально возможное значение стандартной корреляционной длины <^ошт, причем для соединений, включающих несколько подрешеток, величина <^0шт различна для
различных подрешеток. Полученные результаты находятся в согласии с данными по расчету изменения электронной плотности в узлах кристаллической решетки Кронига—Пенни.
БИБЛИОГРАФИЧЕСКИЕ ССЫЛКИ
1. Шриффер Дж. Теория сверхпроводимости. М., 1970.
2. Серегин П. П. Физические основы мессбауэровской спектроскопии. СПб., 2002.
3. Бондаревский С. И., Мастеров В. Ф., Серегин П. П. Центры меди в полупроводниковой и сверхпроводниковой фазах УВа2Си3О7-х // ФТТ. 1990. Т. 32. Вып. 10. С. 3150-3154.
4. Насрединов Ф. С., Серегин Н. П., Серегин П. П., Бондаревский С. И. Мессбауэров-ское исследование двухэлектронной акцепторной примеси цинка в кремнии // ФТП. 2000. Т. 34. Вып. 3. С. 275-277.
5. Бондаревский С. И., Еремин В. В., Серегин Н. П. Новые экологически безопасные и высокоэффективные методы сепарации радионуклидов широкого спектра применения: наука, техника, ядерная медицина // Труды международной конференции «Экологические проблемы и пути их решения в XXI веке: образование, наука, техника». СПб., 2000. С. 82-84.
6. Серегин Н. П., Насрединов Ф. С., Бондаревский С. И. Исследование высокотемпературных сверхпроводников методом эмиссионной мессбауэровской спектроскопии // Труды международной конференции «Экологические проблемы и пути их решения в XXI веке: образование, наука, техника». СПб., 2000. С. 94-96.
CHANGING ELECTRONIC DENSITY IN SITES OF CRYSTALLINE LATTICE UNDER SUPERCONDUCTING PHASE TRANSITION
Presented are results of the calculation of the changes electron density in sites of one dimensional crystalline lattice Kronig-Penney during superconducting phase transition: the transition from normal to superconducting phase is accompanied by electronic density growth in the elementary cell centre that is in agreement with the data of Mossbauer spectroscopy for Nd1.85Ce0,15CuO4, La1.85Sr0,15CuO4, YBa2Cu3O7-x, YBa2Cu4O8, Tl2Ba2Can_1CunO2n+4 (n = 1, 2), Bi2Sr2Can-1CunO2n+4 (n = 1, 2) u HgBa2Can-1CunO2n+2 (n = 1, 2, 3).