Формула Гельфанда-Яглома

Список литературы:

  1. R. Rattazzi, "The Path Integral approach to Quantum Mechanics", глава 3.1.
  2. R. Rosenfelder, "Path Integrals in Quantum Physics", глава 1.3.

Мотивация: квазиклассическое приближение

Рассмотрим амплитуду перехода частицы из точки xi в точку xf за время T. На языке функционального интеграла мы можем записать:

(1)U(xi,xf;T)=DxeiS[x]/.

Сделаем замену переменных – перейдем к интегрированию по флуктуациям y относительно классического пути xcl, то есть x=xcl+y. В этом случае находим:

(2)U(xi,xf;T)=eiS[xcl]/Dyexp(i[12δ2Sδx2y2+16δ3Sδx3y3+]).

Теперь перемасштабируем величину y=y~ и тогда:

(3)U(xi,xf;T)=NeiS[xcl]/Dy~exp(i12δ2Sδx2y~2+O(y~3)).

В пределе 0 мы можем ограничиться квадратичным членом и таким способом должны построить квазиклассическое приближение. Логика этого построения ничем не отличается от логики взятия интегралов методом перевала.

Напомню, что для действия вида S[x]=tatbdt[mx˙2/2V(x)], описывающего движение частицы в произвольном потенциале V(x), мы ранее получали (см. семинар 3, упражнение 1):

(4)S[xcl+y]=S[xcl]+tatbdt[m2y˙2(t)12V(xcl(t))y2(t)]+.

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

(5)ω2(t)=V(xcl(t))/m.

Представляет интерес разобрать вычисление таких функциональных интегралов в общем виде.

Квадратичные лагранжианы

Рассматриваем произвольный квадратичный лагранжиан вида L=m2x˙2+b(t)xx˙c(t)2x2+d(t)x˙e(t)x, который без потери общности интегрированием по частям можно привести к виду:

(6)L=12mx˙212c(t)x2e(t)x.

Классический путь удовлетворяет уравнению движения:

(7)mx¨cl+c(t)xcl+e(t)=0,xcl(ta)=xa,xcl(tb)=xb,

а для амплитуды перехода мы можем записать:

(8)U(xb,tb;xa,ta)=eiScl/y(ta)=0y(tb)=0Dyexp(itatbdt[m2y˙212c(t)y2])=eiScl/F(tb,ta).

Обратите внимание, что F(tb,ta) не зависит от источника e(t).

Метод Гельфанда-Яглома

На лекции было показано, что

(9)F(tb,ta)=m2πi1f(tb,ta),

где функцию f(tb,ta) следует находить путем решения дифференциального уравнения

(10)2f(t,ta)t2=c(t)mf(t,ta),f(ta,ta)=0,f(t,ta)t|t=ta=1.

Упражнение 1. Найдите f(tb,ta) для свободной частицы и гармонического осциллятора (c=mω2).

Решение. Для свободной частицы находим, что f(t,ta)=At+B, где константы определяются начальными условиями. Окончательно получаем f(t,ta)=tta.

Для осциллятора решение дифференциального уравнения имеет вид: f(t,ta)=Asin(ω(tta)+ϕ). Из начальных условий находим константы и окончательно получаем f(t,ta)=sin(ω(tta))/ω. То есть

(11)f(tb,ta)=sin(ωT)ω,T=tbta.

Связь с классическим действием

Напомним, что классическое действие Scl(xb,tb;xa,ta) позволяет легко вычислять начальный импульс частицы, при ее движении по классической траектории, pa=Scl/xa. Действительно,

(12)Sxa=tatbdt(Lxxxa+Lx˙ddtxxa)=Lx˙xxa|tatb=Lx˙|tapa,

где второе слагаемое было проинтегрировано по частям, после чего мы воспользовались уравнением Эйлера-Лагранжа.

Теперь рассмотрим уравнение движения

(13)ddtLx˙Lx=0

и продифференцируем его по импульсу pa при фиксированном xa, но свободном xb. При этом

(14)Lpa=Lx˙J˙+LxJ,J=xclpa,

поэтому находим

(15)ddt(2Lx˙2J˙)+ddt(2Lx˙xJ)2Lxx˙J˙2Lx2J=ddt(2Lx˙2J˙)+[ddt(2Lx˙x)2Lx2]J=0.

Так как классическая траектория вблизи левого края: xclxa+pat/m+, то начальные условия равны:

(16)J(0)=0,J˙(0)=1/m.

Теперь подставляем наш лагранжиан L=12mx˙212c(t)x2e(t)x и получаем, что

(17)mJ¨+c(t)J=0,

то есть уравнение совпадает с уравнением Гельфанда-Яглома и можно отождествить mJ(t)=f(t,ta). Поэтому находим

(18)1f(tb,ta)=1mJ(tb)=1m1xb/pa=1m2Sclxaxb.

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

(19)U(xb,tb;xa,ta)=12πi|2Sclxaxb|eiScl/.

Упражнение 2. Проверить формулу для свободной частицы и гармонического осциллятора.