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

КВАНТЫ билеты / 16. Вариационный метод

.pdf
Скачиваний:
74
Добавлен:
23.02.2015
Размер:
367.36 Кб
Скачать

Вариационный метод

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

ˆ

 

 

 

 

 

 

 

 

 

 

 

ˆ

H . Для упрощения дальнейших рассуждений считаем, что спектр

H не

вырожден и дискретен. Следовательно, необходимо решить уравнение

 

 

 

 

 

ˆ

 

En

= En

 

En

 

(7.59)

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

на собственные векторы

 

En

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

и собственные значения En гамильтониана H .

Если при

отыскании

 

En

и

En

процедура диагонализации

гамильтониана

ˆ

 

 

 

 

причинам

затруднительна или

нельзя

H по каким-либо

 

 

 

 

 

 

 

 

 

 

 

ˆ

ˆ

использовать стационарную теорию возмущений (например, когда V

~ H 0 ),

то представляет интерес применение вариационного метода. Основная суть его заключается в следующем. Вместо задачи, связанной с решением уравнения (7.59), рассматривается эквивалентная задача по поиску условного

экстремума среднего значения E g

 

 

 

 

ˆ

 

g

 

 

оператора H в состоянии

 

 

 

 

 

 

E g = g

 

ˆ

 

g

, где

 

g

 

g 1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

как

функционала

 

от вектора состояния

 

g . Обозначим

 

функционал

 

 

 

 

 

 

 

 

 

посредством E

g

 

. Напомним, что под функционалом понимают некоторое

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

правило, сопоставляющее множеству векторов

 

g числовое множество

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

g .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В вариационном исчислении существует теорема, которую применительно к нашему случаю можно сформулировать так: функционал

 

 

 

 

 

 

E

g

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

E

g

=

E

g g E

,

 

 

 

 

 

 

 

 

 

 

 

обусловленная заданием вариации g вектора g , и вычисленная с

точностью до членов первого порядка малости (ее называют первой вариацией)

 

 

 

 

 

 

 

 

 

E

 

g

E

g

=

 

 

 

 

g ,

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

равна нулю.

 

 

 

 

 

 

В этой связи условие на экстремум функционала

E

g будет

 

 

 

 

 

аналогично вариационному уравнению

 

 

 

 

 

 

E

g = 0 .

 

(7.60)

 

 

 

 

 

Таким образом, суть вариационного метода состоит в поиске решений уравнения (7.60).

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

1. Вначале на основе физических критериев и качественного анализа изучаемой квантовой системы выбирается для n-го состояния некоторый вектор состояния g( , , ...) , зависящий от нескольких непрерывных

вариационных параметров , , ... . Такой вектор называется пробным.

2.С пробным вектором g( , , ...) производят вычисления

 

 

 

 

 

 

 

 

 

 

функционала E

g( , , ...)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

g( , , ...)

=

g( , , ...)

 

ˆ

 

 

g( , , ...) = E( , , ...) .

(7.61)

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При этом

E

g( , , ...)

 

сводится

к

обычной функции

E( , , ...) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

аргументом которой являются вариационные параметры.

3. Затем выполняется поиск значений неизвестных вариационных параметров, при которых имеет место экстремальность функционала

 

 

 

 

 

E

g( , , ...)

. Поиск сводится

к

отысканию

первых производных по

 

 

 

 

 

 

 

 

параметрам , , ... от функции

E( , , ...) ,

т.

е. к решению

системы

алгебраических уравнений

 

 

 

 

 

 

 

E( , , ...)

 

E( , , ...)

…= 0 .

(7.62)

 

 

 

 

 

 

 

 

Набор параметров , , ..., полученный из этой системы уравнений,

определяет

решение

 

 

g( , , ...)

 

уравнения

(7.60) и

величину

 

 

 

E

 

g( , ,

...) . Причем

вектор

 

g( , , ...)

приближенно

отвечает

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

собственному вектору

 

E

 

 

оператора

 

ˆ

 

g( , , ...)

с более

 

 

 

 

 

 

 

 

H , а величина E

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

лучшей точностью – его собственному значению En .

 

 

 

 

 

Для энергии E0 основного состояния выполняется неравенство

 

 

 

 

 

 

 

 

 

E

0

E

 

g( , , ...) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Равенство имеет место

 

только в том

случае,

когда

 

 

g( , , ...) является

 

 

 

собственным вектором

 

 

E0

оператора Гамильтона

 

 

ˆ

 

 

 

 

 

 

 

 

H соответствующем

собственному значению E0

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

E

 

= E

= E

 

ˆ

 

 

E

= E

 

.

 

 

 

 

 

 

H

 

 

 

 

 

0

 

 

Е0

0

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

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

в первом порядке теории возмущения, которые определяются с помощью формулы

 

 

En =

0

ˆ

0

+

0

 

ˆ

0

 

 

 

0

 

 

 

En | H0

| En

En

V

En

= En +Vnn ,

 

эквивалентны

вычислениям

вариационным методом с пробным вектором

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

ˆ

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

En ( , ,...) , имеющим вид собственного вектора

En

оператора H0 .

 

Для иллюстрации применения вариационного метода рассмотрим

одномерный гармонический осциллятор

(см.

п. 4.3). Для него

оператор

 

ˆ

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гамильтона H

= H0 , записанный в координатном представлении, имеет вид

 

 

 

ˆ ( x)

 

2

 

d 2

 

 

 

m 2 xˆ

2

.

 

(7.63)

 

 

 

H

 

2m dx2

 

 

2

 

 

 

 

 

 

Поставим

задачу

по

определению

 

 

 

полной

энергии Е0

основного

состояния осциллятора и 0 (x) -функции, соответствующей этому состоянию. При выборе пробной функции учтем, что она, во-первых, должна

обращаться в нуль при стремлении x

к бесконечности: x ,

во-вторых,

 

 

 

 

 

 

 

 

ˆ

 

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

должна иметь узлов.

 

 

 

 

 

 

 

 

 

Перечисленным требованиям

удовлетворяет следующая

функция

0 (x, ) , зависящая от одного параметра

 

 

0

(x, ) Ce x2 /2 , 0 .

(7.64)

 

 

 

 

 

 

 

 

 

Из условия ее нормировки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 (x, )

 

2 dx C2 e x2 dx 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определяется множитель

 

 

 

 

1/4

 

 

C

.

 

 

 

 

 

 

 

 

 

 

С помощью (7.63) и (7.64) вычисляется функционал E 0 (x, )

 

 

 

1

 

2

m 2

1

 

 

E 0 (x, ) =

ˆ

(x, )dx

 

 

 

 

 

 

 

 

= E( ) . (7.65)

 

 

 

 

 

0 (x, )H 0

 

 

 

 

 

 

 

 

 

 

2

 

2m

2

 

 

 

Следовательно, функционал E 0 (x, ) от функции

0 (x, ) приобретает

вид функции E( ) от параметра . Поэтому задача состоит в определении

минимума этой функции. Первая производная функции E( )

по

обращается в нуль, когда =

 

 

 

 

 

=

m

.

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя это значение в (7.65), находят энергию E0

основного состояния

 

 

 

.

 

(7.66)

 

 

 

E0 = E( ) =

2

 

 

 

 

 

 

Соответствующая основному состоянию 0 (x) -функция, согласно (7.64), имеет вид

 

m 1/ 4

 

m x2

/ 2

 

 

0 (x) 0 (x, )

 

 

 

e

 

 

.

(7.67)

 

 

 

 

 

 

 

 

 

 

 

Отметим, что из-за простоты рассмотренной задачи полученные

вариационным методом значения энергии

E0

основного состояния и вид

функции 0 (x) в точности совпадают с аналитическими решениям стационарного уравнения Шредингера (4.19).