Основы интегрирования по путям в квантовой механике, часть 2.

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

  1. Р. Фейнман и А. Хибс, "Квантовая механика и интегралы по траекториям", глава 3.5
  2. R. Rattazzi, "The Path Integral approach to Quantum Mechanics", глава 1.4.

Классическая траектория и квантовые флуктуации

На прошлом занятии мы установили, что амплитуда перехода частицы из точки xa в другую точку xb за время T может быть представлена в виде функционального интеграла

(1)U(xa,xb;T)=Dx(t)eiS[x(t)]/.

Рассмотрим альтернативный способ вычислить такой объект на примере свободной частицы. Для этого представим все возможные траектории x(t) в виде классического пути и отклонений от него, то есть

(2)x(t)=xcl(t)+δx(t),xcl(t)=xa+xbxatbta(tta).

Поскольку начальная и конечная точки траектории фиксированы, то δx(ta)=δx(tb)=0. Отклонения δx(t) связаны с квантовыми флуктуациями.

Найдем действие для свободной частицы, воспользовавшись введенной декомпозицией:

(3)S=m2tatbdt{x˙cl2(t)+2x˙cl(t)δx˙(t)+[δx˙(t)]2}=Scl+mx˙clδx|tatbmtatbdtx¨clδx+m2tatbdt[δx˙(t)]2=Scl+m2tatbdt[δx˙(t)]2.

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

(4)S=Scl+12tatbdttatbdtδ2Sδx(t)δx(t)δx(t)δx(t)|x(t)=xcl(t)+

Тогда для амплитуды перехода частицы справедливо:

(5)U(xa,xb;T)=Dx(t)eiS[x(t)]/=eiScl/F0(ta,tb).

Для свободной частицы находим

(6)F0(tbta)=Dδx(t)exp[itatbdtm2(δx˙)2].

Отметим, что поскольку δx(t) равна нулю на концах, то фактор F0 не зависит от xa и xb, а определяется исключительно параметрами ta и tb. Трансляционная инвариантность во времени сводит эту зависимость к F0(ta,tb)=F0(tbta). Выражение для классического действия имеет вид:

(7)Scl=m2(xbxa)2tbta.

Фактор F0(tbta) вычислите самостоятельно в качестве упражнения и сравните с результатом, полученным на прошлом семинаре. Правильный ответ имеет вид:

(8)F0(tbta)=m2πi(tbta).

Представление Гейзенберга

В представлении Шредингера операторы физических величин не зависят от времени, зато эта зависимость содержится в волновых функциях, что в конечном счете и приводит к тому, что средние значения, измеряемые экспериментально, зависят от времени. Математический аппарат квантовой механики можно сформулировать альтернативным образом, перенеся всю зависимость от времени на операторы. Такой подход известен как представление Гейзенберга. Если физическая величина описывается оператором O^ в представлении Шредингера, то в представлении Гейзенберга ей соответствует оператор

(9)O^(t)=eiHtO^eiHt.

В частности, для оператора координаты x^ получаем

(10)x^(t)=eiHtx^eiHt.

В представлении Гейзенберга также можно ввести базисные вектора. Пусть задан момент времени t, рассмотрим вектор

(11)|x,teiHt|x,

который является собственным для x^(t) (обратите внимание на знак "+" в экспоненте, который противоположен знаку "-", соответствующему эволюции в представлении Шредингера), поскольку

(12)x^(t)|x,t=x|x,t.

В этих терминах, амплитуда перехода из состояния xi в момент времени ti в состояние xf в момент времени tf может быть записана как

(13)xf,tf|xi,ti=Dx(t)eiS[x(t)]/.

Усредняем с помощью интеграла по траекториям

Теперь рассмотрим произвольную функцию A(x). Она определяет оператор в Гильбертовом пространстве

(14)A^(x^)|xx|A(x)dx.

Рассмотрим функциональный интеграл

(15)D[x]A(x(t1))eiS[x]/=D[xa]D[xb]dx1A(x1)eiS[xa]+S[xb]=dx1xf|eiH(tft1)/|x1x1|eiH(t1ti)/|xiA(x1)=xf|eiH(tft1)/A^(x^)eiH(t1ti)/|xi=xf|eiHtf/A^(x^(t1))eiHti/|xi=xf,tf|A^(x^(t1))|xi,ti.

Здесь мы разбили всю траекторию движения на два участка a и b, которые соединены точкой x1 в момент времени t1. Аналогичным образом можно установить смысл и более общего функционального интеграла, под которым стоит разновременное произведение двух функций. Повторяя аналогичные рассуждения, получаем:

(16)D[x]O1(x(t1))O2(x(t2))eiS[x]/=xf,tf|T[O^1(t1)O^2(t2)]|xi,ti,

где хронологическое упорядочение имеет смысл:

(17)T[O^1(t1)O^2(t2)]=θ(t2t1)O^2(t2)O^1(t1)+θ(t1t2)O^1(t1)O^2(t2).

Упражнение 1. С помощью функционального интеграла вычислите xf,tf|x^(t)|xi,ti для свободной частицы.

Решение. Нам необходимо вычислить функциональный интеграл

(18)Dxx(t)eiS[x]/.

Представим все возможные пути в виде x(t)=xcl(t)+δx(t), где соответствующие величины были определены в первом разделе. Для действия воспользуемся найденным представлением

(19)S=Scl+m2tatbdt[δx˙(t)]2.

В итоге получаем

(20)Dxx(t)eiS[x]/=xcl(t)eiScl/F0(tbta)=xcl(t)U(xa,xb;T).

Источники и дифференцирование по ним

Экспериментально, чтобы определить свойства физической системы (например, атома), мы позволяем ей взаимодействовать с внешним возмущением (источник, например, электромагнитная волна) и изучаем отклик на это воздействие. Подобная схема имеет формальный аналог и в теоретическом описании.

Рассмотрим систему с лагранжианом L(x,x˙) и позволим ей взаимодействовать с произвольным внешним источником, то есть добавим слагаемое

(21)L(x,x˙)L(x,x˙)+f(t)x(t).

Теперь амплитуда перехода из состояния xa в состояние xb становится функционалом от f(t):

(22)U(xa,xb,T;f)=Dxexp{idt[L(x,x˙)+f(t)x(t)]}.

Этот функционал содержит в себе всю физическую информацию, которая нужна для нахождения средних значений x^ (и любых операторов, которые являются функциями x). В частности,

(23)δiδf(t1)δiδf(tn)U(xa,xb,T;f)|f=0=xb,tb|T[x^(t1)x^(tn)]|xa,ta.

Свободная частица с источником

В качестве примера рассмотрим движение свободной частицы. Добавляем слагаемое с источником в действие и получаем:

(24)U(xa,xb,T;f)=Dxexp{idt[mx˙2(t)2+f(t)x(t)]}.

Траектории представим в виде x(t)=xcl(t)+δx(t). Эта подстановка позволяет представить действие S в виде суммы S=Scl+δS, причем флуктуационная часть не будет зависеть от источника f(t), так как она квадратична по отклонениям δx. Поэтому выражение для δS остается прежним и определяется движением для свободной частицы. Формально это следует из равенства:

(25)tatbdt(mx˙clδx˙+fδx)=mx˙clδx|tatbtatbdt(mx¨clf)δx=0.

Таким образом, получаем

(26)U(xa,xb,T;f)=m2πiTeiScl/.

Теперь вычислим классическое действие Scl. Для этого нужно найти классическую траекторию, которая удовлетворяет уравнению движения:

(27)mx¨cl(t)=f(t),x(0)=xa,x(T)=xb.

Ответ удобно представить в виде xcl=xfree+xf, где первое слагаемое соответствует движению в случае равенства нулю внешней силы

(28)xfree=xa+xbxaTt,

а второе удовлетворяет уравнению

(29)mx¨f=f,xf(0)=0,xf(T)=0.

Решение этой граничной задачи запишем через функцию Грина

(30)xf(t)=0TdtG(t,t)f(t).

Упражнение 2. Найдите функцию Грина G(t,t).

Решение. Искомая функция Грина удовлетворяет уравнению d2dt2G(t,t)=1mδ(tt) и граничным условиям G(0,t)=G(T,t)=0. Найдем левое u(t) и правое v(t) решения.

Для левого решения имеем t2u=0 и u(0)=0, откуда следует, что u(t)=t. Аналогично, для правого решения найдем v(t)=tT. Вронскиан решений равен W(t)=uvuv=T. Окончательно находим:

(31)G(t,t)=1mT{t(tT),если t<t(tT)t,если t>t

После нахождения классической траектории мы можем приступить к нахождению классического действия. Получаем

(32)Scl=0Tdt(mx˙cl22+fxcl)=Sfree+0Tdtfxfree+Sf.

Здесь Sfree определяется выражением (7), а для последнего слагаемого находим:

(33)Sf=0Tdt[fxf+m2x˙f2]=120Tdtfxf=120Tdt0Tdtf(t)G(t,t)f(t).

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

(34)U(xa,xb,T;f)=m2πiTexp[im2T(xbxa)2]×exp[i0Tdtxfree(t)f(t)]exp[i20Tdt0Tdtf(t)G(t,t)f(t)].

Упражнение 3. Воспроизведите результат упражнения 1 с помощью соотношения (23).

Упражнение 4. Вычислите парную корреляционную функцию xb,tb|T[x^(t)x^(t)]|xa,ta.

Решение. Необходимо вычислить вторую производную

(35)δiδf(t)δiδf(t)U(xa,xb,T;f)|f=0=U(xa,xb,T)[xfree(t)xfree(t)iG(t,t)].

Флуктуации квадрата скорости и типичная траектория

Упражнение 5. С помощью функционального интеграла вычислите величину

(36)xb,tb|T[(x^(t+τ)x^(t))2τ2]|xa,ta,

где предполагается, что τT. Покажите, что ответ будет иметь вид IU(xa,xb,T)[x˙free2+imτ].

Первое слагаемое соответствует квадрату скорости при движении по классической траектории, а второе описывает квантовые флуктуации. Если время τ</(mx˙free2), то второй член доминирует.

Полученный результат применим, если в качестве τ взять время Δt – элементарный шаг в дискретной версии функционального интеграла. Тогда мы приходим к оценке

(37)ΔxΔtmΔtΔxΔtm.

Это означает, что в пределе Δt0 типичное отклонение Δx намного больше ожидаемого, то есть типичная траектория "скачкообразна".