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

книги / Релаксационные явления в полимерах

..pdf
Скачиваний:
3
Добавлен:
13.11.2023
Размер:
23.13 Mб
Скачать

С учетом вида гамильтониана (VII. 2) получим»

 

 

Ц (D.) =

(р),

£>, (t, г,

р)]

 

L1(Dl) =

J [V (I Г -

Г

I). D f (г, / , р,

р'. Я ,)] dr dp

(VII. 36)

L2(^l) = J

 

 

 

 

 

 

(I

r !),

£>2* (r<r/-

P- РГ>°i)] * Г dPr и т- д>

 

где H(p) — pz/2m.

Далее задача решается методом последовательных приближе­ ний, и в первом приближении получаем кинетическое уравнение для одночастичной функции распределения в виде:

dD' % Г' Р) + £ •— 1

Р) - I J (I г -

г' |), D, (/. *<2>, р<2>) X

 

 

X £>, (/,

Я(2)\ р<2)')1 dr dp'

(VII. 37)

Из этого уравнения можно вывести уравнение Больцмана, что было сделано Боголюбовым [37] и другими авторами для различ­ ных конкретных задач [60—64]. Рассмотрен также вывод кинетиче­ ских уравнений из уравнения Лиувилля для систем частиц с внут­ ренними степенями свободы [65—68]. Из уравнения (VII. 37), или уравнения Больцмана с учетом законов сохранения, выполняю­ щихся при взаимодействиях в системё, можно вывести кинетиче­ ское уравнение вида:

~ д Г = 2 (“’vjiO'v - æ'nvO'n) (VII.38)

где Wpv — отнесенная к единице времени вероятность перехода из состояния р в состояние у, а — вероятность заполнения состоя­ ния р.

Точное решение нелинейного интегро-дифференциального кине­

тического уравнения (VII. 37) в

общем

случае невозможно, хотя

для отдельных частных случаев

[69, 70]

имеются доказательства

существования и единственности решения кинетических уравнений подобного типа.

Практически кинетические уравнения такого типа решаются для малых отклонений от равновесного состояния с использова­ нием линеаризации уравнения. Однако в любом случае анализ полученного уравнения дает многое для понимания особенностей происходящих в системе процессов и позволяет провести хотя бы качественную оценку исследуемых макроскопических характери­ стик.

Метод кинетических уравнений еще не нашел достаточно ши­ рокого применения в физике полимеров. Однако можно .назвать некоторые работы [27—36], в которых методом кинетических урав­ нений решаются различные задачи физики полимеров с исполь­ зованием диффузионных многочастичных уравнений или уравне­ ний типа (VII. 38) применительно к модели полимерной цепи из

жестких или полужестких элементов с использованием модели Изинга. Применимость рассмотренного выше метода кинетиче­ ских уравнений длямикроописания физико-химических свойств, полимеров определяется в основном видом потенциала взаимодей­ ствия частиц системы и степенью ее начального возмущения. Трудности применения метода кинетических уравнений встре­ чаются в том случае, когда характер взаимодействия частиц и их конфигурация таковы, что имеется существенная многочастичная корреляция движения частиц. Однако, если время релаксации та­ ких коллективных, коррелированных движений многих частиц меньше характерного временного масштаба заметной эволюции функции распределения, метод кинетических уравнений применим. Все это позволяет говорить о том, что рассмотренные выше схемы методов кинетических уравнений для классических систем могут быть применены для широкого исследования физических свойств растворов полимеров и отдельных макромолекул,' а также при ре­ шении отдельных задач других, более сложных систем. Эффектив­ ное применение этого метода для широкого исследования свойств полимеров различных классов требует дополнительной разработ­ ки математического аппарата рассматриваемых методов.

Метод кинетических уравнений можно изложить, используя формализм функций распределения в фазовом пространстве, как это сделано выше, или с помощью методов вторичного квантова­ ния в фазовом пространстве [71].

Не останавливаясь на методе кинетических уравнений для квантовых систем, отметим лишь, что применение его в физике по­ лимеров также возможно для решения целого ряда проблем, на­ пример при анализе процессов спин-решеточной релаксации в полимерных системах, при построении микротеории явления термофлюктуационного механизма разрушения и прочности полиме­ ров [6 , 23], для микроописания элементарных актов механизма электропроводности [78, 79], в особенности ее электронной состав­ ляющей, например, в полимерных полупроводниках.

ФЛЮКТУАЦИИ И ФИЗИЧЕСКИЕ СВОЙСТВА ПОЛИМЕРОВ

Представления теории флюктуаций и ее методы получили доста­ точно широкое применение при рассмотрении физических явлений и процессов в электрических цепях, радиотехнике и электронике, теории информации, теории распространения и рассеяния электро­ магнитных и акустических волн в,средах, при анализе условий устойчивости [81—85]. В физике полимеров представления о флюк­ туациях и корреляции случайных величин привлекались в ограни­ ченном числе случаев, например в связи с обсуждением структуры и свойств полимеров [6 , 23]. В ряде работ [86—91] исследованы взаимосвязь и влияние флюктуационных характеристик на меха­ нические и диэлектрические свойства материалов, разрабатывается флюктуационная теория сеток, успешно применяются представле»

ния о корреляциях при разработке теории механических свойств анизотропных армированных полимерных систем. Открывающиеся возможности более полного использования представлений и мето­ дов теории флюктуаций в физике полимеров заслуживают внима­ ния и вызывают необходимость изучения этого круга вопросов фи­ зической статистики и термодинамики.

Переходя к рассмотрению основных корреляционных и спект­ ральных характеристик флюктуаций, прежде всего отметим, что под корреляциями флюктуаций в общем случае следует понимать определенный вид статистической взаимосвязи во времени и про­ странстве флюктуирующих физических величинAi{t,г), Az(t, г), ...

. . ., An(t, г) системы и их флюктуаций x\{t,r) = Ax(t, г)

A\(t, г), . . . . xn(t, r ) = An(t, г) An{t, г). Корреляцию флюктуа­ ций описывают с помощью матрицы корреляций Rx.Xk(t{, t2 г,, г 2),

представляющей собой матрицу центральных моментов второго по­

рядка [97J:

_______________

 

^Х(хк (*г

*2’ Гг гг)e xi (*г rl) хк (*2’ гг)

39)

Черта сверху означает усреднение по всей совокупности зна­

чений случайных функций Ai(t, г) и Ak{t, г). В конкретных зада­

чах такое усреднение

обычно заменяют временным

и простран­

ственным усреднением. Диагональными элементами матрицы RХХк (г = k) будут функции Rx.Xk(tv rv Гг) автокорреляции или

просто корреляции флюктуаций xt(t, г) в различных точках про­ странства в различные моменты времени. Недиагональные эле­ менты Rxtxk (IФ k) определяют функции взаимных корреляций флюктуаций x{(ti, Г|) и xk{tb г2). Матрица RX[xk симметрична по

индексам xt и х к . Флюктуации некоррелированы в точках (/|, Г|) и (^2>гг)« если Rxtxk= 0 , й вообще некоррелированы, если Rxlxk= 0 .

Для стационарных флюктуаций в системах с неизменными во времени внешними условиями, например в состояниях термо­ динамического равновесия систем, средние значения физических

величин At не зависят от времени, а матрица корреляций зави­ сит лишь от разности времен т = / , —*i2, т. е. Rx.Xk— RX[xk(r» rv

Точно так же в пространственно однородных средах Rxtxk зависит

лишь от разности р = r t — г2, т. е. RX[xk= R x{xk(U

h ,

р).

Если

к тому же пространство еще и изотропно, Rxixk зависит от

вре­

мени

и модуля вектора р. Вместо общей матрицы

корреляций

J?

(/„ t2 г (, г2) в зависимости от условий задачи рассматривают

более простые матрицы RXixk(t\» *2) и Rxtxk(rv гг)*

 

 

 

В случае временной автокорреляции стационарной

флюктуа­

ции

величины Ai функция автокорреляции имеет

вид

 

является четной функцией т и удовлетворяет условиям

RXl{0 ) > 0

и | Д * № [< /4 (0 ). Аналогичными свойствами обладает

функция

пространственной автокорреляции в пространственно

однородной

и изотропной системе.

общем виде матрицы

корреляций

Для нахождения в

*2> rv г2) или

*а) и Rx{xk(ri> гг) необходимо мик­

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

В случае стационарных флюктуаций в термодинамически рав­ новесной системе для матрицы временных корреляций флюктуи­ рующих величин имеет место формула Хинчина.— Винера [94]:

 

+°°

(VII. 40)

V f c (T)=

J

 

—ОО

 

Применим обратное преобразование

Фурье к уравнению

(VII. 40) :

+со

 

 

 

sw > - - 5 r

!

(VII. 41)

 

— со

 

Матрица временных корреляций Rx(xk(т) и матрица спектраль­ ной плотности 5x<Xft((o) связаны между собой преобразованиями

Фурье. Матрица Sx.Xk(со) является эрмитовой

(со) =

четной по параметру ш [S*.*fe(©) = Sx.X/l(— ©)],

каждый ее элемент

определяет спектральную интенсивность флюктуаций величин Ai(t) и Ak(t). Вследствие Фурье-сопряженности Rxixk(т) и Sx.Xk(et)

выполняется соотношение Д© -Дт~1, где Дт — интервал времен­ ной корреляции флюктуаций, а Д© — ширина спектра этих флюк­ туаций. Формулы (VII. 40) и (VII. 41) показывают, что корреля­ ционная матрица флюктуаций полностью определяет их спектр, и наоборот. Так, для физической величины А, флюктуации кото­ рой затухают во времени по экспоненциальному закону, корре­ ляционная функция будет иметь вид Rxx (т) = а ехр (— Ьх). Тогда матрица спектральной плотности такой флюктуации в соответ­ ствии с уравнением (VII. 41) будет иметь вид:

+ оо

 

 

s xx (©) = ■— J ехр (-Ьт + i<ùт)

yt

(VII.42)

— ОО

 

 

где а и b — некоторые постоянные.

Для определения вида матрицы спектральной плотности флюк­ туаций в термодинамически равновесных системах можно исполь­

зовать флюктуационно-диссипационную теорему [92, 93]

s *i*k <«> =

К * (ю) “

в» <®>1cth Ш

(VI1,43)

где <Xih{о») — матрица обобщенного

адмитанса системы,

опреде­

ляемая уравнениями

 

 

 

Ai W f =

2k atk (®) fokехР (-fa>T)

(VII. 44)

Здесь Аг (t) *— усредненные величины физических характери­ стик системы (усредненные отклики системы), = /олехр(—mt) — вызывающие эти отклики обобщенные сторонние силы, изменяю­ щиеся во времени по гармоническому закону. Уравнение (VII. 44) в конкретных задачах обычно получают путем линеаризации соот­ ношений, вытекающих из уравнения

Ù------ 2 кM l (VII. 45)

где Ü— усредненное значение внутренней энергии системы, через

fk обозначены сторонние

обобщенные

силы. Если kT h®

и

флюктуации классические, формула

(VII. 43)

принимает вид:

 

V

»

<•> -

■ & “ •<<•>-

 

и

<vn-

 

Это выражение

с

учетом (VII.

44)

можно

преобразовать

к

виду:

 

 

 

 

 

 

 

 

<•> - TF

<•> -

 

H

' th w

(vn-47)

Здесь a^ 1— матрица, обратная матрице ajh адмитанса системы, Sfifk() — спектральная матрица сил /{, если рассматривать спон­ танные флюктуации величины Ai как результат действия на сис­

тему

некоторых фиктивных

сторонних

флюктуационных

сил f,-.

В

случае флюктуации

х

одной

величины А по

формуле

(VII. 47) имеем:

 

 

 

 

 

Sxx (©) = J L

cth {fi(ù/2kT) Ima (о)

(VII. 48)

Отнесенная к единице времени диссипация энергии <3 в систе­ ме при действии внешней силы f0exp(—Ш) определяется форму­ лой:

Q = 4 ‘ to|/o|2/,“a(<û)

(VII. 49)

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

В случае пространственно однородных флюктуаций простран­

ственная

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

матрица

г2)

имеет

вид

Д (р)

и для нее также

справедлива

обобщенная

формула

Ви-

*ixk

 

 

 

 

 

 

нера — Хинчина:

+00

 

 

 

 

 

 

 

 

 

 

 

(Р) =

|

(ft) exp (ftp) Л

(VII. 50)

 

 

00

 

 

 

 

где p =

rj — г2,

(ft) — спектральная матрица.

 

В случае изотропности флюктуаций это уравнение прини­

мает вид;

 

оо

 

 

 

4п

 

 

 

Г

 

 

 

Rx,xk (P) =

J каХ;Хк (k) sin *Р dk

(VII. 51)

Физический смысл и свойства матрицы ax.Xk(k) и матрицы

ах.Хк{к)

аналогичны таковым

для матрицы SXixk(©). Некоторые

важные

свойства

матриц вида Rx.Xk(p) и их использование в кон­

кретных задачах рассмотрены в работе [81]. Для спектральной матрицы oXiXft(k) пространственной корреляционной функции

Rxtxk{р) в термодинамически равновесных системах также выпол­

няется флкжтуационно-диссипационная теорема. Рассмотрим несколько примеров.

Пусть система, например образец полимерного материала, при заданной температуре Т помещена в поле изменяющихся во вре­ мени упругих деформаций и [89, 90]. Внутренним параметром т) в случае полимерной системы может быть число молекул опреде­ ленной конформации, число сегментов молекулярной цепи с опре­ деленной ориентацией, величина проекции сегмента, степень ближ­ него или дальнего порядков и т. п. При малых деформациях уравнение отнесенной к единице объема свободной энергии такой системы имеет вид:

2 F = 2f» + a u 2 + 2v«! + §|* (VII.52)

где | — отклонение внутреннего параметра т) от его. равновесного значения т|о в отсутствие деформаций, а коэффициенты а, р, у яв­ ляются равновесными значениями вторых производных свободной энергии по соответствующим переменным.

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

Щ~ Р (°Р “ Y2)"”1

р2 _ kT .

Ê—Т г Р

ÔF

где

д%

 

и - M L

6 ~ V

(VII. 53)

К= - у -

рdF Ри~ да

и коэффициенты корреляций

 

(VII. 54)

Поскольку термодинамические силы Fu и

в принятом при­

ближении являются линейными функциями деформации и и внут­

реннего параметра £, коэффициенты корреляции

Ru5 и RFU сов­

падают, т. е. Rui = RFlFfi= R.

 

 

В случае циклических деформаций с частотой ш получаем для

комплексного модуля упругости М(о)

следующее

выражение:

М (ш) = а — -^-(1 +

t©T|)

 

где tç = .(Lp)_1, L — кинетический коэффициент в уравнении типа

Î = - I ( P I + Y«)

Используя приведенные выражения и метод расчета, изложен­ ный в работе [8 6, 87], получаем выражения для механических ха­ рактеристик системы через ее флюктуационные характеристики:

feT___1

 

1

 

 

V

' И3 '

1 - R3

 

(VII. 55)

 

 

 

 

 

 

kT

1

 

R2

<0Т£

 

V и3 l — R2 1 + (Û2T |

(VII. 56)

JH г (ш ) =

 

 

 

 

 

V

' и3

l - R 3

kT

(VII. 57)

 

 

 

 

 

 

 

(VII. 58)

fe_____________

*

o3R3 A2 = R3

(VII. 59)

V ' ü3 *

1 -

R3

 

kT

 

 

 

 

R3

 

(DT

(VII. 60)

 

 

 

 

 

 

(VII. 61)

Здесь Mi ( и ) — динамический модуль упругости, М2 (ю) — мни­ мая часть комплексного модуля упругости М(со), Ми— нерелак-

сированный по параметру g модуль упругости, MR — релаксированный модуль упругости, Ai — дефект модуля упругости, Д2— от­ носительный дефект модуля упругости, tg Ô — принимаемый обыч­ но за меру внутреннего трения тангенс угла механических потерь, т\— время релаксации параметра g.

Приведенные формулы показывают, что механические характе­ ристики и внутреннее трение системы в рассматриваемом случае полностью определяются квадратичными флюктуациями выбран­ ных переменных и коэффициентом корреляции этих переменных. Из этих формул также следует, что диссипация упругой энергии в сис­ теме, дефекты упругих модулей равны нулю в случае статистиче­ ской независимости упругой деформации и и внутреннего парамет­ ра g (коэффициент корреляции в этом случае равен нулю). Если статистическая зависимость деформации и внутреннего параметра

вырождается

в линейную связь

между

этими

переменными

(R = zЬ1 ),

диссипация энергии и

дефекты

модулей

достигают

максимума.

Время релаксации

также

зависит

от

параметра

корреляции

R,

возрастая при R -* 1 . В

полимерах

флюктуации

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

Аналогичным методом рассматривается взаимосвязь флюктуа­ ций с диэлектрическими характеристиками материалов, в том чис-

. ле и полимерных [87]. Результатом такого рассмотрения являются следующие выражения, определяющие диэлектрические характе­ ристики полимера через квадратичные флюктуации напряженности электрического поля Е, вектора электрической индукции D, пара­ метра корреляции R между внутренним параметром системы g (количество полярных групп, например, с определенной ориента­ цией) и напряженностью электрического поля Е:

е, (и) =

 

 

(VII. 62)

kT

1

R2

cûTg

V

* Ег

' 1— R2 '

1 + сА|

 

 

 

( V II . 63)

 

 

 

( V II . 64)

 

 

 

(VII. 65)

(VII. 66)

(Vir. 67)

(VII. 68)

где ei (со)— динамическое значение диэлектрической проницае­ мости; еи, ел— релаксированное и иерелаксированное по параметру | значение диэлектрической проницаемости; Д4— дефект диэлек­ трической проницаемости; Д2 — относительный дефект диэлектриче­ ской проницаемости; tgô 8 — тангенс угла диэлектрических потерь; х%— время релаксации параметра |; т — время релаксации, опреде­ ляющее величину диэлектрических потерь.

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

Подобная схема рассмотрения позволяет установить взаимо­ связь механических и электрических характеристик системы [2 0]. В случае одного релаксирующего параметра из расчетов, приве­ денных в работе [2 0], вытекают следующие простые соотношения:

а \мп —

( “ )] 8 г ( й )

=

& [в| (S ) —М2euJ( а )

(VII. 69)

[ М „ -

М , (« .)]

[е , (Я )

-

е4„я] =[0 , ( а ) - р „ ] 2

(VII. 70)

 

М2( а )

е2 ( а )

=

-4 яр| ( а )

(VII. 71)

где (Ù— частота изменения упругих деформаций; Q — частота из­ менения электрического поля; р — пьезоэлектрические коэффи­ циенты.

Остальные обозначения те же, что и в предыдущих выраже­ ниях. Проверка формул (VII. 69) — (VII. 71) показала их хоро­ шую согласованность с имеющимися экспериментальными дан­ ными [2 0].

ЛИТЕРАТУРА

1.В. А. Каргин, Г. Л. Слонимский, Краткие очерки по физико-хнмии полимеров, Изд. «Химия», 1967.

2. А. Т о бо ль с ки й, Свойства и структура полимеров, "Изд. «Химия», 1964.

3.Т. Алфрей, Механические свойства высокополимеров, ИЛ, 1952.

4.Л. Трелоар, Физика упругости каучука, ИЛ, 1953.

5.Г. М. Бартенев, Ю. В. Зеленев, сб. «Релаксационные явления в твер­ дых телах», под ред. В. С. Постникова, Изд. «Металлургия», 1968.

6.Г. М. Бартенев, Ю. С. Зуев, Прочность и разрушение высокоэластич­

7.

ных материалов, Изд. «Химия», 1969.

 

Структура макромоле­

В. Н. Ц в е тк ов, В. Э. Эскин, С. Я. Френкель,

8.

кул в растворах, Изд. «Химия», 1964.

 

 

 

А. А. Т а г е р, Физико-химия полимеров, Изд. «Химия», 1968.

телах», Изд.

9.

Г. П. Михайлов, сб. «Релаксационные

явления в

твердых

 

«Металлургия», 1968: Г П. Михайлов,

Т. И. Б о р ис о ва ,

Диэлектриче­

 

ские потери и поляризация полярных полимеров и

сополимеров, там же,

Г. П. Михайлов, С. П. Кабин, Б. И. Сажин, ЖТФ, 5, 4, 591 (1955).

10.В. Е. Гуль, В. Н. Кузнецов, Структура и механические свойства, поли­ меров, Изд. «Высшая школа», 1966.

11. Ч. Тен форд, Физическая химия полимеров, Изд. «Химия», 1965.

12.Д. Ферр и, сб. «Переходы и релаксационные явления в полимерах», Изд. «Мир», 1968.

13.М. В. В олькенштей н, Конфигурационная статистика полимерных цепей, Изд. АН СССР, 1959.

14.T. М. Бирштейн, О. Б. Птицын, Конформации макромолекул, Изд.

 

«Наука», 1964.

 

Polymer Chemistry, N. Y., 1953.

15. P. J. F 1 о г у , Principles of

16. A. C h a r l e s b y , J. Polymer Sci., 11, 513 (1953);

J. Proc. Roy. Soc., A231,

.

521

(1955);

A241, 495

(1957).

 

 

 

17. Л. С. Присс, Колл, ж., 19, 607 (1958).

 

 

 

18. Г. Л. С ло ним ски й, Доклад на I Всесоюзном сойещанни по релаксацион­

 

ным явлениям* в полимерах, МГПИ, 1967.

 

Бартенев, Укр. физ. ж.,

19. А. П. М о л о т к о в ,

10.

В. Зеленев, Г. М.

 

12, 303

(1967).

 

сб. «Релаксационные

явления в твердых телах»,

20. Л.

А.

С ве р гу н ен к о,

 

Изд. «Металлургия»,

1968; Л. А. Све р гу не нк о ,

Уч. зап. Ереванского гос.

 

ун-та,

2, 16 (1968); Л. А. Све р гу н ен к о,

В. Ф. Ушаков, Доклад

 

на I Всесоюзном совещании по релаксационным явлениям в полимерах,

 

МГПИ, 1967.

 

 

сб. «Релаксационные явления

21. Т. Д. Ш ер мер гор, Б. И. Да ринский,

 

в твердых телах», Изд. «Металлургия», 1968.

Т. Сырников, Основы моле­

 

22. И. Г. Михайлов,

В. А. Со л овь ев , 10.

 

кулярной акустики, Изд. «Наука», 1964.

 

 

 

 

23. И. П. Бородин, Ю. В. Зеленев, ВМС, 10, 1486 (1968).

 

24. С. И. Мешков, Г.

Н.

П а ч е в с к а я, В. С.

П о с тн и ко в , Доклад на

 

I Всесоюзном совещании по релаксационным явлениям в полимерах, МГПИ,

 

1967.

 

 

 

 

 

на

I Всесоюзном совещании

 

25. Ю. В. Зеленев, И. П. Бородин, Доклад

 

по релаксационным явлениям в полимерах, МГПИ, 1967.

26. T. H. X а з а н о в и ч,

Доклад на I Всесоюзном совещании по релаксацион­

 

ным явлениям в полимерах, МГПИ, 1967.

 

 

 

27.Ю. Я. Готлиб, А. А. Даринский, Физика тв. тела, 11, 1717, 1969.

28.P. Е. Rouse, J. Chem. Phys., 21, 1272 (1953).

29.В. H. Zi mm, J. Chem. Phys., 24, 269 (1956).

30.M. Mooney, J. Polymer Sci., 34, 599 (1959).

31.Ю. Я. Готлиб, Доклад на I Всесоюзном совещании по релаксационным

явлениям в полимерах, МГПИ, 1967.

32. Ю. Я. Готлиб, М. В. В олькенштейн, ДАН СССР, 39, 821 (1953).

33.Настоящий сборник, стр. 263—297.

34.Ю. Я- Готлиб, Физика тв. тела, 6, 2938 (1964).

35. Ф. П.

Г ри го рь ев а, 10. Я. Готлиб, ВМС, All, 962 (1962).

36. В. Ф.

Попов, Л. С. Присс, Труды Международной конференции по кау­

чуку и резине, Изд. «Химия», 1971.

37. H. Н.

Б о г о л ю б о в , Проблемы динамической теории в статистической фи­

38.

зике,

Гостехиздат,

1946.

статистическая

механика,

Йзд.

«Мир»,

И. П р и г о ж и н,

Неравновесная

39.

1964.

 

 

 

 

 

 

R. В rout, J. P r i g o g i n e , Physica, 22, 621 (1956).

 

 

40.

J. P r i g о g i n e, R. B a 1e s c u, Physica, 25, 281 (1959).

 

46,-53

41.

R. R e s i b o u i s , J. P ri go gi ne,

Acad. Roy. Belg. Bull. Classe Sci.,

 

 

 

 

Квантовая

теория поля,

Гостехиздат,

1952.

43.P. Б p a y T , Фазовые переходы* Изд. «Мир», 1967.

44.Р. Б а л е с к у, Статистическая механика заряженных частиц, Изд. «Мир»,