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

электронное научно-техническое издание

НАУКА и ОБРАЗОВАНИЕ

Эл № ФС77 - 30569. Государственная регистрация №0421100025. ISSN 1994-0408

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

77-30569/239563

# 10, октябрь 2011 Я. А. Бутко

УДК 517.987.4

МГТУ им. Н.Э. Баумана yanabutko@yandex.ru

1. Введение

В настоящей работе рассматривается новый метод исследования и описания линейной динамики. Метод основан на представлении соответствующих эволюционных полугрупп (или, что то же самое, решений соответствующих эволюционных уравнений) с помощью формул Фейнмана, то есть в виде пределов конечнократных интегралов при стремлении кратности к бесконечности. Как известно, для многих начально-краевых задач функции Грина неизвестны в явном виде. В то же время для некоторых таких задач удается получить формулы Фейнмана, содержащие конечнократные интегралы только от элементарных функций. Такие формулы Фейнмана позволяют проводить непосредственные вычисления решений эволюционных уравнений, пригодны для аппроксимации переходных вероятностей случайных процессов, полезны для компьютерного моделирования стохастической и квантовой динамики. Термин "формула Фейнмана" был введен в статье [22]; в серии работ [22]-[27] был развит метод получения формул Фейнмана для эволюционных уравнений на основе использования теоремы Чернова [10]. В последнее десятилетие этот метод активно применяется для описания различных типов динамики в областях евклидовых пространств и римановых многообразий, в бесконечномерных линейных и нелинейных пространствах, при исследовании -адических аналогов уравнений математической физики (см., например, [1] — [9], [17] — [21]); а также

для построения поверхностных мер на бесконечномерных многообразиях (см., например, [23] — [27]).

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

2. Предварительные сведения

В настоящей работе рассматриваются эволюционные уравнения следующего вида: f (t, q) = Af (t, q), где A некоторый линейный оператор, действующий на функцию f (t, •) переменной q £ Q, Q — некоторое пространство, которое мы будем называть конфигурационным пространством системы, описываемой этим уравнением, t > 0. Рассмотрим задачу Коши

f (t, q ) = Af (t, q), f (0, q) = fo(q).

Если A — это ограниченный оператор на некотором банаховом пространстве (X, У • ||X) функций переменной q и f0 £ X, то решение задачи Коши предста-вимо в виде f (t, q) = (etAf0)(q), где оператор etA определяется как сумма ряда

то n

etA = ^ An, причем ряд сходится в равномерной операторной топологии.

n=0 П'

При этом, ||etA|| < et|A|, то есть оператор etA снова является ограниченным для любого t > 0. Кроме того, как видно из определения etA, справедливы соотношения etA о esA = e(t+s)A и e0A = Id, где Id — тождественный оператор на X.

Как правило, однако, оператор A не является ограниченным. В этом случае приведенная выше схема решения задачи Коши (1) обобщается описанным ниже образом. Пусть символ L(X) обозначает пространство всех непрерывных линейных операторов на X с сильной операторной топологией. Если Dom(A) С X — это линейное подпространство и A : Dom (A) —> X — линейный оператор, то Dom(A) означает область определения A. Однопараме-трическое семейство (Tt)t>0 ограниченных линейных операторов Tt : X — X

называется сильно непрерывной полугруппой, если Т0 = Ы, Т^ = Т ◦ Т для всех я, £ > 0 и Нш^0 \\Т^ф — ф\\х = 0 для всех ф £ X. Если (Т^>0 — сильно непрерывная полугруппа на банаховом пространстве (X, \ \ • \ \ х), то генератором этой полугруппы называется оператор А, определенный по формуле

Т4ф — ф

Аф := lim

t-ю

с областью определения

Г | Тф — ф "1

Dom (А) := < ф G X 1 lim-£- существует как сильный предел >.

Таким образом, если А — ограниченный оператор на X, то Tt = etA — это сильно непрерывная полугруппа на X с генератором А. И в случае, если генератор А — неограниченный оператор, будем иногда использовать обозначение

tA

e для соответствующей полугруппы.

Для корректно поставленной в банаховом пространстве X задачи Коши (1) ее решение представляется в виде f (t, q) = Ttf0(q) [12]. И значит решение задачи Коши (1) равносильно построению полугруппы (Tt)t>0 с заданным генератором А. Как правило, полугруппу (Tt)t>0 не удается получить в явном виде, но удается различными методами ее аппроксимировать. В настоящей работе используется метод приближения, основанный на теореме Чернова (см. [22]).

Теорема 1 (Чернова). Пусть X банахово пространство, F : [0, то) ^ L(X) — (сильно) непрерывное отображение такое, что F(0) = Id и ||F(t)|| < ect для некоторой константы c G [0, то) и всех t > 0. Пусть D — это линейное подпространство D(F'(0)) такое, что сужение оператора F'(0) на D замыкаемо. Пусть (А^(А)) — соответствующее замыкание. Если (А^(А)) является генератором сильно непрерывной полугруппы (Tt)t>0, то для всех t0 > 0 последовательность операторов (F(t/n))n)nGN сходится к (Tt)t>0 при n ^ то в сильной операторной топологии равномерно по t G [0, t0], то есть Tt = limn^TO(F (t/n))n.

Заметим, что производная в нуле функции F : [0, e) ^ L(X), е > 0, — это линейное отображение F'(0) : D(F'(0)) ^ X такое, что

F'(0)g:=lim Ä—,

v >J t\0 t

t

где D(F'(0)) —векторное пространство всех тех элементов g £ X, для которых этот предел существует.

Семейство операторов (F(t))t>0 называется эквивалентным по Чернову полугруппе (Tt)t>0 (будем обозначать это так: F (t) ~Tt), если это семейство удовлетворяет всем требованиям теоремы Чернова по отношению к этой полугруппе, то есть по теореме Чернова в пространстве L(X) локально равномерно по t выполняется равенство

Tt = lim (F(t/n))n. (2)

Это равенство мы будем называть формулой Фейнмана. Мы используем такую терминологию, так как во многих случаях операторы F(t) оказываются интегральными операторами, то есть в правой части формулы Фейнмана стоит предел кратных интегралов при возрастании кратности к бесконечности, а именно Ричард Фейнман ([13], [14]) впервые рассмотрел конструкцию функционального интеграла как предела обыкновенных многократных интегралов по пространствам неограниченно возрастающей размерности. Любое представление решения начальной (или начально-краевой) задачи для эволюционного уравнения (или, эквивалентно, представление полугруппы, разрешающей данную задачу) в виде предела кратных интегралов при возрастании кратности к бесконечности мы будем называть формулой Фейнмана.

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

Замечание 1. Пусть операторы A, B, A + B являются генераторами сильно непрерывных полугрупп etA, etB и et(A+B) на некотором банаховом

пространстве X соответственно, причем операторы А и B не коммутируют.

Тогда etA о etB = et(A+B) = etB 0 etA. Тем не менее, можно показать, что

etA ◦ etB — e*(A+B), etB ◦ etA - e*(A+B),

и значит, из теоремы Чернова следует равенство

et(A+B) = lim reПA о enBГ = lim [enB о enAf.

n—^^o n—то

Последняя формула широко известна как формула Троттера.

Замечание 2. Пусть A — ограниченный оператор. Тогда, можно показать, что Id +tA — etA, и значит, по теореме Чернова

etA = lim fId+ -ä] .

n—^\ n J

что обобщает классическую формулу анализа ex = limn—TO (1 + П)П.

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

Пусть Q — некоторое локально компактное отделимое пространство, (X, | • Ух) — банахово пространство некоторых непрерывных функций на Q с нормой

||/У = sup |/ (q)|. Пусть (Tt)t>0 — сильно непрерывная полугруппа на X с

geQ

генератором A. Пусть функция а( ) : Q — [ci, c2] С (0, то) непрерывна на Q, aX С X. Тогда оператор A, определенный по формуле

Аф(д) = a(q)( A0)(q), Dow (А) = Dow (A),

также является (см. [11], [15], [16] и ссылки в них) генератором сильно непрерывной полугруппы, которую мы обозначим символом (Tt)t>0. Оператор А будем называть мультипликативным возмущением генератора A, а полугруппу (Tt)t>o — полугруппой с мультипликативно возмущенным функцией a( ) генератором.

Пусть операторы Tt при всех t > 0 удовлетворяют условию ||Tt|| < etk при некотором k > 0. Рассмотрим семейство операторов (F(t))t>0 на X, действующих следующим образом:

F(t)0(q) = (Ta(g)t0)(q), Vф e X, Vq e Q.

Заметим, что семейство (F(t))t>o уже не является однопараметрической полугруппой.

Теорема 2. Семейство (F(t))t>0 эквивалентно по Чернову полугруппе (Tt)t>0, то есть локально равномерно по t > 0 в L(X) справедлива формула Фейнмана

T = lim |F(t/n)Г.

n—

Проверим выполнение условий теоремы Чернова. Очевидно, что F(0) = T0 = Id. Кроме того, ||F(t)|| < ec2kt, так как

||F(t)0|x = sup |(Ta(g)t0)(q)|

qeQ

< sup |(Ta(qa)i^)(q)|

< SUp ||(Ta(qo)t|| • ||ф||х

qoGQ

< ec2kt||ф||х.

Покажем, что семейство (F(t))t>0 сильно непрерывно:

lim ||F(t)0 - ф||х = lm0 SUP |(Ta(q)t0)(q) - )|

0

qGQ

< lim SUP |(Ta(qo)i0)(q) - 0(q)| 1—0 q,qo€Q

= lim sup ||ТаЫ1ф - ф||х

qoGQ

< lim sup ЦТСф — ф||х 1—0 ce[ci,c2]

= 0.

Далее, для любого ф G Dom(A) имеем:

lim l0

F (1)Ф—Ф

— Аф

= lim sup

X 1—0 qGQ

< lim sup

1—0 q,qoGQ

< lim sup

1—0 q,qoGQ

C2t

— а(д)Аф(д)

(Ta(qo)t ф)Ы— ф(д)

1

a(qo)i

— а(^0)Аф(д)

< lim 1f ||А(Тяф — ф)||хds

10

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

=0

1

a

Здесь мы воспользовались стандарными свойствами сильно непрерывных полугрупп: Т5ф Е Эош(А) для любого ф Е Эош(А) и любого й > 0, а также

t

Т^ф — ф = А / Т^ф^й для любого ф Е X. о

Таким образом, все условия теоремы Чернова выполнены, а значит, семейство (£)^>0 эквивалентно по Чернову полугруппе (Т )^>0. ►

Следствие 1. (ср. с [9]) Пусть Q = К^, (Т^>0 — феллеровская полугруппа, т.е. сильно непрерывная сжимающая сохраняющая положительность полугруппа на пространстве (М^) непрерывных убывающих на бесконечности к нулю функций. Пусть РДд, ¿у) — переходная вероятность соответствующего феллеровского случайного процесса, то есть Т^ф(д) = / ф(у)Рг(#, ¿у).

Тогда, по теореме 2, семейство операторов (^(£)^>0, действующее по формуле

^ЖзО = ^ ф(у),

м^

эквивалентно по Чернову полугруппе (Т )г>0 с мультипликативно возмущенным функцией а( ) генератором. И значит, справедлива формула Фейнмана

Т;фЫ = Иш / ••• / ф(дп •••

п^^ м^ Ум" (3)

• • • Ра(д„_1)£/п

(Яп—1,^Яп)

Замечание 3. Мультипликативное возмущение генератора марковского процесса равносильно процедуре случайной замены времени (см. [5]). Отметим, что Р(£, х, ¿у) = Ра(ж)г(ж, ¿у) в общем случае НЕ является переходной вероятностью какого-либо случайного процесса. Тем не менее, если переходная вероятность Р^(#, ¿у) исходного процесса известна, то формула (3) позволяет аппроксимировать не выражающуюся в явном виде переходную вероятность модифицированного процесса.

Пример 1 (формула Фейнмана для диффузии с переменным коэффициентом). Рассмотрим процесс броуновского движения в М^. Генератор этого процесса — это оператор Лапласа А = 1 А. Плотность переходной вероятности задается гауссовской экспонентой

р (х) = (2п^/2 ехр( —

2

Пусть (Тг)г>о — полугруппа на С(М^) с мультипликативно возмущенным функцией а() генератором А Эта полугруппа соответствует диффузионному процессу с переменным коэффициентом диффузии л/а. Тогда по Теореме 2 для любого ф Е справедлива формула Фейнмана:

Т(фЫ = Нт / • • • /РМ«,)^2 ехр { _ ^^ } • • •

- ТГ^-^^ ? ф(Тп)2Т1 • • • 2ТП

Пример 2 (Процесс типа Коши с переменным коэффициентом). Рассмотрим процесс Коши в М^. Плотность его переходной вероятности задается формулой

= Г (2 + 2) [п|х|2 '

где Г( ) — это гамма-функция Эйлера. Рассмотрим мультипликативное возмущение генератора процесса Коши функцией а( ). Тогда для соответствующей полугруппы (7^>0 с мультипликативно возмущенным генератором по Теореме 2 имеем:

Т*ф(то) = Иш / • ••

г( 2+2

V 2 2

а(то)£/п

[(а(то)^/п)2 + (п|то _ Т1|)2Р+1)/2

а(Тп_1)^/п Ф(Т )2Т 2т

ф(Тп)2Т1 • • • 2Тп.

[(а(тп_1)^/п)2 + (п|Тп_1 _ тп|)2р+1)/2

п

4. Заключение

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

Работа выполнена при поддержке гранта Президента Российской Федерации МК-943.2010.1. и гранта РФФИ 10-01-00724-а.

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

1. Бутко Я.А. Формулы Фейнмана и функциональные интегралы для диффузии со сносом в области многообразия // Мат. Заметки. 2008. Т. 83, N. 3. С. 333-349.

2. Бутко Я.А., Гротхаус М., Смолянов О.Г. Формула Фейнмана для параболического уравнения второго порядка в области // Доклады РАН. 2008. Т. 421, N. 6. С. 727-732.

3. Бутко Я.А., Смолянов О.Г. Формулы Фейнмана в квантовой и стохастической динамике // Современные проблемы математики и механики. 2011. Т. 6, N. 1. С. 61-75.

4. Бутко Я.А., Смолянов О.Г., Шиллинг Р.Л. Формулы Фейнмана для фелле-ровских полугрупп // Доклады РАН. 2010. T. 434, N. 1. C. 7-11.

5. Портенко Н.И., Скороход А.В., Шуренков В.М. Марковские процессы // Итоги науки и техники. Серия Современные проблемы математики. Фундаментальные направления. Т. 46. Теория вероятностей-4. М.: ВИНИТИ, 1989г. 248с.

6. Butko Ya. A. Function integrals corresponding to a solution of the Cauchy-Dirichlet problem for the heat equation in a domain of a Riemannian manifold // J. of Math. Sci. 2008. Vol. 151, no. 1. Pp. 2629-2638.

7. Butko Ya., Grothaus M., Smolyanov O.G. Lagrangian Feynman Formulae for Second Order Parabolic Equations in Bounded and Unbounded Domains // IDAQP. 2010. Vol. 13, no. 3. Pp. 377-392.

8. Butko Ya. A., Schilling R.L., Smolyanov O.G. Hamiltonian Feynman-Kac and Feynman formulae for dynamics of particles with position-dependent mass // Int. J. Theor. Phys. 2011. Vol. 50. Pp. 2009-2018.

9. Butko Ya. A., Schilling R.L., Smolyanov O.G. Lagrangian and Hamiltonian Feynman formulae for some Feller semigroups and their perturbations // IDAQP. 2011, to appear.

10. Chernoff P. Product formulas, nonlinear semigroups and addition of unbounded operators // Mem. Am. Math. Soc. 1974. Vol. 140.

11. Dorroh J.R. Contraction semi-groups in a function space // Pacific J.Math. 1966. Vol. 19, no. l.Pp. 35-38.

12. Engel K. J., Nagel R. One-parameter semigroups for linear evolution equations. Springer, 1995. 586 Pp.

13. Feynman R. P. Space-time approach to nonrelativistic quantum mechanics // Rev. Mod. Phys. 1948. Vol. 20. Pp. 367-387.

14. Feynman R.P. An Operator Calculus Having Applications in Quantum Electrodynamics //1951. Phys. Rev. Vol. 84. Pp. 108-128.

15. Gustafson K, Lumer G. Multiplicative perturbation of semigroup generators // Pacific J. Math. 1972. Vol. 41, no. 3. Pp. 731-742.

16. Lumer G. Perturbation de generateurs infinitesimaux du type "changement de temps" // 1974. Ann. Inst. Fourier. Vol. 23, N. 4. Pp. 271-279.

17. Obrezkov O.O. The Proof of the Feynman-Kac Formula for Heat Equation on a Compact Riemannian Manifold // IDAQP. 2003. Vol. 6, no. 2. Pp. 311-320.

18. Obrezkov O., Smolyanov O.G., Truman A. The Generalized Chernoff Theorem and Randomized Feynman Formula // Doklady Math. 2005. Vol. 71, no. 1. Pp. 105-110.

19. Sakbaev V. G., Smolyanov O. G. Dynamics of a Quantum Particle with Discontinuous Position-Dependent Mass // Dokl. Math. 2010. Vol. 82, no. 1. Pp. 630-634.

20. Smolyanov O.G. Feynman type formulae for quantum evolution and diffusion on manifolds and graphs // Quant. Bio-Informatics, World Sc. 2010. Vol. 3. Pp. 337-347.

21. Smolyanov O.G., Shamarov N.N. Feynman and Feynman-Kac formulae for evolution equations with Vladimirov operator // Doklady Math. 2008. Vol. 77, no. 3.Pp. 345-349.

22. Smolyanov O.G., Tokarev A.G., Truman A. Hamiltonian Feynman path integrals via the Chernoff formula// J. Math. Phys. 2002. Vol. 43, no. 10. Pp. 5161-5171.

23. Smolyanov O.G., Weizsacker H.v. and Wittich O., Diffusion on compact Riemannian manifolds, and surface measures // Doklady Math. 2000. Vol. 61. Pp. 230-234.

24. Smolyanov O.G., Weizsacker H. v., Wittich O. Brownian Motion on a Manifold as Limit of Stepwise Conditioned Standard Brownian Motions // Stochastic Proceses, Physics and Geometry: New Interplays. II: A Volume in Honor of Sergio Albeverio. Ser. Conference Proceedings. Canadian Math. Society. Providence: AMS. 2000. Vol. 29. Pp. 589-602.

25. Smolyanov O. G., Weizsacker H. v., Wittich O. Chernoff's theorem and the construction of semigroups // Evolution Equations: Applications to Physics, Industry, Life Sciences and Economics. Birkhauser, Prog. Nonlinear Differ. Eq. Appl. 2003. Vol. 55. Pp. 349-358.

26. Smolyanov O.G., Weizsacker H.v., Wittich O. Surface Measures and Initial Boundary Value Problems Generated by Diffusions with Drift // Doklady Math. 2007. Vol. 76, no. 1. Pp. 606-610.

27. Smolyanov O.G., Weizsacker H. v., Wittich O. Chernoff's Theorem and Discrete Time Approximations of Brownian Motion on Manifolds // Potent. Anal. 2007. Vol. 26, no. l.Pp. 1-29.

electronic scientific and t echnical periodical

SCIENCE and EDUCATION

El № FS77 - 30569. №0421100025. ISSN 1994-0408

Feynman formula for semigroups with multiplicatively perturbed generators

77-30569/239563

# 10, October 2011 Ya. A.Butko

Bauman Moscow State Technical University

yanabutko@yandex.ru

A new method to describe linear dynamics is considered. This method is based on representations of corresponding evolution semigroups (or, what is the same, representations of solutions of corresponding evolution equations) by Feynman formulae, i.e. by limits of n-fold iterated integrals when n tends to infinity. Green functions of many initial-boundary value problems are not known explicitly, whereas it is possible to obtain Feynman formulae containing only elementary functions as integrands for some of these problems. Such Feynman formulae allow to calculate solutions of evolution equations directly, to approximate transition probabilities of stochastic processes, are useful for computer modeling of quantum and stochastic dynamics. The notion "Feynman formula" (in this context) and the method to obtain such formulae were introduced in works of Smolyanov and his coauthors in the late nineties. In the last decade it has been actively applied to describe different types of dynamics in domains of Euclidean spaces and Riemannian manifolds, in infinite dimensional linear and non-linear spaces. In the present note the multiplicative perturbations of generators of strongly continuous semigroups on a Banach space of some continuous functions are considered. A Feynman formula is obtained for the semigroups with perturbed generators. Therefore, a new formula is given for the description and the investigation of the perturbed dynamics. Also some particular examples, when the obtained Feynman formula contains only elementary functions as integrands, are considered in this note. Keywords: Feynman formula, approximation of semigroups, approximation of transitional probabilities, multiplicative perturbations

References

1. Butko Ya. A. Feynman Formulas and Functional Integrals for Diffusion with Drift in a Domain on a Manifold // Math. Notes. 2008. Vol. 83, no. 3. Pp. 301-316.

2. Butko Ya.A., Grothaus M., Smolyanov O.G. Feynman Formula for a Class of Second-Order Parabolic Equations in a Bounded Domain // Doklady Math. 2008. Vol. 78, No. l.Pp. 590-595.

3. Butko Ya.A., Smolyanov O.G. Formuly Fejnmana v kvantovoj i stohasticheskoj dinamike // Sovremennye problemy matematiki i mehaniki. 2011. Vol. 6, no. 1. Pp. 61-75.

4. Butko Ya.A., Schilling R.L., Smolyanov O.G., Feynman formulae for Feller semigroups // Doklady Math. 2010. Vol. 82, no. 2. Pp. 679-683.

5. Portenko N.I., Skorohod A.V., Shurenkov V.M. Markovskie processy // Itogi nauki i tehniki. Serija Sovremennye problemy matematiki. Fundamental'nye napravlenija. Vol. 46. Teorija verojatnostej-4. M.: VINITI, 1989. 248 p.