'

К вопросу о фазовом переходе парамагнетик-спиновое стекло в модели Изинга

Понравилась презентация – покажи это...





Слайд 0

К вопросу о фазовом переходе парамагнетик-спиновое стекло в модели Изинга Дальневосточный федеральный университет Школа Естественных наук П.Д.Андрющенко К.Н.Нефедев PitAndMind@gmail.ru knefedev@phys.dvgu.ru 1 Дальневосточный Федеральный Университет, Школа Естественных Наук, кафедра компьютерной безопасности, г. Владивосток, ул. Суханова 8, 690950 Дальневосточный Федеральный Университет, Школа Естественных Наук, кафедра компьютерных систем, г. Владивосток, ул. Суханова 8, 690950


Слайд 1

Актуальность 2


Слайд 2

Используемые определения 3 Конфигурация - магнитное состояние системы спинов, которое характеризуется уникальным распределением магнитных моментов по направлениям в пространстве  Вырождение - существование различных конфигураций системы, для которых некоторая физическая величина принимает одинаковые значения. Фрустрация - (от лат. frustratio – неудовлетворённость) наличие фрустраций в системе характеризуется наличием спинов, для которых основное состояние (ground state) не достижимо даже при Т=0


Слайд 3

Используемые определения 4 Спиновое стекло - магнитная фаза, в которую магнетик переходит ниже определённой критической температуры (температура замерзания Tf). В данной фазе наблюдается термодинамическое неравновесное метастабильное магнитное состояние, характеризующееся «замороженным» пространственным распределением ориентации спиновых магнитных моментов Статистическая сумма Распределение Гиббса – закон распределения вероятностей всех возможных конфигураций статистического ансамбля. Вероятность конфигурации есть функция температуры, спинового избытка, внешнего магнитного поля и обменного интеграла [1] [1]Базаров И. П., Геворкян Э. В., Николаев П. Н. Термодинамика и статистическая физика. Теория равновесных систем. — М.: МГУ, 1986.


Слайд 4

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


Слайд 5

Используемые определения 6 Фрустрация Вырождение


Слайд 6

Модель Изинга +1 -1 Фрустрированная модель спинового стекла на простой квадратной решетке со спинами Изинга Схема Реальная решетка 7 Атомы основной решетки Атомы примеси


Слайд 7

Модель Изинга Система описывается с помощью функционала энергии, называющегося гамильтонианом, который ассоциирует каждую уникальную конфигурацию системы с конкретным значением энергии Jij - обменный интеграл h - внешнее магнитное поле 8 [2]Гинзбург С. Л., Необратимые явления в спиновых стеклах, М., 1989 [2]


Слайд 8

Модель Изинга При этом вероятность любой из возможных конфигураций задается распределением Гиббса. В нашей модели каждый спин взаимодействует только с ближайшими z = 4 соседями посредством прямого ферромагнитного или антиферромагнитного обменного взаимодействия, распределенного случайным образом в узлах решетки, с условием, что 9


Слайд 9

Модель Изинга 10


Слайд 10

Модель Изинга +1 -1 Фрустрированная модель спинового стекла на простой квадратной решетке со спинами Изинга Схема Реальная решетка 11 Атомы основной решетки Атомы примеси


Слайд 11

Спиновое стекло 12 СПИНОВОЕ СТЕКЛО - магнетик, в котором ниже определённой температуры (температура замерзания Tf)возникает термодинамическое неравновесное метастабильное магнитное состояние, которое характеризуется «замороженным» (отсутствуют термодинамические флуктуации) пространственным распределением ориентации спиновых магнитных моментов AuFe CuMn AgMn  CuCo AuMn


Слайд 12

Экспериментальные данные Низкополевая магнитная восприимчивость ?(T) сплавов AuFe с концентрацией железа 1, 2, 5 и 8 ат. % [3] 13 [3]Г. А. Петраковский, Спиновые стекла, Соросовский образовательный журнал, т. 7, №9, 2001 Зависимость удельной теплоемкости от температуры в полях H = 0 и 5Т. На вставке показана низкотемпературная область при H = 0 T [4] [4]Г. А. Петраковский, Состояние спинового стекла в , Физика твердого тела, 2011, том53, вып. 9


Слайд 13

Алгоритм 14 Разработан программный инструментарий для моделирования динамики физических величин, характеризующих систему с заданным (в т.ч. знакопеременным) короткодействующим взаимодействием, а также исследования фазовых переходов в системах с разными значениями обменного интеграла (ферромагнетик, антиферромагнетик, спиновое стекло). Алгоритм реализует схему Монте – Карло.


Слайд 14

Алгоритм 15 1. Создание решетки 1.1 Создание массива с магнитными моментами(спинами) 1.2 Создание массива со связями 1.2 Создание массива энергий Линейный размер = n


Слайд 15

Алгоритм 16 2. Создание фрустраций Задача: Сделать решетку фрустрированной (для спинового стекла) Сложность заключается в том, что необходимо соблюсти граничные условия [5]К.Биндер, Д.Хеерман. Моделирование методом Монте-Карло в статистической физике. Перевод с англ., "Наука", Москва, 1995. Мелкие системы (порядка 5000 элементов) конструируюстя без ошибок В более крупных системах возникают ошибки (порядка 0.5%)


Слайд 16

17 Алгоритм 3. Монте-Карло При моделировании возможно задать количество Монте-Карло проходов Каждый Монте-Карло проход состоит: 1.Случайным образом генерируем число от 1 до .Выбранное число соответствует порядковому номеру атома в решетке. 2. Считаем его энергию Е1, записываем. Затем переворачиваем его магнитный момент и считаем энергию Е2. 3. Сравниваем энергии Е1 и Е2. Если Е1>Е2, то вероятность того, что спин перевернется из исходного состояния P = 1, если же Е1<Е2, то вероятность высчитывается по формуле: 4. Генерируем число S случайным образом от 0 до 1. Считаем, что спин перевернулся, если S > P. Иначе, спин остается в исходном состоянии. 5. Повторяем пункты 1-4 раз.


Слайд 17

18 Алгоритм 4. Расчет параметра порядка Для разделения парамагнитной и спинстекольной фаз предлагается использовать параметр «порядка», который представляет собой отношение количества элементов, находящихся в максимальном кластере с заданной энергией -4 (-4 и -2) к общему числу элементов, т.е. порог протекания или порог перколяции. Легко видеть, что данная задача сходна с задачей обхода графа. Существуют несколько возможных алгоритмов обхода графа.


Слайд 18

19 Алгоритм 4.1 Суть метода обхода в глубину: 2. Рассматривается вершина V2, смежная с V1. Она выбирается. 1. Поиск начинается с некоторой фиксированной вершины V1 3. Процесс повторяется с выбранной вершиной. 4. Если на очередном шаге мы работаем с вершиной Vn, и нет вершин, смежных с ней и не рассмотренных ранее, то возвращаемся из вершины Vn к вершине, которая была до нее. Если эта вершина V1, то процесс просмотра закончен


Слайд 19

20 Алгоритм 4.2 Суть метода обхода в ширину: Такой порядок обхода обеспечивается благодаря тому, что создается очередь. Рассматривается первая вершина, в конец очереди записываются смежные с ней вершины, затем рассматривается следующая в очереди, добавляются ее смежные вершины в конец очереди и так до тех пор, пока есть вершины в очереди. Именно такой подход мы и реализовали в программе. При поиске в ширину сначала рассматриваются все вершины, смежные со стартовой, то есть находящиеся от нее на расстоянии 1, затем вершины, находящиеся от старта на расстоянии 2, и т.д.


Слайд 20

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


Слайд 21

22 Результаты Ферромагнетик 2D ПК 4 соcеда Подпись указать -4, М, (-4)+(-2), число МК проходов, размер системы, 1000х1000 характеристики, число всех конфигураций,


Слайд 22

23 Результаты 2. Спиновое стекло (подпись (-4)+(-2))


Слайд 23

Зависимость пп от числа проходов МК 24


Слайд 24

Cверхмасшабирование и дальнейшее развитие исследований Будет проведено распараллеливание алгоритма с целью увеличения числа частиц в системе и генерации большего числа конфигураций за МК проход. Планируется рассчитать поведение исследуемых систем во внешнем магнитном поле, установить зависимость пп от Т для заданного внешнего поля. Планируется вычислить температурную зависимость теплоемкости и магнитной восприимчивости, в т.ч. в ZFC и FC режимах. 25


Слайд 25

Выводы Разработан алгоритм, написана программа ЭВМ на языке С#. Проведены численные эксперименты моделирующие поведение ферромагнетика, и спинового стекла на решетке Изинга. Предложена схема вычисления параметра порядка. Установлена его температурная зависимость, а также критические температуры фазовых переходов в парамагнетик-спиновое стекло, суперпарамагнетк-парамагнетик. 26


Слайд 26

Заключение Закон температурного спада предлагаемого параметра «порядка» (-4) для ферромагнитной фазы совпадает с законом температурного поведения намагниченности с критическим индексом 1/8 (Онсагер []) Существование неаналитической функции в области перехода спиновое стекло-парамагнетик позволяет утвердительно ответить на вопрос о существовании фазового перехода PM-SG (PM-SPM) 27


Слайд 27

Спасибо за внимание 28


×

HTML:





Ссылка: