Русская Википедия:Статистическая сумма
Статистическая сумма (или статсумма) (обозначается <math>Z</math>, от Шаблон:Lang-de — сумма по состояниям) — это нормировочный коэффициент в знаменателе соответствующего статистического (вероятностного) распределения, при котором интегральная сумма этого вероятностного распределения (т.е. полная вероятность) по всем возможным состояниям равна 1. Статистическая сумма - важная величина в термодинамике и статистической физике, содержащая информацию о статистических свойствах системы в состоянии термодинамического равновесия. Она может являться функцией температуры и других параметров, таких как объём. Многие термодинамические величины системы, такие как энергия, свободная энергия, энтропия и давление, могут быть выражены через статистическую сумму и её производные.
Статистическая сумма в каноническом ансамбле
Определение
Предположим, что имеется подчиняющаяся законам термодинамики система, находящаяся в постоянном тепловом контакте со средой, которая имеет температуру <math>T</math>, а объём системы и количество составляющих её частиц фиксированы. В такой ситуации система относится к каноническому ансамблю. Обозначим точные состояния, в которых может находиться система, через <math>j</math> <math>(j=1,2,3,\ldots)</math>, а полную энергию системы в состоянии <math>j</math> — <math>E_j</math>. Как правило, эти микросостояния можно рассматривать как дискретные квантовые состояния системы.
Каноническая статистическая сумма — это
- <math>Z=\sum_j e^{-\beta E_j},</math>
где обратная температура <math>\beta</math> определена как
- <math>\beta\equiv\frac{1}{k_BT},</math>
а <math>k_B</math> — это постоянная Больцмана. В классической статистической механике было бы некорректно определять статистическую сумму в виде суммы дискретных членов, как в приведённой выше формуле. В классической механике координаты и импульсы частиц могут меняться непрерывно, и множество микросостояний несчётно. В таком случае необходимо провести разбиение фазового пространства на ячейки, то есть два микросостояния считаются одинаковыми, если их различия в координатах и импульсах «не слишком велики». При этом статистическая сумма принимает вид интеграла. Например, статистическая сумма газа из <math>N</math> классических частиц равна
- <math>Z=\frac{1}{N!h^{3N}}\int \exp[-\beta H(p_1,\ldots,p_N,x_1,\ldots,x_N)]\,d^3p_1\ldots d^3p_N\,d^3x_1\ldots d^3x_N,</math>
где <math>h</math> — некоторая величина размерности действия (которая должна быть равна постоянной Планка для соответствия квантовой механике), а <math>H</math> — классический гамильтониан. Причины появления множителя <math>N!</math> объяснены ниже. Для простоты в этой статье будет использоваться дискретный вид статистической суммы, но полученные результаты в равной мере относятся и к непрерывному виду.
В квантовой механике статистическая сумма может быть записана более формально как след по пространству состояний (который не зависит от выбора базиса):
- <math>Z=\mathrm{tr}\,(e^{-\beta H}),</math>
где <math>H</math> — оператор Гамильтона. Экспонента от оператора определяется с помощью разложения в степенной ряд.
Смысл и значимость
Сначала рассмотрим, от чего она зависит. Статистическая сумма является функцией температуры <math>T</math>, а также энергий микросостояний <math>E_1,E_2,E_3</math> и т. д. Энергии микросостояний определяются другими термодинамическими величинами, такими как число частиц и объём, а также микроскопическими свойствами, такими как масса частиц. Эта зависимость от микроскопических свойств является основной в статистической механике. По модели микроскопических составляющих системы можно рассчитать энергии микросостояний, а следовательно, и статистическую сумму, которая позволяет рассчитать все остальные термодинамические свойства системы.
Статистическая сумма может быть использована для расчёта термодинамических величин, поскольку она имеет очень важный статистический смысл. Вероятность <math>P_j</math>, с которой система находится в микросостоянии <math>j</math>, равна
- <math>P_j=\frac{1}{Z}e^{-\beta E_j}.</math>
Статистическая сумма входит в распределение Гиббса в виде нормировочного множителя (она не зависит от <math>j</math>), обеспечивая равенство единице суммы вероятностей:
- <math>\sum_j P_j=\frac{1}{Z}\sum_j e^{-\beta E_j}=\frac{1}{Z}Z=1.</math>
Вычисление термодинамической полной энергии
Чтобы продемонстрировать полезность статистической суммы, рассчитаем термодинамическое значение полной энергии. Это просто математическое ожидание, или среднее по ансамблю значение энергии, равное сумме энергий микросостояний, взятых с весами, равными их вероятностям:
- <math>\langle E\rangle=\sum_j E_jP_j=\frac{1}{Z}\sum_j E_j e^{-\beta E_j}=-\frac{1}{Z}\frac{\partial}{\partial\beta}Z(\beta,\;E_1,\;E_2,\;\ldots)=-\frac{\partial\ln Z}{\partial\beta}</math>
или, что то же самое
- <math>\langle E\rangle=k_B T^2\frac{\partial\ln Z}{\partial T}.</math>
Можно также заметить, что если энергии микросостояний зависят от параметра <math>\lambda</math> как
- <math>E_j=E_j^{(0)}+\lambda A_j</math>
для всех <math>j</math>, то среднее значение <math>A</math> равно
- <math>\langle A\rangle=\sum_j A_jP_j=-\frac{1}{\beta}\frac{\partial}{\partial\lambda}\ln Z(\beta,\;\lambda).</math>
На этом основан приём, позволяющий вычислить средние значения многих микроскопических величин. Нужно искусственно добавить эту величину к энергии микросостояний (или, на языке квантовой механики, к гамильтониану), вычислить новую статистическую сумму и среднее значение, а затем в итоговом выражении положить <math>\lambda</math> равным нулю. Аналогичный метод применяется в квантовой теории поля.
Связь с термодинамическими величинами
В этом разделе приведена связь статистической суммы с различными термодинамическими параметрами системы. Эти результаты могут быть получены с помощью метода, описанного в предыдущем разделе, и различных термодинамических соотношений.
Как мы уже видели, энергия равна
- <math>\langle E\rangle=-\frac{\partial\ln Z}{\partial\beta}.</math>
Флуктуация энергии равна
- <math>\langle\delta E^2\rangle\equiv\langle(E-\langle E\rangle)^2\rangle=\frac{\partial^2\ln Z}{\partial\beta^2}.</math>
Теплоёмкость равна
- <math>c_v=\frac{\partial\langle E\rangle}{\partial T}=\frac{1}{k_B T^2}\langle\delta E^2\rangle.</math>
Энтропия равна
- <math>S\equiv-k_B\sum_j P_j\ln P_j=k_B(\ln Z+\beta\langle E\rangle)=\frac{\partial}{\partial T}(k_B T\ln Z)=-\frac{\partial F}{\partial T},</math>
где <math>F</math> — свободная энергия, определяемая как <math>F=E-TS</math>, где <math>E</math> — полная энергия, а <math>S</math> — энтропия, так что
- <math>F=\langle E\rangle-TS=-k_B T\ln Z.</math>
Статистическая сумма подсистем
Предположим, что система состоит из <math>N</math> подсистем, взаимодействие между которыми пренебрежимо мало. Если статистические суммы подсистем равны <math>\zeta_1,\;\zeta_2,\;\ldots,\;\zeta_N</math>, то статистическая сумма всей системы равна произведению отдельных статистических сумм:
- <math>Z =\prod_{j=1}^N\zeta_j.</math>
Если подсистемы обладают одинаковыми физическими свойствами, то их статистические суммы одинаковы: <math>\zeta_1=\zeta_2=\ldots=\zeta</math>, и в этом случае
- <math>Z=\zeta^N.</math>
Из этого правила, однако, есть одно известное исключение. Если подсистемы — это тождественные частицы, то есть, исходя из принципов квантовой механики, их невозможно различить даже в принципе, общая статистическая сумма должна быть разделена на <math>N!</math>:
- <math>Z=\frac{\zeta^N}{N!}.</math>
Это делается, чтобы не учитывать одно и то же микросостояние несколько раз.
Статистическая сумма большого канонического ансамбля
Определение
Аналогично канонической статистической сумме для канонического ансамбля, можно определить большую каноническую статистическую сумму для большого канонического ансамбля — системы, которая может обмениваться со средой и теплотой, и частицами, и имеет постоянную температуру <math>T</math>, объём <math>V</math> и химический потенциал <math>\mu</math>. Большая каноническая статистическая сумма, хотя и более сложна для понимания, упрощает расчёт квантовых систем. Большая каноническая статистическая сумма <math>\mathcal{Z}</math> для квантового идеального газа записывается как:
- <math>\mathcal{Z}=\sum_{N=0}^\infty\,\sum_{\{n_i\}}\,\prod_i e^{-\beta n_i(\varepsilon_i-\mu)},</math>
где <math>N</math> — общее количество частиц в объёме <math>V</math>, индекс <math>i</math> пробегает все микросостояния системы, <math>n_i</math> — число частиц в состоянии <math>i</math>, а <math>\varepsilon_i</math> — энергия в состоянии <math>i</math>. <math>\{n_i\}</math> — всевозможные наборы чисел заполнения каждого микросостояния, такие что <math>\sum_i n_i=N</math>. Рассмотрим, например, слагаемое, соответствующее <math>N=3</math>. Один из возможных наборов чисел заполнения будет <math>\{n_i\}=0,\;1,\;0,\;2,\;0,\ldots</math>, он даёт вклад в слагаемое с <math>N=3</math>, равный
- <math>\prod_i e^{-\beta n_i(\varepsilon_i-\mu)}=e^{-\beta(\varepsilon_1-\mu)}\,e^{-2\beta(\varepsilon_3-\mu)}.</math>
Для бозонов числа заполнения могут принимать любые целые неотрицательные значения при том, что их сумма равна <math>N</math>. Для фермионов, в соответствии с принципом запрета Паули, числа заполнения могут быть равны только 0 или 1, но их сумма опять же равна <math>N</math>.
Частные случаи
Можно показать, что указанное выражение для большой канонической статистической суммы математически эквивалентно следующему:
- <math>\mathcal{Z}=\prod_i\mathcal{Z}_i.</math>
(Это произведение иногда берётся по всем значениям энергии, а не по отдельным состояниям, и в этом случае каждая отдельная статистическая сумма должна быть возведена в степень <math>g_i</math>, где <math>g_i</math> — число состояний с такой энергией. <math>g_i</math> также называется степенью вырождения.)
Для системы, состоящей из бозонов:
- <math>\mathcal{Z}_i=\sum_{n_i=0}^\infty e^{-\beta n_i(\varepsilon_i-\mu)}=\frac{1}{1-e^{-\beta(\varepsilon_i-\mu)}},</math>
а для системы, состоящей из фермионов:
- <math>\mathcal{Z}_i=\sum_{n_i=0}^1 e^{-\beta n_i(\varepsilon_i-\mu)}=1+e^{-\beta(\varepsilon_i-\mu)}.</math>
В случае максвелловско-больцмановского газа необходимо корректно подсчитывать состояния и делить больцмановский множитель <math>e^{-\beta (\varepsilon_i-\mu)}</math> на <math>n_i!</math>
- <math>\mathcal{Z}_i=\sum_{n_i=0}^\infty\frac{e^{-\beta n_i(\varepsilon_i-\mu)}}{n_i!}=\exp\left(e^{-\beta(\varepsilon_i-\mu)}\right).</math>
Связь с термодинамическими величинами
Так же как и каноническая статистическая сумма, большую каноническую статистическую сумму можно использовать для вычисления термодинамических и статистических величин системы. Как и в каноническом ансамбле, термодинамические величины не фиксированы, а статистически распределены вокруг среднего значения. Обозначая <math>\alpha=-\beta\mu</math>, получаем средние значения чисел заполнения:
- <math>\langle n_i\rangle=-\left(\frac{\partial\ln\mathcal{Z}_i}{\partial\alpha}\right)_{\beta,\;V}=\frac{1}{\beta}\left(\frac{\partial\ln\mathcal{Z}_i}{\partial\mu}\right)_{\beta,\;V}.</math>
Для больцмановских частиц это даёт:
- <math>\langle n_i\rangle=e^{-\beta(\varepsilon_i-\mu)}.</math>
Для бозонов:
- <math>\langle n_i\rangle=\frac{1}{e^{\beta(\varepsilon_i-\mu)}-1}.</math>
Для фермионов:
- <math>\langle n_i\rangle=\frac{1}{e^{\beta(\varepsilon_i-\mu)}+1},</math>
что совпадает с результатами, получаемыми с помощью канонического ансамбля для статистики Максвелла — Больцмана, статистики Бозе — Эйнштейна и статистики Ферми — Дирака соответственно. (Степень вырождения <math>g_i</math> отсутствует в этих уравнениях, поскольку индекс <math>i</math> нумерует отдельные состояния, а не уровни энергии.)
Общее число частиц
- <math>\langle N\rangle=-\left(\frac{\partial\ln\mathcal{Z}}{\partial\alpha}\right)_{\beta,\;V}=\frac{1}{\beta}\left(\frac{\partial\ln\mathcal{Z}}{\partial\mu}\right)_{\beta,\;V}.</math>
Флуктуация общего числа частиц
- <math>\mathrm{var}\,(N)=\left(\frac{\partial^2\ln\mathcal{Z}}{\partial\alpha^2}\right)_{\beta,\;V}.</math>
Внутренняя энергия
- <math>\langle E\rangle=-\left(\frac{\partial\ln\mathcal{Z}}{\partial\beta}\right)_{\mu,\;V}+\mu\langle N\rangle.</math>
Флуктуация внутренней энергии
- <math>\mathrm{var}\,(E)=\left(\frac{\partial^2\ln\mathcal{Z}}{\partial\beta^2}\right)_{\mu,\;V}.</math>
- <math>\langle P\rangle=\frac{1}{\beta}\left(\frac{\partial\ln\mathcal{Z}}{\partial V}\right)_{\mu,\;\beta}.</math>
Механическое уравнение состояния
- <math>\langle PV\rangle=\frac{\ln\mathcal{Z}}{\beta}.</math>
Литература
- Кубо Р. Статистическая механика. — М.: Мир, 1967.
- Хуанг К. Статистическая механика. — М.: Мир, 1966. (Huang, Kerson, «Statistical Mechanics», John Wiley & Sons, New York, 1967.)
- Исихара А. Статистическая физика. — М.: Мир, 1973. (Isihara A. «Statistical Physics». — New York: Academic Press, 1971.)
- Kelly, James J. Lecture notes.
- Книга:Ландау Л. Д., Лифшиц Е. М.: Статистическая физика.