Классическая модель XY - Classical XY model
В классическая модель XY (иногда также называют классический ротор (ротатор) модель или же O (2) модель) это решетчатая модель из статистическая механика. В целом модель XY можно рассматривать как специализацию Стэнли. п-векторная модель [1] за п = 2.
Определение
Учитывая D-размерный решетка Λ, на каждый узел решетки j ∈ Λ есть двухмерный, вектор единичной длины sj = (cos θjгрех θj)
В конфигурация вращения, s = (sj)j ∈ Λ это задание угла −π < θj ≤ π для каждого j ∈ Λ.
Учитывая переводно-инвариантный взаимодействие Jij = J(я − j) и точечно-зависимое внешнее поле , то конфигурационная энергия является
Случай, когда Jij = 0 кроме ij ближайший сосед называется ближайший сосед дело.
В вероятность конфигурации дается Распределение Больцмана с обратной температурой β ≥ 0:
куда Z это нормализация, или же функция распределения.[2] Обозначение указывает на ожидание случайной величины А(s) в пределе бесконечного объема, после периодические граничные условия были наложены.
Строгие результаты
- Существование термодинамический предел для свободная энергия и спиновые корреляции были доказаны Ginibre, распространяя на этот случай Неравенство гриффитса.[3]
- С использованием Неравенство гриффитса в постановке Жинибра, Айзенмана и Саймона[4] доказал, что двухточечная спиновая корреляция ферромагнетики XY модель в размере D, связь J > 0 и обратная температура β является преобладают по (т.е. имеет верхняя граница дается) двухточечной корреляцией ферромагнитный Модель Изинга в измерении D, связь J > 0 и обратная температура β/2
- Следовательно, критический β модели XY не может быть меньше двойной критической температуры модели Изинга
Одно измерение
Как и любой "ближайший сосед" п-векторная модель со свободными (непериодическими) граничными условиями, если внешнее поле равно нулю, существует простое точное решение. В случае свободных граничных условий гамильтониан имеет вид
поэтому статистическая сумма факторизуется при изменении координат
Это дает
куда это модифицированная функция Бесселя первого вида. Статистическая сумма может использоваться для нахождения нескольких важных термодинамических величин. Например, в термодинамическом пределе (), свободная энергия за спин
Используя свойства модифицированных функций Бесселя, удельную теплоемкость (на спин) можно выразить как[5]
куда , и - короткодействующая корреляционная функция,
Даже в термодинамическом пределе нет расхождения в теплоемкости. Действительно, как и одномерная модель Изинга, одномерная XY-модель не имеет фазовых переходов при конечной температуре.
То же самое вычисление для периодического граничного условия (и все еще час = 0) требует формализм матрицы переноса, хотя результат тот же.[6].
Статистическая сумма может быть оценена как
который можно рассматривать как след матрицы, а именно произведение матриц (в данном случае скаляров). След матрицы - это просто сумма ее собственных значений, и в термодинамическом пределе выживет только наибольшее собственное значение, поэтому статистическая сумма может быть записана как повторяющееся произведение этого максимального собственного значения. Это требует решения проблемы собственных значений
Обратите внимание на расширение
Такой подход матрицы переноса также требуется при использовании свободных граничных условий, но с приложенным полем . Если прикладное поле достаточно мал, чтобы его можно было рассматривать как возмущение системы в нулевом поле, то магнитная восприимчивость можно оценить. Это делается с помощью собственных состояний, вычисленных с помощью подхода матрицы переноса, и вычисления сдвига энергии со вторым порядком теория возмущений, то сравнивая с разложением по свободной энергии . Один находит [7]
куда это Постоянная Кюри (значение, обычно связанное с восприимчивостью магнитных материалов). Это выражение справедливо и для одномерной модели Изинга с заменой .
Два измерения
Двумерная XY-модель с взаимодействиями ближайших соседей является примером двумерной системы с непрерывной симметрией, которая не имеет дальнего порядка, как того требует Теорема Мермина – Вагнера. Точно так же нет обычного фазового перехода, который был бы связан с нарушение симметрии. Однако, как будет обсуждаться позже, система действительно показывает признаки перехода из неупорядоченного высокотемпературного состояния в квазиупорядоченное состояние ниже некоторой критической температуры, называемое Переход Костерлица-Таулеса. В случае дискретной решетки спинов двумерная XY-модель может быть оценена с использованием подхода матрицы передачи, сводя модель к задаче на собственные значения и используя наибольшее собственное значение из матрицы передачи. Хотя точное решение трудно найти, можно использовать определенные приближения для получения оценок критической температуры. что происходит при низких температурах. Например, Маттис (1984) использовал приближение к этой модели, чтобы оценить критическую температуру системы как
Двухмерная XY-модель также была детально изучена с использованием Монте-Карло моделирования, например, с Алгоритм мегаполиса. Их можно использовать для вычисления термодинамических величин, таких как энергия системы, удельная теплоемкость, намагниченность и т. Д., В диапазоне температур и временных масштабов. В моделировании Монте-Карло каждый спин связан с непрерывно меняющимся углом (часто его можно разделить на конечное число углов, как в соответствующем Модель Поттса, для простоты вычислений. Однако это не является обязательным требованием.) На каждом временном шаге алгоритм Метрополиса выбирает одно вращение случайным образом и поворачивает его угол на некоторое случайное приращение. . Это изменение угла вызывает изменение энергии системы, которая может быть положительной или отрицательной. Если отрицательный, алгоритм принимает изменение угла; если положительный, конфигурация принимается с вероятностью , то Фактор Больцмана для изменения энергии. Метод Монте-Карло использовался для проверки различными методами критической температуры системы и оценивается как[9] . Метод Монте-Карло также может вычислять средние значения, которые используются для вычисления термодинамических величин, таких как намагниченность, спин-спиновая корреляция, корреляционные длины и удельная теплоемкость. Это важные способы характеристики поведения системы вблизи критической температуры. Намагниченность и квадрат намагниченности, например, можно вычислить как
куда количество вращений. Средняя намагниченность характеризует величину чистого магнитного момента системы; во многих магнитных системах он равен нулю выше критической температуры и самопроизвольно становится ненулевым при низких температурах. Точно так же среднеквадратичная намагниченность характеризует среднее значение квадрата чистых компонентов спинов по решетке. Любой из них обычно используется для характеристики параметра порядка системы. Строгий анализ модели XY показывает, что намагниченность в термодинамическом пределе равна нулю, а квадрат намагниченности приблизительно соответствует[10] , которая обращается в нуль в термодинамическом пределе. Действительно, при высоких температурах эта величина приближается к нулю, поскольку компоненты спинов будут иметь тенденцию быть случайными и, следовательно, их сумма будет равна нулю. Однако при низких температурах для конечной системы среднеквадратичная намагниченность увеличивается, что позволяет предположить, что существуют области спинового пространства, которые выровнены, чтобы вносить ненулевой вклад. Показанная намагниченность (для решетки 25x25) является одним из примеров этого, который, по-видимому, предполагает фазовый переход, в то время как такой переход не существует в термодинамическом пределе.
Кроме того, используя статистическую механику, можно связать средние термодинамические величины с такими величинами, как удельная теплоемкость, вычисляя
Удельная теплоемкость указана при низких температурах вблизи критической температуры. . В удельной теплоемкости нет характеристики, согласующейся с критической характеристикой (например, дивергенцией) при этой прогнозируемой температуре. Действительно, оценка критической температуры исходит из других методов, например, из модуль спиральности, или температурная зависимость дивергенции восприимчивости.[11] Однако в теплоемкости есть особенность в виде пика вблизи . Было показано, что это положение пика и высота пика зависят от размера системы;[12] однако этот элемент остается конечным для всех размеров решетки и, похоже, сходится к конечному значению (хотя не исключено, что элемент является куспидом, это маловероятно).
Природу критических переходов и образования вихрей можно выяснить, рассмотрев непрерывную версию XY-модели. Здесь дискретные спины заменяются полем представляющий угол вращения в любой точке пространства. В этом случае угол вращения должны плавно меняться при смене положения. Разложив исходный косинус в ряд Тейлора, гамильтониан можно выразить в приближении континуума как
Непрерывная версия XY модели часто используется для моделирования систем, которые обладают параметрами порядка с одинаковыми видами симметрии, например сверхтекучий гелий, гексатические жидкие кристаллы. Это отличает их от других фазовых переходов, которые всегда сопровождаются нарушением симметрии. Топологические дефекты в XY-модели приводят к вихревой переход от низкотемпературной фазы к высокотемпературной неупорядоченная фаза. Действительно, тот факт, что при высоких температурах корреляции затухают экспоненциально быстро, а при низких температурах затухают по степенному закону, хотя в обоих режимах M(β) = 0, называется Переход Костерлица – Таулеса. Костерлиц и Таулесс представили простой аргумент, почему это могло быть так: это рассматривает основное состояние, состоящее из всех спинов с одинаковой ориентацией, с добавлением одного вихря. Их присутствие дает энтропию примерно , куда - эффективный масштаб длины (например, размер решетки для дискретной решетки). Между тем, энергия системы увеличивается из-за вихря на величину . Если сложить их вместе, свободная энергия системы изменится из-за спонтанного образования вихря на величину
В термодинамическом пределе система не способствует образованию вихрей при низких температурах, но способствует им при высоких температурах, превышающих критическую температуру. . Это указывает на то, что при низких температурах любые возникающие вихри будут стремиться аннигилировать с антивихрями, чтобы снизить энергию системы. Действительно, качественно это будет так, если наблюдать «снимки» спиновой системы при низких температурах, где вихри и антивихри постепенно объединяются и аннигилируют. Таким образом, низкотемпературное состояние будет состоять из связанных пар вихрь-антивихрь. Между тем, при высоких температурах будет набор несвязанных вихрей и антивихрей, которые могут свободно перемещаться по плоскости.
Чтобы визуализировать модель Изинга, можно использовать стрелку, указывающую вверх или вниз, или представить ее в виде точки, окрашенной в черный / белый цвет, чтобы указать ее состояние. Чтобы визуализировать систему вращения XY, вращения можно представить в виде стрелки, указывающей в каком-либо направлении, или в виде точки с некоторым цветом. Здесь необходимо представить спин со спектром цветов, обусловленным каждой из возможных непрерывных переменных. Это можно сделать, используя, например, непрерывный и периодический красный-зеленый-синий спектр. Как показано на рисунке, голубой соответствует нулевому углу (указывает вправо), а красный соответствует углу 180 градусов (указывает влево). Затем можно изучить снимки спиновых конфигураций при разных температурах, чтобы выяснить, что происходит выше и ниже критической температуры XY-модели. При высоких температурах спины не будут иметь предпочтительной ориентации, и будет происходить непредсказуемое изменение углов между соседними спинами, так как не будет предпочтительной энергетически выгодной конфигурации. В этом случае цветовая карта будет выглядеть сильно пиксельной. Между тем при низких температурах одна возможная конфигурация основного состояния имеет все спины, направленные в одну и ту же ориентацию (один и тот же угол); они будут соответствовать областям (доменам) цветовой карты, где все спины имеют примерно одинаковый цвет.
Чтобы идентифицировать вихри (или антивихри), присутствующие в результате перехода Костерлица – Таулеса, можно определить изменение угла со знаком, пройдя круг точек решетки против часовой стрелки. Если полное изменение угла равно нулю, это соответствует отсутствию вихря; тогда как полное изменение угла соответствует вихрю (или антивихрю). Эти вихри являются топологически нетривиальными объектами, которые входят в пары вихрь-антивихрь, которые могут разделяться или парно-аннигилировать. На палитре эти дефекты можно идентифицировать в областях с большим цветовым градиентом, где все цвета спектра встречаются вокруг точки. Качественно эти дефекты могут выглядеть как источники потока, направленные внутрь или наружу, или как водовороты вращений, которые вместе вращаются по часовой стрелке или против часовой стрелки, или как объекты гиперболического вида с некоторыми вращениями, направленными в сторону дефекта, и некоторыми вращениями, направленными от дефекта. Поскольку конфигурация изучается на больших временных масштабах и при низких температурах, наблюдается, что многие из этих пар вихрь-антивихрь сближаются и в конечном итоге пары аннигилируют. Только при высоких температурах эти вихри и антивихри высвобождаются и развязываются друг с другом.
В непрерывной XY-модели высокотемпературная спонтанная намагниченность исчезает:
Помимо, расширение кластера показывает, что спиновые корреляции группируются экспоненциально быстро: например,
При низких температурах, т.е. β ≫ 1, спонтанная намагниченность остается нулевой (см. Теорема Мермина – Вагнера ),
но затухание корреляций происходит только по степенному закону: Фрёлих и Спенсер[13] нашел нижнюю границу
а Макбрайан и Спенсер нашли верхнюю границу для любых
Три и выше измерения
Независимо от диапазона взаимодействия при достаточно низкой температуре намагниченность положительна.
- При высокой температуре спонтанная намагниченность исчезает: . Помимо, расширение кластера показывает, что спиновые корреляции группируются экспоненциально быстро: например, .
- При низкой температуре инфракрасный порт показывает, что спонтанная намагниченность строго положительна: . Кроме того, существует однопараметрическое семейство экстремальных состояний: , так что но предположительно в каждом из этих экстремальных состояний усеченные корреляции затухают алгебраически.
Фаза перехода
Как упоминалось выше, в одном измерении модель XY не имеет фазового перехода, в то время как в двух измерениях она имеет Переход Березинского-Костерлица-Таулеса между фазами с экспоненциально и степенно убывающими корреляционными функциями.
В трех и более измерениях XY-модель имеет фазовый переход ферромагнетик-парамагнетик. При низких температурах спонтанная намагниченность отлична от нуля: это ферромагнитная фаза. С повышением температуры спонтанная намагниченность постепенно уменьшается и исчезает при критической температуре. Он остается нулевым при всех более высоких температурах: это ферромагнитная фаза.
В четырех измерениях и выше фазовый переход имеет средние критические показатели теории поля (с логарифмическими поправками в четырех измерениях).
Трехмерный случай: критические показатели
Трехмерный случай интересен тем, что критические показатели при фазовом переходе нетривиальны. Многие трехмерные физические системы принадлежат одному и тому же класс универсальности как трехмерная модель XY и имеют одни и те же критические показатели, в первую очередь магниты с легкой плоскостью и жидкие Гелий-4. Ценности этих критические показатели измеряются с помощью экспериментов, моделирования Монте-Карло, а также могут быть рассчитаны теоретическими методами квантовой теории поля, такими как ренормгруппа и конформный бутстрап. Методы ренормгруппы применимы, потому что считается, что критическая точка модели XY описывается фиксированной точкой ренормгруппы. Конформные методы начальной загрузки применимы, потому что они также считаются унитарным трехмерным конформная теория поля.
Самое важное критические показатели трехмерной XY-модели . Все они могут быть выражены всего двумя числами: масштабные размеры и поля комплексного параметра порядка и ведущего синглетного оператора (такой же как в Гинзбург – Ландау описание). Еще одно важное поле - это (такой же как ), размерность которой определяет степень поправки к масштабированию . Согласно вычислению конформного бутстрапа,[14] эти три измерения задаются:
0.519088(22) | |
1.51136(22) | |
3.794(8) |
Это дает следующие значения критических показателей:
общее выражение () | численная величина | |
---|---|---|
α | -0.01526(30) | |
β | 0.34869(7) | |
γ | 1.3179(2) | |
δ | 4.77937(25) | |
η | 0.038176(44) | |
ν | 0.67175(10) | |
ω | 0.794(8) |
Методы Монте-Карло дают совместимые определения:[15] .
Смотрите также
- Классическая модель Гейзенберга
- Бозон Голдстоуна
- Модель Изинга
- Модель Поттса
- п-векторная модель
- Переход Костерлица – Таулеса
- Топологический дефект
- Сверхтекучая пленка
- Сигма модель
Примечания
- ^ Стэнли, Е. (1968). «Зависимость критических свойств от размерности спинов». Phys. Rev. Lett. 20 (12): 589–592. Bibcode:1968ПхРвЛ..20..589С. Дои:10.1103 / PhysRevLett.20.589.
- ^ Чайкин П.М.; Лубенский, Т. (2000). Принципы физики конденсированного состояния. Издательство Кембриджского университета. ISBN 978-0521794503.
- ^ Жинибр Дж. (1970). «Общая формулировка неравенств Гриффитса». Comm. Математика. Phys. 16 (4): 310–328. Bibcode:1970CMaPh..16..310G. Дои:10.1007 / BF01646537.
- ^ Aizenman, M .; Саймон, Б. (1980). «Сравнение ротора самолета и модели Изинга». Phys. Lett. А. 76 (3–4): 281–282. Bibcode:1980ФЛА ... 76..281А. Дои:10.1016/0375-9601(80)90493-4.
- ^ Бадалян, Д. (1996). «К термодинамике классических спинов с изотропным гейзенберговским взаимодействием в одномерных квазипериодических структурах». Physica B. 226: 385–390. Дои:10.1016/0921-4526(96)00283-9.
- ^ Мэттис, округ Колумбия (1984). «Передаточная матрица в модели плосковращателя». Phys. Латыш. 104 А (6–7): 357–360. Bibcode:1984ФЛА..104..357М. Дои:10.1016/0375-9601(84)90816-8.
- ^ Мэттис, Д. К. (1985). Теория магнетизма II. Серия Спрингера по физике твердого тела. ISBN 978-3-642-82405-0.
- ^ Ota, S .; Ota, S.B .; Фанле, М. (1992). «Микроканоническое моделирование методом Монте-Карло для двумерной XY-модели». J. Phys .: Condens. Иметь значение. 4: 5411. Дои:10.1088/0953-8984/4/24/011.
- ^ Hsieh, Y.-D .; Kao, Y.-J .; Sandvik, A.W. (2013). "Метод конечных размеров для перехода Березинского-Костерлица-Таулеса". J. Stat. Механизм .: Theory Exp. 2013. arXiv:1302.2900. Дои:10.1088 / 1742-5468 / 2013/09 / P09001.
- ^ Tobochnik, J .; Честер, Г.В. (1979). "Монте-Карло исследование модели плоского спина". Phys. Ред. B. 20 (9): 3761–3769. Дои:10.1103 / PhysRevB.20.3761.
- ^ Биндер, К. (2013). Применение метода Монте-Карло в статистической физике. Springer Science & Business Media. ISBN 978-3-642-51703-7.
- ^ Van Himbergen, J.E .; Чакраварти, Судип (1981). «Модуль спиральности и теплоемкость классической модели XY в двух измерениях». Phys. Ред. B. 23: 359. Дои:10.1103 / PhysRevB.23.359.
- ^ Fröhlich, J .; Спенсер, Т. (1981). «Переход Костерлица – Таулеса в двумерных абелевых спиновых системах и кулоновский газ». Comm. Математика. Phys. 81 (4): 527–602. Bibcode:1981CMaPh..81..527F. Дои:10.1007 / bf01208273.
- ^ Честер, Шай М .; Лэндри, Уолтер; Лю, Цзюнюй; Польша, Давид; Симмонс-Даффин, Дэвид; Су, Нин; Вичи, Алессандро. «Выделение пространства OPE и точных критических показателей модели O (2)». Журнал физики высоких энергий. 2020 (6): 142. arXiv:1912.03324. Дои:10.1007 / JHEP06 (2020) 142. ISSN 1029-8479.
- ^ Хазенбуш, Мартин (26 декабря 2019 г.). «Монте-Карло исследование улучшенной модели часов в трех измерениях». Физический обзор B. 100 (22): 224517. arXiv:1910.05916. Дои:10.1103 / PhysRevB.100.224517. ISSN 2469-9950.
Рекомендации
- Евгений Демидов, Вихри в XY-модели (2004)
дальнейшее чтение
- Х. Э. Стэнли, Введение в фазовые переходы и критические явления(Издательство Оксфордского университета, Оксфорд и Нью-Йорк, 1971);
- Х. Кляйнерт, Калибровочные поля в конденсированных средах, Vol. I, "СУПЕРПОТОК И ВИХРЕВЫЕ ЛИНИИ", стр. 1–742, Vol. II, «НАПРЯЖЕНИЯ И ДЕФЕКТЫ», стр. 743–1456, World Scientific (Сингапур, 1989 г.); Мягкая обложка ISBN 9971-5-0210-0 (также доступно в Интернете: Vol. я и Vol. II )