Критическая температура и размер решетки с помощью алгоритма Вольфа для двумерной модели Изинга

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

Сначала я подумал, что это ошибка в программе, но поведение действительно имеет смысл. На бесконечной решетке алгоритм Вольфа должен создавать кластеры всех размеров на Т С . Это означает, что большинство кластеров, которые он пытается создать, больше, чем решетка, используемая в моделировании, и большинство кластеров в конечном итоге достигают границ и заполняют большую часть решетки. Это также восходит к точке зрения Каданова, что истинная критичность возможна только в бесконечных системах.

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

Есть ли литература по этому эффекту?

Как люди компенсируют это на практике?

Есть ли формула для регулировки температуры для различных размеров решетки?

Это хорошо понятный эффект конечного размера, возникающий из-за того, что длина корреляции в вашей системе не может стать больше, чем размер системы. Это подразумевает сдвиг (кажущейся) критической точки порядка 1 / н , если вы находитесь на н × н тор. Насколько мне известно, это впервые было изучено Фердинандом и Фишером в 1969 году (Phys. Rev. 185, 832).
@YvanVelenik: Если бы вы могли немного рассказать об этом, это был бы хороший ответ!
@ACuriousMind: Сделано (но это немного схематично, так как у меня сейчас не так много времени; в частности, было бы неплохо представить картину теплоемкости для 2d-модели Изинга на конечном торе, а также как предельное количество).

Ответы (1)

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

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

Теперь, конечно, достаточно большая конечная система будет по-прежнему демонстрировать поведение, «напоминающее» фазовый переход, но со сглаженными сингулярностями. Затем для экстраполяции результатов на бесконечные системы требуется (i) определение аналогов предельных величин конечного объема (в частности, критической температуры), (ii) изучение того, как эти величины конечного объема изменяются при увеличении размера системы. Чтобы помочь с этой процедурой экстраполяции, физики разработали различные теории масштабирования конечного размера.

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

Первая подробная теория масштабирования конечного размера была разработана Фердинандом и Фишером в 1969 году в классической статье, опубликованной в журнале Phys. Откр. 185, 832 . Они использовали точные результаты, доступные для двумерной модели Изинга, для анализа эффектов конечного размера на свободную энергию и удельную теплоемкость.

Удельная теплоемкость модели Изинга конечного объема не расходится. Однако он по-прежнему показывает резкое увеличение в узкой области вокруг «истинной» критической точки. Т с . Фишер и Фердинанд предложили определить аналог конечного объема Т с ( л ) критической температуры как значение температуры, при которой теплоемкость максимальна. Затем они утверждали, что

Т с Т с ( л ) л 1 / ν ,
где ν - критический показатель, связанный с удельной теплоемкостью, которая определяется выражением ν "=" 1 для двухмерной модели.

Спасибо. Это именно то, что я искал. Т с оценки, которые я получил, дают константу пропорциональности а * почти точно в 4 раза от Фердинанда и Фишера -.3603. Я выбирал температуру, глядя на логарифмические графики распределения размеров кластеров и выбирая те, которые давали самую прямую линию возможных размеров. Не очень строго, но я получаю соответствие 3 значащим цифрам с их формулой для размеров решетки 64, 127, 512 после корректировки с коэффициентом 4. Интересно, почему коэффициент 4.