Geum.ru
   Авторефераты по темам  >>  Разные специальности - [часть 1]  [часть 2]

Численное исследование влияния дальнодействующей корреляции дефектов на критическое поведение спиновых систем

Автореферат кандидатской диссертации

 

На правах рукописи

Колесников Вячеслав Юрьевич

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ВЛИЯНИЯ

ДАЛЬНОДЕИСТВУЮЩЕИ КОРРЕЛЯЦИИ ДЕФЕКТОВ НА

КРИТИЧЕСКОЕ ПОВЕДЕНИЕ СПИНОВЫХ СИСТЕМ

01.04.02 Ч теоретическая физика

Автореферат

диссертации на соискание ученой степени

кандидата физико-математических наук

Омск - 2009


Работа выполнена на кафедре теоретической физики ГОУ ВПО Омский государственный университет им. Ф.М. Достоевского.

Научный руководитель:аа доктор физико-математических наук,

профессор Прудников Владимир Васильевич.

Официальные оппоненты:аа доктор физико-математических наук,

профессор Аплеснин Сергей Степанович,

доктор физико-математических наук, ведущий научный сотрудник, профессор Щур Лев Николаевич.

Ведущая организация: Институт физики им. Х.И. Армиханова

ДагН - РАН, г. Махачкала.

Защита состоится л___ > декабря 2009 г. в____ часов на заседании диссертаци

онного совета Д 212.179.04 при ГОУ ВПО Омский государственный университет

им. Ф.М. Достоевскогопо адресу: 644077, г. Омск, пр. Мира, 55а.

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО Омский государственный университет им. Ф.М. Достоевского.

Автореферат разослан л___ ноября 2009 г.

Ученый секретарь

диссертационного совета Д 212.179.04

кандидат физико-математических наукаа Вершинин Г.А.


Общая характеристика работы Актуальность темы

Исследование критического поведения неупорядоченных систем представляет большой теоретический и экспериментальный интерес, поскольку большинство реальных твердых тел содержат замороженные дефекты структуры, присутствие которых влияет на характеристики систем. В частности, при фазовых переходах поведение таких систем может существенно меняться. В большинстве работ исследование ограничивается рассмотрением точечных дефектов низкой концентрации. Согласно критерию Харриса [1], такие дефекты изменяют критическое поведение только систем с расходящейся теплоемкостью (изингоподобные магнетики). В то же время вопрос о влиянии на критическое поведение эффектов корреляции дефектов практически не исследован. В рамках этой же проблемы можно поставить вопрос о влиянии на критическое поведение протяженных дефектов (дислокаций, границ зерен), что еще больше приближает исследователей к описанию реальных материалов.

В работе [2] представлена модель изотропной неупорядоченной системы с даль-нодействующей корреляцией дефектов. Был получен критерий существенности ее влияния на критическое поведение систем и показано, что дефекты, обладающие свойством дальней пространственной корреляции, могут при определенных условиях изменять критическое поведение не только систем с однокомпонентным параметром порядка (модель Изинга), но и систем с двухкомпонентным (XY-модель) и трехкомпонентным (модель Гейзенберга) параметром порядка. Также были получены значения статических критических индексов в однопетлевом приближении с использованием метода двойного , ?-разложения. Однако существует ряд работ [3, 4], посвященных теоретико-полевому описанию однородных и неупорядоченных моделей, которые показывают существенное различие реального критического поведения и предсказаний однопетлевых приближений (особенно результатов -разложения). В работе [5] было осуществлено теоретико-полевое описание критического поведения непосредственно трехмерных систем с дальнодействующей изотропной корреляцией дефектов в двухпетлевом приближении с последовательным применением для анализа рядов разложения методов суммирования и проведен расчет статических и динамического критических индексов для систем с различным числом компонент параметра порядка и различными значениями параметра корреляции. Было выявлено значительное отличие характеристик критического поведения систем с дальнодействующей корреляцией от аналогичных характеристик для однородных систем и систем с некоррелированными дефектами. Удалось установить, что получающаяся картина областей устойчивости различных типов критического поведения для одних и тех же параметров модели, существенно отличается от предсказываемых в работе [2].

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

3


ведение во всем диапазоне концентраций примеси вплоть до порога перколяции. Ренормгрупповое описание не дает ответа на этот вопрос, поскольку применимо лишь в области низких концентраций дефектов.

Цель работы

Целью настоящей диссертации является:

  1. исследование влияния дальнодействующей корреляции дефектов структуры на критическое поведение систем с различным числом компонент параметра порядка посредством численного изучения методами Монте-Карло ферромагнитных трехмерных моделей Изинга и XY.
  2. численное исследование неравновесного критического поведения трехмерных моделей Изинга и XY с линейными дефектами при спиновых концентрациях р = 0.80 и р = 0.60 методом коротковременной динамики при рассмотрении эволюции систем из разных начальных неравновесных состояний. Определение совокупности значений для независимых динамических z: 9 и статических z/, [3 критических индексов с применением метода поправок к скейлпнгу. Сопоставление полученных значений критических индексов для слабо неупорядоченных систем с р = 0.80 с существующими результатами теоретико-полевых расчетов.
  3. численное исследование равновесного критического поведения трехмерной модели Изинга с линейными дефектами со спиновой концентрацией р = 0.80 традиционным методом Монте-Карло и определение совокупности статических критических индексов с применением метода поправок к скейлпнгу. Сопоставление полученных значений критических индексов со значениями аналогичных критических индексов, определенных методом коротковременной динамики.

Научная новизна результатов

  1. Впервые осуществлено компьютерное моделирование неравновесного критического поведения трехмерных моделей Изинга и XY с дальней пространственной корреляцией дефектов в коротковременном режиме. При исследовании критической релаксации модели из различных начальных состояний системы определены значения совокупности динамических и статических критических индексов при применении методики учета поправок к скейлпнгу. Полученные результаты позволяют сделать вывод о существовании различных классов универсального критического поведения для рассматриваемых систем, отвечающих областям слабой и сильной структурной неупорядоченности.
  2. Впервые получено численное подтверждение о существенности влияния дальней пространственной корреляции дефектов на критическое поведение не только изингоподобных систем, как в случае систем с некоррелированным структурным беспорядком, но и систем с многокомпонентным параметром порядка (на примере XY-модели).
  3. Впервые продемонстрировано при сопоставлении результатов компьютерного моделирования неравновесного критического поведения трехмерной модели Изинга с дальней пространственной корреляцией дефектов в коротковременном режиме и ее равновесного критического поведения, что метод коротковременной динамики может служить надежной альтернативой традиционным методам Монте-Карло не только при численных исследованиях однородных систем, но и систем со структур-

4


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

Научная и практическая значимость работы

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

Исследование влияния дефектов структуры и эффектов их корреляции является актуальным направлением современной физики конденсированного состояния, т.к. практически все реальные материалы содержат примеси и другие дефекты структуры. Дальнодействующая корреляция в пространственном распределении дефектов может модифицировать критические свойства неупорядоченных систем. На это указывают эксперименты по рассеянию нейтронов и рентгеновского излучения на различных системах, находящихся в критических точках. В силу этого к моделям систем с дальнодействующей корреляцией дефектов существует несомненный интерес как с общетеоретической точки зрения выявления новых типов критического поведения в неупорядоченных системах, так и с точки зрения реальной возможности проявления дальнодействующей корреляции дефектов в ориента-ционных стеклах [6], полимерах [7] и неупорядоченных твердых телах с дефектами фракталоподобного типа [8].

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

Основные положения, выносимые на защиту

  1. Методика численного исследования неравновесного критического поведения структурно неупорядоченных трехмерных моделей Изинга и XY с дальней пространственной корреляцией дефектов в коротковременном режиме и методика определения значений критических индексов с учетом ведущих поправок к скейлингу.
  2. Наличие нескольких этапов динамического развития слабо неупорядоченных систем после микроскопического временного масштаба: области с характеристиками однородной системы, кроссоверной области и области, характеризующейся влиянием структурного беспорядка.
  3. Подтверждение расширенного критерия Харриса о влиянии дефектов с дальней пространственной корреляцией на критическое поведение не только изингопо-добных систем, но и систем с многокомпонентным параметром порядка (на примере XY-модели).
  4. Возникновение при концентрациях спинов большей порога спиновой перколя-

5


ции двух классов универсального критического поведения, отвечающих слабой и сильной структурной неупорядоченности.

Апробация работы

Основные результаты диссертационной работы докладывались и обсуждались на III Международной конференции Фундаментальные проблемы физики(Казань, 2005), The 3-rd International Workshop Hangzhou 2006 on Simulational Physics (Hangzhou, 2006), Семинаре по вычислительным технологиям в естественных науках (Таруса, 2009), Международной конференции Фазовые переходы, критические и нелинейные явления в конденсированных средах (Махачкала, 2009), а также на научных семинарах кафедры теоретической физики ОмГУ.

Публикации

Список публикаций автора по теме диссертации включает 10 статей и тезисов докладов, опубликованных в российских и иностранных журналах, сборниках трудов и материалах конференций.

Структура и объем диссертации

Диссертация состоит из введения, четырех глав, и заключения. Объем диссертации - 108 страниц машинописного текста, в том числе 28 рисунков, 13 таблиц, и список цитируемой литературы из 108 наименований.

Краткое содержание работы

Во введении обоснована актуальность выбранной темы диссертационной работы и сформулированы основные цели исследований.

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

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

Рассматривается модель неупорядоченной спиновой системы в виде кубической решетки с линейным размером Lи наложенными граничными условиями. Микроскопический гамильтониан неупорядоченной модели Изинга записывается в виде

н = Ч^23з№Рз(!)

где Jij- короткодействующее обменное взаимодействие между закрепленными в узлах решетки спинами <jj, принимающими значения 1. Немагнитные атомы примеси образуют пустые узлы. Числа заполнения piпри этом принимают значения 0 или 1 и описывают распределение замороженных дефектов структуры с дальнодействующей изотропной пространственной корреляцией д(х Ч у) ~ \хЧу\~а7 описываемой моделью Вейнриба-Гальперина [2] с параметром корреляции а. В диссертации

6


была осуществлена численная реализация этой модели в виде системы линейных дефектов, характеризуемых а = 2. Для этого в заполненной спинами трехмерной решетке "вычеркивались" линии, параллельные осям координат, до достижения заданной концентрации примесей Cimp= 0.2. Для обеспечения изотропности в распределении линейных дефектов в решетке удалялось одинаковое число линий в каждом направлении при случайном положении линий в координатных плоскостях. Дополнительно накладывалось условие непересекаемости этих линейных дефектов, что позволяло гарантировать существование в системе единого протекающего спинового кластера. Физическое обоснование этого условия связано с тем фактом, что в реальных кристаллах дислокации, как линейные дефекты, распределяются в образце равномерно с вероятностью их пересечения, близкой к нулю. Условие непересекаемости линейных дефектов соответствует приближениям модели Вейнриба-Гальперина, поскольку их пересечение привело бы к дополнительным вершинам взаимодействия, которые отсутствуют в гамильтониане модели [2].

Для определения критической температуры было осуществлено компьютерное моделирование системы в состоянии равновесия при различных температурах. Для снижения влияния эффектов критического замедления и корреляции различных спиновых конфигураций был применен наиболее эффективный в этом смысле од-нокластерный алгоритм Вольфа. За один шаг Монте-Карло на спин (MCS/s) принималось 5 переворотов кластера Вольфа. Процедуре установления термодинамического равновесия в системе отводилось 104 MCS/s и 105 MCS/s отводилось на статистическое усреднение вычисляемых характеристик системы при заданной примесной конфигурации. Для определения средних значений термодинамических и корреляционных функций наряду со статистическим усреднением применялось усреднение по различным примесным конфигурациям: по 15000 образцам для решеток с линейными размерами L= 16 и 32, а для решеток с L= 64 и 128 - по 10000 образцам. Анализировались зависимости кумулянта Биндера 4-го порядка U{T,L) и отношения г/L(T,L) корреляционной длины ? к Lот температуры

где S= $><*, F= [(< Ф >)]/pL\ Ф = \ Е

- координаты -го узла трехмерной решетки. Скобки (...) означают статистическое усреднение по шагам Монте-Карло, а [...] - усреднение по различным реализациям распределения дефектов структуры в системе при заданной спиновой концентрации р. По точкам пересечения графиков данных величин для различных размеров решетки (рис. 1) была определена критическая температура Тс = 3.9281(1).

При найденной критической температуре было осуществлено численное исследование неравновесной критической динамики в коротковременном режиме для трехмерной слабо неупорядоченной модели Изинга (со спиновой концентрацией

7

3,9282

3,928

и

0,65-

(а)

0.54 0^

^

^V

0,60-

0,535.

\3,9265

3,9270аа 3,92

0,55.

0,50.

X.аа 64

^128 жаа жаа 1аа >

3,922а 3,924а 3,926а 3,928а 3,930 J


Рис. 1: Зависимости кумулянта Биндера 4-го порядка U{T,L) (а) и отношения г/L(T,L) (b) от температуры для различных размеров решетки

р = 0.8) с линейными дефектами. Было осуществлено моделирование критической релаксации системы из полностью упорядоченного начального состояния (с начальной намагниченностью то = 1), а также критической эволюции системы из начальных неупорядоченных состояний с то <С 1.

Для моделирования спиновых конфигураций в системе был применен алгоритм Метрополиса. В работе была исследована релаксационная динамика системы, описываемая моделью А в классификации моделей критической динамики, проведенной Хоэнбергом и Гальпериным [9]. Алгоритм Метрополиса, реализующий динамику односпиновых переворотов, наилучшим образом соответствует релаксационной модели А и позволяет провести сравнение получаемого в результате моделирования критической релаксации системы динамического критического индекса zс результатами теоретико-полевого описания [5] критической динамики модели А для трехмерных систем с дальнодействующей пространственной корреляцией дефектов структуры.

Традиционное моделирование критического поведения системы взаимодействующих частиц методом Монте-Карло наталкивается на трудности, связанные в основном с явлением критического замедления, характеризующимся тем, что время релаксации системы, как и время корреляции состояний, неограниченно растут по мере приближения к критической температуре и степенной характер их асимптотической зависимости от приведенной температуры определяется динамическим индексом z. Для структурно неупорядоченных систем эта проблема еще более существенна, т. к. их неравновесное критическое поведение определяется индексом z, принимающим большие значения, чем для однородных систем. Для уменьшения эффектов влияния критического замедления применяют кластерные алгоритмы Вольфа или Свендсена-Ванга, но эти алгоритмы столь существенно меняют динамику системы по сравнению с алгоритмом Метрополиса, что для получения информации о характеристиках критической динамики их применять нельзя. В связи с этим в работе был применен метод коротковременной динамики (МКД) для получения значений как динамического, так и статических критических индексов. Особенностью МКД является то, что информация об универсальном критическом поведении может быть получена на относительно малых макроскопических промежутках времени (от 1000 до 2000 MCS/s) на ранней стадии развития системы в критической точке или ее окрестности.

МКД был обоснован результатами аналитических и численных исследований,


проведенных в работах [10, 11]. Так, в работе [10] на основе ренормгруппового анализа было показано, что после микроскопически малого промежутка времени tm\cдля к-то момента намагниченности системы реализуется обобщенно-однородная (скейлинговая) форма


M^(t,T,L,m0) = b-kiilvM^ (ь-Ч^т^ЦЬ^то


(4)


где t- время, г = (Т Ч Тс)/Тс - приведенная температура, L- инейный размер решетки, Ъ - произвольный масштабный фактор, /3, v:z- известные критические индексы, хо - новый независимый критический индекс, задающий скейлинговую размерность начального значения намагниченности т$.

Для неупорядоченных систем вычисление осуществляется в виде

M^(t)


(5)


В работе проводилось усреднение вычисляемых величин по 3000 различным примесным конфигурациям.

Для начального состояния системы с то = 1 (спины ориентированы в одном направлении) и для решеток с достаточно большими линейными размерами Lуравнение (4) для намагниченности (к = 1) при выборе фактора Ъ = tl'zпринимает следующий вид:


M(t,T) = t~fj/vzM (l,t1/l/zr\ ~ t~fj/vz (l+ At1/vzT+ 0(t2)) Х


(6)


В пределе г Ч>ж 0

M(t) ~ t~Plvz.(7)

Другой определяемой в работе величиной является кумулянт Биндера, характеризуемый выражением

U(t)

1.

(8)

Mt2\t) (M(t)f

Конечномерный скейлинговый анализ показывает, что в критической точке поведение кумулянта Биндера описывается степенным законом:


U{t) ~ t


d/z


О)


где d- размерность системы.

Продифференцировав обе части уравнения (6) по т, при г = 0 получим еще одно соотношение (в дальнейшем называемое логарифмической производной намагниченности), описываемое степенной зависимостью от времени в критической

точке:

l/vz

д

(10)

дт

т=0

\nM{t,r)а ~t

9


Рис. 2: Усредненные значения намагниченности (а) и ее логарифмической производной (Ь) и кумулянта Биндера (с) на временном интервале до 1000 MCS/s в двойном логарифмическом масштабе

В работе осуществлялось моделирование кубических решеток с размером L= 128 при критической температуре Тс = 3.9281. Временное поведение намагниченности и кумулянта Биндера исследовалось на временах до 1000 MCS/s. С целью численного определения логарифмической производной намагниченности проводился расчет намагниченности для двух близких к Тс температур, т. е. для Т_ = 3.9250 и Т+ = 3.9310. На рис. 2 приведены итоговые зависимости намагниченности, кумулянта Биндера и логарифмической производной намагниченности от времени в двойном логарифмическом масштабе.

Анализ зависимости кумулянта U(t) показал, что во временном интервале [50,150] MCS/s, степенному характеру зависимости U(t) соответствует значение динамического индекса z~ 2.02, описывающее критическое поведение однородной модели Изинга [12], а влияние линейных дефектов начинает проявляться лишь на временах t> 400 MCS/s. Выявленные динамические кроссоверные явления учитывались также при анализе временных зависимостей намагниченности и ее логарифмической производной.

В диссертации был осуществлен учет поправок к асимптотической зависимости измеряемых величин за счет влияния конечности моделируемых систем, т. к. только учет данных поправок к скейлингу позволяет получать корректные значения критических индексов в термодинамическом пределе LЧ> сю. Для этого было использовано следующее выражение для временной зависимости наблюдаемых величин X(t):

X(t)~ts(l + Axt-W'z},(11)

где Ах - неуниверсальные амплитуды, и является хорошо известным критическим индексом поправки к скейлингу, а показатель 8 = Ч(З/vzв случае X= M(t)) 8 = d/zв случае X= U(t) и 8 = l/uzв случае X= dT\nM(t). Теоретико-полевая оценка в двухпетлевом приближении для показателя и модели Изинга дает значение и = 0.8. При анализе полученных кривых показатель и варьировался в интервале от 0.7 до 1.0. Для расчета значений критических индексов fijvz,

Таблица 1: Значения показателей /3/isz, ^jvz^ d/zи минимальные значения среднеквадратичных погрешностей а аппроксимации при различных значениях критического индекса ио

UJ

(3/uz

а

l/uz

а

d/z

а

0.7

0.2112

0.0100

0.556

0.0053

1.183

0.0100

0.8

0.2096

0.0088

0.559

0.0049

1.205

0.0100

0.9

0.2101

0.0093

0.553

0.0070

1.213

0.0099

1.0

0.2090

0.0095

0.558

0.0072

1.227

0.0098

10


d/z: 1/vzи cu/zна временном интервале, соответствующем влиянию структурного беспорядка, был использован метод наименьших квадратов для наилучшей аппроксимации значений M(t): U(t) и dT\nM(t) выражением (11). По минимумам зависимостей погрешности аппроксимации а от 6 и cu/zопределялись значения соответствующих критических показателей. В таблице 1 приведены полученные итоговые значения критических показателей, соответствующие минимальным погрешностям процедуры аппроксимации (11), для различных значений показателя си.

Видно, что значения этих показателей слабо зависят от изменения си в рассмотренном интервале, но значение си = 0.8 оказывается более предпочтительным, т.к. обеспечивает лучшую аппроксимацию полученных данных. На основе значений показателей были определены динамический критический индекс z= 2.489(21), статические критические индексы v= 0.719(22), [3 = 0.375(45). Сравнение этих значений критических индексов со значениями, вычисленными в работе [5] в рамках теоретико-полевого подхода, z= 2.495, v= 0.716, [3 = 0.350, показывает их хорошее согласие в пределах статистических погрешностей моделирования и аппроксимации.

Далее исследовалась критическая эволюция системы из начальных неупорядоченных состояний с то <С 1. Согласно теории МКД, для этого режима можно получить следующие соотношения для временных зависимостей намагниченности M(t), ее второго момента M^2>{t) и автокорреляционной функции A(t):

M(t)~te,e= (x0-г)-;(12)

MW(t)~t<*,C2 = (d-2г)-;(13)

Isа /С

А()~Гс%са = --в.(14)

/С

Используя данные зависимости, были определены показатели 0, С2 и са: а на их основе вычислялись критические индексы (3/v:z7xq. Чтобы вычислить значения критических индексов, для слабо неупорядоченной модели Изинга было реализовано компьютерное моделирование решетки с линейным размером L= 128 и спиновой концентрацией р = 0.80 при двух значениях начальной намагниченности то = 0.02 и то = 0.001 с последующей линейной аппроксимацией результатов к то = 0. Для одной примесной конфигурации проводилось 25 "прогонок" по 700 MCS/s при температуре Тс = 3.9281. Средние значения вычисляемых термодинамических величин получались усреднением по 3000 различным примесным конфигурациям. На рис. 3-4 приведены графики временных зависимостей исследуемых величин в двойном логарифмическом масштабе, что позволяет по наклону линейных участков графиков определять соответствующие показатели.

Из рисунков видно, что на каждом графике могут быть выделены по два линейных участка: для временных интервалов от 10 до 70 MCS/s и от 100 до 650 MCS/s. Это связано с наблюдаемым явлением кроссовера, т. е. перехода от поведения, характерного для однородной системы, к поведению неупорядоченной системы с линейными дефектами. Для каждого линейного участка исследуемых величин были

11

{ 100С

t


Рис. 3: Эволюция намагниченности в двойном логарифмическом масштабе из начальных состояний с т0 = 0.02 (а) и т0 = 0.001 (Ь)

Таблица 2: Значения критических показателей для начального участка эволюции от 10 до 70 MCS/s в сравнении с показателями однородной системы, а также значения для участка от 100 до 650 MCS/s, соответствующие слабо неупорядоченной модели

участки

т0

в

с2

са

t = 10 ч- 70 MCS/s

0.02

0.086(12)

0.964(28)

1.384(26)

0.001

0.099(9)

0.973(19)

1.364(23)

0

0.101(10)

0.975(23)

1.363(26)

однородная, то = 0 [13]

0.108(2)

0.970(11)

1.362(19)

t = 100 -т- 650 MCS/s

0.02

0.152(12)

0.812(21)

1.103(16)

0.001

0.149(10)

0.804(19)

1.047(12)

0

0.149(11)

0.801(20)

1.043(14)

определены степенные показатели при то = 0.02 и то = 0.001 с последующей их линейной аппроксимацией к то = 0. Для выделенных на графиках линейных участков, соответствующих поведению однородной системы, было проведено сравнение получаемых показателей с результатами работы [13]. Из табл. 2 видно, что значения наших показателей сопоставимы с соответствующими показателями однородной системы [13]. Для линейных участков, соответствующих критическому поведению системы с дефектами, также был осуществлен расчет соответствующих показателей, которые представлены в табл. 2.

При использовании соотношений, связывающих показатели 0, с^ и са с критическими индексами, для модели Изинга с линейными дефектами и со спиновой концентрацией р = 0.80 были определены значения (З/v= 0.492(28), z= 2.517(32), xq= 0.867(37). Сопоставление со значениями соответствующих критических индексов Р/ту = 0.489 и z= 2.495, полученными в [5] с применением методов ренорм-группового описания, показывает их хорошее согласие в пределах статистических погрешностей численных исследований.

С целью проверки применимости метода коротковременной динамики для численного описания структурно неупорядоченных систем и полученных с его помо-

(с)

^

т0=о.02 т0 =0.001

А "

,-'ж'

Щ,t

Рис. 4: Эволюция второго момента намагниченности (а, Ь) и автокорреляционной функции (с) в двойном логарифмическом масштабе из начальных состояний с toq = 0.02 и toq = 0.001

12


Рис. 5: Зависимости среднеквадратичной погрешности аппроксимации а кривых теплоемкости (а), намагниченности (Ь), восприимчивости (с) и производной кумулянта по температуре (d) как функций показателей а/и, (З/и, ^ jvи vдля различных значений ио

Таблица 3: Значения критических индексов, полученные при среднем значении uj= 0.76, а также значения, полученные с помощью МКД и теоретико-полевого описания

а

Р

7

V

моделирование в равновесии

-0.078(30)

0.362(20)

1.441(15)

0.710(10)

МКД

-0.157(66)

0.375(45)

1.41(11)

0.719(22)

теоретико-полевое описание [5]

-0.1048

0.3504

1.4453

0.7155

щью результатов в диссертации было реализовано также компьютерное моделирование трехмерной модели Изинга с линейными дефектами традиционными методами Монте-Карло в состоянии равновесия. Для моделирования использовался одно-кластерный алгоритм Вольфа. В состоянии равновесия при критической температуре Тс = 3.9281 был проведен расчет различных термодинамических и корреляционных функций, таких как намагниченность, восприимчивость, корреляционная длина, теплоемкость и кумулянт Биндера U^, для решеток с линейными размерами L= 16-1-128 и спиновой концентрацией р = 0.80. Использование хорошо известных скейлинговых соотношений для рассматриваемых функций с учетом конечнораз-мерных скейлинговых поправок C(L) ~ Lalv{\ + aL^), M(L) ~ L~^v{\ + ЬЬ~Ш), x(L) ~ L7^(l + cL~UJ)) ^ ' г^ Ll^(l+ (1Ь~и) позволяет определить критические индексы а) (3) 7? vи и с помощью статистической обработки данных моделирования. Для этого проводилась линейная аппроксимация зависимостей (XL~A) от L_CJ, где X - исследуемая функция, А - соответствующий критический показатель. Затем исследовались графики зависимостей погрешности аппроксимации а функций XL_A(L_a;) при варьировании значений параметров А и ш. На рис. 5 представлены зависимости а для теплоемкости, намагниченности, восприимчивости и производной кумулянта по температуре как функции соответствующих показателей olJv, (3/v, j/vи vДля различных значений ио.

По минимуму а определялись значения показателей. Минимум погрешности аппроксимации наблюдался при и = 0.90 для olJv, uj= 0.65 для (З/v)uj= 0.70 для 7/^5 oj= 0.80 для г/. Было взято среднее значение и = 0.76(5), при котором и были рассчитаны окончательные значения показателей. Для сравнения в табл. 3 приведены значения критических индексов, полученных с помощью компьютерного моделирования системы в равновесном состоянии, а также с помощью МКД и теоретико-полевых расчетов [5]. Полученные значения находятся в хорошем согласии в пределах погрешностей измерений с результатами применения метода корот-ковременной динамики и теоретико-полевого описания.

13


Рис. 6: Зависимости кумулянтов Биндера С/4 от температуры для слабо неупорядоченной XY-модели (а), сильно неупорядоченных моделей Изинга (Ь) и XY (с) для решеток различных размеров L

В третьей главе осуществлено компьютерное моделирование слабо неупорядоченной трехмерной XY-модели с дальнодействующей корреляцией дефектов. С помощью метода кумулянтов Биндера определена критическая температура. Значения критических индексов получены с помощью МКД.

Трехмерная классическая XY-модель описывает критическое поведение многих физических систем: прежде всего, широкий класс сильно анизотропных магнетиков (с анизотропией типа "легкая плоскость"), а также сверхтекучего гелия АНе и сверхпроводников. Эта модель характеризуется гамильтонианом вида:

Н = -К ^2 Vi^jPiPji(15)

(hi)

где &i= (erf, (jf) - плоский единичный вектор в узле г, сумма берется по всем ближайшим узлам решетки, pi- случайные переменные, задающие распределение в решетке замороженных дефектов структуры с дальнодействующей пространственной корреляцией. Распределение линейных дефектов задается так же, как и для модели Изинга (см. гл. 2). Параметром порядка модели является средний модуль намагниченности:

-

М = sjMj + MlмхЬ = Ч Y, AvVi-(16)

Расчет критической температуры для трехмерной XY-системы с линейными дефектами со спиновой концентрацией р = 0.80 был проведен по методу кумулянтов Биндера. Использовался однокластерный алгоритм моделирования Вольфа. Для размеров решетки L= 32, 64,128 усреднение осуществлялось, соответственно, по 2200, 2200, и 1830 примесным конфигурациям. Для каждой примесной конфигурации нужные величины усреднялись по 1000 MCS/s (за один шаг - три переворота кластера). Для достижения равновесного состояния, проводилась термализация, на которую отводилось 200 MCS/s. На рис. 6(a) приведены графики кумулянтов (2) для различных размеров решетки в диапазоне температур [1.861; 1.8645] с шагом 0.00039. Из пересечения данных кривых было получено значение критической температуры: Тс = 1.8626(5).

Затем при рассчитанной критической температуре было проведено моделирование трехмерной XY-системы с линейными дефектами со спиновой концентрацией р = 0.80 методом коротковременной динамики. Для меньшего влияния конечности размера решетки моделируемой системы был взят максимально возможный с точки зрения временных затрат на вычисления размер решетки L= 128.

14


При исследовании критической релаксации системы для получения зависимостей намагниченности (7), ее логарифмической производной по температуре (10) и кумулянта Биндера 2-го порядка (9) от времени проводилось усреднение по 830 примесным конфигурациям и по 5 прогонкам для каждой примесной конфигурации. Рассматривался временной интервал t= l-i-1000 MCS/s. Для вычисления логарифмической производной намагниченности были дополнительно получены данные по динамике намагниченности для двух температур, близких к критической: Т_ = 1.8502 и Т+ = 1.8749. Усреднение данных при Т проводилось по 110 примесным конфигурациям и по 5 прогонкам для каждой конфигурации. На рис. 7 представлены временные зависимости исследуемых величин. Анализ данных зависимостей с учетом процедуры применения поправки к скейлингу позволил определить следующие значения для критических показателей (З/vz= 0.221(2) 5 d/z= 1.269(3) и \jvz = 0.55(4).

При исследовании критической эволюции системы из начальных неупорядоченных состояний с малой или нулевой начальной намагниченностью ttiqставилось целью определение значений динамических критических индексов в и z, а также отношения статических критических индексов (З/v. Для расчета индекса в были получены зависимости намагниченности (12) от времени для различных малых значений начальной намагниченности то = 0.01, то = 0.0075 и то = 0.005. Затем на основе полученных значений степенных показателей для различных то определялся индекс 0 5 как асимптотическое значение в пределе то Ч> 0. Для случая с то = 0.01 усреднение проводилось по 400 примесным конфигурациям и по 3 прогонкам для каждой примесной конфигурации; для для то = 0.0075 - по 300(3), для то = 0.005 - по 380(3) конфигурациям. При этом для каждой примесной конфигурации при реализации прогонок генерировались различные спиновые конфигурации, соответствующие заданному то. Для получения степенного поведения автокорреляционной функции (14) и второго момента намагниченности (13) расчеты проводились при то = 0.000001. Усреднение данных осуществлялось по 270 примесным конфигурациям и по 3 прогонкам для каждой примесной конфигурации. Исследование эволюции системы проводилось на временном интервале в 700 MCS/s. На рис. 8 показаны полученные зависимости рассматриваемых величин. На основе их анализа были получены значения критических показателей в = 0.374(14), d/z= 1.281(33) и (3/v= 0.534(35).

Рис. 7: Усредненные значения намагниченности (а), ее логарифмической производной (Ь) и кумулянта Биндера (с) для слабо неупорядоченной XY-модели на временном интервале t= 1 Ч 1000 MCS/s в двойном логарифмическом масштабе

15

В табл. 4 представлены итоговые результаты, полученные с помощью МКД для XY-модели с линейными дефектами для спиновой концентрации р = 0.80, а так-


Рис. 8: Эволюция намагниченности (а) из начальных состояний с то = 0.01, то = 0.0075 и то = 0.005, второго момента намагниченности (Ь) и автокорреляционной функции (с) из начального состояния с то = Ю-6 в двойном логарифмическом масштабе

Таблица 4: Значения критических показателей для слабо неупорядоченной XY-модели с линейными дефектами, полученные из МКД, в сравнении с результатами теоретико-полевых вычислений

z

Р

и

в

диссертация

2.358(25)

0.408(55)

0.78(6)

0.374(14)

работа [5]

2.365

0.37

0.76

же приведены значения критических индексов, полученные в [5] при теоретико-полевом описании модели. Видно хорошее согласие значений в пределах погрешностей.

В четвертой главе представлены результаты компьютерного моделирования трехмерных моделей Изинга и XY с линейными дефектами в области их сильной неупорядоченности (со спиновой концентрацией р = 0.60). Ставилось целью проверить, проявляется ли в критическом поведении систем с дальней пространственной корреляцией дефектов зависимость от их концентрации. В моделировании был также применен метод коротковременной динамики, хорошо зарекомендовавший себя при исследовании слабо неупорядоченных систем.

Для получения значений критических температур для сильно неупорядоченных моделей Изинга и XY было осуществлено вычисление равновесных значений кумулянта Биндера U4 (2). Из пересечения температурных зависимостей для решеток с размерами L= 64, 96,128 (рис. 6(Ь,с)) были получены значения критических температур Тс = 3.1956(34) для модели Изинга, иГс = 1.4455(5) для XY-модели. В этом случае для моделирования применялся однокластерный алгоритм Вольфа (25 переворотов кластера за MCS/s). Было использовано 5000 MCS/s для термализа-ции и 50000 MCS/s для статистического усреднения по спиновым конфигурациям. Окончательные результаты были получены усреднением по 500 образцам с различными конфигурациями линейных дефектов.

Далее с помощью МКД была исследована релаксация модели Изинга с линейными дефектами и спиновой концентрацией р = 0.60 из полностью упорядоченного начального состояния с ttiq= 1. Анализировались временные зависимости намагниченности (7), ее логарифмической производной по температуре (10) и кумулянта Биндера 2-го порядка (9). Результирующие кривые на временном интервале t= 1 -г- 2000 MCS/s были получены усреднением по 80 образцам с различными конфигурациями распределения линейных дефектов в решетке и по 25 прогонкам для каждого образца. Зависимость логарифмической производной намагниченности при критической температуре была получена на основе разностной схемы по двум кривым M(t) при температурах Т = 3.1956 и Т = 3.2180. На рис. 9 пред-

16


JVI(tJ)


(a)


Sin M(t,T)|t=0


U(t,T)


(c)



100 t, MCs/s 1000


100 t, MCs/sаа 100


100 t, MCs/sа 100


Рис. 9: Усредненные значения намагниченности (а), ее логарифмической производной (Ь) и кумулянта Биндера (с) для сильно неупорядоченной модели Изинга в двойном логарифмическом масштабе

M(Z)(t), хЮ

0,010-|

0,008

0,007

0,005 -I


m =0.002


t, MCs/s


A(t)



100аа t, MCs/s


100а t, MCs/sа 100


Рис. 10: Эволюция намагниченности из начальных состояний с различными то (а), второго момента намагниченности (Ь) и автокорреляционной функции (с) из начального состояния с то = 0 в двойном логарифмическом масштабе для сильно неупорядоченной модели Изинга

ставлены полученные зависимости в двойном логарифмическом масштабе. Анализ данных зависимостей с применением процедуры учета поправок к скейлингу позволил рассчитать значения критических показателей: (З// = 0.396(52) , z= 2.707(34) , /5 = 0.420(97) и и = 1.061(155).

При исследовании эволюции системы из неупорядоченных состояний с то <С 1 осуществлялся анализ временных зависимостей следующих величин: намагниченности M(t) (12), второго момента намагниченности M^2\t) (13) и автокорреляционной функции A(t) (14). Используя данные зависимости, в диссертации определялись показатели ^, С2 и са, а на их основе вычислялись критические индексы Р/ту и z. При расчете критического индекса в для сильно неупорядоченной модели Изинга было реализовано компьютерное моделирование решетки с линейным размером L= 128 и спиновой концентрацией р = 0.60 при двух значениях начальной намагниченности то = 0.002 и то = 0.01 с последующей линейной аппроксимацией результатов к то = 0. Для расчета значений с^ и са осуществлялось исследование эволюции системы из полностью неупорядоченных состояний с то = 0. Временные зависимости измеряемых величин изучались на временном интервале t= 1 -г-1000 MCS/s при критической температуре Тс = 3.1956. Данные усреднялись по 30 различным примесным конфигурациям при 25 "прогонках" для каждой примесной конфигурации. На рис. 10 представлены итоговые зависимости измеряемых величин в двойном логарифмическом масштабе. Их анализ позволил определить значения показателей в = 0.210(19), с^ = 0.808(27) и са = 0.902(32), а на основе соотношений, связывающих показатели с критическими индексами, были определены значения (З// = 0.409(52), z= 2.699(90) для модели Изинга с линейными дефектами при р = 0.60.

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

17


 

sin M(t,x;

100 t, MCs/s 1000

10а ioot)MCs/sаа 10o


t,MCs/sаа 100


Рис. 11: Усредненные значения намагниченности (а), ее логарифмической производной (Ь) и кумулянта Биндера (с) для сильно неупорядоченной XY-модели в двойном логарифмическом масштабе

ностей. Итоговые усредненные значения критических индексов представлены в табл. 5.

Аналогичные исследования неравновесной критической динамики были осуществлены для сильно неупорядоченной XY-модели с линейными дефектами и со спиновой концентрацией р = 0.60. На первом этапе осуществлялось моделирование релаксации системы из полностью упорядоченного начального состояния с т$ = 1. Результирующие кривые временной зависимости намагниченности (7) и кумулянта Биндера 2-го порядка (9) при критической температуре Тс = 1.4455(5) на временном интервале t= 1 -г- 3000 MCS/s были получены усреднением по 200 образцам с различными конфигурациями распределения линейных дефектов в решетке и 25 прогонкам для каждого образца. Зависимость логарифмической производной намагниченности (10) при критической температуре была получена на основе разностной схемы по двум кривым M(t) при температурах Тс = 1.4455 и Т+ = 1.4600 на временном интервале t= 1ч-2000 MCS/s. На рис. 11 представлены полученные зависимости в двойном логарифмическом масштабе. Итоговые значения критических показателей: (З/г/ = 0.51(1), z= 2.589(52), (3 = 0.511(72) и v= 1.001(143) были получены при анализе временных зависимостей измеряемых величин с применением к ним процедуры учета поправок к скейлингу.

При исследовании эволюции сильно неупорядоченной ХУ-модели из начальных состояний с то <С 1 осуществлялся анализ временных зависимостей следующих величин: M(t) (12), M^2\t) (13) и A(t) (14) при критической температуре Тс = 1.4455. Для определения критического индекса начальной эволюции в было реализовано компьютерное моделирование решетки с линейным размером L= 128 при трех значениях начальной намагниченности то = 0.004, ttiq= 0.003 и rriQ= 0.002 (рис. 12(a)). Для каждой зависимости M(t) определялся показатель 0(то) с последующей линейной аппроксимацией 0(т Ч> 0). Для расчета показателей С2 и са на основе зависимостей для второго момента намагниченности М^2> (t) и автокорреляционной функции A(t) проводилось моделирование системы из начального состояния с то = 0. Временные зависимости получаемых величин изучались на временном интервале t= 1 -г- 2000 MCS/s. Данные усреднялись по 1300 различным примесным конфигурациям при 25 "прогонках" для каждой примесной

Таблица 5: Значения критических показателей для сильно неупорядоченной модели Изинга с линейными дефектами, полученные из МКД

Zav

(P/v)av

Р

и

в

2.703(48)

0.403(38)

0.420(97)

1.061(155)

0.210(19)

18


Рис. 12: Эволюция намагниченности из начальных состояний с различной то (а), второго момента намагниченности (Ь) и автокорреляционной функции (с) из начального состояния с то = Ю-6 в двойном логарифмическом масштабе для сильно неупорядоченной XY-модели

Таблица 6: Значения критических показателей для сильно неупорядоченной XY-модели с линейными дефектами, полученные из МКД

Zav

(P/v)av

0

V

в

2.550(54)

0.537(25)

0.511(72)

1.001(143)

0.393(45)

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

Анализ данных зависимостей позволил определить значения показателей в = 0.393(45), С2 = 0.745(26) и са = 0.8016(26), а затем значения критических индексов P/i/ = 0.564(48) и z= 2.511(95). Сопоставление значений критических индексов, характеризующих неравновесное критическое поведение сильно неупорядоченной XY-модели из ее различных начальных состояний, показывает их согласие в пределах погрешностей. Итоговые усредненные значения критических индексов представлены в табл. 6.

Полученные результаты показывают, что наличие высокой концентрации дефектов с дальнодеиствующеи корреляцией оказывает сильное влияние на критическое поведение как трехмерной модели Изинга, так и трехмерной XY-модели. Анализ временных зависимостей исследуемых величин не выявил при этом временных режимов, соответствующих критическому поведению однородной системы, как это было в случае неравновесного критического поведения слабо неупорядоченных систем. Критические индексы, вычисленные для сильно неупорядоченных систем, демонстрируют значения значительно отличающиеся от значений, полученных для систем при малых концентрациях дефектов.

В заключении сформулированы основные результаты и выводы диссертации.

Основные результаты и выводы

  1. Осуществлено компьютерное моделирование равновесного критического поведения трехмерной слабо неупорядоченной (р = 0.80) ферромагнитной модели Изинга с линейными дефектами. Получены значения набора статических критических индексов а) (3) 7 и vмодели, индекс асимптотической поправки к скейлингу и) а также значение критической температуры Тс модели.
  2. Впервые для исследования влияния структурного беспорядка с дальней пространственной корреляцией дефектов на критическое поведение спиновых систем был применен метод коротковременной динамики. В результате численного исследовании критической релаксации для слабо неупорядоченной (р = 0.80) трехмерной модели Изинга с линейными дефектами были определены значения совокупности динамических и статических критических индексов при применении методики

19


учета поправок к скейлингу. Полученные значения критических индексов находятся в хорошем согласии с результатами теоретико-полевого описания и результатами проведенного моделирования равновесного критического поведения модели. Сделан вывод, что метод коротковременной динамики может служить надежной альтернативой традиционным методам Монте-Карло не только при численных исследованиях однородных систем, но и систем со структурным беспорядком, обеспечивая при меньших машинных затратах получение более полной информации о критическом поведении структурно неупорядоченных систем.

  1. Впервые методом коротковременной динамики осуществлено численное исследование неравновесного критического поведения трехмерной слабо неупорядоченной (р = 0.80) ферромагнитной XY-модели с линейными дефектами. Определены значения совокупности динамических и статических критических индексов при применении методики учета поправок к скейлингу. Полученные значения критических индексов находятся в хорошем согласии с результатами теоретико-полевого описания. На примере исследования XY-модели впервые получено численное подтверждение о существенности влияния дальней пространственной корреляции дефектов на критическое поведение систем с многокомпонентным параметром порядка.
  2. В отличие от результатов исследования неравновесного критического поведения однородных систем впервые выявлено наличие нескольких этапов динамического развития в эволюции слабо неупорядоченных систем после микроскопического временного масштаба: области с характеристиками однородной системы, кроссоверной области и области, характеризующейся влиянием структурного беспорядка.
  3. Впервые осуществлено компьютерное моделирование неравновесной критической динамики трехмерных сильно неупорядоченных ферромагнитных модели Изинга и XY-модели с дальнодействующей корреляцией дефектов (спиновая концентрация р = 0.60) в коротковременном режиме. Полученные значения критических индексов позволяют сделать вывод о существовании различных классов универсального критического поведения для рассматриваемых систем, отвечающих областям слабой и сильной структурной неупорядоченности.

Основные результаты диссертации опубликованы в работах

  1. Prudnikov V.V., Prudnikov P.V., Dorofeev S.V., Kolesnikov V.Y. Monte Carlo studies of critical behaviour of systems with long-range correlated disorder // Condensed Matter Physics, 2005, v.8, N1 (41), p.213-224.
  2. Prudnikov V.V., Prudnikov P.V., Zheng В., Dorofeev S.V., Kolesnikov V.Y. Short-time critical dynamics of the three-dimensional systems with long-range correlated disorder // Progress of Theoretical Physics, 2007, v.117, N6, pp.973-991.
  3. Прудников П.В., Прудников В.В., Колесников В.Ю., Медведева М.А., Желтышев П.А. Численное исследование влияния протяженных дефектов структуры на критическое поведение трехмерных систем методом коротковременной динамики // Труды Семинара по вычислительным технологиям в естественных науках. Вып. 1. Вычислительная физика / Под ред. P.P. Назирова. - М.: Изд-во КДУ, 2009. - 288 с; с.264-278.
  4. Колесников В.Ю., Прудников В.В. Монте-Карло исследования влияния дальнодействующей корреляции дефектов на критическое поведение трехмерной XY-модели // Вестник Омского

20


университета, 2004, вып. 4, с.34-36.

  1. Прудников В.В., Гергертд Е.А., Колесников В.Ю., Прудников П.В. Компьютерное моделирование критического поведения трехмерной XY-модели с линейными дефектами методом коротковременной динамики // Вестник Омского университета, 2007, вып.4, с.28-33.
  2. Prudnikov V.V., Prudnikov P.V., Zheng В., Dorofeev S.V., Kolesnikov V.Y. Short-time critical dynamics of the three-dimensional systems with long-range correlated disorder // E-print arXiv: cond-mat.diss-nn / 0709.0997, 2007, pp. 1-24.
  3. Прудников В.В., Прудников П.В., Колесников В.Ю. Численное исследование неравновесного критического поведения сильно неупорядоченной трехмерной XY-модели с линейными дефектами // Вестник Омского университета, 2009, вып.2, с.57-62.
  4. Дорофеев СВ., Колесников В.Ю., Прудников В.В., Прудников П.В. Исследование влияния дальнодействующей корреляции дефектов на критическое поведение систем методами компьютерного моделирования // Тезисы докладов III Международной конференции Фундаментальные проблемы физики>(13-18 июня 2005, г. Казань). Казань: Изд-во КГУ, 2005. с. 145.
  5. Prudnikov V.V., Prudnikov P.V., Zheng В., Dorofeev S.V., Kolesnikov V.Y. Short-time critical dynamics of systems with long-range correlated disorder // Abstracts of The 3-rd International Workshop Hangzhou 2006 on Simulational Physics, 2006, p. 14.

10.а Колесников В.Ю., Прудников В.В., Прудников П.В. Компьютерное моделирование крити

ческого поведения ферромагнитных систем с линейными дефектами методом коротковременной

динамики // Сборник трудов международной конференции Фазовые переходы, критические и

нелинейные явления в конденсированных средах>(7-10 сентября 2009, г. Махачкала). Махачка

ла: Изд-во Института физики ДагН - РАН, 2009. - 544с; с.35-38.

Список литературы

[1]а Harris А. В. // J. Phys. С. - 1974. - V. 7. - Р. 1671.

[2] Weinrib A., Halperin B.I. // Phys. Rev. В. - 1983. - V. 27. - P. 413.

[3]а Mayer I. O. // Journal of Physics. A. - 1989. - V. 22. - P. 2815.

[4]а Прудников В. В., Белим СВ., Иванов А. В., Осинцев Е.В., Федоренко А. А. // ЖЭТФ. - 1998. - Т. 114. - С. 972.

[5]а Prudnikov V. V., Prudnikov P. V. and Fedorenko A. A. // Phys. Rev. B. - 2000. - V.62. - P. 8777.

[6]а Binder K., Reger J.D. // Adv. Phys. - 1992. - V. 41. - P. 547.

[7]а Blavats'ka V., von Ferber C, Holovatch Yu. // Phys. Rev. B. - 2001. - V. 64. - P. 041102.

[8]а Korzhenevskii A.L., Luzhkov A. A., Schirmacher W. // Phys. Rev. B. - 1994. - V. 50. - P. 3661.

[9]а Hohenberg P.C., Halperin B.I. // Rev. Mod. Phys. - 1977. - V. 49. - P. 435. [10]а Janssen H.K., Schaub В., Schmittmann B. // Z. Phys. B. - 1989. - V. 73. - P. 539. [11]а Huse D. // Phys. Rev. B. - 1989. - V. 40. - P. 304. [12]а Zheng B. // Int. J. Mod. Phys. B. - 1998. - V. 12. - P. 1419.

[13]а Jaster A., Mainville J., Schulke L., Zheng B. // J. Phys. A. - 1999. - V. 32. - P. 1395. [14]а Kim J. K., de Souza A. J. and Landau D.P. // Phys. Rev. E. - 1996. - V. 54. - P. 2291. [15] Wegner F. J. // Phys. Rev. B. - 1972. - V. 5. - P. 4529.

21


Колесников Вячеслав Юрьевич

ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ВЛИЯНИЯ ДАЛЬНОДЕИСТВУЮЩЕИ КОРРЕЛЯЦИИ ДЕФЕКТОВ НА КРИТИЧЕСКОЕ ПОВЕДЕНИЕ СПИНОВЫХ

СИСТЕМ

01.04.02 - Теоретическая физика

Автореферат

диссертации на соискание ученой степени

кандидата физико-математических наук

     Авторефераты по темам  >>  Разные специальности - [часть 1]  [часть 2]