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

Материал из Википедии — свободной энциклопедии

Модель Изинга — математическая модель статистической физики, предназначенная для описания намагничивания материала.

Каждой вершине кристаллической решётки (рассматриваются не только трёхмерные, но и одно- и двумерные случаи) сопоставляется число, называемое спином и равное +1 или −1 («поле вверх»/«поле вниз»). Каждому из 2^N возможных вариантов расположения спинов (где N — число атомов решётки) приписывается энергия, получающаяся из попарного взаимодействия спинов соседних атомов:

где J — энергия взаимодействия (в простейшем случае одна и та же для всех пар соседних атомов). Иногда также рассматривается внешнее поле h (часто полагаемое малым):

Затем, для заданной обратной температуры β = 1 / k_BT на получившихся конфигурациях рассматривается распределение Больцмана: вероятность конфигурации полагается пропорциональной , и исследуется поведение такого распределения при очень большом числе атомов N.

Например, в моделях с размерностью, большей 1, имеет место фазовый переход: при достаточно низких температурах большая часть спинов ферромагнетика (J>0) будет ориентирована (с близкой к 1 вероятностью) одинаково, а при высоких почти наверняка спинов «вверх» и «вниз» будет почти поровну. Температура, при которой происходит этот переход (иными словами, при которой исчезают магнитные свойства материала), называется критической, или точкой Кюри. В окрестности точки фазового перехода(в данном случае второго рода) ряд термодинамических характеристик расходится. Опыт показывает, что расходимость имеет универсальный характер, и определяется лишь симметрией системы. Впервые критические индексы расходимостей были получены для двумерной модели Изинга в 40-х годах Онсагером. Для остальных размерностей исследования проводится с помощью методов компьютерного моделирования, ренормгруппы. Обоснованием применения ренормализационной группы в данном случае являются блочное построение Каданова и термодинамическая гипотеза подобия.

[править] Алгоритм реализации модели Изинга методом Монте-Карло на компьютере

1. Создать решётку спинов (двумерный массив), спины ориентированы произвольно.
2. Выбрать случайно одну из клеток решётки, стереть значение в ней.
3. Вычислить энергии конфигураций при заполнении этой клетки спином вверх и вниз (либо при всех возможных состояниях, если их больше двух).
4. Выбрать один из вариантов для «стёртого» спина случайно, с вероятностью, пропорциональной , где E(S) — энергия в соответствующем состоянии. (Поскольку все слагаемые, не затрагивающие данный спин, одни и те же, на самом деле вычислять нужно только суммы по соседям).
5. Возвращаемся в пункт 2; по выполнении достаточного числа итераций (определение этого — отдельная и непростая задача) цикл прекращается.

Это незавершённая статья по физике. Вы можете помочь проекту, исправив и дополнив её.