Объем и площадь n-мерной сферы

 

На основании размерности для объема n-мерной сферы радиусом r и шарового слоя толщиной dr получаем

 

,

 

.

 

Для нахождения постоянной вычисляем по всему пространству интеграл

.

 

В декартовых координатах

 

,

 

,

тогда

,

где использован интеграл Пуассона

.

 

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

 

,

тогда

,

где учтено

.

 

Гамма-функция является обобщением факториала и удовлетворяет соотношениям:

Г(n + 1) = n!, ,

 

Г(z + 1) = z Г(z),

 

,

 

, , ,

 

при .

 

Сравниваем выражения в декартовых и сферических координатах:

 

, ,

находим

.

 

В результате объем n-мерного шара и шарового слоя

 

, (П.2.1)

 

. (П.2.2)

Площадь сферы

. (П.2.3)

 

Эллипсоид с полуосями удовлетворяет уравнению

.

Сравниваем с уравнением сферы

,

обобщаем (П.2.1)

,

и находим объем n-мерного эллипсоида

 

. (П.2.1а)


Фазовая траектория

 

С течением времени система изменяет свое микросостояние за счет движения частиц, и точка X перемещается по фазовой траектории согласно уравнениям Гамильтона (2.1)

 

,

 

.

 

Фазовый ансамбль

 

Макросостояние системы характеризуется макроскопическими и термодинамическими характеристиками, в частности:числом частиц N, температурой T,объемом V, давлением P, внутренней энергией U, энтропией S. Одному макросостоянию соответствует множество различных микросостояний. Все они находятся в пределах некоторой области фазового пространства, границы которой зависят от макрохарактеристик. Фазовые траектории микросостояний не выходят за пределы указанной области фазового пространства. Фазовый ансамбль есть множество микросостояний с одинаковыми макрохарактеристиками, т. е. относящихся к одному макросостоянию.

 

Функция распределения микросостояний фазового ансамбля

Точка X фазового пространства описывает микросостояние системы. В интервале

 

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

 

. (2.3)

 

Вероятность нахождения системы в интервале

 

(2.3а)

 

удовлетворяет условию нормировки

 

, (2.4)

 

где интегрирование ведется по всему фазовому пространству.

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

Теорема Лиувилля

Равновесный газ описывается стационарным гамильтонианом и постоянными термодинамическими параметрами. Макросостояние реализуется фазовым ансамблем микросостояний. Это множество точек с течением времени движется по фазовому пространству. Закон их перемещения описывает теорема Лиувилля – при движении точек фазового ансамбля плотность микросостояний в каждой точке фазового пространства постоянна и зависит от гамильтониана.

. (2.5)

 

Аналогично течет несжимаемая жидкость, сохраняя свою плотность. Теорему доказал французский математик Лиувилль в 1838 г. Теорема используется для получения функции распределения состояний по фазовому пространству.

 

 

Жозеф Лиувилль (1809–1882)

 

Доказательство теоремы

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

 

 

Микросостояния с плотностью входят в объем и выходят из него.

Для нахождения числа вошедших за 1с микросостояний представим микросостояние в виде жирной точки на рисунке. Число точек в единице объема равно w. Если все точки двигаются со скоростью , то за 1с через сечение пройдут состояния, которые первоначально заполняли цилиндр с длиной образующей, равной скорости. Умножаем объем цилиндра на плотность состояний, получаем число вошедших состояний

 

.

 

 

 

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

 

,

где использовано

.

 

Если с течением времени плотность изменяется, тогда в объеме появляются и исчезают состояния. За 1с в объеме появляется число состояний

.

 

Каждое состояние описывает реальную систему, поэтому число состояний сохраняется и выполняется уравнение баланса

 

«число появившихся состояний» =

= «число вошедших состояний» – «число вышедших состояний»:

 

.

Сокращаем подобные

.

 

Результат обобщаем на случай изменения координат фазового пространства

.

 

Раскрываем круглые скобки

 

.

 

Последняя скобка равна нулю согласно уравнениям Гамильтона (2.1)

 

, .

 

Используем формулу для полной производной

 

,

 

и получаем теорему Лиувилля

 

 

– полная производная по времени от плотности микросостояний равна нулю. Следовательно, плотность микросостояний фазового ансамбля не изменяется при его движении и зависит от гамильтониана.

 

Примечание: Полный дифференциал функции

 

Рассмотрим функцию . Если изменение функции при переходе между точками

,

 

находящимися в противоположенных вершинах параллелепипеда, не зависит от формы пути, то такая функция называется потенциальной. Выбираем путь, состоящий из трех участков, параллельных осям x, y и z. Тогда изменение складывается из изменений на каждом участке. Если участки бесконечно малые, то такое элементарное изменение функции называется полным дифференциалом

 

.

 

Деление результата на , дает полную производную

 

.

 

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

Потенциальными функциями описываются гравитационное поле и электростатическое поле, магнитное поле описывается непотенциальными функциями. В термодинамике потенциальными функциями являются внутренняя и свободная энергии, непотенциальными функциями описываются работа и теплота.

 

Следствия теоремы

 

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

 

,

 

изменяется лишь форма объема. Аналогично ведет себя несжимаемая жидкость. Учитываем

,

где J – якобиан преобразования между начальными и текущими X координатами, и получаем

 

= 1. (2.6)

 

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

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

 

. (2.6а)

 

Б. Для стационарной системы функция распределения не изменяется с течением времени согласно теореме Лиувилля, и может зависеть только от интегралов движения. Если система как целое неподвижна и не вращается, то функция распределения зависит от полной энергии, т. е. от гамильтониана:

. (2.6б)

 

В. Для равновесной изолированной системы

 

.

 

Система с равной вероятностью обнаруживается в любом из доступных микросостояний.

 

Г. Теорема не выполняется для диссипативных систем, т. е. при наличии трения и неупругих соударений. Диссипативная сила, действующая на тело со стороны среды, направлена против скорости движения тела относительно среды. Уравнения Гамильтона в виде (2.1) в этом случае не применимы.

 

ПРИМЕР

 

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

 

1. Энергия системы фиксирована, поэтому микросостояние в фазовом пространстве движется по гиперповерхности. Для одномерной системы из одной частицы координаты фазового пространства (x,p). Гамильтониан осциллятора приравниваем полной энергии

 

где

; ;

 

κ – коэффициент жесткости пружины.

 

2. Получаем уравнение фазовой траектории, по которой движется микросостояние:

.

 

Сравниваем с уравнением эллипса

 

,

находим полуоси

, .

 

Микросостояния отличаются друг от друга начальной фазой.

 

 

3. Находим число микросостояний (2.2а)

 

.

 

При и интеграл равен площади эллипса

 

,

тогда число микросостояний

 

, (П.2.4)

где . Следовательно, энергия осциллятора квантуется

, (П.2.4а)

 

где квант энергии; – число микросостояний равно числу квантов энергии осциллятора.

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

 

 

4. Для получения якобиана

 

найдем функции

, ,

 

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

Используем уравнения Гамильтона (2.1)

 

, .

Подставляем гамильтониан

,

получаем

– связь скорости с импульсом,

 

– 2-й закон Ньютона ,

 

где – коэффициент жесткости упругой силы F;

 

.

 

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

,

 

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

 

.

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

,

 

.

 

Для нахождения свободных параметров A и B накладываем начальные условия

,

 

,

получаем

, .

 

Находим закон изменения координат микросостояния с течением времени

 

,

 

.

 

Микросостояние перемещается по эллипсу по часовой стрелке с круговой частотой ω.

 

5. Вычисляем якобиан

 

.

 

Теорема Лиувилля выполняется.