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

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

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

Функция вероятности ошибок (интеграл вероятности):

 

2

x

 

Φ(x) = erf x =

et2dt .

(П.2.8)

π

 

0

 

Чтобы получить приближенное значение Φ(x), подынтегральную функцию представляютстепеннымрядомспоследующиминтегрированиемслагаемых:

 

2

 

x

3

 

x

5

 

n

x

2n+1

 

 

 

 

 

 

 

(1)

 

 

Φ(x) =

 

x

 

 

+

 

 

... +

 

 

 

+... .

 

 

 

 

 

 

 

 

 

 

 

1! 3

 

2! 5

 

n ! (2n +1)

 

 

π

 

 

 

Через функцию вероятности вычисляются интегралы

et2dt = π(1 ±Φ(x))2 ;

±x

t2et2dt = π(1 Φ(x)) 4 ±xex2 2 .

±x

(П.2.9)

(П.2.10а)

(П.2.10б)

П.2.3. Интеграл Ферми. ИнтеграломФерминазывается интегралвида:

I = g(ε) f ((εµ)kT )dε,

0

где f (εµkT )= (e(ε µ )/k T +1)1 – функция распределения Ферми-Дирака,

g(ε)монотонно возрастающая

функция, равная нулю при

ε = 0 , а

µ kT >> 1. Функция f(z) , где

z = (εµ) kT , при заданных

условиях

близка к единице, если ε < µ, и стремится к нулю (~ ez ) при z → ∞. Ее

производная f

ε– симметричная функция z

с резким

максимумом в точке

ε = µ, а при

других значениях ε близка

к нулю.

После вычисления

 

 

 

ε

интеграла по частям и представления функции F(ε) = g(x)dx рядом

0

Тейлора вблизи точки ε = µ имеем

282

 

 

 

 

 

 

 

 

 

 

 

 

m

F

(m)

(µ)

 

 

m

z

dz

 

 

 

 

 

 

 

 

 

 

I = F(µ) +

(kT)

 

 

 

 

z

e

 

 

.

 

 

 

 

 

 

 

 

m !

 

 

 

 

 

z

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

m=1

 

 

 

 

 

µ kT

(e

+1)

 

 

 

 

 

Так как при больших отрицательных значениях z

 

подынтегральная функция

близка к нулю, а

µ kT >> 1,

то с высокой степенью точности нижний

предел

интеграла

 

можно

 

положить

 

 

равным

 

 

 

(−∞).

Поскольку

ez (ez +1)2

 

=ez (ez +1)2 , имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zmezdz

0

 

zmezdz

 

 

zmezdz

 

 

 

[1 +(1)m ] zmez

 

 

 

=

 

 

 

 

 

+

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

dz ,

(e

z

2

(e

z

2

 

(e

z

+

 

2

 

 

 

(e

z

+

2

 

− ∞

 

+1)

 

− ∞

 

+1)

 

 

0

 

1)

 

 

 

0

 

 

 

 

 

 

1)

 

 

т.е. слагаемые с нечетным

 

значением

m

 

исчезают,

а

интеграл I

преобразуется к виду:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(kT)2m F(2m)(µ) z2mez

 

 

 

 

 

 

 

 

 

 

I = F(µ) +2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz .

 

 

 

 

 

 

 

(2m)!

 

 

 

 

0 (1 +e

z 2

 

 

 

 

 

 

 

 

 

 

 

m=1

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку (1 +ez )2

= (1)n1ne(n1) z , то после замены z = y n имеем

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

2m

z

 

 

 

 

 

n1

 

 

 

 

 

 

 

 

 

 

 

 

n1

 

 

 

 

 

e

 

 

dz =

(1)

 

 

 

 

 

 

y

 

2m

 

 

(2m)!

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

y

d y =

 

 

.

 

 

(1 +e

z 2

 

 

n

2m

 

 

n

2m

 

 

)

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, вычисление фермиевского интеграла сводится к формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

 

 

 

 

 

 

 

 

 

 

 

I = F(µ) +2(kT)2m F(2m)(µ)

(1)

 

.

 

 

(П.2.11)

 

 

 

 

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m=1

 

 

 

 

 

 

 

 

n=1

 

n

 

 

 

 

 

 

Если

ограничиться

 

первым

 

(m = 1)

 

слагаемым

в

 

сумме, учитывая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)n1

n2 = π2 12 , получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(ε)dε

 

 

µ

 

 

 

 

 

(πkT)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

= g(ε)dε+

 

 

 

 

 

 

 

gε(µ) = F(µ) +F

 

(µ)(πkT) 6. (П.2.12а)

 

(εµ) kT

 

 

 

 

 

6

 

 

 

0 e

 

 

+1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для частного случая: g(ε) = εs (s > 0 ), имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

µ

s+1

 

 

 

 

s(s +1)

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε dε

 

 

 

 

 

 

 

 

 

 

 

 

πkT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

1 +

 

 

 

 

.

 

 

(П.2.12б)

 

 

 

 

 

 

 

(εµ) kT

+1

 

 

s +1

 

 

 

6

 

 

µ

 

 

 

 

 

 

 

 

 

 

 

 

 

0 e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

283

П.3. ТЕОРЕМА ЛИУВИЛЛЯ

Выделим малый элемент фазового пространства:

(

)

6N

 

k

1 2

6

 

 

dx

.

(П.3.1)

dx

6N =

 

 

= dx dx

...dx

k=1

Точки внутри него, изображающие различные механические системы, не могут возникать и исчезать, как и сами механические системы, которые они изображают. Изменение их числа в рассматриваемом элементе фазового объема происходит только за счет втекания этих точек из других частей пространства или их вытекания из элемента фазового объема в другие части пространства 6N измерений.

Рассмотрим элемент объема фазового пространства с двумя его поверхностями S1 = 1 и S2 = 1 , ортогональными к некоторому выделенному направлению xi (Рис. П.3.1). Число точек, вошедших вдоль оси xi через поверхность S1 с координатой xi за единицу времени, равно ρ(x,t)xiS1 , где xi – проекция их скорости на направление xi , ρ – плотность точек фазового пространства. Через вторую поверхность S2 с координатой xi +dxi вытечет число точек, равное:

ρxi + (ρxxi i ) S2 ,

где второй член в скобках учитывает изменение величин ρ и xi на участке dxi .

Рис. П.3.1. К выводу теоремы Лиувилля

284

Разность между числом точек, втекающих через поверхность S1 и

вытекающих через – S2 при их одинаковой площади (S2 = S1 = 1 ),

(ρxi )

(П.3.2)

xi

 

 

равна изменению числа точек внутри выбранного объема вследствие их движения вдоль одного направления xi в единицу времени.

Но в фазовом пространстве имеется 6N независимых направлений. Через каждое такое направление так же, как и через рассмотренное, могут втекать и вытекать точки фазового пространства. Поэтому полное изменение числа точек в элементе фазового пространства в единицу времени определяется суммой слагаемых (П.3.2)

ρ

6N

(ρxi )

 

 

+

= 0 .

(П.3.3.)

t

 

i=1

x

i

 

 

 

 

 

 

Для трехмерного пространства это – уравнение непрерывности:

ρt +div(ρu) = 0 .

Значит, полученное уравнение (П.3.3) можно рассматривать как уравнение непрерывности (сохранения), обобщенное на 6N -мерное пространство обобщенных координат и импульсов.

В представлении через обобщенные координаты и импульсы оно принимает вид

ρ

3N

ρ

 

ρ

3N

 

q

i

 

p

 

 

 

 

 

 

+ pi

 

 

 

 

 

+

i

 

= 0 .

(П.3.4)

 

 

 

 

 

 

 

 

t

+ qi

q

 

 

 

+ ρ

q

 

 

 

 

i

 

p

i=1

 

i

 

p

 

 

 

i=1

 

 

i

 

 

 

i

 

 

Первые два члена представляют собой полную производную от ρ по времени

(с учетом того, что ρ есть функция q , p , t ), т.е. dρdt . Поэтому последнее можно записать так:

dρ 3N

 

q

i

 

p

 

 

 

 

 

 

 

i

 

 

 

+ ρ

 

 

+

 

 

= 0 .

(П.3.5)

dt

q

 

 

i=1

 

i

 

p

 

 

 

 

 

 

i

 

 

285

Здесь пока отражены только общие кинематические соотношения. Для консервативной системы, у которой обобщенные координаты и импульсы определяются через функцию Гамильтона ( pi = −∂Hqi и qi = ∂Hpi ),

из (П.3.5) следует

dρ

= 0

(П.3.6)

dt

 

 

Равенство (П.3.6) справедливо только для точек, перемещающихся вдоль фазовой траектории, а вообще в разных местах фазового пространства плотность числа точек может быть различна и может меняться со временем, что видно из зависимости ρ(q, p,t).

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

жидкость. Уравнение dρdt = 0 вдоль фазовой траектории и есть теорема

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

Как следствие теоремы Лиувилля, элементарный объем в фазовом пространстве, перемещаясь с течением времени, будет оставаться постоянным по величине, хотя его форма может меняться, т.е.

(

)0

(

)t

(П.3.7)

dx

6N

= dx

6N

286

ВОПРОСЫ ДЛЯ САМОКОНТРОЛЯ

РАЗДЕЛ 2

1.1.Фазовий простір (ФП).

1.2.Характеристики системи у ФП (елементарний об’єм, траекторія, точка зображення).

1.3.Квазікласична система. Число станів.

1.4.Метод ансамблів Гіббса.

1.5.Теорема Ліувілля. Рівняння для функції розподілу.

1.6.Залежність стаціонарної функції розподілу консервативних систем від інтегралів руху.

1.7.Визначення систем за умов їх взаємодії з навколишнім середовищем.

1.8.Мікроканонічний та канонічний розподіли Гіббса.

1.9.Функція станів.

1.10.Умови використання й особливості розподілу Гіббса.

1.11.Великий канонічний розподіл.

1.12.Хімічний (парціальний) потенціал.

1.13.Фазова траекторія лінійного гармонійного осцилятора.

РАЗДЕЛ 2

2.1.Визначення ідеального газу (ІГ).

2.2.ІГ як ансамбль Гіббса.

2.3.Розподіл Максвелла.

2.4.Молекулярно-кінетична модель тиску та рівняння стану ІГ.

2.5.Особливості розподілу Максвелла.

2.6.Мікрота макрохарактеристики ІГ.

2.7.Розподіл енергії класичних частинок за ступенями вільності.

2.8.Зіткнення молекул ІГ між собою, перетин зіткнення, довжина вільного пробігу.

2.9.Розподіл молекул ІГ у зовнішньому полі.

2.10.Число та функція станів ІГ.

287

2.11. Діелектрична проникність ІГ полярних молекул (діелектрик другого роду).

2.12. Температурна залежність діелектричної проникності ІГ полярних молекул.

2.13. Фізичний образ і визначення густини струму термоелектронної емісії. 2.12. Чим зумовлена природна ширина спектральної лінії випромінювання активного середовища?

РАЗДЕЛ 3

3.1.Основний постулат статистичної термодинаміки. Внутрішня енергія макроскопічної системи.

3.2.Статистичне обґрунтування першого закону термодинаміки.

3.3.Молекулярно-кінетична трактовка ентропії і її властивості.

3.4.Квазістатичний (зворотній) процес і умови його реалізації.

3.5.Формула Больцмана. Другий закон термодинаміки.

3.6.Слідства другого закону термодинаміки. Максимальна робота теплових машин.

3.7.Третій закон термодинаміки та його слідства.

3.8.Рівняння стану ідеального газу в адіабатичному процесі.

3.9.Визначення та фізичний сенс термодинамічних потенціалів.

3.10.Критерій оборотності ізотермо-ізохорного та ізотермо-ізобаричного процесів у незамкнених системах.

3.11.Диференціали термодинамічних потенціалів як функції узагальнених сил і параметрів.

3.12.Залежність термодинамічних потенціалів від функції станів.

3.13.Визначення узагальнених сил.

3.14.Основна термодинамічна рівність системи зі змінним числом частинок.

3.13. Базові співвідношення термодинамічної теорії діелектриків і магнетиків. Визначення електричного та магнітного векторів поляризації.

3.16.Алгоритм отримання диференціальних (перехрестних) співвідношень термодинаміки, їх методичне й практичне значення.

3.17.Рівняння Гіббса-Гельмгольца.

288

3.18.Узагальнююча формула теплоємності та її подання для окремих процесів.

3.19.Теплові та пружні коефіцієнти і зв’язок між ними.

3.20.Необхідні і достатні умови стійкої рівноваги однорідної системи.

3.21.Закон Кюрі для класичних парамагнетиків.

3.22.Отримання формули для ентропії ідеального газу (Сакура-Тетроді).

3.23.Коефіцієнт корисної дії теплових машин оборотного циклу.

РАЗДЕЛ 4

4.1.Метод комірок у фазовому просторі.

4.2.Класифікація частинок у квантовій механіці.

4.3.Розподіли Фермі-Дірака та Бозе-Ейнштейна.

4.4.Принципи побудови класичної і квантової теорії теплоємності ідеального газу двоатомних молекул.

4.5.Енергія коливань і обертань квантових частинок.

4.6.Функція коливальних і обертальних станів двоатомних молекул.

4.7.Молекулярна модель кристала, класична теорія теплоємності.

4.8.Внутрішня енергія та функція станів кристала.

4.9.Основні положення теорії теплоємності Дебая.

4.10.Фізичне обґрунтування правомірності допущення в теорії теплоємності Дебая.

4.11.Асимптотика при високих і низьких температурах.

4.12.Утворення електронного газу у металах. Функція розподілу.

4.13.Властивості електронного газу поблизу абсолютного нуля температури.

4.14.Енергія Фермі.

4.15.Температурна залежність хімічного потенціалу електронного газу.

4.14. Рівноважне випромінювання (РВ) нагрітих тіл.

4.17.Фізичний образ абсолютно чорного тіла.

4.18.Спектральна густина випромінювання (визначення, алгоритм обчислення, асимптотика при високих та низьких температурах).

4.19.Функція розподілу та властивості фотонного газу.

4.20.Термодинамічні функції та параметри РВ.

РАЗДЕЛ 5

289

5.1. Фізичний сенс флуктуаційних процесів у макросистемах і їх характеристики.

5.2.Міра ймовірності флуктуацій внутрішніх параметрів замкнених систем.

5.3.Міра ймовірності флуктуацій у квазізамкнених системах.

5.4.Загальний метод обчислення ймовірності малих флуктуацій термодинамічних параметрів.

5.5.Яким чином обчислюється міра флуктуацій двох незалежних термодинамічних параметрів однорідної макросистеми?

5.6.Кореляційний момент і радіус (коефіцієнт) кореляції.

5.7.Що є підґрунтям узагальнених умов стійкості рівноваги?

5.8. Визначення масштабу флуктуацій термодинамічних параметрів у однорідному середовищі.

5.9.Модель флуктуацій в струмових приладах.

5.10.Принцип визначення дисперсії дробового шуму.

5.11.Фізична модель теплового шуму, дисперсія відповідного струму.

5.12.Шляхи зменшення впливу флуктуацій струму на вимірювальні прилади.

5.13.Залежність дисперсії струму від ширини частотної смуги пропускання вимірювального приладу.

РАЗДЕЛ 6

6.1.Явища релаксації і переносу. Час релаксації.

6.2.Кінетичні коефіцієнти.

6.3.Фізична основа кінетичного рівняння Больцмана.

6.4.Визначення потоків фізичної якості.

6.5.Як обчислюється дрейфова похідна в рівнянні Больцмана?

6.6.Фізичні допущення наближеного рішення кінетичного рівняння Больцмана.

6.7.Базова формула обчислення електропровідності.

6.8.Фізичні допущення та математичний алгоритм отримання коефіцієнта та кута Холу.

290

СПИСОК ЛИТЕРАТУРЫ

1.Анго А. Математика для электро- и радиоинженеров // пер. с франц. под ред. К.С. Шифрина. – М.: Наука, 1967. – 780 с.

2. Ансельм А.И. Основы статистической физики и термодинамики.

М.: Наука, 1973. – 424с.

3.Базаров И.П. Термодинамика. – М.: Высш. Шк., 1976. – 447 с.

4.Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеров и учащихся втузов. – М.: Наука, 1981. – 720 с.

5.Біленко І.І. Фізичний словник. – К. Вища шк., 1976. – 447 с.

6.Варикаш В.М., Болсун А.И., Аксенов В.В. Сборник задач по статистической физике. – Минск: Вышэйшая шк., 1979. – 223 с.

7.Васильев А.М. Введение в статистическую физику. – М.: Высш. шк., 1980. – 272 с.

8.Вентцель Л.Г. Теория вероятностей. – М.: Наука, 1969. – 576 с.

9.Головко Д.Б., Менторський Ю.Л. Загальні основи фізики. Термодинаміка.

Молекулярна фізика. Навч. посібн. – К.: Либідь, 1993. – 112 с.

10. Задачи по термодинамике и статистической физике// под ред. П. Ландсберга. – М.: Мир, 1974. – 640 с.

11.Жирифалько Л. Статистическая физика твердого тела. – М.: Мир, 1975. – 382 с.

12. Казанский В.Б.

Методические

указания

к

решению

задач

по

курсу

«Термодинамика

и

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

физика».

Харьков: ХГУ, 1984. – 50 с.

13.Казанский В.Б. Статистическая физика и термодинамика. Задачи, основные понятия и положения. Метод. пособие. – Харьков: ХНУ, 2004. – 112 с.

14.Киттель Ч. Статистическая термодинамика. – М.: Наука, 1977. – 336 с.

15.Ландау Л.Д., Лифшиц Е.М. Статистическая физика. – М.: Наука, 1964. –

291