Одномерная модель Изинга с подвижными примесями Текст научной статьи по специальности «Математика»

УДК 531.19

с.В. сёмкин1 В.П. смагин2

Владивостокский государственный университет экономики и сервиса Владивосток. Россия

Одномерная модель Изинга с подвижными примесями*

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

Ключевые слова и словосочетания: фазовые переходы, модель Поттса, решетка Бете.

S.V. Semkin V.P. Smagin

Vladivostok State University of Economics and Service Vladivostok. Russia

The one-dimensional Ising model with mobile impurities

We solve the one-dimensional Ising model with movable magnetic impurities. For this model, we found the correlation functions. It is shown that using the selection parameters of interatomic interaction, mobile system with impurities that are in thermodynamic equilibrium, can be approximated to a system with frozen impurities.

Keywords: phase transitions, Potts model, Bethe lattice.

Частным случаем решетки Бете [1, 2] является одномерная цепочка подвижных магнитных и немагнитных атомов, для которой можно построить точное решение с помощью трансфер-матрицы [1]. Это решение обладает тем преимуществом, что оно позволяет сравнительно легко рассчитать корреляционную функцию.

Запишем статистическую сумму одномерной цепочки, состоящей из N узлов, в следующем виде:

Zn = 2м + Lafaf+X + га} + Л^)} , (1)

где г = ц/кТ и используется циклическое граничное условие = &ц+\.

1 Сергей Викторович Сёмкин - канд. физ.-мат. наук, доцент кафедры информационных технологий и систем; e-mail: Li15@rambler.ru.

2 Виктор Павлович Смагин - д-р. физ.-мат. наук, зав. лаб. геофизических полей; e-mail: Li15@rambler.ru.

* Работа выполнена при финансовой поддержке Министерства образования и науки РФ по заданию № 2014/292 на выполнение государственных работ в сфере научной деятельности в рамках базовой части государственного задания.

Вычислив тем или иным способом (1) можно найти химический потенциал и спонтанную намагниченность из следующих соотношений:

,, 1 д1п Zы 1 д1п Zм

рМ = ——11 и р =--т-2-. (2)

и N дке N дг

Для вычисления ZдГ воспользуемся следующим способом [1]. Рассмотрим трансфер-матрицу V:

/ 1 е(Г+ке)/2 е(г-ке)/2 \

V = \ (г+Ле)/2 оК+Ь+Г+ке &-К+Ь+Г | (3)

\ (1—Ле)/2 е-К+Ь+г еК+Ь+г-ке I

Если Ху, А-2 и Я3 и - собственные числа матрицы (3), то статистическая сумма (1) равна = Я^ + А^ + Я^. Пусть Я^ есть максимальное из собственных чисел

матрицы V. Тогда в термодинамическом пределе (Л^ —* оо) формулы (2) перейдут в

рМ = ^ „ р = (4)

Собственные числа (3) находятся из характеристического уравнения:

Я3 + аЛ2 + ЬЛ + с = 0, (5)

где

а = -(1 + 2е^к+гсЪ(ке)) Ъ = 2(е2<Ук+гЬЪ(2К) + - 1)егсЬ(Ле))

с = -4е^+2г(е^сЬ(А') - 1)$Ъ(К).

Дифференцируя (4) по Г и по Не и используя (5), можно записать выражения для нахождения М и Г через Лу и производные коэффициентов а, Ь и с:

1 да I дь -¡2да,1 дь дс

пМ = - ее ъ = - ** ^ (6)

^ ЗЯ?+2аЯ1+й' ^ (ЗЯ|+2аЯ1+Ь)А1' (6)

Решая (5) численно или по формулам Кардано и находя Я^, можно, используя

(6), найти зависимость намагниченности М от р и ]ье.

Определим для модели Изинга с подвижными примесями корреляционные

2 2

функции. Ковариацию величин О^ и оу (средние значения которых равны концентрации магнитных атомов р), рассматриваемую как функция от расстояния между этими узлами, будем называть «позиционной» корреляционной функцией:

Эта функция показывает коррелированность в расположении магнитных атомов (или атомов примеси); впрочем, даже если все д^ равны нулю, это еще не означает, что примеси распределены в решетке полностью хаотично. Ковариацию самих величин о^ и оу

д™р = <о*о)> - р2М2

будем называть «магнитно-позиционной» корреляционной функцией, поскольку она характеризует связь между расположением и магнитными моментами атомов

одновременно. Для того чтобы описать корреляцию только между величинами магнитных моментов, определим «магнитную» корреляционную функцию:

8ц = -м2 = сд™р + р2м2)/(я5 + р2) - и2.

значения этой функции есть разницы условных средних от произведений оуоу и произведений условных средних оу и оу при условии, что обе эти величины не равны нулю.

рассмотрим вычисление этих корреляционных функций для линейной цепочки. Для одномерной цепочки, состоящей из N узлов

где

У(?1.<Г1+1) = ехР {К<Г1<*1+1 + + г(Р1 + о1+1)/2 + + а1+1)/2}

элемент трансфер-матрицы.

'0 о

0 | и запишем (7) в следующем виде:

Введем матрицу 5 = [ 0 1

0 ь

{а} а}) = -^-ТтУ^БУ}-^"-^1-

(8)

Для вычисления следа в (8) рассмотрим ортогональную матрицу Р, приводящую симметричную трансфер-матрицу V к диагональному виду:

Тогда (8) запишется в виде:

ДГ1 0 0

= 7-Тг

яГу+1

где 5. Учитывая, что Zдг = Я^ + Я-У + Я^, и переходя к пределу N -* оо, получим

<^2вт/> = 4гТг

о

(9)

выражая (9) через элементы матрицы 5, получим:

{а}о}) = + §12§211Щ' + §13§31 ^

Я2

Я1

А1 /

(10)

вычисляя аналогичным образом {(Г2), можно показать, что оно равно 5ц. Поэтому позиционная корреляционная функция является суммой двух убывающих геометрических прогрессий:

д* = <<7?<т/> - р2 = ахГ + Ъх>2-1, (11)

где а = 512^21, Ь = 51з5з1, х1 = Я2/Я15 х2 = Яз/А^ Для расчета магнитно-позиционной корреляционной функции нужно использовать вместо матрицы Б матрицу

0 °\ гщ>

$1 = I 0 1 0 1' что приводит для д-у к уравнению, аналогичному (11):

\0 0 -1/

= агхС + М2У"!- (12)

Оказывается, что при 1ге = 0 (в этом случае намагниченность М также обращается в ноль - в одномерной цепочке нет возникновения спонтанной намагниченности при ненулевой температуре) один из коэффициентов а или Ь в формуле (11) обращается в ноль (какой именно, зависит от нумерации корней и Я3 характеристического уравнения), то есть позиционная корреляционная функция в этом случае имеет вид убывающей геометрической прогрессии. Что касается магнитно-позиционной корреляционной функции (12), то она также становится при Ке = 0 геометрической прогрессией, но если в (11) знаменатель прогрессии Ху то в (12) - и наоборот. Если же Не Ф 0, то все четыре коэффициента а, Ь, а^ и в выражениях (11) и (12) не равны нулю

Рис. 1. Концентрационная зависимость корреляционных функций ближайших соседей в одномерной цепочке при К = 1,1ге = 0, у = 1. По горизонтальной оси - концентрация магнитных атомов р, по вертикальной - значение корреляционной функции.

Кривая 1 - д\2 (р), кривая 2 - д7^ (р), кривая 3 - д^. (р)

На рисунках 1 и 2 показаны графики позиционной д±2 (1), магнитно-позиционной д7^ (2) и магнитной д^ (3) корреляционных функций в зависимости от концентрации магнитных атомов р; на рис. 1 - при ке = 0, а на рис. 2 - при Ке = 0,3.

Рис. 2. Концентрационная зависимость корреляционных функций ближайших соседей в одномерной цепочке при К = 1, ке = 0,3, у = 1. По горизонтальной оси -концентрация магнитных атомов р, по вертикальной - значение корреляционной функции. Кривая 1 - д^2 (р) = 0, кривая 2 - д^ (р), кривая 3 - д^ (р)

Позиционная корреляционная функция дзависит, в частности, от значения величины и - эффективного потенциала кулоновского взаимодействия.

Рис. 3. Концентрационная зависимость корреляционных функций ближайших соседей в одномерной цепочке при К = 1, /1е = 0, у = у0. По горизонтальной оси - концентрация магнитных атомов , по вертикальной - значение корреляционной функции. Кривая 1 - д^2 (р) = 0, кривая 2 - д^р (р), кривая 3 - д^ (р)

Подберем теперь величину II с таким расчетом, чтобы при заданных значениях внешнего поля /1е, температурного параметра К и концентрации магнитных атомов

"О

р корреляционная функция д^2 равнялась бы нулю. Согласно изложенному выше при Не = 0 такой выбор и приводит к обращению в ноль всех значений позиционной корреляционной функции д?., то есть распределение примесей в этом случае является не коррелированным на любых расстояниях.

О 8 -1-1-1-1-1-1-1-1-1-

8

0.7 0.6 0.5 0.4 0.3 0.2 0.1 0

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Рис. 4. Концентрационная зависимость корреляционных функций ближайших соседей в одномерной цепочке при К = 1,Ке = 0,3, у = у0. По горизонтальной оси -концентрация магнитных атомов р, по вертикальной - значение корреляционной функции. Кривая 1 - д%2 (р), кривая 2 - д™2 (р), кривая 3 - д^ (р)

Строго говоря, это не означает (хотя и не исключает этого), что распределение примесей полностью случайно, поскольку для полностью случайного распределения необходимо равенство нулю не только парных ковариаций, но и тройных, четверных и т.д. Однако равенство нулю парной корреляционной функции позволяет рассматривать распределение примесей как «приближенно случайное». Если же Не Ф 0, то обращение в ноль только д\2 уже не приводит к обращению

~ V ~

в ноль всех значений о г., но все же можно считать, что расположение примесей

75

более хаотично, чем при д^2 ^0. В дальнейшем будем называть распределение примесей с д^2 = 0 «псевдохаотичным».

Рис. 5. Концентрационная зависимость корреляционных функций для соседей, следующих за ближайшими в одномерной цепочке при К = 1,Не = 0,3, у = у0. По горизонтальной оси - концентрация магнитных атомов р, по вертикальной - значение корреляционной функции. Кривая 1 - ^3(р), кривая 2 - кривая 3 - д£з(р)

На рисунках 3 и 4 показаны графики позиционной (1), магнитно-позиционной д^ (2) и магнитной д^ (3) корреляционных функций в зависимости от концентрации магнитных атомов р, причем для каждого значения величина у = и// подбирается так, чтобы = 0. На рис. 3 /ь = 0, а на рисунке 4 Не = 0,3. На рисунке 5 показаны графики функций д^ д^ и д^ (кривые 1, 2 и 3 соответственно) при ке = 0,3 и для такого у, при котором д^2 = 0. Видно, что хотя д^ отлично от нуля, но ее величина значительно меньше других корреляционных функций.

1. Бэкстер, Р. Точно решаемые модели в статистической механике / Р Бэкстер. - М.: Мир, 1985. - 486 с.

2. Сёмкин, С.В. Модель Поттса на решетке Бете с немагнитными примесями / С.В. Сёмкин, В.П. Смагин // Журнал экспериментальной и теоретической физики. -2015. - Т. 148, Вып. 4 (10). - С. 729-733.

© Сёмкин, С.В., 2016 © Смагин, В.П., 2016

Для цитирования: Сёмкин, С.В. Одномерная модель Изинга с подвижными примесями / С.В. Сёмкин, В.П. Смагин // Территория новых возможностей. Вестник Владивостокского государственного университета экономики и сервиса. - 2016. - № 2. - С. 114-120.

For citation: Semkin, S.V. The one-dimensional Ising model with mobile impurities / S.V. Semkin, V.P. Smagin // The Territory Of New Opportunities. The Herald of Vladivostok State University of Economics and Service. - 2016. - № 2. - P. 114-120.

Дата поступления: 20.02.2016.