Симплектические томограммы в представлении фейнмановскими интегралами по траекториям Текст научной статьи по специальности «Математика»

УДК 530.145.1

А. К. Федоров, С. О. Юрченко

T(e, ß, п) =

2тг|п| J

м 2

x2 ) ф(xi)dxidx2

Показано, что предлагаемый метод позволяет находить решения для эволюционного уравнения Фоккера-Планка в континуальных интегралах. Найдены квантовые томограммы гармонического осциллятора.

E-mail: akfedorov@student.bmstu.ru

Ключевые слова: квантовая томография, уравнение Шредингера, интеграл Фейнмана.

В настоящее время активно развивается новый подход к квантовой механике, основанный на использовании неотрицательных функций распределения вероятности — симплектических томограмм [1-9]. Они представляют интерес, поскольку могут быть измерены экспериментально в квантово-оптических экспериментах [1, 2]. При этом томографическая формулировка эквивалентна другим известным на сегодняшний день формулировкам квантовой механики, а томограмма связана с различными функциями квазираспределений, например: W-функцией Вигнера [4, 5], P-функцией Глаубера-Сударшана [6, 7], Q-функцией Хусими-Кано [8].

Томографическое представление квантовой механики использует линейное каноническое преобразование фазового пространства как действие симплектической группы SP2(R):

где д,р и е, а — обобщенные координаты и импульсы, п, п —

произвольные действительные постоянные с условием — пп = 1.

Преобразование (1) является каноническим. По определению оно сохраняет невырожденные кососимметрические билинейные формы: скобку Пуассона в классическом случае и коммутатор — в квантовом. В классическом случае отсюда следует, что уравнения Гамильтона выполнены как для новых е, а, так и для старых переменных д и р. Уравнения Гамильтона следуют из принципа наименьшего действия,

(1)

поэтому

¿У (pdq - H (q,p)dt) = 8 J (a de - H(e, a)dt) = 0, (2)

где H(e, a) — гамильтониан системы после преобразования (1).

Частным случаем преобразования (1) является действие матриц вращения группы SO(R)

cos 9 sin 9N ч— sin 9 cos 9,

где 9 е R/2nZ — угол поворота.

В таком случае преобразование (1) является поворотом фазового пространства на угол 9, а в общем случае (1) — это поворот фазового пространства с взаимным масштабированием по осям p и q.

Симплектическая томограмма T(e, j, п) наблюдаемой e, которая является линейной комбинацией квадратурных компонент

e = jq + ПР,

определяется через волновую функцию

T (e,j,n) = \FtJ,;n [ф}(e)\2,

где — линейный унитарный оператор. Таким линейным унитарным оператором является интегральный оператор дробного преобразования Фурье

T {e,p,n)= 1

2n|n|

ф{х)ехр^ i^x2--Xjdx

(3)

,2n n

где x — обобщенная координата.

Исходя из определения (3) томограмма T(е, j, n) представляет собой положительную нормированную и однородную (с порядком — 1) функцию:

T(e,j,n) > 0, Jt(e,j,n)de = 1, T(e,j,n) = |Л|T(\e,\j,Xq);

Л — произвольная постоянная.

Свойство однородности позволяет записать для томограммы уравнение Эйлера дГ дГ дГ

+ 8JJj+v=—T ■ (4)

Решение уравнения (4) имеет вид

Tf (e,J,n) = Л f(-,~

\J\ \J J

где f — произвольная функция.

2

Решение уравнения Эйлера для томограммы позволяет снизить ее размерность. Несмотря на то, что томограмма является функцией трех переменных, благодаря нормировке ее представление возможно через функцию двух переменных. Подобные функции называются томограммами Френеля.

Целью данной работы является представление симплектической томограммы через интеграл Фейнмана по траекториям. Предложенный ранее метод классического пропагатора [9] требует знания начального распределения вероятности. Это является существенным усложнением, поскольку такой формализм предполагает измерения как начального, так и конечного состояния системы. Технически, с учетом особенностей экспериментальной реализации [1, 2], это означает ансамблевые измерения начального и конечного состояний. В случае использования некоторой модельной волновой функции и предлагаемого метода становится достаточно однократного измерения.

Пусть эволюция рассматриваемой квантово-механической системы происходит согласно уравнению Шредингера. Рассмотрим задачу Коши для уравнения Шредингера

1 д

—2Аф(х,г) + у = г—ф(х,^);

ф(х,Ь) | 0 = ф0(х), т = К = 1;

волновые функции ф(х, ¿) будем рассматривать с учетом стандартного условия нормировки

/|ф(х,<)|2 Лх = 1-

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

Для произвольных моментов времени ¿1 и ¿2 справедливо соотношение

ф(х2^2) = ! С(х2,Х1,^2 — ¿1)ф(х1, ¿1)^X1, (5)

где С(х2, х1, ¿2 — ¿1) — функция Грина уравнения Шредингера.

Функция Грина С(х2,х1,^2 — ¿1) представляет собой амплитуду перехода из начального состояния с волновой функцией ф(х1, ¿1) в состояние с волновой функцией ф(х2,£2). Таким образом, задача Коши для дифференциального уравнения Шредингера свелась к интегральному уравнению (5), которое решается численно, например, методом

Монте-Карло [10, 11]. При этом функция Грина обладает свойством

G(X2,X1,Í2 - ti) =0 Vil > Í2.

В работе Фейнмана [12] показано, что при переходе (x1, í1)^(x2, t2) функция Грина G(x2, x1, t2 — t1) может рассматриваться как интеграл по траекториям

G(x2,x1,t2 — t1) = J exp(iS(x2,x1 ))D[x(t)], (6)

n

где S(x2, x1) — функционал действия, инвариантный относительно канонических преобразований (1) [13]; под D[x(t)] поздразумевается интегрирование по конфигурационному пространству П функций x(t).

В настоящей работе рассматриваются линейные канонические преобразования фазового пространства. Поэтому будет использована га-мильтонова формулировка интегралов по траекториям — интегралы по фазовому пространству обобщенных координат x = (x1,..., xn) и импульсов p = (p1,... ,pn). Впервые подобная концепция изложена в работе [14].

Предполагается, что задан гамильтониан системы H

H = p2 + V (x).

В силу справедливости вариационного принципа (2) интеграл

I = J (px — H) dt принимает экстремальное значение. Его можно представлять в виде

I = /(px — pT — V(x))dt = £(pAx¿ — ЙAtj — jj V(x)dt).

j=1 t"

Воспользуемся известными формулами [14]

p? д At, ( Аж, )2 lAX

7 - —At, =--- (p, +---x—~.

j 2 7 2 VJ AtJ 2 At,

,, Пусть

e =, A| )■

Вычислим интеграл Эйлера-Пуассона

exp

< - А, , =( - A,3/2

Таким образом, интеграл по траекториям (6) понимается в следующем смысле:

exp

i (px — H)dt d p\.. .d pNd x1.. .d xN_1 =

Ш—j

exp

-i Л Ax2 гj

j=i 3 N

£

At

3 = i

V(x)dt

j Jtj

j-i

d xi . . . d xN i =

—rl 2ni\ 3/2 f

— exp(iS (x))d3x1 ...d3xN _1.

3=i

В соответствие с этим

N

V[x(t)] = lim П

n At,-/

3=i 3

2п% \3/2dzx

Для решения поставленной задачи будем использовать связь (3) между симплектической томограммой Т(е, л, п) и волновой функцией ф(х2^). Для ф(х2,12) справедливо уравнение (5).

Тогда симплектическую томограмму можно представить в виде

Т =

1

2ж\п\

G(x2, х1, t2 — t1)exp( %—x\ — —х2 ) ф(x1)dx1dx2 \2ц п /

. (7)

Или с учетом представления функции Грина через интеграл Фей-нмана по траекториям (6)

Т =

1

2п\п\

V[x(t)] exp(iS) exp( --х2 )^(x1)dx1dx2

2п п

. (8)

Уравнение (8) представляет собой основной результат настоящей работы. Таким образом, можно ввести функцию

P (e,^,n)=j' G(x?,x1,t2 — t^exp^^xX? — dx?

2

2

Функцию Р можно, как и функцию Грина, представить в виде интеграла Фейнмана

P (e, n) = exp( ¿S + —ж2 — —ж2) D[x(t)]dx2 I V 2n n /

и записать результат (8) как

т (e,^,n) =

2n|n|

P (e, n)^(x)dx

Полученный результат предполагает следующую физическую интерпретацию: если известны начальное условие для уравнения Шредингера и гамильтониан системы, то результат гомодинного детектирования (экспериментальное измерение томограммы) в экспериментах типа [1,2] выражается через уравнение (8).

Применяя неравенство Коши-Буняковского-Шварца можно получить следующую оценку:

1

т(e^n) <

2n|nl

|P(e,^,n)| dx

(9)

Благодаря свойствам томограммы Т(е, п) условие нормировки можно сформулировать в общем виде:

„ 2

1

2n|n|

P (e, n)^(x)dx

de = 1.

(10)

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

Н = ^2 + V (ж), (Я = ш = т = 1)

имеет вид обобщенного уравнения Фоккера-Планка [3]

дт _ dT •

dt ^ дп

гд ¿п д гд ¿п д

т, (11)

где /б — оператор интегрирования по переменной е.

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

дт öt

^дТ = 0.

дп

(12)

2

2

1

Для гармонического осциллятора

V й-1 ■ ^ - ^ + = 0. (13)

Для параметрического осциллятора с частотой ш(г)

£ - + -0. (14)

Как уже было отмечено, томограмма представляет собой решение уравнения Фоккера-Планка. Поэтому формула (8) позволяет искать решение уравнения (11) с учетом связи между томограммой и волновой функцией (3). Задачи о представлении решения для данного уравнения рассматривались ранее в работах [15, 16], потому что существует его тесная связь с формулировкой Фейнмана интегралов по траекториям и представлением решения уравнения Шредингера [17, 18].

Решение уравнения Фоккера-Планка для гармонического осциллятора. В качестве примера представим с помощью найденной формулы (8) решение уравнения Фоккера-Планка для гармонического осциллятора (13). Как известно [12], функция Грина имеет вид

1 /г[(ж? + ж^соэг - 2x1x2]\

е(Х2'Х1'г)- 72ПШ ехр(--). (15)

Запишем решение уравнения Фоккера-Планка для функции Грина (15) как

т (б,^,п)= 1

2п|п|

(ф 2 ie \ f dxi^o

dx2 exp— x2--x2 =x

V2n n ' J V2nisint

x ('i[(x1 + x2) cos t — 2x1x2] ^

(16)

2 sin t

Вычислив P-функцию уравнения (16), можно получить волновые функции осциллятора. В общем случае решение имеет вид

2-n ( е2 \ ( е \ = , ¡ . 2 , = ехЫ —н2 ^-=), (17)

nl^n^2 + п2) v + П2/ v vV + П где Hn(x) — полиномы Эрмита.

Например, томограмма основного состояния осциллятора

1 ( е2 \

То (е,^,п) = , , 2 2N exp ———2 . (18)

+ п2) V + п2^

Томограмма Tí первого возбужденного состояния может быть легко вычислена:

2

2е

х

(р2 + п2) \/ п(р2 + п2)

(19)

Томограмма основного (п = 0) и первого возбужденного состояния (п = 1) гармонического осциллятора

(г = + п2)-1/2)

Г1(е,р,п)

х exp

1 р2 + п

Следует отметить, что в случае соотношения параметров системы отсчета р = 1 и п = 0, томограммы состояния вырождаются в обычные плотности вероятности в координатном представлении, что видно на рисунке. Численное решение уравнения (13) совпадает с найденным при помощи формулы (8) представлением (16) томограммы основного и первого возбужденного состояний.

Заключение. Предложено представление симплектических томограмм через интегралы Фейнмана по траекториям (8). Показан метод решения интегро-дифференциального уравнения Фоккера-Планка (11) в томографическом представлении квантовой механики на основе использования связи между симплектической томограммой и волновой функцией, а также метода интегралов по траекториям. В качестве примера найдены томограммы различных состояний гармонического осциллятора.

Томографические интегралы по траекториям (7, 8) закладывают основу для применения современных методов численного интегрирования в квантовой томографии, позволяющих проводить численные эксперименты для различных типов квантовых систем. Поскольку задача сводится к вычислению многомерных интегралов, может быть применим вариационный метод Монте-Карло [10, 11], который имеет имеет существенное преимущество в производительности по сравнению с другими численными методами интегрирования в задачах квантовой механики и физики конденсированного состояния.

Авторы признательны профессору Ю.Е. Лозовику и доценту Я.А. Бутко за плодотворные дискуссии. Работа выполнена при поддержке РФФИ.

СПИСОК ЛИТЕРАТУРЫ

1. Beck M., S m i t h e y D. T., Raymer M. G. Experimental determination of quantum-phase distributions using optical homodyne tomography // Physical Review A. - 1993. - Vol. 48. - P. 890-893.

2. S m i t h e y D. T., Beck M., Raymer M. G., Faridani A. Measurement of the Wigner distribution and the density matrix of a light mode using optical homodyne tomography: Application to squeezed states and the vacuum // Physical Review Letters. - 1993. - Vol. 70. - P. 1244-1247.

3. I b o r t A., Man'ko V. I., M a r m o G., S i m o n i A., Ve n t r i g l i a F. An introduction to the tomographic picture of quantum mechanics, arXiv:0904.4439v1 [quant-ph], (2009).

4. W i g n e r E. P. On the quantum correction for thermodynamic equilibrium // Physical Review. - 1932. - Vol. 40. - P. 749-759.

5. M o y a l J. E. Quantum mechanics as a statistical theory // Proceedings of the Cambridge Philosophical Society. - 1949. - Vol. 45. - P. 99-124.

6. G l a u b e r R. J. Photon Correlation // Physical Review Letters. - 1963. - Vol. 10.

- P. 84-86.

7. S u d a r s h a n E. C. G. Equivalence of semiclassical and quantum mechanical descriptions of statistical light beams // Physical Review Letters. - 1963. - Vol. 10.

- P. 227-279.

8. H u s i m i E. Some formal properties of the density matrix // Proceedings of the Physico-Mathematical Society of Japan. - 1940. - Vol. 23. - P. 264-314.

9. М а н ь к о О. В. Классические пропагаторы квадратичных квантовых систем // Теоретическая и математическая физика - 1999. - T. 121, № 2. - C. 285-296.

10. F o u l k e s W. M. C., M i t a s L., N e e d s R. J., R a j a g o p a l G. Quantum Monte Carlo simulations of solids // Reviews of Modern Physics. - 2001. - Vol 73, no. 33.

11. P o 11 o c k E. L., C e p e r l e y D. M. Path-integral computation of superfluid densities // Physical Review B. - 1987. - Vol. 36, no. 33.

12. F e y n m a n R. P., H i b b s A. R. Quantum mechanics and path integrals. - N.Y.: McGraw-Hill, 1965.

13. С м о л я н о в О. Г., С м о л я н о в а М. О. Преобразования интеграла Фейн-мана при нелинейных преобразованиях фазового пространства // Теоретическая и математическая физика. - 1994. - № 100. - C. 3-12.

14. G a r r o d C. Hamiltonian path-integral methods // Reviews of Modern Physics. -1966. - Vol. 38, no. 3. - P. 483-494.

15. Морозов В. Г. О континуальных интегралах для обобщенного уравнения Фоккера-Планка // Теоретическая и математическая физика. - 1984. - T. 58, №1.

- C. 79-95.

16. D r o z d o v A. N., T a l k n e r P. Path integrals for Fokker-Planck dynamics with singular diffusion: Accurate factorization for the time evolution operator // Journal of Chemical Physics. - 1998. - Vol. 109, no. 8. - P. 2080-2091.

17. B o 11 c h e r B., B u t k o Y a. A., Schilling R. L., S m o l y a n o v O. G. Feynman formulas and path integrals for some evolution semigroups related to т-quantization // Rus. J. Math. Phys.. - 2011. - Vol. 18, no. 4. - P. 387-399.

18. B u t k o Y a. A., Schilling R. L., S m o l y a n o v O. G. Hamiltonian Feynman-Kac and Feynman formulae for dynamics of particles with position-dependent mass // Int. J. Theor. Phys. - 2011. - Vol. 50, no. 7. - P. 2009-2018.

Статья поступила в редакцию 3.11.2011

Алексей Константинович Федоров родился в 1993 г. Студент МГТУ им. Н.Э. Баумана. Автор пяти научных работ в области нерелятивистской квантовой механики.

A.K. Fedorov (b. 1993) — student of the Bauman Moscow State Technical University. Author of 5 publications in the field of non-relativistic quantum mechanics.

Станислав Олегович Юрченко родился в 1985 г., окончил МГТУ им. Н.Э. Баумана в 2009 г. Канд. физ.-мат. наук, доцент кафедры "Физика" МГТУ им. Н.Э. Баумана. Автор более 35 научных работ в области теоретической физики.

S.O. Yurchenko (b. 1985) graduated from the Bauman Moscow State Technical University in 2009. Ph. D. (Phys.-Math.), assoc. professor of "Physics" department of the Bauman Moscow State Technical University. Author of more than 35 publications in the field of theoretical physics.