CC BY
22
структуры и моделирование 2017. №1(41). С. 30-37
УДК 538.91
КРИТИЧЕСКОЕ ПОВЕДЕНИЕ ДВУМЕРНОЙ НЕУПОРЯДОЧЕННОЙ МОДЕЛИ ИЗИНГА С ДЕФОРМИРОВАННОЙ СТАТИСТИКОЙ
В.Н. Бородихин
доцент, к.ф.-м.н., e-mail: [email protected]
Омский государственный университет им. Ф.М. Достоевского
Аннотация. В работе впервые исследована двумерная неупорядоченная модель Изинга с деформированной статистикой Цаллиса, со спиновыми концентрациями 0.95 и 0.8. Получены значения критических температур и критических показателей. Для неупорядоченной модели с деформированной статистикой выявлено возникновение нового типа критического поведения, зависящего от концентрации примесей.
Ключевые слова: фазовые переходы, деформированная статистика, модель Изинга.
Введение
В работе [1] рассматривали двумерную однородную модель Изинга с деформированной статистикой Цаллиса [2]. Были определены критические температуры и критические показатели для систем с различной степенью деформации статистики. Исследование данной модели представляет теоретический и экспериментальный интерес, поскольку критическое поведение данной модели качественно согласуется с критическим поведением реальных систем, в частности, так называемых манганитов, например, La0mY0mCa0.33MnO3 [3,4]. Однако вопрос о связи подобных систем и неэкстенсивной статистики нуждается в уточнении.
Деформированная статистика характеризуется параметром деформации q, связанным со степенью отклонения от статистики Гиббса. В [1] было показано, что в отличие от стандартной однородной двумерной модели Изинга, в которой индекс теплоёмкости а = 0, в двумерной однородной модели Изинга с деформированной статистикой индекс а является положительным. Согласно критерию Харриса [5] немагнитные примеси оказывают существенное влияние на критическое поведение в случае а > 0, и таким образом, представляется интересным исследовать влияние вмороженного беспорядка в рамках двумерной модели Изинга с деформированной статистикой, а также исследовать влияние деформированной статистики на критерий Харриса.
В работе впервые были определены значения критических температур и критических показателей для двумерной неупорядоченной модели Изинга со спиновой концентрацией p = 0.95 с деформированной статистикой Цаллиса.
1. Моделирование систем с деформированной статистикой
Выражение для неаддитивной энтропии Цаллиса имеет вид [2]:
1 -Е Р
Бя = к-, (1)
д - 1
где рг — вероятность нахождения системы в состоянии 1, д — параметр деформации статистики: при 0 < рг < 1 р^ > рг для д < 1, и р^ < рг для д > 1. При д =1 статистика не является деформированной.
В неэкстенстивной статистической теории [2] вводится определение:
п
<Н> = ^ Ргег = и, (2)
г=1
где Н — гамильтониан системы, ег — одно из возможных энергетических состояний, и — внутренняя энергия, Рг так называемое эскортное распределение [6]
г, _ р1 _ 4£г)я (3)
р = ~п = ~п : , (3)
Ер9 Е^-Ч)*
3 3
где
-' = - л---, (4)
Т р + (1 - д )ви
3
здесь - — лагранжев множитель, е-х — деформированная экспонента, удовлетворяющая свойству
[1 - (1 - д)-'вг]1/(1-) = ( 0-X, 1 -(1 ~д)< ^ (5)
10,1 - (1 - д)х < 0.
Следуя [1], в качестве физической температуры выберем -'. При компьютерном моделировании системы использовали алгоритм метро-полиса с модифицированной вероятностью переворота спинов, обусловленной деформированной статистикой:
Р(2) ~42)/т ^ = роУ = , (6)
г 6д
где Рг(2) эскортное распределение (3), соответствующее состоянию после переворота спина, а Рг(1) - состоянию до переворота. Состояние энергии ег определяется гамильтонианом двумерной неупорядоченной модели Изинга с учётом взаимодействия ближайших соседних спинов вг:
Н =2 Е ЗцРгРз , (7)
г,3
где — константа обменного ферромагнитного взаимодействия, рг — случайная переменная, описываемая функцией распределения
РЫ = р6(Рг - 1) + (1 - РгЩрг) (8)
с р = 1 — с, где с — концентрация примесей. Система имеет вид квадратной решётки размера Ь с периодическими граничными условиями. С узлами решётки связаны спины вг, принимающие значения ±1 и немагнитные атомы примеси (пустые узлы с вг = 0). Примесь равномерно распределяется по всей системе, её положение фиксировано для отдельной примесной конфигурации в процессе моделирования системы. Концентрация спинов определяется суммированием абсолютных значений спинов по всем узлам решётки:
1 ь2
р = Ь Е N. (9)
г=1
При моделировании системы вычисляли намагниченность на спин т, энергию на спин е, восприимчивость х
х = (т)т(т2 • <10>
где (...) означает статистическое усреднение по шагам Монте-Карло, а черта — усреднение по примесным конфигурациям. Также вычисляли теплоёмкость
С = (11)
Таблица 1. Критические индексы и критические температуры двумерной однородной и неупорядоченной модели Изинга (р=1, 0.95) с деформированной статистикой
р V Т ± с а/и в/и 1/и
1 1 2.269 0 0.125 1.75
0.8 1.888 ± 0.005 0.223 ± 0.007 0.1247 ± 0.003 1.770 ± 0.042
0.6 1.773 ± 0.007 0.261 ± 0.019 0.121 ± 0.003 1.696 ± 0.01
0.95 1[7] 2.088 0.125(5) 1.75(2)
0.8 1.743 ± 0.006 0.166 ± 0.023 0.1063 ± 0.01 1.605 ± 0.04
0.6 1.662 ± 0.006 0.182 ± 0.022 0.135 ± 0.004 1.650 ± 0.05
и кумулянт Биндера
V =1 - ■ (12)
3(т2)2
Отношение критических индексов теплоёмкости а/у, намагниченности -/V и восприимчивости определяли из зависимости параметров от размеров системы Ь в критической точке Тс бесконечной системы
С(Ь,Т) - Ьа/у; (13)
т(Ь,Т) - Ь~в/и;
х(ь, Т) - Ь^,
где V - критический индекс корреляционной длины. Для двумерной однородной модели Изинга он равен 1, и как показано в работе [1], практически не зависит от степени деформации статистики.
Таблица 2. Критические индексы и критические температуры двумерной неупорядоченной модели Изинга (р=0.8) с деформированной статистикой
р д Т Тс -/V
0.8 1[7] 1.5077 0.11(1) 1.78(2)
0.8 1.2905 ± 0.009 0.126 ± 0.028 1.3474 ± 0.238
0.6 1.2413 ± 0.01 0.1185 ± 0.02 1.332 ± 0.225
2. Результаты моделирования
Было осуществлено моделирование двумерной однородной и неупорядоченной модели Изинга с деформированной статистикой. Спиновая концентрация неупорядоченной системы составляла р = 0.95, р = 0.8, параметр деформации д = 0.8,0.6, размеры системы выбирались Ь = 256,128,64, 32. При моделировании использовали 500-1000 примесных конфигураций для неупорядоченных систем, для однородных систем 100-400 прогонок. Отношение критических индексов а^, -/V и определяли из соотношений (13) в двойном логарифмическом масштабе в критической точке. Критические температуры определяли по пересечению кумулянтов Биндера (12). На рис. 1 для примера приведено пересечение кумулянтов (р=0.95, q=0.6). На рис. 2-5 приведены графики критических показателей при различных концентрациях и параметрах деформации.
Результаты моделирования приведены в табл. 1. Для однородной модели в табл. 1 в отсутствие деформации ^=1) приведены точные значения, для неупорядоченной (р=0.95) и q=1, результаты взяты из [7]. Для спиновой концентрации р = 0.8 результаты моделирования приведены в табл. 2.
Рис. 2. Критический показатель 7, р=1.0, ц=0.6.
Рис. 4. Критический показатель а, p=0.95, q=0.6.
Рис. 5. Критический показатель 7, p=0.8, q=0.8.
3. Заключение
Сравнивая результаты данной работы для однородной деформированной модели Изинга и работы [1], можно выявить некоторые отличия. Критические температуры, полученные в рамках методики пересечения кумулянтов Биндера, в пределах погрешности совпадают с [1], как и критические показатели 7. Критические показатели а почти совпадают в пределах погрешности. Что касается критического показателя в, то в [1] выявлена сильная зависимость данного показателя от параметра деформации статистики q. Так для д = 0.8 в = 0.075, а для д = 0.6 в = 0.025, т.е. с уменьшением параметра деформации критический индекс в практически стремится к нулю. В то же время в данной работе не выявлено подобной зависимости. Критические показатели определяли из логарифмической зависимости соответствующих термодинамических величин в критической точке бесконечной системы.
Как известно, критические показатели двумерной неупорядоченной модели Изинга практически не зависят от концентрации примесей. Однако же в рамках деформированной неупорядоченной модели можно сказать, что и при малых концентрациях примесей (0.05), согласно значениям критических показателей а и 7, появляется новый тип критического поведения. Данный вывод подтверждает критерий Харриса (а > 0) соответствующей однородной модели с деформированной статистикой д = 1. Что касается примесных концентраций 0.2, то здесь не удалось определить показатель а с хорошей точностью. Тем не менее в частности индекс 7 демонстрирует существенное отличие от недефор-мируемой модели, а также от деформируемой модели с малой концентрацией
примесей (0.05). Следовательно, для данной примесной концентрации также выявлен новый тип критического поведения.
Таким образом, влияние деформации статистики Цаллиса для неупорядоченной модели Изинга приводит к появлению нового типа критического поведения, зависящего от концентрации примесей.
Литература
1. Crokidakis N., Soares-Pinto D.O., Reis M.S., Souza A.M., Sarthour R.S., Oliveira I.S. Finite size analysis of a two-dimensional Ising model within a nonextensive approach // Phys. Rev. E. 2009. V. 80. P. 051101.
2. Tsallis C. Introduction to Nonextensive Statistical Mechanics: Approaching a Complex World. New York : Springer, 2009. P. 382.
3. Reis M.S., Amaral V.S., Araujo J.P. and Oliveira I.S. Magnetic phase diagram for a non-extensive system: Experimental connection with manganites // Phys. Rev. B. 2003. V. 68. P. 014404.
4. Amaral V.S., Araujo J.P., Pogorelov Y.P., Tavares P.B., Souza A.M., Vieira J.M. Discontinuous transition effects in manganites: magnetization study in the paramagnetic phase // J. Magn. Magn. Mater. 2002. V. 242. P. 655-658.
5. Harris A.B. Effect of random defects on the critical behaviour of Ising models // J. Phys. C. 1974. V. 7. P. 1671.
6. Beck C., Schlogl F. Thermodynamics of Chaotic Systems: An Introduction. Cambridge, : Cambridge University Press, 1993. P. 308.
7. Heuer H.O. Monte Carlo Simulation of Disordered 2-Dimensional Ising Systems // Europhys. Lett. 1991. V. 16. P. 503-508.
STUDY OF THE SYSTEMS BEHAVIOR IN PHASE TRANSITIONS ON THE COALESCENCE STAGE
V.N. Borodikhin
Associate Professor, Ph.D. (Phys.-Math.), e-mail: [email protected] Dostoevsky Omsk State University
Abstract. For the first time the two-dimensional disordered Ising model with deformed Tsallis statistics, with spin concentrations of 0.95 and 0.8 is investigated. The values of critical temperatures and the critical exponents are obtained. For disordered model with deformed statistics the emergence of a new type of critical behavior depending on the concentration of impurities is revealed.
Keywords: phase transitions, deformed statistics, Ising model.
Дата поступления в редакцию: 11.12.2016
