Вестник Северо-Кавказского государственного технического университета. 2011. № 1 (26)
10
References
1. Levshin, L. V. Luminescence and its measure. Molecular luminescence / L. V. Levshin, A. M. Saletskiy. – M. :
MGU, 1989. – 272 p.
2. Parker, C. Photoluminescence of solutions / C. Parker. – M. : Mir, 1972. – 511 p.
3. Bagnich, V. V. Kinetics of triplet excitation annihilation in a organic glasses / V. V. Bagnich, A. V. Konash // Optic
and spectroscopy. 2002. Vol. 92. № 4. P. 556 – 563.
4. Nonradiative transfer of electronic excitation energy / V. L. Ermolaev, E. N. Bodunov, E. B. Sveshnikova,
T. I. Shakhverdov. – Leningrad : Nauka, 1977. – 311 p.
5. Sapunov, V. V. On triplet-triplet annihilation kinetics in liquid solutions / V. V. Sapunov // Optic and spectroscopy.
1995. Vol. 79. № 4. P. 606 – 612.
6. Bryukhanov, V. V. Temperature effect on triplet-triplet annihilation rate of 1, 2-benzantracene in polymeric matrix /
V. V. Bryukhanov, I. G. Samusev, S. G. Karstina // Journal of Applied Spectroscopy. 2004. Vol. 71. № 1. P. 49 – 53.
7. Kulikova, O. I. Kinetics organic molecules delayed fluorescence in n-paraffin matrices and its mathematical model /
O. I. Kulikova, T. V. Zheludkova, V. V. Solodunov // News of HEI. Phusics. 2007. № 12. P. 3 – 8.
8. Harmon, L. A. Triplet excitation transport kinetics in vapor-deposited naphthalene / L. A. Harmon, R. Kopelman //
J. Phys. Chem. 1990. Vol. 94. № 9. P. 3454 – 3461.
УДК 537.612:519.6
ДВУМЕРНАЯ МОДЕЛЬ ИЗИНГА
С ТРЕМЯ ВОЗМОЖНЫМИ СОСТОЯНИЯМИ УЗЛОВ
Ратнер И. М.
НОУ «Ставропольский институт управления»
Разработан алгоритм и программа на языке Паскаль для вычисления элементов и собственных значений
трансфер-матрицы решетки Изинга с тремя возможными состояниями узлов. Приводятся результаты расчетов
собственных значений решетки без учета и с учетом внешнего магнитного поля. Результаты могут служить для
проверки моделей диагонализации трансфер-матрицы решетки подобного типа.
The algorithm and program on Pascal language to calculate elements and proper values of Easing transfer-matrix
lattice with three possible cross-points conditions have been developed. The results of lattice proper values calculations are
given without taking and with taking into consideration the external magnetic field. The results can be used to check models of
transfer-matrix lattice diagonalization of such type.
Ключевые слова: алгоритм на языке Паскаля, двумерная модель Изинга, трансфер-матрицы решетки.
Key words: algorithm on Pascal language, Easing two-dimensional model, transfer-matrix lattices.
Двумерная модель Изинга без магнитного поля была точно решена в работах Онсагера –
Кауфман [1, 2]. Предполагалось, что каждый узел решетки имеет два возможных состояния,
нумеруемых как +1 и –1. Это соответствует утверждению, что источником ферромагнетизма
является взаимодействие спинов атомов.
Представляет значительный интерес исследование возможности фазового перехода
в системе, узлы которой обладают тремя возможными состояниями. Это
может быть либо система
спаренных электронов, либо система взаимодействующих орбитальных моментов атомов.
В настоящей работе рассматривается двумерная решетка, узлы которой имеют три
возможных состояния
s, нумеруемых как s={+1, 0, –1}. Параметры взаимодействия по горизонтали
J
1
и по вертикали J
2
считаем различными. Для вычисления статистической суммы модели найдем
элементы матрицы взаимодействия (трансфер-матрицы), которую представим в симметричной
форме [3]:
)exp()
2
(exp)
2
exp(,...,,...,
!
2
!
1
!
1
1
1
!
!!
11 kk
kkkk
n
k
kk
nn
ss
ssssss
ssPss
θθβ
++
=
++
=
∏
, (1)
где
T
B
=
β
,
T
J
1
1
=
θ
,
2
2
J
T
θ
= , 2)
и
T абсолютная температура в энергетических единицах: n – число атомов в ряду, B – величина
магнитного поля, последний атом в ряду взаимодействует с первым (циклические граничные
условия Борна – Кармана [4]):
s
n+1
≡
s
1
. (3)
Форма записи матричных элементов в (1) такая же, как и для решетки с с двумя состояниями
узлов, отличие в том, что величины
s
k
могут принимать три значения, так что члены статистической
суммы образуют матрицу третьего порядка: