CC BY
12
УДК 519.6+530.1
DOI: 10.14529/ mmp170413
ДИНАМИКА ВЗАИМОДЕЙСТВИЯ ЕЛОХОВСКИХ ДОМЕННЫХ ГРАНИЦ В ДВУМЕРНОЙ НЕЛИНЕЙНОЙ СИГМА-МОДЕЛИ
Ф.Ш. Шокиров
Физико-технический институт им. С.У. Умарова АН РТ, г. Душанбе, Республика Таджикистан
Проведено численное моделирование процессов взаимодействия 180-градусных доменных границ блоховского типа в фазовом пространстве (2+1)-мерной суперсимметричной 0(3) нелинейной сигма-модели. Метод проведения численных расчетов основан на специальном применении свойств стереографической проекции, где проецированием изосферы на комплексную плоскость устраняется проблема бесконечно больших величин, возникающих в обычной проекции. Таким образом осуществляется необходимая для численного подхода параметризация исследуемой модели в комплексном виде, в которой преодолена сингулярность, возникающая на полюсах изосферы. Использована трехслойная разностная схема второго порядка точности по времени и по координате на пятиточечном шаблоне с весами явного типа. Предложен комплексный программный модуль, реализующий алгоритм численного расчета пространственно-временных топологических структур в трехмерных решетках. Получены модели лобовых столкновений, где в зависимости от динамических параметров наблюдаются процессы образования связанных (бионных) состояний доменных границ, дальнодей-ствующих моделей, прохождения доменных границ магнитных доменов друг сквозь друга, а также формирования радиально-симметричных осциллирующих солитонов.
Ключевые слова: численное моделирование; взаимодействие доменных границ; нелинейная сигма-модель; изотопическое пространство.
Введение
Исследование динамики доменных границ (ДГ) м^гнитоупорядочбнных веществ является актуальной задачей, как с прикладной точки зрения в разработке высокотехнологичных устройств, так и для развития фундаментальных исследований теории магнетизма [1-3]. Теоретические исследования динамики взаимодействующих ДГ в магнитных наноструктурах весьма затруднительны и в данном случае, наиболее эффективным инструментом является численное моделирование. При определенных толщинах тонких пленок энергетически выгодными являются ДГ нееловского типа, где при переходе между доменами вектор намагниченности остается параллельным плоскости поверхности кристалла [3, 4]. Численное исследование динамики взаимодействия нееловских ДГ проведено, в частности, в работе [4] в рамках (2 1)-,мерной анизотропной 0(3) нелинейной сигма-модели (НСМ), где получены дальнодействую-щие ДГ, а также (бионные) модели связанных состояний ДГ.
В настоящей работе приведены результаты численных исследований процессов взаимодействия ДГ блоховского типа, являющиеся энергетически наиболее выгодными в безграничных средах. Плотность функции Лагранжа и гамильтониан исследуемой сигма-модели в стандартной (изоспиновой) параметризации можно записать в следующем виде [5]:
L = g[dßsa3»Sa + (s\ - 1)], (1)
H = g[(dosa)2 + (дкsa)2 + (1 - s¡)], (2)
g =1/2; a = 1, 2, 3; k = 1, 2; ц = 0,1, 2; SaSa - 1 = 0.
Заметим, что 0(3) HCM (n-поле, АЗ-поле) в специальной параметризации (29)
2д^9 + sin(29)(1 - дфдУ) = 0, (3)
2 cos 9д^фд^ф + sin 9д^ф = 0
меридианного сечения ф(х,у,Ь) = 0 (рис. 1 а) изотопического пространства блохов-ской сферы
S2 = SU (2)/U (1) = S0(3)/SO(2) (4)
сводится [4, 5] к вполне интегрируемому уравнению синус-Гордона (УСГ) следующего
ВИДЯ!
2Ш + sin(29) = 0, (5)
где □ - оператор Даламбера. В выражениях (3) и (5) переменные 9(x,y,t) = 0 и ф(х,у,Ь) = 0 описывают эйлеровы углы, связанные с изоспиновыми параметрами НСМ (единичным изовектором S(si, s2, s3) АЗ-поля) следующим образом:
S(si,s2,s3) = (sin9 cos ф, sin9 sin ф, cos 9), (6.1)
* = x + iy = s1+s2 = e* tg2. (6.2)
1 ± s3 2
Для численной схемы применен алгоритм, разработанный в работе [6], где использованы свойства стереографической проекции: точки верхней полусферы (s3 > 0) проецируются на касательную комплексную плоскость, проходящую через «северный полюо точки нижней полусферы (s3 < 0) проецируются на касательную плоскость, проходящую через «южный полюо блоховской сферы S2 (4). В точках «экватора> (s3 = 0) специальным образом производится «прошивка> решения, и, таким образом, осуществляется взаимно-однозначная проекция (компактификация S2 — Щотр) всех точек комплексной плоскости z(x,y) (6.2) та сферу S2 (см. также [5]). Использована трехслойная разностная схема второго порядка точности 0(т2 + h2) на пятиточечном шаблоне с весами явного типа [7]. Аппроксимация проведена на прямоугольной сетке (х,у): hx(3e+3), hy(2e+3) трехмерной решетки (hxy, tT)\ Tmax(5e+4). Устойчивость разностной схемы удовлетворяет требованиям для гиперболических систем уравнений т < min(h/|A|max) [8].
В качестве модельных функций использованы известные топологические решения УСГ (5) в виде ДГ
(ш ш \ (ш ш \ 9 вЛ—х — — хо I +B2¡ ~ у — — уо tg 2 = e \ki ki ) \h¿ k2 J, (7)
n 3n
ф(х,у,Ц = е, е = 2' Y.
На первом этапе в качестве тестовых задач получены модели стационарных и движущихся ДГ (7) блоховского типа УСГ (5) с одномерным распределением намаг-НИЧ6ННОСТИ. Для полученных решений УСГ (5), с учетом (3) и (5), добавлением специально подобранных возмущений ф^0) = е — шт (рис. 1) [4, 5] построены аналогичные
решения для НСМ (1). На основе полученных динамических моделей во втором и третьем пунктах рассмотрены модели двухсолитонных взаимодействий. В частности, получены модели бионных и дальнодействующих взаимодействий ДГ, а также хорошо локализованных, осциллирующих (бризерных) солитонов. В заключительной части обсуждены основные результаты проведенных исследований и подведены итоги.
Рис. 1. Динамика S(si, s2, S3) в изопространстве сферы S2 (4): а) ^(0) = 0; Ь) ^(0) > 0
1. Стационарные и движущиеся доменные границы
В этой части работы приведем результаты численного моделирования стационарных и движущихся ДГ вида (7) блоховского типа с одномерным распределением намагниченности для УСГ (5) и НСМ (1). Соответствующим подбором параметров топологические решения вида (7) можно рассмотреть в более удобной для численного моделирования форме
в(х,уо = 2 tg~1 (е/(ж'*}), (8)
ф(х,уо,г) = £ - ит,
где угловые параметры 9(х, ^ и ф(х, ¿) зависят только от первой пространственной переменной. Заметим, что решения (8) при £ = 0,^н£ = |, 3г описывают соответственно динамику нееловских [4] и блоховских ДГ.
На рис. 2 а приведены результаты численного моделирования эволюции стационарного решения (8) в рамках УСГ (5) для ф = 2 (и = 0). В данном случае изоспи-новая структура поля взаимодействующих ДГ описывается изовектором
S(si, S2, S3) = (1 + е )"1(0, 2ex, 1 - e2x).
(9)
Полученная модель является абсолютно устойчивой: 5Еи ~ 10_6 (Ей, рис. 2 а). Для тестовых моделей, описанных на рис. 2 а, добавлением вращения (и = 0, 5; рис. 1 Ь) вектору АЗ-поля (6.1) изотопического пространства сферы S2 (4) в следующем виде:
т
S(sb s2, s3) = (1 + е ) 1 (2ex si^, 2ex co^, 1 - e2x),
2
2
(10)
получены устойчивые (5Еи е (10"4; 10"3)) ДГ в рамках 0(3) НСМ (1) (рис. 2 Ь).
Заметим, что НСМ обладают свойством формальной лоренц-инвариантности, которое позволяет провести исследование динамики движущихся решений. На рис. 2 с, с!
приведены результаты применения преобразований Лоренца для решений (8) при ш = 0 и ш = 0, 5 соответственно. Интегралы энергии движущихся ДГ сохраняются с высокой точностью: 5Еп(ш = 0) « 10_6; 5Еп(ш = 0, 5) е (10"4; 10"3).
На основе полученных стационарных и движущихся ДГ (8) УСГ (5) и НСМ (1) в следующих пунктах построены модели двухсолитонных взаимодействий с вариацией параметров начальных условий в следующем виде:
п 3п
0 <У! (¿с) = -У2(го) < 1, |Ш12 (¿0 )|< 1, £ =2, у •
В отличие от эволюции ДГ УСГ (5), в рамках 0(3) НСМ (1) наблюдается опреде-
Еп
(рис. 2 с!) ДГ. В данном случае, вследствие введенных возмущений (ш = 0) в динамку вектора АЗ-поля (6.1) происходит излучение лишней энергии в виде линейных волн возмущений, которые поглощаются специально установленными по периметру области моделирования граничными условиями [4 6].
Рис. 2. Плотность энергии (ОН) и проекция изовектора 8(й1 ,в2,в3) на плоскость * (ж, у) (6.2) ДГ (8) УСГ (5) и НСМ (1) при ф = п/2 - шт: а) ш = 0,и (¿о) = 0; Ь) ш = 0, 5,у(г0) = 0; с) ш = 0, V(¿0) « 0, 0995; а) ш = 0, 5,у(¿0) « 0, 0995. Общее время моделирования I е [0; 120]
2. Взаимодействие доменных границ уравнения синус- Гордона
В этой части работы рассмотрены модели лобовых столкновений ДГ вида (8) УСГ (5) (и = 0) при £ = 2, 32е- На рис. 3 приведены результаты лобового (кинк-антикинк) столкновения ДГ (8), движущихся со скоростью ■и(¿о) ~ ±0,0995 при £12 = где изоспиновая структура взаимодействующих полей при Ь = 0 имеет следующий вид:
8(^1, 82, 83) = М(ex)(fl(ex),f2(ex),fз(ex)), (11)
М(ех) = (1 + е2(ж+хо))_1(1 + е-2(ж-ж°))-1, ^(ех) = 0,
¡2(ех) = 2(ех+х° + е-х+х°)(1 - е2х°), fз(ex) = (1 - е2х°)2 - (ех+х° + е-х+х°)2.
В этом случае лобовое взаимодействие ДГ приводит к синхронному изменению ориентации их изоспиновой структуры (рис. 3 а): ф12(|) ^ ф12(Щ-). При этом ДГ свободно проходят друг сквозь друга (рис. 3 Ь). Интеграл энергии системы взаимодействующих ДГ сохраняется с высокой точностью: 8Еи(ш = 0) ~ 10-6 (рис. 3 с). Очевидно, что аналогичные процессы наблюдаются для ДГ с ф12(3т), где топологические солитоны при взаимодействии свободно проходят друг сквозь друга с соответствующим изменением ориентации проекций изовектора 8 ф12( 3П) ^ ф12( 2 )■
В следующем этапе были рассмотрены модели взаимодействия ДГ (8) аналогичные экспериментам, описанным на рис. 3, но с взаимно отличающимися ориентациями изоспиновой структуры (рис. 4 а). В этом случае, во всех экспериментах данной конфигурации наблюдается проявление дальнодействующих сил [4, 5]. На рис. 4 приведен один из примеров вышеуказанных экспериментов для случая - ф1(2) ^^ ф2(Зп) в котором ДГ меняют направление движения без явного взаимодействия (рис. 4 а, Ь). Интеграл энергии системы дальнодействующих ДГ сохраняется с высокой точностью: 8Еи(ш = 0) е (10-6; 10-5) (рис. 4 с).
В следующей части настоящей работы рассмотрены аналогичные взаимодействия ДГ, но в рамках 0(3) НСМ (1), где в динамику изоспиновой структуры (6.1) ДГ добавлены и=0
рис. 1 Ь, 2 Ь, с!).
3. Взаимодействие доменных границ в нелинейной сигма-модели
В этом пункте приведены результаты лобовых взаимодействий ДГ (8) при ш = 0. Отметим, что все эксперименты данной серии приводят к двум формам эволюции системы взаимодействующих ДГ - образованию связанного состояния ДГ и формированию локализованных осциллирующих (бризерных) солитонов.
Результаты первой серии экспериментов показали, что при лобовом взаимодействии ДГ (8) обладающих одинаковыми значениями частоты вращения вектора АЗ-поля (ш1 = ш2 = 0), независимо от других параметров числе иной модели (012(ío), vi2(to)) - наблюдается образование СВЯЗАННОГО (бионного) состояния. На рис. 5 приведен пример образования бионного состояния из двух взаимодействующих ДГ (8) при
S(si, S2, S3) = M(ex)(fi(ex) 81п(шт),¡2(ex) со8(шт), f3(ex)), (12)
¿ = 0.0 £ = 52.8 £ = 75.0
ЧИрт" 1-|- =| 1 ■■■Ч11П11.........*1.
II ...................Ч ...............
Рис. 3. Лобовое столкновение ДГ (8) УСГ (5) при у12(£0) - ±0, 0995, ф12 = п/2 (и = 0): а) проекция изовектора 8(5^ 52, 53) на плоскость £(х, у) (6.2); Ь) эволюция плотности энергии (РН); с) интеграл энергии (Еп). Общее время моделирования £ £ [0; 75]
С = 0.0 С = 50.4 Г = 90.0
Рис. 4. Дальнодействующие ДГ (8) УСГ (5) при -и12(£0) - ±0,0995, ф1 = п/2, ф2 = 3п/2 (и = 0): а) проекция изовектора 8(^,в2,в<3) на плоскость £(х,у) (6.2); Ь) эволюция плотности энергии (ОН)] с) интеграл энергии (Еп). Общее время моделирования £ £ [0; 90]
3П
ф12 = у, ^12(£о) - 0, 0995, и; = и = 0, 5,
где /1 (ех) = /2(ех); М(ех) и /^(ех) (г = 2, 3) определены в (11).
На рис. 5 с, с! приведены иллюстрации осциллирующей динамики плотности энергии (ОН) бионного поля в центральной точке (рис. 5 с), а также интеграл энергии (Еп) модели (рис. 5 с1). Как видно из рис. 5 с период (Р) осцилляций связанных ДГ при £ — 100 приближается к постоянной величине (в данном случае: Р — 10). Потеря энергии бионной системы при £ £ [0; 175] составляет порядка 1,1%; интеграл энергии сохраняется с точностью 5Еп(и = —0, 5) £ (10_3; 10"2) (рис. 5 с1).
Качественно новые результаты получены в следующей серии экспериментов, где при взаимодействии ДГ наблюдается формирование осциллирующих (бризерных) со-литонов (см., например, рис. 6). Численные моделирования показали, что указанное свойство системы взаимодействующих ДГ (8) имеет место при одинаковых значениях углового параметра £: е1 (£0) = е2(£0) и взаимно обратных значений частоты и: и1 = — и2
,52, 53) = М(ех)(±/1 (ех) в1п(ит), /2(ех) сов(ит), /3(ех)), (13)
где М(ех) и /г(ех) (г = 1, 2, 3) определены в (12).
t = Л П . Я6.7 69.Я 79.7 86.4 9S.4
х х
t е [0.0,175.0]
С) d)
Рис. 5. Формирование связанного состояния взаимодействующих ДГ (8) 0(3) НСМ (1) при у12(¿0) « ±0, 0995 £12 = 3п/2 ш12 = —0, 5: а) эволюция плотности энергии (ОН); Ь) проекция изовектора 8(й1 ,в2,в3) на плоскость *(ж, у) (6.2); с) эволюция БН(0, 0); с1) интеграл энергии (Еп). Общее время моделирования е [0; 175]
На рис. 6 а приведен пример образования единого бризерного поля при лобовом
ш
ш12 = ±0, 5
лишняя энергия в виде линейных волн возмущений (рис. 6 а, при I > 90), существенная часть (рис. 6 с!) которой поглощается специальными граничными условиями [4 6]. Тем не менее, часть излучаемой энергии, отражаясь от границ области моделирования Ь[3002 х 1001], оказывает существенное влияние на эволюцию сформированного бризера (рис. 6 а, при > 140).
Для максимального устранения влияния излучаемых волн в следующей серии экспериментов площадь моделирования была увеличена вдвое Ь[3002 х 2001]. Но, в данном случае, увеличение площади пространственной сетки привело к неожиданным результатам, где при взаимодействии ДГ происходит формирование двух бризерных полей (см., например, рис. 7). В примере, приведенном на рис. 7 процесс взаимодействия ДГ можно разделить на два этапа формирование связанного (бионного) состояния ДГ при I ~ [0,140] и формирование двух осциллирующих солитонов (бри-зеров) при > 140.
Динамика проекций 8(^1 ,в2 ,в3) на плоско сть * (ж,у) (6.2) вышеописанной модели приведена на рис. 8 а. Как видно из данной иллюстрации, начальная ротационная динамика проекций изовектора 8(^1, в2, в3) при взаимодействии ДГ (ф12 = 2 ± шт) пере-
Рис. 6. Формирование осциллирующего состояния поля при взаимодействии ДГ (8) 0(3) НСМ (1) - у12(£0) — ±0, 0995, е12 = 3п/2, и12 = ±0, 5: а) эволюция плотности энергии (БН); Ь) проекция изовектора 8(й1, в2, в3) на плоскость £(х,у) (6.2); с) эволюция БН(0, 0); с1) интеграл энергии (Еп). Общее время моделирования £ £ [0; 300]
* I = 0.0 , 62.1 64.8 Л 73.8 . 76.5
■10 .15 -10 -5 0 5 10 15 Ю .,5 !в в 5 « И ",0 '',5 -10 V > Ю И И1> -15 -10 -5 о 5 10 15 ',0 -10 -5 0 5 10 15
Рис. 7. Формирование двойного бризера при взаимодействии ДГ (8) 0(3) НСМ (1) эволюция плотности энергии (БН): ^12(£0) — ±0, 0995, е12 = 3п/2, и12 = ±0, 5. Общее время моделирования £ £ [0; 175]
ходит на осциллирующую динамику (рис. 8 а), ^ Е [62,1; 152,1]). На рис. 8 Ь приведена динамика центральной точки Б И(0,у1) одного из сформированных осциллирующих солитонов (динамика БИ(0, у2) абсолютно идентична БИ(0,у1)). Потеря энергии системы взаимодействующих ДГ при £ Е [62,1; 152,1] составила 5Еп(и = ±0, 5) « 9% (рис. 8 с).
Рис. 8. Формирование двойного бризера при взаимодействии ДГ (8) 0(3) НСМ (1) _ ^12(£о) ~ ±0, 0995, е12 = п/2, и12 = ±0, 5: а) проекция изовектора 8(^,52,53) на плоскость £(х, у) (6.2); Ь) эволюцня БИ(0, у1); с) интеграл энергии (Еп). Общее время моделирования £ Е [0; 175]
4. Обсуждение полученных результатов
Результаты группы экспериментов первой части работы показали, что блохов-ские ДГ (8) также устойчивы в рамках 0(3) НСМ (1) в стационарном (рис. 2 Ь) и движущимся (рис. 2 с!) состояниях. По отношению к аналогичным известным решениям УСГ (рис. 2 а, с) при эволюции блоховских ДГ модели (1) наблюдаются незначительные потери энергии (не более 1, 0% для £ = 120, 0 < ш(Ь0),у(£0) < 1).
На основе численных моделей, полученных в первой части работы, во втором и третьем частях построены модели лобовых взаимодействий ДГ. На первом этапе численная схема и разработанная компьютерная программа была апробирована на тестовых моделях, при моделировании взаимодействий ДГ УСГ (5) (см., например, рис. 3, 4). На основе разработанного метода [6] и предложенного комплекса программ в третьей части работы проведены эксперименты по исследованию взаимодействия ДГ модели (1). В данном случае в динамику изоспиновой структуры ДГ (8) вводится возмущение в виде вращения изовектора 8(й1 ,в2,в3) в пространстве блоховской сферы в2 (4) (см., например, рис. 5 8). Результаты, полученные в предыдущих экспериментах с ДГ (8) УСГ (5) можно разделить на три группы ДГ проходящие друг сквозь друга (рис. 3), дальнодействующие ДГ (4) и связанные состояния дальнодей-ствующих ДГ. Эксперименты четвертой части работы, проведенные в рамках 0(3) НСМ (1) выявили следующие свойства взаимодействующих ДГ (8): а) е1 = е2, и1 = ш2 - формирование связанного состояния ДГ;
b) £1 = е2, и1 = —и2 - формирование связанного состояния ДГ с переходом в радиально-симметричные осциллирующие (бризерные) состояния;
c) £1 = £2, и1 = —и2 - проявление дальнодействующих сил;
с!) прохождение ДГ друг сквозь друга не наблюдается.
Во втором случае (Ь) количество сформированных бризерных состояний, в некоторой степени зависит от размерности области моделирования Ь(х,у) (рис. 6, 7).
Процессы формирования осциллирующего солитона вследствие взаимодействия ДГ в рамках (1+1)-мерной 0(3) НСМ были определены также в работе [9]. В численных экспериментах указанной работы, при взаимодействии одномерных ДГ, обладающих дополнительной динамикой вращения вектора АЗ-поля, прохождения солитонов друг сквозь друга также не наблюдались. В Н аС Т0ЯЩ6 и работе осциллирующие соли-тоны, сформированные при взаимодействии ДГ обладают, в частности, характерной динамикой изоспиновой структуры. Заметим, что аналогичной динамикой изовекто-ра 85,82,833) обладают бризерные решения 0(3) НСМ, Н&ИД6ННЫ6 в работе [10].
Процессы формирования осциллирующего солитона при взаимодействии доменных границ, обнаруженные В Н аС ТОЯТД6 и работе численными моделированиями, были теоретически рассмотрены, в том числе в работах [11, 12].
Заключение
Реальную динамику солитонных решений, где проявляются их особые частицепо-добные свойства можно выявить исследованием процессов их взаимодействия. Статические свойства солитонных решений в большинстве случаев поддаются аналитическому исследованию, но процессы их взаимодействия удается исследовать лишь численными методами [13]. В настоящей работе построением численных эволюционных моделей фронтальных столкновений ДГ осуществлена попытка исследования динамики их взаимодействия в рамках суперсимметричной модели (1). Все солитон-ные решения 0(3) НСМ (1), использованные при построении моделей, были получены численными методами, добавлением специально подобранных возмущений к известным решениям (2+1)-мерного УСГ (5) [4, 5].
Таким образом, в настоящей работе получены численные модели, описывающие эволюцию процессов фронтальных столкновений ДГ магнитных доменов в рамках (2 1)-,мерной 0(3) НСМ. В качестве тестовых задач для апробации алгоритма, численной схемы и компьютерных кодов получены модели стационарных, движущихся и взаимодействующих доменных границ (2+1)-мерного УСГ. Основные результаты настоящей работы приведены в третей части, где описываются модели взаимодействия ДГ, проведенные в рамках 0(3) НСМ. В частности, получены бионные модели взаимодействующих ДГ, модели дальнодействующих ДГ, а также модели, описывающие процессы формирования радиально-симметричных бризерных солитонов. Предложен программный модуль, позволяющий провести комплексный анализ эволюции взаимодействующих решений 0(3) НСМ теории поля, с учетом ее групповых особенностей в двумерном псевдоевклидовом пространстве.
Автор выражает благодарность профессору Х.Х. Муминову за полезные обсуждения и ценные замечания.
Литература
1. Барьяхтар, В.Г. Динамика доменных границ в слабых ферромагнетиках / В.Г. Барьях-тар, В.А. Иванов, М.В. Четкин // Успехи физических наук. - 1985. - Т. 146, № 3. -С. 417-458.
2. Волков, В.В. Динамика доменных стенок в ферромагнетиках / В.В. Волков, В.А. Боков // Физика твердого тела. - 2008. - Т. 50, № 2. - С. 193-221.
3. Филиппов, Б.Н. Статические свойства и нелинейная динамика доменных границ с вих-реподобной внутренней структурой в магнитных пленках (Обзор) / Б.Н. Филиппов // Физика низких температур. - 2002. - Т. 28, № 10. - С. 991-1032.
4. Муминов, Х.Х. Динамика взаимодействия доменных границ в (2+1)-мерной 0(3) нелинейной сигма-модели / Х.Х. Муминов, Ф.Ш. Шокиров // Известия АН РТ. - 2015. -Т. 161, № 4. - С. 57-64.
5. Муминов, Х.Х. Математическое моделирование нелинейных динамических систем квантовой теории поля / Х.Х. Муминов, Ф.Ш. Шокиров. - Новосибирск: Изд-во СО РАН, 2017. - 375 с.
6. Муминов, Х.Х. Многомерные динамические топологические солитоны в нелинейной анизотропной сигма-модели / Х.Х. Муминов // Доклады академии наук РТ. - 2002. - Т. 45, № 10. - С. 28-36.
7. Самарский, A.A. Теория разностных схем / A.A. Самарский. - M.: Наука, 1977.
8. Тихонов, А.Н. Уравнения математической физики / А.Н. Тихонов, A.A. Самарский. -M.: Мир, 1981.
9. Муминов, Х.Х. Взаимодействие динамических и топологических солитонов в 1D нелинейной сигма-модели / Х.Х. Муминов, Ф.Ш. Шокиров // Доклады академии наук РТ. -2016. - Т. 59, № 3-4. - С. 120-126.
10. Шокиров, Ф.Ш. Математическое моделирование бризеров двумерной 0(3) нелинейной сигма-модели / Ф.Ш. Шокиров // Математическое моделирование и численные методы. - 2016. - Т. 12, № 4. - С. 3-16.
11. Белова, Т.И. Солитоны и их взаимодействия в классической теории поля / Т.И. Белова, А.Е. Кудрявцев // Успехи физических наук. - 1997. - Т. 167, № 4. - С. 377-406.
12. Gervais, J.L. Perturbation Expansion Around Extended-Particle States in Quantum Field Theory / J.L. Gervais, A. Jevicki, B. Sakita // Physical Review D. - 1975. - V. 12, № 4. -P. 1038-1051.
13. Маханьков, В.Г. Солитоны и численный эксперимент / В.Г. Маханьков // Физика элементарных частиц и атомного ядра. - 1983. - Т. 14, № 1. - С. 123-180.
Фарход Шамсидинович Шокиров, К^НДИДсХТ физико-математических наук, ведущий научный сотрудник, отдел «Наноматеримы и нанотехнологии>, Физико-технический институт им. С.У. Умарова АН РТ (г. Душанбе, Республика Таджикистан), farhod0475@gmail.com.
Поступила в редакцию 25 сентября 2011 г.
MSC 65Z05 DOI: 10.14529/mmp170413
DYNAMICS OF INTERACTION OF BLOCH TYPE DOMAIN WALLS IN A TWO-DIMENSIONAL NONLINEAR SIGMA MODEL
F.Sh. Shokirov
S.U. Umarov Physical-Technical Institute of AS ET, Dushanbe, Tajikistan
E-mail: farhod0475@gmail.com
Numerical simulation of the interaction of 180-degree Bloch-type domain walls in the phase space of the (2+l)-dimensional supersymmetric 0(3) nonlinear sigma model is carried out. The method of numerical calculations is based on the special application of the properties of stereographic projection, where the projection of the isosphere onto the complex plane eliminates the problem of infinitely large quantities arising in the ordinary projection. Thus, the parametrization of the model under study in a complex form, necessary for the numerical approach, is realized, in which the singularity arising at the poles of the isosphere is overcome. A three-layer explicit difference scheme of the second order of accuracy with respect to time and coordinate on a five-point template is used. A complex programme module is proposed that implements the algorithm for the numerical calculation of space-time topological structures in three-dimensional lattices. The models of frontal collisions are obtained, where, depending on the dynamic parameters, processes of formation of bound (bion) states of domain walls, long-range models, passage of domain walls of magnetic domains through each other, as well as the formations of radially symmetric breathers are observed.
Keywords: numerical simulation; interaction of domain walls; nonlinear sigma model; isotopic space.
References
1. Bar'yakhtar V.G., Ivanov B.A., Chetkin M.V. Dynamics of Domain Walls in Weak Ferromagnets. Physics-Uspekhi, 1985, vol. 28, no. 7, pp. 563-588. DOI: 10.1070/PU1985v028n07ABEH003871.
2. Volkov V.V., Bokov V.A. Domain Wall Dynamics in Ferromagnets. Physics of the Solid State, 2008, vol. 50, no. 2, pp. 199-228. DOI: 10.1134/S1063783408020017.
3. Filippov B.N. Static Properties and Nonlinear Dynamics of Domain Walls with a Vortexlike Internal Structure in Magnetic Films (Review). Low Temperature Physics, 2002, vol. 28, no. 10, pp. 707-738. DOI: 10.1063/1.1521291.
4. Muminov Kh.Kh., Shokirov F.Sh. [Dynamics of Interaction of Domain Walls in (2+1)-Dimensional Non-Linear Sigma Model]. Izvestiya AN RT [News of the Academy of Sciences of the Republic of Tajikistan], 2015, vol. 161, no. 4, pp. 57-64. (in Russian)
5. Muminov Kh.Kh., Shokirov F.Sh. Matematicheskoye modelirovaniye nelineynykh dinamicheskikh sistem kvantovoy teorii polya [Mathematical Modeling of Nonlinear Dynamical Systems of Quantum Field Theory]. Novosibirsk: Publishing House SB RAS, 2017. (in Russian)
6. Muminov Kh.Kh. [Multidimensional Dynamic Topological Solitons in a Nonlinear Anisotropic Sigma Model]. Doklady Akademii nauk Respubliki Tadzhihistán [Reports of the Academy of Sciences of the Republic of Tajikistan], 2002, vol. 45, no. 10, pp. 28-36. (in Russian)
7. Samarskiy A.A. The Theory of Difference Schemes. N.Y., Marcel Dekker, 2001.
8. Tikhonov A.N., Samarskiy A.A. Equations of Mathematical Physics. N.Y., Dover Publications, 2011.
9. Muminov Kh.Kh., Shokirov F.Sh. [Interactions of Dynamical and Topological Solitons in ID Nonlinear Sigma Model]. Doklady Akademii nauk Respubliki Tadzhikistan [Reports of the Academy of Sciences of the Republic of Tajikistan], 2016, vol. 59, no. 3-4, pp. 120-126. (in Russian)
10. Shokirov F.Sh. [Mathematical Modeling of Breathers of Two-Dimensional 0(3) Nonlinear Sigma Model]. Mathematical Modeling and Computational Methods, 2016, vol. 12, no. 4. pp. 3-16. (in Russian)
11. Belova T.I., Kudryavtsev A.E. Solitons and Their Interactions in Classical Field Theory. Physics-Uspekhi, 1997, vol. 40, no. 4, pp. 359-386. DOI: 10.1070/PU1997v040n04ABEH000227.
12. Gervais J.L, Jevicki A., Sakita B. Perturbation Expansion Around Extended-Particle States in Quantum Field Theory. Physical Review D, 1975, vol. 12, no 4, pp. 1038-1051.
13. Makhankov V.G. [Solitons and Numerical Experiments]. Fizika elementamykh chastits i atomnogo yadra [Physics of Elementary Particles and Atomic Nuclei], 1983, vol. 14, no. 1, pp. 123-180. (in Russian)
Received September 25, 2017
