
КВАНТЫ билеты / 16. Вариационный метод
.pdf
Вариационный метод
Для консервативных квантово-механических систем приближенные решения стационарного уравнения Шредингера в некоторых случаях удобно искать с помощью вариационного метода, к изучению которого сейчас приступим. Предположим, что имеется некая квантовая система, физическую модель которой характеризует независящий от времени оператор Гамильтона
ˆ |
|
|
|
|
|
|
|
|
|
|
|
ˆ |
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).