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

Метода термодинамика(лекции)

.pdf
Скачиваний:
37
Добавлен:
23.02.2015
Размер:
2.65 Mб
Скачать

ней излучение. Нетрудно найти температуру, при которой плотность энергии

излучения

σT 4

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

энергии

теплового движения

частиц 3kTn 2 .

Из равенства σT 4

= 3kTn 2

следует

T = 3 3kn σ . Для

достаточно

разреженной плазмы,

когда

плотность

частиц

n 1021 м-3,

T 3 104

К. При температурах полной

ионизации

(T 105

К) плотность

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

Электростатический потенциал плазмы. Рассмотрим модель плазмы в предположении ее полной ионизации, т.е. состоящей из одинакового числа положительно и отрицательно заряженных частиц (N+ = N= N / 2 ). При этом свободной энергией излучения пренебрегаем. Внутренняя энергия плазмы E складывается из кинетической энергии хаотического движения ее

частиц Eид = 3NkT 2 =CVT (внутренняя энергия идеального

газа) и

 

N

 

энергии электростатического взаимодействия Eq

= ej ϕ(rj ) 2

 

 

j =1

 

N

 

 

E = Eид + ej ϕ(rj )

2 ,

(4.64)

j =1

 

 

где ϕ(rj ) – потенциал поля всех остальных зарядов в точке нахождения j -го заряда. Энергия одного положительного (отрицательного) заряда в произвольной точке (например, r0 = 0 ), равна eϕ+(0) (eϕ(0) ). Тогда за суммарную энергию кулоновского взаимодействия можно принять

E

q

=

1 Neϕ

(0)

1 Neϕ

(0) .

 

(4.65)

 

 

2

 

+

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно методу самосогласованного поля, распределение плотности

положительных и отрицательных зарядов определяется формулой

n (r )

= n exp ±e ϕ[n(r )] kT

,

(4.66)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где n0 = N /V , ϕ(r ) –

потенциал в

точке

r ,

который

порождается

остальными зарядами. После подстановки (4.66) в уравнение Пуассона имеем

212

 

 

 

 

 

 

ϕ(r ) = −4πρ(r ) = 4πen0

eϕ kT

e

eϕ kT

 

(4.67)

 

 

 

 

 

 

e

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

 

достаточно

разреженной

изотропной

 

плазмы

(eϕ << kT ,

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∆ =

 

 

r

 

 

) уравнение упрощается

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ = χ2ϕ,

 

χ2

=

 

8πe2n0

.

 

 

 

(4.68)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kT

 

 

 

 

 

 

Общее решение

ϕ(r) =C eχr r +C eχr

r

с учетом его ограниченности во

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

всем пространстве (ϕ(r → ∞) < M ) принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

ϕ(r) =C

1

eχr

 

r .

 

 

 

 

 

(4.69)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

вблизи

точки

r = 0 ,

где

находится положительный или

отрицательный заряд, потенциал ϕ(r 0) = ±e r ,

то C1

= ±e . Потенциал,

создаваемый всеми остальными зарядами в окрестности

точки r = 0 ,

определяется разностью между ϕ(r) и ±e / r

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ (0)

= lim e ±eχr 1 =

χe .

 

(4.70)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

±

r 0 r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя (4.70) в (4.64) и (4.65), находим энергию кулоновского взаимодействия и внутреннюю энергию плазмы

E

q

= −Ne2χ,

 

 

 

 

 

 

 

 

 

 

 

 

E =

3

NkT Ne

2 8πe2

 

N

(4.71)

2

kT

V

 

 

 

 

 

 

 

Отрицательное значение Eq обусловлено преимуществом сил притяжения,

поскольку каждый заряд находится в окружении противоположно заряженной среды (дебаевская экранировка).

4.10.3. Термодинамические потенциалы и параметры. При независимых переменных температуре T и объема V термодинамическим потенциалом плазмы есть свободная энергия F(T,V ) (3.40)

213

 

 

F = −T

 

E

dT = −T

E

+E

 

 

 

 

 

 

 

ид

 

q

dT = Fид +Fq ,

(4.72)

 

 

 

2

 

2

 

 

 

 

 

 

 

 

T

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3/2

 

 

 

 

 

 

 

 

 

 

 

 

 

2πmkT

 

 

 

 

 

 

 

 

где F

= −NkT ln

eV

– свободная энергия идеального газа, а

 

 

 

 

 

 

 

 

 

ид

 

 

 

N

 

 

h

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F = −2 Ne2

8πe2N

 

свободная

энергия

кулоновского взаимодействия.

q

3

kTV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При температурах T 105 К необходимо к (4.72) добавить свободную энергию излучения (4.58)

F ′ = F αVT3 .

Уравнение состояния ( p = −∂FV ), энтропия (S = −(FT )V ) и

теплоемкость (CV = −T (2FT 2 )V ) плазмы определяются формулами

 

 

p = NkT

1 Ne2

8πNe2

,

 

 

 

 

(4.73)

 

 

 

V

3

kTV 3

 

 

 

 

 

 

S = (C

)

lnT +Nk lnV

1 Ne2

 

8πNe2

,

(4.74)

 

V ид

 

 

3

 

 

kT 3V

 

 

C

= (C

)

+ 1 Ne2

8πNe2

, (C )

 

=

3 Nk .

(4.75)

V

V ид

2

kT 3V

V

ид

 

2

 

 

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

 

 

A

 

 

Ne2

8πNe2

 

 

 

+

 

 

,

A = −

 

 

.

(4.76)

 

 

 

pV = NkT 1

 

 

 

 

 

 

 

 

 

 

3NkT

kT

 

 

 

 

V

 

 

 

 

Здесь уже нельзя выделить, как в уравнении Ван-дер-Ваальса, разную

степень взаимодействия между частицами при изменении их плотности. 214

С помощью формул (4.73) и (4.75) находим следующие условия устойчивости разреженной плазмы

 

 

 

 

 

 

 

 

 

 

NkT

 

 

1

 

 

2

8πNe

2

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

+

 

Ne

 

 

 

 

 

 

 

< 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

V T

 

 

 

 

V

 

 

 

2

 

 

 

kTV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.77)

 

 

T

=

 

 

 

 

 

 

 

T

 

 

 

 

 

 

> 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CV

 

 

(C

)

 

 

+

1

Ne

2 8πNe2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

kT 3V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V ид

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

NkT

 

 

 

 

ид

 

 

 

T

 

 

 

 

 

T

 

 

p

 

 

 

 

 

p

 

,

 

 

 

 

 

 

(4.78)

 

 

 

 

 

> −

 

 

 

 

=

 

 

 

 

 

 

 

 

<

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V T

 

 

 

V

 

 

 

V T

 

 

CV

 

 

 

 

(CV )ид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

Выделим мысленно в плазме с концентрацией электронов n прямоугольный параллелепипед длиной dx и сечением S (объем dV = Sdx ) (рис. 4.12). Вследствие большой массы скорость движения ионов много меньше, чем у электронов, поэтому их можно считать неподвижными. Пусть

внекоторый момент времени электроны выделенного объема сместились по отношению к ионам на величину ξ(x). Тогда возникающий здесь объемный

заряд и его плотность становятся равными

215

 

Рис. 4.12

 

dq =enS ξ(x) enS ξ(x +dx)

 

≈ −en dξ

Sdx = −en dξ dV;

 

dx

 

 

dx

 

ρ =

dq

=en dξ .

(4.79)

 

 

dV

dx

 

Учитывая связь локального заряда с напряженностью электрического поля E

уравнением div E = = 4πρ, или эквивалентным ему для одномерного случая

dΕ dx = −4πen (dξ dx) получаем

 

 

E = −4πenξ .

 

(4.80)

Из уравнения движения электрона (mξ = eE ) следует

 

 

ξ + 4πne2

ξ = 0 ,

или ξ + ω2ξ = 0 .

(4.81)

 

m

 

 

 

 

0

 

 

 

 

 

 

 

 

Эти формулы и определяют продольные ( E

ξ ) электрические колебания в

плазме с циклической частотой

 

 

 

 

 

 

 

 

ω

0

=

4πne2 .

 

(4.82)

 

 

 

 

m

 

 

 

 

 

 

 

 

 

Так, при

концентрации плазмы

n ~ 1016

м–3 её собственная

частота

ω0 = 5 109

с– 1, что соответствует дециметровым волнам. Впервые наличие

колебаний в плазме было установлено Дж. Релеем (1906) и независимо И.

Ленгмюром (1929 г.), получившим формулу для ω0 . Последняя была названа ленгмюровской частотой колебаний плазмы.

216

5.ТЕОРИЯ МАЛЫХ ФЛУКТУАЦИЙ

5.1.Определение и значение флуктуаций

Флуктуации – случайные отклонения наблюдаемых значений физических величин от их средних значений. Для макроскопических систем наблюдаемые значения физических величин с очень большой точностью совпадают с их статистическими средними, а сколько-нибудь значительные флуктуации встречаются редко. Так, относительная флуктуация аддитивной величины в системе из N частиц обратно пропорционально N . Этот результат является

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

позволяет установить границы применимости термодинамических понятий и закономерностей. Кроме того, флуктуации являются существенной особенностью многих физических явлений. Например, с точки зрения термодинамики (рассматриваются только средние значения физических величин!), импульс макроскопического тела, находящегося в среде с температурой T , в состоянии равновесия равен нулю. Согласно же статистической физике, в соответствии с законом о равнораспределении энергии по степеням свободы средняя кинетическая энергия такого тела

должна быть равна 3NkT / 2 , а его третья часть совпадает с Nmvx2 / 2 . Таким образом, средний квадрат каждой декартовой компоненты скорости тела равен vx2 (vx )2 = kT / m . Это объясняет броуновское движение

взвешенных частиц, которое послужило одним из первых объектов изучения флуктуаций. Хаотические отклонения плотности частиц – причина таких явлений как критическая опалесценция, голубая окраска неба, шумы радиоаппаратуры и т.п. В работах А. Эйнштейна и М. Смолуховского (1905-

1906 гг.) построена количественная теория флуктуаций, которая оказалась в

217

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

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

вытекающим

из соотношения неопределенностей, является

неравенство

τT << 1,

где τ = x xt – время, характеризующее скорость (xt

) изменения

физической величины около среднего значения (~ время релаксации), T – температура (см. Л.Д. Ландау, Е.М. Лифшиц. Статистическая физика. М.: Наука.

– 1964, 568 с, гл. 12, § 112.) При очень низких температурах или очень быстром изменении величины x (очень мало τ) флуктуации нельзя рассматривать термодинамически, инапервыйпланвыступают чистоквантовыефлуктуации.

Явление флуктуации практически может наблюдаться в двух случаях: 1) размеры системы достаточно малы; флуктуации происходят часто и масштаб их относительно велик; 2) размеры системы большие; флуктуации происходят также часто, но отклонение системы от состояния равновесия весьма малые. Наглядным примером флуктуаций в малых системах может служить броуновское движение.

В настоящем курсе основное внимание уделяется флуктуациям физических величин в больших классических системах и не рассматриваются квантовые флуктуации.

5.2. Мера вероятности и масштаб флуктуаций

5.2.1. Флуктуации в замкнутых системах . Пусть замкнутая система в состоянии статистического равновесия имеет энтропию

S0 . При переходе в неравновесное состояние энтропия принимает значение

218

S . Считаем, что изменение состояния системы обусловлено изменением некоторого внутреннего параметра ξ (например, плотности, давлении),

значение которого зависит от состояния всей системы, а в состоянии равновесия ξ = ξ0 . Энтропия системы S(ξ) – функция параметра ξ. Согласно формуле Больцмана (S = k ln w +const ), вероятность величины ξ иметь значение в интервале между ξ и ξ +dξ равна

 

S(ξ)S(ξ0 )

 

w(ξ)dξ = const e

k

dξ.

(5.1)

 

 

 

Эта формула была впервые применена А. Эйнштейном в 1910 г. к исследованию флуктуаций. Здесь в качестве меры вероятности флуктуации рассматривается изменение энтропии. В состоянии равновесия энтропия имеет максимальное значение S0 = S(ξ0 ) , поэтому разность энтропии S(ξ) S(ξ0 ) отрицательная.

Считаем отклонение параметра ξ от его равновесного значения малым.

Представляя S(ξ) рядом по степеням ξ и ограничиваясь членом второго порядка малости, получимсучетом S(ξ0 )ξ = 0

S(ξ) = S(ξ0 ) a2 (ξ ξ0 )2 ,

где a = −∂2S(ξ0 )ξ2 > 0 . Нормировочнуюпостояннуюопределимизусловия

a

2

conste

 

(ξξ0 ) dξ = 1.

2k

−∞

Хотя здесь использовалось выражение для w(ξ), относящееся к малым отклонениям от ξ0 , но ввиду быстрого убывания подынтегральной функции с

увеличением

ξ ξ0

область интегрирования можно распространить на все

значения

от −∞ до

+∞.

После интегрирования

 

получим

const = 1

2πk a .

Таким

образом,

вероятность отклонения

ξ

от ξ0

определяется распределением Гаусса:

219

 

 

 

 

 

1

 

e

(ξξ0 )2

 

 

w(ξ)dξ =

 

 

2k a dξ,

(4.2)

 

 

 

 

 

 

 

2πk a

 

 

 

 

 

 

 

 

 

 

а величина k a = (ξ ξ0 )2 (ξ)2

 

есть

дисперсия

(средний квадрат

отклонения); ее еще называют центральным моментом второго порядка. Если в разложении S(ξ) учитывать последующие члены, то распределение вероятностей флуктуации становится негауссовым и может быть нессиметричным. В этом случае требуется, помимо дисперсии, задание моментов более высокого порядка.

 

Из

формулы следует, что отклонения

физического параметра ξ от

 

 

гораздо чаще встречается в интервале

 

 

 

 

, чем вне его.

ξ = ξ

 

ξ ξ

(ξ)2

0

 

 

0

 

 

 

Величина (ξ)2 дает представление о масштабе флуктуаций.

Знание среднего квадрата отклонения (ξ ξ0 )2 позволяет найти дисперсию для любой функции случайной величины f(ξ). Учитывая малую величину отклонения, имеем f = f(ξ) f(ξ0 ) (f(ξ0 )ξ)(ξ ξ0 ). Откуда следует

(f )2 = ∂f(ξ0 ) ξ 2 (ξ)2 .

Заметим, что приведенная формула (5.2) применима к флуктуациям в системах с постоянной энергией.

5.2.2. Флуктуации в квазизамкнутых системах . Произвольную квазизамкнутую систему можно рассматривать как малую часть замкнутой системы или как подсистему, погруженную в термостат с постоянной температурой T0 . Считаем, что флуктуации происходят только в подсистеме, тогда как термостат остается равновесной системой. Состояние

подсистемы определяется некоторым внешним параметром λ. При переходе

220

из равновесного в неравновесное состояние он меняется от λ0 до λ, при этом изменяются и термодинамические параметры, характеризующие подсистему. Предполагаем изменения λ настолько медленным, что в каждый данный момент времени подсистема находится в равновесном состоянии, и ее термодинамические параметры связаны между собой равновесными соотношениями. Причиной перехода подсистемы из равновесного в неравновесное состояние будем считать действие некоторого внешнего теплоизолированного источника работы. При изменении параметра λ на величину λ = λ λ0 источник совершает работу A(λ) .

Поскольку термостат, подсистема и источник работы составляют замкнутую систему, к ним применима формула (5.1), где общее изменение энтропии состоит из изменения энтропии подсистемы Sи термостата S0 :

S = ∆S0 +∆S.

Поэтому вероятность перехода подсистемы в состояние с λ под действием внешнего источника роботы определяется формулой

S0 +∆S

(5.3)

dw(λ) = const e k dλ.

Запишем основные термодинамические равенства для подсистемы и термостата

E′ = ∆A p0V ′ +T0S′ ∆E0 = p0V0 +T0S .

Здесь T0 и p0 – равновесные значения температуры и давления подсистемы и термостата, Eи V – энергия и объем подсистемы, A – работа внешнего источника (но не термостата!). Поскольку внутренняя энергия и полный объем замкнутой системы остаются постоянными (V0 = ∆V , E′ +∆E0 = 0 ), то из указанных равенств следует

S +∆S′ = −∆A(λ)T0 .

Подстановка этих величин в (5.3) дает

221