ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«БЕЛГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ
ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ»
(Н И У
«Б е л Г У»)
ИНСТИТУТ ИНЖЕНЕРНЫХ ТЕХНОЛОГИЙ И ЕСТЕСТВЕННЫХ НАУК
КАФЕДРА ТЕОРЕТИЧЕСКОЙ И МАТЕМАТИЧЕСКОЙ ФИЗИКИ
СФЕРИЧЕСКИ СИММЕТРИЧНЫЕ ЭВОЛЮЦИОННЫЕ УРАВНЕНИЯ
ДИВЕРГЕНТНОГО ТИПА ДЛЯ ПСЕВДОВЕКТОРНОГО
СОЛЕНОИДАЛЬНОГО УНИМОДАЛЬНОГО ПОЛЯ
Магистерская диссертация
обучающейся по направлению подготовки 03.04.02. Физика,
Программа физика конденсированного состояния
очной формы обучения,
группы 07001637
Понамарева Алина Эдуардовна
Научный руководитель:
доктор физико-математических
наук, Вирченко Ю.П.
Рецензент:
кандидат физико-математических
наук, доцент Москаленко Н.И.
Белгород 2018
Аннотация
Найдены все возможные эволюционные во времени t пространственно однородные по
x ∈ R3 уравнения для псевдовекторного соленоидального поля M(x, t), x ∈ R3 , имеющие
дивергентный тип, ковариантные относительно вращений пространства R3 и такие, что «генератор» эволюции представляет собой дифференциальный оператор не выше второго порядка.
Весь класс таких уравнений описывается двумя феноменологическими постоянными γ и γ 0 , на
относительную величину которых, накладывается одно условие. Вычислен символ линейного
дифференциального оператора (дисперсионное уравнение), полученного линеаризацией каждого уравнения из описанного в работе класса, около равновесных состояний. Соотношения
между частотой и волновым вектором распространяющихся плоских волн поля M(x, t) – решения дисперсионного уравнения показывают, что все уравнения найденного в работе класса
не обладают диссипацией. Решение задачи основано на перечислении всех возможных алгебраически независимых тензоров второго порядка, которые представляют собой ковариантные
относительно вращений мономы тензорной алгебры с одним порождающим элементом – псевдотензором первого ранга Mj , j = 1, 2, 3. На основе этих мономов найден общий вид потока
плотности поля M(x), который содержит пространственные производные не более чем первого
порядка. Эти дифференциальные выражения плотностей потоков поля определяют линейно
независимый набор термодинамических сил. Из всего набора допустимых термодинамических
сил выбрано 7 таких, которые гарантируют наличие инварианта M2 (x, t) = M 2 движения.
Наконец, на заключительном этапе, из этих семи термодинамических сил отобраны две таких,
которые сохраняют соленоидальность поля M2 (x, t) в процессе эволюции.
Ключевые слова
алгебраическая независимость
дифференциальный оператор
инвариант
линейная независимость
плотность потока поля
плотность магнитного момента
псевдовектор
псевдоскаляр
псевдотензор
символ Леви-Чивита
соленоидальность
тензор
тензорная алгебра
уравнение Ландау-Лифшица
унимодальность
эволюционное уравнение
Оглавление
Список обозначений
Предметный список обозначений
Введение
Глава 1. Описание класса эволюционных уравнений
псевдовекторного поля
4
5
6
11
1.1. Уравнение Ландау-Лифшица
1.2. Стационарные решения уравнения Ландау-Лифшица
1.3. Проблема диссипации
11
14
17
Глава 2. Класс эволюционных уравнений псевдовекторного поля
2.1. Постановка задачи об описании класса плотностей
потоков псевдовекторного поля
2.2. Построение множества линейно независимых тензоров Tijkl
2.3. Общий вид эволюционного уравнения дивергентного типа
для псевдовекторного поля
21
21
Глава 3. Эволюционные уравнения псевдовекторного
унимодального поля
3.1. Построение множества линейно независимых
термодинамических сил
3.2. Эволюционные уравнения с законом сохранения
плотности магнитного момента
Заключение
Литература
24
28
43
43
44
53
54
Список обозначений
В работе мы придерживаемся следующих правил при употреблении
шрифтов для обозначения математических объектов и операций над ними.
• Для обозначения математических операторов (функционалов), для которых в математике имеются устоявшиеся аббревиатуры на основе букв
латинского алфавита, мы употребляем шрифт «roman» – A, B, C, ...; a, b,
c, ... . Например, Re и Im – реальная и мнимая части комплексного числа.
Если таковых устоявшихся аббревиатур не имеется, то мы используем для
обозначения математических объектов различные шрифты, в зависимости
от природы объекта, перечисленные ниже.
• Для обозначения стандартных математических структур используется
ажурный шрифт – A, B, C, ..., например, R – множество действительных
чисел, Z – множество целых чисел, N – множество натуральных чисел.
• Операторы, отображения, функционалы обозначаются прописными
буквами шрифта «sanserif» – A, B, C, ....
• Для обозначения числовых величин (параметров, функций и их аргументов) используются буквы латинского в шрифте «italic» – a, b, c, ... и
греческого алфавитов. При этом латинские буквы i, j, k, l – обозначают целые числа.
• Для обозначения векторов жирные буквы латинского алфавита. Их
компоненты нумеруются индексами i, j, k, l, m, n. При этом принимается
тензорное соглашение о суммировании по повторяющимся парным индексам.
• Для обозначения множеств различных математических объектов используется готический шрифт A, B, C, , ....
• Для обозначения классов множеств различных математических объектов используется шрифт «calligraphic» – A, B, C, , ....
4
Предметный указатель
δij – символ Кронекера
εijk – символ Леви-Чивита
M, либо Mj – псевдовекторное поле плотности магнитного момента
H, либо Hj – вектор внешнего магнитного поля
x – радиус-вектор пространственной точки в R3
t – временной параметр
F, либо Fj – векторное поле термодинамических сил
Sij – псевдотензорное поле плотности потока магнитного момента
Tijkl – тензор коэффициентов плотности потока магнитного момента
q – волновой вектор плоской волны
A – вектор-амплитуда плоской волны
n – единичный вектор
5
Введение
Фундаментальной проблемой механики сплошных сред является формулировка для каждого типа среды адекватного теоретически обоснованного
эволюционного уравнения (см., например, [1]). Эта проблема тесно связана
с основной проблемой неравновесной термодинамики – описанием имеющихся внутри каждой физической среды процессов передачи тепла, которые возникают вследствие протекающих в ней диссипативных процессов.
Действительно, развитие механики и электродинамики сплошных сред, в
последние десятилетия, оказалось сосредоточенным на проблеме конструирования адекватных эволюционных уравнений для описания динамики
сред со сложной структурой. Решение этой проблемы стало очень востребованным в связи с тем, что развитие физики конденсированного состояния,
естественным образом, привело к изучению сред, претерпевающих многочисленных фазовые превращения, в результате которых среды изменяют
свои качественные свойства. При этом для характеризации локального термодинамического состояния каждой такой среды приходится использовать
дополнительные интенсивные термодинамические параметры, которые называются параметрами порядка изучаемой среды и которые связаны, как
правило, со спонтанным нарушением какой-либо симметрии. Кроме того,
естественным следствием указанной линии развития механики и электродинамики сплошных сред, а также следствием изучения термодинамики
сред, характеризуемых, наряду с традиционными интенсивными термодинамическими параметрами такими как плотность, температура, концентрации различных химических компонентов и т.д., различными дополнительными параметрами порядка различной математической природы явилось
то, что появилась необходимость построения неравновесной термодинамики таких сред.
Одним из путей решения указанной проблемы является так называемый
гамильтонов подход построения эволюционных уравнений сред со спонтанно нарушенной симметрией (см., например, [2, [3]). Этот подход оказался
очень плодотворным при решении многих задач из указанного выше круга, а, с другой стороны, оказалось, что он обладает рядом недостатков и,
по видимому, его основные положения являются слишком стеснительны6
ми при построении таких динамических уравнений, предсказания которых
полностью согласовывались бы с экспериментальными данными. По этой
причине, оказалось желательным создание альтернативного подхода при
построении эволюционных уравнений для твердотельных сред с нарушенной симметрией. Причем такой подход должен быть более общим, чем указанный гамильтонов формализм для того, чтобы иметь возможность конструировать на его основе эволюционные уравнения для описания всех известных в настоящее время типов конденсированных сред. Ясно, что такую
наиболее общую теоретическую схему такого подхода можно создать положив в ее основу только наиболее общие физические теоретические принципы, нерушимость которых не вызывает сомнений в любых мыслимых в
настоящее время физических ситуациях.
По нашему мнению, решение общей проблемы конструирования эволюционных уравнений должно состоять из нижеследующих положений.
Пусть макроскопическое термодинамически равновесное состояние изучаемой сплошной среды, с физической точки зрения, вполне характеризуется фиксированным набором интенсивных термодинамических параметров {ξl ; l = 1 ÷ N }. Согласно определению интенсивных термодинамических параметров, они являются плотностями соответствующих экстенсивных термодинамических функций изменяющихся пропорционально объему
области пространства, занимаемой средой. Полнота набора {ξl ; l = 1 ÷ N }
означает, что любая измеряемая физическая характеристика среды может
быть определена на основе известных значений термодинамических параметров из указанного списка. Более того, подразумевается, что этот полный
набор термодинамических характеристик состояния среды является минимальным, то есть удаление из него какого-либо из параметров приведет к
тому, что получаемый таким образом более узкий список параметров уже
не описывает состояние среды однозначным образом.
Одним из основных положений механики сплошных сред является предположение о том, что локально в каждой пространственной точке, характеризуемой радиус-вектором x, состояние среды вполне описывается набором
полей {ξl (x); l = 1 ÷ N } – функций, зависящих от пространственной точки
x. Значения этих полей в каждой точке x представляют наборы значений
плотностей {ξl ; l = 1÷N } в том случае, если бы весь объем среды находился в том же термодинамическом состоянии, что и ее объем расположенный
в физически малой окрестности точки x.
В пространственно неравновесном состоянии значения полей {ξl (x); l =
7
1 ÷ N } в каждой точке x изменяются со временем, то есть в рамках механики сплошных сред необходимо рассматривать эти поля зависящими от
времени t.
Предполагается, что список полей ξl (x, t), l = 1 ÷ N допускает разбиение на группы таким образом, что набор полей в каждой из них преобразуется посредством тензорных (спин-тензорных) представлений группы
вращений O3 евклидового пространства. Это положение является тезисом
А.Эйнштейна о тензорном характере фундаментальных физических законов.
Таким образом, объектом изучения механики фиксированной сплошной
среды являются наборы зависящих от времени полей {ξl (x, t); l = 1 ÷ N },
x ∈ Ω ⊂ R3 плотностей полного набора термодинамических характеристик
среды, которые представляются наборами тензорных плотностей. Тогда
возникает проблема формулировки тех эволюционных дифференциальных
уравнений
ξ˙m (x, t) = L[ξ] (x, t) , m = 1 ÷ N ,
(0.1)
m
определяемых операторами Lm [·], m = 1 ÷ N – «генераторами» эволюции, где здесь и далее ξ = hξ1 , ..., ξN i. Решения такой системы уравнений
должны описывать динамику рассматриваемой пространственно распределенной термодинамической системы в представлении ее в виде сплошной
среды.
Можно сформулировать набор естественных с общетеоретической точки
зрения требований, предъявляемых к конструированию этих уравнений,
как «фундаментальных уравнений» механики сплошных сред. Таковыми
являются следующие.
1. Операторы Lm [·], m = 1 ÷ N не должны явно зависеть от времени t,
что является следствием однородности времени – независимости динамики
системы от выбора начала отсчета времени.
2. Операторы Lm [·], m = 1 ÷ N не должны явно зависеть от от пространственной точки x, что является следствием пространственной однородности (трансляционной инвариантности) – независимости динамики системы
от выбора начала отсчета пространственных координат.
3. Совокупность эволюционных уравнений (1) преобразуется ковариантным образом при преобразованиях вращений группы O3 , что является
следствием изотропии пространства – независимости динамики системы
от выбора направления координатных осей для ее описания.
8
К этой совокупностей требований, следующих из общих принципов теоретической физики, добавим еще требование локальности эволюционных
уравнений. Математическая формулировка этого требования состоит в следующем.
4. Операторы Ll [ξ], l = 1÷N (вообще говоря, нелинейные), действующие
на ξ, являются дифференциальными и, более того, мы будем считать, что
они имеют не выше второго.
Физическим источником этого требования является то, что динамика
сплошных сред основана на описании эволюции физических полей, обладающих малыми градиентами, а каждая частная производная по пространственным переменным, неявным образом, пропорциональна с малым параметром – отношением пространственного масштаба r0 порядка линейных
размеров области, физическое состояние которой уже с большой точностью
можно считать термодинамическим, к пространственному масштабу L, характеризующему неоднородности в системе. Здесь r0 – линейный размер
пространственной области, начиная с которого можно говорить, что в ее
объеме содержится достаточно много частиц так, что для состояния этой
системы частиц уже применимо локальное термодинамическое описание.
В естественных окружающих нас физических средах в земных условиях r0
имеет порядок ∼ 10−6 см. В то же время масштаб L ∼ 10−4 ÷ 1 см и более.
В связи с требованием 4, предъявляемым к эволюционным операторам
Ll [ξ], l = 1 ÷ N , заметим, что их всегда можно представить в следующем
каноническом виде
(L[ξ])l = ∇m S[ξ]
+ fl (ξ) ,
(0.2)
lm
где fl , l = 1 ÷ N – функции от переменных набора ξ, которые называем
источниками, и каждый дифференциальный оператор Slm [ξ] при фиксированном значении l = 1 ÷ N представляет вектор c компонентами, нумеруемыми индексом m = 1, 2, 3. Каждый такой оператор представляет
плотность потока физического поля ξl (x, t) в пространственно временной
точке hx, ti.
В частном случае, если при некотором значении l источник fl (ξ) = 0,
то уравнение (2) с соответствующим номером является с математической
точки зрения уравнением дивергентного типа. В этом случае мы говорим,
что это уравнение представляет собой локальный закон сохранения.
Перечисленные положения представляют собой совокупность довольно
общих положений, сформулированных на основе общих физических прин9
ципов, которые имеют универсальных характер для всех физических сред.
Однако, они определяют только лишь общую математическую структуру
эволюционных уравнений. Они совершенно недостаточны для формулировки окончательного вида совокупности динамических уравнений, описывающих эволюцию состояний фиксированной изучаемой среды. Для получения
окончательного вида системы уравнений требуются еще какие-то физические требования, связанные с соотношением между величинами тех или
иных физических параметров – локальных характеристик среды. Однако,
по нашему мнению, очень важной задачей является описание классов систем уравнений, характеризуемых определенным типом поведения наборов
полей ξ(x) локальных термодинамических характеристик по отношению к
преобразованиям вращения пространства. В настоящей роботе такая задача решается в том случае, когда набор ξ состоит из одного псевдовектора,
и при этом в описании термодинамического состояния среды нет никаких
дополнительных (выделенных) векторных и тензорных величин, которые
не являются пространственно распределенными характеристиками системы. На базе полученного описания класса эволюционных уравнений для
указанного типа полей описаны все типы эволюционных уравнений, обладающие специальными свойствами их решений: унимодальностью и соленоидальностью. Решение этих задач может иметь непосредственное отношение к построению адекватных уравнений ферродинамики сферически
симметричных ферромагнетиков в отсутствие внешнего магнитного поля.
10
Глава 1. Описание класса эволюционных уравнений
псевдовекторного поля
В этой главе мы кратко проанализируем современное состояние основ
ферродинамики и поставим задачу о нахождении допустимого расширения класса уравнений для описания эволюции ферромагнитным образом
упорядоченных сплошных сред.
1.1. Уравнение Ландау-Лифшица
Рассмотрим описание динамики ферромагнитным образом магнитоупорядоченной среды. При этом мы будем пренебрегать изменениями других
ее локальных термодинамических характеристик. Мгновенное локальное
состояние такой среды описывается зависящем от времени полем M(x, t)
плотности магнитного момента (плотностью намагниченности) ферромагнитной среды в пространственной точке, описываемой радиус-вектором x,
в момент времени t. Поле M(x, t) принимает псевдовекторые значения, то
есть оно преобразуется как векторное поле при непрерывных преобразованиях группы O3 поворотов пространства R3 , но, в отличие от векторного
поля, оно не изменяется при отражениях пространства R3 (см., например,
[4], [5]).
В настоящее время является общепринятым, что эволюция магнитной
структуры в магнитоупорядоченной твердотельной среде описывается так
называемым уравнение Ландау-Лифшица (см., например, [6], [7]). В том
случае, когда среда является ферромагнитным образом магнитоупорядоченной, это уравнение записывается в следующем общем виде
δ W [M]
Ṁ = M,
,
(1.1.1)
δ M(x, t)
Здесь и далее точкой мы обозначаем частную производную по времени t. В
правой части уравнения стоит вариационная производная δW[M]/δM(x, t)
функционала W[M] по полю M(x, t) в пространственно-временной точке
(x, t). Этот функционал представляет собой плотность энергии поля намагниченности в среде, описываемой полем M(x, t). Cама производная представляет эффективное «среднее» магнитное поле в пространственной точке
11
x среды. Энергия, которой обладает магнитная структура, связана с имеющимися в ней внутренними напряжениями, возникающими в том случае,
если в среде реализуется неоднородное распределение намагниченности.
В простейшем случае сферически симметричного ферромагнетика в
присутствии внешнего магнитного поля вектором его напряженности
H(x, t), в общем случае зависящей от пространственной точки с радиусвектором x и времени t, функционал W [M] имеет следующий вид
Z h
i
∇j Mk (x, t) ∇j Mk (x, t) − Hj (x, t)Mj (x, t) dx , (1.1.2)
W[M] = γ
R3
γ > 0 – постоянная, которая называется гиромагнитным отношением. В
этой формуле ∇ – векторный оператор дифференцирования (градиент).
Кроме того, здесь и далее, использованы принятые в тензорной алгебре
индексные обозначения векторов и тензоров. По повторяющимся индексам
подразумевается, как это принято, суммирование по их значениям 1, 2, 3
(см. [4], [5]).
Для ферромагнетика с плотностью энергии W [M] вида (2) уравнение
(1) принимает следующий вид
Ṁ = γ[M, ∆M − H] .
(1.1.3)
где ∆ – дифференциальный оператор Лапласа в R3 .
Используя тензорные обозначения уравнение (3) представляется в виде
Ṁj (x, t) = γεjkl Mk (x, t) ∆Ml (x, t) − Hl (x, t) , j = 1, 2, 3.
(1.1.4)
В частности, при Hj = 0, имеем
Ṁj (x, t) = γεjkl Mk (x, t)∆Ml (x, t) ,
j = 1, 2, 3.
(1.1.5)
Здесь Mj (x, t) и Hj (x, t) – декартовы компоненты, соответственно, полей M(x, t) и H(x, t), M(x, t) = hM1 (x, t), M2 (x, t), M3 (x, t)i, H(x, t) =
hH1 (x, t), H2 (x, t), H3 (x, t)i, εjkl – универсальный антисимметричный псевдотензор третьего ранга (символ Леви-Чивита). Его значения остаются
неизменными при поворотах пространства R3 и равными ±1 Кроме того, здесь и далее принимается, что, во всех формулах, свободные («говорящие») нижние индексы i, j, k, l, m, n принимают значения {1, 2, 3}. Тот
факт, что описывается динамика именно сферически симметричного ферромагнетика, отражается в том, что уравнение (3) инвариантно относительно пространственных вращений вектора M в том случае, когда H(x, t) = 0.
12
Обычно, этот факт в теоретической физике выражают утверждением о
том, что в уравнении (3) c H(x, t) = 0 не используются никакие иные векторные (тензорные) физические поля, кроме самого поля M(x, t).
Уравнение (3) обладает тем замечательным свойством, что оно сохраняет во времени унимодальность поля M(x, t) (постоянство M2 (x, t)), в чем
мы немедленно убеждаемся скалярно умножив обе части уравнения (3) на
M(x, t),
Mj (x, t)Ṁj (x, t) = γεjkl Mj (x, t)Mk (x, t) ∆Ml (x, t) − Hl (x, t) = 0 ,
где учтено тензорное тождество εjkl Ml Mk = 0. Таким образом,
d 2
M = 2Mj Ṁj = 0
dt j
и Ml2 (x, t) = M 2 является инвариантом движения.
Свойство уравнения сохранять абсолютную величину намагниченности
в каждом физически малом объеме среды важно с физической точки зрения, так как в реальных (относительно медленных) эволюционных процессах, происходящих в магнитной структуре, практически не происходит
изменения температуры и, как следствие, не происходит изменения зависящей только от локальной температуры T (x, t) величины |M(x, t)|.
Заметим, что уравнение (4) представимо в виде (0.2). В самом деле,
положив N = 3 и выбрав в качестве набора параметров ξ три компоненты
псевдовектора Mj , уравнение (4) представим в виде (0.2)
Ṁj (x, t) = ∇k Sjk (x, t) + fj ,
(1.1.6)
Sjk (x, t) = γεjml Mm (x, t)∇k Ml (x, t) ,
(1.1.7)
где
так как
∇k Sjk (x, t) = γεjml ∇k Mm (x, t)∇k Ml (x, t) =
= γεjml ∇k Mm (x, t) ∇k Ml (x, t) + γεjml Mm (x, t)∆Ml (x, t) ,
где первое слагаемое тождественно равно нулю, как свертка по индексам m, l антисимметричного псевдотензора εjml и симметричного тензора
[∇k Mm (x, t)][∇k Ml (x, t)], и, кроме того,
fj (x, t) = −γεjml Mm (x, t)Hl (x, t) .
13
(1.1.8)
В частности, уравнение (5) имеет дивергентную форму, если внешнее магнитное поле равно нулю. Таким образом, уравнение Ландау-Лифшица в
рассматриваемом случае вкладывается в класс эволюционных уравнений,
описанный пп.1-4 во Введении.
Заметим, что значения тензорного поля Sjk (x, t) в каждой пространственно-временной точке hx, ti с радиус вектором x составляют псевдотензор
второго ранга, то есть он, в отличие от тензора второго ранга, изменяет
знак при отражениях пространства. Это связано с тем, что символ ЛевиЧивита εjkl – алгебраический третьего ранга (см. [7]), который является
пседотензором, то есть он не изменяет знак при преобразованиях отражения R3 . Точно также значения поля Mj (x, t) в каждой пространственновременной точке представляют псевдовектор, то есть они при отражении
пространства R3 не изменяются. Наконец, векторный дифференциальный
оператор ∇j , то есть он преобразуется как вектор при ортогональных преобразованиях пространства R3 , он изменяет знак на обратный при отражении. Тогда, так как поле Sjk (x, t), которое можно назвать плотностью
потока магнитного момента, представляет собой тензорное произведение
с последующей сверткой по одному из индексов тензора второго ранга
γεjml Mm (x, t) и псевдотензора ∇k Ml (x, t), при непрерывных преобразованиях пространства преобразуется как тензор второго ранга, а при отражениях изменяет знак.
1.2. Стационарные решения уравнения Ландау-Лифшица
Проанализируем кратко равновесные унимодальные M2 (x) = M 2 решения уравнения Ландау-Лифшица в отсутствии внешнего магнитного поля. При этом мы ограничимся рассмотрением только одномерных соленоидальных решений с целью сравнения их со стационарными решениями получаемыми в этой работе более общих эволюционных уравнений для
псевдовекторного поля. В общем случае стационарные решения уравнения
Ландау-Лифшица, согласно (5), определяются уравнением
[M, ∆M] = 0 .
(1.2.1)
∆M = λM ,
(1.2.2)
Тогда
где λ = λ[M] – локальный функционал
1)
1)
от M. Рассмотрим одномерные
То есть его значения зависят только от значений поля M и его производных в точке x.
14
решения M(x), которые зависят только от x1 , т.е.
M00 (x1 ) = λ(x1 )M(x1 ) .
(1.2.3)
Условие соленоидальности поля M(x1 ), то есть условие (∇, M) = 0 в этом
случае сводится к M10 = 0, где здесь и далее штрихом обозначена частная
производная по x1 . Таким образом, M1 (x1 )= const. Подставляя такого типа решение в дифференциальное уравнение (3) для компоненты M1 (x1 ),
находим 0 = M100 (x1 ) = λM1 (x1 ), то есть M1 λ = 0. Тогда либо λ = 0, либо
M1 = 0.
В первом случае M00 (x1 ) = 0. Тогда M0 (x1 ) = const. Для соблюдения
условия унимодальности M2 (x1 ) = M 2 необходимо, чтобы эта постоянная
интегрирования была равна нулю. Следовательно, в рассматриваемом случае имеется только один тип равновесных решений M(x1 ) = M, где M
–произвольный постоянный вектор.
Рассмотрим второй случай M1 = 0. Условие унимодальности поля
M(x1 ) требует, чтобы M2 (x1 ) = M 2 . Достаточно рассмотреть случай, когда
M = 1. Это приводит к условию
(M00 (x1 ))2 = λ2 (x1 ) .
(1.2.4)
Таким образом имеем следующую систему уравнений для компонент
M2 (x1 ) и M3 (x1 ),
M200 (x1 ) = ±λ(x1 )M2 (x1 ) ,
M300 (x1 ) = ±λ(x1 )M3 (x1 ) ,
λ(x1 ) = ±|(M200 (x1 ))2 + (M300 (x1 ))2 |1/2
(1.2.5)
(1.2.6)
Ввиду M1 = 0 и условия M = 1, положим теперь
M2 (x1 ) = cos ϕ(x1 ) ,
M3 (x1 ) = sin ϕ(x1 ) .
Продифференцировав по x1 , получаем
M20 (x1 ) = −ϕ0 (x1 ) sin ϕ(x1 ) ,
M3 (x1 ) = ϕ0 (x1 ) cos ϕ(x1 ) .
M200 (x1 ) = −(ϕ0 (x1 ))2 cos ϕ(x1 ) + ϕ00 (x1 ) sin ϕ(x1 ) ,
(1.2.7)
M300 (x1 ) = −(ϕ0 (x1 ))2 sin ϕ(x1 ) + ϕ00 (x1 ) cos ϕ(x1 ) .
(1.2.8)
Тогда
(M200 (x1 ))2 + (M300 (x1 ))2 =
= [(ϕ0 (x1 ))2 cos ϕ(x1 ) + ϕ00 (x1 ) sin ϕ(x1 )]2 +
15
= [(ϕ0 (x1 ))2 sin ϕ(x1 ) − ϕ00 (x1 ) cos ϕ(x1 )]2 =
0
4
00
2
= (ϕ (x1 )) + (ϕ (x1 )) .
(1.2.9)
Подстановка выражений (7) - (9) в систему уравнений (5), (6) дает
−[(ϕ0 (x1 ))2 cos ϕ(x1 ) + ϕ00 (x1 ) sin ϕ(x1 )] =
1/2
0
4
00
2
= ± (ϕ (x1 )) + (ϕ (x1 ))
cos ϕ(x1 ) , (1.2.10)
−[(ϕ0 (x1 ))2 sin ϕ(x1 ) − ϕ00 (x1 ) cos ϕ(x1 )] =
1/2
0
4
00
2
= ± (ϕ (x1 )) + (ϕ (x1 ))
sin ϕ(x1 ) . (1.2.11)
Умножением уравнения (10) на sin ϕ(x1 ) и уравнения (11) – на sin ϕ(x1 ) и
вычитанием одного из другого получаем
ϕ00 (x1 ) = 0 .
(1.2.12)
Тогда из (10) и (11) получаются два уравнения
−(ϕ0 (x1 ))2 cos ϕ(x1 ) = ±(ϕ0 (x1 ))2 cos ϕ(x1 ) ,
−(ϕ0 (x1 ))2 sin ϕ(x1 ) = ±(ϕ0 (x1 ))2 sin ϕ(x1 ) ,
из которых, ввиду того, что sin ϕ и cos ϕ не могут обратиться в нуль, следует,
(ϕ0 (x1 ))2 = ∓(ϕ0 (x1 ))2 .
Здесь знак + в правой части не может реализоваться, а в случае знака −
получается тождество 2) .
Условие (12) приводит к решениям в виде ϕ(x1 ) = kx1 + ϕ0 , k = const
и, соответственно, получаются решения исходной системы уравнений (5) в
виде «поперечной» спирали с осью вдоль x1 ,
M = h0, cos ϕ, sin ϕi ,
ϕ = ϕ0 + kx1 .
(1.2.13)
Следует заметить, что такие равновесные решения реализуются и до
перехода к пределу сплошной среды, в случае «микроскопической» т.н.
векторной модели на кристаллической решетке (см., например, [8],[9]).
2)
Либо ϕ(x1 ) = const, что дает уже полученные нами ранее решения в виде постоянного вектора M.
16
1.3. Проблема диссипации
Эволюционные уравнения вида (1.2) являются основой ферродинамики
и широко используются для описания процессов, происходящих в магнитных структурах, протекающих вследствие имеющихся внутри них неравновесных пространственных распределений плотности магнитного момента,
возникших, например, в результате каких-либо технологических процессов.
Однако, несмотря на признанную в теории магнетизма адекватность динамического уравнения (1.2) и его общеупотребительность при решении различных практических задач ферродинамики, оно, как известно, обладает
плохим, с физической точки зрения, свойством, которое не соответствуют
физическим представлениям. Это уравнение не описывает диссипативных
процессов внутри магнитной структуры, то есть оно не описывает наличие
«магнитного трения», действие которого, в конце концов, приводит магнитоупорядоченную среду к состоянию равновесия. На это обстоятельство
неоднократно обращалось внимание и предпринимались попытки как-то
«исправить» это уравнение, чтобы учесть диссипацию (см., например, [8]).
Наличие магнитного трения устанавливается на уровне микроскопической
квантовой теории магнетизма [7], в рамках которой вычисляется декремент
затухания спиновых волн как мнимая часть массового оператора магнонов.
Так как ясно, что магнитное трение нельзя рассматривать как квантовый
эффект, то оказывается, что между макроскопической теорией эволюции
неравновесных магнитных структур и их микроскопической теорией имеется противоречие, которое желательно преодолеть.
Продемонстрируем указанное противоречие на примере возмущения
стационарного решения уравнения (1.3) –пространственно однородного
распределения поля M(x) = M. С этой целью Рассмотрим возмущение
M(x, t) = M + m(x, t)
(1.3.1)
стационарного решения в виде плоской монохроматической волны
m(x, t) = m exp(i(k, x)−iωt), где m – комплекснозначный вектор. Подстановка этого решения в уравнение (1.3) приводит к следующему уравнению
для векторной амплитуды m:
iωm = γk2 [M, m] ≡ Am ,
(1.3.2)
где матрица A, определяемая правой частью этого уравнения дается следующими равенствами
(Am)j = γk2 εjkl Mk ml ,
17
то есть
(A)jk = γk2 εjkl Mk .
(1.3.3)
Согласно (2), для того чтобы имелось решение указанного типа (1) нужно, чтобы вектор m был собственным для матрицы A. Заметим, что матрица A антисимметрична. Тогда она обладает только чисто мнимыми собственными числами. Поэтому частота ω в (2) Должна быть вещественной,
что означает отсутствие затухания возмущений рассматриваемого типа.
Так как матрица A вещественна, то каждому ее собственному числу λ соответствует собственное число λ∗ . Для матрицы 3-го порядка это означает,
что она обязательно имеет нулевое собственное число.
В самом деле, матрица A, очевидно, обладает нулевым собственным значением с собственным вектором M. Следовательно, det A = 0. Однако, этот
собственный вектор нефизический, так как необходимо учесть условие унимодальности для возмущенного решения (1),
(M + m(x, t), M + m(x, t)) = M 2 .
(1.3.4)
Отсюда, удерживая члены первого порядка по m, находим
(m, M) = 0 .
(1.3.5)
Сужение действия матрицы A в плоскость, ортогональную вектору M, дает
2 × 2-матрицу
(A)jk = −γk2 M εjk3 .
Так как собственные числа 2 × 2-матрицы
0 1
−1 0
равны ±i, собственные частоты возмущения m(x, t) равны ω = ±γk2 M .
Соответствующими собственными векторами в ортобазисе he1 , e2 , e3 i,
e3 k M для этих частот являются e1 ± ie2 .
В связи с описанной ситуацией относительно бездиссипативности уравнения (1.2), его нельзя рассматривать как определяющее, в полной мере,
динамику пространственного распределения плотности намагниченности в
ферромагнитной среде на больших интервалах времени, когда диссипативные процессы начинают играть существенную роль и, вследствие неравновесного распределения плотности магнитного момента в среде, происходит
18
перекачка части энергии взаимодействия между различным образом намагниченных малых пространственных областей среды в тепловую энергию.
Таким образом, возникает вопрос, каким образом нужно изменить уравнение (1.2), в частности, уравнение (1.3), чтобы устранить указанный дефект макроскопической теории. Трудность решения этой задачи состоит в
том, что простая модификация уравнения (1.2), как это делается в других
разделах механики сплошных сред, невозможна. Обычно это осуществляется посредством добавления в правую его часть слагаемого в виде дифференциального оператора второго порядка с отрицательным символом, которое бы описывало, на феноменологическом уровне, процессы взаимного
«трения магнитных моментов» среды, находящихся в различных пространственных точках среды. Такому примитивному преобразованию основного
уравнения ферродинамики (1.2) препятствует важное с физической точки зрения требование, которому должны удовлетворять любые слагаемые,
добавляемые в правую часть уравнения. Любое видоизменение уравнения
(1.2) должно быть такими, чтобы, в течение эволюции системы, сохранялась величина M2 (x, t) = M 2 = const. В настоящей работе предлагается
подход к решению указанной проблемы в рамках общего подхода к составлению макроскопических уравнений, описывающих эволюцию состояния
сплошной среды, которые сформулированы в пп. 1-4 введения.
Заметим также, что уравнение (1.2), в частности (1.3), не сохраняет дивергенцию поля M(x, t). Это обстоятельство приводит к тому, что
не сохраняется соленоидальность поля суммарной магнитной индукции
B(x, t) = H(x, t) + 4πM(x, t) в ферромагнитной среде, то есть невозможно
удовлетворить равенству (∇, B)(x, t) = 0, выполнение которого требуется
системой уравнений Максвелла, так как даже при наличии внешнего магнитного поля H(x, t), которое является соленоидальным (∇, H)(x, t) = 0.
В связи со сделанным замечанием, по-видимому, нужно также потребовать, чтобы искомое «правильное» уравнение для плотности магнитного момента должно обладать, наряду с указанным выше инвариантом
M2 (x, t) = M 2 , также и инвариантом (∇, M)(x, t).
Продемонстрируем указанный факт несохранения соленоидальности поля M(x, t). С этой целью, как и выше, рассмотрим возмущение m(x, t) =
m(t) exp(i(k, x)) постоянного поля M – стационарного решения уравнения
(1.3). Подстановка этой функции в (1.3) приводит к уравнению для ампли-
19
туды m(t)
ṁ = −γk2 [M, m] .
(1.3.6)
Подстановка же этой функции в условие соленоидальности (∇, m(x, t)) = 0
приводит к условию поперечности
(k, m(t)) = 0 .
(1.3.7)
Из этого условия в сочетании с (6) получаем
(k, [, M, m(t)] = 0 .
(1.3.8
Так как, кроме того, ввиду сохранения унимодальности суммарного поля
должно, аналогично (5), быть удовлетворено условие (m(t), M) = 0, то
отсюда и из (8) следует, что совместимость всех условий, налагаемых на
возмущение, возможна только при k k M.
В свете сформулированных нами общих проблем ферродинамики, в следующих главах настоящей работы мы частично ответим на поставленные
выше вопросы. Нами полностью будет описан класс сферически симметричных эволюционных уравнений для псевдовекторого поля M(x, t), которые обладают инвариантом M2 (x, t) = M 2 движения, в том случае, когда
эти уравнения содержат пространственные производные не выше второго
порядка. Далее, из этого класса уравнений выделен класс уравнений, которые, в дополнение к указанному инварианту, обладают также инвариантом
(∇, M)(x, t). Оказалось, что этот подкласс довольно узок он описывается
двумя вещественными параметрами (физическими константами). Показано, однако, все уравнения в этом подклассе не описывают диссипативную
динамику псевдовекторного поля.
20
Глава 2. Класс эволюционных уравнений
псевдовекторного поля
В этой главе мы поставим и решим задачу об описании класса всех
эволюционных уравнений дивергентного типа для псевдовекторного поля,
удовлетворяющих сформулированным во введении положениям 1-4, в предположении об отсутствии у них дополнительных выделенных векторных
(тензорных) характеристик.
2.1. Постановка задачи об описании класса плотностей
потоков псевдовекторного поля
В этом разделе будет сформулирована постановка задачи об описании
класса всех сферически симметричных эволюционных уравнений дивергентного типа для псевдовекторного поля M(x, t). Такая постановка задачи основана на простой идеологической схеме о том, что эволюционные
уравнения, описывающие динамику термодинамического состояния среды,
должны формулироваться в виде в локальных законах сохранения. Именно, на такой схеме, в частности, основан вывод уравнений гидродинамики
(см. [11]), где в основу кладутся локальные законы сохранения плотностей
энергии, импульса и массы вещества.
Поставим задачу об описании класса всех дифференциальных уравнений, которые имеют вид
Ṁj = ∇k Sjk ,
(2.1.1)
в которых плотность потока Sjk (x, t) не зависит явно от x и t, а является локальным функционалом поля M(x, t), то есть Sjk является функцией значений поля Mj (x, t) и его производных не выше первого порядка в
пространственно-временной точке hx, ti. Причем будем требовать, чтобы
зависимость от производных была линейной. Эти ограничения являются
следствием требования 4. Таким образом, мы будем нужно описать класс
всех уравнений (1), в которых плотность потока поля Mj имеет вид
Sjk = Tjkln ∇l Mn ,
(2.1.2)
где Sj,k является псевдотензором второго ранга и Tj,k,l,n является функцией
от значений поля M в hx, ti. Следовательно, необходимо описать класс всех
выражений вида (2), который является линейным многообраз ием.
21
Здесь и далее, так как значения функционалов всех тензорных величин
вычисляются в фиксированной пространственно-временной точке, то, если это не вызывает недоразумений, мы будем опускать пространственно
временные аргументы.
Заметим, что, в общем положении в пространстве дифференцируемых
полей, функции Mj (x, t) и их частные производные ∇l Mn функционально
независимы. В противном случае, если существует какая-то функциональная зависимость f (Mj , ∇k Ml ) = 0, то она представляет собой дифференциальное уравнение в частных производных для поля Mj . Следовательно,
допустимые поля Mj (x, t),которые удовлетворяют этому дифференциальному уравнению. Многообразие решений этого уравнения имеет коразмерность, не меньшую 1, что означает, что такая ситуация невозможна в общем
положении.
Ввиду функциональной независимости поля Mj и его производных, на
основании (2), можно записать
∂Sjk
= Tjkln .
∂(∇l Mn )
(2.1.2)
Левая часть этого равенства представляет собой тензор четвертого ранга.
Следовательно, совокупность коэффициентов Tjkln в представлении (2) образует тензор четвертого ранга. Тогда задача описания класса всех уравнений вида (1) сводится к задаче описания всех функций Tjkln [M], которые
принимают значения в виде тензоров четвертого ранга. Каждый тензор
четвертого ранга представим в виде разложения
Tjkln =
81
X
(a)
f (a) Tjkln .
(2.1.3)
a=1
по базису каких-то стандартных линейно независимых тензоров четвертого
(a)
ранга Tjkln со скалярными коэффициентами (всего имеется 34 = 81 элементов базиса). В тензорной алгебре в трехмерном пространстве имеется два
исключительных универсальных тензорных объектов: тензор второго ранга
δij и псевдотензор εjkl . Таким образом в тензорной алгебре всегда имеется
2 образующих. Тензоры четвертого ранга можно строить только из первой из них, так как вторая является антисимметричным псевдотензором
тензором. Получающийся на его основе тензор четвертого ранга εjkm εlnm
алгебраически выражается через δij , в виде δjl δkn − δjn δkl . Таким образом,
всего на основе тензора δij можно построить 3 алгебраически независимых
22
тензора четвертого ранга:
(1)
Tjkln = δjk δln ,
(2)
Tjkln = δjl δkn ,
(3)
Tjkln = δjn δkl
и взять их в качестве элементов базиса. Если тензор четвертого ранга можно разложить имеет в своем разложении только эти элементы базиса (в (3)
имеется три слагаемых), а скалярные функции зависят только от M2 , то
в этом случае мы будем говорить, что Tjkln является сферически симметричным в пространстве векторов M.
Мы в настоящей работе при описании класса уравнений вида (1) для
псевдовекторного поля будем использовать только такие элементы бази(a)
са – тензоры четвертого ранга Tjkln , которые являются мономами (при их
конструировании используется только операция тензорного умножения и
оператор свертки) в тензорной алгебре с образующими δij , εjkl , Mj , а ска(a)
лярные функции f (a) зависят только от M2 . В связи с этим мономы Tjkln
рассматриваются мы рассматриваем как эквивалентные, если они отличаются на множитель, возможно зависящий от M2 . 3)
Опишем кратко план дальнейшего изложения в этой главе.
1). На первом этапе строится наиболее широкое множество T линейно
(a)
независимых тензоров Tijkl , a = 1÷N четвертого ранга в тензорной алгебре
с порождающими элементами δij , εijk и Mj . Эти тензоры составляют базис
разложения тензора Tijkl общего вида. Такое построение решает задачу об
описании класса всех допустимых тензоров Tijkl в виде
Tijkl =
N
X
(a)
f (a) Tijkl .
(2.1.4)
a=1
Заметим при этом, что перестановка любой пары свободных индексов в
(a)
фиксированном элементе Tijkl приводит к другому ее элементу, если только
оба индекса не принадлежат тензору δ, либо псевдотензору ε.
2). После этого, из множества S всех псевдотензоров второго ранга
(a)
(a)
(a)
Sij = Tijkl ∇k Ml , Tijkl ∈ T, полученных посредством свертки по паре индексов выбирается базис линейно независимых пседотензоров второго ранга. Число элементов этого базиса будем обозначать посредством N̄ . Построение базиса S решает задачу об описании класса всех плотностей потоков
3)
В физике такое положение часто формулируется на интуитивном уровне в виде требования, чтобы
в физической системе не было никаких других выделенных векторов, кроме вектора Mj .
23
поля Mj на основе формулы
Sij =
N̄
X
(a)
f (a) Sij .
(2.1.5)
a=1
3) К каждому псевдотензору Sij вида (5) применяется операция свертки
с вектором ∇j . При этом оператор ∇j записывается слева от Sij . В результате, получается описание класса всех псевдовекторов Fi = ∇j Sij , которые
представляют поля Fi (x, t) «термодинамических» сил. Они получаются в
виде суммы
¯
N̄
X
(a)
Fi =
f (a) Fi .
(2.1.6)
a=1
(a)
Так как множество векторов Fi , которое получается посредством применения операции свертки к каждому элементу совокупности линейно независимых тензоров может приводить к совокупности линейно зависимых
(a)
векторов, то необходимо выделить из множества ∇j Sij линейно независимый набор F всех возможных термодинамических сил.
Описание класса всех термодинамических сил завершает описание класса всех эволюционных уравнений для псевдовекторного поля Mj , которые
имеют вид
Ṁi = Fi .
(2.1.7)
2.2. Построение множества линейно независимых тензоров Tijkl
Необходимо решить следующую задачу (см. [12]). Построить базис линейно независимых тензоров четвертого ранга в тензорной алгебре c образующими Mj , δij , εijk . Очевидно, что все элементы такого базиса могут
быть получены применением операций тензорного умножения и свертки.
Иными словами они являются элементами полугруппы с указанными порождающими элементами.
Прежде чем дать решение этой задачи, необходимо отметить следующее.
Символ Кронекера δ не должен использоваться в операции свертки, так
как, после такой операции по одному из индексов, он исчезает из конструируемого монома. По этой причине, множество всех тензоров Tijkl четвертого
ранга распадается на три класса: K0 состоит из мономов Tjkkl , в составе
которых символ Кронекера δ отсутствует, K1 состоит из тензоров, имеющих вид тензорного произведения символа Кронекера на тензор второго
24
ранга Tij , в состав которого этот символ не входит и класс K2 , который
состоит из тензоров четвертого ранга, которые представляются тензорным
произведением двух символов Кронекера.
Второе замечание заключается в том, что символ Леви-Чивита ε при
конструировании линейно независимых мономов не может быть использован в составе монома Tijkl более одного раза и как сомножитель в тензорном
произведении и как сомножитель в составе свертки. Это связано с тем, что
для тензорного произведения символов ε справедлива формула (см., например, [4])
δil δim δin
εijk εlmn = det δjl δjm δjn ,
δkl δkm δkn
которая сводит их к линейным комбинациям из символов δ. Следовательно,
множество всех искомых мономов разбивается на два класса: L0 состоит из
мономов, при построении которых отсутствует символ ε и мономов класса
L1 , которые содержат в своем составе символ ε.
Третье замечание касается того, что степень вхождения образующей Mj
в состав мономов класса L0 должна быть четной, а в состав мономов класса
L1 – нечетной, так как Mj является псевдовектором, а εklm – псевдотензором и только при четности суммарной степени их вхождения результирующий моном Tijkl тензором.
Наконец, укажем, что имеется только один линейно независимый скаляр Mj Mj в рассматриваемой нами полугруппе, а псевдоскаляры в ней отсутствуют. Это означает, что свертки Mj Mj учитываются при построении
базиса. При появлении таких скаляров в результате проведения операций
свертки онимогут быть включены в состав коэффициентов f (a) .
Прежде чем строить линейно независимые тензоры Tijkl сделаем еще
одно замечание. При описании программы решения задачи в пп.1-3. мы совсем не уделили исключили из рассмотрения возможные типы плотностей
потоков Sij , которые не содержат пространственных производных от Mj .
Это связано с тем, что имеется только два линейно независимых элемента, с точностью до принятого нами понятия эквивалентности, которые
являются результатами тензорного умножения и свертки из образующих
Mj и εijk , а именно, Mi Mj и εijk Mk . Но эти элементы тензорной алгебры
являются тензорами, а не псевдотензорами. Поэтому они не могут входить
в состав множества S линейно независимых плотностей потоков Sij .
Перейдем к построению линейно независимых тензоров Tijkl . Для то25
го чтобы полным образом перечислить все такие тензоры, необходимо их
убдобным образом проклассифицировать. Так как символ Леви-Чивитта
может входить в состав тензора Tijkl , то есть быть использован в операциях умноженяи и свертки, не более одного раза, то распределим, прежде
всего, все тензоры на два класса L0 и L1 . К классу L0 отнесем тензоры
в которых символ Леви-Чивитта εijk не используется, а к классу L1 – те
тензоры, в которых он используется один раз. Определим также классы
тензоров, которые классифицируют все перечисляемые тензоры согласно
признаку, сколько раз в них входит символ Кронекера. Число вхождений
k символа δij может принимать значения k = 0, 1, 2. Соответственно, согласно каждому из этих значений, введем классы K0 , K1 , K2 . Тогда легко
видеть, что класс всех всех тензоров представляется объединением
[ [
Lj ∩ Kl .
(2.2.1)
j=0,1 l=0,1,2
Таким образом, для перечисления класса всех тензоров, обладающих нужными свойствами нужно, описать каждый из классов, входящих в это объединение.
Очевидно, что L1 ∩ K2 = ∅, как как при тензорном умножении двух
символов Кронекера уже не остается индексов для сомножитель ε.
Рассмотрим класс K2 ∩ L0 . Как и выше, при тензорном умножении двух
символов Кронекера уже нет возможности для тензорного умножения на
псевдовектор Mj . Поэтому имеется 3 линейно независимых тензора этого
класса:
δij δkl , δik δjl , δil δjk .
(2.2.2)
Имеется только два линейно независимых тензора второго ранга, построенных на основе образующих M и ε, а именно, Mi Mj и εijk Mk , которые
пригодны для построения тензоров классов K1 ∩ L0 и K1 ∩ L1 .
Рассмотрим тензоры класса K1 ∩L0 . Эти тензоры перечисляются по признаку: какая пара из четырех индексов {i, j, k, l} связывается с символом
δ. После выбора такой пары остальные два индекса относятся к тензору,
который является тензорным произведением M M . Так как символ ε имеет
ранг 3, то он не может участвовать в построении тензоров этого класса,
так как для этого необходимо три индекса. В результате, класс K1 ∩ L0
26
состоит из шести элементов,
δij Mk Ml ,
δik Mj Ml ,
δil Mk Mj ,
δjk Mi Ml ,
δjl Mi Mk ,
δkl Mi Mj .
(2.2.3)
Точно также перечисляются шесть линейно независимых тензоров класса K1 ∩ L1 :
δij εklm Mm ,
δik εjlm Mm ,
δil εjkm Mm ,
δjk εilm Mm ,
δjl εikm Mm ,
δkl εijm Mm .
(2.2.4)
Рассмотрим теперь тензоры классов K0 ∩ L0 и K0 ∩ L1 .
При построении тензоров четвертого ранга класса K0 ∩ L0 нельзя использовать символы δ и ε. Поэтому имеется один тензор четвертого ранга
в этом классе, который порождается образующей M и он имеет вид тензорного произведения M M M M , т.е.
Mi Mj Mk Ml .
(2.2.5)
Тензоры класса K0 ∩ L1 распадаются на две группы. К одной из них
мы отнесем все тензоры, у которых все индексы символа ε являются свободными, и поэтому имеется только один индекс, который присваивается
псевдовектору M. В этой группе имеется следующие четыре тензора, согласно тому, что имеется только четыре возможности выбрать этот исключительный индекс:
Mi εjkl ,
Mj εikl ,
Mk εijl ,
Ml εijk .
(2.2.6)
Ко второй группе отнесем те тензоры, у которых один из индексов у
символа ε является связанным сверткой с псевдовектором M. Заметим,
что символ ε не может иметь два связанных индекса, так как εijk Mj mk ≡
0 (другой образующей первого ранга кромеM не имеется). В результате,
имеется шесть линейно независимых тензоров этой группы, согласно тому,
что имеется шесть возможностей выбора двух индексов из групы четырех
индексов {i, j, k, l} тех, по которым производится свертка с символом ε:
Mi Mj εklm Mm ,
Mi Mk εjlm Mm ,
Mj Mk εilm Mm ,
Mi Ml εjkm Mm ,
Mj Ml εikm Mm ,
27
Mk Ml εijm Mm .
(2.2.7)
Суммируем результаты проведенного в этом разделе анализа в виде отдельного утверждения.
Теорема 1. Множество T линейно независимых тензоров Tijkl четвертого ранга в тензорной алгебре с образующими M , ε, δ состоит из 26 тензоров,
которые перечисляются списками (2)-(7).
Согласно формуле (1.4), из этой теоремы следует что общий вид тензора
Tijkl :
26
X
(a)
Tijkl =
f (a) Tijkl ,
(2.2.8)
a=1
(a)
где тензоры Tijkl , a = 1 ÷ 26 даются списками (2)-(7).
2.3. Общий вид эволюционного уравнения дивергентного типа
для псевдовекторного поля
В этом разделе реализуется п. 2) программы исследования, описанной в
(a)
разд.2.1. Мы выделим из всего набора плотностей потоков Sij , a = 1 ÷ 26
максимально широкий набор линейно независимых плотностей. Необходимость решения такой задачи возникает в связи с тем, что при сворачи(a)
(a)
вании тензоров Tijkl ∇k Ml с линейно независимыми коэффициентами Tijkl ,
a = 1÷26 могут возникать линейно зависимые комбинации. Наличие таких
линейных комбинаций приводит к тому, что слагаемые в сумме
26
X
(a)
f (a) (M2 )Tijkl ∇k Ml
a=1
становятся линейно зависимыми, несмотря на произвольные скалярные
функции f (a) (M2 ). Решение задачи о выделении максимально широкого
набора линейно независиымх плотностей потоков решает задачу о выделении максимально широкого набора термодинамических сил в общем случае, когда коэффициенты f (a) (M2 ), a = 1 ÷ 26 зависят от M2 .
Заметим, что из тензоров класса K2 , сверткой с тензором ∇k Ml , получаются плотности потоков
δij (∇, M) ,
∇i Mj ,
28
∇j Mi .
Заменим этот набор следующим эквивалентным ему набором, в котором
тензоры второго ранга либо симметричны, либо антисимметричны (про(a)
цедура симметризации) и которые мы обозначим Sij , a = 1, 2, 3, соответственно с (2.2), в порядке их представления в приведенном списке:
∇i Mj + ∇j Mi ,
δij (∇, M) ,
∇i Mj − ∇j Mi .
(2.3.1)
Соответственно, из тензоров класса K1 ∩L0 получается набор плотностей
потоков (см. (2.3))
δij (M, ∇)M2 ,
Mj ∇i M2 ,
Mi ∇j M2 ,
Mj (M, ∇)Mi ,
Mi (M, ∇)Mj ,
Mi Mj (∇, M)
или эквивалентный ему набор, получаемый процедурой симметризации
δij (M, ∇)M2 ,
Mj ∇i M2 + Mi ∇j M2 ,
Mi ∇j M2 − Mj ∇i M2 ,
Mi (M, ∇)Mj + Mj (M, ∇)Mi ,
Mi (M, ∇)Mj − Mj (M, ∇)Mi ,
Mi Mj (∇, M) ,
(2.3.2)
(a)
которые мы обозначим Sij , a = 4 ÷ 9. Здесь и далее, как и приведенном
выше случае, номера для плотностей потоков даются в порядке их представления в приведенном списке.
Из тензоров класса K1 ∩ L1 получаются плотности потоков следующего
вида (см. (2.4)):
δij (M, [∇, M]) ,
εjkl Mk ∇i Ml ,
[M, ∇]j Mi ,
εikl Mk ∇j Ml ,
[M, ∇]i Mj ,
εijm Mm (∇, M)
или, после процедуры симметризации,
δij (M, [∇, M]) ,
εikl Mk ∇j Ml + εjkl Mk ∇i Ml ,
εikl Mk ∇j Ml − εjkl Mk ∇i Ml ,
[M, ∇]i Mj + [M, ∇]j Mi ,
[M, ∇]i Mj − [M, ∇]j Mi ,
(a)
εijm Mm (∇, M) ,
(2.3.3)
которые обозначаются Sij , a = 10 ÷ 15
Тензору класса K0 ∩ L0 (см. (5)) соответствует симметричная плотность
потока
(16)
Sij = Mi Mj (M, ∇)M2 .
(2.3.5)
29
Из тензоров класса K0 ∩ L1 получаются следующие плотности потоков
(см. (2.6), (2.7)) соответственно из первой и второй групп выражений для
тензоров Tijkl :
Mi [∇, M]j ,
Mi Mj (M, [∇, M]) ,
εijl (M, ∇)Ml ,
Mj [∇, M]i ,
Mi [M, (M, ∇)M]j ,
Mj [M, (M, ∇)M]i ,
εijk ∇k M2
Mi [M, ∇]j M2 ,
Mj [M, ∇]i M2 ,
εijm Mm (M, ∇)M2
или, после процедуры симметризации, которые мы обозначаем согласно
(a)
порядку следования Sij , a = 17 ÷ 26,
Mi [∇, M]j +Mj [∇, M]i , Mi [∇, M]j −Mj [∇, M]i , εijl (M, ∇)Ml , εijk ∇k M2 .
(2.3.6)
εijm Mm (M, ∇)M2 ,
Mi Mj (M, [∇, M]) ,
Mi [M, (M, ∇)M]j + Mj [M, (M, ∇)M]i ,
Mi [M, (M, ∇)M]j − Mj [M, (M, ∇)M]i ,
Mi [M, ∇]j M2 + Mj [M, ∇]i M2 ,
Mi [M, ∇]j M2 − Mj [M, ∇]i M2 . (2.3.7)
Лемма 1. Справедливо тождество
1 (20)
(18)
(19)
Sij = Sij + Sij .
2
(2.3.8)
Равенство (8) записывается в виде
1
εijl ∇l M2 = εijl (M, ∇)Ml + εijl [M, [∇, M]]l ,
2
а это последнее равенство проверяется сверткой его с εkij , что является
переходом к эквивалентному равенству ввиду антисимметрии по i, j обеих
его частей. Так как εkij εijl = 2δkl , то в результате свертки получаем
∇k M2 = 2(M, ∇)Mk + 2[M, [∇, M]]k .
30
(2.3.9)
Последнее слагаемое преобразуется следующим образом
[M, [∇, M]]k = εkij Mi εjlm ∇l Mm = (δkl δim − δkm δil )Mi ∇l Mm =
= Mm ∇k Mm − (M, ∇)Mk .
Подстановкой этого выражения в (9) убеждаемся в его тождественном выполнении.
Лемма 2. Справедливы тождества
1 (13)
(19)
(15)
Sij = Sij − Sij .
2
1 (15)
(14)
(18)
(19)
Sij = Sij + Sij − Sij .
2
(2.3.9)
(2.3.10)
Доказательство такое же как и доказательство Леммы 1. От равенства
(9), ввиду антисимметрии тензоров, входящих в его состав, переходим к
эквивалентному ему равенству сверткой εnij ,
εnij (εikl Mk ∇j Ml − εjkl Mk ∇i Ml ) = 2(δnl δjk − δjl δnk )Mk ∇i Ml =
= 2(M, ∇)Mn − 2Mn (∇, M) = εnij (εijl (M, ∇)Ml − εijl Ml (∇, M)) =
(19)
= εnij (Sij
(15)
− Sij ) .
Аналогично доказывается равенство (10),
εnij ([M, ∇]i Mj − [M, ∇]j Mi ) = 2[[M, ∇], M]n =
= ∇n M2 − 2Mn (∇, M) = εnij (εijl ∇l M2 /2 − εijl Ml (∇, M)) =
1 (20)
(15)
= Sij − Sij .
2
Следовательно,
1 (20)
(14)
(15)
Sij = Sij − Sij .
2
(20)
Подставляя в правую часть равенства выражение (8) для Sij получаем
(10).
Доказанные леммы 1 и 2 позволяют нам исключить из числа линейно
(13)
(14)
(20)
независимых плотностей потоков, плотности Sij , Sij , Sij . Проанализируем теперь плотности потоков, пропорциональные M 4 . Имеет место
31
Лемма 3. Справедливы тождества
(24)
Sij
=
(26)
Sij
1 (22)
(19)
Sij − M2 Sij .
2
(22)
= Sij
(2.3.11)
(20)
− M2 Sij .
(24)
(2.3.12)
(24)
Рассмотрим плотность Sij и вычислим εnij Sij ,
εnij Mi [M, (M, ∇)M]j − Mj [M, (M, ∇)M]i =
= εnij Mi εjkl Mk (M, ∇)Ml − Mj εikl Mk (M, ∇)Ml =
= (δnk δil − δnl δik )Mi Mk (M, ∇)Ml + (δnk δjl − δnl δjk )Mj Mk (M, ∇)Ml =
= Mn (M, ∇)M2 − 2M2 (M, ∇)Mn .
Сворачивая теперь
(24)
εijn εnkl Skl
(24)
= (δik δjl − δil δjk )Skl
(24)
= 2Sij
=
= εijn Mn (M, ∇)M2 − 2εijn M2 (M, ∇)Mn ,
Откуда следует тождество (11).
Точно также доказывается тождество (12). Вычисляя свертку
(26)
2
2
εnij Sij = εnij Mi [M, ∇]j M − Mj [M, ∇]i M =
2
= εnij Mi εjkl Mk ∇l M − Mj εikl Mk ∇l M
2
=
= (δnk δil − δnl δik )Mi Mk ∇l M2 + (δnk δjl − δnl δjk )Mj Mk ∇l M2 =
= 2Mn (M, ∇)M2 − 2M2 ∇n M2 .
Тогда
(26)
Sij
=
1
(26)
εijn εnkl Skl = εijn Mn (M, ∇)M2 − εijn M2 ∇n M2 =
2
(22)
= Sij
(20)
− M2 Sij .
Принимая во внимание доказанную лемму и Лемму 2, можно исклю(24)
(26)
чить из линейно независимых плотностей потоков плотности Sij и Sij .
32
Таким образом, нам остается проверить, что имеет место линейная незави(a)
симость плотностей потоков Sij с номерами a = 1 ÷ 12, 15 ÷ 19, 21 ÷ 23,
25. Обозначим это множество номеров посредством S.
(a)
Покажем, теперь, что все перечисленные псевдотензорные поля Sij (x),
a ∈ S линейно независимы. Для этого нужно проанализировать равенство
X
(a)
ca (M2 )Sij = 0
(2.3.13)
a∈S
при любых значениях индексов i, j = 1, 2, 3 и для любой пространственной
(a)
точки с радиус-вектором x, от которой зависят значения полей Sij , и показать, что оно возможно только при равенстве нулю всех коэффициентов ca ,
a ∈ S, которые понимаются как функции от M2 . Допустим, что (13) имеет
место. Тогда, заменив все векторное поле M, которое определяет значения
(a)
каждой из плотностей Sij , на λM, где λ – произвольный вещественный
множитель, из равенства (13), ввиду произвольности поля M, получаем
X
X
(a)
(a)
2
2
2
λ
ca (λ M )Sij + λ
ca (λ2 M2 )Sij +
a=10÷19, 6=13,14,16
a=1÷3
3
+λ
X
2
ca (λ M
2
(a)
)Sij
+λ
a=4÷9
(a)
X
4
ca (λ2 M2 )Sij = 0 .
(2.3.14)
a=16,21÷23,25
Ввиду линейной независимости функций λs с различными показателями
s ∈ N, из (14) следуют четыре равенства
X
(a)
ca (0)Sij = 0 ,
(2.3.15)
a=1÷3
(a)
X
ca (0)Sij = 0 ,
(2.3.16)
a=10÷19, 6=13,14,16
X
(a)
c0a M2 Sij +
a=1÷3
X
a=10÷19, 6=13,14,16
X
(a)
ca (0)Sij = 0 ,
(2.3.17)
a=4÷9
(a)
X
c0a M2 Sij +
(a)
ca (0)Sij = 0 .
(2.3.18)
a=16,21÷23,25
Проанализируем каждое из полученных равенств.
Лемма 4. Единственным решением уравнения (15) является c1 (0) =
c2 (0) = c3 (0) = 0.
(a)
Ввиду произвольности полей Sjk , подставим в равенство (15) их выражения в том случае, когда их значения вычисляются на основе фиксированного псевдовекторного поля Mj = Aj exp((n, x)), где Aj – произвольный
33
постоянный псевдовектор, а n – произвольный вектор. После элементарных
преобразований получаем равенство
c1 (0)δij (n, A) + c2 (0)Aj ni + c3 (0)Ai nj = 0 .
(2.3.19)
Ввиду произвольности векторов n и A, положим, что n ⊥ A. Тогда
(n, A) = 0 и из (19) получаем
c2 (0)Aj ni + c3 (0)Ai nj = 0 .
Сворачивая обе части этого равенства с ni , находим, что c2 (0)Aj n2 , что возможно только при c2 (0). Поэтому получаем, что имеет место c3 (0)Ai nj = 0,
что возможно только при c3 (0) = 0. Следовательно, равенство (19) возможно только при c1 (0)δij (n, A) = 0. Полагая, ввиду произвольности векторов,
(n, A) 6= 0, находим, что c1 (0) = 0. Таким образом, равенство (15) возможно при всех направлениях только при c1 (0) = c2 (0) = c3 (0) = 0.
Рассмотрим равенство (16), которое в подробной записи выглядит следующим образом:
c10 (0)δij (M, [∇, M]) + c11 (0)(εjkl Mk ∇i Ml + εikl Mk ∇j Ml )+
+c12 (0)([M, ∇]j Mi + [M, ∇]i Mj ) + c15 (0)εijm Mm (∇, M)+
+c17 (0)(Mi [∇, M]j + Mj [∇, M]i ) + c18 (0)(Mi [∇, M]j − Mj [∇, M]i )+
+c19 (0)εijl (M, ∇)Ml = 0 .
(2.3.20)
Это уравнение генерирует два независимых уравнения для симметричной
и антисимметричной части тензора, стоящего в его левой части. Запишем
симметризованное уравнение
c10 (0)δij (M, [∇, M]) + c11 (0)(εjkl Mk ∇i Ml + εikl Mk ∇j Ml )+
+c12 (0)([M, ∇]j Mi + [M, ∇]i Mj ) + c17 (0)(Mi [∇, M]j + Mj [∇, M]i ) = 0 .
(2.3.21)
Подстановка, аналогичная той, которая указана выше, приводит к равенству
(c12 (0) − c17 (0))(Ai [A, n]j + Aj [A, n]i ) = 0 ,
(2.3.22)
так как первые два слагаемых при такой подстановке тождественно обращаются в нуль.
34
Лемма 5. Имеют место равенства
c12 (0) = c17 (0) .
(2.3.23)
При условии [A, n] 6= 0 равенства (23) получаются последовательными свертками сначала с Aj , а затем – с [A, n]i .
Следствием этой леммы является вытекающее из (21) равенство
c10 (0)δij (M, [∇, M]) + c11 (0)(εjkl Mk ∇i Ml + εikl Mk ∇j Ml )+
+ c12 (0)([M, ∇]j Mi + [M, ∇]i Mj + Mi [∇, M]j + Mj [∇, M]i ) = 0 .
(2.3.24)
Лемма 6. Имеют место равенства
c10 (0) = c12 (0) = c17 (0) .
(25)
Вычислим след тензора в обеих частях равенства (24). В результате,
так как
εjkl Mk ∇j Ml + εjkl Mk ∇j Ml = εj ∇j Mk Ml = 0 ,
и (M[∇, M]) = ([M, ∇]M), получим
3c10 (0) + 4c12 (0) (M, [∇, M]) = 0 .
(2.3.26)
В общем положении (M, [∇, M]) 6= 0. Следовательно,
3c10 (0) + 4c12 (0) = 0 .
(2.3.27)
Свернем теперь обе части (24) с Mi Mj . Тогда, так как оператор
(M, [M, ∇]) ≡ 0, то получаем равенство
c10 (0) + 2c12 (0) M2 (M, [∇, M]) = 0 .
(2.3.28)
Следовательно, при тех же условиях ([M, ∇], M) 6= 0, имеем
c10 (0) + 2c12 (0) = 0 .
(2.3.29)
Из уравнений (27), (29) следует, что c10 (0) = c12 (0) = 0. Принимая во
внимание (23), получаем (25).
Следствие. c11 (0) = 0.
35
Используя равенства (25), уравнение (24) запишем в виде
c11 (0)(εjkl Mk ∇i Ml + εikl Mk ∇j Ml ) = 0 .
(2.3.30)
Выражение в скобках не равно тождественно нулю. В этом можно убедиться подстановкой суммы A(1) exp((n(1) , x)) + A(2) exp((n(2) , x)).
Рассмотрим, теперь, антисимметричную часть уравнения (20)
c15 (0)εijm Mm (∇, M) + c18 (0)(Mi [∇, M]j − Mj [∇, M]i )+
+c19 (0)εijl (M, ∇)Ml = 0 .
(2.3.31)
Лемма 7. Имеет место равенство
c15 (0) + c19 (0) = 0 ,
c18 = 0 .
(2.3.32)
Свернем обе части равенства (31) с εkij . Воспользовавшись тождеством εkij εijl = 2δkl , находим
c15 (0)Mk (∇, M) + c18 (0)εkij Mi [∇, M]j + c19 (0)(M, ∇)Mk = 0 .
(2.3.33)
Заметим, что εkij Mi [∇, M]j = ∇k M2 /2 − (M, ∇)Mk . Тогда
c15 (0)Mk (∇, M) + c18 (0)∇k M2 /2 + (c19 (0) − c18 (0))(M, ∇)Mk = 0 . (2.3.34)
Подставляя в это равенство M = A exp((n, x)) находим
c15 (0)Ak (n, A) + c18 (0)nk A2 + (c19 (0) − c18 (0))Ak (n, A) = 0 .
(2.3.35)
Ввиду произвольности векторов A и n, положим (A, n) = 0. Тогда
c18 (0)nk A2 = 0
(2.3.36)
и, следовательно, c18 (0) = 0. Поэтому из (34), в условиях (A, n) 6= 0, получаем c15 (0) + c19 (0) = 0.
Подстановка равенств (32) в (34) получаем уравнение
c15 (0)(Mk (∇, M) − (M, ∇)Mk ) = 0 .
Из него следует
Лемма 8. c15 (0) = c19 (0) = 0.
36
(2.3.37)
Существует поле M такое, что Mk (∇, M) 6= (M, ∇)Mk . Для этого
достаточно положить M = A(1) exp((n(1) , x)) + A(2) exp((n(2) , x)). Тогда из
(37) следует c15 (0) = 0.
Теорема 2. Единственным решением уравнения (16) является
ca (0) = 0 ,
a = 10 ÷ 19, 6= 13, 14, 16 .
Доказательство следует из Лемм 5-8.
Проанализируем равенство (21)
δij (c01 M2 (∇, M) + c4 (0)(M, ∇)M2 )+
+c02 M2 (∇i Mj + ∇j Mi ) + c03 M2 (∇i Mj − ∇j Mi )+
+c5 (0)(Mj ∇i M2 + Mi ∇j M2 ) + c6 (0)(Mi (M, ∇)Mj + Mj (M, ∇)Mi )+
+c7 (0)(Mi ∇j M2 − Mj ∇i M2 ) + c8 (0)(Mi (M, ∇)Mj − Mj (M, ∇)Mi )+
+c9 (0)Mi Mj (∇, M) = 0 .
(2.3.38)
Рассмотрим его симметричную часть
δij (c01 M2 (∇, M) + c4 (0)(M, ∇)M2 ) + c02 M2 (∇i Mj + ∇j Mi )+
+c5 (0)(Mj ∇i M2 + Mi ∇j M2 ) + c6 (0)(Mi (M, ∇)Mj + Mj (M, ∇)Mi )+
+c9 (0)Mi Mj (∇, M) = 0 .
(2.3.39)
Подставим в нее поле M(x) = A exp[i(n, x)]. В результате, получим
(c01 + 2c4 (0))δij A2 (A, n) + (c02 + 2c5 (0))A2 (Aj ni + Ai nj )+
+(c6 (0) + c9 (0))(A, n)Ai Aj = 0 .
(2.3.40)
Лемма 9. Справедливы равенства
c02 + 2c5 (0) = 0 ,
c01 + 2c4 (0) = 0 ,
c6 (0) + c9 (0) = 0 .
(2.3.41)
Ввиду произвольности векторов A и n, положим (A, n) = 0. Откуда
получаем первое равенство.
Далее, из (41) при (A, n) 6= 0 следует равенство
(c01 + 2c4 (0))δij A2 + (c6 (0) + c9 (0))Ai Aj = 0 ,
37
(2.3.42)
которое возможно только при выполнении других равенств (41).
Таким образом, равенство (39) превращается в следующее
0
2
2
c1 δij M (∇, M) − (M, ∇)M /2 +
0
2
2
2
+c2 M (∇i Mj + ∇j Mi ) − (Mj ∇i M + Mi ∇j M )/2 +
+c6 (0) Mi (M, ∇)Mj +Mj (M, ∇)Mi )−Mi Mj (∇, M) = 0 . (2.3.43)
Допустим, что тождество имеет место при каких-то значениях постоянных
c01 , c02 и c6 (0). Зафиксируем поле M в какой-либо пространственной точке
и добавим к этому значению малое поле
m = A exp((n, x)) .
(2.3.44)
Подставим в поле M + m в равенство (43). Тогда удерживая слагаемые
линейные по m имеем
0
2
c1 δij M (∇, m) − (M, ∇)(M, m) +
0
2
+c2 M (∇i mj + ∇j mi ) − (Mj ∇i (M, m) + Mi ∇j (M, m) +
+ c6 (0) Mi (M, ∇)mj + Mj (M, ∇)mi ) − Mi Mj (∇, m) = 0 ,
или, после подстановки (44)
0
2
c1 δij M (A, n) − (M, n)(M, A) +
0
2
+c2 M (Aj ni + Ai ni ) − (M, A)(Mj ni + Mi nj ) +
+ c6 (0) (Mi Aj + Mj Ai )(M, n) − Mi Mj (n, A) = 0 .
(2.3.45)
Полагая в этом равенстве (M, n) = (A, n) = (M, A) = 0, находим
c02 M2 (Aj ni + Ai ni ) = 0 .
Откуда имеем c02 = 0. Полагая в полученном, таким образом, равенстве
0
2
c1 δij M (A, n) − (M, n)(M, A) +
+ c6 (0) (Mi Aj + Mj Ai )(M, n) − Mi Mj (n, A) = 0 .
(2.3.46)
38
(M, n) = 0, (A, n) 6= 0, получим
c01 δij M2 − c6 (0)Mi Mj = 0 .
(2.3.46)
Ввиду произвольности вектора M, получаем c01 = c6 (0) = 0.
Принимая во внимание Лемму 9, находим, что нами доказана следующая
Лемма 10. Единственным решением уравнения (39) является
c01 = c02 = c03 = 0 ,
ca = 0 ,
a = 4 ÷ 9.
Рассмотрим антисимметричную часть уравнения (38),
c03 M2 (∇i Mj − ∇j Mi ) + c7 (0)(Mi ∇j M2 − Mj ∇i M2 )+
+c8 (0)(Mi (M, ∇)Mj − Mj (M, ∇)Mi ) = 0 .
(2.3.47)
Подставим в этом уравнение M = A exp((n, x)).
c03 (ni Aj − nj Ai ) + 2c7 (0)(Ai nj − Aj ni ) = 0 .
(2.3.48)
Ввиду произвольности векторов A и n, получаем c03 + 2c7 = 0. Тогда (47)
представимо в виде
0
2
2
2
c3 M (∇i Mj − ∇j Mi ) + (Mi ∇j M − Mj ∇i M )/2 +
+c8 (0)(Mi (M, ∇)Mj − Mj (M, ∇)Mi ) = 0 .
(2.3.49)
Применим тот же прием, который нами был использован при доказательстве Леммы 9. Зафиксировав вектор M, линеаризуем около него это уравнение.
0
2
c3 M (ni Aj − nj Ai ) + (M, A)(Mi nj − Mj ni ) +
+c8 (0)(M, n)(Mi Aj − Mj Ai ) = 0 .
(2.3.50)
Полагая, сначала, что (M, A) = (M, n) = 0, получим, что c03 = 0. Затем
при условии, что (M, n) 6= 0, получаем c8 (0) = 0. Таким образом доказана
Лемма 11. Единственным решением уравнения (47) является
c03 = c7 (0) = c8 (0) = 0 .
39
Из Лемм 10 и 11 следует следующая
Теорема 3. Единственным решением уравнения (17) является
c01 = c02 = c03 = 0 ,
ca = 0 ,
a = 4 ÷ 9.
Наконец, проанализируем уравнение (18). Выделим из него симметричную и антисимметричную части
c010 δij M2 (M, [∇, M]) + c011 M2 (εjkl Mk ∇i Ml + εikl Mk ∇j Ml )+
+c012 M2 ([M, ∇]j Mi + [M, ∇]i Mj ) + c017 M2 (Mi [∇, M]j + Mj [∇, M]i )+
+c16 (0)Mi Mj (M, ∇)M2 + c21 (0)Mi Mj (M, [∇, M])+
+c23 (0) Mi [M, (M, ∇)M]j + Mj [M, (M, ∇)M]i +
2
2
+ c25 (0) Mi [M, ∇]j M + Mj [M, ∇]i M = 0 . (2.3.51)
c015 εijm Mm (∇, M) + c018 (Mi [∇, M]j − Mj [∇, M]i )+
+ c019 εijl (M, ∇)Ml + c22 (0)εijm Mm (M, ∇)M2 = 0 .
(2.3.52)
Проанализируем симметричную часть (51) уравнения, снова применяя метод линаризации этого уравнения около постоянного поля M. В результате
линеаризации (51) получаем
c010 δij M2 (M, [∇, m]) + c011 M2 (εjkl Mk ∇i ml + εikl Mk ∇j ml )+
+c012 M2 ([M, ∇]j mi + [M, ∇]i mj ) + c017 M2 (Mi [∇, m]j + Mj [∇, m]i )+
+2c16 (0)Mi Mj (M, ∇)(M, m) + c21 (0)Mi Mj (M, [∇, m])+
+c23 (0) Mi [M, (M, ∇)m]j + Mj [M, (M, ∇)m]i +
+ 2c25 (0) Mi [M, ∇]j (M, m) + Mj [M, ∇]i (M, m) = 0 ,
(2.3.53)
а после подстановки m = A exp((n, x)) –
c010 δij M2 (M, [n, A]) + c011 M2 ([M, A]ni + [M, A]i nj )+
+c012 M2 ([M, n]j Ai + [M, n]i Aj ) + c017 M2 (Mi [n, A]j + Mj [n, A]i )+
+2c16 (0)Mi Mj (M, n)(M, A) + c21 (0)Mi Mj (M, [n, A])+
40
+c23 (0)(M, n) Mi [M, A]j + Mj [M, A]i +
+ 2c25 (0)(M, A) Mi [M, n]j + Mj [M, n]i = 0 , (2.3.54)
Пусть (M, n) = 0, A k n, то (53) превращается в уравнение
c011 + c012 = 0 .
(2.3.55)
Если же (M, n) = 0, A k M, то
c012 M2 − c017 M2 + 2c25 (0)(M, A) = 0 .
(2.3.56)
Полагая здесь M = ±A, найдем, что
c012 = c017 ,
c25 (0) = 0 .
(2.3.57)
Заметим также, что слагаемое с коэффициентом c16 имеет другой порядок
по |M|, по сравнению с другими слагаемыми. Поэтому нужно положить
c16 = 0. После подстановки полученных равенств в (54) находим
0
2
0
2
c10 δij M (M, [n, A]) + c11 M ([M, A]ni + [M, A]i nj )−
−([M, n]j Ai + [M, n]i Aj ) − (Mi [n, A]j + Mj [n, A]i ) +
+c21 (0)Mi Mj (M, [n, A])+
+c23 (0)(M, n) Mi [M, A]j + Mj [M, A]i = 0 ,
(2.3.58)
Если теперь положить, что все три вектора лежат в одной плоскости и
(M, n) = 0, но вектор A расположен произвольно по отношению к этим
векторам. Тогда получаем уравнение
0
2
c11 M ([M, A]ni + [M, A]i nj )−
−([M, n]j Ai + [M, n]i Aj ) − (Mi [n, A]j + Mj [n, A]i ) = 0 .
(2.3.59)
Полагая теперь в этом равенстве, что три вектора могут принимать произвольные значения, ввиду того, что векторное выражение в скобках не
равно тождественно нулю при любых векторах M, A, n, находим c011 = 0.
Наконец, рассмотри остальные слагаемые в (58),
c010 δij M2 (M, [n, A]) + c21 (0)Mi Mj (M, [n, A])+
41
+c23 (0)(M, n) Mi [M, A]j + Mj [M, A]i = 0 .
(2.3.60)
Полагая, что смешанное произведение (M, [n, A]) не равно нулю, вычислим
след от обеих частей уравнения (60). В результате, получим уравнение
3c010 + c21 (0) = 0 .
(2.3.61)
Если же свернуть уравнение (60) c Mi Mj , то получим уравнение
c010 + c21 (0) = 0 .
(2.3.62)
Из этих двух уравнений следует, что c010 = c21 (0) = 0. Следовательно, на
основании уравнения (60), получаем
c23 (0)(M, n) Mi [M, A]j + Mj [M, A]i = 0 .
(2.3.63)
Отсюда следует c23 (0) = 0.
Анализ уравнения антисимметричной части (52) уравнения (18) осуществляется точно таким же образом. В результате мы приходим к заключению, что справедлива следующая
Теорема 4. Единственным решением уравнения (18) является
c0a = 0 ,
a0 = 10 ÷ 19, 6= 13, 14, 16
ca = 0 ,
a = 16, 21 ÷ 23, 25 .
На основании Леммы 4 и Теорем 1-3, приходим к заключению о том, что
общий вид плотности Sij (см. (14)) потока магнитного момента состоит из
21 слагаемых и имеет вид
X
(a)
f (a) (M2 )Tijkl ∇k Ml .
(2.3.64)
a=1÷25;6=13,14,20,24
Все слагаемые являются в общем положении (при произвольных коэффциентных функциях f (a) (M2 ) являются линейно независимыми.
Следовательно, по этой же причине и на основании (2.1.1), справедливо
Следствие. Общее эволюционное уравнение для псевдовекторного поля M имеет вид
X
(a)
Ṁi =
∇j f (a) (M2 )Tijkl ∇k Ml .
(2.3.65)
a=1÷25;6=13,14,20,24
42
Глава 3. Эволюционное уравнение
для псевдовекторного унимодального поля
3.1. Построение множества линейно независимых
термодинамических сил
В этой главе мы завершаем решение той задачи, которая была поставлена в Главе 1. В этом разделе мы отберем линейно независимые термоди(a)
(a)
намические силы Fi = ∂Sij /∂xj соответствуют отобранным плотностям
(a)
Sij , a = 1 ÷ 25, a 6= 13, 14, 20, 24 описанным в Главе 2. Такая процедура
необходима в связи с тем, что вычисление сверток с каждым из линейно
независимых тензорных полей, которые представляют термодинамические
силы, может привести к тому, что список полученных таким образом полей
F (a) (x) окажется линейно зависимым набором.
Вычисление дивергенций всех возможных линейно независимых плот(a)
ностей потоков Sij , a = 1 ÷ 25, a 6= 13, 14, 20, 24, которые перечислены
формулами в разд. 2.3, приводит к следующим термодинамическим силам.
При этом дивергенции, соответствующие псевдотензорам Sij из формулы
(13), приводят к двум совпадающим значениям.
Псевдотензорам, представленным в (2.3.1), соответствуют силы F[x; M]
следующего вида:
∇i (∇, M) , ∆Mi .
(3.1.1)
Псевдотензорам, представленным формулой (2.3.2), соответствуют силы
∇i (M, ∇M2 ) ,
(∇, M∇i M2 ) ,
(∇, Mi ∇M2 ) ,
(∇, M(M, ∇)Mi ) ,
(∇, Mi (M, ∇)M) ,
(∇, Mi M(∇, M)) .
(3.1.2)
Псевдотензорам, представленным формулой (2.3.3), соответствуют силы
∇i (M, [∇, M]) ,
([∇, M], ∇i M) ,
εikl ∇j Mk ∇j Ml ,
(∇, [M, ∇]Mi ) ,
(∇, [M, ∇]i M) ,
(16)
Поле Fi , соответствующее псевдотензору Sij
(3.1.3)
(2.3.5) имеет вид
(∇, Mi M(M, ∇)M2 ) .
43
[∇, M(∇, M)]i .
(3.1.4)
Псевдотензорам, представленным формулой (2.3.6), соответствуют силы
(∇, Mi [∇, M]) ,
(∇, M[∇, M]i ) ,
[∇, (M, ∇)M]i .
(3.1.5)
Наконец, псевдотензорам, представленным формулой (2.3.7), соответствуют силы
(∇, Mi M(M, [∇, M])) ,
(∇, M[M, (M, ∇)M]i ) ,
(∇, Mi [M, (M, ∇)M]) ,
(∇, M[M, ∇]i M2 ) ,
(∇, Mi [M, ∇]M2 ) ,
[∇, M(M, ∇M2 )]i . (3.1.6)
Таким образом в окончательном списке линейно независимых термоди(a)
намических сил имеется 21 поле Fi [x; M], нумерация которых соответствует тому порядку, в котором они перечислены последовательно в формулах (1-6). В соответствии с этим общее эволюционное уравнение дивергентного типа для псевдовекторного поля M(x, t), сферически симметричного, и у которого термодинамические силы представляют собой локальные функционалы от поля M и содержат пространственные производные
не выше второго порядка, имеет вид
X
Ṁ =
F(a) [x; M] .
(3.1.7)
a=1÷25;6=13,14,20,24
3.2. Эволюционные уравнения с законом сохранения
плотности магнитного момента
Как известно (см., например, [1]), локальное термодинамическое состояние диэлектрической ферромагнитной среды, в пренебрежении ее механическими деформациями, изменениями со временем t концентраций имеющихся в ней примесных атомов, а также при пренебрежении проявлением в
ней магнитоэлектрических эффектов, полностью характеризуется локальной температурой T (x, t) и плотностью магнитного момента M(x, t) в каждой пространственной точке с радиус-вектором x. Таким образом, изменение со временем локального термодинамического состояния описывается
изменением физических полей T (x, t) и M(x, t). Нашей задачей является
изучение эволюции поля M(x, t) в пренебрежении тепловыми эффектами.
Мы должны среди всех термодинамических сил – слагаемых в правой
части эволюционного уравнения (1.7) отобрать только те, которые обеспечивают постоянство скалярного поля M2 (x, t) = M 2 = const, когда поле
если M(x, t) является его решением.
44
Проанализируем последовательно каждое из слагаемых в правой части
уравнения (37) и выделим из всего их списка те, которые удовлетворяют
нужному требованию. Прежде всего заметим, что поставленному требованию – обеспечению сохранения квадрата псевдовекторного поля M2 (x, t)
не отвечают силы с номерами a = 3, 4, 6, 15, 21, 23 из списка предыдущего раздела. Это связано с тем, что в выражения этих сил содержится
дифференцирование поля M2 (x, t) = M 2 .
Заметим, что термодинамическая сила с номером a = 12 в формуле (33)
соответствует уравнению (3) Главы 1. Она, заведомо обладает требуемым
свойством, так как Mi εikl ∇j Mk ∇j Ml ≡ 0. Поэтому мы сразу же удалим из
проводимого далее анализа.
Точно также удалим из списка термодинамическую силу с номером
a = 22, так как она также обладает требуемым свойством при условии
сохранения M2 (x, t) = M 2 . В самом деле, справедливы следующие тождественные преобразования
Mi (∇, M[M, (M, ∇)M]i ) =
= (∇, M)(M, [M, (M, ∇)M]) + Mi (M, ∇)[M, (M, ∇)M]i =
= Mi [(M, ∇)M, (M, ∇)M]i + Mi [M, (M, ∇)2 M]i = 0
с обращением в нуль первого слагаемого во второй строке и обоих слагаемых в первой.
Таким образом, нам нужно проанализировать следующие термодинамические силы F на предмет существования таких универсальных линейных
зависимостей между соответствующими им функциями от радиус-вектора
x, которые представляются скалярными произведениями (M, F):
∇i (∇, M) ,
∆Mi ,
∇i (M, [∇, M]) , ([∇, M], ∇i M) , (∇, [M, ∇]Mi ) , (∇, [M, ∇]i M) ,
[∇, M(∇, M)]i , (∇, Mi [∇, M]) , (∇, M[∇, M]i ) , [∇, (M, ∇)M]i .
(∇, M(M, ∇)Mi ) , (∇, Mi (M, ∇)M) , (∇, Mi M(∇, M)) .
(∇, Mi M(M, [∇, M])) , (∇, Mi [M, (M, ∇)M]) .
(3.1.8)
Эти термодинамические силы разбиты на четыре группы согласно тому, в
какой степени в каждую из них входит поле M. Наличие линейной зависимости, точно также как и в разделе 2.2, нужно проверять внутри каждой
45
из групп. В последующем анализе наличия линейной зависимости нумерацию выражений (M, F(a) ) и соответствующим из них коэффициентов ca
при принимаем согласно представленному здесь списку.
Итак, для квадратичных по полю M выражений проанализируем равенство
c1 (M, ∇(∇, M)) + c2 (M, ∆M) = 0 .
Подстановка в него поля M(x) = A exp[i(n, x)] дает равенство
c1 (A, n)2 + c2 n2 A2 = 0 ,
которое возможно только при c1 = c2 = 0, что усматривается полагая,
ввиду произвольности векторов A и n, сначала n ⊥ A, а затем n k A.
Проанализируем теперь возможность линейной зависимости между выражениями третьего порядка по M,
c11 Mi (∇, M(M, ∇)Mi ) + c12 Mi (∇, Mi (M, ∇)M)+
+c13 Mi (∇, Mi M(∇, M)) = 0 .
(3.1.9)
Аналогично, подстановка выражения для поля M(x) в виде плоской волны
с произвольной амплитудой A и волновым вектором n приводит к равненству
c11 + c12 + c13 = 0 .
(3.1.10)
Подставим теперь в (9) поле = [x, n]. Тогда (∇, M) = 0 и
(M, ∇)Mi = ni (n, x) − xi n2 .
Тогда, ввиду (∇, M) = 0, (M, n) = 0
(∇, M(M, ∇)Mi ) = (∇, M(ni (n, x) − xi n2 )) =
= (M, ∇(ni (n, x) − xi n2 )) = −Mi n2 .
Наряду с этим, по той же причине,
(∇, Mi (M, ∇)M) = (∇k Mi ) · (M, ∇)Mk + Mi (∇k M, ∇)Mk =
= εikl nl [M, n]k + εikl nl εkmn Mm nn =
= −Mi n2 + (δin δlm − δim δln )nl nn Mm = −2Mi n2 .
Подстановка этих соотношений дает
c11 + 2c12 = 0 .
46
(3.1.11)
Наконец, подставим в (39) поле M = x. Тогда
(∇, M(M, ∇))Mi = (∇, Mi (M, ∇)M) = (∇, Mi M) = 4xi
и поэтому из (9) следует
c11 + c12 + 3c13 = 0 .
Из этого равенства, вместе с (10) и (11) следует c11 = c12 = c13 = 0.
Проанализируем возможность линейной зависимости между двумя выражениями четвертого порядка порядка по M,
c14 Mi (∇, Mi M(M, [∇, M])) + c15 Mi (∇, Mi [M, (M, ∇)M]) = 0 .
(3.1.12)
Подставим в это равенство поле Mj (x) = Aj + ReBj exp[i(n, x)] с произвольными вещественными векторами A и n, а также комплексным вектором B, которого (A, B) = 0 и |ReB| = |ImB|, (ReB, ImB) = 0 так, что
B2 = 0. Такой выбор поля Mj (x) обеспечивает равенство Mj (x)Mj (x) =
const. Так как в этом случае
Mi (∇, Mi M(M, [∇, M])) = εjkl Mi ∇j Mi Mk Mm ∇m Ml =
= εjkl [M 2 ∇j Mk Mm ∇m Ml − (Mk Mm ∇m Ml )Mi ∇j Mi ] ,
где второе слагаемое обращается в нуль в предположении, что M2 = const,
и, кроме того,
Mi (∇, Mi [M, (M, ∇)M]) = εklm Mi ∇j Mi Mj Mk ∇l Mm =
= εklm [M 2 ∇j Mj Mk ∇l Mm − (Mj Mk ∇l Mm )Mi ∇j Mi ] ,
где второе слагаемое также обращается в нуль в том же предположении.
Тогда равенство (12) сводится к следующему
c14 εjkl ∇j Mk Mm ∇m Ml + c15 εklm ∇j Mj Mk ∇l Mm = 0 .
(3.1.13)
Подстановка указанного явного выражения для поля в это равенство дает
c14 εjkl nj nm [ImBm ei(n,x) ][ImBl ei(n,x) ]Mk +
i(n,x)
+ (Am + [ReBm e
i(n,x)
])[ReBl e
+c15 εklm nj nl [ImBm ei(n,x) ][ImBj ei(n,x) ]Mk +
47
]Mk +
i(n,x)
+ (Aj + [ReBj e
])[ReBm e
i(n,x)
]Mk = 0 .
Меняя во втором слагаемом индексы суммирования j ⇔ l во втором слагаемом и, сравнивая оба получившихся слагаемых, имеем
(c14 − c15 )εjkl nj nm (Am + [ReBm ei(n,x) ])[ReBl ei(n,x) ]Mk +
+ (c14 εjkl nm − c15 εjkm nl )nj [ImBm ei(n,x) ][ImBl ei(n,x) ]Mk = 0 .
(3.1.14)
Дифференцируя дважды по Ak и по Am получим равенство
(c14 − c15 )(εjkl nj nm + εjml nj nk ) = 0 .
Сверткой с nm находим (c14 − c15 )εjkl nj = 0. В силу произвольности вектора nj получаем c14 = c15 . В этом случае второе слагаемое в (14) также
обращается в нуль, в силу равенства нулю свертки симметричного и антисимметричного тензоров.
Воспользуемся полученным равенством коэффициентов в уравнении
(13) и покажем, что они должны быть равны нулю. Для этого нужно показать, что равенство
εjkl ∇j Mk Mm ∇m Ml + εklm ∇j Mj Mk ∇l Mm = 0 .
(3.1.15)
не может иметь место для любых гладких полей Mj (x), удовлетворяющих
условию Mj Mj = M 2 = const.
Совершим во втором слагаемом замены индексов суммирования: сначала m ⇔ j, а затем l ⇒ j. В результате, получим
h
i
εjkl ∇j Mk Mm ∇m Ml − ∇m Mk Mm ∇j Ml = 0 .
Подставим в это равенство поле Mj в виде Mj = Ajk xk с произвольной
невырожденной матрицей A, (A)jk = Ajk det A = 0. Такое поле не удовлетворяет поставленному условию Mj Mj = const, однако в случае, если приведенное равенство выполняется, то оно обязано выполняться и для линейных полей Mj = Ajk xk , у которых матрица A является антисимметричной,
так как в любой точке x, например в нулевой, где оно выполняется, оно также должно выполняться и при малых отклонений от этой точки. При этом
в силу сохранения Mj Mj для матрицы должно иметь место A + AT = 0.
В результате указанной подстановки получаем
h
i
2
εjkl Akj (AM)l + Mk (A )lj − (AM)k Alj − Mk Alj SpA = 0 .
48
Положим матрица A антисимметрична. Тогда A2 – симметричная матрица
и SpA = 0. Поэтому
i
h
εjkl Akj (AM)l − (AM)k Alj = 0
вследствие равенства нулю свертки симметричного и антисимметричного
тензоров. Далее, ввиду того, что xj – произвольный вектор и det A 6= 0,
то Mj – произвольный вектор. Тогда, по той же причине, AM – произвольный вектор, и потому дифференцированием по этому вектору, получаем равенство εjkn Akj = εjnl Alj = 0, которое должно быть верным для
любой антисимметричной матрицы A. Из него следует, что εjkn Akj = 0.
Откуда сверткой с εnlm получаем Alm = Am l, что противоречит антисимметрии матрицы A. Полученное противоречие показывает, что равенство
(15) невозможно для любых полей Mj , и поэтому c14 = c15 = 0.
Перейдем теперь к анализу возможности наличия тождеств в списке
(8) выражений (M, F), имеющих степень два по полю Mj . Прежде всего
заметим, что в списке таких выражения имеются тождественно равные
нулю при выполнении условия M2 = const. А именно, таковыми являются
выражения с номерами a = 4, 5, 8. В самом деле, для выражения с силой
под номером a = 4 имеем
Mi ([∇, M], ∇i M) = Mi εjkl ∇k Ml ∇j Mi =
= εjkl ∇k Ml Mi ∇j Mi − Ml εjkl (∇k Mi )(∇j Mi ) = 0 ,
где оба слагаемых, полученных последовательным применением оператора ∇k равны нулю; первое вследствие предположения Mi Mi = const,
Mi ∇j Mi = 0, а второе – вследствие свертки симметричного тензора с антисимметричным.
Для выражения с силой под номером a = 5, последовательно применяя
оператор ∇j к каждому сомножителю, имеем
Mi (∇, [M, ∇]Mi ) = Mi ∇j εjkl Mk ∇l Mi =
= εjkl (∇j Mk )Mi ∇l Mi + Mi εjkl Mk ∇j ∇l Mi = 0
где первое слагаемое равно нулю в силу Mi ∇j Mi = 0, а второе – ввиду
свертки симметричного и антисимметричного тензоров. По тем же причинам, равно нулю выражение, соответствующее силе под номером a = 8,
Mi (∇, Mi [∇, M]) = Mi ∇j Mi εjkl ∇k Ml =
49
= Mi (∇j Mi )εjkl ∇k Ml + Mi Mi εjkl ∇j ∇k Ml = 0 .
Таким образом следующие термодинамические силы Fi из списка (8)
обеспечивают сохранение квадрата поля M:
([∇, M], ∇i M) ,
(∇, [M, ∇]Mi ) ,
(∇, Mi [∇, M]) .
(3.1.16)
Проанализируем, теперь, возможность существования линейной зависимости между скалярными функциями, представленными оставшимися
выражениями (M, F):
Mi ∇i (M, [∇, M]) , Mi (∇, [M, ∇]i M) , Mi [∇, M(∇, M)]i ,
Mi (∇, M[∇, M]i ) , Mi [∇, (M, ∇)M]i .
(3.1.17)
Подставим в возможное соотношение линейной зависимости между ними
α1 Mi ∇i Mj εjkl ∇k Ml + α2 Mi ∇j εikl Mk ∇l Mj + α3 Mi εijk ∇j Mk ∇l Ml +
+α4 Mi ∇j Mj εikl ∇k Ml + α5 Mi εijk ∇j Ml ∇l Mk = 0 .
(3.1.18)
линейное выражение для поля Mj = Ajk xk с произвольной невырожденной
матрицей A. Тогда имеем равенство
α1 (A2 )jn εjkl Alk + α2 Ain εikl (A2 )kl + α3 Ain εijk Akj SpA+
+α4 Ain εikl Alk SpA + α5 Ain εijk (A2 )kj = 0 .
Так как матрица A антисимметрична, то SpA = 0 и A2 симметрична. Тогда это равенство превращается в α1 (A2 )jn εjkl Alk = 0. Предположим, что
α1 6= 0. Следовательно, так как матрица A невырождена (все ее собственные числа чисто мнимые), то у матрицы A2 также нет нулевых собственных чисел и поэтому εjkl Akl не может быть ее собственным вектором с
нулевым собственным числом. Отсюда следует, что εjkl Alk = 0 и сверткой
обеих частей этого равенства с εjmn находим, что Amn = Anm . Полученное
противоречие доказывает, что α1 = 0.
(1)
Учитывая равенство α1 = 0 подставим в (18) поле Mj (x) = Aj + Mj ,
(1)
Mj = ReBj exp[i(n, x)] с теми же что и ранее использованными свойствами коэффициентов A и B. Такая подстановка приводит к следующему
(2)
(2)
Mi εikl α2 nl ∇j Mk Mj − α3 nj ∇l Mk Mj +
+
(2)
α4 nk ∇j Mj Ml
50
−
(2)
α5 nj ∇l Mj Mk
= 0,
(1)
где ∇k Mj
(2)
= nk Mj , и далее,
(2)
(1)
(2)
Mi εikl (α2 − α3 )nl nj (Mk Mj − Mk Mj )+
(2)
(2)
(1)
(2)
(2)
(1)
+α4 nk nj (Ml Mj − Mj Ml ) − α5 nj nl (Ml Mk − Mk Ml ) = 0 . (3.1.19)
Положим в полученном равенстве, сначала, n ⊥ M(1) при любом значении радиус-вектора x, от которого зависит M(1) и, следовательно, n ⊥
M(2) . Тогда из него следует, что
(1)
α4 Ai εikl nk Ml nj Aj = 0 .
Выбрав вектор n так, чтобы (n, A) 6= 0 и при этом n 6 k A, получим, что
(1)
Ai εikl nk Ml 6= 0, так как A ⊥ M(1) . Отсюда следует, что α4 = 0. Учитывая это равенство в (19), подставим в него вектор n, перпендикулярный
A. Тогда слагаемое с α5 исчезает, и поэтому из получившегося равенства
следует, что
(2)
(2)
(1)
Mi εikl (α2 − α3 )nl nj (Mk Mj − Mk Mj ) = 0 .
(3.1.20)
Выбрав теперь вектор n, получим из этого равенства
h
i
(2)
(2)
(1)
(1)
(1)
(2)
(2)
(1)
(α2 − α3 )εikl nl nj Ai (Mk Mj − Mk Mj ) + Mi (Mk Mj − Ak Mj ) =
h
i
(2)
(2)
(1)
(2)
(2)
= (α2 − α3 )εikl nl nj Ai (Mk Mj + Mi (Mk Mj =
h
i
(2)
(2)
(1)
(2)
= (α2 − α3 )(n, M ) (A, [M , n]) + (M , [M , n]) = 0 .
Это равенство должно выполняться при любом векторе n, который мы выберем так, чтобы (n, M(2) ) 6= 0, более того, мы его выберем так, чтобы он
был перпендикулярен ни ReB, ни ImB. Но тогда (M(1) , [M(2) , n]) 6= 0 и
равенство не может выполняться при α2 6= α3 для любого радиус-вектора
x, так как это слагаемое в скобках содержит квадратичную комбинацию
относительно exp[i(n, x)], а первое слагаемое – имеет первую степень относительно этой экспоненты. Отсюда следует, что α2 = α3 .
Покажем, наконец, что
Mi εikl ∇j Mk ∇l Mj − ∇l Mk ∇j Mj = 0
(3.1.21)
для полей, удовлетворяющих условию Mj Mj = M 2 . Во первых, учитывая
равенство нулю свертки εikl Mi Mk , левая часть (51) записывается в виде
Mi εikl (∇j Mk ) · (∇l Mj ) − (∇l Mk ) · (∇j Mj ) .
51
Введем матрицу Akl = ∂Mk /∂l . Тогда, в терминах этой матрицы, последнее
выражение имеет вид
2
Mi εikl (A )kl − Akl SpA .
Дифференцируя же по xk условие Mj Mj = M 2 имеем дополнительную
связь матрицы A с вектором Mj , Mj Ajk = 0, k = 1, 2, 3. Запишем это
равенство в системе координат, где Mj = 0, j = 2, 3, M1 6= 0,
A1k = 0 , k = 1, 2, 3 .
(3.1.22)
Тогда, в силу этих равенств, исследуемое выражение, действительно, обращается в нуль,
2
ε1kl (A )kl − Akl SpA =
= (A21 A13 − A31 A12 ) + (A22 A23 − A32 A22 ) + (A23 A33 − A33 A32 ) −
− (A23 − A32 )(A11 + A22 + A33 ) =
= (A21 A13 − A31 A12 ) − (A23 − A32 )A11 = 0 .
Следовательно, (51) имеет место в выбранной нами системе координат, и
поэтому оно верно, в следствие его инвариантности, в любой системе координат.
Теорема 1. Весь класс термодинамических сил F, который обеспечивают выполнение равенства (M, M) = 0 дается следующим списком
εikl ∇j Mk ∇j Ml ,
εikl ∇j Mj Mk Mm ∇m Ml ,
∇j Mi εjkl ∇k Ml ,
εjkl ∇k Ml ∇i Mj , ∇j εjkl Mk ∇l Mi ,
εikl ∇j Mk ∇l Mj − ∇l Mk ∇j Mj .
52
Заключение
В работе решена задача об описании класса всех эволюционных уравнений с дифференциальным эволюционным оператором по пространственным производным дивергентного типа для псевдовекторного поля M(x).
При этом на возможный выбор эволюционного оператора наложены дополнительные условия: он должен содержать пространственные производные
не выше второго порядка, должен быть ковариантным при преобразованиях группы O3 (повороты и отражения евклидового пространства) и при
этом должен быть сферически симметричным. Далее, из этого класса эволюционных уравнений выделен класс таких из них, которые обладают инвариантом M2 (x, t) = M 2 .
В результате проведенного исследования оказалось, что все уравнения
исследуемого класса имеют следующий вид
Ṁi = γ1 εikl ∇j Mk ∇j Ml + γ2 εikl ∇j Mj Mk Mm ∇m Ml +
+γ3 εjkl ∇k Ml ∇i Mj + γ4 ∇j εjkl Mk ∇l Mi + γ5 ∇j Mi εjkl ∇k Ml +
+γ6 εikl ∇j Mk ∇l Mj − ∇l Mk ∇j Mj .
(1)
с произвольными постоянными γa и a = 1 ÷ 6. В частности, если все постоянные γa = 0 при a = 2÷6 из этого общего уравнения получается сферически симметричное уравнение Ландау-Лифшица для плотности магнитного
момента ферромагнетика в отсутствии внешнего магнитного поля.
В том случае, если накладывается дополнительное условие на термодинамические силы, которое состоит в том, что требуется сохранение соленоидальности поля M, правая часть последнего уравнения упрощается
таким образом, что остаются только два линейно независимых слагаемых
Ṁi = γ(εikl ∇j Mk ∇j Ml −εjkl ∇k Ml ∇i Mj )+γ 0 (εjkl ∇j Mk ∇l Mi −εikl ∇j Mk ∇l Mj ) .
(2)
Заметим, что найденные эволюционные уравнения, в общем случае,
не обладают свойством инвариантности относительно замены t ⇒ −t,
M ⇒ −M. А именно, такой тип инвариантности, которым обладает уравнение Ландау-Лифшица, нарушается в том случае, когда хотя бы одна из
постоянных γa , a = 3÷6 отлична от нуля. Тогда можно ожидать, что наличие таких слагаемых приводит к эволюционным уравнениям, которые не
обладают обратимостью движения.
53
Литература
1. Трусделл К. Первоначальный курс рациональной механики сплошных
сред / М.: Мир, 1975.– 592 с.
Truesdell C. A first Course in Rational Continuum Mechanics / Baltimore,
Maryland: The John Hopkins University, 1972.
2. Исаев А.А., Ковалевский М.Ю., Пелетминский С.В. Гамильтонов подход к теории антиферромагнитных систем // ТМФ. – 1993. – 95:1. –
C.58-73.
3. Исаев А.А., Ковалевский М.Ю., Пелетминский С.В. Гамильтонов подход в теории конденсированных сред со спонтанно нарушенной симметрией / Физика элементарных частиц и атомного ядра.– 1996. – 27. –
2. – C.431-492.
4. Рашевский П.К. Риманова геометрия и тензорный анализ / М.: Наука,
1967. – 664 c.
5. Мак-Коннел А.Дж. Введение в тензорный анализ с приложениями к
геометрии, механике и физике / М.: Гос. изд. ФМЛ, 1963. – 412 c.
6. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред / Теоретическая физика т.8 / М.: Наука, 1982. – 620 c.
7. Ахиезер А.И., Барьяхтар В.Г., Пелетминский С.В. Спиновые волны /
М.: Наука, 1967. – 368 c.
8. Академик НАН Украины Виктор Григорьевич Барьяхтар. Жизнь в
науке / Нац. акад. наук Украины, Нац. науч. центр «Харьк. физ.-техн.
ин-т» / К. : Наукова думка, 2010. - 328 с.
9. Вирченко Ю.П. Основное состояние векторной решеточной модели //
Belgorod State University Scientific Bulletin. Mathematics & Physics. –
2012. – 23(142);29. – С.54-66.
10. Клюев А.С., Вирченко Ю.П. Основное состояние векторной решеточной модели с парным взаимодействием. Cлучай вырожденного
обменного интеграла // Belgorod State University Scientific. Bulletin
Mathematics & Physics. – 2014. – 5(176);34. – С.126-133.
54
11. Ландау Л.Д., Лифшиц Е.М. Гидродинамика / Теоретическая физика
т.6 / М.: Наука, 1986.– 732c.
12. Вирченко Ю.П., Д.А. Чурсин Плотность потока магнитного момента сферически симметричного магнетика // Belgorod State University
Scientific Bulletin. Mathematics & Physics. – 2015. – No11(208); 39. –
С.191-196.
13. Гантмахер Ф.Р. Теория матриц / М.: М.: Наука, 1966. – 576 c.
14. Барьяхтар В.Г., Белых В.Г., Соболева Т.К. Макроскопическая теория
релаксации коллективных возбуждений в неупорядоченных и неколлинеарных магнетиках / Теор. и мат. физика.– 1988.– 77;2.– C.311Џ318.
55
Отзывы:
Авторизуйтесь, чтобы оставить отзыв