Содержание к диссертации
Стр,
Введение г 4
ГЛАВА I. МОДЕЛЬ ДИСКРЕТНЫХ ОРИЕНТАЦИИ 25
1. Классификация собственных функций
оператора Фоккера-Планка по представ
лениям группы симметрии куба 25
2. Качественный анализ младших собст
венных функций оператора 31
3. Вывод модели дискретных ориентации 37
ГЛАВА П. ПОГРАНСЛОЙНЫЕ СОБСТВЕННЫЕ ФУНКЦИИ
ОПЕРАТОРА "МАЛОЙ ДИФФУЗИИ" 51
1. Главный член погранслойной асимпто
тики собственных функций: вывод и реше
ние "параболического уравнения" 51
2. Построение полного асимптотического
разложения типа внутреннего погра
ничного слоя в окрестности пере
вальной точки 56
3. Асимптотика главных членов разложе
ния в окрестности вершины 67
ГЛАВА Ш. ЛОКАЛЬНЫЕ РАЗЛОЖЕНИЯ В ОКРЕСТНОСТИ
ВЕРШИНЫ 77
1. Постановка задачи 77
2. Пересечение двух пограничных слоев 81
3. Сшивание асимптотик 88
4. Угловой пограничный слой 93
5. Однозначная разрешимость уравнений
3.
углового пограничного слоя ГЛАВА ІУ. АОШШОТШШСШ. СПЕКТРАЛЬНЫЙ АНАЛИЗ ^ ' 1. Вариационные оценки низших собст в "среднем" 115 3. "Почти" инвариантность собственного 4. Форма линии поглощения. Асимптотика билинейной формы резольвенты оператора 132 ПРИЛОЖЕНИЕ 145 ЛИТЕРАТУРА 158 4. Введение к работе Асимптотический спектральный анализ сингулярно возмущенного оператора Фоккера-Нланка или соответствующего диффузионного оператора с малым параметром при старших производных представляет собой важную задачу математической физики [Я, Интерес к эллиптическим и параболическим дифференциальным уравнениям с малым параметром тесно связан с задачей изучения малых случайных возмущений динамических систем Й. Впервые такая задача ставилась в статье Понтрягина, Андронова и Витта[з}, результаты которой относились к одномерным и, частично, к двумерным динамическим системам. В сущности, речь идет об асимптотических задачах, возникающих при стремлении к нулю параметра, характеризующего малость возмущения. Малые возмущения на больших интервалах времени, вообще говоря, оказывают существенное влияние на поведение системы. Так, например, даже если соответствующая детерминированная система имеет асимптотическое положение равновесия, эффект малых случайных возмущений состоит в том, что все траектории системы покидают область притяжения соответствующего положения равновесия с вероятностью І. В случае, когда возмущение динамической системы X - fecx)представляет собой гауссовский процесс белого шума, то возмущенное движение становится случайным процессом 'X-t , который удовлетворяет стохастическому дифференциальному уравнению w'X.fc- *>(**)* 6W dvJj ф Здесь W-t - Ч, -мерный винеровский процесс (броуновское движение), %W - векторное поле сноса, - малый вещественный параметр. К задачам, которые позволяют учесть влияние малых случай- ных возмущений, относятся задачи о предельном поведении стационарного распределения диффузионного процесса и задача о выходе процесса из области. Чтобы проиллюстрировать природу явления, удобно представить себе случайные возмущения как диффузию частиц в детерминированном поле сноса бСх) . Следуя Вентцелю и Фрейдлину [4^, различают три случая согласно тому, как диффундируют частицы а) вместе с полем в) поперек поля и с) против поля а) в) с) Для диффузии типа а) первые результаты были получены Левинсо-ном \.б\ Для диффузии типа в) Хасьминский И вычислил асимптотику распределения вероятности тех точек на границе области, в которых траектории процесса впервые ее покидают (решение проблемы Колмогорова). Для наиболее сложного случая с) первые строгие результаты были получены вероятностными методами Вент-целем и Фрейдлиным М, которые показали, что время выхода растет экспоненциально с уменьшением коэффициента диффузии. Более полные результаты были получены Людвигом ВД. Используя формальные разложения, он вывел асимптотическую формулу для среднего времени выхода и минимального положительного собственного значения соответствующего эллиптического оператора Zg, . Теоретически движение вектора намагниченности суперпара- магнитных частиц рассмотрели Неель и Броун [_II-I3JS. Они ввели следующую модель движения вектора М в отсутствие поля: движению И сопоставляется диффузионный процесс на сфере постоянного радиуса \ И I . Наиболее важное предположение, сделанное Броуном, состоит в том, что случайные тепловые воздействия имеют времена корреляции гораздо меньшие, чем временные характеристики рассматриваемой динамической системы. Это позволяет считать случайные воздействия на вектор И случайным процессом со свойствами белого шума и использовать для дальнейшего расчета уравнение Фоккера-Планка. Броун [і2,із}, а вслед за ним Аарони [j5^ рассматривали уравнение диффузии на всей сфере, которое допускает разделение переменных. Броун [із"] проанализировал решение уравнения диффузии для различных предельных случаев и получил численные оценки для первого собственного значения и его асимптотическое поведение при малых коэффициентах диффузии по сравнению с величиной поля сноса. В работах [16-18] Аарони впервые рассмотрел релаксацию в частицах с кубической магнитной анизотропией (в соответствующем уравнении Фоккера-Еланка переменные уже не делятся). Путем расчета на ЭВМ он получил зависимость первых трех собственных чисел от коэффициента диффузии и обнаружил, что в случае отрицательной константы анизотропии выполняется соотношение 1:2:3 между ними при малых коэффициентах диффузии. й Экспериментальное изучение магнитных свойств мелких частиц получило широкое развитие после работ Вина [І4І, который назвал такие частицы суперпарамагнитными (СП), а термин СП используется для описания поведения магнитных свойств этих При изучении явления суперпарамагнетизма методом ЯП? наблюдаемая форма спектра возникает в результате усреднения некоторых величин (функционалов) по всевозможным траекториям процесса Ml"t) . Информация об этих траекториях в принципе имеется в уравнении диффузии, но она там содержится в "зашифрованном" виде, неудобном для применения. На практике обычно используют более простую модель, в которой предполагается, что М vt) может принимать лишь конечное множество значений (состояния намагниченности вдоль одной из "легких" осей), а процесс движения представляет собой марковский процесс с конечным числом состояний, так называемую модель дискретных ориентации (ЩО). Этот процесс описывается с помощью времен релаксации, обратные величины которых пропорциональны вероятностям перехода вектора М из одного состояния (от одной "легкой" оси) в другое (к другой "легкой" оси). В работах Белозерского и Павлова JI9-2I] была предпринята попытка обоснования ВДО на основе анализа времен выхода диффузионного процесса из специальных областей и их комбинаций на единичной сфере, которые они называют фундаментальными. В работе Белозерского, Павлова и автора \22~] было показано, что вопрос обоснования ВДО - чисто спектральный вопрос, связанный в первую очередь с наличием серии экспоненциально малых собственных чисел оператора Фоккера-Шганка. Там было показано, что ВДО естественным образом возникает, если следить за эволюцией функции распределения ориентации магнитного момента в определенном масштабе времени и ограничиться несколькими частиц. членами ее разложения в ряд Фурье по собственным функциям оператора Фоккера-Бланка. При этом следует сохранять те ее компоненты, которые отвечают экспоненциальной серии собственных чисел. Тогда величина первого неучтенного собственного значения по отношению к учтенным оказывается удобным критерием эффективности замены диффузионной модели на ВДО. Примерно в то же время Матковский и Шусс L23-25] предприняли попытку исследования собственных функций диффузионного оператора с малым параметром при старших производных, отвечающих экспоненциально малым собственным значениям. Они свели проблему разыскания собственных функций к построению погран-слойных решений однородного диффузионного уравнения. Методом параболического уравнения они получили главный член асимптотики решения. В последней из работ 125 j Матковский и Щусс интересовались лишь главным собственным значением и соответствующей собственной функцией, которые определяют лишь скорость стремления решения уравнения Фоккера-Планка к равновесному распределению. В теории сунерпарамагнетизма этот этап еще не является завершающим. Дело в том, что вектор намагниченности суцерпара-магнитной частицы имеет несколько равноценных положений равновесия. Такая симметрия задачи приводит к появлению целой серии экспоненциально малых собственных чисел оператора Фоккера -Бланка, имеющих один порядок малости. Именно это обстоятельство позволяет смотреть на релаксацию функции распределения как на марковский процесс с конечным числом состояний в определенном масштабе времени, т.е. приводит к ВДО. Настоящая работа посвящена обоснованию перехода от диффузионной модели к МДО на основе подродного анализа асимптотики низших собственных значений диффузионного оператора и соответствующих собственных функции, а также оценке погрешности в форме линии поглощения при замене диффузионной модели на МДО. Первая глава начинается с исследования основного оператора диффузионной модели 2е с малым параметром при старших производных ч5~=-ЄД<\іЧ ?4>?\ґ, (В.і) который действует в пространстве L^(jc.) квадратично суммируемых функций на единичной сфере Sd . Соответствующее уравнение для функции распределения магнитного момента имеет вид (В.2) -=* Здесь Ч> - потенциал поля сноса ^ - - yvp Ч>= Vi-C^+xV^V), «Ч^г*, і (в.з) Координатные плоскости X = О, U = О, ^ = 0 разбивают единичную сферу на восемь сферических треугольников, которые называются фундаментальными областями Ьс.^ , I = 1,2,...,8. Симметрия задачи позволяет свести изучение оператора 2. на всей сфере к исследованию краевых задач в какой-либо из фундаментальных областей , скажем, в области 2L *г= А \г (B.4) *v = 0 либо^^ = 0 на каждой из сторон сферического треугольника ->Сі. Оказывается, что собственные функции таких краевых задач, отвечающие экспоненциально малым собственным числам -Д<0Ч\^Р\Г-А\Гр А=0(І"^), W>0 (В.5) практически постоянны во внутренней части фундаментальной области 0^=. oov4sir (В.6) Асимптотика младших собственных функций подробно исследована в гл. ПиШ. В гл.I мы строим ВДО, опираясь на эти результаты. Именно, заметим, что линейные комбинации младших собственных функций оператора *& ^-X feutf* (в-7) образуют базис в пространстве N$ , натянутом на множество собственных функций этого оператора, отвечающих экспоненциально малым собственным числам. На основании анализа главного члена асимптотики младших собственных функций оператора 21 матрица О в (В.7) может быть выбрана так, что функции иГ; в главном порядке своей асимптотики удовлетворяют соотношениям ^1 = tu і Я_х< (В.8) Матрица М сужения оператора -* ]K/g в новом базисе ^U^Y, задает систему дифференциальных уравнений - К (В.9) Эти уравнения и определяют МДО при малом . Вычисленная в гл.ІУ асимптотика низших собственных значений оператора 2! позволяет заключить, что перенормированная матрица и- \% U) имеет в качестве предела при -> 0 матрицу, которая оказывается марковской ?iw Н V Ю = L. (B.I0) Это ее свойство связано с тем, что низшие собственные числа /* » v > Р оператора 5С г удовлетворяют так называемым соотношениям Аарони: м : v : Р =1:2:3. Аарони установил эти со- отношения путем численного эксперимента на ЭВМ Ll7j. В нашем случае соотношения Аарони получаются простым подсчетом числа точек минимума и точек перевала потенциала у , дающих основной вклад в соответствующие вариационные оценки, полученные в гл.ІУ. Возможность замены диффузионной модели на ВДО основана на следующем основном предположении о разделении спектра оператора 5-е Вычисление собственные чисел, проведенное в гл.ІУ на основании вариационного принципа, показывает, что это соотношение выполняется VoXVubOU^). Рассмотрим наряду с процессом И 1т) соответствующий процесс nvf) с "укрупненными состояниями". Теорема \М . Пусть с\ W) - функция распределения процес-са МС+) где UA* )4:)- функция распределения процесса Ml"t) , удовлетворяющая уравнению Фоккера-Планка (В.2). Пусть, далее,~tU) -функция, задающая временной масштаб, такая, что Для всех начальных распределений U0 «ионного процесса К (А) , не сильно уклоняющихся от стационарного распределения процесс Hfo) является асимптотически марковским (при -* 0), и его функция распределения эволюционирует по закону ііууч ^Guo) = с ^ Со). Следующие главы посвящены получению результатов, на которых основана эта теорема. Во второй главе строится формальная погранслойная асимптотика младших собственных функций оператора * . Ввиду экспоненциальной малости соответствующих собственных чисел, дело сводится к построению решения уравнения (B.I2) - о г имеющего погранслойный характер вблизи стороны Z сферического треугольника Sla . Погранслойное решение уравнения изучается в трубчатой окрестности седловой точки потенциала ч » расположенной вдоль границы фундаментальной области -^М . Вводя на естественные координаты ( 5 , К ) (S =0в седловой точке) и переходя к новым переменным N * tfL Ь' > (В.ІЗ) получаем следующую задачу для выделения погранслоиного решения 7 іГІ * О tf-* 1 (B.I4) где с і 7 t Сґ - формальные ряды по степеням с полиномиальными по ^ коэффициентами, причем ряд в главном порядке нигде не обращается в нуль на контуре , а главный член ряда & обращается в нуль лишь в критических точках потенциала ч* ? = &(*)+ОСО *= сс*н осо , cw\ -о,C'(4)Uо. (ВЛ5) S-=.o $so (т = 1 л ОСО Погранслойное решение уравнения (В. 14) получается в виде Лемма 2,.1. Главный член асимптотики ^*(s, vO удовлетворяет нелинейному уравнению в частных производных Единственное нетривиальное полиномиальное решение этого уравнения, регулярное в окрестности седловой точки S =0, имеет вид где У Is) - специальное решение уравнения Бернулли & + &&>*- »о As CCsH с to (Б.І9) -\ -г а г г - л^А Здесь Rs) - ненулевое решение однородного уравнения F'(*W fcU)C"'w Fte) Старшие члены разложения |к(^^) удовлетворяют уже линейным уравнениям, разрешимость которых в полиномах гарантирует Лемма 2..2. . Старшие члены асимптотики ^k(S,v>) являются решениями рекуррентной системы неоднородных линейных дифференциальных уравнений Ч>Уг U ^V IS (В.20) Правые части т ус зависят от уже построенных решений ^е » ^ \С-1, ив предположении полиномиальности последних сами являются полиномами по V : степени не выше Я/К+І . Коэффициенты полиномиального решения удовлетворяют рекуррентной системе обыкновенных сингулярных дифференциальных уравнений сщ[Ам^-Ы\Щ^ \ (B.2D Требование регулярности решений ^\cte>W в окрестности перевальной точки S = 0 устраняет произвол при решении уравнений (В.20). В частности, следующий член разложения $1 -полином третьей степени по V где коэффициенты $Аг , ^4Л задаются явными формулами. Уравнение Бернулли (В.19) и все его высшие аналоги (В.21) имеют сингулярности в седловой точке потенциала S = 0 и в вершинах ( ^ = S> )- точках максимума потенциала *f на границе стороны , сферического треугольника Э_ ^ . При этом регулярные решения этих уравнений в окрестности перевальной точки S = 0 перестают быть регулярными в окрестности точек вершины S = S» . Вычисляя асимптотику главных членов в окрестности вершины, например * (*)» ІЇіА\ ~ fe (S-S*)^*-S) * OCS-S^ ) (в.22) получаем при К > 2 (В.23) что указывает на непригодность полученного разложения для описания пограничного слоя в \fiL - окрестности вершины <> = S В третьей главе строятся локальные асимптотические разложения, описывающие поведение пограничного слоя в окрестности вершин. Здесь существует два типа погранслоиных решений. К первому типу относится решение, описывающее сшивание двух пограничных слоев, приходящих в вершину, ко второму - решение, описывающее сшивание четырех пограничных слоев (угловой пограничный слой). Пусть Z и Ь - стороны сферического треугольника -И-4 , приходящие в вершину S = S . Тогда погранслойные решения первого и второго типов определяются соответственно как решения следующих задач е 'ъъ\* е г* Сделав масштабное преобразование где х , Ч- - декартовые локальные координаты на единичной сфере, получим уравнение I 1 Г—7 " &)[ ^ *^*я * ^ ^ = (в.*) Здесь функции іП» и /± - формальные ряды по степеням Є , с полиномиальными по \ и 3 коэффициентами. Решение, описы-ыввдее сшивание двух погранслоев, определяется как решение уравнения (В.26), удовлетворяющее условиям (В.27) (В.28) соответствуют угловому пограничному слою. Решение уравнения (В.26) с условиями (В.24) строим в виде: где ОХ - ряд по степеням Є и we с полиномиальными по 3 коэффициентами A = 21 ^ L t о«кЛС|л) - --- ^-3} Необходимость введения логарифмических членов по связана с наличием логарифмического слагаемого в главном члене ряда V в окрестности вершины, построенного в гл.П, и делает впоследствии возможным сшивание разложения (В.29) с разложением (B.I8). Подстановка (В.30) в уравнение (В.26) дает рекуррентную систему для определения членов ряда иг . Выпишем уравнения для главных членов ряда А^оо-1^оо^^0>-^о-^| +1341-^-1^. в с !TJ 2Ъг Здесь функции А {,1 определяются через потенциал Ч> следующим образом 4*=-xa,U^ , ^з=-^ги^) (В.32) Уравнения старших порядков имеют следующий вид М -^^+^^-<л=^С|Л) (в.зз) также разрешимы в полиномах по "J Доказана однозначная разрешимость полученных уравнений в главных порядках, исходя из вида граничных условий (В.27) и условий сшивания с разложением гл.П. При этом оказывается, что <*oZ = О при Z^> I. Если \}q*" - аналогичное решение, сосредоточенное в окре-стности стороны Ь то решение ЖоеЛ , описывающее угловой пограничный слой, ищем в виде + &E<-V(l)Er&) + ГЕЛ'(|)Е*)\ , (В"35) где ОТ и Jt> - уже описанные ряды с полиномиальными по "Т и "J коэффициентами соответственно, а ряд \ подлежит определению. Уравнения для коэффициентов ряда 1 определяются уже известными рядами <П и Л> и однозначно разрешимы в пространстве функций, представимых в виде суммы полиномов по ^ с коэффициентами, зависящими от 3 > и наоборот, если потребовать выполнимости граничных условий на ряд р : П = И = О и ограничить рост коэффициентов ряда \ при "|/5-» усло- r-.U*4**) (В.36) Из этого условия следует, что решение ЯГ^*- ^ ; 1^* (В.37) автоматически сшивается с разложениями 1)^ и lie* ЧР11 -* и "3 "^ соответственно, которые, в свою очередь, хорошо сшиты с разложениями гл.П. Следовательно, решение 1/^е*- описывает угловой пограничный слой. Вычислены главные члены раз- ложения для ряда і В четвертой главе проводится спектральный анализ диффузионного оператора Ссе . Используя вариационный принцип для определения его собственных значений ( 2С =<t В Lj^Wo ) \^ упік j \уц\г WoAx) )Ul\WeA%\ (B#38) мы доказываем наличие экспоненциальной серии собственных чисел \к-\ \(г))- > *ь(*-) и получаем оценку сверху для девятого собственного значения Теорема АЛ. Справедливы вариационные оценки для первых собственных чисел оператора Я.^ ( м = Х^ = Аэ = Ач , V = \* = ^ = \^_ , j> = Дх ) V ^ Q.J* ( і -V ОС*')) р і Ъ/а С і* ОСО) При доказательстве теоремы пробные функции конструируются из главных членов формальной погранслойнои асимптотики собственных функций оператора , изученной в главах П и Ш. Критерием близости пробных функций к младшим собственным функциям оператора 2. служит отношение где А* - вариационная оценка для восьмого собственного числа А X . С другой стороны, явный вид пробных функций позволяет доказать их близость в Ц wo -норме к кусочно постоян- ным функциям tffcl-fcu. (В.40) где d І хі - коэффициенты матрицы асимптотических значений собственных функций Я/ к в фундаментальных областях Я. ( , введенной в (В.7). Доказана следующая аппроксимационная теорема Теорема н Ъ . Кусочно-постоянные функции хорошо аппроксимируют младшие собственные функции оператора ^g. в том смысле, что 4'w- мл;' Здесь S" - любое положительное число. Полученная оценка становится эффективной, если спектр оператора S^ разделяется (В.II). Все дальнейшие теоремы основаны на предположении разделения спектра (В.II), и поэтому их следует считать условными. Аппроксимационная теорема позволяет изучить геометрию подпространства младших собственных функций. Оказывается, что это подпространство - "почти" инвариантно относительно оператора Теорема цМ . Пусть i^^Ux - стохастические состояния, введенные в гл.1, г - проектор на пространство первых восьми младших собственных функций оператора 2е » тогда матрица оператора СЗС, в базисе lW;\\^\ почти диагональна а пространство младших собственных функций "почти" инвариантно относительно оператора умножения tL : (здесь через ^С С обозначены точки минимума потенциала vp в фундаментальных областях Si I ). Теорема ^.Ч обобщается на случай оператора умножения на операторно-значную функцию во вспомогательном конечномерном пространстве. В последнем параграфе главы полученные результаты применяются для оценки погрешности в форме линии поглощения при замене диффузионной модели на МДО. В диффузионной модели форма линии W(.%) определяется билинейной формой резольвенты R^. оператора-зєД^'Н^ в точке Ч R^U^^-^)' (В.4І) вычисленной на элементах специального вида, которые, так же как и параметры Эс и "Z. , определяются физикой задачи. В ВДО форма линии Whao {.ъ ) определяется выражением, в котором диффузионный оператор - Ъ-У* g. заменяется на марковскую матрицу ъсН , а оператор умножения на операторно-значную функцию С5С в пространстве «- ^w0 - на оператор-функцию, принимающую конечное число значений, по одному в каждой фундаментальной области. Оценена погрешность при замене диффузионной модели на ЩО: \ \tf Ы - WM(Bg C%) \ < С (О г. (в.42) Таким образом, в предлагаемой работе дается последовательное оправдание перехода диффузионного процесса в симметричном поле сноса в марковский процесс с конечным числом состояний, основанное на асимптотике спектральных характеристик диффузионного оператора. В качестве приложений этих результатов дана оценка погрешности формы линии поглощения, полученной с помощью модели дискретных ориентации, тем самым определена область применимости этой модели в теории суперпарамагнетизма. В приложении вычислена по методу Лапласа асимптотика некоторых интегралов, возникающих при доказательстве вариационных оценок.
ОПЕРАТОРА "МАЛОЙ ДШУЗИИ" И ФОРМА ЛИНИИ
ПОГЛОЩЕНИЯ В ТЕОРИИ МЕССБАУЭРОВСКИХ
СПЕКТРОВ 106
венных значений 106
2. Асимптотика собственных функций
подпространства младших собственных
функций относительно оператора
умножения 123
Заключение 143Похожие диссертации на Математическое оправдание модели дискретных ориентаций
