Токофазовое соотношение джозефсоновских SFS-контактов с учетом s-d-рассеяния Текст научной статьи по специальности «Физика»

Токофазовое соотношение джозефсоновских SFS-контактов

с учетом s-d-рассеяния

Т.С. Румянцева1, Н. Г. Пугач2,а, Н. В. Кленов3,6

1 РНЦ «Курчатовский институт». Россия, 123182, Москва, пл. Академика Курчатова, д. 1.

2-3 Московский государственный университет имени М. В. Ломоносова, физический факультет,

2 кафедра магнетизма; 3 кафедра атомной физики, физики плазмы и микроэлектроники.

Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.

E-mail: " pugach@magn.ru, ь nvklenov@nm.ru Статья поступила 01.06.2010, подписана в печать 18.10.2010

Получены аналитические выражения коэффициентов токофазовой зависимости джозефсоновских переходов с ферромагнитной прослойкой в модели, в которой рассеяние из s- в d-зону считается основным типом рассеяния носителей заряда в ферромагнетике. Показано, что величина коэффициентов осциллирует и затухает с увеличением толщины прослойки. Как период, так и длина затухания для первой и второй гармоник отличаются в два раза. Получен вид температурной зависимости критического тока, демонстрирующий 0-7Г переходы по температуре. Показана возможность создания потокового кубита на основе таких контактов.

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

УДК: 538.945. PACS: 74.50,+г, 74.45,+с.

Введение. Джозефсоновские контакты и возможность их применения для создания квантового компьютера

Джозефсоновский контакт представляет собой два сверхпроводника (5), разделенных некоторой слабой связью. В настоящей работе рассматривались контакты, разделенные ферромагнитной прослойкой (Т7) (рис. 1). За счет взаимодействия волновых функций левого и правого сверхпроводящих берегов возникает стационарный эффект Джозефсона: через прослойку протекает сверхпроводящий ток . Вид зависимости тока /5 от разности фаз параметра сверхпроводящего порядка на берегах перехода можно получить из следующих соображений: во-первых, в отсутствие разности фаз ток должен равняться нулю, во-вторых, смена знака разности фаз должна приводить к смене направления тока, и, в-третьих, изменение разности фаз на 27т не должно влиять на ток [1]. Наиболее простым и часто встречающимся соотношением, удовлетворяющим этим условиям, является синусоидальное соотношение

=к эш

где ¡с — критический ток перехода, максимально возможный сверхпроводящий ток через прослойку. Возможны и другие случаи, в частности в разложении могут присутствовать гармоники высшего порядка:

hW) = Л sin ip + В sin 2ip + С sin 3ip ■

(1)

Js

Рис. 1. Джозефсоновский контакт с ферромагнитной прослойкой: 5 — сверхпроводник, Р — ферромагнетик, — сверхпроводящий ток

В контактах с ферромагнитной прослойкой величина критического тока может быть отрицательна /с < 0. Это состояние называется 7г-состоянием, так как разность фаз в отсутствие тока, т.е. в основном состоянии, становится равна не нулю, а 7г [2]. Причина возникновения этого состояния заключается в воздействии магнитного поля на куперовскую пару. Так как ее составляют электроны с противоположным направлением спина, то вследствие зеемановского расщепления обменное поле дает разный вклад в энергии электронов пары. Это приводит к осцилляциям параметра порядка. Устанавливается сверхпроводящее состояние типа ЛОФФ фазы (Ларкина, Овчинникова, Фулде, Ферелла) [3, 4]. Тогда параметр порядка приближенно выглядит следующим образом:

Ф/ = Ф/о ехр = ф/о ехР ехр

1 1 1

~г = 7— +1-^,

Ч Ч1 Ч 2

где Фд) — параметр порядка на границе сверхпроводник-ферромагнетик, — комплексная длина когерентности, — длина затухания, 2тг^2 — период осцилляций параметра порядка. Для несинусоидальной токофазовой зависимости (1), при изменении разности фаз на 7г, меняется знак нечетных гармоник. Если основное состояние достигается при разности фаз, равной нулю, тогда это 0-состояние и первый коэффициент в разложении (1) А> 0. Если основное состояние достигается при разности фаз, равной тг, тогда это 7г-состояние и А <0. Важное значение имеют 7г-контакты в связи с возможностью их использования для создания квантового бита.

Для создания квантового бита необходимо построить систему с двухъямным потенциалом [5]. Между ямами потенциала возможно туннелирование, снима-

ющее вырождение, тогда двум базисным состояниям кубита отвечают антисимметричная и симметричная суперпозиции состояний в левом и правом минимумах [6]. Для реализации таких кубитов необходимы контакты с нетривиальной токофазовой зависимостью [7, 8]. В качестве схемы для реализации так называемого «потокового» кубита рассмотрим двухконтактный низкоиндуктивный сверхпроводящий интерферометр. Он представляет собой два включенных параллельно джо-зефсоновских контакта. Такой кубит называется потоковым, так как величиной, кодирующей информацию, является «захваченный» магнитный поток. Мы рассматривали случай с одинаковыми параметрами контактов и с одинаковыми токофазовыми характеристиками. Условия для реализации кубита выглядят следующим образом [7-9]:

\2(В1+В2)\>\А1+А2\, (2)

Из этих условий следует, что необходимо найти контакты, токофазовая зависимость которых достаточно сильно отличается от синусоидальной.

1. Модель

1.1. Плотность состояний s - и с?-электронов

В работе [10] группой Рязанова были представлены экспериментальные данные, показывающие, что в ферромагнетике превышает длину затухания параметра порядка приблизительно в три раза, что вступало в противоречие с решением, следующим из уравнения Узаделя [2]. Так как это требовало объяснения, было предложено учесть следующие механизмы рассеяния в ферромагнетике: профессор А. И. Буздин [11] предложил объяснить это рассеянием с переворотом спина, и группа A.B. Ведяева [12] предложила модель, учитывающую s — d-рассеяние в ферромагнетике. В настоящей работе использована модель для расчета токофазовой зависимости, разработанная группой Ведяева.

В работе [13] приведен модельный расчет плотности состояний для никеля и для сплава меди с никелем (рис. 3 и 4 в работе [13]). Авторы работы показали, что на уровне Ферми плотность ¿-электронов значительно больше, и, так как вероятность рассеяния пропорциональна плотности состояний, вероятность рассеяния в ¿-зону должна быть больше. Мы рассматривали контакты со сверхпроводниками из алюминия, в нем нет электронов в ¿-зоне, следовательно, ¿-электроны отражаются зеркально на границе сверхпроводника и ферромагнетика и не образуют куперовских пар, поэтому рассеяние в ¿-зону приводит к развалу пар.

1.2. Выражение для тока в модели, учитывающей s-d -рассеяние

С учетом описанного типа рассеяния группой Ведяева были решены уравнения Горькова и получено следующее выражение для сверхпроводящего тока [12]:

7Г ^

k dk Re

(1 - R2) sin p

(3)

где T

температура,

Den(</?,a,...)'

квазиимпульс Ферми

с, с\, С2 — действительные части импульсов квазичастиц р, pf, pi в сверхпроводнике и ферромагнетике соответственно, р — разность фаз параметра порядка, а — толщина ферромагнитной прослойки, ш = 7гГ(2я + 1) — мацубаровские частоты, п — целое число, Den — осциллирующая функция, зависящая от параметров перехода. Так как фаза р является аргументом функции Den, токофазовая зависимость отлична от синусоидальной.

1.3. Выражения для гармоник токофазовой зависимости

Рассмотрим токофазовую зависимость, описываемую выражением (3). Введем следующие обозначения: вектор-строка Ai = (cos(apf) sin(apf)), вектор-столбец

( cos(ap4_)\

А2 = , \ , где pt, pi — импульсы электронов \ып(ар4.) J

в состояниях со спином против и вдоль намагниченности в ферромагнетике, которые являются решением уравнений Горькова [12]. Введем матрицу Q:

Q =

'^■(l-R2)+^¡P2(l-r2) + rR)

где Д -воднике,

г

Д2 1 V1 ' Д2 V1 1 1V ) 1 Д2

ширина энергетической щели в сверхпро-

= Cl~C2

Ci+C2

коэффициент отражения: г =

р _ с,с+с2с ^ 14СПОЛЬЗуЯ введенные обозначения, выраже-

С +C¡ С>2

ние для тока (3) можно записать в более удобной форме

(1 - R2) sin р

TT ^

Re

(1 — R2) cosíp-{-\\Q\2

dE.

Теперь можно найти коэффициенты разложения этого выражения в ряд Фурье по р. Используя формулы для производящей функции полиномов Чебышева второго рода Un(co$p) [14], получим

z\, 2 =

(AIQA2)

(1 ^R2)

(A1QA2)Y (1 -&))

1.

Получены следующие выражения для коэффициентов А и В:

z dE,

¿:J dE.

(4)

R — коэффициент андреевского отражения: R =

— С--С2С,

где в качестве г выбирается то из значений г\ и г2, которое удовлетворяет условию < 1.

2. Исследование токофазовой зависимости 5^5-контактов

2.1. Поиск областей, в которых происходит переход из 0- в 7г-состояние

С помощью полученных выражений (4) была исследована токофазовая зависимость для джозефсоновских контактов с прослойкой из ферромагнитного металла. Модельные параметры сверхпроводника соответствуют алюминию: критическая температура Тс = 1.196 К, квазиимпульс Ферми Щ = 1.75 А^1 [15]. Параметры ферромагнитной прослойки приблизительно соответствуют значениям для никеля: квазиимпульс Ферми

15 ВМУ. Физика. Астрономия. № 1

и

В

lis

IS

3

2.0

1.5

5 Э

Я» 10 § £

а U

0.5

1 1 1 1 1 б

1 1 1 1 1 1 а-28А

1 \ 1 \ 1 \ 1 \ 1 \ 1 ------а ~ 27,8 А

2 3 4 5 Фаза ф, рад

0.2 0.4 0.6 0.8 1.0 Температура, К

1.2

Рис. 2. а — Токофазовая зависимость. Пунктиром показано 7г-состояние, толщина прослойки а = 32 А; сплошной линией показано 0-состояние, а = '23 А; штриховой линией показан переход из 0- в 7г-состояние, а = '27.8 А. 7=0.5 К, /=28 А. б — Зависимость сверхпроводящего тока от температуры. Сплошной линией показано 0-состояние (толщина прослойки а = '28 А). Штриховой линией показан переход из 0в 7г-состояние (толщина прослойки а = '27.8 А). Температура перехода из 0- в 7г-состояние 7 = 0.12 К,

I ='28 К

^• = 1.18 А -1, И = 0.2 эВ [16]. На рис. 2, а представлена зависимость сверхпроводящего тока от фазы при различных значениях толщины прослойки. Сплошной линией показана зависимость вдали 0-7г-перехода, это обычное 0-состояние с синусоидальной зависимостью. Пунктирной линией также показана зависимость вдали от 0-7г-перехода, но при другом значении толщины прослойки. По отношению к предыдущей кривой она имеет сдвиг по фазе на тг, направление тока противоположно направлению тока в 0-состоянии — это синусоидальное 7г-состояние. Промежуточное состояние перехода из 0в 7г-состояние показано штриховой линией. Ее вид отличен от синусоидального, в токофазовой зависимости должна присутствовать, по крайней мере, вторая гармоника. Переход из 0- в 7г-состояние можно также наблюдать при изучении зависимости критического тока от температуры. Такой переход наблюдался в эксперименте группой Рязанова [17]. На рис. 2,6 представлены модельные кривые температурного перехода. Если в 0-и 7Г-СОСТОЯНИЯХ абсолютная величина критического тока просто убывает с температурой, то в промежуточном состоянии есть температура, при которой критический ток меняет свое направление и контакт переходит из 0в 7г-состояние, в этой точке токофазовая зависимость должна быть несинусоидальна. На приведенном графике это температура 0.12 К.

2.2. Зависимость критического тока и гармоник от толщины прослойки

Для более внимательного исследования гармоник высшего порядка были рассчитаны кривые, показывающие зависимость коэффициентов от толщины прослойки. Результаты расчетов показаны на рис. 3. Минимумы на кривых означают переход через нуль и смену знака. Области, где первая гармоника положительна, это 0-состояние, области, где она отрицательна, это 7г-состояние. Графики для тока и первого коэффициента совпадают практически при всех значениях толщины прослойки, и, следовательно, ток определяется в основном первой гармоникой, т.е. синусоидален: = Аэт^. Области, где первая гармоника меняет свой знак, соответствуют промежуточному состоянию

1(Г

В

10 10 10

- ° 10 § и

а £ ю-5

о.

5Г я О

10" 10" 10" 10"

V\,\A\

-\в\

0 50 100 150 200 Толщина прослойки, А

Рис. 3. Зависимость критического тока и гармоники А (сплошная линия) и гармоники В (штрихпунктирная линия) токофазового разложения от толщины прослойки. Область первого минимума первой гармоники (25-31 А) показана отдельно на врезке. Штрихпунк-тирной линией на врезке показана удвоенная гармоника В. Область выполнения условий: 27.7-30.5 А. 7 = 0.5 К, / = 28 А

перехода из 0- в 7г-состояние. В этих областях коэффициент А проходит через нуль, в то время как коэффициент В еще отличен от нуля, поэтому в таких областях и реализуется несинусоидальное состояние, показанное на предыдущих рисунках /¿-(р) = Вз1п'2р. Периоды осцилляций первого и второго коэффициентов отличаются приблизительно в два раза. Благодаря этому появляются области, где первый коэффициент А находится в своем минимуме, в то время как второй коэффициент В близок к максимуму и токофазовая зависимость несинусоидальна. С увеличением толщины прослойки коэффициенты не только осциллируют, но и затухают, причем второй коэффициент В затухает в два раза быстрее первого коэффициента F, длина затухания которого в данном случае приблизительно равна заданной длине свободного пробега. Значение тока в областях с несинусоидальной зависимостью тем меньше, чем больше толщина прослойки, но даже при

небольших толщинах значения тока в этих областях малы по сравнению со значением тока в 0- и 7г-со-стояниях, поэтому обнаружить несинусоидальную зависимость на практике довольно трудно и имеет смысл рассматривать области с небольшой толщиной прослойки, где критический ток еще сравнительно велик, а различия в гармониках еще не так существенны. В работах [18, 19] были представлены экспериментальные свидетельства существования второй гармоники, что позволяет нам и дальше производить расчеты для поиска таких контактов.

2.3. Влияние длины свободного пробега на затухание критического тока

В используемой модели коэффициент затухания параметра порядка пропорционален длине свободного пробега [12], поэтому имеет смысл рассматривать контакты с большей длиной свободного пробега. На рис. 3 на врезке изображены результаты более подробного расчета области, показанной на основном графике овалом. Предположительно, в этой области токофазовая зависимость не синусоидальна, и эту область можно рассмотреть в качестве возможной для создания ку-бита. Исследование условий (2) в этой области показывает, что они выполняются при 27.7-30.5 А. Также были проведены расчеты для контактов с меньшей длиной свободного пробега I. Результаты представлены на рис. 4. Видно, что при меньшем значении I ток затухает быстрее, так как чем меньше I, тем больше эффект от рассеяния. Для меньшей длины свободного пробега область выполнения условий очень узкая. Поэтому более чистые контакты являются более перспективными для создания кубита. Тем не менее такую узкую область также можно реализовать, если изготовить контакт с прослойкой, близкой по толщине к требуемому значению, затем, меняя температуру, «ловить» искомую область.

Для исследования влияния в-с!-рассеяния на протекание тока был произведен расчет зависимости критического тока от толщины прослойки при очень большой

101

1(Г 10"1 10":

10"

« N

« ч

15 ю-4

S8 _5 & * 10 5

о.

5Г я U

10~° ю-7 10~8 10"9

чу/ ^ч. \ :

У \v' \ s \ V Г \ ! \ \ У

\ -----/= 12.2 Â < Ч

- / = 28Â \

........ 1= 1000Â \ " ^ \

О 50 100 150 200

Толщина прослойки, А

Рис. 4. Зависимость критического тока от толщины прослойки при разных значениях длины свободного пробега /: сплошной линией при / = 28 А, штрих-пунктиром — при / = 12.2 А, пунктирной линией — / = 1000 А. Область выполнения условий при / = 12.2 А: 30.55-31.1 А. Г =0.5 К

длине свободного пробега. Результаты представлены на рис. 4. Для «чистого» контакта затухание сверхпроводящего тока происходит медленнее: не по экспоненциальному, а по степенному закону [19].

Заключение

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

Работа выполнена при финансовой поддержке РФФИ (гранты 09-02-12176, 10-02-00569, 10-02-00870).

Список литературы

1. Шмидт В.В. Введение в теорию сверхпроводимости. М., 1982.

2. Buzdin A.I. И Rev. Mod. Phys. 2005. 77. P. 935.

3. Larkin A.I.. Ovchinnikov Yu.N. //JETP. 1964. 47. P. 1136.

4. Fulde PFerrell R.A. 11 Phys. Rev. A. 1964. 135. P. 550.

5. van der Wal Caspar H., ter Haar A.C.J., Wilhelm F.K. et al. // Science. 2000. 290. P. 773.

6. Валиев К.A. I ! УФН. 2005. 175, № 1. С. 3.

7. Klenov N.. Kornev V., Vedyayev A. et al. 11 J. Phys. Conference Series. 2008. 97. P. 0"l2037.

8. Кленов H.В., Корпев В.К., Соловьев И.И. и др. // Нелинейный мир. 2005. 3, № 1-2. С. 75.

9. Amin M.H.S., Smirnov A.Yu., Zagoskin A.M. et al. // Phys. Rev. В. 2005. 73. P. 064516-1-5.

10. Ryazanov V.V., Oboznov V.A., Prokofiev Л.S. et al. // J. Low Temp. Phys. 2004. 136. P. 385.

11. Obobznov V.A., Bolginov V.V., Feofanov A.K. et al. // Phys. Rev. Lett. 2006. 96. P. 197003.

12. Vedyayev A., Ryzhanova N.. Pugach N. // J. Magnetism ancf Magnetic Mat. 2006. 305. P. 53.

13. Breuers F., Vedyayev A.V., Giorgino M. // Phys. Rev. В. 1973. 7. P. 380.

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

15. Ашкрофт Н., Мермип Н. Физика твердого тела. Т. 1, 2. М„ 1979.

16. Shelukhin V., Tsukernik A., Karpovski M. et al. // Phys. Rev. В. 2006. 73. P. 174506.

17. Frolov S.Al, Harlingen DJ. van, Oboznov V.A. et al. // Phys. Rev. В. 2004. 70. P. 144505.

18. Sellier H., Baraduc C., Lefloch F., Calemczuk R. // Phys. Rev. Lett. 2004. 92. P. 257005.

19. Рязанов В.В., Обозное В.А, Больгинов В.В., Россолен-ко А.Н. // Матер, симпоз. «Нанофизика и наноэлектро-ника». Нижн. Новгород, 2008. 1. С. 42.

20. Buzdin А.1., Bulaevski L.N., Panyukov S.V. 11 JETP Lett. 1982. 35. P. 147.

16 ВМУ. Физика. Астрономия. M 1

Current-phase relation of SFS junctions in the presence of s-d scattering T.S. Rumyantseva1, N. G. Pugach1' , N.V. Klenov , f

1 Russian Research Center "Kurchatov Institute". Kurchatov sq. 1, Moscow 123182, Russia. 2Department of Magnetism; 3Department of Atomic Physics, Plasma Physics, and Microelectronics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia. E-mail: apugach@magn.ru, b nvklenov@nm.ru.

An analytic form of the CPR coefficients in superconductor-ferrornagnet-superconductor junctions is presented taking into account s-d scattering as a basic pair-breaking mechanism. It is shown that magnitude of the coefficients oscillates and damps as a function of the ferromagnetic layer thickness, the oscillation period and damping length are differ in two times. Critical current has been calculated as a function of temperature. It is demonstrated the possibility of a design flux qubit based on such kind of Josephson junction.

Keywords: superconductivity, Josephson junctions, superconductor-ferromagnet structure, flux qubit. PACS: 74.50,+r, 74.45,+c. Received 1 June 2010.

English version: Moscow University Physics Bulletin 1(2011).

Сведения об авторах

1. Румянцева Татьяна Сергеевна — инженер; тел.: (906) 051-84-94, e-mail: rumyanceva-tatyana@yandex.ru.

2. Пугач Наталия Григорьевна — канд. физ.-мат. наук, ассистент; тел.: (495) 939-47-87, e-mail: pugach@magn.ru.

3. Кленов Николай Викторович — канд. физ.-мат. наук, ассистент; тел.: (495) 939-43-51, e-mail: nvklenov@nm.ru.