Линейным гармоническим осциллятором называется система, потенциальная энергия которой квадратично зависит от координаты:

Здесь m — масса частицы, а ω — собственная частота осциллятора. На рис. 11.1 зависимость (1)

изображена графически. Кривая U(x) своей крутизной и бесконечно большой высотой напоминает потенциальную яму. Ниже мы увидим, что линейный осциллятор, действительно, проявляет некоторые свойства частицы в бесконечно высокой потенциальной яме. Например, он имеет бесконечное число дискретных уровней. Но в отличие от отвесных стенок ямы, потенциал осциллятора растёт плавно, и, как следствие, появляется некоторая вероятность обнаружить частицу достаточно далеко от начала координат. Плавная форма потенциала позволяет осциллятору при определённых условиях проявить свойства классической (не квантовой) частицы. Для этого достаточно, чтобы длина волны де Бройля была меньше характерных размеров области изменения потенциала. В случае потенциальной ямы, либо потенциального барьера, такая возможность полностью исключена, так как там потенциал меняется скачком в одной точке. Перейдём к количественному решению задачи.

Напишем одномерное уравнение Шредингера с потенциальной энергией (1):

У него нет естественных граничных условий. Дискретные уровни энергии получаются как следствие ограниченности волновой функции.

Преобразуем уравнение (2): вместо координаты x введём безразмерный аргумент

и вместо E — безразмерную энергию осциллятора

Легко убедиться, что обратная величина подкоренного выражения в (3) равна произведению комптоновской длины волны электрона  (1.3.1) и длины волны D = l/(2p) =c/w, соответствующей собственной частоте осциллятора. С принятыми обозначениями уравнение Шредингера принимает вид:

Здесь штрихами обозначено дифференцирование по координате y.

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

11.1 Решение волнового уравнения

Уравнение (5) решаем методом разложения в ряд с предварительным выделением множителя, быстро убывающего на больших расстояниях от начала координат. Обозначим посредством ψ_∞ волновую функцию при больших значениях аргумента: y ? ε. Ей соответствует асимптотическая форма уравнения Шредингера

Множитель  перед функцией появляется при двукратном дифференцировании экспоненты , а именно:

.

В области больших значений аргумента можно пренебречь единицей по сравнению с y² в множителе . Следовательно, рассматриваемая экспонента может служить асимптотическим приближением точного решения:

Решение уравнения (5) при произвольных значениях аргумента будем искать в виде

Согласно (5), функция f(y) удовлетворяет уравнению

Ищем решение в виде ряда:

Дважды продифференцируем искомую функцию

Подставив эти разложения в (1.3), получим бесконечную цепь линейных уравнений для коэффициентов разложения A_ν:

Соберем степени с одинаковыми показателями, для чего индекс первой суммы увеличиваем на два:

Сумма тождественно равна нулю, когда исчезает каждый коэффициент при y^ν. Отсюда вытекает рекуррентное соотношение для коэффициентов разложения (1.4):

Задав A₀ и A₁, мы получим все коэффициенты A_n, соответственно, с чётными и с нечётными номерами. Формальное решение задачи получено.

Теперь покажем, что условию ограниченности волновой функции удовлетворяет только конечная сумма, но не бесконечный ряд. Для этого покажем, что ряд, коэффициенты которого подчиняются условию (1.5), растёт быстрее, чем exp( y²/2). Разложим экспоненту в ряд Тейлора:

Здесь ν принимает только чётные значения. Мы ввели обозначение

Отношение

при больших значениях ν стремится к 1/ν. В то же время из (1.5) следует:

Теперь ясно, что функция, описываемая соотношением (1.5), растёт быстрее, чем экспонента (1.6). Следовательно, произведение (1.4) при больших значениях аргумента неограниченно возрастает. Поэтому физический смысл имеет только такое решение (1.3), в котором сумма (1.4) содержит лишь конечное число слагаемых.

Покажем, как требование конечности числа слагаемых приводит к дискретному спектру энергетических уровней осциллятора. Пусть n — номер последнего члена ряда (1.4), не равного нулю:

Из (1.5) вытекает связь между величиной энергетического уровня и его номером:

Итак, мы снова получили дискретные уровни энергии. Каждому уровню с номером n соответствует ровно одна волновая функция:

где

причём нижний предел равен нулю при чётном n и единице — при нечётном. Уровни энергии, согласно (4) и (1.7), принимают дискретный ряд значений:

Обратим внимание на то, что наименьшее из возможных значений энергии, равное

отлично от нуля. Это состояние соответствует так называемым нулевым колебаниям.

Рекуррентное соотношение (1.5) с учётом (1.7) принимает вид

Оно позволяет восстановить все коэффициенты суммы (1.8а) с точностью до общего множителя; последний может быть найден из условия нормировки волновой функции (1.8). Поскольку ν < n, то сумма (1.8а) — знакопеременная. Функция f_n является чётной при чётных значениях её номера n, и нечётной — при нечётных n. Она пропорциональна известным в математике полиномам Эрмита H_n(y). Выпишем для справки нормированную волновую функцию состояния осциллятора с номером n:

Приведём несколько первых полиномов Эрмита:

В общем случае

Исследуем полученное решение.

11.2 Свойства квантового осциллятора

Линейный осциллятор имеет эквидистантную систему уровней: разность энергий двух соседних уровней постоянна и равна ħω. Именно такой квант энергии излучается или поглощается при переходе между соседними уровнями. На рис.11.2.1 слева приведены графики волновой функции для трёх первых значений n. По горизонтальной оси отложены значения безразмерного аргумента y.

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

Волновые функции линейного осциллятора и рассмотренной выше задачи о прямоугольной потенциальной яме имеют некоторые сходные черты. Во–первых, у них одинаково количество узлов. Здесь мы ещё раз встречаемся с проявлением осцилляционной теоремы, упомянутой в девятой главе. Так, волновая функция основного состояния ψ₀ не обращается в нуль ни в одной точке на прямой, а графики функций ψ₁ и ψ₂ пересекают горизонтальную ость, соответственно, один и два раза. Для остальных состояний имеет место тот же самый результат. Действительно, из теории специальных функций известно, что полином Эрмита n–го порядка H_n имеет ровно n корней. Во–вторых, квантовую частицу можно обнаружить и в области, запрещённой для движения в классической механике, что иллюстрирует график на рис.11.2.1 справа.

В состояниях с сильным возбуждением квантовый осциллятор приобретает свойства классической частицы. На рис. 11.2.2 схематически изображён график вероятности в

пределе n » 1. Пунктиром, как и выше, отмечена потенциальная кривая, а синим цветом — вероятность обнаружения. Функция W(x) быстро осциллирует, причём её огибающая (она обозначена на рисунке красным цветом) монотонно растёт от центра к периферии. Такой рост имеет аналогию в классической механике, когда вероятность обнаружения частицы на отрезке длиной Δx обратно пропорциональна её скорости V:

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

x = x₀ sin ωt.

Классическую вероятность обнаружения частицы на отрезке от x до x + dx определим как отношение времени прохождения отрезка

к половине периода колебаний:

Вероятность значительно увеличивается вблизи точек поворота x = ± x₀. Для достаточно больших значений n в апертуру прибора с конечной разрешающей способностью попадает много горбов, поэтому осцилляции волновой функции в классическом пределе незаметны. Но для основного состояния квантовая теория и классическая механика дают принципиально разные ответы, что иллюстрирует рис.11.2.3. Зелёная линия на нём

обозначает уровень нулевых колебаний E₀ из (1.10). Классическая вероятность увеличивается по мере приближения к точке x₀, определяемой условием U(x₀) = E₀, но не может перейти через эту границу. Квантовая теория предсказывает уменьшение вероятности при приближении к границе, причём частица может быть обнаружена в классически недоступной области x > x₀.

11.3. Нулевые колебания

В самом низком состоянии осциллятор имеет отличную от нуля энергию (1.10), определяемую его собственной частотой. Нулевые колебания осциллятора имеют чисто квантовую природу и находят своё объяснение в соотношении неопределённостей. Полная энергия осциллятора равна сумме кинетической и потенциальной

^.

Если частица локализована внутри области размером x, то, согласно принципу Гайзенберга, её импульс не может быть меньше, чем

^.

Таким образом, первое слагаемое в (3.1) уменьшается по мере увеличения x, а второе — растёт:

^.

Полная энергия как функция x имеет минимум в точке

^,

равный

^.

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

^,
Формально она выводится следующим образом. Пусть состояние частицы описывается функцией y(x), причём средние значения импульса и координаты для простоты вывода предполагаются равными нулю. Напишем очевидное неравенство



где a — действительная постоянная. Далее вычислим три полезных выражения:

Они позволяют преобразовать левую часть (3.3) в квадратичный по a трёхчлен:



Этот трёхчлен не меняет знака ни при каких значениях a, если его дискриминант отрицателен, либо равен нулю. Отсюда вытекает условие

Легко видеть, что оно тождественно (3.2).