НАИБОЛЕЕ ТИПИЧНЫЕ ОШИБКИ, ДОПУСКАЕМЫЕ В КВАЗИСОВРЕМЕННОЙ АБСТРАКТНОЙ ФИЗИКЕ – 5
РЕШЕНИЕ НЕКОТОРЫХ КЛЮЧЕВЫХ ЗАДАЧ АТОМНОЙ ФИЗИКИ В РАМКАХ КЛАССИЧЕСКОЙ СТАТИСТИЧЕСКОЙ ФИЗИКИ
§ 29. Электромагнитный механизм дифракции микрочастиц на
монокристаллах.
На первый взгляд, в рамки обычных
классических представлений не укладывались закономерности, которые проявлялись
при отражении любых микрочастиц от граней совершенных монокристаллов. Частицы
проявляли себя так же, как и рентгеновские лучи с длиной волны, равной длине
волны де Бройля
h/mv, для которых выполнялись условия Вульфа-Брэгга при
отражении от кристаллических плоскостей. Особенно знаменательным было то, что
существование гипотетических волн микрочастиц (волн материи) было предсказано де
Бройлем за два года до экспериментов по дифракции микрочастиц на монокристаллах.
Попытаемся произвести анализ этого явления, ставшего поначалу сенсационным и
необычным, а в дальнейшем просто роковым для судеб физики.
При рассмотрении каких бы то ни было моделей дифракции микрочастиц
в результате их взаимодействия с внешними макрообъектами следует учитывать,
прежде всего, те экспериментальные данные, которые можно отнести к
разряду твердо
установленных фактов. К настоящему времени с
высокой степенью точности и воспроизводимости результатов констатируется
следующее:
1.
Явления дифракции характерны для микрочастиц любой природы – электронов,
протонов, нейтронов, а также для атомов и молекул, за что их и нарекли
своеобразными волнами материи. Наличие у частиц
заряда или его отсутствие может сказаться на коэффициенте отражения, но не на
характере дифракционной картины. Здесь, пожалуй, можно опустить из рассмотрения,
например, эффекты рассеяния p-мезонов на протонах,
которые также предполагается интерпретировать с точки зрения дифракционных
механизмов.
2.
Дифракция микрочастиц имеет в общем
случае не поверхностный, а скорее объемный характер, обнаруживаясь при
прохождении через монокристаллы, облете препятствий. В случае же отражения от
поверхности монокристаллов картина дифракции в большей степени определяется
физической природой монокристалла и в меньшей степени – состоянием его
поверхности, в частности, процессами адсорбции или концентрацией дефектов на
поверхности. Последние можно рассматривать как малые возмущения к основной
картине дифракции на монокристалле, обусловленной его структурой.
3.
Доминирующим фактором дифракции является величина относительной скорости между
микрочастицей и макрообъектом. Если же говорить точнее, то для системы
координат, связанной с монокристаллом, главным является импульс микрочастицы. Но
эксперимент можно поставить так, что монокристалл будет двигаться с некоторой
скоростью навстречу частицам. В том случае, когда будут двигаться навстречу друг
другу и микрочастица, и монокристалл, не совсем ясно, что понимать под длиной
волны де Бройля в разных системах отсчета, не говоря уже о механизме
возникновения такой волны.
4. В
данных экспериментах отмечается поразительная корреляция положения дифракционных
максимумов от кристаллографических характеристик макрообъектов-мишеней, от
взаимной ориентации векторов импульса микрочастицы и кристаллической
решетки.
Последнее обстоятельство традиционно является самым сильным
доводом в пользу сложившейся концепции корпускулярно-волнового
дуализма, а также обоснованием де-бройлевского
формализма. Оставляя вне критики исходные пункты квантовой механики, основанные
Борном, де Бройлем, Бором, Шредингером, Гейзенбергом и др. и приведшие, в конце
концов, к возникновению целых направлений в фундаментальных исследованиях (от
квантовой химии до прикладной математики), можно попытаться дать альтернативные
интерпретации дифракционных механизмов и моделей.
Стало уже почти общим утверждение о нематериальности волн де
Бройля, т.е. их существование оценивается почти
исключительно с позиций эффективности математической модели. Можно согласиться,
что это действительно фантастично и плодотворно, но, с точки зрения материальной
физики, выглядит не вполне достаточной и непротиворечивой моделью. Отметим лишь
некоторые принципиальные моменты:
1.
Волнам де Бройля присущ широкий спектр необычных свойств – вакуумная
дисперсия, пространственная локализация,
необычные соотношения групповой и фазовой скоростей. Неясна даже гипотетически
их природа, и, как следствие, "волны де Бройля не имеют ничего общего с волнами,
рассматриваемыми в классической физике" [1]. Если волны де Бройля – это лишь
статистическое (общепринятое) описание движения микрочастиц, то ведь должен,
кроме количественного статистического описания, присутствовать реальный
“природный” механизм рассеяния, который хотя бы приблизительно способен был
объяснять наличие дифракции.
2. Мягко
говоря, дискуссионным подходом к проблеме можно было бы назвать решительное
нерассмотрение
электромагнитной природы взаимодействия
микрочастиц с макрообъектом. При условиях, сопутствующих явлениям дифракции,
вряд ли могут быть отменены (и не отменяются) законы классической механики и
электродинамики. Тем не менее, в рамках де-бройлевских моделей
электростатическое взаимодействие с периодическими полями кристаллической
решетки никак не отражается. Это не может не вызвать вопросов, потому что,
например, при различных твердотельных эффектах, связанных с миграцией электронов
по решетке, такие взаимодействия в основном только и рассматриваются. При этом
вполне естественным выглядит требование того, чтобы распределение электронных
облаков имело ту же самую пространственную периодичность, которая характерна в
целом для атомов кристалла.[1] Здесь достаточно
вспомнить соответствующие теоремы Блоха и Ванье [2].
3. Среди
экспериментальных данных, наглядно отражающих кристаллографию мишеней, обычно
привлекают внимание эффекты каналирования электронов и протонов вдоль атомных
плоскостей, а также эффекты теней от плоскостей кристалла. Получающиеся при этом
картинки, сходные с дифракционными, находят, как известно, совершенно ясные
объяснения вне де-бройлевского формализма. Наконец, стоит заметить, что
микрочастице вовсе не обязательно мигрировать сквозь кристалл, чтобы
провзаимодействовать с его периодическими полями, поскольку структура таких
полей посредством электромагнитных полей может быть "прощупана" еще на подлете к
поверхности монокристалла. Ни для кого же, в конце концов, не является секретом,
что силовое электромагнитное взаимодействие между частицами происходит
непрерывно и на любых расстояниях, хотя, возможно, и с разной интенсивностью в
зависимости от конкретных условий.
4.
Поскольку мерой механического движения, определяющей динамические свойства тел,
является импульс – вполне закономерно связывать дифракцию, явно или неявно, с
величиной . Каким бы образом падающий пучок не коллимировался и не
монохроматизировался, все равно приходится иметь дело со статистической физикой,
т.к. различные частицы испытывают рассеяние в существенно отличающихся случайных
условиях. Даже если принимать за основу некую локальную микроскопическую модель
рассеяния, то пришлось бы учитывать, например, возможности рассеяния вблизи
атомов и в междоузлиях или иные факторы влияния. Если же рассматривать процесс
рассеяния на всем кристалле сразу (модель большой молекулы), тогда обмен импульсом
с кристаллом может иметь разброс в широком диапазоне значений, определяемых
спектральным распределением электронов кристалла по импульсам. Таким образом,
приходится иметь дело с функциями распределения либо по координатам или углам,
либо по импульсам.
Из корреляционной теории, получившей существенные результаты
спустя 10-20 лет после работ де Бройля, Дэвиссона и Джермера, можно было бы
позаимствовать ряд полезных идей. Математический аппарат стохастических
процессов, основанный на Фурье-преобразованиях, это - канонические разложения случайных
функций в ряды Фурье, формулы Хинчина-Винера, спектральные разложения
корреляционных функций с дискретным или непрерывным спектрами [3].
Итак, попытаемся все-таки проанализировать сложный механизм,
управляющий дифракцией микрочастиц на монокристаллах. Как известно, явление
дифракции микрочастиц заключается в появлении резких дифракционных максимумов
при отражении пучка частиц от монокристаллов в случае выполнения условий
Вульфа-Брэгга (рис. 29.1):
2d
sin α = nλ D,
(29.1)
где d - постоянная
кристаллической решетки, α - угол скользящего падения, n - порядок дифракционного
максимума, λ D
= h/mv - длина волны де Бройля налетающих
частиц.
Ситуация может несколько проясниться, если приведенные условия
Вульфа-Брэгга переписать в несколько иной форме, а именно:
2p
sin a
=2pn = nh/d =
nћk,
(29.2)
где 2
pn - величина, на которую изменяется нормальная по отношению к
поверхности компонента импульса микрочастицы при зеркальном отражении от
поверхности монокристалла, nћk - дискретный спектр
нормальной компоненты импульса электронов кристаллической решетки,
k = 2π/d - модуль волнового вектора
обратной решетки. В таком виде
записи формулы Вульфа-Брэгга просматривается определенный резонанс при
обмене нормальной компонентой импульса между налетающей частицей и электронной
плазмой кристаллической решетки, т.е. действует закон сохранения проекции
суммарного импульса на нормаль к поверхности монокристалла. При этом проекция
импульса налетающей частицы на плоскость, совпадающую с гранью монокристалла,
может сохраниться автоматически при равенстве углов падения и отражения.
Однако угол падения может быть и не равен углу дифракции (рис.
29.2). При этом условия дифракции (29.1) несколько видоизменяются:
2
a (p2 cos a2 –
p1
cos a1 )
= nλ
D,
(29.3)
где
а – постоянная
кристаллической решетки в направлении, перпендикулярном поверхности кристалла.
При переходе к импульсам это же соотношение запишется в виде
p2
cos a2 –
p1
cos a1 =
nћk
(29.4)
Как мы видим,
это соотношение уже не похоже на закон Вульфа-Брэгга, а именно: здесь
отсутствует коэффициент 2, и угол падения не равен углу отражения.
В завершение сказанного следует отметить, что спектр компонент
импульсов электронов кристаллической решетки может быть найден методом простого
Фурье-анализа функций распределения электронов решетки по координатам и по
импульсам. Ранее Альфредом Ланде [4] также обращалось внимание на то, что
дифракция микрочастиц может быть объяснена без искусственной и неэкономной
гипотезы о волновой природе частиц, а определенными механическими свойствами
всей кристаллической решетки в целом, дискретной передачей импульсов телом,
имеющим периодическую структуру. Однако многие работы А. Ланде остались без должного
внимания.
§ 34. Стационарное уравнение Шредингера. Метод Фурье
Рассмотрим изолированную систему
атомов, которая не обменивается энергией с окружающей средой. Из классической
механики известно, что при движении замкнутой (консервативной) системы ее полная
энергия Е не меняется,
поэтому все точки в фазовом пространстве, изображающие состояние системы в
разные моменты времени, должны лежать на некоторой гиперповерхности, соответствующей
начальному значению энергии
Е. Уравнение этой поверхности в переменных p и q имеет
вид:
H(p,q)
= K(p) + U(q) = E,
(34.1)
где
H(p,q) – функция Гамильтона (или
гамильтониан), K(p) – кинетическая энергия,
зависящая от обобщенных импульсов, U(q) – потенциальная энергия,
зависящая от обобщенных координат.
В декартовых координатах закон сохранения полной энергии
E для отдельного электрона с
потенциальной энергией
U выглядит так:
p2/2m + U(x,y,z) = E,
(34.2)
где
p – импульс электрона, m – масса электрона. Полная
энергия E в среднем имеет одно и то
же значение в каждой точке траектории электрона. Воспользуемся этим
замечательным свойством энергии E для определения средней
электронной плотности в атомах.
Как уже отмечалось, в процессе длительного движения за счет
флуктуаций импульсов и координат электрон может побывать в самых неожиданных
точках пространства и в широком диапазоне
значений импульсов и кинетической энергии. Исходя из статистических
закономерностей, можно заранее сказать, что чем дальше точка находится от ядра,
особенно если речь идет о расстояниях r, значительно превышающих
средний радиус атома, тем с меньшей вероятностью можно встретить там электрон.
Другими словами, плотность вероятности w(x,y,z) пребывания электрона в
различных точках пространства или
функция распределения электронной плотности должна стремиться к нулю при r → ¥. Отсюда следует, что
функция распределения
w(x,y,z) для атома должна быть
абсолютно интегрируемой во всем пространстве, и для нее может быть введена
нормировка в виде (31.5).
Попробуем составить дифференциальное уравнение, из которого можно
было бы определить функцию w(x,y,z) с использованием всей
известной нам информации об атомах? в том числе и граничных условий для w(x,y,z). При этом мы учтем тот факт,
что импульс электрона в различных точках в атоме или молекуле может принимать не
произвольные значения, а на него накладывается ограничение при помощи
соотношения (34.2). Следовательно, при статистическом подходе можно рассматривать некоторое
пространственное распределение электронов по импульсам в соответствии с
выражением (33.2). Учитывая, что
импульс является вектором, в дальнейшем будем исследовать векторное поле p(x,y,z). При этом сразу отметим, что
функции распределения электронов по координатам и импульсам в атомах и молекулах
будут существенно отличаться от функций распределения, полученных Максвеллом и
Больцманом в молекулярной физике.
Характерно, что в статистике
Максвелла функция распределения по скоростям не зависит от координат, а зависит
от средней температуры газа, которая считается постоянной во всем объеме.
Примерно также обстоит дело с функцией распределения Больцмана, которая зависит
от координат, т.е. от потенциальной энергии и от температуры. Таким образом, обе
функции распределения считаются независимыми и входят в общую функцию
распределения по скоростям и по координатам частиц в виде произведения.
Естественно, что это является определенным приближением, поскольку средняя
кинетическая энергия частиц в потенциальном поле в различных точках пространства
обычно не является постоянной. Соотношение (34.2) накладывает ограничения на
допустимые значения импульсов частиц в потенциальном поле при заданной полной
энергии E и, следовательно, вносит
определенное уточнение в статистику электронов по сравнению со статистикой
Максвелла-Больцмана.
Таким образом, в рассматриваемой нами статистике электронов мы не
используем такого понятия, как температура частиц, которая была бы постоянной во
всем рассматриваемом объеме, а
учитываем тот факт, что кинетическая энергия электрона при его заданной полной
энергии E является функцией координат в
соответствии с уравнением (34.2).
Данная статистика более пригодна к внутриатомным движениям, где в
малых областях пространства с относительно малым количеством электронов имеются
очень сильные и неоднородные электромагнитные поля, и где становится невозможно
представлять распределения электронов по скоростям (или импульсам) и координатам
раздельно. Кроме этого, для отдельного электрона в атоме можно указать
определенные интегралы движения, такие, как полная энергия E, модуль полного момента
количества движения L и проекция этого момента на
ось симметрии Lz, чего нет в статистике
Максвелла - Больцмана за исключением полной энергии E.
Теперь рассмотрим более подробно некоторые специфические
особенности, присущие спектральному методу Фурье, в приложении к нашей задаче.
Состояние системы частиц описывается совокупностью 6N канонических переменных
q и p, подчиняющихся уравнениям
Гамильтона. Полагая совокупность канонических переменных случайной величиной,
как это имеет место при статистическом рассмотрении динамической микромодели,
каждому микросостоянию (q, p) сопоставляем некоторую
плотность вероятности, или функцию распределения w(q,p,t), которая статистически описывает движение
фазовых точек в фазовом пространстве.
Как уже было отмечено, согласно теореме Лиувилля движение фазового
ансамбля в фазовом пространстве можно по аналогии с гидродинамикой рассматривать
как движение несжимаемой фазовой жидкости. Это означает следующее: если
распределение случайной величины по координате w(x) сжимается в несколько раз
по оси x, то во столько же раз
распределение случайной величины по импульсам v(px) расширяется по оси
импульсов px,
сохраняя тем самым фазовый объем
Г. В силу независимости уравнений движения в проекциях то же
самое, очевидно, можно сказать и про остальные оси координат:
D
x D
px = const.
(34.3)
Аналогичный
результат можно получить и из теории адиабатических инвариантов. Смысл этой
теории заключается в том, что при адиабатических процессах в системе,
совершающей финитное движение, т.е. при медленном изменении некоторого
параметра l, характеризующего систему
или внешнее воздействие, сохраняется в среднем некоторая величина I,
называемая адиабатическим инвариантом. На языке формул это запишется
как
<dI/dt> = 0,
(34.4)
где I
обозначает
интеграл
(34.5)
вычисленный по траектории движения при заданных Е
и l.
Интегралу (34.5) может быть придан наглядный геометрический смысл,
если воспользоваться понятием фазовой траектории системы. В случае одной степени
свободы фазовое пространство сводится к двумерной системе координат (q, p), и
фазовая траектория системы, совершающей периодическое движение, представляет
собой замкнутую кривую в этой плоскости. Интеграл (34.5), взятый вдоль этой
кривой, представляет собой заключенную внутри нее площадь. Он может быть записан
и как двойной интеграл по двумерной плоской области
I =
(1/2p)
òò
dpdq,
(34.6)
т.е. по
аналогии с теоремой Лиувилля (31.1).
Отсюда можно заключить, что переменные q и p при движении системы в
фазовом пространстве ведут себя по аналогии с сопряженными переменными в
преобразованиях Фурье и к ним можно применить спектральный метод Фурье. Другими
словами, функция распределения системы частиц по импульсам будет вести себя
подобно спектру Фурье по пространственным частотам. При этом импульс рх с точностью до некоторого
размерного коэффициента может служить аналогом пространственной частоты fx для сопряженной переменной х.
По аналогии с функцией распределения электронов по координатам
w(x,y,z) вида (31.2) рассмотрим
функцию распределения электронов по импульсам
v(px, py, pz ), для которой справедливо
соотношение
v(px) = çj(px)ç2,
(34.7)
где j (px) – комплексная амплитуда
функции распределения по импульсам.
В соответствии с методом Фурье частотный спектр c(fx) функции Ф(x) можно найти с помощью
интеграла
¥
c(fx) =
òФ(x)
exp(-2p i
fx x)dx.
(34.8)
-¥
Для обратного преобразования Фурье, как известно,
имеем
¥
Ф(x) =
ò
c(fx) exp(2p i
fx x) d fx.
(34.9)
-¥
при
этом выполняется равенство Парсеваля
ò
çФ(x)ç2
dx =
ò
çc(fx)ç2
d fx.
(34.10)
-¥
-¥
Для функций распределения
w(x) и v(px) имеется аналогичное
соотношение нормировки
¥
¥
ò
çФ(x)ç2dx =
ò
çj(рx)ç2
dрx = 1.
(34.11)
-¥
-¥
Продолжая аналогию для сопряженных переменных х и
рх, мы можем
записать
¥
j(рx) =
òФ(x)
exp(-2p i
рxx/h)dx/Ö
h
и
(34.12)
-¥
¥
Ф(х) = ò
j(рx)
exp(2p i
рxx/h)dрх/Ö
h,
(34.13)
-¥
где h = 2p ћ – размерный коэффициент, определяемый из
эксперимента. В соотношениях (34.12) и (34.13) хорошо видно, что переменная
px / h играет такую же роль, что и
пространственная частота fx в (34.8) и (34.9), т.е.
px аналогично hfx = ћ kx. Теперь становятся понятными
соотношения, очень часто используемые в атомной физике как
постулат,
px =
ћ
kx , p = ћ k.
(34.14)
Согласно методу Фурье
динамическую переменную в уравнении можно заменить производной от комплексной
амплитуды по сопряженной переменной. Например, всем хорошо знакомы
дифференциальные уравнения для частных производных по времени в методе
Фурье
¶
Ф/ ¶
t = iw
Ф =
2p
if Ф; ¶
2
Ф/ ¶
t 2 =( iw)
2 Ф =
- w
2
Ф,
(34.15)
где
частота w выступает в роли сопряженной переменной
для t.
По аналогии с (34.15) для переменной x
в методе Фурье справедливы уравнения
¶
Ф/ ¶ x
= 2p i
f x Ф
;
¶
2Ф/¶
x 2 = - (2p
f x) 2 Ф
,
(34.16)
Теперь в соответствии с (34.14) заменяем пространственную частоту
fx в (34.16) через
импульс px по формуле
px = ћ kx
= ћ 2 π fx .
(34.17)
В итоге
получаем
¶
Ф/¶ x
= i p xФ/
ћ;
-i ћ
¶
Ф/¶ x
= p xФ;
- ћ
2¶
2Ф/¶
x 2 = px 2 Ф
,
(34.18)
С помощью соотношений (34.18) задача на нахождение средних
значений импульсов может быть решена полностью в конфигурационном пространстве,
т.е. минуя переменные импульсов рх.
С использованием уравнений (34.18) закон сохранения полной энергии
электрона в атоме
p2/2m +
U(x,y,z,) = E
(34.19)
запишется для трех координат в виде
-
(ћ2/2m)
DФ(x,y,z)
+ U(x,y,z) Ф(x,y,z)
= EФ(x,y,z).
(34.20)
Таким образом, мы получаем
стационарное уравнение Шредингера,
опираясь только на теорему Лиувилля
для функций распределения электронов по координатам и по импульсам,
спектральный метод Фурье и закон
сохранения полной энергии в атомных системах.[2]
При выводе статистического уравнения Шредингера (34.20) нами были
использованы наиболее общие свойства движения электронов в атомных системах,
т.е. движения в центральных кулоновских полях ядер, поэтому оно может быть
справедливо для общего анализа различных атомных, молекулярных и других систем,
где не учитываются разного рода релятивистские эффекты. С помощью данного
уравнения мы можем отыскать функцию Ф(x,y,z) как комплексную амплитуду
функции распределения электронов w(x,y,z) для конкретного вида
функции потенциальной энергии
U(x,y,z) и заданной полной
энергии E электрона.
Необходимо подчеркнуть еще одну особенность данной задачи. Условие
нормировки (31.5) требует, чтобы функция Ф на бесконечности обращалась в
нуль. Подобные ограничения на функцию называются граничными условиями, а задача
с такими условиями называется краевой. Как правило, решение подобных задач
сводится к отысканию собственных функций и собственных значений некоторого
параметра уравнения, удовлетворяющего частному решению линейного
дифференциального уравнения (задача Штурма - Лиувилля). В качестве собственных
значений в данной задаче могут выступать различные интегралы движения электронов
в центральных полях (полная энергия E, модуль момента количества
движения L, а также проекция
механического момента на ось симметрии системы Lx).
§ 36. Связь
энергии с частотой. Динамическое уравнение Шредингера
Рассмотрим
систему, состоящую из большого
числа атомов, т.е. вещество. Распределение электронной плотности в веществе при
стационарном движении электронов описывается стационарной функцией распределения
w(x,y,z), уравнение для которой мы
вывели в § 31. Отметим, что под стационарным распределением следует понимать
усредненную за бесконечное время наблюдения плотность электронной плазмы в веществе. Однако из-за
статистического характера движения электронов электронная плотность не остается
постоянной, а непрерывно флуктуирует во времени за счет естественного
орбитального движения электронов в атомах и молекулах. Таким образом, для
установления более полной динамической картины поведения электронной плазмы
необходимо учитывать временную зависимость функции распределения. При этом
функция распределения будет иметь вид: w(x,y,z,t). Требуется определить
характер этой зависимости и, в частности, зависимость от времени комплексной
амплитуды функции распределения
Ф(x,y,z,t).
С этой целью воспользуемся законом сохранения полного заряда или
числа частиц в веществе. Этот закон в дифференциальной форме известен как
уравнение непрерывности для плотности заряда r и плотности электрического тока j = vr, где v - средняя скорость электронов,
¶r /¶ t + div j = 0. (36.1)
Используя статистический метод описания, выразим плотность
электрического заряда через функцию распределения
электронной плотности w(x,y,z,t) и заряд электрона e
r(x,y,z,t) = ew(x,y,z,t) = e |Ф(x,y,z,t)|2.
(36.2)
Подставив это выражение в
(36.1) и сократив заряд e, получим уравнение
непрерывности для функции
распределения электронной плотности
¶
/¶
t |Ф|2
+ div |Ф|2
v =
0.
(36.3)
Данное уравнение на языке
статистической механики может быть
интерпретировано следующим образом. Первый член уравнения означает изменение
функции распределения или электронной плотности во времени в данной точке
пространства. Второй член имеет смысл потока функции распределения или потока
плотности вероятности через малую сферу, окружающую данную точку, в соответствии
с определением дивергенции вектора. Вполне естественно, что от функции
распределения некоторой физической
величины (заряда, массы, энергии и
т.д.) можно всегда перейти к описанию поведения во времени самой
физической величины в терминах механики сплошных сред.
Таким образом, использование той или иной функции распределения
является мостиком или связующим звеном между описанием движения дискретных
объектов в статистической механике и в механике непрерывных сред, которые всегда
являются некоторой идеализацией реального вещества, состоящего из атомов и
молекул. Недооценка этого подхода породила в квантовой теории представление об
отдельной частице как о протяженном (размытом) в пространстве объекте: например, волна де Бройля, волновой
пакет или электрон в виде облака, хотя речь идет, как правило, всего лишь о
функции распределения, то есть плотности вероятности местонахождения
частицы в заданном объеме.
С целью определения временной зависимости Ф(x,y,z,t) приведем уравнение (36.3) к виду
¶
/¶
t (Ф*Ф) = -
div(Ф*Ф v) .
(36.4)
Преобразуем правую часть (36.4) к симметричной форме с заменой
скорости v на импульс электрона p = mv
¶ /¶ t(ФФ*) = - (1/2m)div(Ф*рФ + ФрФ*).
(36.5)
Теперь воспользуемся методом
Фурье, применявшимся нами для вывода стационарного уравнения Шредингера (34.18),
и произведем замену импульса p в (36.5) дифференциальным
оператором[3] в соответствии с
выражением
рФ =
- iћ
ÑФ.
(36.6)
Подставив (36.6) в (36.5) и выполнив
дифференцирование, получаем
¶ /¶ t(ФФ*) = (iћ/2m)(Ф*DФ + ФDФ*).
(36.7)
Для
приведения уравнения к удобному виду с целью совместного его решения со
стационарным уравнением Шредингера (34.18) умножим все члены (36.7) на
i ћ, раскроем скобки в левой
части, а в правую часть прибавим и вычтем слагаемое
Ф*UФ. После соответствующей
перегруппировки слагаемых в получившемся уравнении имеем
iћФ*¶
Ф/¶ t
– [iћФ*¶
Ф/¶
t]* = - (ћ2/2m)Ф*DФ +
Ф*UФ –
[- (ћ2/2m)Ф*DФ + Ф*UФ]*.
(36.8)
Данное сложное
уравнение состоит из двух комплексно сопряженных частей более простого
уравнения
iћ
¶Ф/¶
t = - (ћ2/2m)DФ +
UФ,
(36.9)
где
произведено сокращение на
Ф*.
Мы получили полностью классическим путем динамическое уравнение
Шредингера, которое было введено им также в виде постулата, исходя из волновых
соображений в 1926 году. Если потенциальная энергия U(x,y,z) не зависит от времени, то в
уравнении (36.9) можно произвести разделение переменных, записав функцию
Ф в виде произведения временной и пространственной
функций
Ф(x,y,z,t) =
В(t)П(x,y,z)
.
(36.10)
Подставим это выражение в уравнение (36.9)
iћ
П
dB/dt = -( ћ2/2m)
ВDП +
U ВП
(36.11)
и
поделим обе части уравнения (36.11) на ВП
iћ
1/В
dB/dt = - ћ2/2m
(DП)/П +
U.
(36.12)
Хорошо видно, что левая часть уравнения (36.12) зависит только от
времени, а правая часть - только от координат. Для любых значений t и координат это возможно только в том
случае, если обе части уравнения равны некоторой константе Е, называемой параметром разделения. В
результате мы получаем два независимых уравнения
iћ
dB/dt = E B
(36.13)
и
(- ћ2/2m) DП + U П = E
П.
(36.14)
Второе
уравнение является стационарным уравнением Шредингера, которое было выведено
нами ранее методом Фурье. Уравнение (36.13) легко решается и приводит к
следующей временной зависимости функции В(t):
В(t)
= exp(-iEt/ћ)
= exp(-iw
t),
(36.15)
где
w = E/ћ -
боровская частота. Данное уравнение могло быть получено и на основе
гармонического Фурье-анализа из соотношения
¶ Ф/¶ t = iw Ф.
(36.16)
Проанализируем
различные варианты решения уравнений (36.13) и (36.14). Если параметр разделения
Е является комплексным, то
функция В(t) будет изменяться во времени
по модулю, что приведет к изменению электронной плотности в данной области
пространства, т.е. соответствует нестационарному процессу, а полная энергия не
является интегралом движения системы. При этом будет происходить либо
поглощение, либо излучение электромагнитной энергии в зависимости от
характеристик системы. В том же случае, если параметр разделения является
действительной величиной, модуль функции В(t) становится постоянным и зависимость В(t) является гармонической с частотой
колебаний w. Это
соответствует стационарному движению электронов в атомах и молекулах, т.е.
параметры движения в системе остаются постоянными.
§ 36. Функции распределения электронов по импульсам в
периодических структурах
Рассмотрим поведение электронов в периодической атомной структуре.
Предположим, что атомы расположены в упорядоченной линейной цепочке с
бесконечной протяженностью и расстояния между центрами атомов равны
a (рис. 36.1). Подобные
цепочки атомов можно наблюдать в достаточно крупных совершенных монокристаллах,
состоящих из однородных атомов.
Допустим также, что электроны в атомах цепочки движутся по
некоторым стационарным орбитам, в результате чего устанавливается равновесное
распределение электронов в пространстве. Естественно ожидать, что в данном
случае распределение электронных облаков имеет пространственную периодичность с
периодом, равным
а.
Таким образом, как функция распределения электронов по координатам
w(x,y,z), так и комплексная
амплитуда функции распределения Ф(x,y,z) для цепочки атомов имеют
пространственную периодичность с периодом
а.
Разложим периодическую функцию Ф(x,y,z) в ряд Фурье в комплексной
форме по пространственным гармоникам по оси x, полагая при этом, что по
осям y и z предварительно произведено усреднение
электронной плотности,
Ф(x)
= Cn
exp (inkx),
(36.1)
где
k = 2π/a – пространственная частота
первой гармоники, Cn – коэффициенты разложения
или амплитуды Фурье-гармоник. Cn ряда (36.1) можно определить посредством
интегрирования на этом отрезке
a/2
Сn
=
1/a ò
Ф(x)exp(-inkx)dx.
(36.2)
-a/2
Достаточным условием применимости
разложения Фурье являются некоторые математические ограничения, сводящиеся к
тому, что функция Ф(x) должна быть непрерывной
или иметь конечное число разрывов, а также в пределах одного периода иметь
конечное число максимумов и минимумов (условие Дирихле). Вполне естественно, что
большинство функций распределения удовлетворяют этим условиям.
Можно предположить, что распределение электронной плотности в
реальных периодических структурах имеет сравнительно гладкий характер, т.е. не
имеет очень резких перепадов, как и большинство других распределений в
статистической физике. В таком случае ряд Фурье (36.1) будет быстро спадать с
увеличением n, т.е. не будет иметь в своем
составе очень высокочастотных пространственных гармоник.
Хорошо известно, что слагаемые ряда Фурье являются взаимно
ортогональными функциями при интегрировании их по оси x. Таким образом, ряд (36.1)
можно рассматривать как разложение функции Ф(x) по собственным
ортогональным функциям Фn(x), причем квадрат модуля
коэффициента |Сn|2
означает спектральный вклад n-й гармоники в общий частотный спектр.
Существенно отметить, что согласно (36.1) функция распределения
электронов по координатам в периодической структуре имеет дискретный частотный
спектр. Можно ожидать, что и функция распределения электронов по импульсам также
будет дискретной.
Подставив значение функции Ф(x) из ряда Фурье (36.1) в
выражение функции j (px) (33.12), мы можем вычислить
ее явный вид
(36.3)
Воспользовавшись свойствами дельта – функции
d(k)
= (1/2p
) exp(ikx)dx
(36.4)
и
вычислив интеграл (36.3) по х,
получаем
(36.5)
Мы
действительно получили дискретный спектр функции j (px), причем импульсы
электронов px могут принимать здесь лишь значения n ћ kx
pxn = nћk, (n = 1,2,…,¥).
(36.6)
Отсюда можно заключить, что
такая электронная система в периодической структуре может обмениваться лишь
дискретными импульсами, кратными
ћ kx, что и объясняет, в
конечном итоге, дифракцию любых частиц на монокристаллах, а также условия Вульфа
– Брэгга (29.1)-(29.2) и соотношения (29.3) – (29.4). Более подробно обо всем
этом можно прочесть в монографии [5].
ЛИТЕРАТУРА
1.
Блохинцев Д.И. Основы квантовой механики. – М.: Наука, 1983.
2.
Киттель Ч. Квантовая теория твердых тел. – М.: Наука, 1967.
3.
Вентцель Е.С. Теория вероятностей. – М.: Наука, 1962.
4. Lande A. The case against quantum dualiti / "Philos.
Sci.",1962. V.29, N.1, p.1-6.
5. Шаляпин А.Л., Стукалов В.И.
Введение в классическую электродинамику и атомную физику. Второе издание,
переработанное и дополненное. Екатеринбург, Изд-во Учебно-метод. Центр УПИ, 2006, 490
с.
http://s1836.narod.ru
[1] По Интернету как-то приключился ехидный
запрос: а нейтронография, а нейтроны? На встречный вопрос: а магнитный момент, а
внутренняя зарядовая структура нейтронов? – больше вопросов не было.
Действительно, если заряд мы еще только начинаем постигать, то что мы можем
толковать о структуре элементарных частиц. Складывается такое впечатление, что
де Бройль поначалу пытался представить себе электрон в виде бильярдного шара,
который с треском отскакивает от кристалла. Не получилось. Ну, тогда – волна.
Если есть желание выстраивать образные ассоциации, то электрон – это, скорее,
шарик для пинг-понга, налетающий на небоскребы Манхеттена. Обмен импульсом имеет
место, но только посредством полей или волн. Существенным фактом является то,
что обмен осуществляется дискретным образом, и в этом “виноват” гармонический
анализ Фурье, а точнее – периодическая структура кристалла, которую “видит”
электрон. Мы часто употребляем словосочетание – акустическая модель физического
вакуума (эфира). Так что электрон, скорее всего, “слышит” все эти строго
регулярные “кварталы большого города”.
[2] Это уравнение используется для вычисления
стационарной электронной плотности, а также целого ряда других характеристик в
атомах и молекулах. Уравнение (34.20) было постулировано Шредингером в 1926 г.
на основе
корпускулярно-волновых представлений де Бройля и до последнего времени не
могло быть выведено в рамках квантовой теории.
[3] В принципе, по поводу линейных операторов
можно было бы вообще воздержаться от комментария, т.к. данный математический
объект знаком любому студенту, начиная с первого семестра, по линейной алгебре,
а со второго семестра – по теории дифференциальных уравнений. Более того, можно
было бы записать все уравнения, так ни разу и не прибегнув к алгебре операторов,
если бы не одно обстоятельство… Мы имеем в виду тщательно насаждаемое мнение,
что «математический аппарат квантовой механики резко отличается от
математического аппарата классической механики». И даже более того, оказывается,
«что все взаимоотношения между механическими величинами в квантовой области
могут быть выражены на языке операторов определенного класса. В этом заключается
фундаментальное значение введения операторов в квантовую механику».
Поразительно, но в глазах (или умах?) студентов прививается четкий стереотип:
нет операторов – это “классика”, есть операторы – это “кванты”.