УДК 548:537.611.44
Способ определения по экспериментальной кривой намагничивания энергетических вкладов одноосной и кубической анизотропии в полную энергию магнитной анизотропии наночастиц
Сергей В. Комогорцев*
Институт физики СО РАН, Академгородок, 50/38, Красноярск, 660036,
Россия
Сибирский государственный технологический университет,
Мира, 82, Красноярск, 660049,
Россия
Сергей И. Смирнов^
Красноярский государственный педагогический университет,
Лебедевой, 89, Красноярск, 660049,
Россия
Надежда А. Момот*
Сибирский федеральный университет, Свободный, 79, Красноярск, 660041,
Россия
Рауф С.Исхаков§
Институт физики СО РАН, Академгородок, 50/38, Красноярск, 660036,
Россия
Получена 18.08.2010, окончательный вариант 05.09.2010, принята к печати 25.09.2010 На основе результатов численного моделирования в рамках композиционной модели Стоне-ра-Вольфарта предложены способы определения по экспериментальной кривой намагничивания энергетических вкладов одноосной и кубической анизотропии в полную энергию магнитной анизотропии частиц. Предложенные способы апробированы на примере экспериментальных кривых намагничивания наночастиц Оо-И''.
Ключевые слова: численное моделирование, ферромагнитные наночастицы, магнитная анизотропия, модель Стонера-Вольфарта.
Кривые намагничивания ансамблей наночастиц с размерами, близкими к критическому размеру однодоменности, характеризуются наибольшими значениями величины коэрцитивной силы и остаточной намагниченности, а также временной и температурной стабильностью этих параметров. Перемагничивание однодоменных частиц осуществляется посредством процессов вращения намагниченности и описывается на основе модели Стонера-Вольфарта [1]. Задача данной работы — интерпретация экспериментальных кривых намагничивания наночастиц Ni-Co, капсулированных в углеродных нанотрубках с помощью
§ [email protected] © Siberian Federal University. All rights reserved
численного моделирования этой системы на основе модели Стонера-Вольфарта. Наночасти-цы N1-00, капсулированные в углеродных нанотрубках, предварительно были исследованы в работе [2]. Здесь установлены особенности морфологии, структурного состояния, а также определены величины константы обменного взаимодействия данных частиц. Было показано, что размеры исследуемых систем близки к критическому размеру однодоменности, а величины коэрцитивной силы и остаточной намагниченности незначительно меняются в области температур (78 ^ 300 К). Эти данные указывают на возможность интерпретации основных параметров кривых намагничивания исследуемых наночастиц №-Со в рамках модели Стонера-Вольфарта. Методом рентгеновской дифракции установлено, что магнитные наночастицы Со1-ж№ж характеризуются ГЦК структурой. По величине уширения дифракционных рефлексов на рентгенограммах выполнена оценка размеров наночастиц В, результаты которой представлены в табл. 1.
Таблица 1. Магнитные параметры наночастиц Coi_x
D, nm Hc, Oe 300 K 77 K Mr /Ms 300 K 77 K ku 300 K 77 K Keff, 105 300 K erg/cm3 77 K
Co 15 630 680 0.62 0.62 0.26 0.26 36.1 38.4
CO80NÍ20 22 600 670 0.69 0.76 0.23 0.20 27.1 28.1
C050NÍ50 29 480 620 0.66 0.66 0.25 0.25 24.3 27.4
C020NÍ80 37 300 390 0.64 0.67 0.25 0.24 15.2 16.1
Ni 31 140 290 0.59 0.73 0.28 0.21 8.5 10.5
Обратимая часть кривой намагничивания в области приближения намагниченности к насыщению для всех образцов изучаемых частиц Со1_ж№ж следует степенной зависимости М ~ Н-2 (закон Акулова, см.[7]):
М(Н)/Ма = 1 - (2аКе///МН)2 , (1)
где Ма — намагниченность насыщения, Кец — константа эффективной анизотропии, а — симметрийный коэффициент. Выполнение этой зависимости (рис. 1) свидетельствует как об отсутствии обменных взаимодействий между исследуемыми частицами, так и о ведущей роли процессов неоднородного вращения намагниченности при намагничивании системы частиц [4].
Форма экспериментальных кривых намагничивания не соответствует по отдельности форме петель с одноосной либо с кубической магнитной анизотропией, рассчитанных в стандартной модели Стонера-Вольфарта. Величины остаточной намагниченности тг = Мг/Ма для всех образцов наночастиц Со1-ж№ж (табл. 1, пример на рис. 2) оказываются большими, чем тг =0,5 (положительная одноосная магнитная анизотропия) и меньшими, чем тг = 0, 86 (отрицательная кубическая магнитная анизотропия), характерными для монокристаллов Со1-ж№ж [5]. Экспериментальные значения величин остаточной намагниченности для наночастиц могут быть получены в рамках комбинированной модели Стонера-Вольфарта, в которой полная энергия магнитной анизотропии однодоменной частицы имеет одновременно и кубическую, и одноосную составляющие [6]. Необходимость использования данной модификации модели Стонера-Вольфарта при интерпретации экспериментальных кривых намагничивания исследуемых наночастиц становится очевидной в связи с установленными ранее особенностями морфологии и структуры наночастиц Со-№, а именно: несферичность формы частиц, а также наличие дефектов упаковки, вблизи которых атомы характеризуются ГПУ ближайшим окружением.
В композиционном варианте модели Стонера-Вольфарта, рассмотренном в [6], энергия магнитной анизотропии каждой наночастицы складывается как из энергии случайно ориен-
Рис. 1. Экспериментальные кривые приближения намагниченности к насыщению наночастиц Co1-xNix при температуре 77 К
Рис. 2. Типичная петля гистерезиса, поправленная на размагничивающее поле, на примере наночастиц Co80Ni20 при температуре 300 К
тированной кубической магнитокристаллической анизотропии, так и из энергии случайно ориентированной одноосной магнитной анизотропии. Полная магнитная энергия одной частицы в композиционном варианте модели Стонера-Вольфарта запишется в виде
Е = (К 1 + Ки) [киа2и3 + (1 - ки)(а\а\ + а\а| + а^а?) - Ъ^афх + а^ + азвз)] (2)
где ки = Ки/(К 1 + Ки) и Н = ШаН/ (К1 + Ки); М3 — намагниченность насыщения одной частицы; Н — внешнее магнитное поле; Ки и К1 — константы одноосной и кубической магнитной анизотропии частицы; аиз — направляющий косинус орта намагниченности по отношению к легкой оси одноосной анизотропии; а1, а2, аз — направляющие косинусы орта локальной намагниченности по отношению к одной из осей выбранной системы координат, жестко связанной с одной из осей кубической анизотропии; във2,вз — направляющие косинусы вектора внешнего магнитного поля.
1.0000
■ 0,1
• 0.2,0.8
▲ 0.4,0.6
★ 0.6,0.4
□ 0.8,0.2
о 1,0
0.000 0.002 0.004 0.006 0.008 0.010 h-2
Рис. 3. Изменение формы петли гистерезиса в модели с комбинированной симметрией магнитной анизотропии
Рис. 4. Приближение намагниченности к насыщению в композиционном варианте модели Стонера-Вольфарта
Пример расчетных кривых намагничивания, полученных с использованием энергии в
форме (2), приведен на рис. 3. При ки = 0 либо 1 получаются классические кривые модели Стонера-Вольфарта, характеризующиеся хорошо известными параметрами. В композиционном варианте модели Стонера-Вольфарта (ки = 0, 5) видна модификация формы кривой намагничивания, а также основных параметров петли гистерезиса. В работе [6], посвященной численному моделированию, были рассчитаны кривые намагничивания для случая Ки > 0, К1 < 0, адекватного нашей экспериментальной ситуации. Расчетные графики приближения намагниченности к насыщению в композиционном варианте модели Стонера-Вольфарта для различных вкладов одноосной и кубической анизотропии в осях (т, Л.-2), где т = М/Ма соответствуют закону Акулова (рис. 4):
m(h)
aHa 2
1- =1-
H
2aKeff L(K„ + Ki)hJ
1 -
ao L h J
2
(3)
Симметрийный коэффициент а в законе Акулова для системы однодоменных частиц со случайной одноосной анизотропией равен у/1/15, а для частиц с кубической анизотропией —^ 2/105. При наличии комбинированной симметрии магнитной анизотропии коэффициент ао должен зависеть от относительной доли вклада энергии одноосной анизотропии в результирующую энергию анизотропии. Как оказалось, зависимость а0 (ки) имеет нелинейный характер и может быть хорошо аппроксимирована квадратичным полиномом [6]:
а2(ки)=0, 37ки - 0,17ки + 0,08. (4)
Для предельных случаев одноосной и кубической анизотропии симметрийный коэффициент a2 принимает значения 4/15 и 8/105, впервые вычисленные для одноосных и кубических кристаллов в [3]. Используя результаты численного моделирования в рамках композиционной модели Стонера-Вольфарта, по экспериментальной кривой намагничивания, предварительно поправленной на размагничивающее поле (рис. 2), можно определить как полную энергию магнитной анизотропии исследуемых частиц Ni-Co, так и энергетические вклады одноосной и кубической анизотропии в полную энергию магнитной анизотропии частиц. Сопоставляя экспериментальные значения остаточной намагниченности частиц (табл. 1) с величинами остаточной намагниченности, рассчитанной в рамках композиционной модели Стонера-Вольфарта, мы определили величину ku = Ku/(K1 + Ku) — вклад энергии одноосной магнитной анизотропии в полную энергию магнитной анизотропии (табл. 1). Далее, используя полученное значение ku, мы вычислили величину симметрийного коэффициента a в законе Акулова для каждого исследуемого образца. По тангенсам наклона прямолинейных участков экспериментальных кривых приближения намагниченности к насыщению (рис. 1) определили произведение величины полной энергии анизотропии на симметрийный коэффициент, т.е. величину aKeff. Мы использовали вычисленную величину коэффициента a для определения величины Keff (табл. 1).
Зависимость коэрцитивной силы (приведенной к полю анизотропии Hc/Ha) от размера частицы (приведенного к величине критического размера однодоменности D/Ro) в двойных логарифмических координатах характеризуется максимумом (вблизи D/Ro = 1, 3 — критический размер однодоменности). Величина критического размера однодоменности для
i a
каждого образца вычислялась по формуле R0 = 1,44*/——г для сферической частицы [7], с
V M2
использованием величин магнитных констант, определенных из магнитометрических измерений [2]. Величина поля магнитной анизотропии частицы рассчитана как Ha = 2Kef f /Ms. На рис. 5 также дано сопоставление наших данных с результатами работы [8] по магнитным наночастицам CosoNi2o. Наблюдаемое смещение позиции максимума зависимости Hc/Ha от 1 (D = R0) на рис. 5, на наш взгляд, обусловлено несферической формой исследуемых частиц.
2
0.16-,
0.14-
■ 300 к
• СоЛИ
▲ 77 К
0.12-
0.10-
0.08-
0.06-
0.04-
0.02-
Рис. 5. Зависимость коэрцитивной силы (приведенной к полю анизотропии) от размера частиц (приведенного к величине критического размера однодоменности)
Таким образом, предложен новый способ определения энергетических вкладов одноосной и кубической анизотропии в полную энергию магнитной анизотропии частиц на основе результатов численного моделирования экспериментальной кривой намагничивания. Применимость метода продемонстрирована на примере набора наночастиц N1-00 переменного состава.
Работа выполнена в рамках федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009--2013 годы.
Список литературы
[1] E.C.Stoner, E.P.Wholfarth, A mechanism of magnetic hysteresis in heterogeneous alloys, Phil. Trans. Roy. Soc. A, 240(1948), 559-564.
[2] S.V.Komogortsev, R.S.Iskhakov, K.A.Shayhutdinov, et.al., Magnetic properties of NixCo1-x Nanoparticles in Carbon Nanotubes, The Phys. Met. Met., 102(2006), Suppl. 1, S67--S70.
[3] Н.С.Акулов, Л.В.Киренский, Магнитный метод исследования внутренних напряжений в ферромагнитных металлах, ЖТФ, 9(1939), №13, 1145-1150.
[4] С.В.Комогорцев, Р.С.Исхаков, А.Д.Балаев и др., Магнитные свойства ферромагнитных наночастиц Fe3C, капсулированных в углеродных нанотрубках, ФТТ, 49(2007), вып.49, 700-703.
[5] Magnetic properties of metals, Edited by H.P.J.Wijn, Springer-Verland Berlin Heidelberg, 1991.
[6] S.I.Smirnov, S.V.Komogortsev, Magnetization curves of randomly oriented ferromagnetic singledomain nanoparticles with combined symmetry of magnetic anisotropy, J. Magn. Magn. Mat., 320(2008), 1123--1127.
[7] С.В.Вонсовский, Магнетизм, Наука, Москва, 1971.
[8] G.V.Kurlyandskaya, S.M.Bhagat, C.Luna, M.Vazquez, Microwave absorption of nanoscale CoNi powders, Jour. Appl. Phys, 99(2006), 104308.
Technique for Determination the Energy of Uniaxial and Cubic Magnetic Anisotropy in Magnetic Nanoparticles from Experimental Magnetization Curve
Sergey V. Komogortsev Sergey I. Smirnov Nadezhda A. Momot Rauf S. Iskhakov
Technique for determination the energy of uniaxial and cubic magnetic anisotropy in magnetic nanoparticles from experimental magnetization curve is proposed. The technique is based on the numerical simulation of experimental magnetization curve in the terms of the composite mode of the Stoner- Wolfarth model. Investigation of magnetic anisotropy for Co-Ni nanoparticles is carried out by proposed technique.
Keywords: numerical modeling, ferromagnetic nanoparticles, magnetic anisotropy, Stoner-Wolfart model.