Рівняння Ландау — Ліфшиця (магнетизм)

Рівня́ння Ланда́у — Лі́фшиця — рівняння, що описує рух намагніченості в наближенні континуальної моделі у твердих тілах. Вперше введене Л. Д. Ландау та Є. М. Ліфшицем у 1935 році.

Формулювання

Для бездисипативного середовища та за відсутності спін-поляризованого струму рівняння Ландау-Ліфшиця зазвичай записується у вигляді

${\frac {\partial \mathbf {M} }{\partial t}}=-|\gamma |[\mathbf {M} \times \mathbf {H} ^{\mathrm {eff} }],\qquad (1)$

где $\mathbf {M} \equiv \mathbf {M} (\mathbf {r} ,t)$ — щільність магнітного моменту (намагніченість), $\gamma$ — деяка феноменологічна стала, $\mathbf {H} ^{\mathrm {eff} }\equiv \mathbf {H} ^{\mathrm {eff} }(\mathbf {r} ,t)$ — так зване ефективне магнітне поле.

Рівняння в основному використовується для феро- та феримагнетиків. У загальному випадку стала $\gamma$ не дорівнює гіромагнітному співвідношенню і в рамках феноменологічної теорії має розглядатись як величина, що визначається з експерименту. Їхня відмінність зумовлена вкладом орбітальних моментів. Тому за умови, що магнітні іони знаходяться в $S$-стані (тобто орбітальні моменти відсутні), можна вважати, що $\gamma$ дорівнює гіромагнітному відношенню з високим степенем точності.. Це виконується для CdCr₂Se₄, залізо-ітрієвого гранату Y₃Fe₅O₁₂, пермалою Fe_20+xNi_80-x та більшості інших феро- та феримагнітних матеріалів.

Ефективне магнітне поле визначається як варіаційна похідна вільної енергії за магнітним моментом

$\mathbf {H} ^{\mathrm {eff} }(\mathbf {r} ,t)=-{\frac {\delta F}{\delta \mathbf {M} }}.\qquad (2)$

У випадку, коли розглядається магнетик далеко від температури Кюрі або за нульової температури, то вільна енергія $F$ дорівнює внутрішній $E$.

В формулюванні (1) зберігається довжина вектора намагніченості. Це легко показати, домноживши обидві частини (1) скалярно на $\mathbf {M}$, що дасть

${\frac {\partial \mathbf {M} ^{2}}{\partial t}}=0.\qquad (3)$

Цей факт дає підставу казати про прецесію намагніченості.

Строге виведення рівняння руху намагніченості в континуальному наближенні неможливий, тому часто постулюється можливість формального переходу від рівняння руху оператора спіну $\mathbf {S} _{n}$

$i\hbar {\frac {\partial \mathbf {S} _{n}}{\partial t}}=[{\mathcal {H}},\mathbf {S} _{n}],\qquad (4)$

до рівняння (1) шляхом заміни $\mathbf {S} _{n}\to -{\frac {a^{3}}{2\mu _{B}}}\mathbf {M} (\mathbf {r} _{n})$ і розкладу поля намагніченості $\mathbf {M} (\mathbf {r} _{n+n_{0}})$ поблизу точки $\mathbf {r} _{n}$ в ряд Тейлора. Тут $[\bullet ,\bullet ]$ — комутатор, ${\mathcal {H}}$ — гамільтоніан, $\mathbf {S} _{n}$ — оператор спіну для n-го вузла ґратки, а $\mathbf {r} _{n}$ — його радіус-вектор, $a$ — стала ґратки, $\mu _{B}$ — магнетон Бора.

Модифікації

Врахування дисипації, впливу температури чи спін-поляризованих струмів потребує модифікації вихідного рівняння (1), яка зазвичай зводиться до появи додаткових доданків в правій частині (1). Релаксаційні члени можуть мати різну розмірність і різну кількість параметрів. Але для наближеного опису процесів в феромагнетиках за невеликої дисипації може використовуватись рівняння в будь-якій з наведених нижче форм . Кожне з них можна перетворити з одного в інше.

Релаксаційний член в формі Ландау — Ліфшиця

Ландау та Ліфшиць запропонували наступну модифікацію:

${\frac {\partial \mathbf {M} }{\partial t}}=-|\gamma |[\mathbf {M} \times \mathbf {H} ^{\mathrm {eff} }]-{\frac {|\gamma |\lambda }{M^{2}}}\left[\mathbf {M} \times [\mathbf {M} \times \mathbf {H} ^{\mathrm {eff} }]\right],\qquad (5)$

де $\lambda$ — парметри дисипації. Інколи за параметр дисипації приймають величину $\lambda _{1}=|\gamma |\lambda$.

Рівняння Ландау — Ліфшиця — Гільберта

Часто використовується релаксаційний член в формі Гільберта:

${\frac {\partial \mathbf {M} }{\partial t}}=-|\gamma |[\mathbf {M} \times \mathbf {H} ^{\mathrm {eff} }]+{\frac {\alpha }{M}}\left[\mathbf {M} \times {\frac {\partial \mathbf {M} }{\partial t}}\right],\qquad (6)$

де $\alpha$ — параметр дисипації. Формальний перехід між рівняннями (5) та (6) можна здійснити заміною

$\gamma \to {\frac {\gamma }{1+\alpha ^{2}}},\quad \lambda \to {\frac {\alpha M}{1+\alpha ^{2}}}.\qquad (7)$

В зв'язку з від'ємним значенням гіромагнітного відношення зустрічаються визначення параметрів релаксації з протилежними знаками в (5) та (6) .

Рівняння Блоха — Бломергена

Прикладом рівняння з дисипацією, що допускає зміну довжини вектора намагніченості, може слугувати модифіковане рівняння Блоха чи рівняння Блоха — Бломергена:

${\frac {\partial \mathbf {M} }{\partial t}}=-|\gamma |[\mathbf {M} \times \mathbf {H} ^{\mathrm {eff} }]-\omega _{r}(\mathbf {M} -\chi _{0}\mathbf {H} ^{\mathrm {eff} }),\qquad (8)$

де $\chi _{0}$ — так звана статистична сприйнятливість, що визначається як відношення намагніченості насичення до абсолютної величини ефективного поля, а $\omega _{r}$ — частота релаксації.

Вплив спін-поляризованого струму

Спін-поляризований струм зазвичай описують додатковим доданком в правій частині (1) вигляду $|\gamma |\mathbf {T}$. Один з підходів до його конкретизації полягає в розкладі вектора $|\gamma |\mathbf {T}$ за осями, направленими вздовж $\mathbf {M}$, $[\mathbf {M} \times \mathbf {m} _{\mathrm {ref} }]$ та $\mathbf {M} \times [\mathbf {M} \times \mathbf {m} _{\mathrm {ref} }]$. Тут $\mathbf {m} _{\mathrm {ref} }$ — одиничний вектор вздовж намагніченості опорного шару. В припущенні, що довжина вектора намагніченості не змінюється, перша проєкція буде дорівнювати нулю, а дві інші

$\mathbf {T} _{\parallel }=-{\frac {|\gamma |a_{J}}{M_{s}}}\mathbf {M} \times [\mathbf {M} \times \mathbf {m} _{\mathrm {ref} }],\quad \mathbf {T} _{\perp }=|\gamma |b_{J}[\mathbf {M} \times \mathbf {m} _{\mathrm {ref} }],\qquad (9)$

де коефіцієнти $a_{J}$ та $b_{J}$ пропорційні густині струму, залежать від параметрів структури, що поляризує, та кута між $\mathbf {M}$ и $\mathbf {m} _{\mathrm {ref} }$.

Інші форми запису

Для аналітичного аналізу частіше за все рівняння Ландау-Ліфшиця записується в кутових змінних сферичної системи координат $\theta$ та $\phi$. В такому випадку вектор намагніченості можна представити як

$M_{x}+M_{y}=M_{s}\sin \theta e^{i\phi },\quad M_{z}=M_{s}\cos \theta ,$

де $M_{s}$ — намагніченість насичення. Щоб перейти в (1) до кутових змінних, домножимо рівняння на варіацію намагніченості $\delta \mathbf {M}$, виразивши в кутових змінних проєкцію лівої частини на вісь аплікат. Далі, після запису варіації енергії та намагніченості через варіації кутів, отримаємо

$\sin \theta {\frac {\partial \theta }{\partial t}}=-{\frac {|\gamma |}{M_{s}}}{\dfrac {\delta E}{\delta \phi }},\quad \sin \theta {\frac {\partial \phi }{\partial t}}={\frac {|\gamma |}{M_{s}}}{\dfrac {\delta E}{\delta \theta }}.\qquad (10)$

Отримання рівнянь в кутових змінних, що містять додаткові члени, відбувається аналогічно. Так, для запису в формі Ландау — Ліфшиця — Гільберта маємо

$\sin \theta {\frac {\partial \theta }{\partial t}}=-{\frac {|\gamma |}{M_{s}}}{\dfrac {\delta E}{\delta \phi }}-\alpha \sin ^{2}\theta {\frac {\partial \phi }{\partial t}},\quad \sin \theta {\frac {\partial \phi }{\partial t}}={\frac {|\gamma |}{M_{s}}}{\dfrac {\delta E}{\delta \theta }}+\alpha {\frac {\partial \theta }{\partial t}}.\qquad (11)$

Література

• Ахиезер, А. И., Барьяхтар, В. Г. Пелетминский, С. В. Спиновые волны., М.: Наука, 1967, — 368 с.
• Гуревич, А. Г., Мелков, Г. А. Магнитные колебания и волны. М.: Физматлит, 1994. — 464 с.,  ISBN 5-02-014366-9.
• Зависляк, И. В., Тычинский, А. В., Физические основы функциональной микроэлектроники. К.: УМК ВО, 1989, — 105. с.
• Звездин, А. К, Звездин, К. А, Хвальковский, А. В. Обобщенное уравнение Ландау — Лифшица и процессы переноса спинового момента в магнитных наноструктурах. УФН, 178 436—442, (2008) http://dx.doi.org/10.3367/UFNr.0178.200804i.0436
• Ландау, Л. Д., Лифшиц, Е. М. К теории дисперсии магнитной проницаемости ферромагнитных тел. Phys. Zs. Sowjet., 1935, 8, С. 153-169.
• Скроцкий, Г. В. Еще раз об уравнении Ландау — Лифшица. УФН
• Gilbert, T. A phenomenological theory of damping in ferromagnetic materials. IEEE Transactions on Magnetics, 2004, 40, pp. 3443-3449. http://dx.doi.org/10.1109/TMAG.2004.836740
• Hubert, Alex; Rudolf Schäfer (1998). Magnetic domains: the analysis of magnetic microstructures. Springer. с. 557. ISBN 3540641084.