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

ВведениеПравить

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

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

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

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

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

Распределение энергии Максвелла может быть выражено как дискретное распределение энергии: N i N = exp ( E i / k T ) j exp ( E j / k T ) ( 1 ) , \frac {N_i} {N} = \frac {\exp\left (-E_i/kT \right)} {\sum _ {j} ^ {} {\exp\left (-E_j/kT\right)}} \qquad\qquad (1) ,

где Ni является числом молекул имеющих энергию Ei при температуре системы T, N является общим числом молекул в системе и kпостоянная Больцмана. (Отметьте, что иногда вышеупомянутое уравнение записывается с множителем gi, обозначающим степень вырождения энергетических уровней. В этом случае сумма будет по всем энергиям, а не всем состояниям системы). Поскольку скорость связана с энергией, уравнение (1) может использоваться для получения связи между температурой и скоростями молекул в газе. Знаменатель в уравнении (1) известен как каноническая статистическая сумма.

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

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

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

В случае идеального газа, состоящего из не взаимодействующих атомов в основном состоянии, вся энергия находится в форме кинетической энергии. Кинетическая энергия соотносится с импульсом частицы следующим образом E = p 2 2 m , E =\frac {p^2} {2m} ,

где p2 — квадрат вектора импульса p = [ p x , p y , p z ] \,\mathbf{p}=[p_x,p_y,p_z] .

Мы можем поэтому переписать уравнение (1) как: N i N = 1 Z exp [ ( p x 2 + p y 2 + p z 2 ) 2 m k T ] ( 3 ) , \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постоянная Больцмана. Это распределение Ni / N пропорционально функции плотности вероятности f p \,f_\mathbf{p} нахождения молекулы в состоянии с этими значениями компонентов импульса. Таким образом: f p ( p x , p y , p z ) = C Z exp [ ( p x 2 + p y 2 + p z 2 ) 2 m k T ] ( 4 ) 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) по всем значениям px, py и pz должен быть равен единице. Можно показать, что: + 1 Z exp [ ( p x 2 + p y 2 + p z 2 ) 2 m k T ] d p x d p y d p z = 1 Z ( 2 π m k T ) 3 / 2 ( 5 ) . \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 m kT \right) ^ {3/2} \qquad\qquad (5) .

Таким образом, чтобы интеграл в уравнении (4) имел значение 1 необходимо, чтобы c = Z ( 2 π m k T ) 3 ( 6 ) . c = \frac {Z} {(\sqrt {2 \pi mkT}) ^ 3} \qquad\qquad (6) .

Подставляя выражение (6) в уравнение (4) и используя тот факт, что pi = m vi, мы получим f p ( p x , p y , p z ) = ( 1 2 π m k T ) 3 exp [ ( p x 2 + p y 2 + p z 2 ) 2 m k T ] ( 7 ) . 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 v \,f_\mathbf{v} пропорциональна плотности распределения по импульсам: f v d 3 v = f p ( d p d v ) 3 d 3 v f_\mathbf {v} d^3v = f_\mathbf {p} \left (\frac {dp} {dv} \right) ^3 d^3v

и используя p = m v \,\mathbf{p}=m\mathbf{v} мы получим: f v ( v x , v y , v z ) = ( m 2 π k T ) 3 exp [ m ( v x 2 + v y 2 + v z 2 ) 2 k T ] ( 8 ) , 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) ,

что является распределением Максвелла по скоростям. Вероятность обнаружения частицы в бесконечно малом элементе dvx, dvy, dvz около скорости v = [ v x , v y , v z ] \,\mathbf{v}=[v_x,v_y,v_z] равна f v ( v x , v y , v z ) d v x d v y d v z f_\mathbf {v} \left (v_x, v_y, v_z\right) dv_x dv_y dv_z

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

Интегрируя, мы можем найти распределение по абсолютной величине импульса f p = θ = 0 π ϕ = 0 2 π   f p p 2 sin  Синус  ( θ ) d θ d ϕ = 4 π ( 1 2 π m k T ) 3   p 2 exp [ p 2 2 m k T ] 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]

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

Наконец, используя соотношения p2 = 2 m E и fE dE = fp dp, мы получаем распределение по кинетической энергии: f E = f p d p d E = 2 E π ( k T ) 3   exp [ E k T ] f_E=f_p \frac{dp} {dE} =2\sqrt {\frac {E} {\pi (kT) ^3}} ~ \exp\left [\frac {-E} {kT} \right]

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

Распределение Максвелла для вектора скорости [vx,vy,vz] — является произведением распределений для каждого из трех направлений: f v ( v x , v y , v z ) = f v ( v x ) f v ( v y ) f v ( 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 ) = m 2 π k T exp [ m v i 2 2 k T ] ( 9 ) f_v (v_i) = \sqrt {\frac {m} {2 \pi kT}} \exp \left [\frac {-mv_i^2} {2kT} \right] \qquad\qquad (9)

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

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

Обычно, более интересно распределение по абсолютному значению, а не по проекциям скоростей молекул. Модуль скорости, v определяется как: v = v x 2 + v y 2 + v z 2 ( 10 ) v = \sqrt {v_x^2 + v_y^2 + v_z^2} \qquad\qquad (10)

поэтому модуль скорости всегда будет больше или равен нулю. Так как все vi распределены нормально, то v2 будет иметь хи-квадрат распределение с тремя степенями свободы. Если f ( v ) \,f(\mathbf{v}) функция плотности вероятности для модуля скорости, то: f ( v ) d v = P ( χ 2 | 3 ) d χ 2 , f\left (v\right) dv = P (\chi^2|3) d\chi^2 , где χ 2 = m v 2 k T \chi^2 = \frac {mv^2} {kT}

таким образом, функция плотности вероятности для модуля скорости равна f ( v ) d v = 4 π v 2 ( m 2 π k T ) 3 / 2 exp ( m v 2 2 k T ) d v ( 11 ) 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) дает распределение скоростей, или, другими словами, долю молекул, имеющих специфическую скорость, часто более интересны другие величины, такие как средние скорости частиц. В следующих подразделах мы определим и получим наиболее вероятную скорость, среднюю скорость и среднеквадратичную скорость.

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

наиболее вероятная скорость, vp — вероятность обладания которой любой молекулой системы максимальна, и которая соответствует максимальному значению f(v). Чтобы найти её, необходимо вычислить df / dv, приравнять её нулю и решить относительно v: d f ( v ) d v = ( m 2 π k T ) 3 / 2 exp ( m v 2 / 2 k T ) [ 8 π v + 4 π v 2 ( m v / k T ) ] = 0 ( 12 ) \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 = 2 k T m = 2 R T μ ( 13 ) v_p = \sqrt {\frac {2kT} {m}} = \sqrt {\frac {2RT} {\mu}} \qquad\qquad (13)

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

Undefined control sequence \infin \langle v \rangle = \int\limits_0 ^ {\infin} v \, f (v) \, dv \qquad\qquad (14)

Подставляя f(v) и интегрируя, мы получим v = 8 k T π m = 8 R T π μ ( 15 ) \langle v \rangle = \sqrt {\frac {8kT} {\pi m}} = \sqrt {\frac {8RT} {\pi \mu}} \qquad\qquad (15)

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

Undefined control sequence \infin \langle v^2 \rangle = \int\limits_0 ^ {\infin} v^2 \, f (v) \, dv \qquad\qquad (16)

Подставляя f(v) и интегрируя, мы получим v 2 = 3 k T m = 3 R T μ ( 17 ) \sqrt {\langle v^2 \rangle} = \sqrt {\frac {3kT} {m}} = \sqrt {\frac {3RT} { \mu}} \qquad\qquad (17)

Страница: 0

en Maxwell–Boltzmann distribution

ПримечанияПравить

См. такжеПравить

СсылкиПравить

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

ЛитератураПравить

Для статьиПравить