Ричард Фейнман - 4. Кинетика. Теплота. Звук
Теперь вернемся к вопросу о том, можно ли пренебрегать столкновениями. Почему это не имеет значения? Мы могли бы повторить все наши доводы, используя не конечную высоту h, а бесконечно малую высоту h, столь малую, что для столкновений между высотами 0 и h было бы слишком мало места. Но это не обязательно: наши доводы, очевидно, основаны лишь на анализе значений энергий и на сохранении энергии; при столкновениях же происходит обмен энергиями среди молекул. Но нам довольно безразлично, следим ли мы за одной и той же молекулой, раз происходит лишь обмен энергиями с другой молекулой. И получается, что если мы даже сделаем это достаточно тщательно (а такую работу тщательно проделать, конечно, труднее), то результат будет тот же.
Интересно, что найденное нами распределение по скоростям имеет вид
n>u~e-к.э./kT. (40.4)
Этот способ описания распределения по скоростям —когда подсчитывается число молекул, проходящих через выделенную площадку с заданной минимальной z-составляющей скорости,— отнюдь не самый удобный. Например, чаще хотят знать, сколько молекул в заданном объеме газа движется, имея z-составляющую скорости между двумя заданными значениями, а это, конечно, из (40.4) сразу не получишь. Поэтому придадим нашей формуле удобную форму, хотя то, что мы получили, — это весьма общий результат. Заметим, что невозможно утверждать, что любая молекула в точности обладает той или иной наперед заданной скоростью; ни одна из них не движется со скоростью, в точности равной 1,7962899173 м/сек. Итак, чтобы придать нашему утверждению какой-то смысл, мы должны спросить, сколько молекул можно найти в заданном интервале скоростей. Нам придется говорить о том, как часто встречаются скорости в интервале между 1,796 и 1,797 и т. п. Выражаясь математически, пусть f(u)du будет долей всех молекул, чьи скорости заключены в промежутке u и u+du, или, что то же самое (если du бесконечно мало), долей всех молекул, имеющих скорость и с точностью до du. На фиг. 40.5 представлена возможная форма функции f(u), а заштрихованная часть ширины du и средней высоты f(u) — это доля молекул f(u)du. Таким образом, отношение площади заштрихованного участка ко всей площади под кривой равно относительному числу молекул со скоростью и внутри отрезка du.
Фиг. 40.5. Функция, распределения скоростей.
Заштрихованная площадь равна f(u)du — это относительное число частиц, скорости которых заключены внутри отрезка du около точки u.
Если определить f(u) так, что относительное число молекул будет просто равно площади заштрихованного участка, то полная площадь под кривой — это все 100% молекул, т. е.
Теперь остается только найти это распределение, сравнив его с результатом доказанной ранее теоремы. Сначала надо выяснить, как выразить через f(u) число молекул, проходящих за 1 сек через заданную площадку со скоростью, превышающей u?
Это число не равно интегралу (хотя это первое, что приходит в голову), ведь нас интересует число молекул, проходящих через площадку за секунду. Более быстрые молекулы будут пересекать площадку, так сказать, чаще, чем более медленные, поэтому, чтобы найти число проходящих молекул, надо умножить плотность молекул на скорость. (Мы уже обсуждали это в предыдущей главе, когда подсчитывали число столкновений.)
Полное число молекул, проходящих через поверхность за время t, равно числу молекул, способных достигнуть поверхности, а это молекулы, проходящие к поверхности с расстояния ut. Таким образом, число молекул, достигающих площадки, определяется не просто числом молекул, движущихся с данной скоростью, а равно этому числу, отнесенному к единице объема, и умноженному на расстояние, которое они пройдут, прежде чем достигнут площадки, сквозь которую они, по-видимому, должны пройти, а это расстояние пропорционально и. Значит, нам предстоит вычислить интеграл от произведения и на f(u)du, взятый от и до бесконечности, причем мы уже знаем, что этот интеграл обязательно должен быть пропорционален ехр(-mu2/2kT), а постоянную пропорциональности еще надо определить:
Если теперь продифференцировать интеграл по и, то мы получим подынтегральное выражение (со знаком минус, потому что и — это нижний предел интегрирования), а дифференцируя правую часть равенства, мы получим произведение и на экспоненту (и на некоторую постоянную). Сократим в обеих частях и, и тогда
Мы оставили в обеих частях равенства du, чтобы помнить, что это распределение; оно говорит нам об относительном числе молекул, имеющих скорость между u и u+du.
Постоянная С должна определиться из условия равенства интеграла единице в согласии с уравнением (40.5). Можно доказать, что
Используя это обстоятельство, легко найти С=Ц(m/2pkT).
Поскольку скорость и импульс пропорциональны, можно утверждать, что распределение молекул по импульсам, отнесенное к единице импульсной шкалы, также пропорционально ехр(-к.э./kT). Оказывается, что эта теорема верна также в теории относительности, если только формулировать ее в терминах импульсов, тогда как в терминах скоростей это уже не так; поэтому сформулируем все в терминах импульсов:
f(p)dp=ce-к.э./kTdp. (40.8)
Это значит, что мы установили, что вероятности, определяемые энергиями разного происхождения (и кинетической и потенциальной), в обоих случаях выражаются одинаково: ехр(-энергия/kT); таким образом, наша замечательная теорема приобрела форму, весьма удобную для запоминания.
Однако пока мы говорили только о «вертикальном» распределении скоростей. Но мы можем спросить, какова вероятность того, что молекула движется в другую сторону? Конечно, эти распределения связаны друг с другом и можно получить полное распределение, исходя из какого-то одного, ведь полное распределение зависит только от квадрата величины скорости, а не от ее z-составляющей. Распределение по скоростям не должно зависеть от направления и определяться только функцией u2 — вероятностью величины скорости. Нам известно распределение z-составляющей, и мы хотим получить отсюда распределение других составляющих. В результате полное распределение по-прежнему пропорционально ехр(-к.э./kT), только теперь кинетическая энергия состоит из трех частей: mv2x/2, mv2y/2 и mv2z/2, суммируемых в показателе экспоненты. А можно записать это и в виде произведения:
f(vx,, vy, vz) dvx dvy dvz~
Вы можете убедиться в том, что эта формула верна, ибо, во-первых, распределение зависит только от v2 и, во-вторых, вероятности данных vгполучаются после интегрирования по всем vxи vyи это должно привести к (40.7). Но обоим этим требованиям удовлетворяет только функция (40.9).
§ 5. Удельные теплоемкости газов
Посмотрим теперь, как можно проверить теорию и оценить, насколько хороша классическая теория газов. Мы уже говорили, что если U—внутренняя энергия N молекул, то формула pV=NkT=(g-1)U иногда и для некоторых газов может оказаться правильной. Мы знаем, что для одноатомного газа правая часть равна 2/3 кинетической энергии движения центров масс атомов. В случае одноатомного газа кинетическая энергия равна внутренней энергии, поэтому g-1 == 2/з.
Но предположим, что мы столкнулись с более сложной молекулой, которая может вращаться и колебаться, и предположим (в классической механике это так), что энергии внутренних движений также пропорциональны kT. Поэтому при заданной температуре молекула, кроме кинетической энергии kT, имеет внутреннюю энергию колебания и вращения. Тогда полная энергия U включает не только кинетическую энергию, но и вращательную энергию и мы получаем другие значения у. Наилучший способ измерения g — это измерение удельной теплоемкости, характеризующей изменение энергии при изменении температуры. К этому способу мы еще вернемся, а пока предположим, что нам удалось экспериментально определить g с помощью кривой PVg , соответствующей адиабатическому сжатию.