поляризацияволнылинейнаяивектор Eω0 направленвдольоси Х,тосоставляющаясилы Fx = eEω0 cos(ωt −2πz /λ) = −∂∂Ux
и U (1) = −exEω0 cosωt , а значит Hˆ (1)(t) = −exEˆ ω0 cosωt .
Подсчитаем матричные элементы оператора Hˆ (1)(t). По определению:
Hmn(1)(t) = ∫ψm* Hˆ (1)ψndV = −eEω0 cosωt ×∫ψm* xˆψndV = = −exmn Eω0 cosωt ,
если ввести обозначение: exmn = e∫ψm* xψndV .
Чтобы уяснить смысл этого выражения рассмотрим случай
равных индексов n = m. Тогда exmm = e∫ψm* xψmdV = ex , т.е. эта величина есть средний дипольный момент атома в состоянии m в направлении х, если один из зарядов, например, ядро
поместить в начале координат. В формуле для Hmn(1) индексы смешанные, и мы можем интерпретировать это как средний «дипольный момент перехода» межу состояниями с номерами n и m. Совокупность величин xmn дает матрицу:
|
x11 |
x12 |
x13......x1n ... |
∫ |
x |
x |
x |
.......x ... |
|
22 |
23 |
2n |
|
xjk = |
ψ*j xψk dx 21 |
. |
|
x31 |
x32 |
x33........x33... |
|
|
xn1 |
xn2 |
xn3 |
........xnn .. |
|
Таким образом, производные коэффициентов am(1)(t) могут быть представлены в виде:
dadtm(1) = − i (−exmn Eω0 cosωt)eiωmnt =
= 2i exmn Eω0[ei(ωmn +ω)t +ei(ωmn −ω)t ] .
Чтобы получить am(t) нужно проинтегрировать от 0 до t:
a |
(1) |
(t)= |
|
i |
ex |
|
E |
0 |
|
t |
ei(ωmn +ω)τ |
+ |
ei(ωmn −ω)τ |
|
= |
|
|
|
m |
|
|
|
mn |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
ω |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
i(ωmn +ω) |
|
i(ωmn −ω) |
|
= |
|
1 |
ex |
|
E |
0 ei(ωmn +ω)t −1 |
+ |
ei(ωmn −ω)t −1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
mn |
ω |
(ωmn +ω) |
|
(ωmn |
|
|
|
|
|
|
|
|
|
|
−ω) . |
|
|
Квадратыэтихкоэффициентовописываютвероятностипереходов между уровнями с номерами n и m. Если Em>En, то переход n m происходит с поглощением энергии и ωmn>0; если Em<En, то переход n m происходит с излучением энергии и ωmn<0. Поскольку, как мы увидим, в физически содержательных случаях |ωmn| ω, в случае поглощения энергии первый член в квадратных скобках много меньше второго, в случае излучения – наоборот.
Рассмотрим случай излучения. Тогда квадрат величины am(t) определяет, очевидно, вероятность излучения кванта за время t:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
(t) |
= a |
|
|
(t)a |
* |
(t) = |
1 e2 |
|
xmn |
|
2 Eω2 |
( |
i(ωmn +ω)t |
)( |
−i(ωmn +ω)t |
) |
|
|
|
|
|
|
|
|
|
|
m |
m |
m |
4 2 |
|
(ωmn +ω)2 |
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
|
−1 e |
|
|
−1 = |
= |
|
|
2 |
|
e2 |
|
xmn |
|
2 Eω2 |
[1−cos(ω |
mn |
+ω)t]= |
Eω2e2 |
|
xmn |
|
2 |
|
sin2 (ωmn +ω)t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 2 (ωmn +ω)2 |
2 (ωmn +ω)2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
Таким образом, вероятность перехода пропорциональна квадрату напряженности поля и квадрату дипольного момента перехода (поэтому оно и называется дипольным). Из формулы
видно, что вероятность для любого момента времени имеет острый пик при |ωmn| ω. Чтобы убедиться в этом, введем за-
мену: |
(ωmn +ω)t |
|
|
=ξ , тогда: |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
2 |
|
(t) = |
Eω2e2 |
|
xmn |
|
t |
2 sin2 |
ξ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
4 2 |
|
|
|
ξ2 |
|
|
|
|
|
|
|
|
|
|
|
Исследуем поведение множителя f (ξ)= sin2 ξ .
ξ 2
(sin(x)**2)/x**2
|
Рис. 24 |
|
Приξ 0,(применяяправилоЛопиталя) f (ξ)= |
2sinξcosξ →1; |
при ξ=π/2, |
|
2ξ |
f (ξ) =(2/π)2 0,4. Следующий максимум при |
ξ=3π/2, |
f (ξ) =(2/3π)2=0,045; третий |
при ξ=5π/2, |
f (ξ) =(2/5π)2=0,016 и т.д.
Максимум очень узкий, т.к. при ω~1015 и t 10-8 сек (время “жизни”) ξуже больше 107, а f (ξ)<10−14 !
Поэтому заметная вероятность перехода под действием вынуждающей силы имеет место в очень узком интервале частот вблизи ωmn, и атомные колебательные системы имеют очень высокую (~107) добротность. По своему смыслу вероятность должна быть пропорциональна времени, а не квадрату синуса! Для получения правильного результата нужно учесть еще немонохроматичность реальной волны и ширину (пусть малую!) спектральной линии. С этой целью проинтегрируем полученную формулу в пределах ширины спектральной линии по интервалу частот. Поскольку ширина очень узка, интегрировать можно в пределах от 0 до . Воспользуемся уже введенной заменой переменных и соотношением их дифференциалов
dω = 2dξt :
∞ |
|
|
|
2 |
|
|
|
|
Eω2e2 |
|
xmn |
|
2 ∞ |
2 sin2 ξ |
|
2dξ |
|
|
|
|
|
|
|
|
|
|
|
|
∫ |
am (t) |
|
dω = |
|
|
∫t |
|
= |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
4 |
2 |
|
|
|
ξ |
|
|
|
|
t |
0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
Eω2e2 |
|
xmn |
|
2 |
|
∞ sin2 |
ξ |
|
|
|
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
t(∫0 ξ2 |
|
dξ ={2}). |
|
|
|
|
|
|
|
|
2 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
πE2e2 |
|
x |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Окончательно |
|
w |
= |
a (t) |
|
= |
|
ω |
|
|
|
mn |
|
|
t . Вероят- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
mn |
|
m |
|
|
|
|
|
|
4 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ность можно выразить и через плотность энергии электромаг-
|
E2 |
|
H 2 |
|
E2 |
|
(E0 )2 |
нитногополя.ВсистемеСГС |
|
+ |
|
= |
|
= |
ω |
= ρ . |
|
|
|
|
8π |
|
8π |
|
4π |
|
8π |
ω |
|
|
|
|
|
Откуда wmn = 2π2e2 xmn 2 t ρω , где ρω – плотность энергии в
2
плоской электромагнитной волне на частоте ω.
Дипольное излучение и поглощение
Лекция 21
§59 Матрица дипольных переходов. Правила отбора их физический смысл.
«Разрешенные» и «запрещенные» переходы. Характер поляризации излучения
Итак, вероятность перехода связана с матричным элементом матрицы дипольного момента. Например, для радиальной сим-
метрии: erjk = e∫ψ*j rˆψk dV . Для многих состояний эти элементы оказываются равными нулю. Рассмотрим уже известные
системы.
∞
а). Гармонический осциллятор. xmn = ∫ψm* xψndx .
−∞
Полученные нами ранее волновые функции име-
|
|
|
|
ψn (x) = Hn (x)e− |
αx2 |
|
|
|
где Hn (x) |
|
ют |
вид: |
|
2 |
, |
|
|
полино- |
мы |
Чебышева, |
|
подчиняющиеся |
|
рекуррентной |
фор- |
муле: |
|
|
|
|
xHn (x) = nHn−1(x) |
+ |
1 Hn+1(x). |
Поэто- |
|
|
|
|
|
|
|
|
|
|
|
1 H |
|
|
2 |
|
|
му |
H |
m |
xH |
n |
= H |
m |
nH |
n−1 |
+ |
m |
H |
n+1 |
. В силу |
орто- |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
гональности |
|
|
полиномов, |
|
|
|
матричный |
элемент |
xmn = ∞∫ Hm xHndx = ∞∫ |
HmnHn−1dx + |
1 ∞∫ Hm Hn+1dx ≠ 0, |
−∞ |
−∞ |
|
2 −∞ |
только при условии, что m=n+1 или m=n–1, значит, возможны переходы только при n n 1, т.е. n= 1, откуда
ωmn = Em − En = ∆n ω0 =ω0 . Поэтому в линейном осцилляторе возможны переходы только между соседними уровнями
с излучением одной единственной частоты. |
|
б). Центральное поле. |
|
|
Матричные |
элементы |
координат: |
xjk = ∫ψ*j xψk dV , yjk = ∫ψ*j yψk dV , z jk = ∫ψ*j zψk dV .
Общий вид собственных функций оператора момента ко-
личества движения мы получили ранее: |
ψ = Plm (θ)eimϕ . |
Ввиду независимости |
радиальной |
функции от |
углов можно положить |
r=1. Тогда |
x = cosϕcos−θ1θ, |
y = sinϕsinθ, z = cosθ, dV = sinθdθdϕ.
Пусть в состоянии j и k квантовые числа имеют значения l’, m’ и l, m, соответственно. Тогда, подставляя функции в интегралы, имеем:
xjk = |
12∫π [ei(m−m'+1) +ei(m−m'−1)]dϕ×π∫Plm (θ)Pl 'm' (θ)sin2θdθ |
|
2 0 |
0 |
yjk = |
|
1 |
2∫π [ei(m−m'+1) −ei(m−m'−1)]dϕ×π∫Plm (θ)Pl 'm' (θ)sin2θdθ |
|
|
|
|
2i 0 |
0 |
z jk = |
2∫π ei(m−m')dϕ × |
π∫Plm (θ)Pl 'm' (θ)sinθ cosθdθ |
|
0 |
|
0 |
Для переходов с ∆m=0 матричные элементы
Последний интеграл по φ не равен нулю только при m=m’, откуда ∆m=0. Первые два не равны нулю только при ∆m= 1.
xjk = yjk = 0,
т.е. при таких переходах колебания происходят «вдоль z” и возникающая волна поляризована линейно.
Для∆m= 1матричныйэлементz jk =0,а xjk и yjk «сдвинуты по фазе» на π/2 (из-за множителя 1/i), возникающая волна поляризована по кругу. Для рассмотрения второго интеграла ограничимся матрицей z jk , не равной нулю при m=m’, которую можно записать в виде:
z jk = 2ππ∫Plm (cosθ) Pl 'm (cosθ)sinθ cosθdθ .
|
0 |
z jk |
|
Тогда, если ввести η = cosθ , формула для |
примет вид: |
z jk = 2π ∫1 |
ηPlm (η) Pl 'm (η)dη. Напомним, |
что |
Plm (η) суть, |
−1 |
|
|
так называемые, шаровые функции, для которых в математике доказана рекуррентная формула:
ηPlm (η) = 2l l++m1Pl−1,m (η)+ l −2lm++11Pl+1,m (η).
Опять в силу ортогональности:
|
1 |
l +m |
|
l −m +1 |
|
|
|
zjk = 2π |
|
|
|
Pl−1,m + |
|
Pl+1,m |
Pl 'mdη ≠ 0 |
|
2l +1 |
2l +1 |
|
−∫1 |
|
|
|
|
только при l’=l 1, т.е. l= 1.
Все, что мы получили для момента количества движения верно и для атома водорода и водородоподобных ио-
нов, но нужно еще учесть радиальную функцию, поскольку
ψnlm = Rnl (r)Plm (θ)eimϕ и dV = r2 sinθdθdϕdr . к рассмо-
тренным уже интегралам добавится множитель вида:
∞∫Rn'l ' (r) rRnl (r)r2dr .
0
Таким образом, для l и m правила отбора сохранятся, а при их выполнении последний интеграл не обращается в нуль при любых n и n’. Следовательно, все переходы по n разрешены.
ЛЕКЦИЯ 22
§60 Эмпирические законы теплового излучения. Закон Вина. Формула Релея-Джинса. «Ультрафиолетовая катастрофа». Формула Планка для плотности излучения в термодинамическом равновесии
Одной из проблем, которая не подавалась решению в рамках классической физики в «допланковское» время, была проблема излучения нагретых тел, или теплового излучения (не инфракрасного, а излучения, возбуждаемого тепловым движением атомов). Еще в XIX-м веке экспериментально были установлены внешние закономерности этого явления. Распределение интенсивности излучения по частотам (длинам волн) имело вид характерной кривой с максимумом, с убыванием как в сторону малых, так и в сторону больших длин волн или частот. Высота максимума пропорциональна пятой степени абсолютной температуры, а его положение на шкале длин волн определяется
формулой: λ(max) = C /T (закон смещения Вина). Было установлено, что полная энергия, излучаемая единицей поверхно-
сти нагретого тела за секунду на всех длинах волн пропорциональна четвертой степени абсолютной температуры u =σT 4 (закон Стефана-Больцмана).
Чтобы рассмотреть равновесие тела с окружающим его излучением, введем понятие поглощательной и излучательной способности тела. Первая показывает, какая доля (Aν) падающей в интервале частот ν, ν+Δν энергии остается внутри тела, превращаясь в тепло. Тело у которого Aν=const стали называть «серым», а если Aν=1, – «абсолютно черным» или просто «черным».Вэкспериментахчерноетелообычномоделируютзачер-
ненной полостью, часто сферической, с малым отверстием, в которое свет попадает почти по касательной к поверхности и полностью поглощается при многократных отражениях. Излучательная способность (Eν) – энергия излучаемая 1 см2 нагретого тела за секунду в том же интервале частот.
Еще в середине XIX столетия Кирхгофом с использованием второго закона термодинамики было показано, что при равновесии рассматриваемого тела с окружающим его излучением выполняется соотношение:
Eν = 8cπ ρν ,
Aν
где ρν – спектральная плотность излучения в объеме. Это указывает, что рассматриваемое отношение не зависит от тела и спектральная плотность излучения, находящегося в равновесии с телом, есть универсальная функция температуры, т.е.:
ρ(ν)dν = F(ν,T )dν .
Виномбылополучено,чтоформула,удовлетворяющаязаконам термодинамики должна иметь вид (функция Вина):
ρ(ν)dν =ν3F ν dν .
T
Эта формула оказалась очень полезной, т.к. подтвердила, по крайней мере, две экспериментальные закономерности. ИнтегрированиеэтойформулыповсемчастотамдаетзаконСтефана Больцмана:
∞ |
∞ |
|
|
ν |
|
∞ |
|
|
|
|
u = ∫ |
ρ(ν)dν = ∫ν |
3 |
4 |
∫ξ |
3 |
F(ξ)dξ =σT |
4 |
, |
|
F |
dν =T |
|
|
|
0 |
0 |
|
T |
|
0 |
|
|
|
|