CC BY
60
Вестник Челябинского государственного университета. 2011. № 39 (254). Физика. Вып. 12. С. 15-18.
Л. Н. Котов, Л. С. Носов, А. В. Голов, В. А. Устюгов
О МАГНИТНОЙ ПЕРЕОРИЕНТАЦИИ СВЧ-ПОЛЕМ В АНТИФЕРРОМАГНИТНЫХ НАНОЧАСТИЦАХ1
Получено численное решение уравнения Ландау—Лифшица в форме Г ильберта, описывающее движение вектора антиферромагнетизма в однодоменной антиферромагнитной частице под действием СВЧ-полей при нулевом постоянном магнитном поле. Определены зависимости пороговой амплитуды поля, при которой возникает магнитная переориентация, от частоты, параметра затухания, угла между переменным полем и направлением намагниченностей подрешёток. На частотной зависимости пороговой амплитуды поля выявлен глубокий минимум, который может быть использован для уменьшения энергии поля, необходимой для изменения магнитной структуры в антиферромагнетиках.
Ключевые слова: антиферромагнетик, СВЧ, перемагничивание.
Интерес к нелинейной магнитной динамике в магнитоупорядоченных средах обусловлен разнообразием нелинейных эффектов, возникающих при воздействии на диссипативную магнитную систему ВЧ- или СВЧ-полем [1-8]. Нелинейные эффекты, которые возникают при больших углах прецессии вектора намагниченности, можно использовать, например, для модуляции электромагнитного излучения, эффективность которой определяется величиной угла прецессии [6]. В данной работе методом численного моделирования исследовано явление переориентации векторов намагниченностей подрешёток в антиферромагнитных однодоменных наночастицах, возникающей во время действия на них СВЧ-магнитного поля. При моделировании не учитывались процессы параметрического распада и ограничения на максимальные углы прецессии магнитных моментов. Прецессия неоднородных колебаний, например, суловские неустойчивости, обусловленных трёх- и четырёх-магнонными взаимодействиями, возникать не будет, если частота прецессии совпадает с минимально возможной частотой спектра спиновых волн юк, где к — волновой вектор спиновой волны, т. е. соответствует дну спин-волновой «зоны» [6-7]. Для антиферромагнетиков, в соответствии с [7-8], частота неоднородных колебаний по порядку величины юк ~ у • q • М • k2 , где у — гиромагнитное отношение; q — параметр неоднородного внутриподрешёточного обмена; М — намагниченность подрешёток. Если волновой вектор оценить как k ~Ш, где d < 100 нм — размер однодоменной антиферромагнит-ной частицы, то частота неоднородных колеба-
1 Работа выполнена при поддержке гранта РФФИ
№ 10-02-01327.
ний может быть выражена как юк « у • q • М / d2 [7]. Если считать, что значения параметров внутри- и межподрешёточного обмена одного порядка, то получается, что юк примерно в 104 раз больше частоты внешнего поля, которое близко частоте антиферромагнитного резонанса
(АФМР) *у>/НА (М + НА) [^ где НА —
поле кубической анизотропии; X — параметр обмена между подрешётками. Величина амплитуды переменного магнитного поля, используемого в численном моделировании в 100 раз больше величины поля анизотропии. Следовательно, для антиферромагнитых однодомённых наночастиц, основная (однородная) мода спин-волнового спектра, соответствующая дну «зоны», оказывается достаточно далеко удалена (по частоте в 104 раз) от первой (неоднородной) спин-волновой моды. В этом случае амплитуда поля возбуждения неоднородных колебаний является наивысшей, и поэтому в решаемой задаче с ростом амплитуды СВЧ-поля не будет происходить явления распада однородной прецессии на неоднородные колебания вектора намагниченности [5].
Для переориентации векторов намагниченностей подрешёток в частице из одного в другое устойчивое положение необходимо изменить угол между начальным и конечным положением вектора, например, на 90° [1-3]. Ранее магнитная переориентация в антиферромагнит-ных частицах изучалась только в импульсных и постоянных полях большой амплитуды [6-7]. Явление АФМР при магнитной переориентации, в условиях которого возможно значительное уменьшение амплитуды магнитных полей, при этом не учитывалось, поэтому возможно значительное уменьшение амплитуды внешних полей, поэтому необходимо выявить влияние
АФМР на наименьшую (пороговую) амплитуду поля различной частоты, при которой наблюдается переориентация векторов намагниченностей подрешёток в антиферромагнитных частицах с разным значением параметра затухания. В работах [4-5] исследовалась нелинейная динамика двухслойной обменно-связанной магнитной структуры с антиферромагнитной связью, где было показано, что в результате воздействия переменного магнитного поля возможно наблюдение динамического магнитного фазового перехода [4] или нелинейной стохастической динамики [5]. В данной работе также учитывается ан-тиферромагнитная связь между слоями, поэтому здесь использован подобный подход к решению задачи переориентации векторов намагниченностей подрешёток, как и в [4-5]. Также, для упрощения расчётов, использована разработанная ранее численная методика решения уравнений магнитной динамики в кубических ферромагнетиках [2-3]. Поэтому в данной работе рассмотрена антиферромагнитная однодоменная сферическая частица кубической структуры. Плотность магнитной энергии частицы представим в виде суммы плотностей энергий обменного поля, кубической анизотропии (в предположении однои-онного источника анизотропии) и зеемановской энергии намагниченностей подрешёток частицы в переменном магнитном поле [8]:
U (mi, m2 ) = -ШДО2 + Uan (mi ) +
+Uan (m 2 ) M (ml + m 2 )• H -, (1)
где
Uan (m) = K1 ' (m2 + m2mZ + m)mZ ) -
+K2m2x m2ym2z — плотность энергии кубической анизотропии (подрешётки считаем эквивалентными); К1 > 0 и К2 < 0 — первая и вторая константы кубической анизотропии; H~ = h • sin (со - f) — переменное магнитное поле линейной поляризации с частотой ю; m. = = M JM — вектор направляющих косинусов намагниченности Mi 7-й подрешётки, M = |M.|; оси x, у и z совпадают с кристаллографическими осями частицы [100], [010] и [001]. Пусть вектора ферромагнетизма и антиферромагнетизма равны m = m1 + m2 и l = = m, - m„ соответственно [8]. Зависимость U =
1 2 L J an
= Uan(m) от направления вектора m для рассматриваемого случая аналогична зависимости для ферромагнитных частиц, описанных в работе [2], приводящей к появлению взаимно перпендикулярных осей лёгкого намагничивания [8]. Особенностями зависимости плотности энер-
гии кубической анизотропии от направления вектора намагниченности является наличие трёх типов особых точек. Точки минимума и максимума магнитной энергии характеризуют устойчивое и неустойчивое равновесие векторов m и 1 соответственно. В этих точках эффективное поле анизотропии и вектора намагниченности параллельны. Если же вектора m и 1 параллельны осям [110], то энергия соответствует неустойчивому равновесию и особым точкам типа «седло» [9]. Во время колебаний вектора 1 около этого положения возможна переориентация вектора антиферромагнетизма 1. Эти утверждения для рассматриваемого антиферромагнетика справедливы и в том случае, когда рассматриваются обе подрешётки одновременно, а вектор ферромагнетизма т = 0. При отсутствии переменного поля вектор 1 занимает одно из трёх устойчивых положений, соответствующих минимумам энергии: параллельно осям х, у или г. Уравнение движения векторов намагниченностей подрешёток может быть записано в форме Гильберта с учётом модификаций [8]:
^т,
dt
:-[mi Х H(eff ]■
+a
m x
dmi
dt*
i = 1,2,
(2)
где а — безразмерный параметр диссипации;
НеЦ =Ьи*/дт1 — приведённое эффектив-
ное магнитное поле, действующее на намагниченность 7-й подрешётки М где и = и(т1, т2)/2К — приведённая плотность свободной энергии; t = t • юА — приведённое время; ю = = ю/юА — приведённая частота, где юА = уНА = = 2уК1/М — частота колебаний вектора намагниченности в условиях естественного ферромагнитного резонанса (ФМР). Приведённая амплитуда переменного магнитного поля Н* = Ь/НА = = Ь • М/2К1. Такие переменные удобны для численного моделирования и сравнения их с экспериментом. Например, частота естественного АФМР в таких переменных при а = 0 равна
+ 2ХМ2 /2КХ . При всех расчётах считалось, что постоянное магнитное поле равно 0, и использовались следующие параметры: время наблюдения = 2, время действия СВЧ-поля т* = 1, К2/К1 = -0,5, М2/2К1 = 9, у = 103. Численные расчёты проводились по методу Рунге—Кутты 7-8-го порядка [10]. Пусть до включения поля Н, направленного параллельно плоскости ху, вектор
1 направлен вдоль оси у. При определённых значениях амплитуды и частоты поля Н~ возникает переориентация вектора 1 в положение, перпендикулярное начальному положению (рис. 1).
Это явление носит пороговый характер. Пусть пороговая амплитуда переориентации И* соответствует минимальной амплитуде Н-поля на фиксированной частоте, при которой происходит переориентация вектора 1 и он ориентируется вдоль оси г или х. Зависимость пороговой
7 * .. ^ *
амплитуды ^ от приведённой частоты ю для антиферромагнитной частицы, которая подобна такой же зависимости для ферромагнитной частицы [2-3], приведена на рис. 2.
Интересно отметить, что зависимость ^ (га ) близка к зависимости пороговой амплиту-
ды возбуждения нелинейного ФМР от частоты [6]. Особенностью поведения зависимости h* (<в*^ является наличие низкочастотных минимумов, обусловленных возбуждением нелинейных мод из-за связи колебаний подрешёток. Наблюдающийся широкий минимум на зависимости h*b (a>*j с ростом параметра затухания а смещается в область низких частот. Амплитуда h* в области минимума увеличивается с ростом а, как и для ферромагнитных частиц [2]. Как видно из рис. 3, значение h*b зависит от угла в между направлениями поля h* и вектора 1. Начиная с малых углов и до в ~ п/4 амплитуда h*b медленно уменьшается, а затем возрастает для всех частот переменного поля (рис. 3). При значении угла в ~ п/2 переориентация
Рис. 1. Траектория движения вектора антиферромагнетизма I в сферической частице под действием переменного магнитного поля (переориентация из положения 1 в перпендикулярное положение 2). т* = 67, Ь* = 141, в = к/4
ю
7 *
Рис. 2. Зависимость пороговой амплитуды \ переориентации вектора I в антиферромагнитной однодомённой частице от приведённой частоты т* при параметрах диссипации а: 1 — 0,01, 2 — 0,005, 3 — 0,001, в = к/4
Р
Рис. 3. Зависимость пороговой амплитуды h* переориентации вектора l в антиферромагнитной однодоменной частице от угла в при а = 0,01 на частотах т*:
1 — 40, 2 — 50, 3 — 60, 4 — 70
не наблюдается. Это связанно с тем, что некол-линеарное внешнему полю положение вектора намагниченности для антиферромагнетика является энергетически более выгодным, чем коллинеарное, в отличие от ферромагнетика.
Таким образом, в данной работе показана переориентация вектора l в антиферромагнит-ных однодоменных наночастицах при условии, что при больших амплитудах поля не возникает параметрического распада. Показано, что на зависимости h* для кубических анти-
ферромагнетиков имеется глубокий минимум, который может быть использован для уменьшения амплитуды CBЧ магнитного поля при переориентации вектора l в антиферромаг-нитных частицах. Поэтому необходимо синтезировать такие антиферромагнетики, для которых характерны области минимумов пороговых амплитуд CBЧ-полей. Явление переориентации позволяет изменять магнитную структуру в антиферромагнитных наночастицах и в ансамблях наночастиц и может быть использовано для создания CBЧ переключающих устройств. Кроме того, процесс CBЧ-переориентации вектора намагниченности в антиферромагнитных частицах характеризуется значительно меньшим временем, чем в ферромагнитных частицах.
Список литературы
1. Котов, Л. Н. Переориентация вектора намагниченности в однодоменной частице импульсом высокочастотного поля / Л. Н. Котов, Л. C. Носов
// Письма в Журн. техн. физики. 2003. Т. 29, № 20. С. 38-42.
2. Котов, Л. Н. Переориентация намагниченности в однодоменных частицах и отклик на импульс поля / Л. Н. Котов, Л. С. Носов // Журн. техн. физики. 2005. Т. 75, № 10. С. 55-60.
3. Котов, Л. Н. Изменение магнитной структуры ансамблей однодоменных частиц и их отклик на радиоимпульс поля / Л. Н. Котов, Л. С. Носов // Журн. техн. физики. 2008. Т. 78, № 5. С. 60-65.
4. Шутый, А. М. Динамический ориентационный фазовый переход в двухслойной магнитосвязанной структуре / А. М. Шутый, Д. И. Семенцов // Письма в Журн. эксперимент. и теорет. физики. 2002. Т. 75, № 5-6. С. 287-290.
5. Шутый, А. М. Стохастическая динамика намагниченности в обменносвязанной слоистой структуре / А. М. Шутый, Д. И. Семенцов // Письма в Журн. эксперимент. и теорет. физики. 2003. Т. 78, № 8. С. 952-956.
6. Моносов, Я. А. Нелинейный ферромагнитный резонанс. М. : Наука, 1971. 210 с.
7. Львов, В. С. Нелинейные спиновые волны. М. : Наука, 1987. 272 с.
8. Гуревич, А. Г. Магнитный резонанс в ферритах и антиферромагнетиках. М. : Наука, 1973. 464 с.
9. Мигулин, В. В. Основы теории колебаний / В. В. Мигулин, В. И. Медведев, Е. Р. Мустель, В. Н. Парыгин. М. : Наука, 1973. 392 с.
10. Бордовицина, Т. В. Современные численные методы в задачах небесной механики. М. : Наука, 1984. 136 с.
