Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

nonequilibrium_statistical_operator

.pdf
Скачиваний:
13
Добавлен:
10.02.2015
Размер:
697.54 Кб
Скачать

Для решения полученной системы разностных уравнений, удобно ввести в

рассмотрение оператор трансляций pˆ , следующим образом действующий на функцию ϕm ( j, A)

pˆ nϕm ( j, A) = ϕm ( j + n, A).

(П.46)

С помощью этого оператора система уравнений может быть записана в виде

ϕ1

 

 

 

η0 +

1

 

ϕ0 (A) 1 (pˆ

1

+ pˆ )ϕ1 ( j, A) + 2η2ϕ2 ( j, A),

(П.47)

( j, A) =

2

η2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pˆ 1ϕ

( j, A) = 1

η ϕ

(A) +

 

1 η

pˆ 1 + p

pˆ

1 ϕ ( j, A) ϕ

 

( j, A)

, (П.48)

2

 

 

 

4

 

1 0

 

 

 

 

2

2

 

 

 

 

0

 

1

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

( j, A) = 1 η ϕ

(A) +

 

1 η

pˆ 1

+ pˆ

pˆ

ϕ

( j, A) ϕ

( j, A). (П.49)

 

2

 

 

4

1 0

 

 

 

 

2

 

2

 

 

 

 

0

 

1

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исключим из этих уравнений функции ϕ2 ( j, A )

и ϕ2 ( j, A ) . Для этого из

уравнения (П.49) вычитаем (П.48), получаем

 

 

 

 

 

 

 

 

 

(1pˆ 1 )ϕ

2

( j, A) = η ( pˆ

pˆ 1 )ϕ ( j, A) ( pˆ

1)ϕ

( j, A)

= η ( pˆ

2 1)ϕ

( j, A) .

 

 

 

 

 

0

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(П.50)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее из уравнения (П.49) выражаем функцию ϕ2 ( j, A ) и подставляем

результат в уравнение (П.47), находим

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ1 ( j, A) = η0ϕ0 (A) +

(η12 −η22 )(pˆ + pˆ 1 )2η0η2 pˆ

ϕ1 ( j, A) +

η

 

( j, A)

. (П.51)

 

 

 

 

 

 

 

 

η

 

 

 

 

 

2 η ϕ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

Домножаем полученное уравнение на оператор p€-1 слева

 

 

 

 

 

 

( pˆ 1)ϕ ( j, A) ( pˆ

1)ϕ (A) +( pˆ 1)

12 −η22 )(pˆ + pˆ

1 )2η0η2 pˆ

ϕ ( j, A) +

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

η1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

81

+2 η2 ( pˆ 1)ϕ2 ( j, A)

η1

и исключаем отсюда функцию ϕ2 ( j, A) с помощью уравнения (П.50). В результате получаем следующее уравнение, содержащее только одну неизвестную функцию ϕ1( j, A)

( pˆ 1)ϕ1( j, A) 11 (η12 −η22 )( pˆ2 pˆ +1pˆ1) 2η0η2 ( pˆ2 pˆ) +2η0η2 ( pˆ2 1) ϕ1( j, A),

которое удобно представить в виде

pˆ 2

pˆ

1 + g(1pˆ )

ϕ ( j, A) = 0

,

(П.52)

 

 

 

1

 

 

где

g =1+

η12 2η0η2

= 2ch2 (β ω

/ 2)(1+cth(β ω ))1. (П.53)

 

η2

−η2

0

e

 

1

2

 

 

Уравнение (П.52) представляет собой разностное уравнение

ϕ1( j + 2, A ) −ϕ1

( j 1, A ) = g(ϕ1

( j +1, A ) −ϕ1

( j, A )). (П.54)

Решение полученного одномерного разностного уравнения ищем методом Эйлера, согласно которому решение ищется в виде ϕ1( j, A) = λj . Подстановка этого выражения в уравнение (П.54) приводит к следующему характеристическому уравнению на λ

λ3 g(λ2 −λ) 1 = 0 .

(П.55)

Корни уравнения (П.55) имеют вид

λ1 =1, λ2,3

≡ λ± =

1

(g 1±

(g 1)

2

4 )=

1±

z02 +e−β ωe ch2 (β ω0

/ 2)

, (П.56)

2

 

 

 

 

 

 

 

 

 

 

1

2

+e

−β ω

2

(β ω0

/ 2)

 

 

 

 

 

 

z0

e ch

 

 

82

и обладают следующими

свойствами

λ+ 1, λ1, λ+λ=1.

В

результате

общее решение разностного уравнения (П.54) запишется в форме

 

z

+

j

+

j

(П.57)

ϕ1( j, A ) =< S j

A >= C1( A ) +C1

( A )λ+ +C1

( A )λ,

где C1 (A), C1± (A) - произвольные постоянные, значения которых зависят от вида оператора Aи должны определяться из дополнительных соотношений, играющих роль граничных условий. Выражение для остальных неизвестных функций имеет аналогичный вид

+

j +

j

(П.58)

ϕm ( j, A ) = Cm ( A ) +Cm ( A )λ+ +Cm ( A )λ,

где Cm (A), Cm± (A) являются некоторыми линейными комбинациями констант

C1 (A), C1± (A) и ϕ0 ( A).

В бесконечной цепочке Изинга, величины корреляционных функций

 

при

j → ∞

должны быть

ограниченными, в

связи

с чем

ϕm ( j, A )

 

коэффициент

+

равен нулю.

Способ

определения

постоянных

Cm

(A )

проиллюстрируем на

примере

нахождения

констант

Cm (A ), Cm (A )

для двухчастичного коррелятора

z

z

z

).

В силу

C1 (A ), C1 (A )

< S j

Si >= ϕ1 ( j, Si

трансляционной симметрии модели функцию ϕ ( j, S z )

можно представить в

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

i

 

 

 

 

 

 

 

 

 

виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ ( j, S z ) = ϕ (

 

j i

 

, S z ) =< S z

 

S z >= C (S z ) +C

(S z )λ| ji| .

(П.59)

 

 

 

1

i

1

 

 

 

0

 

 

ji

0

1

0

 

 

1

 

0

 

 

Полагая далее в уравнениях (П.48) и (П.49)

 

 

 

 

 

 

 

 

 

 

 

A =1, получим два соотношения,

связывающих три неизвестные функции

 

C (S z ), C(S z ) и< S z

>

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

1

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

12η1

 

 

 

1

 

 

 

η1

 

 

 

 

 

 

C1 (S0z ) +C1(S0z )λ2

< S0z >= −

+

 

 

 

 

,

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2η2

 

 

 

 

 

 

2η2

 

 

 

 

C (S z ) +C(S z )

λ

 

−λ2

 

 

η +η

< S z

 

 

 

 

η

 

 

 

 

0

2

>=

 

 

 

 

1

.

(П.60)

 

1−η

 

4(1−η )

 

1

0

1 0

 

1

−η

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

Третье соотношение получим из тождества

83

 

 

< S z S z >= ϕ (0, S z ) = C (S z ) +C(S z ) = 1 .

 

 

(П.61)

 

 

0

 

0

 

1

0

1

0

1

0

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Система

уравнений

(П.60), (П.61) позволяет определить намагниченность

< S z > и константы C (S z ), C(S z ) двухчастичного коррелятора (П.59):

 

0

 

 

1

0

 

1

0

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

sh(β ω0

/ 2)

 

 

 

 

 

 

 

 

< S0

>=

 

 

 

 

 

 

 

 

,

 

 

(П.62)

 

 

 

sh2 (β ω

/ 2) +exp(−β ω )

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

e

 

 

 

 

C (S z ) =< S z >2

,

 

 

C

(S z ) =

1

 

exp(−β ωe )

 

.

(П.63)

 

 

4 sh

 

 

1

0

0

 

 

 

1

0

2 (β ω

 

/ 2) +exp(−β ω )

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

e

 

 

В результате выражение для двухчастичного коррелятора примет вид

 

 

< S0 Sm >=< S0

>

2

+

− < S0

>

2

λ.

(П.64)

 

 

z z

z

 

1

z

 

 

m

 

 

 

 

 

 

 

4

 

 

 

 

 

 

Найдем

также

выражения

для

 

 

трехчастичных

корреляторов

вида

< S0z Smz 1Smz

+1 > и

< S0z Smz Smz +1 >.

Для

этого

полагаем

в

уравнениях (П.47) и

(П.51)

A = S0z , в результате получаем

 

 

 

 

z

z

 

η0

+

1

 

 

z

> +η1 (p

1

z

z

> +2η2

< S0

Sm >=

2

η2

< S0

 

+ p)< S0

Sm

 

 

 

 

 

 

 

 

 

 

 

 

 

< S z S z

>= η < S z

> +

(η12

−η22 )(p + p1 )2η0η2 p

< S z S z

 

 

 

 

 

 

 

0

m

 

0

 

0

 

 

 

 

η1

 

 

0 m

 

 

 

 

 

 

 

 

 

 

 

 

 

< S z S z S z + >, (П.65)

0 m 1 m 1

> +2

η2

< S0z Smz Smz +1 >.

 

η1

 

 

 

(П.66)

После некоторых простых преобразований отсюда непосредственно следуют следующие выражения для трехчастичных корреляционных функций

< S z S z

S z

>=< S z >< S z S z > +

 

1

− < S z >2

sh(β ω )(1

+cth(β ω

/ 2))λm ,

0 m1

m+1

0

0 2

 

4

0

 

0

e

 

 

 

 

 

 

 

 

 

 

(П.67)

84

< S z S z S z

>=< S z >< S z S z > +

 

1

− < S z >2

th(β ω

/ 2)(1

)λm . (П.68)

0 m m+1

0

0 1

 

4

0

 

0

 

 

 

 

 

 

 

 

 

 

 

Аналогичным образом из системы (П.47-П.49) можно найти корреляционные функции как с большим числом спинов, так и с различной пространственной структурой. Анализ, проведенный М.П. Желифоновым [23], показал, что в

общем

случае

n-частичная

корреляционная функция

< S jz S jz

...S jz > может

 

 

 

 

 

 

 

 

 

 

 

1 2

n

быть записана в виде

 

 

 

 

 

 

 

 

 

 

 

 

< S j1 S j2

...S jn >=< S j

>

 

m

 

λi

,

(П.69)

 

 

n

+∑∑Cn

 

 

z

z

z

z

 

 

α

qα

 

 

 

 

 

 

 

 

 

α=1 i

 

 

 

 

где величины

qα представляют собой всевозможные линейные комбинации

 

 

i

 

 

 

 

 

 

 

 

 

 

вида

| j1 + j2 +... + jα jα+1 ... j2α |

координат

частиц,

входящих в

коррелятор; Cnα - произвольные постоянные, m=n/2 при n четном и m=(n+1)/2 при n нечетном. В качестве граничных значений для определения Cnα более сложных чем ϕm ( j, S0z ) , корреляторов можно использовать уже найденные значения Cmα .

85

Литература

1.Зубарев Д. Н., Статистический оператор для неравновесных систем, ДАН

СССР 140, 92 (1965).

2.Зубарев Д. Н., Локально равновесный ансамбль Гиббса и его связь с теорией флуктуаций и явлениями переноса, ДАН СССР 162, 532 (1965).

3.Зубарев Д. Н., Процессы переноса в системах с внутренними степенями свободы, ДАН СССР 162, 794 (1965).

4Зубарев Д. Н., Статистический оператор для нестационарных процессов,

ДАН СССР 164, 537 (1965).

5.Зубарев Д. Н., Калашников В. П., Экстремальные свойства неравновесного статистического оператора, ТМФ 1, 137 (1969).

6.Зубарев Д. Н., Граничные условия для статистических операторов в теории неравновесных процессов и квазисредние, ТМФ 3, 276 (1970).

7.Зубарев Д. Н., Калашников В. П., Построение статистических операторов для неравновесных процессов, ТМФ 3, 126 (1970).

8.Зубарев Д. Н., Калашников В. П., Теория возмущений и интегральные уравнения для неравновесных статистических операторов, ТМФ 5, 406 (1970).

9.Зубарев Д. Н., Калашников В. П., Эквивалентность некоторых методов в статистической механике необратимых процессов, ТМФ 7, 372 (1971).

10.Зубарев Д. Н., Калашников В. П., Вывод необратимого во времени обобщенного основного кинетического уравнения. Дубна: ОИЯИ, Препринт

P4-P5658, 1971, с. 3-31.

11.Kuzemsky A. L., Theory of transport processes and the method of the nonequilibrium statistical operator, International Journal of Modern Physics B 21, 2821 (2007)

12.Markiv B. B., Tokarchuk R. M., Kostrobij P.P., Tokarchuk M.B., Nonequilibrium statistical operator method in Renyi statistics, Physica A 390, 785 (2011).

86

13.Ryazanov V. V., Lifetime distributions in the methods of non-equilibrium statistical operator and superstatistics, The European Physical Journal B 72, 629 (2009).

14.Adams J., Reinholz H., Redmer R. and French M., Linear response treatment of the Hall effect within the Zubarev formalism, J. Phys. A 39, 4723 (2006).

15.Huanga X.-G., Sedrakiana A., Rischkea D.H., Kubo formulas for relativistic fluids in strong magnetic fields, Annals of Physics 326, 3075 (2011).

16.Kochelaev B.I. Spin temperature and non-equilibrium phonons, NMR and more (a volume in honor of Anatole Abragam). – Paris: Les Editions de Physique les Ulis, 1994. – P. 279-292.

17.Альтшулер С.А., Валишев Р.М., Кочелаев Б.И., Хасанов А.Х., Исследование фононной лавины по Мандельштам-бриллюэновскому рассеянию света в условиях насыщения магнитного резонанса, ЖЭ'ТФ 62, 239 (1972).

18.Фазлеев Н.Г. ФНТ 5, 380 (1979).

19.Фазлеев Н.Г. ФНТ 6, 1422 (1980).

20.Nigmatullin R. and Tayurskii D. A straightforward memory function calculation with the non-equilibrium statistical operator method, Physica A 175,

275(1991).

21.Боголюбов Н.Н. Избранные труды. В 3 т. – Киев: Наук. Думка, 1970.

22.Кессель А.Р., Берим Г.О. Магнитный резонанс изинговских магнетиков. –

М.: Наука, 1982. – 147 с.

23. Желифонов М.П. Высшие корреляционные функции изинговского ферромагнетика. Случай линейной модели для спина S=1/2, ТМФ 8, 401 (1971).

24.Bariev R.Z., Zhelifonov M.P., The new formulation of Ising problem, Physica Letters 50A, 105 (1974).

25.Bariev R.Z., Higher-order correlation functions of the planar Ising model, Physica 83A, 388 (1976).

26.Bariev R.Z., On the rotational symmetry of the spin-spin correlation function of the two-dimensional Ising model., Physica Letters 55A, 456 (1976).

27.Зубарев Д.Н. Неравновесная статистическая термодинамика. – М.: Изд. Наука, Гл. ред. физ.-мат. лит., 1971. – 415 с.

28.Зубарев Д.Н., Морозов В.Г., Репке Г. Статистическая механика неравновесных процессов. Т. 1.. – М.: Физ.-мат. лит., 2002.- 432 с.

29.Зубарев Д.Н., Морозов В.Г., Репке Г. Статистическая механика неравновесных процессов. Т. 2.. – М.: Физ.-мат. лит., 2002.- 296 с.

30.Ляпилин И.И. Введение в теорию кинетических уравнений. Учебное пособие. - Екатеринбург: УГТУ-УПИ, 2004. - 332 с.

87

31.Ляпилин И.И., Калашников В.П. Неравновесный статистический оператор и его приложения к кинетике парамагнитных явлений в проводящих кристаллах. – Екатеринбург: УрО РАН, 2008.- 366 с.

32.Биккин Х.М., Ляпилин И.И. Неравновесная термодинамика и физическая кинетика. – Екатеринбург: УрО РАН, 2009. – 500 с.

33.Пригожин И. От существующего к возникающему: время и сложность в физических науках. - М.: Наука, 1985.

34.Гуров К. П. Феноменологическая термодинамика необратимых процессов. - М.: Наука, 1978.

35.Кubо R., Statistical-Mechanical Theory of Irreversible Processes: I. General Theory and Simple Applications to Magnetic and Conduction Problems, J. Phys. Soc. Japan. 12, 570 (1957).

36.Пелетминский С.В., Яценко А.А., Метод производящего функционала и вириальные разложения в неравновесной статистической механике, ТМФ,

3, 287 (1970).

37.Покровский Л.А., Получение обобщенных кинетических уравнений с помощью неравновесного статистического оператора, ДАН СССР 183, 806 (1968).

38.Glauber R.J. Time Dependent Statistics of the Ising Model, J. Math. Phys. 5,

294(1963).

39.Matsudaira N. Some dynamical properties of the Ising ferromagnet, Canad. J. Phys. 45, 2091 (1967).

40.Matsudaira N. Some Dynamical Properties of the Ising Ferromagnet. II. Cubic Lattices, J. Phys. Soc. Jap. 23, 232 (1967).

41.Tanaka M., Takahashi K., Kinetic Ising Model with the Bilinear and Biquadratic Interactions, J. Phys. Soc. Jap. 43, 1832 (1977).

42.Буишвили Л.Л., Звиададзе Н.Д. Препринт ин-та теор. физики, №82, Р.

Киев, 1970.

43.Провоторов Б.Н. О магнитном резонансном насыщении в кристаллах,

ЖЭТФ 41, 1582 (1961).

44.Fisher M. E., Perpendicular susceptibility of the Ising model, J. Math. Phys. 4,

124(1963).

45.Allan G. A., Betts D. D., The frequency-dependent initial perpendicular susceptibility of the Ising model, Can. J. Phys. 46, 15 (1968).

46.Halperin B. L, Hohenberg P. C, Ma S., Renormalization-group methods for critical dynamics: I. Recursion relations and effects of energy conservation, Phys. Rev. B 10, 139 (1974).

47.Houtappel R. M. F., Order-disorder in hexagonal lattices, Physica 16, 425 (1950).

88

48.Callen H. B. A note on Green functions and the Ising model, Phys. Lett. 4, 161 (1963).

49.Suzuki M. Generalized exact formula for the correlations of the Ising model and the other classical systems, Phys. Lett. 19, 26 (1965).

50.Гантмахер Ф.Р. Теория матриц. – М.: Наука, 1966.- 576 с.

89

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]