Исследование модели Изинга на гексагональной решетке методом Ванга-Ландау Текст научной статьи по специальности «Физика»

УДК 537.9

Б01: 10.21779/2542-0321- 2016-31-4-71-76 М.А. Магомедов1'2, А.К. Муртазаев1'2, Л.К. Магомедова1'2 Исследование модели Изинга на гексагональной решетке методом Ванга-Ландау

1 Институт физики Дагестанского научного центра РАН; Россия, 367003, г. Махачкала, ул. М. Ярагского, 94;

2Дагестанский государственный университет; Россия, 367001, г. Махачкала, ул. М. Гаджиева 43а; magomedov_ma@mail.ru

Антиферромагнитная модель Изинга на гексагональной решетке исследована высокоэффективным алгоритмом Ванга-Ландау метода Монте-Карло. Вычислена плотность состояний системы и рассчитаны температурные зависимости внутренней энергии и, свободной энергии ¥, энтропии Б, теплоемкости С и кумулянтов Биндера икт. Показано отсутствие вырождения основного состояния и наличие фазового перехода второго рода. Определена температура фазового перехода Тм = 1.5.

Ключевые слова: модель Изинга, плотность состояний, энтропия, фазовые переходы, метод Монте-Карло, алгоритм Ванга-Ландау.

В последние десятилетия значительное внимание уделяется экспериментальному и теоретическому исследованию различных низкоразмерных систем, в том числе наноструктур. Изучение различных свойств этих объектов открывает широкие перспективы для экспериментальных приложений. В ближайшие годы нанотехнологии позволят совершить поистине гигантский технологический скачок в самых различных областях науки и техники [1-6].

Одними из интересных низкоразмерных систем являются графен и графенопо-добные наноструктуры, обладающие целым рядом уникальных свойств. Данные нано-системы обладают структурой, изображенной на рисунке 1. Нами проведены высокоточные исследования двумерной модели Изинга на гексагональной решетке. Данная модель достаточно хорошо описывает графен и подобные ему наноситемы [6, 7].

Схематически модель представлена на рисунке 1. Как видно из рисунка, спины в узлах решетки располога-ются в двух подрешетках. Атомы, входящие в одну подрешетку А, обозначены на рисунке одним цветом, а атомы, входящие во вторую подрешетку В, -другим цветом. Каждый узел имеет трех ближайших соседей [7].

Гамильтониан модели может быть представлен в следующем виде:

Н = - 7 X 8 ,8 , ,

(,, 7> 8,. =±1.

(1)

Рис. 1. Модель Изинга на гексагональной решетке с учетом взаимодействия между первыми и вторыми ближайшими соседями.

При положительном обменном взаимодействии / все атомы в основном состоянии упорядочиваются одинаково, а при отрицательном значении / атомы образуют две подрешетки, спины в ко-

торых направлены в основном состоянии антипараллельно. В данной модели нет вырождения и отсутствует фрустрация. Энтропия такой системы в основном состоянии должна быть равна нулю. Мы исследовали в данной работе антиферромагниную модель Изинга (/ < 0).

В последние годы низкоразмерные структуры интенсивно исследуются с помощью энтропийных методов Монте-Карло, одним из вариантов которого является алгоритм Ванга-Ландау [7-10]. Данный алгоритм позволяет найти функцию плотности состояний системы, зная которую можно легко рассчитать все остальные характеристики системы. Особенно эффективным алгоритм Ванга-Ландау оказался при моделировании различных наноструктур.

Алгоритм Ванга-Ландау основан на том, что, совершая случайное блуждание в пространстве энергий с вероятностями, обратно пропорциональными плотности состояний g(Е), мы получаем равномерное распределение по энергиям. Подобрав такие вероятности перехода, при которых посещение всех энергетических состояний стало бы равномерным, можно получить изначально неизвестную плотность состояний g(E), зная которую, можно вычислить значения необходимых термодинамических параметров при любой температуре. Так как плотность состояний g (Е) очень быстро растет с увеличением размеров исследуемых систем, для удобства хранения и обработки больших чисел пользуются величиной 1п g(E).

Важным обстоятельством является то, что плотность состояний g (Е) не зависит от температуры, следовательно, рассчитав ее однократно, мы можем вычислить значения любых термодинамических параметров системы при любой температуре.

В данной работе алгоритм Ванга-Ландау был использован нами в следующем виде [7].

• Задается произвольная начальная конфигурация спинов. Стартовые значения плотности состояний g(E) = 1 , гистограммы распределений по энергиям Н(Е) = 0 ,

стартовый модификационный фактор / = /0 = е1 « 2.71828.

• Многократно совершаем шаги в фазовом пространстве пока не получим относительно плоскую гистограмму Н (Е) (т. е. пока не будут посещены примерно одинаковое количество раз все возможные энергетические состояния системы). В качестве критерия "плоскости" гистограммы нами принималось условие отклонения числа посещений всех возможных (с ненулевой плотностью g(E') ^ 1) энергетических состояний не более чем на 10 % от среднего значения по системе.

• При этом вероятность перехода из состояния с энергией Е1 в состояние с энергией Е2 определяется по формуле р = g(E1 )/g(E2). Если переход в состояние с энергией Е2 состоялся, то для энергии Е2 проводится модификация плотности состояния g(E2 / х g(E2), и гистограммы Н (Е2 Н (Е2)+1 иначе меняем параметры для

дификационный фактор / > /тп . В качестве минимального значения модификационно-го фактора нами принималось /т-п = 1.0000000001.

• После расчета плотности состояний системы для любой интересующей нас температуры расчитываются различные термодинамические параметры, такие, как энтропия, внутренняя энергия, свободная энергия, теплоемкость, намагниченность, вос-

приимчивость и т. д. Некоторые формулы (2-6) для расчета термодинамических параметров приведены ниже.

Более подробно алгоритм Ванга-Ландау изложен в работах [8-10].

400

300

UJ

Ö) с

200

100

—■— L = 6

—•— L = 12

—'»— L = 18

—т— L = 24

Рис. 2. Плотность состояний g(E) для модели Изинга на гексагональной решетке при различных линейных размерах системы

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

Зная плотность состояний системы, можно рассчитать любой интересующий нас термодинамический параметр A(T) для любой температуры по следующей формуле:

X Ag (E )e - E 7 kBT

A(T) = -

X g (E )e -

E

Например, внутрення энергия U, свободняя энергия F, энтропия S и теплоемкость C системы могут быть рассчитаны следующим образом:

X Eg (E )e - E 7

U (T) = E

(2)

I g (E )e -

=< E >

E

F (T )=-kBT Infi g (E >

- ЩквТ

\

E

(3)

У

5 (T )_

U (T )- F (t )

T

C (T )_

_ eu (T ) _ (E 2) -( E

dT

kBT2

(5)

(6)

Рассчитанные из плотности состояний температурные зависимости различных термодинамических параметров приведены на рисунке 3. Отметим важную особенность алгоритма Ванга-Ландау: значения любых термодинамических параметров можно определить для любой температуры, с любым шагом, при этом объем необходимых вычислений, в отличие от других классических алгоритмов метода Монте-Карло, вырастает незначительно.

Температурные зависимости кумулянтов Биндера по параметру порядка и энтропии приведены на рисунке 4. Как видно из рисунка, при TN = 1.5 в системе происходит фазовый переход второго рода. Энтропия при высоких температурах стремится к значению ln2, при низких температурах - к нулю (отсутствует вырождение), а в точке фазового перехода испытывает перегиб и не зависит от линейных размеров системы.

ÙJ

-1.2 -

—— 1 = 6

—•— L = 12

L = 18

L = 24

—— L = 6

—•— L = 12

—*— L = 18

—.— L = 24

—L = 6

—■— L = 12

L = 18

—T— L = 24

О

Т Т

Рис. 3. Температурные зависисмости внутренней энергии Е, теплоемкости С, параметра

порядка q и восприимчивости %, рассчитанные из плотности состояний g(E)

2

1-■-1-'-1-■-1-■- 0 8 -.-1-1-1-1-.-.-1-

Т т

Рис. 4. Температурные зависисмости кумулянтов Биндера Ukm и энтропии S

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

Работа выполнена при финансовой поддержке проектной части госзадания Мино-брнауки России в научной деятельности (проект 3.1262.2014К).

Литература

1. Nakano A., Bachlechner M.E., Kalia R.K. et al. Multiscale Simulation of Nanosystems // Computing in Science and Engg. - 2001. - V. 3 (4). - P. 56-66.

2. Cheng C., Bukkapatnam S.T.S., Raff L.M., et al. Monte Carlo simulation of carbon nanotube nucleation and growth using nonlinear dynamic predictions // Chem. Phys. Lett. -2012. - V. 530 (0). - P. 81-85.

3. De S., Wen X., Bordas S.P.A. et al. Defect engineering of 2d monatomic-layer materials // Modern Physics Letters B. - 2013. - V. 27 (23). - P. 1330017.

4. Liu Z.S., Sechovsky V.Divis M. A new combined quantum simulation approach for nanomagnets // Physica E. - 2013. - V. 47 (0). - P. 128-133.

5. Fabritius T., Laflorencie N. Wessel S. Finite-temperature ordering of dilute graphene antiferromagnets // Phys. Rev. B. - 2010. - V. 82 (3). - P. 035402.

6. Masrour R., Bahmad L.Benyoussef A. Size effect on magnetic properties of a nano-graphene bilayer structure: A Monte Carlo study // JMMM. - 2012. - V. 324 (23). - P. 3991-3996.

7. Магомедов М.А., Муртазаев А.К., Магомедова Л.К. Плотность состояния и энтропия модели Изинга на гексагональной решетке с учетом взаимодействия вторых ближайших соседей // Вестник ДГУ. - 2016. - Т. 31, вып. 1. - С. 43-50.

8. Shell M.S., Debenedetti P.G.Panagiotopoulos A.Z. Generalization of the WangLandau method for off-lattice simulations // Phys. Rev. E. - 2002. - V. 66 (5). - P. 56703.

9. Landau D.P., Tsai S.-H.Exler M. A new approach to Monte Carlo simulations in statistical physics: Wang-Landau sampling // Am. J. Phys. - 2004. - V. 72 (10). - P. 1294-1302.

10. Zhou C.Bhatt R.N. Understanding and improving the Wang-Landau algorithm // Physical Review E. - 2005. - V. 72 (2). - P. 025701.

Поступила в редакцию 14 ноября 2016 г.

UDC 537.9

DOI: 10.21779/2542-0321-2016-31-4-71-76

The study of the Ising model on the hexagonal lattice by Wang-Landau method

M.A. Magomedov1,2, A.K. Murtazayev1,2, L.K. Magomedova1,2

institute of Physics DSC RAS; Russia, 367003, Makhachkala, Yaragskу st., 94; 2Dagestan State University; Russia, 367025, Makhachkala, M. Gadzhiyev st., 43 a; magomedov_ma@mail.ru

The Antiferromagnetic Ising model on the hexagonal lattice is investigated by high-performance Wang-Landau algorithm of the Monte Carlo method. The density of states of the system are calculated and the temperature dependences of internal energy U, free energy F, entropy S, the heat capacity C and Binder cummulant Ukm. The absence of the ground state degeneracy is shown. A second order phase transition On TN = 1.5 is also determined.

Keywords: Ising model, density of states, entropy, phase transition, Monte Carlo method, Wang-Landau algorithm.

Received 14 November, 2016