Сетевое издание

Современные проблемы науки и образования

ISSN 2070-7428

"Перечень" ВАК

ИФ РИНЦ = 1,006

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Губайдуллин Р.Р. 1 Гареева М.Я. 1 Максимочкин В.И. 2

---

1 ФГБОУ ВПО «Стерлитамакский филиал Башкирского государственного университета»

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

Большинство естественных минералов содержат дефекты структуры, вакансии, примесные атомы и прочие включения, присутствие которых влияет на магнитные характеристики систем. В работе представлено теоретическое исследование влияния неоднородности состава на магнитные свойства естественных ферримагнетиков. Предложена модель магнетика, как аналога минерала магнетита, примесные включения в котором размещались равномерно или по нормальному распределению, создавая тем самым однородную и неоднородную по составу конфигурации. Гетерогенность, как показано в работе, расширяет диапазон составов с аномальными кривыми термонамагниченности. Установлено, что неоднородность состава увеличивает значение намагниченности и температуры Кюри, что находит подтверждение в экспериментальных работах с горными породами. В частности, окисление приводит к нарушению стехиометрии, вызванное диффузией ионов металлов в неподвижном кислородном каркасе. Возможно также и применение модели к восстановительным процессам с магнетиками, находящимся в природной углеводородной среде.

Статья в формате PDF

248 KB

ферримагнитные минералы

термонамагничивание

модель Изинга

1. Белов Н.В. Очерки структурной минералогии. - М. : Недра, 1976. - 323 с.

2. Кравченко С.Н., Орлюк М.И., Русаков О.М. Новый подход к интерпретации региональной Западно-Черноморской магнитной аномалии // Геофизический журнал. - 2003. - Т. 25. - № 2. - С. 135-143.

3. Néel L. Propriétésmagnétiques des ferrites: ferrimagnétisme et antiferromagnétisme // Ann. Phys. - 1948. - V. 3. - № 2. - P. 137-198.

4. Ozima M., Larson E.E. Study of Self-reversal of TRM in Some Submarine Basalts // Journal of geomagnetism and geoelectricity. - 1968. - V. 20. - № 4. - P. 337-351.

5. Ozima M., Ozima M. Self-reversal of remanent magnetization in some dredged submarine basalts // Earth and Planetary Science Letters. - 1967. - V. 3. - P. 213-215.

6. Uyeda S. Magnetic Interaction between Ferromagnetic Materials Contained in Rocks // Journal of geomagnetism and geoelectricity. - 1955. - V. 7. - № 1-2. - P. 9-36.

Неоднородность по составу, наличие примесей и различных дефектов являются характерными особенностями естественных минералов и горных пород различной природы. Эти особенности могут серьезно влиять на магнитные свойства минералов и создавать различные аномалии в процессах намагничивания.

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

В настоящей работе представлено исследование влияния неоднородности состава магнитных минералов на примере трехмерной модели Изинга. Была создана модель двух-подрешеточного ферримагнетика, как аналог минерала магнетита, в которой примесные включения занимали места в одной из подрешеток.

Описание модели

Модель представляет собой кубическую решетку размером L×L×L, содержащую N=L³ узлов. Все узлы разделены на два типа A и B таким образом, что ближайшими соседями узлов одного типа являлись узлы другого типа. Магнитные спины S_A и S_B, принимающие значения ±1 для S_A и ±2 для S_B, занимают соответственно места A и B в решетке. Примесные включения C, не имеющие магнитного спина, занимают места B и заменяют, таким образом, спины S_B. Общая «молекулярная формула» для модели выглядит следующим образом:

S_A∙S_B(1-x)∙C_x,

где x - задает концентрацию примесных включений в подрешетке B.

Для задания однородной и неоднородной по составу конфигурации в модели применяются два вида распределения примесей по решетке -равномерное и нормальное.

Полная магнитная энергия системы определяется гамильтонианом:

где суммирование происходит для различных пар ближайших соседей, J_AB-константа обменного взаимодействия между магнитными спинами, которую можно также назвать обменным интегралом между подрешетками A и B. Взаимодействия с константами J_AA и J_BB отсутствуют ввиду отсутствия ближайшего соседства спинов одного типа, а дальний порядок в моделировании не принимался во внимание. Для устранения влияния поверхностных эффектов на систему накладывались периодические граничные условия.

Суммарная намагниченность определяется:

где S - среднее значение магнитного момента на спин, µ - магнетон Бора.

Для достижения термодинамического равновесия системы при заданной температуре T применялся следующий алгоритм (алгоритм Метрополиса).

1. Формируется начальная конфигурация с заданными значениями L, x и T, а также видом распределения примесей.
2. Производится случайное пробное изменение в начальной конфигурации путем переворота какого-нибудь спина.
3. Вычисляется изменение энергии системы, обусловленное произведенным изменением конфигурации.
4. Если ΔE ≤ 0, то конфигурация принимается. Если ΔE ≥ 0, то конфигурация принимается с вероятностью exp(-ΔE/kT).

Шагом моделирования методом Монте-Карло (МК) называется последовательность пунктов 2-4 алгоритма. Равновесное состояние начальной конфигурации достигалось путем повторения 1000 раз на спин шага МК. После чего определялось среднее значение модуля магнитного момента на спин по 10 шагам МК.

Конечное абсолютное значение магнитного момента на спин (S=M/µN) для заданных L, x,T и вида распределения примесей определялось как среднее по 10 случайно сформированным конфигурациям.

Симуляция производилась при размерах стороны кубических решеток L=10, 20, 30, в интервале температур T от 0 до 10. Концентрация x принимала значения 5, 10, 15 ... 75.

Модель реализована в программе на языке C++ с использованием компилятора g++ в среде Linux. Приводимые данные обезразмерены.

Результаты исследования и их обсуждение

На рисунке 1 показано сравнение результатов моделирования температурных зависимостей магнитного момента на спин для нормального и равномерного распределения примесей при L=10 и L=30.

Рис. 1. Результаты моделирования при x = 0.35 для решеток с L=30 и L=10 при нормальном и равномерном распределениях примесных включений.

Видно, что кривые L=10 (норм.) и L=10 (равн.) практически совпадают, в то время как кривые L=30 (норм.) и L=30 (равн.) различны. Кроме того, кривые при L=30 имеют гораздо большую точность. Другими словами, на решетках малых размеров нельзя создать различия в нормальном и равномерном распределении, что говорит о предпочтении моделирования на решетке с L=30. Поэтому дальнейшие результаты показаны для решеток размером 30×30×30.

Результаты моделирования для двух распределений и различных концентраций примесей показаны в таблице 1, на рисунках 2 и 3.

Таблица 1 - Результаты моделирования для решетки L=30 при равномерном (Равн.) и нормальном (Норм.) распределении примесей (в скобках значения температур соответствующих S_max)

Температурные зависимости магнитного момента принимают различные типы, описанные в теории ферримагнетизма Нееля [3]. Максимальные значения магнитный момент для кривых P-типа и N-типа достигает практически при одной и той же температуре (T=5) для обоих распределений, что характеризует S_max как универсальный индекс, практически не зависящий от x. Наличие точки компенсации T_k, напротив, зависит как от степени разбавления примесями, так и от вида распределения немагнитных включений. Таким образом, самообращение намагниченности минерала может присутствовать и в более широком диапазоне составов в зависимости от вида беспорядка немагнитных атомов.

При увеличении х наблюдается уменьшение температуры Кюри и начального магнитного момента. Для равномерного распределения примесных включений стоит отметить стабильное уменьшение точки Кюри, в то время как при нормальном распределении температура Кюри изменяется незначительно, а также начальный магнитный момент S₀ (T=0) имеет более высокие значения в сравнении с равномерным порядком примесей.

Рис. 2. Результаты моделирования для решетки L=30 при равномерном распределении примесей для различных значений концентрации x.

Рис. 3. Результаты моделирования для решетки L=30 при нормальном распределении примесей для различных значений концентрации x.

На рисунке 4 показаны термозависимости абсолютного значения магнитного момента при x=35 для нормального и равномерного распределений примесей. Для нормального распределения наблюдается более высокая точка Кюри и больший магнитный момент по сравнению с равномерным порядком. Такие результаты можно наблюдать в экспериментальных работах на магнитных минералах [4-6], где при окислении образцов в процессе нагрева происходит увеличение температуры Кюри и увеличение удельной намагниченности. Окисление осуществляется путем диффузии и только диффузии катионов в неподвижном анионном каркасе [1], анионы кислорода в образованной ими плотнейшей упаковке диффундировать не могут. Как следствие, в процессе таких превращений должен существовать градиент степени окисления и нарушение стехиометрии. Чаще всего окисление железосодержащих природных магнетиков приводит к образованию магнетитовой фазы, и можно говорить о неоднородности состава минерала после окисления. То есть применительно к рисунку 4 кривая равномерного размещения примесей - есть кривая нагрева минерала, а кривая нормального распределения - кривая охлаждения после окисления.

Рис. 4. Результаты моделирования дляx = 0.35 решетки L=30 при нормальном и равномерном распределениях примесных включений.

Помимо окислительных процессов возможно и восстановление, также с выделением магнетитовой фазы [2], что тоже приведет к неоднородности и изменениям свойств минерала, описание которых возможно с помощью приведенной модели.

Заключение

В настоящей работе рассмотрены особенности влияния фазовой неоднородности минералов на магнитные свойства на примере построенной модели. Приведенные теоретические исследования могут найти применения для объяснения процессов термонамагничивания в экспериментах на природных и синтезированных ферримагнетиках. Представленная модель носит описательный характер явлений и не может иметь применимость для оценочных расчетов.

Рецензенты

• Мустафина Светлана Анатольевна, доктор физико-математических наук, профессор, заведующая кафедрой математического моделирования Стерлитамакского филиала БашГУ, г. Стерлитамак.
• Гималтдинов Ильяс Кадирович, доктор физико-математических наук, профессор, заведующий кафедрой прикладной математики и механики Стерлитамакского филиала БашГУ, г. Стелитамак.

Библиографическая ссылка