ЗАДАЧА 15.4. Двумерная модель Изинга

а.         Обобщите программу Ising_demon на случай двумерной моделиИзинга на квадратной решетке с динамикой опрокидывания спинов.Возьмите решетку со стороной,   равной  L.  Полное число спинов N

о

равняется N = L . Используйте периодические краевые условия, как показано на рис. 15.2, так чтобы спины в левом столбце взаимодействовали со спинами в правом столбце и т. д.

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

о

среднюю энергию демона и <М > как функции от Е. Начальные параметры удобно выбрать равными L = 5 и h = 0. С помощью формулы (15.8) определите зависимость температуры от энергии системы.

в.         Модифицируйте свою программу, чтобы сделать несколько «сним-ков» конфигураций спинов. Опишите качественно вид этих конфигу-раций при различных энергиях и температурах. Упорядочены они илинет? Имеются ли домены спинов, ориентированных только вверх иливниз?