ИЗВЕСТИЯ АКАДЕМИИ НАУК СССР
Т. XXXVI, № 7 |
СЕРИЯ ФИЗИЧЕСКАЯ
Р. Г. ХЛЕБОПРОС и Л. В. МИХАЙЛОВСКАЯ
СПИНОВЫЕ КОЛЕБАНИЯ В ФЕРРОМАГНИТНОМ СЛОЕ
Возросшие требования к интерпретации очектра спин-волнового резонанса тонких магнитных пленок из-за увеличения диапазона частот и чувствительности эксперимента обусловили появление в последнее время значительного количества теоретических работ, посвященных исследованию энергетического спектра ферродиэлектрического слоя. При этом обнаружен целый ряд особенностей спектра спиновых колебаний; среди них наиболее интересны поверхностные магнито статические волны [1—3] и поверхностные волны, обусловленные закреплением?линов на границе тонкого слоя [4—7].
В настоящей работе рассмотрим в феноменологической постановка •спектр магнитных колебаний изотропного однородно-намагниченного под произвольным углом ферромагнитного слоя*. Равновесное положение намагниченности М0 (рис. 1) определяется уравнением магнитостатики:
coo sin (0о — 9) = 0,5&>м sin 260, - (1)
где (во = у#о, 1(°дг = 4яуЛ/0, постоянное магнитное поле Н0 направлено под углом 9 к нормали (оси z). Малые колебания намагниченности относительно равновесного положения удобно представить в компонентах ту и тв = тож/соз6о = — mz / sin 60. В представлении плоских волн m, h, e ~ ~ exp [i(kr — a>t) ] из уравнений Ландау — Лифшица и Максвелла можно получить в пренебрежении эффектами релаксации соотношение между.со и k в виде
Здесь |
——, G = — - cos (00 — Э) — cos2
обменный параметр
т] = А I 2пМ02, (А — постоянная обменного взаимодействия). &2 = к1 -f- &Д х = у&я2 -f- ky2, ki2 = kf,2 -j- /Cy2, ke = kK cos 00 — kz sin 60, величина волново-
го вектора электромагнитной волны fte = -^-Ve, e — диэлектрическая по^
^стоянная, с — скорость света.
Это дисперсионное соотношение представляет при заданных частоте Q, поле ^о, угле 6 и неоднородности в плоскости пленки % уравнение шестой -степени относительно kz. Колебания магнитного момента при этом являются результатом сложения шести плоских волн:
] \ А„ exp(ikznz).
* На Всесоюзной конференции по магнетизму (Красноярс-к, 1971 г.) нам стало известно, что подобная задача независимо решалась Б. Н. Филипповым. Статья с близкими к пашей работе результатами опубликована в журнале «Физика металлов и металловедение». 32, 911 (1971).
колебания в ферромагнитном слое
Учет граничных условий для электромагнитного доля и для магнитного момента позволяет определить допустимые собственные значения /ё„, ш и с точностью до постоянного множителя соответствующие им амплитуды А,,.
Если выполняется равенство k,sin2fl« = 0, то уравнение (2) переходит в бикубическое, которое мы и будем рассматривать в дальнейшем*.
На рис. 2 изображена симметричная по k, дисперсионная кривая (2) при е = 0 и х = kt/. Собственные колебания при этом представляют сумму трех стоячих волн, две из kotopbix правополяризованны, а одна — ле-вополяризованна п является всегда поверхностной, так как имеет k,s < 0.
Рис. L Равновесное положение намагниченности М.
Рис. 2. Дисперсионное соотношение (2) при 6 =
хг = 0, 2— V'- = 2.°10-"!,? — Г]х2 = 2.10-', 4-их2 = 2). Координаты пересечения с осью Q = 0: ' — в, 2r|*22., ==^(C + sin во) ±
_ ± У(С + I)2 -f- 4т)^ж2. Здесь и в дальнейшем индексы «+» и «—» означают соответственно право- и левополяризованные колебания |
r,kl-2,0 |
•0,2 yjrf!J/!JI/ 0,1 /jit2 Рис. 2 |
7*1 - |
я/2,
Обращает на себя внимание наличие полосы частот, где правогюдяризован-ные колебания существуют с комплексным значением kz. На эту особенность спектра в частном случае перпендикулярно намагниченной пленки со свободными спинами на границах впервые обратили внимание авторы [9].
Влияние условия, в которых находятся слшиы на поверхности слоя, учтем, пренебрегая электромагнитным запаздыванием в двух предельных случаях: \} а = 0, х = 0, вэ =t= 0; 2) е -= О, и =И= О, О» = О.
1. В случае однородных в плоскости колебаний (х = 0) решения диг,-
* Асимметрия дисперсионной кривой, обусловленная наличием в (2) членов с ае-четными степенями fen существенна в магнитостатической области, проявляется
дри о <С 8й *С —— [1 исследовалась, в [3\. •...
Р. Г. Хлебопрос и Л. В. Михайловская
персионного уравнения (2) записываются так:
(3) |
ГДе |
Воспользуемся следующими обменными граничными условиями для намагниченной под произвольным углом пленки, полученными в [10]:
dm* Bi
(4) |
dmy р2 —— -_-my = 0 dzn d
при 2 = ±d. Здесь б! = p cos 28о, р2 = р cos2 Go, 3 = К.4/А, Ks — постоянная поверхностной анизотропии, 2d — толщина слоя. В этой записи, если Кг > 0, то легкая ось поверхностной анизотропии лежит в плоскости слоя. При другом знаке — ортогонально плоскости.
Рис. 3. Графическое решение -уравнений (3) и (5) в областях: а — G >
> ti, б — G < ti. Кривые А соответствуют уравнению (5). кривые В —
уравнению (3), индексы 0. 1, 2, 3. со соответствуют частотам Q -> О, qi =
" ' " - -• ~. ~ -. - -^ анизотропия
Нетрудно убедиться, что задача содержит только два типа решений — симметричные и антисимметричные относительно замены z на —z. Подставляя симметричные решения
тв = exp (~i(ut) (Aicosktiz -j-^acos kziz], my = exp (— mt) (AiaicoskziZ-t-AzazCOBkrtZ)
в (4), получим следующее уравнение для нахождения собственных значений kti и k2z *:
«1 (Pi cos «i + mi sin mi) (pacos iiz + Ma sin K2) —
(5)
— a2 (pz cos в, + b! sin в,) (p, cos B2 + и* sin в,) = О. Здесь
Здесь и в дальнейшем длн антисимметричных решений соответствующие выражения получаются заменой sin и на —cos и и cos и на sin и.
ебания в ферромпзпитном слое
Форма собственных колебаний магнитного момента определяется отношением амплитуд:
,, |
Ui Sill «i + pi COS И-i |
при
Амплитуда магнитостатической ветви с волновым вектором этом равна нулю.
В предельном случае р2 — Pi-^О правополяризованные колебания с волновым вектором kzi и левололяризованные, которым в дисперсионном уравнении соответствует решение kzZ2 < 0, становятся независимыми.
Из исследования выражения (5) видно, что граничные условия запре-
щают решения в области А*Д < 0, k^ < 0, если р\ < 0, р2 < 0. Выражения (3) для данной частоты Q можно рассматривать как параметриче-
ское задание дисперсионной кривой на плоскости (&г1, А;га). В качестве параметра можно взять эффективное поле G, В области k* < 0 делается замена: &г = г А;/. Тогда собственные значения kzi и /с;2 находятся пак пересечения дисперсионной кривой (3) с кривой (5) на плоскости (йя, &*/), если G <С 7\ и на плоскости {^г!7, Агг/), если G > tl
На рис. 3 показано графическое решение уравнений (3) и (5) для угла 90 <.я/4. Видно, что в зависимости от соотношения между Q, 60 л {} существует одно или два пересечения дисперсионной кривой (3) с кривой (5), т. е. существует одно или два собственных чисто гиперболических колебаний в системе. Собственные значения и/, и2', лежащие в области /, удовлетворяют соотношениям:
pi < р!вф< и/Шв/ < р2, г Tt + Tz
Рис. 4. Дисперсионные кривые (10) для различных значений р: 1 ~ р = 0; 2~ о = = 4-10-2;,? —р = 4-10-'; ' 4—р=1; 5 — р = 2 |
Из выражения (6) легко убедиться, что при этом \Ai\ >>!>-|Л2|, т. е. ведущей в колебании будет правополяризованная волна. Одновременно с этим собственным колебанием для частот Q < Q1[p (QKP определяется из уравнения |/7\ + Та. 11 ~= ву, где в20ЧЬ ц20'= р,ра / р, вф) существует собственное колебание (В1,в3), лежащее в области // и удовлетворяющее соотношениям:
•где рг=ф(У, + Тг) = r^z + J'jp,. В этом случае при достаточно низких частотах или при pz-p,^0 главной в собственном колебании становится левополяризованная волна, т. е. |42|>|Л,|.
Спектр собственных колебаний с &„2 > 0 может проявиться ник видно из рис. d, 6, начиная с пекоторой частоты, когда радиус дисперсионной окружности Т/(Т, + T2)/t> в2„'. Собственными значениями являются не-
Р. Г. Хлебопрос и Л. В. Михайловская
ресечения (3) с (5), попадающие в сектор, ограниченный углом ф = arctg ут\ / Tz. Колебания намагниченности вокруг положения равновесия.представляют теперь суперпозицию тригонометрических правополя-ризовапных волп с поверхностными левополяризоваппыми. В области малых Q, когда У(/1 + tz} I t^Uzo', могут существовать решения, в которых ведущей является гиперболическая волна (|Л2| > |.4i|). Число собственных колебаний по мере увеличения частоты возрастает.
Подобное рассмотрение может быть проведено для анизотропии другого знака и распространено на углы я / 4 ^ 6 и ^ л; / 2. Аналогичная, но более сложная для оценок ситуация имеет место в спектре антисимметричных колебаний.
2. Рассмотрим теперь влияние состояния поверхностных спинов на спектр неоднородных в плоскости (к =^ 0) колебаний намагниченности в частном случае 60 = 0 и (3i = 02 = р.
На рис. 4 изображено дисперсионное соотношение G = G(k?} при заданных Ц и к. Характер дисперсионных кривых определяется величиной параметра р ~ цк2 / Q2. Если р < 1, то в зависимости ог значения G имеем три вида решений для &n2 — А„2(С, Й,хв), (га = 1,2,3). Если G^Gi, то &Д kzz > 0, ksz <С 0. В области gi < G <i Gz решения для &Д k22 — комплексно-сопряженные, &з2 < 0. При G ^ G2 вес решения ftn2 <C 0. Переходные значения полей GtJ& определяются выражениями
(7) |
(8) |
Соответстиующие им значения квадратов волновых векторов равны
/ 1 п/~1 /1~~\
1 + — p=F I/ —р(1 +.—р)
С увеличением р все болео расширяется область значений G, в которой ki2 и k22 — комплексно-сопряженные. Если р ^ 1, то она простирается от О до G2, определяемого формулой (7). В области G > G2 результирующее колебание, как и прежде, представляет суперпозицию трех поверхностных (гиперболических) волн. Ситуация р ~ 1 может возникнуть в радиочастотном диапазоне, если неоднородности в плоскости пленки определяются ее магнитной подструктурой.
Интересующее нас распределение памарштчезкости по толщине слоя определяется проекцией волнового век-тора* Jc на ось z: kin = ^k,,z — -/с (п — 1, 2, 3). Видно, что по мере увеличения к уменьшается область тригонометрических решений (fez2 > 0) 'и, еоли выполняется неравенство т]кг(1 4~ "П^2) > ^2, то в спектре полностью отсутствуют тригонометрические решения.
Как и прежде, решения задачи, разбиваются на "симметричные и антисимметричные..Подставляя в обиенпа'ей'Траничноез. условие (4) решения симметричного класса
-А.п cos kzn |
= exp[i(k,x + kyy-№
/ —(oi)lV ^^„cos
г ферро. |
и выполняя требования непрерывности тангенциальной составляющей
4nk(m -k)
магнитного поля U = — ———-———— и пормальноп составляющей маг-
к
нитяой индукции 1) = h -f- 4лш на границах слоя г = ±й, получим уравнение
которое совместно с дисперсионным соотношением
* = (G |
(10) |
определяет спектр собственных колебаний намагниченности. В (9) введены следующие обозначения:
v cos и — и cos и |
(для простоты взято А'* = х; если н-„2 <С 0, то делается замена: ип = т/). Отношения между амплитудами стоячих волн для найденных собственных колебаний можно записать как
(И) |
Проанализируем предельные случаи. Пусть £^0. Тогда из (10) имеем: G
k\z ->- -}~og/ hz" = к2 —^rz—~^ i k3z~^- — oo, а уравнение (9) переходит в
l£ — Li'
уравнение Vz = 0. Таким образом, в нулевом по 'С, приближении (£—>-0) уравнения (9) и (10) определяют спектр уокеровских магнитостатических колебаний ферр о диэлектрического слоя.
При у. —>• 0 решениями уравнения (10) являются: " '
Q-G |
Q + G
= 0,
Ц Т]
а (9) принимает вид: UiU$ = 0. Амплитуда Аг = 0, ветви право- и лево-поляризованных волн оказываются независимыми.
Рассмотрим теперь случай р<С1. Тогда решения уравнения (10) можно записать как:
(т|хг)(Й-С) 1
A.»^h-
2t] L 2
1 (Й-G)2
(t)x!)2G(5Q + G) |
1 (Q-G)'r _tpr
~>i rf L (q+
В этом приближении (7) и (8) упрощаются:
Gi,2 ^ £3 1=РУ2р + -—р
/'. Г. Хлебопро |
Л. В. Михайловская |
Подставляя эти выражения в (9) и (11), можно убедиться, что для собственных решений, попадающих в область полей G<^Gi или G^>GE, ведущей является одна из волн ai, Az или А3. Так, например, в нулевом приближении (р —>- 0) собственные значения обменного киттелевского спектра и.1 определялись из уравнения
«1 lgui + Р = 0 (12)
н в двух предельных случаях |3 —»- 0 н | ния
оо имели соответственно реше-
В данной ситуации при р*^1 (v^l), как и следовало ожидать, мы имеем в этих предельных случаях малые изменения для собственных значений Ад. Для р -*• 0 вместо уравнения (12) получаем уравнение:
«I tg ifi — - |
Если ||3| л# у-^ 1, то имеем '"_ ' -••"-:.1 - • ". •
[ |
у(6 + v) I.,-
1 ———^ • • —• I = 0.
Во всех этих трех случаях, как видно из (И), имеем: \АГ\ ~> \AZ\ > М8|. Для гиперболических колебаний, которые в нулевом приближении существовали как независимые kt\ и kz3 и области G > Q и одно Ai3 в области G -< Q н определялись из уравнения Mlp3th u\,3 — р, теперь имеем смешанные колебания, которые удовлетворяют уравнениям
^-^й
В результирующем колебании ведущую роль играет амплитуда поверхностной волны Ai или А3.
В тонком слое при одновременном учете обменных и диполь-диполь-ных членов в уравнении движения проявляется при некоторых условиях модифицированный магнитостатичеекий спектр. Он теперь определяется уравнением
1 + 2^- |
щ tg tit \
i + uagai/'--*-.1.'----».
V / 2£. Отношение алгплитуд
^ii л, I |
Л2 COSWi P + WitgB,
COS И-2 |
P+D
Отсюда видно, что когда |
•^ 1, то ведущей в колебании является
ые колебания в ферромагнитном слое
магнитостатическая ветвь \\ д | J' При обратном неравенстве у<^ р,
о <s; 1 ведущими в колебании всегда является обменная волна (j^^/^ij <C < 1).
Существует область полей Gi^G^ G2, в основном включающая в себя область комплексных решений для и± и иг. в которой собственные значения \и\\ £z \u2\ и поптому соответствующие им волны представлены в ре-
зультирующем колебании равноправными амплитудами \А\\ ~ \A-t\ (\А3\ <^ mij)- Уравнение (9) в области комплексных значений Ui и иг принимает вид
= А
=-asin2a-&sh26, X"(cos 2a + ch 26) =
В приближении р ^ 1 эти колебания являются промежуточными по «длине волны» между магнитостатическими и обменными.
3. Как уже, отмечалось, влияние эффектов электромагнитного запаздывания может оказаться существенным только в магнитостатической длинноволновой области спектра. В пренебрежении обменными полями для ортогонально намагниченного слоя дисперсионное соотношение приобретает вид:
На рис. 5 представлено графическое решение этого уравнения нри заданной частоте Q и различных к2. Видно, что если % < &„, то при любом G решения представляют смешанные колебания, состоящие из право- и ле-вополяризованных стоячих волн. Они становятся независимыми, когда и = 0. Если и > ke, то решения представляют сумму двух стоячих волн
Р, Г. Хлебопрос и Л. В. Михайлоеск.
правой поляризации; объемной магнитостатической и поверхностной электромагнитной, Кроме того, в спектре появляется область значений G < <С gi, в которой решения для k\ 2, а значит и для kt iiZ являются комплексно-сопряженными. Поле gi определяется выражением
2 = 4Й2«- |
Соответствующие нения |
находятся из урав-
Собственные значения k%i и k:2 должны удовлетворять уравнению, полученному при наложении обычных граничных условий на магнитные и электрические поля. В предельном случае &i -^ е (е, 8i — диэлектрические-постоянные внутри и снаружи слоя) это уравнение имеет вид (для симметричных по z решений):
где |
<l>i [ (к2 — ke2) cos и, — х&л sin it* ] cos u2 — Фг[ (^2 — kc&) cos u2 — xkzs sin uz} cos Ui = 0,
Отношение амплитуд равно
А-. cos hi G (kz2 — k,,2) — kfz
Наибольшее влияние электромагнитное запаздывание оказывает на магнигосгатический спектр, когда неоднородность в плоскости пленки х того же порядка, что и величина волнолого вектора электромагнитной волны ke.
Институт физики
им. Л. В. Киренского "'
Сибирского отделения
Академии наук СССР
Литература
Damon R. W., Fshbach J. R, J. Phys. Cbem. Solids, 19, 308 (1961). Б у л а е в с к и и Л. Н., Физ. твердого тела, 12, 799 (1970).
Михайловская Л. В., Хлебопрос Р. Г., Физ. твердого тела, 13, 2786 (1971).
Puszkarski Н., Acta phys. polon., A38. 899 (1970). С о к о л о в В. М., Т а в г е р Б. А., Физ. твердого тела, 10, 1791 (1968). Хлебопрос Р. Г., Михайловская Л. В.. Физ. твердого тела 12 2477 (1970) Филиппов Б. Н., Физ. твердого тела, 9, 1048 (1967).
Михайловская Л. В., X л е б о л р о с Р. Г., Физ. твердого тела 11, 2854 (1969) De Wames R. Е., Wolfram Т., J. Appl. Phys. Letters, 15, 269 (1969) Сух у Р., Магнитные тонкие пленки, «Мир», М., 1967, стр. 310