Вариационный метод оценки угловых зависимостей верхнего критического поля многослойных сверхпроводящих наноструктур

_Доклады БГУИР_

2004 октябрь-ДЕКАБРЬ № 4

ЭЛЕКТРОНИКА

УДК 612.382

ВАРИАЦИОННЫЙ МЕТОД ОЦЕНКИ УГЛОВЫХ ЗАВИСИМОСТЕЙ ВЕРХНЕГО КРИТИЧЕСКОГО ПОЛЯ МНОГОСЛОЙНЫХ СВЕРХПРОВОДЯЩИХ НАНОСТРУКТУР

В Н. КУШНИР

Белорусский национальный технический университет пр. Ф. Скорины, 65, Минск, 220027, Беларусь

Поступила в редакцию 9 апреля 2004

Дано вариационное решение граничной задачи для уравнений Гинзбурга-Ландау, описывающих критическое состояние сверхпроводящих многослойных наноструктур. Построены угловые зависимости верхнего критического магнитного поля при различных температурах. Определена эффективная размерность и параметры зародыша сверхпроводящей фазы в зависимости от ориентации вектора напряженности магнитного поля.

Ключевые слова: уравнения Гинзбурга-Ландау, сверхпроводимость, многослойные наноструктуры, критическое магнитное поле.

Введение

Многослойные тонкопленочные структуры на основе чередующихся слоев сверхпроводника (S) и нормального метала (N) (СНС) являются предметом постоянных и интенсивных исследований. Интерес к подобным структурам обусловлен следующими причинами. Во-первых, достаточная простота приготовления дает свободу в выборе материальных параметров СНС — сверхпроводящего и нормального металлов, толщин сверхпроводящих и нормальных слоев, а также количества последних. Это влечет за собой и разнообразие в физических характеристиках СНС — величинах критической температуры, критического тока, параметра анизотропии, сил пиннинга на границе S/N, температурных зависимостей верхних критических магнитных полей и т.п. [1,2]. Во-вторых, простые транспортные измерения указанных величин позволяют исследовать фундаментальные характеристики сверхпроводящего состояния — вид волновой функции сверхпроводящей фазы, эффективную размерность структуры, коэффициент прозрачности на границе S/N, свойства симметрии параметра порядка, зародышеобразование сверхпроводящей фазы и т.п. В-третьих, многослойные анизотропные структуры на основе обычных сверхпроводников (Nb, V, Та) являются и тестовыми системами для моделирования процессов, происходящих в высокотемпературных сверхпроводниках (ВТСП), имеющих естественную слоистую структуру [3, 4].

Одна из фундаментальных характеристик сверхпроводников второго рода — верхнее критическое магнитное поле Нс2.

Для СНС достаточно подробно изучены температурные зависимости Нс2(Т) [5-16]. Основной результат, который был при этом получен, — явление размерного кроссовера. Оно состоит в том, что при температурах ниже температуры кроссовера Т* сверхрешетка СНС в параллельном внешнем магнитном поле ведет себя подобно сверхпроводящей тонкой пленке (20-поведение) с зависимостью Hc2\\(T)~(l-T/Tc*)V2, а при Т>Т* — подобно однородному бес-

конечному образцу (ЗБ-поведение) с зависимостью Яс2 (7 )~( 1-777'6) (Тс — критическая температура).

Полной критической магнитной характеристикой СНС является зависимость Нс2 (Т.\ 9) [21, 22], где 9— угол между вектором напряженности внешнего магнитного поля и поверхностью плоской структуры. Однако моделирование экспериментальных зависимостей Нс2(9) [20, 23, 24] для СНС оказывается гораздо более сложной задачей [29, 30], чем для однородных анизотропных сверхпроводников [25-28].

Для однородных анизотропных сверхпроводников были легко определены и экспериментально полностью подтверждены предельные характеристики Нс2(9): предел тонкой пленки [21, 25]; предел полубесконечной сверхпроводящей среды [21, 26]; предел бесконечной сверхпроводящей среды [21,4]. В первом и втором случаях зависимость Нс2(9) характеризуется пиком (или разрывом производной дНс2(9)/д9) при 9=0. В третьем случае функция Нс2(в) в точке 9=0 является гладкой. Существенно, что 20- или ЗБ-характер угловых зависимостей находится во "взаимно однозначном" соответствии с 20- или ЗБ-характером температурных зависимостей верхнего критического поля. В [27] для малых углов 9 был рассмотрен общий случай сверхпроводящего образца конечной толщины. Оказалось, что вид и аналитические свойства зависимости Но2 (9) связаны с положением зародыша сверхпроводящей фазы, определяемым его взаимодействием с границами образца.

Исследование характеристик Нс2(9) СНС осуществлялось в [29, 30] посредством численного решения соответствующего уравнения в частных производных для волновой функции (ВФ), описывающей сверхпроводящую фазу бесконечной сверхрешетки. Оказалось, что в данном случае 2Б или ЗБ-характер угловых зависимостей соответствует таковому для температурных зависимостей только при 9=0, что проявляется в наличии или отсутствии разрыва производной 8Нс2(9=0)/д9. Данные численные расчеты мало пригодны для моделирования экспериментальных данных в силу их громоздкости, а также наличия большого количества подгоночных параметров.

Цель данной работы — исследовать с помощью наиболее простого вариационного расчета [20, 21] угловые зависимости Нс2(9) СНС конечной толщины. В качестве вариационного используется функционал Гинзбурга-Ландау (ГЛ). Предполагается, что данный вариационный расчет может описывать экспериментальные данные для значений 9< 1 и для температур, не очень близких к Т0.

Формулировка метода

Рассматривается сверхпроводящая структура, состоящая из m плоских чередующихся слоев сверхпроводящего (например, Nb) и нормального (например, Си, Pd) металлов, находящаяся во внешнем магнитном поле. Выберем систему координат следующим образом: ось OZ направлена перпендикулярно плоскостям слоев, координатная плоскость XOY параллельна слоям и совпадает с плоскостью симметрии СНС. Данная сверхпроводящая структура предполагается неограниченной в направлениях ОХ и OY и имеющей толщину L в направлении OZ.

Для внешнего магнитного поля H(0, Н cos в. Н sin#), близкого к критическому, с соответствующим ему вектор-потенциалом А(г)(Н z cos9-Н у sin ft 0, 0), волновая функция (ВФ) Гинзбурга-Ландау (ГЛ) может быть представлена в виде , z). В этом случае

функционал ГЛ может быть записан в форме [21, 22, 19]

F = JJdydz-

ТГ +МгЧ ~4

Здесь z0=k/H0; r/{z; Т) — ступенчатая функция:

Г 1, ге /„,

ц

где ^ЛТ) — "корреляционные длины", характеризующие быстроту изменения ВФ внутри сверхпроводящего и нормального слоев; 15, 1;\< — области, занятые сверхпроводящими и нормальными слоями соответственно. Пространственные переменные масштабированы на корре-

2

ляционную длину магнитное поле масштабировано на величину Нс23=Фо/2 к^ ,

где Фо — квант магнитного потока.

Варьированием (1) получим следующее уравнение в частных производных для цАут)'-

Г я2 я2 1

1 -л

ду2 дг2

При этом функция \\г

т V" /| з,_>+00 ' ^^

= 0,

2

а также условиям сшивания на границах раздела г, (/=1,..., т-1) сверхпроводящих и нормальных слоев [1]

& & '

= % -0),

где рг=р, ЦгЩ, если слева от границы находится нормальный слой, и р, = р 1 , д, щ 1 — в противном случае; параметры р, <у даются микроскопической теорией.

Параметр г0 определяется варьированием (1), что приводит к интегральному соотношению

Цф^^о -

или, в кратких обозначениях, г0 = ^СО80 — ^¡пб^. То же соотношение (6) следует и из граничных условий для магнитного поля Н(г) [22].

Задача (2)-(4) вместе с (6) есть задача на собственные значения Н, наибольшее из которых и есть верхнее критическое поле Нс2(в). Очевидно, численное решение данной задачи очень громоздко и никак не годится для быстрой подгонки экспериментальных данных. В то же время простейший вариационный принцип, заключающийся в выборе пробной ВФ с разделенными переменными:

МНу,*) = /

оказывается вполне эффективным для описания зависимости Нс2(в) хотя бы для достаточно малых значений в. Для его обоснования приведем следующие соображения. При малых значениях в перекрестный член

Г

в

роль перекрестного члена уменьшается и из-за наличия плоскостей раздела, создающих анизотропию. Это легко увидеть формально, включив условия сшивания (5) по методу Лагранжа в (1). При этом в уравнении (3) (которое можно рассматривать как уравнение для квантового осциллятора) появляется дополнительный "потенциал"

Гг

02 г-

где С\, С2 — множители Лагранжа. Следовательно, с увеличением количества границ вклад (8) в (3) можно все с большим основанием считать малым.

Таким образом, подставляя (7) в (1) и варьируя, получим следующие уравнения для /(у),

«(г):

dy

-+#2sin2е у2-2Н2 (zcosG - z0 )sinQ у + в" /(у) = О

с граничными условиями

/(±оо) = 0

dz'

- T|(z) + Н2 (zcosG - z0 )2 - 2H2J7Sin9 (zcosG - z0)+s' g(z) = 0

с граничными условиями

dz

(±1/2) = 0

и условиями сшивания g'(zi+0) = pig'(zi+0), s(zi + О) = q¡ g(zi - О).

Здесь введены обозначения: у = ах \dyyf(y), z = Ъ1 jdzzg2(z), a = jdyf2(y), b = ^dzg2(z),

в' = — \dy a J

df(y) . dy

\2

+ H2 y2sm2Qf2

'dg(z/ . dz ,

- Л(z) g2(z) + H2 (zcosG - zj g2(z) .

В силу

Выражение

Решение задачи (10), (11), соответствующее минимальному значению é, известно [31]:

и

V /1^

<к-

Учитывая (18)—(20), вместо (15) имеем:

л

-у\{г;Т)+ Н2со/0

Таким образом, вместо граничной задачи (4) для уравнения в частных производных (3) имеем граничную задачу (13) для обыкновенного дифференциального уравнения (21).

« юл 1 ¿Нс 2(0) Во введении было сказано о важности величины а(0) =--—--крутизны

Нс2Ф)

характеристики Нс2(Э ). В зависимости от того, является данная величина разрывной или нет при 6=0, можно сделать заключение об эффективной размерности сверхпроводящей фазы в окрестности данного значения 9. В рамках модели ГЛ (1) легко получить следующее общее выражение для «(8):

(9) ( ^ ~~ г)соз9 - (у - у^пе] - г^тб + (у - .у)со89] ^ ^[(г - г)со80 - (у - ^^пб]2^

При использовании пробной ВФ (7) выражение (22) значительно упрощается:

сое©

а

Яс2соэ20 ((г - + |зт0|

Здесь

В частности, при 9=0 имеем

Из (23), (25) видно, что в описанном приближении относительное падение характеристики Нс2( 9) определяется исключительно параметром локализации а = — г)2^

ВФ по переменной ъ. Кроме того, очевидно, что в рамках данного приближения ни при каких температурах не может быть получено ЗБ-поведение сверхпроводящей фазы, за исключением непосредственной окрестности Тс, где <7—>оо.

Пример и обсуждение результатов

В качестве примера рассчитывались критические характеристики Нс2(в) для 3-слойной (1 би-слой) и 17-слойной (8 би-слоев) СНС. 3-слойная структура взята как аналог тонкой пленки однородного сверхпроводника.

Приняты следующие значения толщин ¿4 и сверхпроводящего и нормального слоев: с1\:=0.5. с/<= 1.05. 2,л-(())=22, ;(7'с). Температурные зависимости корреляционных длин взяты в соответствии с приближением ГЛ в виде (1 7/7'к) "2, (Г/Гс)~1/2 [17] с их экст-

9

раполяцией на область низких температур. Здесь Т5 — критическая температура сверхпроводящего материала.

На рис. 1 представлены температурные зависимости #с2ц(7)= =Но2(Т,в=0) для 1-й 8-бислойной структур. Видно, для 1-бислойной структуры зависимость Нс2\\(Т) имеет 2Б-характер практически для всех температур. Для 8-бислойной структуры имеем типичный для СНС график Нс2ц(7) с 20-30 кроссовером при температуре Г*=0,86ГС.

На рис. 2,а,б представлены графики зависимостей Нс2(в) при низких температурах и температурах, близких к Т*, соответственно. Видно, что для 3-слойной структуры зависимость Нс2(в) точно воспроизводится формулой Тинкхама, полученной для пленок однородного сверхпроводника, при любых углах и температурах. Этот результат в рамках предложенной модели очевиден, поскольку верхнему критическому полю соответствует единственное для всех в и Г значение параметра гп =0, соответствующее симметричной волновой функции.

Т/Т

Рис. 1. Зависимости верхнего параллельного критического поля от температуры для 1- и 8-бислойной СНС

-г-

"Т"

0.7

0.6

1

X" 0.5

0.4

0.3 -0.2

1-бислойная СНС

- О - 8-бислойная СНС, гО = -5.75

- • - 8-бислойная СНС, г0 = 0 -Приближение тонкой пленки

0,0 0,2 0,4 0,6

1,0 1,2 1,4 1,е

в

I"

1-бислойная СНС

- О - 8-бислойная СНС, 20 = -5.75

- • - 8-бислойная СНС, 20 = 0 -1' бёаеёает ёЗ ом а ё Т ёэ\ её

0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4

в (гаф

а б

Рис. 2. Графики зависимости Нс2(6) для 1-бислойной СНС (светлые квадраты) и для 8-бислойной СНС (штриховая линия с кружками): а — '/'=0,24 Т3; б — 7,=0,425Г3. Светлыми кружками обозначен участок кривой 20-поведения, черными — участок ЗЭ-поведения сверхпроводящей фазы. Сплошная кривая построена в соответствии с приближением Тинкхама

Для 8-бислойной структуры, как видно из рис. 2,а, 2Б-поведение характерно при низких температурах для значений углов вплоть до 1 рад; при 0>1 кривая Яс2(0) существенно отклоняется от соответствующей кривой для тонкой пленки. Данный факт легко объясняется в терминах параметра локализации О ВФ сверхпроводящей фазы и положения сверхпроводящего зародыша, определяемого параметром г0. Для указанной температуры и при 9<1 размер зародыша вдоль оси ОХ, определяемый параметром а, меньше или порядка с15: наибольшему значению магнитного поля соответствует значение параметра г„=±{Ь/2-с1: -с13/2). т.е. зародыш

локализуется в крайнем сверхпроводящем слое структуры (рис. 3,а). При этой же температуре и при в>\ энергетически выгодным оказывается симметричное решение для функции g(z), соответствующее значению параметра г„=() (рис. 3,6, в). При этом сверхпроводящая фаза образуется, по меньшей мере, в двух центральных 8-слоях, т.е. 2Б-поведения здесь нет.

9И1дт,

1\

ш

Ж

И1 " 1К1

-6 -5 -4 -3 -2 -1 0 1 2 3 4 5

г/Л

9(^)19т

Г

г

>

Л

ч

-5-4-3-2-10 1 2 3 4 5

214,

Рис. 3. Графики 2-составляющсй ^(г) ВФ сверхпроводящей фазы при 7-0,242 Г8: а-9=0,3; б — 0=1,2; в — 0=1,5. Вертикальные линии обозначают границы раздела слоев

б

а

в

При движении к более высоким температурам интервал 2Б-поведения по углам в сужается для 8-слойной СНС. При этом аппроксимация соответствующего участка кривой формулой Тинкхама становится хуже. Однако физическая трактовка, приведенная для низких температур, остается справедливой. Более того, изменение симметрии ^-составляющей ВФ в данном случае проявляется более отчетливо: на рис. 2,6 видим излом кривой Нс2(&) при 6^0,3.

В заключение автор выражает признательность проф. С.Л. Прищепе за полезные обсуждения.

Работа финансировалась в рамках межреспубликанской межвузовской программы "На-ноэлектроника".

THE VARIATIONAL METHOD OF THE UPPER CRITICAL FIELD ANGULAR DEPENDENCES ESTIMATIONS FOR A SUPERCONDUCTING MULTILAYER NANOSTRUCTURES

V.N. KUSHNIR Abstract

The variational solution of a boundary problem for the Ginzburg - Landau equations describing the superconducting multilayer nanostructures critical state is presented. The upper critical magnetic field angular dependences are calculated at different temperatures. The effective dimensionality and the superconducting phase nucleation parameters versus magnetic field strength orientation are explained.

Литература

1. JinB.Y., Ketterson J.B. //Adv. Phys. 1989. Vol. 38. P. 189.

2. Лыков A.H. IIУФН. 1992. T. 162, № 10. C. 1-62.

3. Blatter G., Feigelman M.V., Geshkenbein V.B., etal. //Rev. Mod. Phys. 1994. Vol. 66. P. 1125.

4. Cava R.J. II Science. 1990. Vol. 247. P. 656.

5. Ruggiero S.T., BarbeeJr. T.W., BeasleyM.R. //Phys. Rev. B. 1982. Vol. 26. P. 4894.

6. BanerjeeI., Yang Q.S., Falco C.M., Schuller I.K. //Phys. Rev. B. 1983. Vol. 28. P. 5037.

7. Guimpel J., de la CruzM.E., etal. Ill Low Temp. Phys. 1986. Vol. 63. P. 151.

8. Kanoda K, Mazaki H„ Yamada Т., et al. И Phys. Rev. В 1986. Vol. 33. P. 2052.

9. Locquet J.P., Sevenhans W„ Bruynseraede Y, etal. II IEEE Trans. Magn. 1987. Vol. 23, P. 1393.

10. Jin B.Y., Ketterson J.B., McNiffE.J., et al. И J. Low Temp. Phys. 1987. Vol. 69. P. 39.

11. Coccorese C., Attanasio C„ Mercaldo L.V., etal. //Phys. Rev. В 1998. Vol. 57. P. 7922.

12. Takahashi S., TachihiM. //Phys. Rev. В 1986. Vol. 33. P. 4620.

13. BroussardP.R., Geballe Т.Н. //Phys. Rev. В 1988. Vol. 37. P. 60.

14. Dediu V.I., Kabanov V.V., Sidorenko A.A. //Phys. Rev. В 1994. Vol. 49. P. 4027.

15. KoperdraadR.T. W., LodderA. //Phys. Rev. В 1996. Vol. 54. P. 515.

16. Ciuhu C„ Lodder A. //Phys. Rev. B. 2001. Vol. 64. P. 224526.

17. Kushnir V.N., Prischepa S.L., DeliaRoccaM.L., etal. //Phys. Rev. В 2003. Vol. 68. P. 212505.

18. Kushnir V.N., Prischepa S.L., Cirillo C., etal. IIto be published.

19. Кушнир B.H., ПетровА.Ю., Прищепа C.JI. II Физика низких температур. 1999. Т. 25. С. 1265 [Low Temp. Phys. 1999. Vol.25. Р.948].

20. Cirillo С., Attanasio C„ Maritato L„ et al. II J. Low Temp. Phys. 2003. Vol. 130. P.509.

21. TinkhamM. Introduction to Superconductivity. Krieger, Malbar, FL, 1980.

22. Абрикосов A.A. Основы теории металлов. M., 1987.

23. Chun C.S.L., Zheng G.G., VicentJ.L., Schuller I.K. //Phys. Rev. B. 1984. Vol.29 P.4915.

24. Banerjee /., Schuller I.K. II J. Low Temp. Phys. 1984. Vol.54. P.501.

25. TinkhamM. //Phys. Rev. 1963. Vol.129. P. 2413. //Phys. Lett. 1964. Vol.9. P.217.

26. Yamafuji K, Kusayanagi /<'., Irie F. И Phys. Lett. 1966. Vol.21. P.l 1.

27. Saint-JamesD. //Phys. Lett. 1965. Vol.16. P. 218.

28. Saint-JamesD„ De GennesP.G. //Phys. Lett. 1963. Vol. 7. P. 306.

29. TakanakaK. Ill Phys. Soc. Jpn. 1989.Vol.58. P. 668. ibid 1992. Vol.61. P.2865.

30. TakezawaN., Koyama Т., TachikiM. //PhysicaС 1993. Vol. 207. P. 231.

31. Справочник по специальным функциям / Под ред. М. Абрамовича и И. Стиган. М., 1979.