Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Lec, Beskin / beskin. book

.pdf
Скачиваний:
38
Добавлен:
03.06.2015
Размер:
677.94 Кб
Скачать

131

(a + b)3 = (a + b)(a + b)(a + b) два раза выпадет a, и один раз b. Выпадение a будет

соответствовать событию, в котором частица попадает в левую половину комнаты (а выпадение b в правую). При этом сами коэффициенты задаются следующей формулой

(доказывается по индукции)

CNK =

N !

(190)

 

,

 

 

K!(N K)!

 

где по определению факториала n! = 1 ·2 ·. . . ·n, но и 0! = 1. Эту величину еще называют статистическим весом и обозначают буквой G. Полное же число событий есть сумма коэффициентов при произведениях aK bN −K . Положив a = b = 1, получаем (1 + 1)N = 2N .

Для определения же собственно вероятности нам потребуется еще одна гипотеза. А именно гипотеза о том, что выпадение всех 2N возможных состояний имеют одну и ту

же вероятность. Именно в этом и состоит основная аксиома статистической физики. Как уже отмечалось, она является в значительной степени эвристической. Однако ее справедливось (т.е. справедливость выводов и следствий, основанных на этой гипотезе) была многократно проверена. Поэтому, как и для закона сохранения энергии, принцип равно-

распределения был возведен в ранг универсального закона.

Если же все варианты действительно одинаковы, то тогда мы можем написать для вероятности события, в котором в левой части находится K частиц, а в правой (N K)

частиц

 

 

 

PNK =

N !

.

(191)

 

2N K!(N K)!

 

 

 

Обсудим вкратце свойства этого распределения. Поскольку нас прежде всего интересуют большие значения K и N , то можно воспользоваться соотношением Муавра-Стирлинга

 

 

 

n

n

 

n! 2πn

 

 

 

.

(192)

e

Для дальнейшего нам будет удобнее использовать другую форму

 

ln n! n ln n n + ln(2πn)/2.

(193)

Эти соотношения, как мы видим, справедливы не только для целых, но и для любых значений n.

132

Далее, как видно уже из биноминальных коэффициентов, наибольшие значения вероятности должны достигаться при K = N/2. Если теперь ввести новую переменную

δ = K N/2 и воспользоваться соотношением ln(1 + δ) δ δ2/2, то можно получить

после несложных преобразований

 

 

 

 

 

 

 

 

P (δ) = s

 

 

 

 

 

2

 

 

 

2

 

exp

 

.

(194)

 

πN

N

 

 

 

!

 

Это соотношение и было использовано нами в разделе 1.1.1.

6.2Закон Релея-Джинса

Прежде чем переходить к вычислениям, необходимо точно определить все используемые величины. Закон Рэлея-Джинса это закон для равновесной плотности излучения т.н.”абсолютно чёрного тела”. Здесь спектральная плотность излучения P(ω) определяется так, чтобы поток энергии (измеряется в эрг см−2 с−1 Гц−1) в диапазоне частот от

ω до ω + δω был равен P(ω)δω. Далее, понятно, что абсолютно чёрных тел в природе не

существует. ”Абсолютно чёрное тело” есть удобная модель, которая обычно используется для анализа экспериментов. Обычно ее изображают как замкнутую полость с небольшим отверстием. Но на самом деле отверстие для этой модели вообще не важно, поскольку оно нужно лишь для того, чтобы подчеркнуть принципиальную возможность наблюдать излучение, находящегося внутри полости. В результате свет, попадающий внутрь сквозь отверстие, после многократных отражений будет полностью поглощён, и отверстие снаружи будет выглядеть совершенно чёрным.

При этом, однако, важно, что при нагревании этой полости у неё должно появится собственное излучение. Разумно предположить, что излучение, испущенное внутренними стенками полости, прежде, чем оно сможет выйти за пределы полости, должно неизбежно претерпеть большое количество новых поглощений и излучений. Поэтому излучение внутри полости будет находится в термодинамическом равновесии со стенками. Это и есть основное свойство”абсолютно чёрного тела”, которое и приводит к его, как мы сейчас увидим, достаточно универсальным свойствам.

133

Еще одним важнейшим элементом, лежащим в основе вывода закона Релея-Джинса, стала вновь идея о равнораспределении энергии по степеням свободы. Напомним, что как было показано в разделе 1.2.2, количество различных волн в малом объеме δV = δxδyδz

в диапазоне волновых векторов между kx и kx + δkx, ky и ky + δky , kz и kz + δkz есть (64)10

δN = 2

δV δkxδky δkz

.

(195)

 

 

(2π)3

 

Иными словами, число δN и определяет число степеней свободы для электромагнитного поля в элементе объема δV и в диапазоне волновых векторов δkxδky δkz . Согласно же

принципу равнораспределения, все они являются равновероятными.

Вспомним теперь, что мы обсуждаем излучение ”абсолютно черного тела”, для которого мы в качестве модели выбрали излучение в замкнутой полости. Понятно, что для такого излучения естественно предположить, что оно будет изотропным, т.е. что поток излучения через небольшую поверхность не будет зависеть от ориентации этой поверхности. В этом случае будет иметь смысл лишь зависимость излучения от величины волнового вектора k, а не от его направления. Поэтому элемент ”объема” в трехмерном пространстве волновых векторов δkxδky δkz удобно, как и в нашем обычном пространстве, заменить на 4πk2δk. Вспоминая теперь, что частота ω и волновой вектор k для фотонов связаны соотношением ω = ck, получаем окончательно

δN =

ω2

δV δω.

(196)

π2c3

 

 

 

Предположим теперь, как это сделали Релей и Джинс, что каждому колебанию можно приписать энергию

ε = T

(197)

(напоминаем, что мы везде измеряем температуру в энергетических единицах). Помножив теперь плотность вероятности δN (196) на T , и поделив на объем δV , мы получим

10Здесь дополнительный фактор 2 появился за счет того, что электромигнитные волны (в отличие, например, от звуковых) описываются векторными величинами напряженностями полей. В результате, электромагнитная волна имеет две независимые поляризации.

134

энергию δǫ, заключенную в единице объема в диапазоне частот между ω и ω + δω. Определяя теперь эту энергию как δǫ = w(ω, T )δω, получаем окончательно для плотности энергии w(ω, T ), которая приходится на интервал частот δω:

w(ω, T ) =

ω2T

.

(198)

π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. В итоге, получаем

P(ω, T ) =

ω2

(199)

2c2 T.

Как уже отмечалось, это соотношение с хорошей точностью описывает излучение черного тела при hω¯ T .

6.3Формула Планка

Для вывода формулы Планка (14) нам потребуется результат, основанный на полной версии принципа равнораспределения энергии по степеням свободы. А именно распределение Гиббса одно из ключевых соотношений статистической физики. Как мы уже отмечали выше, в Приложении 6.1 был рассмотрен случай, когда все состояния имеют одинаковую энергию. Если же различные состояния имеют разную энергию, то уже интуитивно ясно, что частицам будет выгодно занимать нижние энергетические уровни. Тем не менее, принцип равнораспределения может быть сформулирован и этом случае.

Дело в том, что физически правильной постановкой задачи будет вопрос о том, каким будет распредение N частиц по уровням энергии εn при заданной величине полной

 

135

энергии

 

K

 

X

(200)

E = Nnεn.

n=1

Здесь Nn число частиц, занимающих энергетический уровень εn, а K количество

различных уровней энергии. Естественно, что при этом должно быть выполнено еще и соотношение

 

K

 

 

 

X

 

(201)

N = Nn.

 

 

n=1

 

 

Ответ на этот вопрос дает формула

 

 

 

G =

N !

,

(202)

 

N1!N2! . . . NK !

которая обобщает соотношение (190). При этом величина G (т.н. статистический вес) показывает, сколькими способами можно разместить N1 частиц в первой ячейке, N2 частиц

во второй, и т.д. Здесь и используется фундаментальное предположение о том, что все состояния с данной энергией равновероятны.

На Рис. 18 для иллюстрации показано, как можно разместить шесть одинаковых молекул в трех ячейках, имеющих энергии 0, ε и , так, чтобы полная энергия системы составляла . Таких комбинаций, как мы видим, всего три. Но количество вариантов,

которыми можно расположить молекулы (т.е. статистические веса соответствующих состояний) при этом различны.

Используя теперь, как и в Приложении 6.1, формулу Муавра-Стирлинга (192), а также (что очень важно!) соотношения (200) и (201), можно показать, что

εn

,

 

Nn = N0 exp T

(203)

где N0 некоторая постоянная. Это и есть соотношение, найденное Дж. Гиббсом (1839-

1903) и носящее его имя. Как и в случае распределения Гаусса (3), оно получается в пределе очень больших значений N . Действительно, при малых N , как показано на Рис. 18,

вероятности различных состояний хотя и отличаются друг от друга, но все же могут совпадать по порядку величины. При больших же N , как и для распределения Гаусса,

практически реализуется лишь наиболее вероятное состояние. Впрочем, и на нашем, чрезвычайно упрощенном примере видно, что наиболее вероятным является состояние, при

136

.

.

Рис. 18: Различные варианты размещения шести одинаковых молекул в трех ячейках, имеющих энергии 0, ε и , так, чтобы полная энергия системы составляла . Показаны также статистические веса G и вероятности выпадения соответстенных состояний

137

котором числа заполнения энергетических уровней плавно уменьшаются с увеличением их энергий.

Здесь, правда, неообходимо отметить еще одно обстоятельство. На первый взгляд, непонятно, в каком месте в наших рассуждениях могла возникнуть температура T . Дело

втом, что на самом деле при анализе соотношений (200)–(202) величина 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.

 

 

(204)

 

 

 

 

 

 

dh

T

 

 

 

Решение этого дифференциального уравнения есть

! .

 

n(h) = n0 exp

T

(205)

 

 

 

 

 

 

mpgh

 

 

Как мы видим, форма уранения (204) совершенно аналогична распределению Гиббса, причем в числителе стоит как раз потенциальная энергия mpgh, а в знаменателе температура T .

Теперь, наконец, мы можем вернуться к идее Планка. Как мы помним, она основывалась на предположении, согласно которому энергия осциллятора может принимать лишь значения, кратные кванту энергии εq, где εq = hω¯ . В этом случае величины Nn (203)

образуют геометрическую прогрессию, сумма которой можно легко определить. Напомним, что сам Планк говорил лишь о квантовании энергии осцилляторов, составляющих вещество стенок полости, в котором находится излучение; при этом, конечно же, предпо-

138

лагалось, что излучение и вещество находятся в термодинамическом равновесии11. Нам

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

K фотонов может быть записана в виде Khω¯ . Как легко проверить, точная формула

(87), содержащая фактор (K + 1/2), приводит к тем же результатам.

Итак, нам вновь нужно посчитать энергию, которую содержат фотоны в диапазоне частот от ω до ω + δω. При этом, как и в пределе Релея-Джинса, число состояний в

элементе объема δV будет задаваться формулой δN = (ω22c3)δV δω (196). Отличие

будет состоять лишь в том, что энергия разных колебаний будет отличаться от T .

Действительно, давайте применим общий принцип, найденный Гиббсом, для фотонов, имеющих одинаковую энергию hω¯ . C учетом общего соотношения (203), вероятность того, что на уровне n с энергией hω¯ будет находиться K фотонов, где K = 0, 1, 2, . . ., может

тогда быть записана в виде

 

T

!

,

 

(206)

PK = P0 exp

 

 

 

 

Khω¯

 

 

 

 

где P0 некоторая константа. При этом сумма геометрической прогрессии PK=0 PK легко

вычисляется, и мы получаем

 

 

 

 

 

 

 

 

 

 

hω¯

!#

−1

(207)

PK = P0

"1 exp T

.

X

 

 

 

 

 

 

 

 

K=0

Поскольку же сумма всех вероятностей должна быть равна единице, получаем окончательно

PK =

"1 exp T !# exp

T

! ,

 

 

hω¯

 

Khω¯

 

Для дальнейшего нам будет удобно записать это выражение в виде

1

PK = Z(β) exp (βKhω¯ ) .

Здесь по определению β = 1/T , а

X

Z(β) = exp(βεK )

K=0

(208)

(209)

(210)

11Излучение в комнате, в которой Вы находитесь, этому свойству не удовлетворяет!

139

есть т.н. статистическая сумма, играющая ключевую роль в статистической физике. В нашем случае, когда εK = Khω¯ , получаем

 

1

(211)

Z(β) =

1 exp(βhω¯ ) .

Теперь нам остается сделать заключительный шаг определить среднюю энергию фотонов на уровне n, соответствующей энергии hω¯ . Зная вероятность PK (209), с которой реализуется каждое состояние с полной энергией Khω¯ , имеем в результате

 

=

Khω¯

exp (βKhω¯ ) .

 

(212)

En

X

 

 

 

 

 

 

 

 

 

 

K=0

 

 

 

Z(β)

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь и становится понятным удобство введения статистической суммы Z(β) (210). Дей-

ствительно, прямым вычислением можно показать, что

 

 

 

 

 

 

1 dZ(β)

 

(213)

 

En =

 

 

 

 

.

 

 

Z(β)

 

В итоге, имеем для плотности энергии

 

 

 

 

 

 

 

 

 

 

w(ω, T ) =

ω2

 

 

 

hω¯

,

(214)

π2c3

exp (¯hω/T ) 1

 

 

 

 

которая и соответствует планковскому распределению (14).

6.4Атом Бора против теории Максвелла

На примере определения времени жизни электрона в атоме Бора, которое следует из классической теории электромагнетизма, мы заодно покажем, как удобно проводить вычисления ”по порядку величины”, т.е. не обращая внимание на численные коэффициенты порядка единицы. Иногда подобный расчет совершенно достаточен, чтобы сказать, например, нужно или нет учитывать тот или иной эффект.

Итак, согласно теории Максвелла, любой заряд e, движущийся с ускорением a, будет терять свою энергию E со скоростью12

δE

 

e2a2

.

(215)

 

δt

c3

 

12Эта формула выводится в книге ”Гравитация и астрофизика”.

140

Поскольку в атоме Бора кинетическая энергия по порядку величины совпадает с потенциальной, то мы можем записать E e2/aB, где aB (112) радиус Бора. С другой стороны, ускорение электрона можно записать как a = e2/a2Bme. В результате, имеем для оценки времени жизни τ

τ

E

 

me3aB3

.

(216)

 

 

E

/δt)

e4

 

 

 

 

 

 

 

 

Конечно, для количественной оценки величины τ можно было бы сразу подставить

числовые значения массы электрона, радиуса Бора и скорости света. Однако иногда бывает удобнее несколько перегруппировать члены, объединив их в известные комбинации.

Действительно, воспользовавшись определением радиуса Бора aB (112), получаем

 

 

hc¯

!

5 h¯ 1

(217)

τ =

 

 

 

 

 

.

e2

mec

c

Как мы видим, для оценки времени жизни нам нужно знать лишь значения постоянной тонкой структуры α = e2/hc¯ 1/137, комптоновской длины волны λC и скорости света c. Имеем в итоге

τ 10−11s.

(218)

6.5Cоотношение неопределенностей точная формулировка

Внимательный читатель должен был бы обратить внимание на то, что соотношение

неопределенностей Гейзенберга x px h/¯ 2 (53) (и, соответственно, соотношение x kx

1/2 (52) для волн) безусловно, требуют уточнения. Действительно, справа стоит четкая

граница h/¯ 2 (или 1/2). Следовательно, и стоящие слева величины x и kx

требуют

точного определения.

 

Как уже говорилось, значения x и kx должны быть определены как среднеквад-

ратичные отклонения от их средних значений x0 и k0. Формальное определение

x для

произвольной волны A = A(x), достаточно быстро спадающей при x → −∞ и x → ∞,

имеет вид

 

 

 

 

 

 

 

(x x0)2

A(x)

2dx

 

(Δx)2 =

R−∞

A(x|)

2dx |

,

(219)

 

 

R−∞ |

 

|

 

 

 

Соседние файлы в папке Lec, Beskin