научное издание мгту им. н. э. баумана
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Формулы Фейнмана для эволюционных полугрупп
# 03, март 2014
В01: 10.7463/0314.0701581
Бутко Я. А.
УДК 517.987.4
Россия, МГТУ им. Н.Э. Баумана [email protected]
1. Введение
В работе описывается подход к решению начальных и начально-краевых задач для эволюционных уравнений, основанный на представлении решений таких задач в виде пределов кратных интегралов при стремлении кратности к бесконечности. Такие представления называются формулами Фейнмана. Формулы Фейнмана позволяют проводить непосредственные вычисления решений эволюционных уравнений, пригодны для аппроксимации переходных вероятностей случайных процессов, могут быть использованы для компьютерного моделирования классической, квантовой и стохастической динамики. Настоящий подход был впервые предложен в работах [48, 49] и в настоящее время активно применяется для описания различных типов динамики в областях евклидовых пространств и римановых многообразий, в бесконечномерных линейных и нелинейных пространствах, в пространствах над полем р-адических чисел (см., например, работы [20, 50]).
Формулы Фейнмана позволяют также связать решения эволюционных уравнений с функциональными интегралами. Допредельные выражения в формулах Фейнмана дают аппроксимации функциональных интегралов либо по вероятностным мерам, порожденным некоторыми случайными процессами, либо по псевдомерам Фейнмана, сосредоточенным на траекториях в конфигурационном или в фазовом пространстве некоторой квантовой системы. Представления решений эволюционных уравнений с помощью функциональных интегралов по вероятностным мерам позволяют применить для описания и исследования рассматриваемой эволюции методы стохастического анализа, в частности, метод Монте-Карло. Функциональные интегралы по псевдомерам Фейнмана играют одну из центральных ролей в аппарате современной квантовой физики, представления решений эволюционных уравнений с помощью таких интегралов полезны при исследовании квазиклассических асимптотик, при построении рядов теории возмущений и т.д. Существует множество различных
подходов к строгому определению интегралов Фейнмана, причем различные подходы приводят к различным классам интегрируемых функций. В одном из подходов, восходящем к самому Р. Фейнману [29], рассматривается интеграл Фейнмана по (бесконечному) пространству траекторий как предел интегралов по конечномерным пространствам при стремлении размерности пространств к бесконечности. Формализация такого подхода и приводит к формулам Фейнмана. При этом метод формул Фейнмана позволяет не только строго определить некоторые функциональные интегралы, но и исследовать такие вопросы теоретической физики, как различение процедуры квантования на языке интегралов Фейнмана [6, 13, 19]. Кроме того, метод формул Фейнмана позволяет установить связи между интегралами по псевдомерам Фейнмана и интегралами по вероятностным мерам. Эта связь помогает исследовать и вычислять интегралы Фейнмана. Более того, соответствующие таким интегралам формулы Фейнмана также доставляют средство для их вычисления [22, 39, 42].
В настоящей работе обобщаются и развиваются некоторые из результатов статей [5, 21, 23]. В частности, получены новые формулы Фейнмана для эволюционных полугрупп, порожденных мультипликативными возмущениями генераторов некоторых исходных полугрупп; при этом рассматриваются полугруппы на некоторых банаховых пространствах непрерывных функций, определенных на произвольном метрическом пространстве. Подход к построению формул Фейнмана для полугрупп с мультипликативно и/или аддитивно возмущенными генераторами иллюстрируется на примерах задачи Коши для уравнения Шредингера, аппроксимации переходных вероятностей некоторых марковских случайных процессов. Кроме того, описан подход к построению формул Фейнмана для решения начально-краевой задачи Коши — Дирихле для дифференциального уравнения параболического типа; получены формулы Фейнмана для решения задачи Коши и Коши — Дирихле для параболического уравнения второго порядка с неограниченными переменными коэффициентами. Обсуждаются функциональные интегралы по вероятностным мерам и псевдомерам Фейнмана, которые могут быть вычислены с помощью полученных формул Фейнмана.
В настоящей работе рассматриваются эволюционные уравнения вида — (Ь^) = Lf (Ь, q), где L — некоторый линейный оператор, действующий на функцию Д (Ь, ■) переменной q е Q; Q — некоторое пространство, которое мы будем называть конфигурационным пространством системы, описываемой этим уравнением; Ь > 0. Рассмотрим задачу Коши
Если L — ограниченный оператор на некотором банаховом пространстве (X, || ■ ||х) функций переменной q и е X, то решение задачи Коши представимо в виде f (Ь^) = (е4ЬДо)^), где
2. Предварительные сведения, обозначения
(1)
oo j-ra
оператор определяется как сумма ряда = ^ —¡^га, причем ряд сходится в равномер-
n=0
П!
ной операторной топологии. При этом, ||etL|| < то есть оператор etL снова является
ограниченным для любого t > 0. Кроме того, как видно из определения etL, справедливы соотношения etL о esL = e(t+s)L и e0L = Id, где Id — тождественный оператор на X.
Как правило, однако, оператор L не является ограниченным. В этом случае приведенная выше схема решения задачи Коши (1) обобщается описанным ниже образом. Пусть символ L(X) обозначает пространство всех непрерывных линейных операторов на X с сильной топологией (последовательность операторов Tn е L(X) сходится к оператору T е L(X) в сильной топологии, если для любого x е X выполняется равенство lim ||Tnx — Tx||X = 0). Если Dom(L) С X — линейное подпространство и L: Dom(L) ^ X — линейный оператор, то Dom(L) означает область определения L. Однопараметрическое
семейство (Tt)t>o ограниченных линейных операторов Tt: X ^ X называется сильно непрерывной полугруппой, если T0 = Id, Ts+t = Ts о Tt для всех s, t > 0 и lim ||Tt<£ — <^||X = 0 для всех ^ е X. Если (Tt)t>0 — сильно непрерывная полугруппа на банаховом пространстве (X, || ■ ||X), то
генератором этой полугруппы называется оператор L, определенный по формуле
— ^
Lw := lim-
r t^o t
с областью определения
TtV — <p
Dom(L) := |^ G X: lim-
t
существует как сильныи предел
Таким образом, если L — ограниченный оператор на X, то Tt = e — сильно непрерывная полугруппа на X с генератором L. И в случае, если генератор L — неограниченный оператор, будем иногда использовать обозначение etL для соответствующей полугруппы.
Теорема 1 (ср. [40], стр. 102). Пусть (L, Dom(L)) —плотно определенный на банаховом пространстве X оператор с непустым резольвентным множеством Задача Коши (1) имеет единственное непрерывно дифференцируемое по t е [0, то) решение f для каждого начального условия f0 е Dom(L) тогда и только тогда, когда (L, Dom(L)) является генератором некоторой сильно непрерывной полугруппы (Tt)t>0 на X. При этом решение такой задачи Коши (1) для каждой f0 е Dom(L) представляется в виде f (t, q) = Ttf0(q).
Замечание 1. Резольвентным множеством оператора L называется множество p(L) всех таких комплексных чисел А, для которых оператор R(A, L) = (A Id — L)-1 является ограниченным оператором в X. Семейство {R(A, L)}Aep(L) называется резольвентой оператора L.
Замечание 2. Для всех начальных условий f0 е X функцию f(t,q) = Ttf0(q) по-прежнему будем называть решением задачи Коши (1) (в англоязычной литературе в этом случае используется термин «mild solution», см., например, [40]).
Замечание 3. Если (Tt)t>0 — сильно непрерывная полугруппа на банаховом пространстве X, то существуют константы а > 0, M > 1, такие, что ||Tt|| < Meat для всех t > 0 [40]. Кроме того, если сильно непрерывная полугруппа (Tt)t>0 удовлетворяет неравен-
ству ||Т4|| < Ме"4, то семейство ^)4>0, где St = е также является сильно непрерывной полугруппой на X, причем |^|| < М для всех £ > 0. Введя на X новую норму ||| • |||,
1 1 1 / 1 1 1 = 8ИР4>о ||Stf ^ поЛуЧим, ЧТо Ц < 1 1 1 / 1 1 1 < М ^ ТО есть н°рмЫ ^ • ^ и 1 1 1 • 1 1 1 эКВивалентны. Причем по отношению к норме ||| • ||| полугруппа ^)4>0 является сжимающей полугруппой [27, 28], т.е. |^|| < 1 для всех £ > 0.
Замечание 4. Если (Т4)4>0 — сильно непрерывная полугруппа на банаховом пространстве X с генератором (£, Бош(Ь)), то справедливы также следующие утверждения ([28], предложение 1.5, следствие 1.6):
а) если ^ G X и t > 0, то /0 Tsf ds G Dom(L) и Ttf - f = L / Tfds;
d 0
б) если f G Dom(L) и t > 0, то Ttf G Dom(L) и -Ttf = LTtf = TtLf;
в) если f G Dom(L) и t > 0, то Ttf - f = / LTfds = / TsLf ds;
00
г) генератор (L, Dom(L)) — замкнутый оператор, множество Dom(L) всюду плотно в X.
Итак, решение задачи Коши (1) равносильно построению полугруппы (Tt)t>0 на X с
заданным генератором L. Как правило, полугруппу (Tt)t>o не удается получить в явном виде, но удается различными методами ее аппроксимировать. В настоящей работе используется метод приближения, основанный на теореме Чернова [24, 48]).
Теорема 2 (Чернова). Пусть X — банахово пространство, F: [0, то) ^ L(X) — (сильно) непрерывное отображение, такое, что F(0) = Id и ||F(t) || < eat для некоторой константы а е G [0, то) ивсехt > 0. ПустьD — линейное подпространство Dom(F'(0)), такое, что сужение оператора F'(0) на D замыкаемо. Пусть (L, Dom(L)) — соответствующее замыкание. Если (L, Dom(L)) является генератором сильно непрерывной полугруппы (Tt)t>0, то для всех t0 > 0 последовательность операторов (F(t/n))U)ueN сходится к (Tt)t>0 при n ^ то в сильной топологии равномерно по t G [0, t0], т.е. Tt = lim (F(t/n))U.
Заметим, что производная в нуле функции F: [0,е) ^ L(X), е > 0, — линейное отображение F'(0): Dom(F'(0)) ^ X, такое, что
F'(0)g := lim F- F(0>g,
v t\0 t
где Dom(F'(0)) — векторное пространство всех тех элементов g G X, для которых этот предел существует (в сильной топологии).
Замечание 5. Условие ||F(t)|| < eat, Vt > 0, в теореме 2 можно заменить условием ||Fm(t)|| < Meamt, Vm G N и t > 0 [27].
Замечание 6. Комбинируя условия теоремы 2 с условиями теоремы Хилле — Ио-сиды, можно получить версию теоремы Чернова, утверждающую существование искомой полугруппы [27].
Замечание 7. Существуют различные обобщения теоремы 2. В частности, в книге [28] теорема Чернова обобщена на случай многозначных генераторов (следствие 6.6), в работе [36] теорема Чернова обобщена на случай полугрупп в банаховом пространтсве, непрерывных в более слабом смысле (в смысле локально выпуклой топологии на X более слабой,
чем топология, порожденная нормой), в работе [17] теорема Чернова обобщена на случай полугрупп нелинейных сжимающих отображений на замкнутых выпуклых подмножествах гильбертова пространства.
Семейство операторов (F(t))t>0 называется эквивалентным по Чернову полугруппе (Tt)t>0 (будем обозначать это так: F(t) ~ Tt), если это семейство удовлетворяет всем требованиям теоремы Чернова по отношению к этой полугруппе, то есть по теореме Чернова 2 в пространстве L(X) локально равномерно по t выполняется равенство
Tt = lim (F(t/n))n (2)
Это равенство мы будем называть формулой Фейнмана. Мы используем такую терминологию, так как во многих случаях операторы F(t) оказываются интегральными операторами, то есть в правой части формулы Фейнмана стоит предел кратных интегралов при возрастании кратности к бесконечности. Именно Ричард Фейнман (см. [29]) впервые рассмотрел конструкцию функционального интеграла как предела обыкновенных многократных интегралов по пространствам неограниченно возрастающей размерности. Любое представление решения начальной (или начально-краевой) задачи для эволюционного уравнения (или, эквивалентно, представление полугруппы, разрешающей данную задачу) в виде предела кратных интегралов при возрастании кратности к бесконечности мы будем называть формулой Фейнмана.
Пределы в формулах Фейнмана совпадают с некоторыми функциональными интегралами по вероятностным мерам или по фейнмановским псевдомерам на множестве траекторий некоторой физической системы. Представление решения начальной (или начально-краевой) задачи для эволюционного уравнения (или, эквивалентно, представление поугруппы, разрешающей данную задачу) в виде функционального интеграла обычно называется формулой Фейнмана — Каца. Таким образом, кратные интегралы в формуле Фейнмана для некоторой задачи аппроксимируют функциональный интеграл в формуле Фейнмана — Каца, представляющей решение этой же задачи. Такие аппроксимации во многих случаях представляют собой кратные интегралы только от элементарных функций и, следовательно, могут быть использованы для непосредственных вычислений и моделирования рассматриваемой динамики.
Пример 1. Пусть операторы Li, L2, Li + L2 являются генераторами сильно непрерывных полугрупп etLl, etL2 и et(Ll+L2) на некотором банаховом пространстве X соответственно, причем операторы L1 и L2 не коммутируют. Тогда etLl о etL2 = et(Ll +L2) = etL2 о etLl. Тем не менее, можно показать, что
etLl о etL2 ^ et(Ll+L2), etL2 о etLl ^ et(Ll+L2),
и, значит, из теоремы Чернова следует равенство
et(Ll+L2) = lim ГеnL о enL2ln = lim TenL2 о enli]n. n^^ n—
Последняя формула широко известна как формула Троттера.
Пример 2. Следствием теоремы Чернова является также формула Поста — Виддера [27, 40], обобщающая известную теорему анализа о втором замечательном пределе:
здесь Я(Х, Ь) — резольвента оператора Ь.
Замечание 8. Конечно же, для одной и той же полугруппы еьь существует множество эквивалентны« ей по Чернову семейств операторов. Например, если Т(£) ~ е4Ь, то
с одной стороны, в множестве семейств, эквивалентный по Чернову искомой полугруппе, естественно выбирать оптимальные для проведения непосредственный вычислений. С другой стороны, различные семейства позволяют установить связь между функциональными интегралами по разным мерам и псевдомерам фейнмановского типа.
Равенство (2) будем называть лагранжевой формулой Фейнмана, если Т(£), £ > 0, — интегральные операторы, ядра которых представляются элементарными функциями; если Т(£) — псевдо-дифференциальные операторы, то говорим о гамилътоновой формуле Фейнмана (см., например, работы [21, 22, 23]). Такая терминология связана с тем, что ла-гранжевы формулы Фейнмана дают аппроксимации для функциональны« интегралов по множеству траекторий в конфигурационном пространстве системы (чья эволюция описывается полугруппой (Т)г>°), в то время как гамильтоновы формулы Фейнмана соответствуют интегралам Фейнмана по траекториям в фазовом пространстве некоторой эволюционной системы. Конечно же, существуют и другие типы формул Фейнмана, в частности, формулы Фейнмана, соответствующие функциональным интегралам по траекториям в пространстве импульсов.
В дальнейшем используются следующие обозначения: х • у — скалярное произведение для х, у € Е^; |х| — норма вектора х € Е^; С = С и дС — замыкание области С с Е^;
Ст(Е^) — пространство т раз непрерывно дифференцируемых функций;
С°'Л(Е^) — пространство функций, удовлетворяющих условию Гельдера с параметром
Ст(Е^) — пространство функций / € Ст(Е^), удовлетворяющих условию да/ € СЬ(М^), |а| < т;
Ст'Л(Е^) —пространство функций / € Ст(Е^), удовлетворяющих условию да / €
€ С°'Л(Е^), |а| = т};
Сьт'Л(Ей) = Сьт(^) П Ст'Л(Е^);
Сс(Мй) — пространство непрерытнык на Е функций с компактным носителем;
Ст(^) = Сс(Ей) П Ст(Е^); Сст'Л(Е^) = Сс(Ей) П Ст'Л(Е^);
для любого ограниченного оператора В имеем также [Т(£) + ¿2В] ~ е4Ь. Таким образом,
cc
oo /in>d
) — пространство бесконечно дифференцируемых на Rd функций с компактным
носителем;
Со.
pd
) — пространство всех непрерывных на Rd функций, убывающих на бесконечно-
сти к нулю (это банахово пространство с нормой
sup |f(x)|).
3. Формулы Фейнмана для полугрупп, порожденных аддитивными возмущениями генераторов
Теорема 3. Пусть X — банахово просранство с нормой || ■ ||х. Пусть (Тк(¿))*>о, к = 1, ...,т, — сильно непрерывные полугруппы на X с генераторами Ьк соответственно. Предположим, что оператор Ь = Ь + ■ ■ ■ + Ьт с областью определения Бош(Ь) =
т
= П Бош(Ьк) замыкаем и это замыкание является генератором сильно непрерывной по-
к=1
лугруппы (Т(¿))г>0 на X. Пусть (^(¿))г>0, к = 1, ..., т, — семейства операторов на X, эквивалентные по Чернову полугруппам (Тк(¿))г>0 соответственно, т.е. для каждого к € € {1, ..., т} выполнены соотношения (0) = Ы, (¿)|| < 1 для некоторого ак > 0 и существует множество С Бош(Ьк), являющееся существенной областью определения
оператора Lk, такое, что lim
t-> о
Fk(t)y - ^
t
- Lk P
X
0 при всех p G Dk. Если существует
множество D С Di П ... П Dm, являющееся существенной областью определения для оператора L, то семейство (F(t))t>0, где F(t) = F1 (t) о • • • о Fm(t), эквивалентно по Чернову полугруппе (T(t))t>0 и, следовательно, формула Фейнмана
T = lim [F(t/n)f,
где предел понимается как предел в сильной операторной топологии, справедлива локально равномерно относительно t > 0.
Доказательство. Очевидно, что семейство (F(t))t>0 сильно непрерывно, F(0) = = Id и
||F(t)|| < ||F1(t>n • ... •|Fm(t)|< e(a1
Пусть множество D С D1 П ... П Dm является существенной областью определения для оператора L. Тогда для любого P G D имеем
lim
F (t)P - P
— LP
X
lim
Fi(t) ◦•••◦ Fm(t)^ - P
- Lip - • • • - Lmp
X
lim
t-> 0
р/.ч p Fm(t)P - P r , F1(t) ◦•••◦ Fm-1(t)P - P Fi(t) О • • • о Fm-1(t) о ----LmP +--
- L1P - • • • - Lm-1P
X
lim
t-> 0
F1 (t) О • • • О Fm-1(t)P - P
. . . lim
t-> 0
t
- L1P - • • • - Lm-1P Fl(t)P - P
<
X
t
- L1P
0.
X
Таким образом, все требования теоремы Чернова выполнены и ^ (¿) ~ Т4. Теорема доказана.
оо
x
t
t
Отметим, что если некоторые из полугрупп (Тд(¿))4>0 известны явно и если ||Тд (£) || < еак4 для некоторого > 0, то в формулах Фейнмана можно взять ТД (£) = Тд (£).
Пример 3 (возмущения шредингеровского типа). Пусть X = С^(К^); (Т4)4>0 — сильно непрерывная полугруппа с генератором (Ь, Бош(Ь)); (Т(¿))4>0 — семейство операторов, эквивалентное по Чернову полугруппе (Т4)4>0; с(^): К ^ К — непрерывная ограниченная сверху функция. Тогда функция е4с( ) является ограниченной при всех £ > 0. При этом оператор с умножения на функцию с(-) плотно определен в пространстве X (так как корректно определен, например, на множестве непрерывных финитных функций) и порождает сильно непрерывную полугруппу е4с операторов умножения на функцию е4с( ). Рассмотрим оператор Ь + с, где Бош(Ь + с) = Бош(Ь) П Бош(с) и (Ь + с)^(д) = Ь^^) + с(д)^(д) для всех € Бош(Ь + с). Предположим, что этот оператор порождает сильно непрерывную полугруппу (Т4Ь+С)4>0 на (такое предположение справедливо, например, если функция с(-) является ограниченной или если оператор умножения на функцию с(-) является Ь-ограниченным оператором). Тогда по теореме 3 справедлива формула Фейнмана
TtL+c = lim lenc о F(t/n)ln
В частности,
= lim [eПc о Tt/„]n,
если ||Tt || < eat для некоторого а > 0. Аналогичным образом показывается справедливость этих формул Фейнмана и в других банаховых пространствах.
Пример 4 (возмущения градиентного типа). Пусть X = Co(Rd); b(-): Rd ^ Rd — ограниченное непрерывное векторное поле; V = ^-, ..., — оператор Гамильтона. Рассмотрим оператор bV, где bV^(q) = b(q) ■ V^(q) для всех ^ G D(bV), и семейство (S(t))t>0 операторов на Co(Rd), определенных по формуле
S (%(q) = <^(q + tb(q)). (3)
Отметим, что в случае постоянного векторного поля b(-) = b G Rd семейство (S(t))t>0 совпадает с полугруппой сдвигов (eibv)t>0. Если векторное поле b(-) уже не является постоянным, то семейство (S(t))t>0 не является полугруппой, но эквивалентно по Чернову полугруппе etbv, порожденной оператором bV. Действительно, S(0) = Id, ||S(t)|| = 1 и для всех ^ G D(bV) П C2(Rd)
lim
t-> o
S(t)y - ^ - bV^
lim sup
q€Rd
<^(q + tb(q)) -^(q) - b(q) .V^(q)
t
<
< lim t( sup |b(q) ■ Hess <^(q + sb(q)) b(q)| VqeRd,se[0,i]
oo
0
d
где Hess p — гессиан p. Таким образом, — S(t) P = bVp для любой p g D(bV) П C2(Rd).
t=0
Отметим, что множество ^ € Д (ЬУ) П С2 (К^) является существенной областью определения оператора ЬУ.
Пусть теперь (Т4)4>0 — сильно непрерывная полугруппа на X = с генератором
(Ь, Бош(Ь)) и семейство (^(¿))г>0 эквивалентно по Чернову полугруппе (Т4)4>0. Рассмотрим оператор Ь + ЬУ, определенный соотношением (Ь + ЬУ)^(д) = Ь^(д) + Ь(д)для всех ^ € Бош(Ь + ЬУ) = Бош(Ь) ПБош(ЬУ). Предположим, что оператор Ь + ЬУ порождает сильно непрерывную полугруппу (Т4Ь+ЬУ)4>0 на С^(К^). Отметим, что предположение выполняется, например, если оператор ЬУ является Ь-ограниченным, т.е. Бош(Ь) С Бош(6У) и для всех € Бош(Ь) и некоторых Л € [0,1), 7 > 0 выполняется оценка
||bVp||x < Ä||Lp||x + Y||р||
X
В частности, bV является L-ограниченным для L = - (-А)а/2, а G (1, 2]. По теореме 3 семейство (S(t) о F(t))t>0 эквивалентно по Чернову полугруппе (TtL+6v)t>0 и справедлива формула Фейнмана
TtL+6v = lim [S(t/n) о F(t/n)!"
Пример 5 (лагранжева формула Фейнмана для возмущений полугруппы, порожденной лапласианом). Пусть X = Сте(Rd); &(•): Rd ^ Rd
— ограниченное непрерывное векторное поле; c(^) G Cb(Rd). Рассмотрим оператор Лапласа L = 1 А, Dom(L) = C,re(Rd). Замыкание этого оператора порождает полугруппу (Tt)t>0 [45, теорема 31.5], разрешающую задачу
df = 1 dt 2А
Коши для уравнения теплопроводности д- = 1 А/, где
Ttp(q) = (2nt)(-d/2) / e-^р(у) dy.
Тогда для полугруппы (Т4Ь)4>0, порожденной оператором Ь := ^ А + ЬУ + с, по теореме 3 справедлива формула Фейнмана
TLp(q0) = lim (^П") ^ 2 ^ /•••/" e" ^ ^e ^ P(?n) dq1... dq„.
Так как |x + b(x)t - y|2 = |x - y|2 + 2tb(x) • (x - y) + t2|b(x)|2, то
г^т , , л. f f n E c(qfc-1) - E b(9fc-i)^(9fc-i-9fc)
TtLp(q0) = lim •• en *=1 e *=i x
X e 2n fc?11 (9к 1)l p^/f (q0 - q1) . . .p^f (qra-1 - qra)p(qra) dq1 ... dq
Г 1x12 ^
где (ж) = (2п£)(-й/2) ехр| —— переходная плотность вероятности броуновского движения. Можно показать, что предел в правой части последней формулы совпадает с функциональным интегралом
Eq0
f c(es)ds f b(5s). -2 f |b(5s)|2ds
e0 e0 e 0 f (&)
(4)
по мере Винера, сосредоточенной на траекториях, начинающихся в точке д0. Таким образом, аппарат формул Фейнмана дает еще один способ доказательства формулы Фейнмана — Каца (4) и формулы Гирсанова. Кроме того, с помощью формул Фейнмана удается обобщить их на случай переменного коэффициента диффузии. Этому посвящен параграф 5 настоящей работы.
Пример 6 (лагранжева формула Фейнмана для динамики квантовой частицы в потенциальном и магнитном полях). Пусть X = Ь2 (М^). Рассмотрим оператор
Ь = -А,
Бош(Ь) = С£°(М^). Замыкание этого оператора порождает (по теореме 1Х.27 из [43] и
it
теореме Стоуна) (полу)группу (Т = е 2 Л)4>0, разрешающую задачу Коши для уравнения Шредингера г д = - 2 А/:
где правая часть понимается в регуляризованном смысле [43, формула 1Х.31]. Пусть Ь € М^, с(-) € СЬ(Е^). Тогда замыкание определенного на С(?°(МЙ) оператора Ь = ^ А — ЬУ — гс порождает (по теореме Х.22 из [43] и теореме Стоуна) сильно непрерывную (полу)группу (Ть)4>0, для которой по теореме 3 справедлива формула Фейнмана
TtV(qo) = lim [e-iПc о e-Пbv о Tt/ra]>(qo) =
Г Г -ic(qk i) i У 1 "k-qk-1 + ПЬ| 2 = lim (2nit)(-dn/2) / ..Je n e k=i 2t/n ^ ... dqn)
J J
причем интегралы снова понимаются в регуляризованном смысле. Отметим, что выражение в последней строке может быть интерпретировано как интеграл Фейнмана по траекториям в конфигурационном пространстве системы, квантование которой описывается гамильтонианом H = — - А — ibV + c:
r Г -i J[c(S(s))-1 |b|2]ds if b-d?(s)
TtL^(qo) = e-itH^(qo) = J e 0[ ( ( )) 211 ] e 0 ( V(£(¿))Ф*0(d£),
причем интеграл берется по псевдомере Фейнмана Ф®, сосредоточенной на траекториях,
t
начинающихся в точке q0, и эвристический «стохастический интеграл» f b ■ d£(s) требует
o
отдельного определения (ср. [15, теорема 5.4], [12]). Аналогичная формула Фейнмана остается справедливой и для более широкого класса магнитных и потенциальным полей b(-) и c(-) (ср. [12, 42]).
4. Формулы Фейнмана для полугрупп, порожденных мультипликативными возмущениями генераторов
Пусть Q — метрическое пространство; X = Cb(Q) — банахово пространство ограниченных непрерывных (числовых) функций на Q с нормой ||/||те = sup |/(q)|. (В качестве
qeQ
Q можно взять, например, метрический граф. Мы не предполагаем локальной компактности Q, так что в качестве Q можно рассматривать также и гильбертово пространство.) Пусть (Tt)t>0 — сильно непрерывная полугруппа на X с генератором (L, Dom(L)); функция а(-) : Q ^ [c1, c2] С (0, то) непрерывна на Q. Тогда aX С X. При этом оператор L, определенный по формуле Lp(q) = a(q)(Lp)(q), Dom(L) = Dom(L), также является генератором сильно непрерывной полугруппы [38, теорема 1], которую мы обозначим (Tt)t>0. Оператор L будем называть мультипликативным возмущением генератора L, а полугруппу (Tt)t>0 — полугруппой с мультипликативно возмущенным функцией a(-) генератором.
Теорема 4. Пусть (F (t))t>0 — сильно непрерывное семейство ограниченных операторов на банаховом пространстве X, эквивалентное по Чернову полугруппе (Tt)t>0. Тогда семейство операторов (F(t))t>0, определенное по формуле
F(t)p(q) = (F (a(q)t)p)(q) (5)
и действующее в пространстве X, эквивалентно по Чернову полугруппе (Tt)t>0 с мультипликативно возмущенным функцией a(-) генератором, т.е. локально равномерно по t > 0 в L(X) справедлива формула Фейнмана
Tt = lim [F(t/n)r
га^те
Доказательство. Если (F(t))t>0 — сильно непрерывное семейство (F(t))t>0 ограниченных операторов на банаховом пространстве X, эквивалентное по Чернову полугруппе (Tt)t>0, то ||F(t)|| < ekt для некоторого k > 0 и всех t > 0, а также F'(0)p = Lp для всех р G D, где D — некоторая существенная область определения генератора L. Покажем сначала, что это семейство действует из X в X. Пусть р G X = Cb(Q); проверим, что для каждого t > 0 функция F(t)p является ограниченной и непрерывной на Q функцией. Обозначим расстояние в Q между точками q и q0 через p(q, q0). Фиксируем t > 0 и q0 G Q. Тогда
lim |^(t)p(q) - ^(t)p(q0)| = lim |F(a(q)t)p(q) - F(a(q0)t)p(q0)| < < lim (|F(a(q)t)p(q) - F(a(q0)t)p(q)| + |F(a(q0)t)p(q) - F(a(q0)t)p(q0)^ < < lim ||[F(a(q)t) - F(a(q0)t)]p|U + lim |F(a(q0)t)p(q) - F(a(q0)t)p(q0)|.
p(?,9o)^0 p(q,qo)^0
Второе слагаемое в последней строке равно нулю, так как для каждого q0 £ Q оператор F(a(q0)t) действует из X в X, т.е. F(a(q0)t)p является непрерывной на Q (а значит и в точке q0) функцией. Кроме того,
lim \\[F(a(q)t) - F(a(qo)t)]^H« = lim, J|F(t- F(tqHL = 0
P(q,qo)—>0 т—To=a(qo )t
в силу сильной непрерывности отображения t ^ F(t). Таким образом, функция F(t)p является непрерывной. При этом,
||F(tHU = sup |F(t)p(q)| = sup |F(a(q)t)p(q)| < sup |F(a(qo)t)^(q)| <
< sup ||F(a(qo)tHU < sup ||F(a(qo)t)|| ■ |MU = eC2fct|MU < то.
Таким образом, функция F(t)p является также ограниченной, а значит, принадлежит пространству X. Кроме того, мы показали, что для любого t > 0 справедлива оценка ||F(t) || < < ec2kt. Заметим, что условие FF(0) = Id также выполнено.
Проверим сильную непрерывность семейства (F(t))t>0. Имеем:
lim\\F(t)p - F(to)pL = lim sup|F(a(q)t)p(q) - F(a(q)to)p(q)| <
qeQ
t—to 11 4 " 4 ' ' 1100 t—to
< ЦтэирН^(а(д)Ь)р — F(а(д)£0)р|| = Ит sup (сЬ)р — F. ^о дед ее[с1,с2Г
В силу сильной непрерывности семейства ^(¿))г>0 на [0, то) для любого фиксированного Т > 0 и любой фиксированной € X выполняется следующее: для любых ¿0 € [0,Т] и е > 0 существует такое $£(Ь0) > 0, что при Ь € [0,Т], удовлетворяющем условию — ¿0| < $£(Ь0), справедливо неравенство 11^^(¿)<р — F(Ь0)р|оо < е. Фиксируем ¿0, Т, при котором ¿0 € [0, Т], и е > 0. Рассмотрим Ц£о = (¿0 — ^^(¿0), ¿0 + ^^(¿0)). Тогда (Ц£о)4ое[0,т] — открытое покрытие компакта [0, Т]. Выберем из него конечное подпокрытие (Ц£о^П=г Пусть = шт{#£(Ь0(1)),..., ^е(Ь0(п))}. Тогда для любых ¿1,Ь2 € [0,Т], |Ь1 — ¿2| < 2, существует такое ¿0(к), что ¿1, ¿2 € (¿0(к) — £е(Ь0(к)), ¿0(к) + ^(¿0(к))). При этом
11^^(¿1)р — F(¿2)И|оо < 11^^(¿1)р — F(¿0(к))р|оо + (¿2)р — F(¿0(к))р|оо < 2е.
Теперь возьмем такое ^ > 0, что с2^ € (0, £е/2). Пусть Ь € [0, Т] таково, что — ¿0| < ^. Тогда для любого с € [с1,с2] выполнено неравенство |сЬ — сЬ0| < $е/2. Следовательно,
sup 11^^(сЬ)р — F(сЬ0)р|ос < 2е и, значит,
се[с1,с2]
Ит sup 11 ^^ (сЬ)р — F (с^0)р|оо = 0. се[с1,с2]
Тем самым сильная непрерывность семества (-Р(Ь))4>0 показана.
Найдем производную в нуле семейства ^(¿))4>0 на множестве О. Для любой р € О
lim
t-> о
F(t)p - p
— Lp
= lim sup
oo t—0 qeQ
F (a(q)t)p(q) — p(q)
— a(q)LP(q)
<
< lim sup
't-> о
qo ,q€Q
F (a(qo)t)P(q) —
— a(qo)LP(q)
lim sup
t-> о
qo€Q
F (a(qo)t)p — p
— a(qo)Lp
lim sup
t—0 ce[ci,c2]
F (ct)p — p
— cLp
< co • lim
^ о
F (z)P — P
— Lp
Таким образом, все условия теоремы Чернова выполнены, а значит, семейство (_Р(^))4>0 эквивалентно по Чернову полугруппе (Т4)4>0. Теорема доказана.
Замечание 9. Пусть операторы Т4 при всех £ > 0 удовлетворяют условию ||Т41| < при некотором к > 0. Тогда в качестве семейства операторов ^ (¿))г>0 на X, можно взять саму полугруппу (Т4)4>о:
= (Т^р)^, р € X, 5 € д.
Заметим, что такое семейство ^(¿))4>0 уже не является однопараметрической полугруппой. Однако по теореме 4 это семейство ^(¿))4>0 эквивалентно по Чернову полугруппе (Т)4>0.
Замечание 10. Утверждение теоремы 4 можно перенести на случай, когда банахово пространство X есть не Сь(д), а некоторое линейное пространство функций на д, полное относительно супремум-нормы
, = sup |p(q)|. При этом надо дополнительно
q€Q _
накладывать/проверять условие, что семейство операторов (F(t))t>0 действует из X в X. Предложение 1. Теорема 4 верна и в случае, если
X = CUQ):={p е Cb(Q): lim p(q) = о),
^ p(q,qo)—J
где q0 — произвольная фиксированная точка неограниченного множества Q.
Доказательство. Достаточно проверить, что lim (F(t)p)(q) = 0 для любой
p(q,qo)—^
p е X. Фиксируем p е X. Тогда
lim |F(t)p(q)| = lim |(F(a(q)t)p)(q)| <
< lim sup |F(a(q')t)p(q)| = lim sup |F(r)p(q)|.
p(q,qo)^^ q'eQ p(q,qo)^^ r e[cii,c2i]
Так как семейство (F(t))t>0 действует из X в X, то для любого t > 0 и любого е > 0 существует Re,t > 0, такое, что для любого q G Q, p(q0, q) > Re,t, выполняется неравенство
|F(t)p(q)| <е/2.
t
t
t
t
DO
0
t
z
DO
OO
Так как семейство (Т (¿))г>0 сильно непрерывно, то для любого е > 0 существует > 0, такое, что для любых г, г' € [с!^, с2£], удовлетворяющих условию |т — г'| < выполняется неравенство ||Т(г)р — Т(г')р||те < е/2 (это было проверено в ходе доказательства теоремы 4).
Фиксируем е > 0. Разобьем отрезок [с^, с2£] точками т0 = с^ < Т1 < ... < т^ = с2£ на части так, что тах |гк — гк-1| < При этом для любого г € [с1£, с2£] существует такое гк,
что |г — гк| < Пусть теперь Я£ = тах Д£ Т,. Тогда для любых ^ О, р(д, д0) > Д£, и
г € [с1£, с2£] имеем:
(г)р(д)| < (г)р(д) — Т(г*)р(д)| + (г*)р(д)| <
г г
< ||Т(г)р — Т(г*)р||те + (г*)р(д)| < 2 + 2 = е.
Таким образом, Ит вир (г)р(д)| = 0, что и доказывает предложение. Теорема доказана.
Следствие 1. Пусть дан марковский случайный процесс (Х4)4>0 с пространством состояний О и переходной вероятностью РЦд, ¿у). Предположим, что соответствующая процессу полугруппа (Т4)4>0, заданная соотношением
Т4р(д) = Ее9[р(Х4)] = / р(у)Р4(д, ¿у),
Я
является сильно непрерывной полугруппой на банаховом пространстве X, где X = Сь(О) или X = С^(О). Тогда семейство операторов (_Р(£))4>0, действующее по формуле
^Р^ = У р(у)ра(д)*(Му),
Я
эквивалентно по Чернову полугруппе (Т4)4>0 с генератором, мультипликативно возмущенным функцией а(-) и справедлива лагранжева формула Фейнмана
Т4р(д0) = ^Нт J ■■^ р(5п)Ра(до)*/п(?0, ¿51 )Ра(«1)*/п(91, ¿92 ) . . . Ра(«„_1)*/п(9п-1, ¿9п). (6)
Я Я
Это следствие вытекает из теоремы 4.
Замечание 11. Следствие 1 выполняется, в частности, в случае, если ^^^ — феллеровский процесс, т.е. О = Е^, (Т4)4>0 — сильно непрерывная, сохраняющая положительность сжимающая полугруппа на пространстве С^(М^).
Замечание 12. Мультипликативное возмущение генератора марковского процесса равносильно процедуре случайной замены времени [11, 28]. Отметим, что Р(¿,х,^у) = = Ра(х)4 (х, ¿у) в общем случае не является переходной вероятностью какого-либо случайного процесса. Тем не менее, если переходная вероятность Р4(д, ¿у) исходного процесса известна,
то формула (6) позволяет аппроксимировать не выражающуюся в явном виде переходную вероятность модифицированного процесса.
Пример 7 (формула Фейнмана для диффузии с переменным коэффициентом). Рассмотрим процесс броуновского движения в Генератором этого процесса является оператор Лапласа L = -Д. Плотность переходной вероятности задается гауссовской экспонентой
pfM (x) = (2nt)-d/2 expj-М2
Пусть (Tt)t>0 — полугруппа на C^ (Rd) с генератором L, мультипликативно возмущенным функцией а(-). Эта полугруппа соответствует диффузионному процессу с переменным коэффициентом диффузии i/a. По теореме 4 для любого р е C^ (Rd) справедлива формула Фейнмана
7M»)=ta / ... / (2na(qo)i/n)"d/2 exp{-^¡-g } X
Rd Rd
X (2na(qra-i)t/n)-d/2 expj-j^^-1—ГТ/Г p(qra) dqi... dqra.
Пример 8 (процесс типа Коши с переменным коэффициентом). Рассмотрим процесс Коши в Плотность его переходной вероятности задается формулой
pt(x) = r( d + 1 ^
,2 2У [п|х|2 + ¿2р+1)/2'
где Г(^) — гамма-функция Эйлера. Рассмотрим мультипликативное возмущение генератора процесса Коши функцией а(-). Тогда для соответствующей полугруппы (Т;);>о с мультипликативно возмущенным генератором по теореме 4 имеем:
Ttp(qo) = lim
гd + 1
2 2
a(qo)t/n
X
[(a(qo)t/n)2 + (n|qo - qi|)2](d+1)/2
_а(дп-1)*/п_ , ,
Х [(а^-^/п)2 + (п|^-1 - дга|)2](^+1)/2^ • • • ^
5. Формулы Фейнмана для параболического уравнения второго порядка с переменными коэффициентами
В этом параграфе мы обобщим результаты примеров 5 и 7 на более широкий класс переменных коэффициентов.
Пусть А(-): ^ ) — такое непрерывное отображение, что для любого х Е Е оператор А(х) является симметричным и положительно определенным. Пусть Да — дифференциальный оператор, действующий на дважды дифференцируемую на функцию р следующим образом:
(Да р)(х) :=1г(А(х) р(2)(х)), (7)
i>
где p(2)(x) := (Hess p)(x) — гессиан p в точке x G Rd. Пусть отображения &(•): Rd ^ Rd и c(-): Rd ^ R непрерывны. В этом параграфе мы будем рассматривать оператор L, определенный на множестве дважды дифференцируемых на Rd функций следующим образом:
Lp(x) := 2(AAp)(x) + b(x) ■ Vp(x) + c(x)p(x). (8)
Рассмотрим задачу Коши
df
■7^(t,x) = Lf(t,x), t> 0, x G Rd;
dtv ; v ; (9)
f (0, x) = fo(x), x G Rd.
При f0 G C^(Rd) будем искать такую функцию f: [0, то) х Rd ^ R, что f (t, ■) дважды дифференцируема на Rd при всех t > 0, f (■, x) дифференцируема на (0, то) при всех x G Rd и f (t, ■) G C^(Rd) при всех t > 0.
Рассмотрим семейство операторов (Fi(t))t>0 (на пространстве C^(Rd), как будет показано в дальнейшем), для которого F1(0) := Id и при всех t > 0 оператор F1 (t) определяется по формуле
Fi<"p(x) := ТЩ^Щр /exp(-)p(y) dy. „0)
Если A(-) — постоянная матрица, то семейство операторов (F1(t))t>0 совпадает с полугруппой (e2 аа )t>0. Если A(-) не является постоянной матрицей, то семейство (F1(t))t>0 уже не является полугруппой. При этом, если A(x) при всех x G Rd является диагональной матрицей, у которой все диагональные элементы равны a(x), где а(-) — некоторая непрерывная, ограниченная, отделенная от нуля функция, то F1(t)p(x) = (Ta(x)tp)(x), где Tt = e2А. Значит, по теореме 4 семейство (F1(t))t>0 эквивалентно по Чернову полугруппе (Tt = e2 аа)t>0. Как будет показано далее, эквивалентность по Чернову сохраняется и в общем случае.
Лемма 1. Семейство (F1(t))t>0 является непрерывным семейством сжатий на пространстве C^(Rd).
Доказательство. Линейность операторов F1 (t) очевидна. Кроме того, для любой функции p G Cc(Rd) С C^(Rd) функция F1(t)p G C^(Rd) в силу свойств гауссовской экспоненты (в случае как ограниченного, так и неограниченного отображения A(-)).
При этом, для любой p G Cc(Rd)
sup |F1 (t)p(x)| = sup
< ||p||^ sup
I -A 1(x)(x - y) ■ (x - y) I , s,
exp --- p(y) dy
л/det A(x)(2nt)dJ V 2t
Rd
1 '' /^-A~1(x)(x - y) ■ (x - y) ,
exp -ö- dy
<
Vdet A(x)(2nt)dJ V 2t
(11)
1
DO
Так как Сс(Ей) плотно в СО (К"), то по теореме об ограниченном линейном отображении [43] продолжается до сжатия на пространстве Со(К"). Это продолжение также действует по формуле (10).
Покажем теперь, что семейство (^\(£))4>0 является непрерывным. Достаточно показать, что Иш Ц^1 (£)р — ^(¿0)р||оо = 0 для каждой р е Сс(К") и любого £0 > 0, а далее воспользоваться (3е) -приемом. Для £0 > 0 равенство верно, так как все функции в формуле (10)
непрерывны по £ е (0, +то), р имеет компактный носитель и отображение
( — (А-1(х)(х — у),х — у)
у) ^ ехр^---
является непрерывным на [£1, £2] х К" х К", где £0 е (£1, £2) и [£1, С (0, то). Для р е Сс(К") имеем также
1 [ (— А-1(х)(х — у) ■ (х — у)\ , . ,, .
!-Ш0уаеПЩ^У ехр( —(— У) ( )р(у)= р<х)'
равномерно по х е К", так как для каждого фиксированного х е К функция
1 /— А-1(х)(х — у) ■ (х — у)
у ^ — ехр
V^det A(x)(2nt)d V
стремится к дельта-функции Дирака при t \ 0, а функция р ограничена и равномерно
непрерывна. Таким образом, limF(t)p = р в Coo(Rd) для всех р е Cc(Rd). Теорема
t\0
доказана.
Лемма 2. Для р е Cc2'7(Rd), 0 <7 < 1, имеем
Fi(t)p(x) = p(x) + ^Aa p(x) + o(t)
при t \ 0 равномерно по x е
Доказательство. Пусть c(t, x) := ((2nt)ddetA(x))-1/2, t > 0, x е Тогда
1 (Fi(t)p - p)(x) = ^ / exp(- A-1(x)(x - У) )(p(y) - p(x)) dy.
Так как p е C^7 (Rd), то ее разложение Тейлора с центром в точке x влечет для некоторого $ е [0,1] равенство
p(y) - p(x) = p(1)(x)(y - x) + 2p(2)($x + (1 - $)y)(y - x)2,
где p(k)(x) понимается как k-линейный функционал на (Rd)k, а обозначение p(k) (x)(y - x)k указывает на применение функционала p(k)(x) к вектору (y - x) е взятому k раз.
Таким образом, 1 - p)(x)
C(t'X) f e(p(1)(x)(y - x) + 1 + (1 - %)(y - x)2) dy
c(t,x) f f-A 1(x)(x - y) ■ (x - y)
exP ( -^——--- ) P(2) (x)(y - x)2dy +
2t J L \ 2t
Rd
+ ^ / exp( -A-1(x)(x ~y) ' (x - y)) (p(2)(^x + (1 - •%) - p(2)(x))(y - x)2dy.
Rd
Так как р(2) удовлетворяет условию Гельдера, второе слагаемое последнего равенства стремится к нулю при t \ 0 равномерно по x G Rd. И, как и в случае постоянной матрицы A, первое слагаемое совпадает с - AAp(x) для всех x G Rd. Теорема доказана.
Теорема5. Пусть A(-): Rd ^ L(Rd) — такое непрерывное отображение, что для любого x G Rd оператор A(x) является симметричным и положительно определенным. Предположим, что существует а, 0 < а < 1, для которого оператор (Aa, Cf'^R^)), действующий по формуле (7), замыкаем, причем это замыкание является генератором сильно непрерывной полугруппы (etAA)t>0 на C^(Rd). Тогда семейство (F1(t))t>0, заданное формулой (10), эквивалентно по Чернову этой полугруппе, т.е.
etAA = lim [F\(t/n)]n
в L(C^(Rd)) локально равномерно по t > 0. Кроме того, для любой р G C^(Rd) справедлива лагранжева формула Фейнмана
etAAp(qo) = lim / ... PA(t/n, qo, .. .pA(t/n, q„-b q„)p(q„) dq1... dq„, (12)
n^J J
Rd Rd
где
/. ч _1_ ( A-1(x)(x - y) ■ (x - y)
PA(t'x'y) := Vdet A(x)(2nt)d eX4--2t-
при x, y G Rd, t > 0.
Доказательство. Как следует из лемм 1 и 2, все условия теоремы Чернова выполнены. Поэтому утверждение данной теоремы вытекает из теоремы Чернова. Теорема доказана.
Замечание 13. При некоторых дополнительных предположениях [31, 32, 35] для полугруппы etAA справедлива следующая формула Фейнмана — Каца:
etAAp(qo) = Eq0 [p(Xt)] = i p(y)PtA(qo,dy), (13)
где Р!А(д0, ^у) = е ^у | Х0 = д0} — переходная вероятность диффузионного процесса X с переменной матрицей диффузии А(-). Функция рА(£,х,у) из теоремы 5 отнюдь не является плотностью переходной вероятности этого случайного процесса; его плотность не выражается в элементарных функциях, однако ее можно аппроксиммировать с использованием формулы Фейнмана (12) и формулы Фейнмана — Каца (13).
Из результатов теорем 3 и 5, а также примеров 3 и 4 вытекает следующее утверждение.
Теорема6. Пусть А(-): К" ^ — такое непрерывное отображение, что для любого
х е К" оператор А(х) является симметричным и положительно определенным. Предположим, что существует а, 0 < а < 1, для которого оператор (Да, С2'а(К")), действующий по формуле (7), замыкаем, причем это замыкание является генератором сильно непрерывной полугруппы (в!Лл)*>0 на Со(К"). Пусть &(•): К" ^ К" — ограниченное непрерывное векторное поле; с(-): К" ^ К ^ — непрерывная ограниченная сверху функция; Допустим, что существует 7, а < 7 < 1, при котором оператор (£, С^'7(К")), заданный формулой (8), замыкаем, причем это замыкание (£, Оош(Ь)) является генератором сильно непрерывной полугруппы (Т!)!>0 на Со(К"), разрешающей задачу Коши (9) для всех /0 е Бош(Ь). Рассмотрим семейство (^1(£))!>0 операторов, заданное формулой (10) и семейство (5(£))!>0, заданное формулой (3). Тогда верны следующие утверждения:
1) 5(£) о е!с о ЗД) - Т*, 3) е!с о 5(£) о ЗД) - Т, 5) ЗД) о 5(£) о е!с - Т,
2) 5(£) о ЗД) о е!с - Т*, 4) ЗД) о е!с о 5(£) - Т, 6) е!с о ЗД) о 5(£) - Т*.
Замечание 14. Если с(-): К" ^ К ^ — непрерывная ограниченная сверху функция, оператор умножения на эту функцию плотно определен в пространстве Со (К") (так он как корректно определен, например, на множестве непрерывных финитных функций) и порождает сильно непрерывную полугруппу е*с операторов умножения на функцию е^-.
Согласно теореме 6, построено шесть различных семейств, эквивалентных по Чернову полугруппе (Т)^. Построим еще одно семейство, пригодное и в случае, если векторное поле &(•): К" ^ К" не является ограниченным.
Предположение 1. Для заданного непрерывного отображения А(-) : К" ^ £(К"), такого, что для любого х е К" оператор А(х) симметричный и положительно определенный, непрерывного векторного поля Ь(-) и непрерывной ограниченной сверху функции с(-) предполагаем, что существует а, 0 < а < 1, при котором оператор (£, С2'"^")), заданный формулой (8), замыкаем, причем это замыкание (£, Бош(Ь)) является генератором сильно непрерывной полугруппы (Т)^ на Со (К"), разрешающей задачу Коши (9) для всех /0 е Бош(Ь).
Замечание 15. Предположение 1 выполнено, например, при с(-) е С1(К"), &(•) е е С2(К") и А(-) е ), если оператор Да является равномерно эллиптическим (т.е.
существует такая константа 7 > 0, что для всех х е К" и £ е К" \ {0} выполнено неравенство £ ■ А(х)£ > 7|£|2). Этот факт следует из доказательства теоремы 2.1.43 (см. [34, том II]), основанного на теореме Хилле — Иосиды — Рэя. В случае неограниченных коэффициентов
вопросы существования соответствующей (феллеровской) полугруппы рассматриваются, например, в работе [18].
Определим семейство операторов (Р2(Ь))4>0 на множестве Сс(Е^) следующим образом: Т2(0) := Ы и при Ь > 0 оператор Т2(Ь) задается формулой
^2(Ь)р(ж) := . 1 х
л/ёе1 А(ж)(2пЬ)^
х У ехр(-А-1(Х)(Х — У) ' (Х - У)) ехр(—А-1Ь(ж) ■ (ж - у))р(у) ¿у, (14)
где А-1Ь(ж) := А-1(ж)Ь(ж), ж е №*.
Лемма 3. Для всех р е (М^), 0 <7 < 1, при Ь \ 0 равномерно по ж е К выполняется равенство
^2(Ь)р(ж) = р(ж) + 2Д^р(ж) + Ьб(ж) ■ Ур(ж) + 2(А-1Ь(ж) ■ 6(ж))р(ж) + о(Ь).
Доказательство. Пусть р е С?'7(К^) и выполнено предположение 1. Тогда для любого ж е К функция = ехр (—А-1Ь(ж) ■ (ж — у)) р(у) переменной у принадлежит
множеству Сс2'7 (М^). Кроме того,
^2(Ь)р(ж) = (ЗД^Хж). Таким образом, по лемме 2 при Ь \ 0 равномерно по ж е К
(Ь) р(ж) = ^Х'^(ж) + ^ (Да^Х'дж) + о(Ь) =
(ж) • (ж •) )р |(ж) + -Да Гехр(-А-1Ь(ж) • (ж - •) )Р 1(ж) + о(
(exp(A-1b(x) ■ (x - ^))p) (x) + - exp(-A-1b(x) ■ (x - ^))p) (ж) + o(t) = p(x) + -tr(A(p(2) + 2A-1b 0 Vp + (A-1 b 0 A-1b)p)) (x) + o(-) =
= p(x) + 2Aap(x) + tb(x) ■ p(x) + 2(A-1b(x) ■ b(x))p(x) + o(t).
Теорема доказана.
Рассмотрим теперь семейство операторов (F(i))t>0, для которых для любой p G (Rd) F (t)p(x) = et[c(x)-1(x)b(x)^(x)]F2(t)p(x) =
_ exp(-c(x)) f expf-A-1(x)(x - y +-b(x)) ■(x - y + -b(x)A p(y) (15)
л/ёе1 А(ж)(2пЬ)^У \ 2Ь
Отметим, что операторы Т(Ь) отличаются от операторов в формуле 3) теоремы 6 тем, что у последних коэффициенты А(-) и А-1() вычисляются в точке ж + ЬЬ(ж), а не в точке ж.
Теорема 7. Пусть выполнено предположение 1. Тогда семейство (F(t))t>o операторов, определенных по формуле (15), эквивалентно по Чернову полугруппе (Tt)t>o на пространстве X = C^(Rd), разрешающей задачу Коши (9) для всех /o G Dom(L), т.е.
Tt = lim (F(t/n))n
в L(C^(Rd)) локально равномерно по t > 0. Таким образом, для любой р G (Rd) справедлива лагранжева формула Фейнмана
Ttp(qo) = Jim J ••• J exp^ - ^^ A-1(qfc-i)b(qfc-i) • (qfc-i - qfc x
X exp( П ^ c(qfc-iH ex^ -^ A 1(qfc-i)6(qfc-i) • b(qfc-iH x ^ fc=i ' ^ fc=i '
* ^ с(?к-1ЛехР *
/г —'
х <р(дп)рл(*/га,до,д!) ■ ■ ■ дга-Ь д„) ¿91.. (16)
где при х, у е Е^, * > 0
/. ч _1_ / А-1(х)(х - у) ■ (х - у)'
Замечание 16. При некоторых дополнительных ограничениях на коэфициенты А(-), &(•), с(-) [31, 32, 35] для полугруппы (Тг)г>0 справедлива формула Фейнмана — Каца
Ttp(qo) =
(17)
о
где Е90 — математическое ожидание выпущенного из точки 90 диффузионного процесса (6)*>0 с переменной матрицей диффузии А(-) и сносом &(•). Таким образом, формула Фейнмана (16) дает пригодные для непосредственных вычислений аппроксимации функционального интеграла в формуле Фейнмана — Каца (17). Кроме того, с учетом формул (12) и (13) можно показать, что выражение в правой части формулы (16) совпадает со следующим функциональным интегралом:
t t
Ttp(qo) = Eq0
exp / c(Xr) dr exrn / A-i(Xr) b(Xr) • dXr ) I x
o
t
1
x exp( — I A-(Xr)b(Xr) • b(Xr) dr )/o(Xt)
2 _
o
(18)
где Е90 — математическое ожидание диффузионного процесса (Хг)г>0 с переменной матри-
г
цей диффузии А(-) и без сноса, а стохастический интеграл / А-1(ХТ)Ь(ХТ) ■ ¿ХТ есть стоха-
t
,-i/ ^) o
стический интеграл Ито. Теперь, приравнивая функциональные интегралы в формулах (17)
t
и (18), получаем аналог формулы Гирсанова — Камерона — Мартина — Реймера — Мару-ямы (о преобразовании меры, порожденной случайным процесом броуновского движения, при добавлении сноса) для случая диффузионных процессов с переменным коэффициентом диффузии.
Перейдем к доказательству теоремы.
Во -первых, Т(Ь)р е С^(М^) для любой функции р е Сс(М.^). Действительно, если
||А(ж)||ма1(^} ^ ТО при |ж| ^ ТО,
ТО
exp( tmaxc(x)
|F(t)p(x)| <-, x€Rd — I dy ^ 0, Ы^то.
| Л < ^det A(x)(2nt)d J У 11
supp (
Если существует такой путь Г, что при x Е Г, |x| ^ то имеем ||A(x)||Mat(dxd) = O(1), то существует такая константа k > 0, что для любых x Е Г, y Е t > 0 выполняется неравенство
( A-1(x)(x + tb(x) - y) ■ (x + tb(x) - y)\ ( k|x + tb(x) - y|2
ехш---- < ехш —
2Ь 2Ь
При этом существует такое Ь* > 0, что для любого Ь е [0,Ь*] имеем |ж + ЬЬ(ж)| ^ то при |ж| ^ то. Тогда по свойствам гауссовской экспоненты для всех Ь е [0,Ь*] при ж е Г, |ж| ^ то выполняется оценка
|F(f)p(x)|< ^det A(x)(2?rt)d J ex4--2t-I dy ^ 0.
supp (
В этом случае положим F(t) = F(t*) для всех t > t*. Во-вторых, для любой p Е Cc(Rd)
sup |F(t)p(x)| < ex^l t max(c(x))
xP»r V x€Rr
x sup
1 f /-A-1(x)(x - y + tb(x))(x - y + tb(x))
exp i dy
Л(ж)(2пЬ)<*У V 2Ь
= ехш Ь тах(с(ж))
V жек^
Значит, по теореме об ограниченном линейном отображении операторы Т (Ь) можно продолжить с множества Сс(К^) до линейных ограниченных операторов из С^(Е^) в причем
||Т(Ь)|| < ехр(Ь тах(с(ж))).
V жек^ /
Далее, по лемме 3 для любой р е С,?'7(Е^), 0 <7 < 1, при * \ 0 равномерно по х е Е выполняются соотношения:
^ (*)р(х) = ег[с(х)-1 А-1(ж)ь(ж)-ь(ж)]^2(%(х) =
= ег[с(х)-2 л-^ж.нн] ^р(х) + * Да^(х) + *Ь(х) ■ Ур(х) + * (А-1Ь(х) ■ Ь(х))р(х) + о(*)) =
= р(х) + 22Дар(х) + *Ь(х) ■ Ур(х) + с(х)р(х) + о(*).
Таким образом, все условия теоремы Чернова выполнены. Теорема доказана.
Замечание 17. Можно показать, что при надлежащих условиях на коэффициенты А(-), &(•), с(-) оператор Ь = —¿Я, где
Яр(х) := — ^(ДАр)(х) — ¿Ь(х) ■ Ур(х) + с(х)р(х),
с подходящей областью определения порождает сильно непрерывную полугруппу (и даже группу) е-гШ на пространстве Ь2(Е^) и семейство (^(*))г>0 операторов, таких, что для любой р е Се(Е^)
ехР(—¿¿с(х)) [ ^.Л-1(х)(х — у — ¿6(х)) ■ (х — у — ¿6(х))\
^(г)р(х) = ТШЩЗрУ ехЧ'-2г-)р(у)<гу
(ср. с формулой примера 6 и формулой (15)) эквивалентно по Чернову этой полугруппе (ср. [12, 42]).
6. Формулы Фейнмана для решения начально-краевых задач
Существуют различные подходы к получению формул Фейнмана для начально-краевых задач. С одной стороны, начиная с решения начально-краевой задачи для более простого уравнения, можно применять теоремы 3 и 4 об аддитивных и мультипликативных возмущениях для получения решения исходной задачи (см., например, [7]). С другой стороны, начиная с решения задачи Коши для исходного уравнения, можно «подправлять» соответствующее семейство эквивалентных по Чернову операторов таким образом, чтобы выполнялись краевые условия (см., например, [21, 10]). Опишем второй подход для решения краевой задачи Коши — Дирихле для некоторого класса уравнений.
Пусть С — ограниченная область в Е^, X = С^(Е^) и дифференциальный оператор (Ь, Бош(Ь)) является генератором сильно непрерывной полугруппы (Тг)г>0 на X. Пусть также дано семейство операторов (^(*))г>0 на X, эквивалентное по Чернову полугруппе (Тг)г>0, причем Ц^1 (*)|| < е"г для некоторого а е Е и всех * > 0, ^'(0)р = Ьр для всех р е О с Бош(Ь), где О — некоторая существенная область определения оператора Ь. Рассмотрим следующую задачу Коши — Дирихле:
М) = Ь/(*,9), *> 0, 9 е С;
/(0,9) = /0(9), 9 е С; (19)
/(*,9) = 0, 9 е дС.
Пусть Y = C0(G) := (p G C(G), lim p(q) = 0} — пространство непрерытнык функций, обращающихся в нуль на границе области. Это банахово пространство с нормой
|у = вирхес |/(х)|.
Предположение 2. Предполагаем, что существует сильно непрерывная полугруппа (Т°)4>0 на пространстве У, разрешающая задачу Коши — Дирихле (19), т.е. если (Ь0, Бош(Ь0)) — генератор полугруппы (Т°)4>0, то для любой р е Бош(Ь0) функция /(*, •) = р е У является решением задачи (19) с начальным условием /0 := р.
Условия существования полугруппы, разрешающей начально-краевую задачу (19), могут быть найдены, например в [37, теорема 3.2.5 и следствие 3.1.4].
Множество С0(С) можно различными способами вложить в пространство Се(М^) непрерывных на Е функций с компактными носителями.
Предположение 3. Предполагаем, что существует вложение Е: С0(С) ^ Сс(Ей), для которого:
(Е0) каждая функция р е С0(С) продолжается до функции Е(р) е Се(Е^), т.е. Е(р) =
(Е1) вложение линейно, т.е. для всех функций р1, р2 е С0(С) и чисел а, Ь е Е справедливо равенство Е(ар1 + Ьр2) = аЕ(р1) + ЬЕ(р2);
(Е2) вложение сохраняет норму, т.е. для любой р е У = С0(С) справедливо равенство
||р||у =вир |Е(р)(9)| = |Е(р)||х;
(Е3) вложение сохраняет существенную область, т.е. если О | с — множество сужений на С функций из О с Бош(Ь), то для любой р е Бош(Ь0) П О|с ее продолжение Е(р) е О.
Замечание 18. Пусть, например, Ь = Д, О = С^Е^). Тогда требование (Е 3) равносильно требованию сохранения при вложении гладкости до второго порядка включительно.
Предположение 4. Предполагаем, что множество
D0 := Dom(L0) П [C1 (Rd) П D]
G
является существенной областью определения генератора Ь0 полугруппы (Т°)¿>0, разрешающей поставленную задачу Коши — Дирихле (19).
Замечание 19. Пусть оператор Ь задан формулой (8), причем коэффициенты являются функциями класса С2'а(Е^) для некоторого а е (0,1), матрица А(9) симметрична при всех 9 е Е и оператор Да равномерно эллиптичен. Рассмотрим в качестве О множество С2'«(^) с СО(Е^) = X, тогда О0 = {р е С2'а(С) : р, Ьр е У}. Если граница области С является гладкой класса С4'а, то предположения 2, 3 и 4 верны (см. [16]).
Пусть ^ : (0, то) ^ (0, то) гладкая функция, монотонно убывающая к нулю при * \ 0, причем = о(*) (например з(*) = к аг^ *2, где к е (0, diaш(С)/п)). Рассмотрим множество С8(4) С С, определенное по формуле С^ = {9 е С: dist(9,дС) > з(*)}. Рассмотрим теперь семейство (^(¿))¿>0 бесконечно дифференцируемых функций Е ^ [0,1], та-
ких, что
1 q Е Gs(t); 0, q Е Rd \ G,
и lim ||^s(i) - = 0 для любого t* > 0. Отметим, что при t \ 0 функции ^s(t)
поточечно сходятся к индикатору области G.
Рассмотрим семейство (Fo (t));>0 операторов, действующих по формуле
F0(t)p(q) = { Id, t = 0; (20)
I ^s(t)(q)[F(t)E(p)](q), t> 0.
где семейство (F(t));>0 эквивалентно по Чернову полугруппе (T;);>0, разрешающей задачу Коши в пространстве X, соответствующую нашей начально-краевой задаче (19). Это семейство операторов действует в пространстве Y. Действительно, если p Е Y = G0(G), то E(p) Е Gc(Rd) С G^(Rd) = X, F(t)E(p) Е X, ^(;)[F(t)E(p)] Е Y.
Теорема 8. Пусть выполнены предположения 2, 3 и 4. Тогда семейство (Fo(t));>0, заданное формулой (20), эквивалентно по Чернову полугруппе (T;o);>0, разрешающей задачу Коши — Дирихле (19).
Доказательство. Проверим выполнение условий теоремы Чернова. Так как Fo (0) = Id по определению, то необходимо проверить сильную непрерывность семейства (Fo(t));>0, оценить норму операторов Fo(t) и найти (сильную) производную Fo(t) в нуле на существенной области Do определения генератора Lo. Итак,
||F„(i)|| = sup "F;<<)p"y = sup SUP'^c|*-(;^<q)'F(t)E(p)](q)| < sup "F< «-.
||Е (р)||х ||р||х
Теперь проверим сильную непрерывность семейства (Т0(Ь))^>о. Во-первых, для любой
р е Y
|То(Ь)р - р^у = Ит8ир|^)(9)[Р(Ь)Е(р)](9) - р(9)| =
= Иглвир|^в(4)(д)([Т(Ь)Е(р)](д) - Е(р)(д)) + р(д)[^(*)(д) - 1]| <
деС
< 11п° ||Т(Ь)Е(р) - Е(р)||х + Игл вир |р(?)| = 0
д е с\св(4)
в силу сильной непрерывности семейства (Е(Ь))1>° на пространстве X и равномерной непрерывности р на компакте С \ С^). Во-вторых, для любого Ь* > 0 и для любой р е Y
Ит ||Е0(Ь)р - Е0(Ь*)р||у = Ит вир^М^(Ь)Е(р)](д) - ^(г)(д)[Т(Ь*)Е(р)](д)| =
дес
= Ит8ир|^в(4)(д)([Т(Ь)Е(р)](д) - [Е(Ь*)Е(р)](д)) + (^(9) -^(г)ШЕ(Ь*)Е(р)](д)|< дес
< Ит ||^в(4) ||у ■ || [Е(Ь)Е(р)] - [Е(Ь*)Е(р)] ||х + ||Е(Ь*)Е(р) ||х ■ ||^(1) - ^ *) ||у = 0.
Пусть теперь р G Do. Найдем
lim
t-> о
F0 - р
- L0p
= lim sup
Y q€G
F0(t)p(q) - p(q)
- LoP(q)
Отметим, что
F0(t)p(q) - p(q)
- LoP(q)
[F(t)E(p)j(q) - p(q)
- LP(q)
q G Gs(t);
^s(t)(q)[F(t)E(p)](q) - p(q)
- LP(q) , q G G \ Gs(t);
q G dG.
При этом
sup
?eGs(t)
[F(t)E(p)](q) - E(p)(q)
- LP(q)
<
F(t)E(р) - E(р)
- L(E(р))
-»■ 0
X
при t ^ 0. Кроме того,
sup «eG\Ge(t)
^e(t)(q)[F(t)E(р)](q) - р(q)
= sup
q€G\Gs(t)
< sup
9eG\G,(t)
- ¿р^
^s(t)[F(t)E(р)] - + ^ «(4)р р
^s(t)(q)
F(t)E(р)(q) - E(р)(q)
(q) -
- L(E(р))(q)
<
+
<
+ sup
q€G\Gs(t)
[F(t)E(р)] - E(р)
L(E(р))(^(^) - 1]
+ sup
q€G\Gs(t)
р(q)
^s(i)(q) - 1
t
t
- ОДр))
+ sup |L(E (р))^)| + sup
X q€G\Gs(t) ?eG\Gs(t)
<
|р(д)! t
при t ^ 0.
Действительно,
[F(t)E(р)] - E(р)
- L(E(р))
0
X
при * ^ 0, так как Е(О0) С О и при всех р е О в предположениях теоремы
lim
t-> о
F(t)р - р
- Ьр
0.
X
Учитывая, что L(E(р))(q) = 0 при q G ÖG для любой р G Do С Dom(Lo), в силу равномерной непрерывности функции L(E(р)) на компакте G \ Gs(t), получаем
lim sup |L(E (р))(q)| = 0.
Далее, для любого q G G \ Gs(t) существует хотя бы одна точка xq G ÖG, такая, что
dist(q,xq) = dist(q,dG) < s(t) = o(t) при t ^ 0. Так как р G D0, то E(р) G C(Rd) и по формуле Тейлора р(q) = E(р)(q) = E(р)(х) + E^)(1) (z)(q - ), где z = + (1 -
t
t
t
t
t
0
t
t
t
t
t
0
t
t
для некоторого ê G [0,1]. Следовательно,
|p(q)b s(t)
q€G\Gä(t)
sup ||E(p)(1)||x - 0
при t — 0. Таким образом, lim
t-> о
F0(i)p - p
--L0p
t
= 0, и, значит, выполнены все
Y
требования теоремы Чернова для семейства (Fo(t))t>o. Следовательно, по теореме Чернова это семейство эквивалентно полугруппе (Tto)t>0, разрешающей задачу Коши — Дирихле (19). Теорема доказана.
Пример 9. Пусть оператор L задан формулой (8), причем коэффициенты являются функциями класса C2'a(Rd) для некоторого а £ (0,1), матрица A(q) симметрична при всех q £ и оператор Да равномерно эллиптичен. Пусть также выполнено предположение 1. В качестве семейства (F(t))t>o операторов на пространстве X рассмотрим семейство, заданное формулой (15), а в качестве D — множество C2'a(Rd). Следовательно, Do = = (р £ G2'a(G): р, Lp £ Y}. Пусть граница области G является гладкой класса С4'а (тогда предположения 2, 3 и 4 верны, см. [16]). Следовательно, по теоремам 7 и 8 решение задачи Коши — Дирихле (19) может быто получено по следующей формуле:
f (t, qo) = TOfo(qo) = lim [F0(t/n)]nfo(qo) =
^(t/n)(qfc-i) I exp ( - У2 A-1b(qfc-i) ■ (qfc-i
Jim/ f (П ^(i/n)(qfc-i)^ exp^ - A 1b(qfc-i) • (qfc-i - qk x
x exp( n ^ c(qfc-iM exp( -Ê A-1(qfc-i)b(qfc-i) • b^-in x ^ fc=i ' ^ fc=i '
x E(p)(qra)pA(i/n,qo,qi) ■ ■ ■ p^(t/n, q„-b q„) dqi... dq„, (21)
Упростим последнее выражение. Пусть xG индикатор области G, т.е. xG(x) = 1 при x £ G, XG(x) = 0 при x £ G. Ни скорость сходимости s(t) ^ 0 (если эта скорость не медленнее, чем o(t)), ни выбор семейства )t>0, аппроксимирующего функцию xg при t ^ 0, не меняют предел в формуле (21). Так как функции ^s(t) гладкие и имеют компактные носители в G, а функция E(р) непрерывна и также имеет компактный носитель K С Rd, то в формуле (21) фактически производится интегрирование непрерывной функции по компакту K x G 1 С (Rd)n. Пусть теперь t £ (0,T] для некоторого фиксированного 0 < T < то. Так как ^s(t) ^ xG при t ^ 0, то можно выбрать s(t), t £ (0,T], так, что при всех n £ N выполняется неравенство
n i t
ys(t/n)(xo) / , exp ^A-1b(qfc-i) ■ (qfc-i - qfcM exp ( -2^c(qfc-i)l x
k—1 / \n K Gn-1 i i
^s(t/n)(xo^ У exp^- ^ A 1 b(qfc-i) • (qfc-i - qfc)^exp^ ^^ c(qfc-i)^
K G"-1 k=1 k=1
x exp( - 2~ ^ A-1(qfc-i)b(qfc-i) • b(qfc-in (П^/«)^-^ - xg^-OH x ^ fc=1 ' V=2 '
fc-i)b(qfc-i) • b(qfc-i ) I I J_J_('^s(t/n)(qfc-i) - XG(qfc-i
' V=2
x E(p)(qra)pA(t/n,qo,qi) • • • ^¿(t/n, q„-i, qra) dqi . . . dq„
1
< -. n
Тогда из формулы (21) вытекает, что
Ttp(qo) = J еХР(- S A 1b(qfc-i) ■ (qfc_i - qfcX
K G"-1 k=1
X exp( - 2~ ^ A-1(qfc-i)b(qfc-i) ■ b(9k-iH exp( -П ^ c(qfc-i) j ^ n fc=i ' fc=i '
X E(p)(qra)pA(t/n, qo, qi) ■ ■ ■ PA(t/n, q„_i, q„) dqi... dqra (22) равномерно по (x0, t) G G X [0, T] для любого T > 0. Кроме того,
exp I - У" A-ib(qfc_i) ■ (qfc-i - qfc ) ) X
G" K G"-1
fc=i
X ex^(-2" ^ A i(qfc_i)b(qfc_i) ■ 6(qfc_i)jexp( n. ^ c(qfc_i)j X ^ n fc=i ' fc=i '
X E(p)(qra)pA(t/n,qo,qi) ■ ■ ■ PA(t/n,qra_i,qra) dqi . . . dq,
<
< ex^ tmax{c(y)} n sup exp(A i(y)b(y) ■ (z - y)) I X
v yeG / VyeG,zeK\g /
X sup pA(t/n,y, z)|E(p)(z)|dz — 0 yeG J
K\G
при n —► то. Действительно, пусть e > 0 и G<s — ^-окрестность области G в 6 > 0. Так как функция E(p) непрерывна на и обращается в ноль на dG, то существует 6 > 0, такое,
что |E(р)| < е/2 на G5 \ G. Тогда
sup pA(t/n,y,z)|E(p)(z)| dz <
yeG ./ K\G
< sup pA(t/n,y,z)|E(p)(z)| dz + sup pA(t/n,y,z)|E(p)(z)| dz <
yeG J yeG ./
g5\g k\g5
< e + sup |E(p)(z)| sup pA(t/n,y,z)dz. (23)
2 zeK yeG ./
к\Gi
В силу быстрого убывания гауссовской экспоненты ра, для выбранного выше 6 > 0 существует N G N такое, что при всех n > N последнее слагаемое в формуле (23) мажорируется константой е/2. Таким образом, доказано следующее утверждение.
Предложение 2. Пусть оператор L задан формулой (8), причем коэффициенты являются функциями класса G2'a(Rd) для некоторого a G (0,1), матрица A(q) симметрична при всех q G и оператор Да равномерно эллиптичен. Пусть выполнено предположение 1. Кроме того, пусть граница области G является гладкой класса G4'а. Тогда решение задачи Коши —
Дирихле (19) может быто получено по следующей формуле Фейнмана:
4°/оЫ = Ит I • • I ехр( - У^ • - 9кЛ х
/& 9о) = УоЫ = ехр(- 5 А 1&(9к-0 • (9к-1 - 9к
С С к=1
х ехр( п 5 ехр(-5 А-1(?к-1)ь(?к-1) • х
X 9о, 91) • • • 9«-1,9«) ¿91... ¿9«, (24)
Замечание 20. Пусть /о е Бош(Ь°) и 0 < Т < то. Тогда решение задачи Коши — Дирихле (19) при £ е [0, Т] может быть представлено с помощью формулы Фейнмана — Каца (см. [52], лемма 3.1, теорема 3.3):
4
/ (£,9о) = Е90
£ < тс
9о е С, £ е (0,Т), (25)
о
где Е90 — математическое ожидание начинающегося в точке 9 е С диффузионного процесса (6)о«т с переменной матрицей диффузии А(), сносом &(•) и поглощением на границе области; тс — время первого выхода (^)о«т из С. Таким образом, конечномерные интегралы в формуле Фейнмана (24) дают аппроксимации функционального интеграла в формуле Фейн-мана — Каца (25), причем эти аппроксимации содержат интегралы только от элементарных функций и не содержат переходные плотности данного диффузионного процесса.
Замечание 21. Аналогичным образом, решение поставленной задачи Коши — Дирихле (19) может быть найденно по формулам Фейнмана, построенным с помощью шести семейств операторов, введенных в теореме 6.
Замечание 22. Теорема 8 может быть обобщена в следующих направлениях. Во-первых, можно рассматривать неограниченные области в К и по-прежнему использовать вложение Со(С) в Сс(К^) со свойствами из предположения 3. Во-вторых, в качестве объемлющего банахова пространства X можно брать не пространство (К^), а какое-либо другое пространство функций (с нормой || • — такое, в которое пространство У = (Со(С), || • вкладывается с сохранением свойств, указанных в предположении 3, причем оператор Ь с надлежащей областью определения является генератором сильно непрерывной полугруппы. Более того, можно рассматривать другие пары банаховых пространств X и У, например У = Ь2(С) С X = Ь2(М^). В-третьих, аналогичным образом можно рассматривать и начально-краевые задачи в области риманова многообразия. В-четвертых, теорема 8 может быть обобщена и на случай нелокальных псевдо-дифференциальных операторов Ь.
Заключение
В настоящей работе описывается подход к решению начальных и начально-краевых задач для эволюционных уравнений, основанный на представлении соответствующих эволюционных полугрупп с помощью формул Фейнмана. В статье обсуждаются некоторые методы
построения формул Фейнмана для различных эволюционных полугрупп, приведены конкретные примеры решения эволюционных уравнений. В частности, получены формулы Фейнмана для эволюционных полугрупп, порожденных мультипликативными возмущениями генераторов некоторых исходных полугрупп. При этом рассматриваются полугруппы на некотором банаховом пространстве непрерывных функций, определенных на произвольном метрическом пространстве; формулы Фейнмана строятся с помощью семейств операторов, эквивалентных по Чернову исходным, невозмущенным полугруппам. Настоящий результат обобщает некоторые результаты работ [4] и [23]. Подход к построению формул Фейнмана для полугрупп с мультипликативно и аддитивно возмущенными генераторами иллюстрируется на примерах задачи Коши для уравнения Шредингера, аппроксимации переходных вероятностей некоторых марковских случайных процессов. Далее в работе рассматривается более широкий класс аддитивных и мультипликативных возмущений конкретного генератора — оператора Лапласа. При этом выводятся формулы Фейнмана для решения задачи Коши для параболического уравнения второго порядка с неограниченными переменными коэффициентами. Кроме того, в статье описывается метод построения формул Фейнмана для решения начально-краевой задачи Коши — Дирихле для дифференциального уравнения параболического типа. Метод также иллюстрируется на примере параболического уравнения второго порядка с переменными коэффициентами. Настоящие результаты обобщают некоторые из результатов работы [21]. В статье также обсуждаются некоторые формулы Фейнмана — Каца и интегралы Фейнмана, совпадающие с полученными формулами Фейнмана.
Работа выполнена при поддержке гранта Министерства образования и науки Российской Федерации (грант № 14.B37.21.0370) и гранта Президента Российской Федерации (грант МК-4255.2012.1).
Список литературы
1. Богачев В.И., Смолянов О.Г. Действительный и функциональный анализ: университетский курс. Москва-Ижевск: НИЦ «Регулярная и хаотическая динамика», 2011. 728 с.
2. Бутко Я.А. Представления эволюционных полугрупп с помощью формул Фейнмана и интегралов Фейнмана по траекториям в фазовом пространстве // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2012. №2. Режим доступа: http://technomag.bmstu.ru/doc/315838.html (дата обращения 01.02.2014)
3. Бутко Я.А. Формулы Фейнмана и функциональные интегралы для диффузии со сносом в области многообразия // Математические заметки. 2008. Т. 83, №3. С. 333-349. Б01: 10.4213/т2т3772
4. Бутко Я.А. Формула Фейнмана для полугрупп с мультипликативно возмущенными генераторами // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2011. № 10. Режим доступа: http://technomag.bmstu.ru/doc/239563.html (дата обращения 01.02.2014).
5. Бутко Я.А., ГротхаусМ., СмоляновО.Г. Формула Фейнмана для параболического уравнения второго порядка в области // Доклады Академии Наук. 2008. Т. 421, № 6. С. 727-732.
6. Бутко Я.А., Дурягин А.В. Формулы Фейнмана для семейства параболических уравнений, соответствующих тау-квантованию квадратичной функции Гамильтона // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2011. №11. Режим доступа: http://technomag.bmstu.ru/doc/251251.html (дата обращения 01.02.2014).
7. Бутко Я.А., Морозов А.В. Представление решения задачи Коши — Неймана для параболического уравнения на полупрямой с помощью лагранжевой формулы Фейнмана // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2011. №11. Режим доступа: http://technomag.bmstu.ru/doc/246219.html (дата обращения 01.02.2014).
8. Бутко Я.А., Смолянов О.Г. Формулы Фейнмана в квантовой и стохастической динамике // Современные проблемы математики и механики. 2011. Т. 6, № 1. С. 61-75.
9. Бутко Я.А., Смолянов О.Г., Шиллинг Р.Л. Формулы Фейнмана для феллеровских полугрупп // Доклады Академии Наук. 2010. T. 434, № 1. C. 7-11.
10. Смолянов О.Г., Толстыга Д.С., Вайцзеккер X. Фейнмановское описание одномерной динамики частиц с кусочно-непрерывной зависимостью массы от координаты // Доклады Академии Наук. 2011. T. 441, № 3. C. 295-298.
11. Волконский В.А. Аддитивные функционалы от марковских процессов // Труды Московского Математического Общества. 1960. T. 9. C. 143-189.
12. Обрезков О.О. Формула Фейнмана для задачи Коши — Дирихле в ограниченной области // Математические заметки. 2005. T. 77, №2. C. 316-320. DOI: 10.4213/mzm2493
13. Орлов Ю.Н., Сакбаев В.Ж., Смолянов О.Г. Скорость сходимости Фейнмановских аппроксимаций полугрупп, порождаемых гамильтонианом осциллятора // Теоретическая и математическая физика. 2012. T. 172, № 1. C. 122-137.
14. Портенко Н.И., Скороход А.В., Шуренков В.М. Марковские процессы // Итоги науки и техники. Серия «Современные проблемы математики». Фундаментальные направления. Т. 46. Теория вероятностей-4. М.: ВИНИТИ, 1989. 248 с.
15. Albeverio S., Brzezniak Z. Oscillatory integrals on Hilbert spaces and Schrodinger equation with magnetic fields // J. Math. Phys. 1995. Vol. 36, no. 5. P. 2135-2156.
16. Baur B., Conrad F., Grothaus M. Smooth contractive embeddings and application to Feynman formula for parabolic equations on smooth bounded domains // Communications in Statistics: Theory and Methods. 2011. Vol. 40, no. 19-20. P. 3452-3464.
17. Brezis H., Pazy A. Semigroups of Nonlinear Contractions on Convex Sets // J. Func. Anal. 1970. Vol.6. P. 237-281.
18. Bottcher B. On the construction of Feller processes with unbounded coefficients // Electronic Communications in Probability. 2011. Vol. 16. P. 545-555.
19. Bottcher B., Butko Ya.A., Schilling R.L., Smolyanov O.G. Feynman formulae and path integrals for some evolutionary semigroups related to tau-quantization // Rus. J. Math. Phys. 2011. Vol. 18, no. 4. P. 387-399.
20. 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. P. 2629-2638.
21. Butko Ya., Grothaus M., Smolyanov O.G. Lagrangian Feynman Formulae for Second Order Parabolic Equations in Bounded and Unbounded Domains // Inf. Dim. Anal. Quant. Probab. 2010. Vol. 13, no. 3. P. 377-392.
22. 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. P. 2009-2018.
23. Butko Ya.A., Schilling R.L., Smolyanov O.G. Lagrangian and Hamiltonian Feynman formulae for some Feller semigroups and their perturbations // Inf. Dim. Anal. Quant. Probab. Rel. Top. 2012. Vol. 15, no. 3. 26 p. DOI: 10.1142/S0219025712500154
24. Chernoff P. Note on product formulas for operator semigroups // J. Func. Anal. 1968. Vol. 2. P. 238-242.
25. Chernoff P. Product formulas, nonlinear semigroups and addition of unbounded operators // Mem. Am. Math. Soc. 1974. Vol. 140. 121 p.
26. Dorroh J.R. Contraction semi-groups in a function space // Pacific J.Math. 1966. Vol. 19, no. 1. P. 35-38.
27. Engel K.J., Nagel R. One-parameter semigroups for linear evolution equations. Springer, 1995. 586 p.
28. Ethier S.E., Kurtz T.G. Markov Processes: Characterization and Convergence. New York: Wiley, 1986. 534 p. (Wiley Ser. Probab. Math. Stat.).
29. Feynman R.P. Space-time approach to nonrelativistic quantum mechanics // Rev. Mod. Phys. 1948. Vol.20. P. 367-387.
30. Feynman R.P. An Operator Calculus Having Applications in Quantum Electrodynamics // Phys. Rev. 1951. Vol. 84. P. 108-128.
31. Freidlin M. Functional integration and partial differential equations. Princeton, New Jersey: Princeton University Press, 1985. 545 p.
32. Friedman A. Partial differential equations of parabolic type. Prentice-Hall, 1964. 427 c.
33. Gustafson K., Lumer G. Multiplicative perturbation of semigroup generators // Pacific J. Math. 1972. Vol.41, no. 3. P. 731-742.
34. Jacob N. Pseudo-differential operators and Markov processes. Vol. I-II. Imperial College Press, 2001. 946 p.
35. Karatzas I., Shreve S. Brownian Motion and Stochastic Calculus. Springer, 1991. 470 p.
36. Kühnemund F. Bi-continuous semigroups on spaces with two topologies: theory and applications. Dissertation der Mathematischen Fakultat der Eberhard-Karls-Universitat Tübingen zur Erlangung des Grades eines Doktors der Naturwissenschaften, 2001. 104 p.
37. Lunardi A. Analytic semigroups and optimal regularity in parabolic problems. Birkhauser, 1995. 424 p.
38. Lumer G. Perturbation de generateurs infinitesimaux du type "changement de temps" // Ann. Inst. Fourier. 1974. Vol. 23, no. 4. P. 271-279.
39. Obrezkov O., Smolyanov O.G., Truman A. The Generalized Chernoff Theorem and Randomized Feynman Formula // Doklady Math. 2005. Vol. 71, no. 1. P. 105-110.
40. Pazy A. Semigroups of Linear Operators and Applications to Partial Differential Equations. New York: Springer-Verlag, 1983. 276 p.
41. Plyashechnik A.S. Feynman formulas for second-order parabolic equations with variable coefficients // Rus. J. Math. Phys. 2013. Vol. 20, no. 3. P. 377-379.
42. Plyashechnik A.S. Feynman Formula for Schrodinger-type equations with time- and space-dependent coefficients //Rus. J. Math. Phys. 2012. Vol. 19, no. 3. P. 340-359.
43. Reed M., Simon B. Methods of Modern Mathematical Physics. Vol. II. Academic Press, 1980. 361 p.
44. Sakbaev V.G., Smolyanov O.G. Dynamics of a Quantum Particle with Discontinuous Position-Dependent Mass // Doklady Math. 2010. Vol. 82, no. 1. P. 630-634.
45. Sato K. Levy Processes and Infinitely Divisible Distributions. Cambridge Univ. Press, 1999. 486 p.
46. Smolyanov O.G. Feynman type formulae for quantum evolution and diffusion on manifolds and graphs // Quant. Bio-Informatics, World Sc. 2010. Vol. 3. P. 337-347.
47. Smolyanov O.G., Shamarov N.N. Feynman and Feynman-Kac formulae for evolution equations with Vladimirov operator // Doklady Math. 2008. Vol. 77, no. 3. P. 345-349.
48. Smolyanov O.G., Tokarev A.G., Truman A. Hamiltonian Feynman path integrals via the Chernoff formula//J. Math. Phys. 2002. Vol.43, no. 10. P. 5161-5171.
49. Smolyanov O.G., Weizsacker H., 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. P. 349-358.
50. Smolyanov O.G., Weizsacker H., Wittich O. Chernoff's Theorem and Discrete Time Approximations of Brownian Motion on Manifolds // Potent. Anal. 2007. Vol. 26, no. 1. P. 1-29.
51. Weizsaecker H., Winkler G. Stochastic Integrals: an Introduction. Vieweg, 1990. 332 p.
52. Zhang G., Jiang M. Parabolic equations and Feynman-Kac formula on general bounded domains // Sci. in China. 2001. Vol. 44, no. 3. P. 311-329.
scientific periodical of the bauman mstu
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Feynman formulae for evolution semigroups
# 03, March 2014
DOI: 10.7463/0314.0701581
Butko Ya. A.
Bauman Moscow State Technical University 105005, Moscow, Russian Federation [email protected]
The paper systematically describes an approach to solution of initial and initial-boundary value problems for evolution equations based on the representation of the corresponding evolution semigroups with the help of Feynman formulae. The article discusses some of the methods of constructing Feynman formulae for different evolution semigroups, presents specific examples of solutions of evolution equations. In particular, Feynman formula is obtained for evolution semigroups generated by multiplicative perturbations of generators of some initial semigroups. In this case semigroups on a Banach space of continuous functions defined on an arbitrary metric space are considered; Feynman formulae are constructed with the help of operator families, which are Chernoff equivalent to the initial unperturbed semigroups. The present result generalizes the author's paper "Feynman formula for semigroups with multiplicative perturbed generators" and some of the results of the joint with O.G. Smolyanov and R.L. Schilling paper "Lagrangian and Hamiltonian Feynman formulae for some Feller processes and their perturbations". The approach to the construction of Feynman formulae for semigroups with multiplicative and additive perturbed generators is illustrated with examples of the Cauchy problem for the Schrodinger equation, the approximation of transition probabilities of some Markov processes.
Further, a wider class of additive and multiplicative perturbations of a particular generator — the Laplace operator — is considered in the paper. And Feynman formula for the solution of the Cauchy problem for a second order parabolic equation with unbounded variable coefficients is proved. In addition, the article describes a method for constructing Feynman formulae for solutions of the Cauchy — Dirichlet problem for parabolic differential equations. The method is also illustrated by a second order parabolic equation with variable coefficients. These results generalize some of the results of the work by Butko, Grothaus and Smolyanov "Lagrangian Feynman formulae for Second Order Parabolic Equations in Bounded and Unbounded Domains". The article also discusses some of the Feynman — Kac formulae and Feynman integrals related to the obtained Feynman formulae.
Publications with keywords: Feynman formula, multiplicative perturbations, equation of evolution, evolution semigroups, Feynman — Kac formulae, additive perturbation, approximation of transition probabilities, Schrodinger equation
Publications with words: Feynman formula, multiplicative perturbations, equation of evolution, evolution semigroups, Feynman — Kac formulae, additive perturbation, approximation of transition probabilities, Schrodinger equation
References
1. Bogachev V.I., Smolyanov O.G. Deystvitel'nyyifunktsional'nyyanaliz: universitetskiy kurs [Real and functional analysis: university course]. Moscow-Izhevsk, Publ. of NITs "Regul-yarnaya i khaoticheskaya dinamika", 2011. 728 p. (in Russian).
2. Butko Ya.A. [Representations of evolution semigroups by Feynman formulae and phase space Feynman path integrals]. Nauka i obrazovanie MGTU im. N.E. Baumana — Science and Education of the Bauman MSTU, 2012, no. 2. Available at: http://technomag.bmstu.ru/ doc/315838.html, accessed 01.02.2014. (in Russian).
3. Butko Ya.A. [Feynman Formulas and Functional Integrals for Diffusion with Drift in a Domain on a Manifold]. Matematicheskie zametki, 2008, vol.83, no. 3, pp. 333-349. DOI: 10.4213/mzm3772 (English translation: Mathematical Notes, 2008, vol. 83, no. 3-4, pp. 301316. DOI: 10.1134/S0001434608030024).
4. Butko Ya.A. [Feynman formula for semigroups with multiplicatively perturbed Generators]. Nauka i obrazovanie MGTU im. N.E. Baumana—Science and Education of the Bauman MSTU, 2011, no. 10. Available at: http://technomag.bmstu.ru/doc/239563.html, accessed 01.02.2014. (in Russian).
5. Butko Ya.A., Grotkhaus M., Smolyanov O.G. [Feynman formula for a class of second order parabolic equations in a bounded domain]. Doklady Akademii Nauk, 2008, vol.421, no. 6, pp. 727-732. (English translation: Doklady Mathematics, 2008, vol.78, no. 1, pp. 590-595. DOI: 10.1134/S1064562408040327).
6. Butko Ya.A., Duryagin A.V. [Feynman formulae for a family of parabolic equations corresponding to the tau-quantization of a quadratic Hamilton function]. Nauka i obrazovanie MGTU im. N.E. Baumana — Science and Education of the Bauman MSTU, 2011, no. 11. Available at: http://technomag.bmstu.ru/doc/251251.html, accessed 01.02.2014. (in Russian).
7. Butko Ya.A., Morozov A.V. [Representation of the solution of the Cauchy-Neumann problem for a parabolic equation on a ray by a Lagrangian Feynman formula]. Nauka i obrazovanie MGTU im. N.E. Baumana — Science and Education of the Bauman MSTU, 2011, no. 11. Available at: http://technomag.bmstu.ru/doc/246219.html, accessed 01.02.2014. (in Russian).
8. Butko Ya.A., Smolyanov O.G. [Formuly Feynmana v stokhasticheskoy i kvantovoy dinamike]. Sovremennyeproblemy matematiki i mekhaniki, 2011, vol. 6, no. 1, pp. 61-75. (in Russian).
9. Butko Ya.A., Smolyanov O.G., Shilling R.L. [Feynman formulae for Feller semigroups]. DokladyAkademii Nauk, 2010, vol. 434, no. 1, pp. 7-11. (English translation: DokladyMath-ematics, 2010, vol. 82, no. 2, pp. 679-683. DOI: 10.1134/S1064562410050017).
10. Smolyanov O.G., Tolstyga D.S., Weizsacker H. [Feynman description of the one-dimensional dynamics of particles whose masses piecewise continuously depend on coordinates]. Doklady Akademii Nauk, 2011, vol. 441, no. 3, pp. 295-298. (English translation: Doklady Mathematics, 2011, vol. 84, no. 3, pp. 804-807. DOI: 10.1134/S1064562411070209).
11. Volkonskiy V.A. [Additive functionals of Markov processes]. Trudy Moskovskogo Matem-aticheskogo Obshchestva, 1960, vol. 9, pp. 143-189. (in Russian).
12. Obrezkov O.O. [Feynman's formula for the Cauchy-Dirichlet problem in bounded domains]. Matematicheskie zametki, 2005, vol. 77, no. 2. pp. 316-320. DOI: 10.4213/mzm2493 (English translation: Mathematical Notes, 2005, vol.77, no. 1-2, pp. 288-293. DOI: 10.1007/s11006-005-0029-8).
13. Orlov Yu.N., Sakbaev V.Zh., Smolyanov O.G. [Rate of convergence of Feynman approximations of semigroups generated by theoscillator Hamiltonian]. Teoreticheskaya i matematich-eskaya fizika, 2012, vol. 172, no. 1, pp. 122-137. (English translation: Theoretical andMath-ematicalPhysics, 2012, vol. 172, no. 1, pp. 987-1000. DOI: 10.1007/s11232-012-0090-x).
14. Portenko N.I., Skorokhod A.V., Shurenkov V.M. [Markov processes]. Itogi nauki i tekhniki. Seriya Sovremennye problemy matematiki. Fundamental'nye napravleniya. T. 46. Teoriya veroyatnostey-4 [Results of science and technology. Series Modern Problems of Mathematics. Fundamental Directions. Vol. 46. Probability Theory-4]. Moscow, VINITIPubl., 1989. 248 p. (in Russian).
15. Albeverio S., Brzezniak Z. Oscillatory integrals on Hilbert spaces and Schrodinger equation with magnetic fields. J.Math. Phys., 1995, vol. 36, no. 5, pp. 2135-2156.
16. Baur B., Conrad F., Grothaus M. Smooth contractive embeddings and application to Feynman formula for parabolic equations on smooth bounded domains. Communications in Statistics: Theory andMethods, 2011, vol. 40, no. 19-20, pp. 3452-3464.
17. Brezis H., Pazy A. Semigroups of Nonlinear Contractions on Convex Sets. J. Func. Anal., 1970, vol.6, pp. 237-281.
18. Böttcher B. On the construction of Feller processes with unbounded coefficients. Electronic Communications in Probability, 2011, vol. 16, pp. 545-555.
19. Bottcher B., Butko Ya.A., Schilling R.L., Smolyanov O.G. Feynman formulae and path integrals for some evolutionary semigroups related to tau-quantization. Rus. J.Math. Phys., 2011, vol. 18, no. 4, pp. 387-399.
20. 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.
21. Butko Ya., Grothaus M., Smolyanov O.G. Lagrangian Feynman Formulae for Second Order Parabolic Equations in Bounded and Unbounded Domains. Inf. Dim. Anal. Quant. Probab., 2010, vol. 13, no. 3, pp. 377-392.
22. 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.
23. Butko Ya.A., Schilling R.L., Smolyanov O.G. Lagrangian and Hamiltonian Feynman formulae for some Feller semigroups and their perturbations. Inf. Dim. Anal. Quant. Probab. Rel. Top., 2012, vol. 15, no. 3. 26 p. DOI: 10.1142/S0219025712500154
24. Chernoff P. Note on product formulas for operator semigroups. J. Func. Anal., 1968, vol.2, pp. 238-242.
25. Chernoff P. Product formulas, nonlinear semigroups and addition of unbounded operators. Mem. Am. Math. Soc., no. 140. American Mathematical Society, Providence, R. I., 1974. 121 p.
26. Dorroh J.R. Contraction semi-groups in a function space. Pacific J. Math., 1966, vol. 19, no. 1, pp. 35-38.
27. Engel K.J., Nagel R. One-parameter semigroups for linear evolution equations. Springer, 1995. 586 p.
28. Ethier S.E., Kurtz T.G. Markov Processes: Characterization and Convergence. New York, Wiley, 1986. 534 p. (Wiley Ser. Probab. Math. Stat.).
29. Feynman R.P. Space-time approach to nonrelativistic quantum mechanics. Rev. Mod. Phys., 1948, vol.20, pp. 367-387.
30. Feynman R.P. An Operator Calculus Having Applications in Quantum Electrodynamics. Phys. Rev., 1951, vol. 84, pp. 108-128.
31. Freidlin M. Functional integration and partial differential equations. Princeton, New Jersey, Princeton University Press, 1985. 545 p.
32. Friedman A. Partial differential equations of parabolic type. Prentice-Hall, 1964. 427 p.
33. Gustafson K., Lumer G. Multiplicative perturbation of semigroup generators. Pacific J. Math., 1972, vol. 41, no. 3, pp. 731-742.
34. Jacob N. Pseudo-differential operators and Markov processes. Vol. 1-2. Imperial College Press, 2001. 946 p.
35. Karatzas I., Shreve S. Brownian Motion and Stochastic Calculus. Springer, 1991. 470 p.
36. Kühnemund F. Bi-continuous semigroups on spaces with two topologies: theory and applications. Dissertation der Mathematischen Fakultat der Eberhard-Karls-Universitat Tübingen zur Erlangung des Grades eines Doktors der Naturwissenschaften, 2001. 104 p.
37. Lunardi A. Analytic semigroups and optimal regularity in parabolic problems. Birkhauser, 1995. 424 p.
38. Lumer G. Perturbation de generateurs infinitesimaux du type "changement de temps". Ann. Inst. Fourier, 1974, vol. 23, no. 4, pp. 271-279.
39. 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.
40. Pazy A. Semigroups of Linear Operators and Applications to Partial Differential Equations. New York, Springer-Verlag, 1983. 276 p.
41. Plyashechnik A.S. Feynman formulas for second-order parabolic equations with variable coefficients. Rus. J.Math. Phys., 2013, vol. 20, no. 3, pp. 377-379.
42. Plyashechnik A.S. Feynman Formula for Schrodinger-type equations with time- and space-dependent coefficients. Rus. J.Math. Phys., 2012, vol. 19, no. 3, pp. 340-359.
43. ReedM., Simon B. Methods of Modern Mathematical Physics. Vol. 2. Academic Press, 1980. 361 p.
44. Sakbaev V.G., Smolyanov O.G. Dynamics of a Quantum Particle with Discontinuous Position-Dependent Mass. Doklady Math., 2010, vol. 82, no. 1, pp. 630-634.
45. Sato K. Levy Processes and Infinitely Divisible Distributions. Cambridge Univ. Press, 1999. 486 p.
46. 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.
47. 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.
48. 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.
49. Smolyanov O.G., Weizsacker H., Wittich O. Chernoff's theorem and the construction of semigroups. In: Evolution Equations: Applications to Physics, Industry, Life sciences and Economics. Birkhauser. Prog. Nonlinear Differ. Eq. Appl., 2003, vol. 55, pp. 349-358.
50. Smolyanov O.G., Weizsacker H., Wittich O. Chernoff's Theorem and Discrete Time Approximations of Brownian Motion on Manifolds. Potent. Anal., 2007, vol. 26, no. 1, pp. 1-29.
51. Weizsaecker H., Winkler G. Stochastic Integrals: an Introduction. Vieweg, 1990. 332 p.
52. Zhang G., Jiang M. Parabolic equations and Feynman-Kac formula on general bounded domains. Sci. in China, 2001, vol. 44, no. 3, pp. 311-329.