ЗАДАЧА 15.6. Профиль намагниченности

а.         Модифицируйте программу conduct, удалив подпрограмму heat иограничив спины 1 и N значениями +1 н -1 соответственно. Оцени-те профиль намагниченности, построив графики среднего значенияспииа в каждом узле как функции номера узла. Возьмите параметрыN = 22, nmcs = 1000 и У = 1. Как меняются ваши результаты с рос-том N?

б.         Вычислите среднюю энергию демона и по ней локальную темпера-туру в каждом узле. Однородна ли температура системы, даже еслинамагниченность неоднородна? Находится ли система в тепловомравновесии?

в.         Влияние этого ограничения легче наблюдать в двумерном и трех-мерном случаях, чем в одномерном. Составьте программу для дву-мерной модели Изинга на квадратной решетке Lx L. Ограничьтеспины в узле (»,/) значениями +1 и -1 для i = 1 и i = L соответ-ственно. Используйте периодические краевые условия по у. Как со-относятся ваши результаты с одномерным случаем?

г.         Откажитесь от периодических краевых условий по у и зафикси-руйте все граничные спины от i = 1 до L/2 равными +1, а все ос-тальные граничные спины равными -1. Возьмите начальную конфигу-рацию, где все спины в левой половине системы равны +1, а всеостальные -1. Проведите моделирование и изобразите конфигурациюспинов после достижения системой равновесия. Соедините на схемекаждую пару спинов противоположного знака. Опишите кривую, раз-деляющую спины +1 и спины -1. Опишите вид этой кривой. Начните сL = 20 и установите, что будет происходить при увеличении L.