beskin
.pdf201
Еще одним важнейшим элементом, лежащим в основе вывода закона Релея-Джинса, стала вновь идея о равнораспределении энергии по степеням свободы. Напомним, что как было показано в разделе 3.2.2, количество различных волн в малом объеме δV = δxδyδz
в диапазоне волновых векторов между kx и kx + δkx, ky и ky + δky , kz и kz + δkz есть (71)31
δN = 2 |
δV δkxδky δkz |
. |
(298) |
|
|||
|
(2π)3 |
|
Иными словами, число δN и определяет число степеней свободы для электромагнитного поля в элементе объема δV и в диапазоне волновых векторов δkxδky δkz . Согласно же
принципу равнораспределения, все они являются равновероятными.
Вспомним теперь, что мы обсуждаем излучение ”абсолютно черного тела”, для которого мы в качестве модели выбрали излучение в замкнутой полости. Понятно, что для такого излучения естественно предположить, что оно будет изотропным, т.е. что поток излучения через небольшую поверхность не будет зависеть от ориентации этой поверхности. В этом случае будет иметь смысл лишь зависимость излучения от величины волнового вектора k, а не от его направления. Поэтому элемент ”объема” в трехмерном пространстве волновых векторов δkxδky δkz удобно, как и в нашем обычном пространстве, заменить на 4πk2δk. Вспоминая теперь, что частота ω и волновой вектор k для фотонов связаны соотношением ω = ck, получаем окончательно для числа волн в объеме δV в
интервале частот от ω до ω + δω |
|
|
|
|
δN = |
ω2 |
δV δω. |
(299) |
|
π2c3 |
||||
|
|
|
Предположим теперь, как это сделали Релей и Джинс, что каждому колебанию можно приписать энергию
ε = T |
(300) |
(напоминаем, что мы везде измеряем температуру в энергетических единицах). Помножив теперь плотность вероятности δN (299) на T , и поделив на объем δV , мы получим
31Здесь дополнительный фактор 2 появился за счет того, что электромагнитные волны (в отличие, например, от звуковых) описываются векторными величинами напряженностями полей. В результате, электромагнитная волна имеет две независимые поляризации.
202
энергию δǫ, заключенную в единице объема в диапазоне частот между ω и ω + δω. Определяя теперь эту энергию как δǫ = w(ω, T )δω, получаем окончательно для плотности энергии w(ω, T ), которая приходится на интервал частот δω:
w(ω, T ) = |
ω2T |
. |
(301) |
|
π2c3 |
||||
|
|
|
Зная же плотность излучения w(ω, T ) (размерность эрг см−3 Гц−1), нетрудно опреде-
лить и спектральную плотность потока P (размерность эрг см−2 c−1Гц−1). Если бы все
фотоны распространялись перпендикулярно элементу площади δS, то тогда можно было бы записать просто P = cw. Действительно, через площадку δS за время δt прошла бы вся энергия, заключенная в объеме δScδt (фотоны распространяются со скоростью света).
В случае же изотропного излучения мы должны учесть то, что фотоны, распространяющийся под углом θ к поверхности δS, будут фактически приходить из объема δScδt cos θ
и, кроме того, половина фотонов вообще будет распространяться в противоположную сторону. Поэтому в этом случае (проверьте!) P = (1/4)cw. В итоге, получаем
|
ω2 |
(302) |
P(ω, T ) = |
4π2c2 T. |
Как уже отмечалось, это соотношение с хорошей точностью описывает излучение абсолютно черного тела при hω¯ T .
Формула (302) является одной из основных расчетных формул в наблюдательной радиоастрономии. Действительно, если из наблюдений известна величина P(ω), то тогда
мы можем определить величину
|
4π2c2 |
(303) |
Tb = |
ω2 P(ω), |
которая называется яркостной температурой источника на частоте ω. Понятно, что величина Tb будет совпадать с реальной температурой излучающей среды лишь в том случае,
если излучение имеет тепловую природу и если частота наблюдения достаточно мала, так что мы можем пользоваться приближением Релея-Джинса. Для источников же, у которых Tb составляет 1015–1020 K, можно с уверенностью сказать, что их излучение не
связано с тепловыми процессами.
203
8.3Формула Планка
Для вывода формулы Планка (20) нам потребуется результат, основанный на полной версии принципа равнораспределения энергии по степеням свободы. А именно распределение Гиббса одно из ключевых соотношений статистической физики. Как мы уже отмечали выше, в Приложении 8.1 был рассмотрен случай, когда все состояния имеют одинаковую энергию. Если же различные состояния имеют разную энергию, то уже интуитивно ясно, что частицам будет выгодно занимать нижние энергетические уровни. Тем не менее, принцип равнораспределения может быть сформулирован и этом случае.
Дело в том, что физически правильной постановкой задачи будет вопрос о том, каким будет распредение N частиц по уровням энергии εn при заданной величине полной
энергии
K |
|
X |
(304) |
E = Nnεn. |
n=1
Здесь Nn число частиц, занимающих энергетический уровень εn, а K количество
различных уровней энергии. При этом заданным должно быть еще и полное число частиц
|
K |
|
|
|
X |
|
(305) |
N = Nn. |
|
||
|
n=1 |
|
|
В результате, ответ на посталенный вопрос дает формула |
|
||
G = |
N ! |
, |
(306) |
|
|||
N1!N2! . . . NK ! |
которая обобщает соотношение (293). При этом величина G (т.е. статистический вес) показывает, сколькими способами можно разместить N1 частиц в первой ячейке, N2 частиц
во второй, и т.д. при выполнении условий (304) и (305). Здесь и используется фундаментальное предположение о том, что все состояния с данной энергией равновероятны.
На Рис. 32 для иллюстрации показано, как можно разместить шесть одинаковых молекул в трех ячейках, имеющих энергии 0, ε и 2ε, так, чтобы полная энергия системы составляла 4ε. Таких комбинаций, как мы видим, всего три. Но количество вариантов,
которыми можно расположить молекулы (т.е. статистические веса соответствующих состояний) при этом различны.
204
.
Рис. 32: Различные варианты размещения шести одинаковых молекул в трех ячейках, имеющих энергии 0, ε и 2ε, так, чтобы полная энергия системы составляла 4ε. Показаны также статистические веса G и вероятности выпадения соответстенных состояний
Используя теперь, как и в Приложении 8.1, формулу Муавра-Стирлинга (295), а также (что очень важно!) соотношения (304) и (305), можно показать, что
εn |
, |
|
Nn = N0 exp − T |
(307) |
где N0 некоторая постоянная. Это и есть соотношение, найденное Дж. Гиббсом (1839-
1903) и носящее его имя. Как и в случае распределения Гаусса (8), оно получается в пределе очень больших значений N . Действительно, при малых N , как показано на Рис. 32,
вероятности различных состояний хотя и отличаются друг от друга, но все же могут совпадать по порядку величины. При больших же N , как и для распределения Гаусса,
практически реализуется лишь наиболее вероятное состояние. Впрочем, и на нашем, чрезвычайно упрощенном примере видно, что наиболее вероятным является состояние, при котором числа заполнения энергетических уровней плавно уменьшаются с увеличением их энергий.
Здесь, правда, неообходимо отметить еще одно обстоятельство. На первый взгляд, непонятно, в каком месте в наших рассуждениях могла возникнуть температура T . Дело
205
втом, что на самом деле при анализе соотношений (304)–(306) величина T возникает
лишь как некоторая константа, природа которой может быть понята лишь при рассмотрении различных приложений. Простейшим примером здесь может служить распределение концентрации изотермической атмосферы n(h) в однородном гравитационном поле
взависимости от высоты h.
Действительно, рассмотрим силы, действующие на элемент объема δSδh, где δS элемент площади, а δh высота. Вниз будет действовать сила гравитации M g = (nmpδSδh)g, где mp масса частиц газа, а вверх сила давления газа, связанная с разностью давлений на верхнюю и нижнюю границы объема [P (h + δh) −P (h)]δS, где P = nT давление газа. Поскольку для изотермической атмосферы P (h + δh) − P (h) = (δn/δh)T δh, имеем
после сокращения на δSδh |
|
|
|
|
|
|
|||
|
dn |
= |
mpg |
n. |
|
|
(308) |
||
|
|
|
|
|
|||||
|
dh |
T |
|
|
|
||||
Решение этого дифференциального уравнения есть |
! . |
|
|||||||
n(h) = n0 exp |
− |
T |
(309) |
||||||
|
|
|
|
|
|
mpgh |
|
|
Как мы видим, форма уранения (308) совершенно аналогична распределению Гиббса, причем в числителе стоит как раз потенциальная энергия mpgh, а в знаменателе температура T .
Теперь, наконец, мы можем вернуться к идее Планка. Как мы помним, она основывалась на предположении, согласно которому энергия осциллятора может принимать лишь значения, кратные кванту энергии εq, где εq = hω¯ . В этом случае величины Nn (307)
образуют геометрическую прогрессию, сумма которой можно легко определить. Напомним, что сам Планк говорил лишь о квантовании энергии осцилляторов, составляющих вещество стенок полости, в котором находится излучение; при этом, конечно же, предполагалось, что излучение и вещество находятся в термодинамическом равновесии32. Нам
же будет удобно воспользоваться тем, что квантование можно ввести и для фотонов. С другой стороны, мы, следуя Планку, рассмотрим упрощенную ситуацию, когда энергия
32Излучение в комнате, в которой Вы находитесь, этому свойству не удовлетворяет!
206
K фотонов может быть записана в виде Khω¯ . Как легко проверить, точная формула (94), содержащая фактор (K + 1/2), приводит к тем же результатам.
Итак, нам вновь нужно посчитать энергию, которую содержат фотоны в диапазоне частот от ω до ω + δω. При этом, как и в пределе Релея-Джинса, число состояний в элементе объема δV будет задаваться формулой δN = (ω2/π2c3)δV δω (299). Отличие будет состоять лишь в том, что энергия разных колебаний будет отличаться от T .
Действительно, давайте применим общий принцип, найденный Гиббсом, для фотонов, имеющих одинаковую энергию hω¯ . Поскольку при эквидистантности уровней осцилляторов энергии излученных фотонов могут принимать лишь значения Khω¯ , где K целое
число, получаем с учетом общего соотношения (307) для вероятности фотону иметь энергию Khω¯
PK = P0 exp − |
T |
! |
, |
|
(310) |
||
|
|
Khω¯ |
|
|
|
|
|
где P0 некоторая константа. При этом сумма геометрической прогрессии PK∞=0 PK |
|||||||
легко вычисляется, и мы получаем |
|
|
|
|
|
|
|
∞ |
"1 − exp |
hω¯ |
!# |
−1 |
(311) |
||
PK = P0 |
− T |
. |
|||||
X |
|
|
|
|
|
|
|
K=0
Поскольку же сумма всех вероятностей должна быть равна единице, получаем оконча-
тельно |
− T |
!# exp − |
T |
! . |
(312) |
|||
PK = "1 − exp |
||||||||
|
|
|
|
hω¯ |
|
Khω¯ |
|
|
Для дальнейшего нам будет удобно записать это выражение в виде |
|
|||||||
1 |
|
|
|
|
|
|
(313) |
|
PK = |
|
exp (−βKhω¯ ) . |
|
|||||
Z(β) |
|
|||||||
Здесь по определению β = 1/T , а |
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
||
|
|
X |
|
|
|
(314) |
||
Z(β) = |
|
|
exp(−βεK ) |
|
|
K=0
есть т.н. статистическая сумма, играющая ключевую роль в статистической физике. В
нашем случае, когда εK = Khω¯ , получаем |
|
||
Z(β) = |
1 |
. |
(315) |
|
|||
1 − exp(−βhω¯ ) |
207
Теперь нам остается сделать заключительный шаг определить среднюю энергию фотонов на уровне n, соответствующей энергии hω¯ . Зная вероятность PK (313), с которой реализуется каждое состояние с полной энергией Khω¯ , имеем в результате
|
= ∞ |
Khω¯ |
exp (−βKhω¯ ) . |
(316) |
|||||||
Eω |
X |
|
|
|
|
|
|
|
|
|
|
K=0 |
|
|
|
Z(β) |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
Здесь и становится понятным удобство введения статистической суммы Z(β) (314). Дей- |
|||||||||||
ствительно, прямым вычислением можно показать, что |
|
||||||||||
|
Eω = − |
|
1 dZ(β) |
. |
|
|
(317) |
||||
|
|
|
|
|
|
|
|||||
|
Z(β) |
dβ |
|
|
|||||||
В итоге, умножая энергию Eω на число состояний δN и деля его на объем δV , получем |
|||||||||||
для спектральной плотности энергии |
|
|
|
|
|
|
|
|
|
|
|
w(ω, T ) = |
ω2 |
|
|
|
hω¯ |
|
|
, |
(318) |
||
π2c3 |
|
|
|
|
|||||||
|
|
exp (¯hω/T ) − 1 |
|
которая и соответствует планковскому распределению (20).
8.4Атом Бора против теории Максвелла
На примере определения времени жизни электрона в атоме Бора, которое следует из классической теории электромагнетизма, мы заодно покажем, как удобно проводить вычисления ”по порядку величины”, т.е. не обращая внимание на численные коэффициенты порядка единицы. Иногда подобный расчет совершенно достаточен, чтобы сказать, например, нужно или нет учитывать тот или иной эффект.
Итак, согласно теории Максвелла, любой заряд e, движущийся с ускорением r¨, будет терять свою энергию E со скоростью33
δE |
|
e2r¨2 |
. |
(319) |
|
||||
δt |
c3 |
|
Поскольку в атоме Бора кинетическая энергия по порядку величины совпадает с потенциальной, то мы можем записать E e2/aB, где aB (119) радиус Бора. С другой
33Эта формула выводится в книге ”Гравитация и астрофизика”.
208
стороны, ускорение электрона можно записать как r¨ = e2/a2Bme. В результате, имеем для оценки времени жизни τ
τ |
(δ |
E |
|
me2c3aB3 |
. |
(320) |
|
|
|||||||
|
E |
/δt) |
e4 |
|
|||
|
|
|
|
|
|
|
Конечно, для количественной оценки величины τ можно было бы сразу подставить
числовые значения заряда и массы электрона, радиуса Бора и скорости света. Однако иногда бывает удобнее несколько перегруппировать члены, объединив их в известные комбинации. Действительно, воспользовавшись определением радиуса Бора aB (119), по-
лучаем |
|
|
|
τ α−4 |
aB |
(321) |
|
|
. |
||
c |
Как мы видим, для оценки времени жизни нам нужно знать лишь значения постоянной тонкой структуры α = e2/hc¯ ≈ 1/137, радиуса Бора aB и скорости света c. Имеем в итоге
τ |
10−11 s. |
(322) |
||||
Формула (321) может быть переписана и в другой удобной форме |
|
|||||
|
1 |
|
c |
|
||
τ = |
|
|
|
|
, |
(323) |
ωc |
ωcre |
где ωc есть частота вращения, а re классический радиус электрона. Эта формула оказы-
вается более универсальной, поскольку она справедлива для любого кругового движения с частотой ωc.
8.5Cоотношение неопределенностей точная формулировка
Внимательный читатель должен был бы обратить внимание на то, что соотношение
неопределенностей Гейзенберга x px ≥ h/¯ 2 (60) (и, соответственно, соотношение для
волн x kx ≥ 1/2 (59)) безусловно, требуют уточнения. Действительно, справа стоит
четкая граница h/¯ 2 (или 1/2). Следовательно, и стоящие слева величины x и kx тре-
буют точного определения.
Как уже говорилось, значения x и kx должны быть определены как среднеквадратичные отклонения от их средних значений x0 и k0. Формальное определение x для
209
произвольной волны A = A(x), достаточно быстро спадающей при x → −∞ и x → ∞,
имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
(x |
x0)2 |
A(x) |
2dx |
|
||||||
(Δx)2 = |
R−∞ |
∞− |
|
A(x|) |
2dx |
| |
, |
(324) |
|||||
где среднее значение x0 есть |
|
|
|
R−∞ | |
|
| |
|
|
|
|
|
||
|
|
∞ x A(x) |
|
|
|
|
|
|
|||||
|
|
|
2dx |
|
|
|
|||||||
x0 |
= |
|
R−∞∞ |
|A(x) |
|2dx |
. |
|
|
(325) |
||||
|
|
|
R−∞ |
| |
|
| |
|
|
|
|
|
|
Эти соотношения на самом деле суть просто обобщения определения среднего арифметического для непрерывно распределенных величин.
Здесь нужно отметить два обстоятельства. Прежде всего, как мы уже отмечали, ”мерой” распределения должен быть именно квадрат амплитуды волны. В противном случае для знакопеременного сигнала можно получить совершенно неверный результат. В частности, для любого сигнала, описываемого нечетной функцией A(−x) = −A(x), вообще
невозможно определить величины, которые должны были стоять в знаменателях (324)– (325). Появление же модуля функции A(x) позволяет использовать это соотношение и
для комплексных сигналов. При этом, как мы видим, в определениях (324)–(325) вместо
A(x) может быть использована и волновая функция ψ(x).
Что же касается величин kx и px, то для них определние среднеквадратичных
отклонений на первый взгляд должно проводиться по совершенно различным формулам.
Если говорить о волновых векторах, то для определения характерного разброса |
kx |
|||||
нужно прежде всего сделать преобразование Фурье |
|
|
|
|||
1 |
∞ |
x |
|
|
||
Aˆ(kx) = |
√ |
|
Z−∞ A(x)e−ik |
xdx, |
(326) |
|
2π |
ˆ |
|
|
|
|
|
|
|
|
|
|
(величина A(kx) и определяет область и относительную роль волновых векторов, задаю- |
||||||||||
щих сигнал A(x)), а затем выполнить операцию, аналогичную (324) |
|
|||||||||
|
∞ |
(kx |
k0) |
2 |
ˆ |
|
2 |
dkx |
|
|
(Δkx)2 = |
R−∞ |
∞− |
ˆ |
| |
2 |
| |
|
|
. |
(327) |
|
|
R−∞ |A(kx)| |
dkx |
|
|
|
|
Здесь k0 среднее значение kx, определяемое по аналогии с (325). С другой стороны, по
канонам квантовой механики квадрат импульса, благодаря соотношению (92), следует
|
|
|
|
|
|
210 |
определять как34 |
|
|
|
|
|
|
2 |
2 |
∞ |
|
d2ψ |
|
|
px |
= h¯ |
Z−∞ |
ψ (x) |
|
dx. |
(328) |
dx2 |
Однако, как можно показать, оба этих определения будут приводить к одному и тому же значению px = h¯ kx.
В качестве примера рассмотрим пространственную волну A(x), задаваемую уже зна-
комой нам функцией Гаусса |
|
|
|
|
|
|
|
A(x) = |
1 |
|
exp |
− |
x2 |
!. |
(329) |
|
|
|
|||||
π1/4√ |
|
2x02 |
|||||
x0 |
Коэффициент перед экспонентой выбран таким образом, чтобы знаменатели в соотношениях (324)–(325) были равны единице. Поскольку же ”весом” в интеграле (324) будет квадрат амплитуды волны A(x), т.е. величина e−x2/x20 , не содержащая двойку в знамена-
теле экспоненты, то, как легко проверить,
|
x2 |
|
||
(Δx)2 = |
0 |
. |
(330) |
|
2 |
||||
|
|
|
С другой стороны, замечательное свойство функции Гаусса g(x) e−x2/2 состоит в том, что она совпадает по форме со своим преобразованием Фурье (326): gˆ(k) e−k2/2. Поэто-
му |
|
√ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x2k2 |
!. |
|
|
|
Aˆ(kx) = |
x0 |
|
|
|
|
|
||||
|
|
|
exp |
− |
0 x |
(331) |
|||||
|
π1/4 |
2 |
|||||||||
Отсюда уже нетрудно получить, что |
|
|
|
|
|
|
|
||||
|
(Δkx)2 = |
|
1 |
, |
|
|
(332) |
||||
|
|
2 |
|
|
|||||||
|
|
|
|
|
|
2x |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
и, следовательно, x px = 1/2. Читателю самому предоставляется возможность убе-
диться, что формула (328) приводит к тому же результату. Таким образом, можно сделать вывод, что для гауссова профиля реализуется наименьшее из значений, разрешенных соотношением неопределенностей Гейзенберга.
34Мы для порядка поставили звездочку, которая означает, что для комплексных волновых функций
здесь должна стоять величина, комплексно сопряженная к ψ(x).