Распределение Максвелла

Функция плотности распределения

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

Распределение Максвелла может быть получено при помощи статистической механики (см. происхождение статсуммы). Как распределение энергии, оно соответствует самому вероятному распределению энергии, в столкновительно-доминируемой системе, состоящей из большого количества невзаимодействующих частиц, в которой квантовые эффекты являются незначительными. Так как взаимодействие между молекулами в газе является обычно весьма небольшим, распределение Максвелла даёт довольно хорошее приближение ситуации, существующей в газе.

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

Распределение энергии Максвелла может быть выражено как дискретное распределение энергии:

${\frac {N_{i}}{N}}={\frac {\exp \left(-E_{i}/kT\right)}{\sum _{j}^{}{\exp \left(-E_{j}/kT\right)}}}\qquad \qquad (1)$,

где $\,N_{i}$ является числом молекул имеющих энергию $\,E_{i}$ при температуре системы $\,T$, $\,N$ является общим числом молекул в системе и $\,k$ — постоянная Больцмана. (Отметьте, что иногда вышеупомянутое уравнение записывается с множителем $\,g_{i}$, обозначающим степень вырождения энергетических уровней. В этом случае сумма будет по всем энергиям, а не всем состояниям системы). Поскольку скорость связана с энергией, уравнение (1) может использоваться для получения связи между температурой и скоростями молекул в газе. Знаменатель в уравнении (1) известен как каноническая статистическая сумма.

Распределение Максвелла

Распределение по вектору импульса

Представленное ниже очень сильно отличается от вывода, предложенного Джеймсом Клерком Максвеллом и позже описанного с меньшим количеством предположений Людвигом Больцманом.

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

$E={\frac {p^{2}}{2m}}$,

где $\,p^{2}$ — квадрат вектора импульса $\,\mathbf {p} =[p_{x},p_{y},p_{z}]$.

Мы можем поэтому переписать уравнение (1) как:

${\frac {N_{i}}{N}}={\frac {1}{Z}}\exp \left[{\frac {-(p_{x}^{2}+p_{y}^{2}+p_{z}^{2})}{2mkT}}\right]\qquad \qquad (3)$,

где $\,Z$ — статсумма, соответствующая знаменателю в уравнении (1), $\,m$ — молекулярная масса газа, $\,T$ — термодинамическая температура, и $\,k$ — постоянная Больцмана. Это распределение $\,N_{i}/N$ пропорционально функции плотности вероятности $\,f_{\mathbf {p} }$ нахождения молекулы в состоянии с этими значениями компонентов импульса. Таким образом:

$f_{\mathbf {p} }(p_{x},p_{y},p_{z})={\frac {C}{Z}}\exp \left[{\frac {-(p_{x}^{2}+p_{y}^{2}+p_{z}^{2})}{2mkT}}\right]\qquad \qquad (4)$

Постоянная нормировки C, определяется из условия, в соответствии с которым вероятность того, что молекулы имеют какой-либо вообще импульс, должна быть равна единице. Поэтому интеграл уравнения (4) по всем значениям $\,p_{x}\,,p_{y}$ и $\,p_{z}$ должен быть равен единице. Можно показать, что:

$\iiint \limits _{-\infty }^{+\infty }{\frac {1}{Z}}\exp \left[{\frac {-(p_{x}^{2}+p_{y}^{2}+p_{z}^{2})}{2mkT}}\right]\,dp_{x}\,dp_{y}\,dp_{z}={\frac {1}{Z}}\left(2\pi mkT\right)^{3/2}\qquad \qquad (5)$.

Таким образом, чтобы интеграл в уравнении (4) имел значение 1 необходимо, чтобы

$c={\frac {Z}{({\sqrt {2\pi mkT}})^{3}}}\qquad \qquad (6)$.

Подставляя выражение (6) в уравнение (4) и используя тот факт, что $\,p_{i}=mv_{i}$, мы получим

$f_{\mathbf {p} }(p_{x},p_{y},p_{z})={\sqrt {\left({\frac {1}{2\pi mkT}}\right)^{3}}}\exp \left[{\frac {-(p_{x}^{2}+p_{y}^{2}+p_{z}^{2})}{2mkT}}\right]\qquad \qquad (7)$.

Распределение по вектору скорости

Учитывая, что плотность распределения по скоростям $\,f_{\mathbf {v} }$ пропорциональна плотности распределения по импульсам:

$f_{\mathbf {v} }d^{3}v=f_{\mathbf {p} }\left({\frac {dp}{dv}}\right)^{3}d^{3}v$

и используя $\,\mathbf {p} =m\mathbf {v}$ мы получим:

$f_{\mathbf {v} }(v_{x},v_{y},v_{z})={\sqrt {\left({\frac {m}{2\pi kT}}\right)^{3}}}\exp \left[{\frac {-m(v_{x}^{2}+v_{y}^{2}+v_{z}^{2})}{2kT}}\right]\qquad \qquad (8)$,

что является распределением Максвелла по скоростям. Вероятность обнаружения частицы в бесконечно малом элементе $\,dv_{x}\,,dv_{y}\,,dv_{z}$ около скорости $\,\mathbf {v} =[v_{x},v_{y},v_{z}]$ равна

$f_{\mathbf {v} }\left(v_{x},v_{y},v_{z}\right)dv_{x}dv_{y}dv_{z}$

Распределение по абсолютной величине импульса

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

$f_{p}=\int _{\theta =0}^{\pi }\int _{\phi =0}^{2\pi }~f_{\mathbf {p} }p^{2}\sin(\theta )\,d\theta \,d\phi =4\pi {\sqrt {\left({\frac {1}{2\pi mkT}}\right)^{3}}}~p^{2}\exp \left[{\frac {-p^{2}}{2mkT}}\right]$

Распределение по энергии

Наконец, используя соотношения $\,p^{2}=2\,mE$ и $\,f_{E}\,dE=f_{p}\,dp$, мы получаем распределение по кинетической энергии:

$f_{E}=f_{p}{\frac {dp}{dE}}=2{\sqrt {\frac {E}{\pi (kT)^{3}}}}~\exp \left[{\frac {-E}{kT}}\right]$

Распределение по проекции скорости

Распределение Максвелла для вектора скорости $\,[v_{x},v_{y},v_{z}]$ — является произведением распределений для каждого из трех направлений:

$f_{v}\left(v_{x},v_{y},v_{z}\right)=f_{v}(v_{x})f_{v}(v_{y})f_{v}(v_{z})$,

где распределение по одному направлению:

$f_{v}(v_{i})={\sqrt {\frac {m}{2\pi kT}}}\exp \left[{\frac {-mv_{i}^{2}}{2kT}}\right]\qquad \qquad (9)$

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

Распределение по модулю скоростей

Обычно, более интересно распределение по абсолютному значению, а не по проекциям скоростей молекул. Модуль скорости, v определяется как:

$v={\sqrt {v_{x}^{2}+v_{y}^{2}+v_{z}^{2}}}\qquad \qquad (10)$

поэтому модуль скорости всегда будет больше или равен нулю. Так как все $\,v_{i}$ распределены нормально, то $\,v^{2}$ будет иметь хи-квадрат распределение с тремя степенями свободы. Если $\,f(\mathbf {v} )$ — функция плотности вероятности для модуля скорости, то:

$f\left(v\right)dv=P(\chi ^{2}|3)d\chi ^{2}$,

где

$\chi ^{2}={\frac {mv^{2}}{kT}}$

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

$f(v)dv=4\pi v^{2}\left({\frac {m}{2\pi kT}}\right)^{3/2}\exp \left({\frac {-mv^{2}}{2kT}}\right)dv\qquad \qquad (11)$

Характерная скорость

Хотя Уравнение (11) дает распределение скоростей, или, другими словами, долю молекул, имеющих специфическую скорость, часто более интересны другие величины, такие как средние скорости частиц. В следующих подразделах мы определим и получим наиболее вероятную скорость, среднюю скорость и среднеквадратичную скорость.

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

наиболее вероятная скорость, $\,v_{p}$ — вероятность обладания которой любой молекулой системы максимальна, и которая соответствует максимальному значению Невозможно разобрать выражение (SVG (MathML можно включить с помощью плагина для браузера): Недопустимый ответ («Math extension cannot connect to Restbase.») от сервера «http://localhost:6011/ru.wikipedia.org/v1/»:): {\displaystyle \,f(v)} . Чтобы найти её, необходимо вычислить $\,df/dv$, приравнять её нулю и решить относительно $\,v$:

${\frac {df(v)}{dv}}=\left({\frac {m}{2\pi kT}}\right)^{3/2}\exp \left(-mv^{2}/2kT\right)\left[8\pi v+4\pi v^{2}(-mv/kT)\right]=0\qquad \qquad (12)$

$v_{p}={\sqrt {\frac {2kT}{m}}}={\sqrt {\frac {2RT}{\mu }}}\qquad \qquad (13)$

Средняя скорость

$\langle v\rangle =\int \limits _{0}^{\infty }v\,f(v)\,dv\qquad \qquad (14)$

Подставляя Невозможно разобрать выражение (Ошибка преобразования. Сервер («https://wikimedia.org/ruwiki/api/rest_») сообщил: «Cannot get mml. upstream connect error or disconnect/reset before headers. reset reason: connection termination»): {\displaystyle \,f(v)} и интегрируя, мы получим

$\langle v\rangle ={\sqrt {\frac {8kT}{\pi m}}}={\sqrt {\frac {8RT}{\pi \mu }}}\qquad \qquad (15)$

Среднеквадратичная скорость

$\langle v^{2}\rangle =\int \limits _{0}^{\infty }v^{2}\,f(v)\,dv\qquad \qquad (16)$

Подставляя Невозможно разобрать выражение (Ошибка преобразования. Сервер («https://wikimedia.org/ruwiki/api/rest_») сообщил: «Cannot get mml. upstream connect error or disconnect/reset before headers. reset reason: connection termination»): {\displaystyle \,f(v)} и интегрируя, мы получим

${\sqrt {\langle v^{2}\rangle }}={\sqrt {\frac {3kT}{m}}}={\sqrt {\frac {3RT}{\mu }}}\qquad \qquad (17)$

Вывод распределения по Максвеллу

Получим теперь формулу распределения так, как это делал сам Джеймс Клерк Максвелл.
Рассмотрим пространство скоростных точек (каждую молекулу представляем как точку в системе координат $Ov_{x}v_{y}v_{z}\!$) в стационарном состоянии газа. Выберем бесконечно малый элемент объема $dv_{x}dv_{y}dv_{z}\!$. Так как газ стационарный, количество скоростных точек в $dv_{x}dv_{y}dv_{z}\!$ остается неизменным с течением времени. Пространство скоростей изотропно, поэтому функции плотности вероятности для всех направлений одинаковы.

$dP(v_{x})=\varphi (v_{x})dv_{x}\qquad dP(v_{y})=\varphi (v_{y})dv_{y}\qquad dP(v_{z})=\varphi (v_{z})dv_{z}\!$

Максвелл предположил, что распределения скоростей по направлениям статистически независимы, то есть компонента $v_{x}\!$ скорости молекулы не зависит от $y-\!$ и $z-\!$ компонент.

$dP(v_{x},v_{y},v_{z})=\underbrace {\varphi (v_{x})\varphi (v_{y})\varphi (v_{z})} _{f(v)}dv_{x}dv_{y}dv_{z}\!$ - фактически вероятность нахождения скоростной точки в объеме $dv_{x}dv_{y}dv_{z}\!$.

$f(v)=\varphi (v_{x})\varphi (v_{y})\varphi (v_{z})\!$

$\ln f(v)=\ln \varphi (v_{x})+\ln \varphi (v_{y})+\ln \varphi (v_{z})\quad {\bigg |}\quad {\frac {\partial }{\partial v_{x}}}\!$

${\frac {f'(v)}{f(v)}}{\frac {\delta v}{\delta v_{x}}}={\frac {\varphi '(v_{x})}{\varphi (v_{x})}}\!$

${\frac {1}{v}}{\frac {f'(v)}{f(v)}}={\frac {1}{v_{x}}}{\frac {\varphi '(v_{x})}{\varphi (v_{x})}}\!$

Правая часть не зависит от $v_{y}\!$ и $v_{z}\!$, значит и левая от $v_{y}\!$ и $v_{z}\!$ не зависит. Однако $v_{x}\!$ и $v_{y}\!$ равноправны, следовательно левая часть не зависит также и от $v_{x}\!$. Значит данное выражение может лишь равняться некоторой константе.

${\frac {1}{v}}{\frac {f'(v)}{f(v)}}=-\alpha \!$

${\frac {\varphi '(v_{x})}{\varphi (v_{x})}}=-\alpha v_{x}\!$

$\varphi (v_{x})=Ae^{-{\frac {\alpha {v_{x}}^{2}}{2}}}\!$

$\int \limits _{-\infty }^{\infty }\varphi (v_{x})\,dv_{x}=1\qquad \Rightarrow \qquad A\int \limits _{-\infty }^{\infty }e^{-{\frac {\alpha {v_{x}}^{2}}{2}}}\,dv_{x}=A\,{\sqrt {\frac {2}{\alpha }}}\,{\sqrt {\pi }}=1\qquad \Rightarrow \qquad A={\sqrt {\frac {\alpha }{2\pi }}}\!$

$\varphi (v_{x})={\sqrt {\frac {\alpha }{2\pi }}}e^{-{\frac {\alpha {v_{x}}^{2}}{2}}}\!$

Теперь нужно сделать принципиальный шаг — ввести температуру. Кинетическое определение температуры (как меры средней кинетической энергии движения молекул):

$<{\frac {mv^{2}}{2}}>={\frac {3}{2}}kT,\!$

где $k=1.38\cdot 10^{-23}\!$ Дж/К - постоянная Больцмана.

$v={\sqrt {v_{x}^{2}+v_{y}^{2}+v_{z}^{2}}}\!$

Ввиду равноправия всех направлений:

$<v_{x}^{2}>=<v_{y}^{2}>=<v_{z}^{2}>={\frac {1}{3}}<v^{2}>={\frac {kT}{m}}\!$

Чтобы найти среднее значение $v_{x}^{2}\!$, проинтегрируем её вместе с функцией плотности вероятности от минус до плюс бесконечности:

${\frac {kT}{m}}=\int \limits _{-\infty }^{\infty }{v_{x}}^{2}\,\varphi (v_{x})\,dv_{x}={\sqrt {\frac {\alpha }{2\pi }}}\,\int \limits _{-\infty }^{\infty }v_{x}^{2}\,e^{-{\frac {\alpha {v_{x}}^{2}}{2}}}\,dv_{x}={\sqrt {\frac {\alpha }{2\pi }}}\left[-2\,{\frac {\delta }{\delta \alpha }}\,\int \limits _{-\infty }^{\infty }e^{-{\frac {\alpha {v_{x}}^{2}}{2}}}\,dv_{x}\right]=-2\,{\sqrt {\frac {\alpha }{2\pi }}}\,{\frac {\delta }{\delta \alpha }}{\sqrt {\frac {2\pi }{\alpha }}}=-2{\sqrt {\alpha }}(-{\frac {1}{2}}\alpha ^{-{\frac {3}{2}}})={\frac {1}{\alpha }}\!$

Отсюда найдём $\alpha \!$:

$\alpha ={\frac {m}{kT}}\!$

Функция распределения плотности вероятности для $v_{x}\!$ (для $v_{y}\!$ и $v_{z}\!$ аналогично):

$\varphi (v_{x})=\left({\frac {m}{2\pi kT}}\right)^{\frac {1}{2}}\,e^{-{\frac {mv_{x}^{2}}{2kT}}}\!$

Теперь рассмотрим распределение по величине скорости. Вернемся в пространство скоростных точек. Все точки с модулем скорости $v\subset [v;v+dv]\!$ лежат в шаровом слое радиуса $v\!$ и толщины $dv\!$, и $dv_{x}dv_{y}dv_{z}\!$ - объем этого шарового слоя.

$dP(v_{x},v_{y},v_{z})=\varphi (v_{x})\varphi (v_{y})\varphi (v_{z})dv_{x}dv_{y}dv_{z}\!$

$\underbrace {dP(v_{x},v_{y},v_{z})} _{dP(v)}=\left({\frac {m}{2\pi kT}}\right)^{\frac {3}{2}}\,e^{-{\frac {mv_{x}^{2}+mv_{y}^{2}+mv_{z}^{2}}{2kT}}}\underbrace {dv_{x}dv_{y}dv_{z}} _{4\pi v^{2}dv}\!$

$dP(v)=\underbrace {4\pi \,\left({\frac {m}{2\pi kT}}\right)^{\frac {3}{2}}\,v^{2}\,e^{-{\frac {mv^{2}}{2kT}}}} _{F(v)}dv\!$

Таким образом, мы получили функцию плотности вероятности $F(v)\!$, которая и является распределением Максвелла.

Границы применимости

Условия применимости распределения Максвелла:

1. Равновесное состояние системы, состоящей из большого числа частиц.

2. Изотропная система.

3. Классическая система. Это значит, что система должна быть не релятивистской и не квантовой (взаимодействие частиц допускается, но только зависящее от относительного положения частиц).

Условия классического рассмотрения

Рассматриваем объем xyz в газе, на который в среднем приходится 1 частица. Чтобы неопределенности в координате и импульсе не играли роли и применялась бы классическая, а не квантовая механика, должны выполняться соотношения:

$xp_{x}\gg h\qquad yp_{y}\gg h\qquad zp_{z}\gg h,\qquad \!$ где $h\!$ - постоянная Планка.

$Vp^{3}\gg h^{3};\qquad V={\frac {1}{n}}\qquad \!$ - объем, приходящийся на частицу - это полный (единичный) объем, поделенный на количество частиц.

$n\,{\left({\frac {h}{p}}\right)}^{3}\ll 1\!$

$n^{\frac {1}{3}}\,{\frac {h}{m\,{\sqrt {\frac {3kT}{m}}}}}\ll 1\!$

${\frac {n^{\frac {2}{3}}h^{2}}{3mkT}}\ll 1\!$

$T\gg {\frac {n^{\frac {2}{3}}h^{2}}{3mk}}=T_{deg}\qquad \!$ - температура вырождения.

При температурах ниже $T_{deg}\!$ газ становится вырожденным, и распределение Максвелла к нему применять нельзя.

См. также

• Статистика Бозе — Эйнштейна
• Статистика Ферми — Дирака
• Статистика Максвелла — Больцмана
• Опыт Ламмерта
• Опыт Штерна
• Распределение Гиббса

Ссылки

http://www.falstad.com/gas/