Из Википедии, бесплатной энциклопедии
  (Перенаправлено из модели XY )
Перейти к навигации Перейти к поиску

Классическая модель XY (иногда также называемый классический ротор ( вращатель ) модель или O (2) модель ) является решеткой модели из статистической механики . В общем, XY-модель можно рассматривать как специализацию n -векторной модели Стэнли [1] для n = 2 .

Определение [ править ]

Для D -мерной решетки Λ на каждый узел j ∈ Λ приходится двумерный вектор единичной длины s j = (cos θ j , sin θ j )

Конфигурации спины , ев = ( ев J ) J ∈ Л является заданием угла - л < θ J & le ; л для каждого J ∈ Л .

Учитывая трансляционно-инвариантное взаимодействие J ij = J ( i - j ) и точечно-зависимое внешнее поле , конфигурационная энергия равна

Случай, когда J ij = 0, за исключением ближайшего соседа ij , называется случаем ближайшего соседа .

Вероятность конфигурации задается распределением Больцмана с обратной температурой р ≥ 0 :

где Z - нормализация или статистическая сумма . [2] Обозначение указывает на ожидание случайной величины A ( s ) в пределе бесконечного объема после наложения периодических граничных условий .

Строгие результаты [ править ]

Следовательно, критическое β модели XY не может быть меньше, чем удвоенная критическая температура модели Изинга.

Одно измерение [ править ]

Как и в любой n -векторной модели «ближайшего соседа» со свободными (непериодическими) граничными условиями, если внешнее поле равно нулю, существует простое точное решение. В случае свободных граничных условий гамильтониан имеет вид

поэтому статистическая сумма факторизуется при изменении координат

Это дает

где - модифицированная функция Бесселя первого рода. Статистическая сумма может использоваться для нахождения нескольких важных термодинамических величин. Например, в термодинамическом пределе ( ) свободная энергия на спин равна

Используя свойства модифицированных функций Бесселя, удельную теплоемкость (на спин) можно выразить как [5]

где , - короткодействующая корреляционная функция,

Точная удельная теплоемкость на спин в одномерной XY-модели

Даже в термодинамическом пределе нет расхождения в теплоемкости. Действительно, как и одномерная модель Изинга, одномерная XY-модель не имеет фазовых переходов при конечной температуре.

То же самое вычисление для периодического граничного условия (и все еще h = 0 ) требует формализма матрицы переноса , хотя результат тот же. [6] .

(Щелкните "показать" справа, чтобы просмотреть подробные сведения о формализме матрицы передачи.)

Статистическая сумма может быть оценена как

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

Обратите внимание на расширение

которая представляет собой диагональное матричное представление в основе собственных функций плоских волн . Собственные значения матрицы просто являются модифицированными функциями Бесселя, вычисленными в , а именно . Для любого конкретного значения эти модифицированные функции Бесселя удовлетворяют и . Следовательно, в термодинамическом пределе собственное значение будет преобладать над следом, и так далее .

Такой подход матрицы переноса также требуется при использовании свободных граничных условий, но с приложенным полем . Если приложенное поле достаточно мало, чтобы его можно было рассматривать как возмущение системы в нулевом поле, то можно оценить магнитную восприимчивость . Это делается путем использования собственных состояний, вычисленных с помощью подхода матрицы переноса, и вычисления сдвига энергии с помощью теории возмущений второго порядка , а затем сравнения с разложением по свободной энергии . Находят [7]

где - постоянная Кюри (значение, обычно связанное с восприимчивостью магнитных материалов). Это выражение справедливо и для одномерной модели Изинга с заменой .

Два измерения [ править ]

Среднеквадратичная магнетизация для решетки 25x25 (Ota: [8] 30x30), предполагающая увеличение магнитного момента, которого нет в термодинамическом пределе.

Двумерная XY-модель с взаимодействиями ближайших соседей является примером двумерной системы с непрерывной симметрией, которая не имеет дальнего порядка, как того требует теорема Мермина – Вагнера . Точно так же нет обычного фазового перехода, который был бы связан с нарушением симметрии . Однако, как будет обсуждаться позже, система действительно показывает признаки перехода из неупорядоченного высокотемпературного состояния в квазиупорядоченное состояние ниже некоторой критической температуры, называемого переходом Костерлица-Таулеса.. В случае дискретной решетки спинов, двумерная XY-модель может быть оценена с использованием подхода матрицы переноса, сводя модель к задаче на собственные значения и используя наибольшее собственное значение из матрицы переноса. Хотя точное решение трудно найти, можно использовать определенные приближения, чтобы получить оценки критической температуры, которая возникает при низких температурах. Например, Маттис (1984) использовал приближение к этой модели, чтобы оценить критическую температуру системы как

Двухмерная XY-модель также была детально изучена с использованием моделирования Монте-Карло , например, с помощью алгоритма Метрополиса . Их можно использовать для вычисления термодинамических величин, таких как энергия системы, удельная теплоемкость, намагниченность и т. Д., В диапазоне температур и временных масштабов. В моделировании Монте-Карло каждый спин связан с непрерывно изменяющимся углом (часто его можно дискретизировать на конечное число углов, как в связанной модели Поттса , для простоты вычислений. Однако это не является обязательным требованием). На каждом временном шаге алгоритм Метрополиса выбирает одно вращение случайным образом и поворачивает его угол на некоторый случайный шаг . Это изменение угла вызывает изменение энергиисистемы, которая может быть положительной или отрицательной. Если отрицательный, алгоритм принимает изменение угла; если положительный, конфигурация принимается с вероятностью , фактор Больцмана для изменения энергии. Для проверки критической температуры системы различными методами использовался метод Монте-Карло, который оценивается в [9] . Метод Монте-Карло также может вычислять средние значения, которые используются для вычисления термодинамических величин, таких как намагниченность, спин-спиновая корреляция, длины корреляции и удельная теплоемкость. Это важные способы охарактеризовать поведение системы вблизи критической температуры. Например, намагниченность и квадрат намагниченности можно вычислить как

Удельная теплоемкость двумерного XY с использованием моделирования Монте-Карло, показывающая особенность в точке над переходом KT

где количество вращений. Средняя намагниченность характеризует величину чистого магнитного момента системы; во многих магнитных системах она равна нулю выше критической температуры и самопроизвольно становится ненулевой при низких температурах. Точно так же среднеквадратичная намагниченность характеризует среднее значение квадрата чистых компонентов спинов по решетке. Любой из них обычно используется для характеристики параметра порядка системы. Строгий анализ модели XY показывает, что намагниченность в термодинамическом пределе равна нулю, и что квадрат намагниченности приблизительно следует [10] , которая обращается в нуль в термодинамическом пределе. Действительно, при высоких температурах эта величина приближается к нулю, поскольку компоненты спинов будут иметь тенденцию быть случайными и, следовательно, их сумма будет равна нулю. Однако при низких температурах для конечной системы среднеквадратичная намагниченность увеличивается, что свидетельствует о наличии областей спинового пространства, которые выровнены, чтобы вносить ненулевой вклад. Показанная намагниченность (для решетки 25x25) является одним из примеров этого, который, по-видимому, предполагает фазовый переход, в то время как такой переход не существует в термодинамическом пределе.

Кроме того, используя статистическую механику, можно связать средние термодинамические величины с такими величинами, как удельная теплоемкость, вычисляя

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

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

Цветовая карта (дискретной) двумерной XY-модели в решетке 250x250 в . Каждое вращение представлено цветом, соответствующим углу между

Непрерывная версия XY-модели часто используется для моделирования систем, которые обладают параметрами порядка с одинаковыми видами симметрии, например, сверхтекучий гелий , гексатические жидкие кристаллы . Это отличает их от других фазовых переходов, которые всегда сопровождаются нарушением симметрии. Топологические дефекты в XY-модели приводят к развязывающему вихрь переходу из низкотемпературной фазы в высокотемпературную неупорядоченную фазу . Действительно, тот факт, что при высоких температурах корреляции затухают экспоненциально быстро, а при низких температурах затухают по степенному закону, даже если в обоих режимах M ( β ) = 0 , называется переходом Костерлица – Таулеса.. Костерлиц и Таулесс представили простой аргумент, почему это могло быть так: это рассматривает основное состояние, состоящее из всех спинов с одинаковой ориентацией, с добавлением затем одного вихря. Их присутствие дает примерно энтропию , где - эффективный масштаб длины (например, размер решетки для дискретной решетки). Между тем, энергия системы увеличивается из-за вихря на определенную величину . Если сложить все вместе, свободная энергия системы изменится из-за спонтанного образования вихря на величину

В термодинамическом пределе система не способствует образованию вихрей при низких температурах, но поддерживает их при высоких температурах, превышающих критическую температуру . Это указывает на то, что при низких температурах любые возникающие вихри будут стремиться аннигилировать с антивихрями, чтобы снизить энергию системы. Действительно, качественно это будет иметь место, если наблюдать «снимки» спиновой системы при низких температурах, когда вихри и антивихри постепенно объединяются и аннигилируют. Таким образом, низкотемпературное состояние будет состоять из связанных пар вихрь-антивихрь. Между тем при высоких температурах будет набор несвязанных вихрей и антивихрей, которые могут свободно перемещаться по плоскости.

Для визуализации модели Изинга можно использовать стрелку, указывающую вверх или вниз, или представить ее в виде точки, окрашенной в черный / белый цвет, чтобы указать ее состояние. Чтобы визуализировать систему вращения XY, вращения можно представить в виде стрелки, указывающей в каком-либо направлении, или в виде точки с некоторым цветом. Здесь необходимо представить спин со спектром цветов, обусловленным каждой из возможных непрерывных переменных. Это можно сделать, используя, например, непрерывный и периодический красный-зеленый-синий спектр. Как показано на рисунке, голубой соответствует нулевому углу (указывает вправо), а красный соответствует углу 180 градусов (указывает влево). Затем можно изучить снимки спиновых конфигураций при разных температурах, чтобы выяснить, что происходит выше и ниже критической температуры XY-модели. При высоких температурахспины не будут иметь предпочтительной ориентации, и будет происходить непредсказуемое изменение углов между соседними спинами, так как не будет предпочтительной энергетически выгодной конфигурации. В этом случае цветовая карта будет выглядеть сильно пиксельной. Между тем при низких температурах одна возможная конфигурация основного состояния имеет все спины, направленные в одну и ту же ориентацию (один и тот же угол); они будут соответствовать областям (доменам) цветовой карты, где все спины имеют примерно одинаковый цвет.они будут соответствовать областям (доменам) цветовой карты, где все спины имеют примерно одинаковый цвет.они будут соответствовать областям (доменам) цветовой карты, где все спины имеют примерно одинаковый цвет.

Различные формы вихрей и антивихрей, показанные в моделировании Монте-Карло на

Чтобы идентифицировать вихри (или антивихри), присутствующие в результате перехода Костерлица – Таулеса, можно определить изменение угла со знаком, пройдя круг из точек решетки против часовой стрелки. Если полное изменение угла равно нулю, это соответствует отсутствию вихря; тогда как полное изменение угласоответствует вихрю (или антивихрю). Эти вихри являются топологически нетривиальными объектами, которые входят в пары вихрь-антивихрь, которые могут разделяться или парно-аннигилировать. На палитре эти дефекты можно идентифицировать в областях с большим цветовым градиентом, где все цвета спектра встречаются вокруг точки. Качественно эти дефекты могут выглядеть как источники потока, направленные внутрь или наружу, или как водовороты вращений, которые все вместе вращаются по часовой стрелке или против часовой стрелки, или как объекты гиперболического вида с некоторыми вращениями, направленными в сторону дефекта, и некоторыми вращениями, направленными от него. Поскольку конфигурация изучается на больших временных масштабах и при низких температурах, наблюдается, что многие из этих пар вихрь-антивихрь сближаются и в конечном итоге пары аннигилируют.Только при высоких температурах эти вихри и антивихри высвобождаются и развязываются друг с другом.

В непрерывной XY-модели высокотемпературная спонтанная намагниченность исчезает:

Кроме того, расширение кластера показывает, что спиновые корреляции кластеризуются экспоненциально быстро: например

При низких температурах, т.е. β ≫ 1 , спонтанная намагниченность остается нулевой (см. Теорему Мермина – Вагнера ),

но затухание корреляций происходит только по степенному закону: Фрёлих и Спенсер [13] нашли нижнюю границу

в то время как Макбрайан и Спенсер нашли верхнюю границу для любых

Три и выше измерения [ править ]

Независимо от диапазона взаимодействия при достаточно низкой температуре намагниченность положительна.

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

Фазовый переход [ править ]

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

В трех и более измерениях XY-модель имеет фазовый переход ферромагнетик-парамагнетик. При низких температурах спонтанная намагниченность отлична от нуля: это ферромагнитная фаза. С повышением температуры спонтанная намагниченность постепенно уменьшается и исчезает при критической температуре. Он остается нулевым при всех более высоких температурах: это парамагнитная фаза.

В четырех и более измерениях фазовый переход имеет критические показатели среднего поля (с логарифмическими поправками в четырех измерениях).

Трехмерный случай: критические показатели [ править ]

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

Наиболее важными критическими показателями трехмерной модели XY являются . Все они выражаются всего двумя числами: масштабными размерностями и комплексным полем параметра порядка и ведущим синглетным оператором (как в описании Гинзбурга – Ландау ). Еще одно важное поле (такое же, как ), размерность которого определяет показатель масштабирования поправки . Согласно расчетам конформного бутстрапа, [14] эти три измерения задаются следующим образом:

Это дает следующие значения критических показателей:

Методы Монте-Карло дают совместимые определения: [15] .

См. Также [ править ]

  • Классическая модель Гейзенберга
  • Бозон Голдстоуна
  • Модель Изинга
  • Модель Поттса
  • n- векторная модель
  • Переход Костерлица – Таулеса
  • Топологический дефект
  • Сверхтекучая пленка
  • Сигма модель

Заметки [ править ]

  1. ^ Стэнли, HE (1968). «Зависимость критических свойств от размерности спинов». Phys. Rev. Lett . 20 (12): 589–592. Bibcode : 1968PhRvL..20..589S . DOI : 10.1103 / PhysRevLett.20.589 .
  2. ^ Чайкин, PM; Лубенский, ТК (2000). Основы физики конденсированного состояния . Издательство Кембриджского университета. ISBN 978-0521794503.
  3. ^ Жинибра, J. (1970). «Общая формулировка неравенств Гриффитса». Comm. Математика. Phys . 16 (4): 310–328. Bibcode : 1970CMaPh..16..310G . DOI : 10.1007 / BF01646537 . S2CID 120649586 . 
  4. ^ Aizenman, M .; Саймон Б. (1980). «Сравнение ротора самолета и моделей Изинга». Phys. Lett. . 76 (3–4): 281–282. Bibcode : 1980PhLA ... 76..281A . DOI : 10.1016 / 0375-9601 (80) 90493-4 .
  5. ^ Бадалян, Д. (1996). «К термодинамике классических спинов с изотропным гейзенберговским взаимодействием в одномерных квазипериодических структурах». Physica B . 226 (4): 385–390. Bibcode : 1996PhyB..226..385B . DOI : 10.1016 / 0921-4526 (96) 00283-9 .
  6. Перейти ↑ Mattis, DC (1984). «Передаточная матрица в модели плосковращателя». Phys. Lett . 104 А (6–7): 357–360. Bibcode : 1984PhLA..104..357M . DOI : 10.1016 / 0375-9601 (84) 90816-8 .
  7. Перейти ↑ Mattis, DC (1985). Теория магнетизма II . Серия Спрингера по физике твердого тела. ISBN 978-3-642-82405-0.
  8. ^ Ota, S .; Ота, SB; Фанле, М. (1992). «Микроканоническое моделирование методом Монте-Карло для двумерной XY-модели». J. Phys .: Condens. Материя . 4 (24): 5411. Bibcode : 1992JPCM .... 4.5411O . DOI : 10.1088 / 0953-8984 / 4/24/011 .
  9. ^ Hsieh, Y.-D .; Kao, Y.-J .; Сандвик, AW (2013). "Метод конечных размеров для перехода Березинского-Костерлица-Таулеса". J. Stat. Механизм .: Теория Эксп . 2013 (9): P09001. arXiv : 1302.2900 . Bibcode : 2013JSMTE..09..001H . DOI : 10.1088 / 1742-5468 / 2013/09 / P09001 . S2CID 118609225 . 
  10. ^ Tobochnik, J .; Честер, GV (1979). "Монте-Карло исследование модели плоского спина". Phys. Rev. B . 20 (9): 3761–3769. Bibcode : 1979PhRvB..20.3761T . DOI : 10.1103 / PhysRevB.20.3761 .
  11. Перейти ↑ Binder, K. (2013). Применение метода Монте-Карло в статистической физике . Springer Science & Business Media. ISBN 978-3-642-51703-7.
  12. ^ Ван Химберген, JE; Чакраварти, Судип (1981). «Модуль спиральности и удельная теплоемкость классической модели XY в двух измерениях». Phys. Rev. B . 23 (1): 359–361. Bibcode : 1981PhRvB..23..359V . DOI : 10.1103 / PhysRevB.23.359 .
  13. ^ Fröhlich, J .; Спенсер, Т. (1981). «Переход Костерлица – Таулеса в двумерных абелевых спиновых системах и кулоновском газе» . Comm. Математика. Phys . 81 (4): 527–602. Bibcode : 1981CMaPh..81..527F . DOI : 10.1007 / bf01208273 . S2CID 73555642 . 
  14. ^ Честер, Шай М .; Лэндри, Уолтер; Лю, Цзюнюй; Польша, Давид; Симмонс-Даффин, Дэвид; Су, Нин; Вичи, Алессандро (2020). «Вырезание пространства OPE и точных критических показателей модели O (2)» . Журнал физики высоких энергий . 2020 (6): 142. arXiv : 1912.03324 . Bibcode : 2020JHEP ... 06..142C . DOI : 10.1007 / JHEP06 (2020) 142 . ISSN 1029-8479 . S2CID 208910721 .  
  15. ^ Hasenbusch, Martin (2019-12-26). «Исследование Монте-Карло усовершенствованной модели часов в трех измерениях» . Physical Review B . 100 (22): 224517. arXiv : 1910.05916 . Bibcode : 2019PhRvB.100v4517H . DOI : 10.1103 / PhysRevB.100.224517 . ISSN 2469-9950 . S2CID 204509042 .  

Ссылки [ править ]

  • Евгений Демидов, Вихри в XY-модели (2004)

Дальнейшее чтение [ править ]

  • Его Превосходительство Стэнли, Введение в фазовые переходы и критические явления , (Oxford University Press, Oxford and New York, 1971);
  • Х. Кляйнерт , Калибровочные поля в конденсированных средах , Vol. I, «СУПЕРПОТОК И ВИХРЕВЫЕ ЛИНИИ», стр. 1–742, Vol. II, «НАПРЯЖЕНИЯ И ДЕФЕКТЫ», стр. 743–1456, World Scientific (Сингапур, 1989) ; ISBN в мягкой обложке 9971-5-0210-0 (также доступно в Интернете: том I и том II ) 

Внешние ссылки [ править ]

  • XY-модель в реальном времени WebGL-симуляция
  • Интерактивное моделирование методом Монте-Карло моделей Изинга, XY и Гейзенберга с трехмерной графикой (требуется браузер, совместимый с WebGL)