электронное научно-техническое издание
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 30569. Государственная регистрация №0421100025. ISSN 1994-0408
Формулы Фейнмана для семейства параболических уравнений, соответствующих т-квантованию квадратичной функции Гамильтона
77-30569/251251
# 10, октябрь 2011 Я.А. Бутко, A.B. Дурягин
УДК 517.987.4
МГТУ им. Н.Э. Баумана [email protected], [email protected]
1. Введение
Формула Фейнмана это представление полугруппы, соответствующей некоторой начально-краевой задаче для эволюционного уравнения, в виде предела интегралов от элементарных функций по декартовым степеням некоторого пространства Е при стремлении степени к бесконечности. Конечнократные интегралы в формулах Фейнмана аппроксимируют функциональные интегралы по некоторым вероятностным мерам или фейнмановским псевдомерам на множестве функций, определенных на интервале [0,£] и принимающих значения в Е. Таким образом, формулы Фейнмана могут быть использованы как для непосредственного вычисления решений начально-краевых задач, так и для представления решений таких задач с помощью функциональных интегралов (такие представления называются формулами Фейнмана-Каца), аппроксимации переходных вероятностей некоторых случайных процессов и переходных амплитуд в квантовой динамике, компьютерного симулирования случайных процессов и моделирования квантовой и стохастической динамики. Формула Фейнмана называется гамилътоновой формулой Фейнмана, если пространство Е, декартовы степени которого рассматриваются, является фазовым пространством некоторой классической гамильтоновой системы. Такие формулы Фей-нмана дают аппроксимации функциональным интегралам по траекториям в
этом самом фазовом пространстве. Аналогично определяются лагранжевы формулы Фейнмана с помощью замены фазового пространства конфигурационным. Лагранжевы формулы Фейнмана и, соответственно, лагранжевы формулы Фейнмана-Каца появились впервые в работе Фейнмана [10] 1948 года, в которой рассматривалось простейшее уравнение Шредингера с потенциалом. Первое строгое математическое доказательство результатов работы [10] было дано Нельсоном в 1964 г. и основывалось на формуле Троттера. Соответствующие гамильтоновы формулы Фейнмана и Фейнмана-Каца были предложены в работе Фейнмана [11], опубликованной в 1951 году. Однако первое строгое доказательство таких формул появилось только в 2002 г. в работе [23]. Возможно, такая длительная задержка связана с тем, что формулы Троттера не достаточно для доказательства гамильтоновых формул Фейнмана и необходимо вместо нее использовать теорему Чернова [9]. Теорема Чернова является естественным обобщением формулы Троттера и может быть использована для доказательства формул Фейнмана и Фейнмана-Каца для широкого класса начальных и начально-краевых задач для разнообразных эволюционных уравнений. Основанный на теореме Чернова метод получения формул Фейнмана и Фейнмана-Каца был разработан в серии работ [23], [24]-[28] и применен для для описания классической, квантовой и стохастической динамики в областях евклидовых пространств и римановых многообразий, в бесконечномерных линейных и нелинейных пространствах, р-адических пространствах и т.д. (см., например, [1] - [8], [12], [14], [17], [18], [19], [20], [23], [24], [26], [28]). В настоящей работе данный метод используется для представлений решения задачи Коши-Дирихле для семейства параболических уравнений второго порядка с помощью гамильтоновых формул Фейнмана. Рассмотренное в работе семейство параболических уравнений отвечает различным видам квантования квадратичной функции Гамильтона классической гамильтоновой системы, в том числе др-, рд-квантованиям и квантованию Вейля.
2. Теорема Чернова и формула Фейнмана
Пусть (X, У • ) — банахово пространство, С(Х) — пространство ограниченных линейных операторов на X, наделенное сильной операторной топологией, || • || — операторная норма на С(Х), Н — тождественный оператор
в X. Если Dom (L) С X линейное подпространство и L : Dom(L) ^ X линейный оператор, то символом Dom (L) обозначим область определения L. Однопараметрическое семейство (Tt)t>о ограниченных линейных операторов Tt : X ^ X называется сильно непрерывной полугруппой, если T0 = Id, Ts+t = Ts о Tt для всех s,t > 0 и limt^0 \\^ф — ф\\х = 0 для всех ф £ X. Пусть (Tt)t>0 — сильно непрерывная полугруппа на банаховом пространстве (X, \\ • \\х). Тогда генератор L полугруппы (Tt)t>0 определяется как Lф := lim Т1ф—ф с областью определения
Dom(L) := {ф £ X| lim Ф X.
I ' t^0 t J
Рассмотрим эволюционное уравнение dt = Lf. Если L — генератор сильно непрерывной полугруппы (Tt)t>0 на банаховом пространстве (X, \\ • \\х), то решение задачи Коши для этого уравнения с начальным условием f (0) = f0 £ X задается как f (t) = Ttf0 для всех f0 £ X. Следовательно, решение задачи Коши для эволюционного уравнения ^ = Lf равносильно построению полугруппы (Tt)t>0 с заданным генератором L. Как правило, действие полугруппы не удается получить в явном, аналитическом виде. Однако, полугруппу с заданным генератором возможно аппроксимировать. Одно из средств аппроксимации полугрупп основывается на теореме Чернова [9] (далее представлена версия теоремы Чернова, введенная в [23]).
Теорема 1 (Чернова). Пусть X — банахово пространство, F : [0, ж) ^ L(q) — сильно непрерывное отображение, такое что F(0) = Id и \\F(t) \\ < eAt для некоторого Л £ [0, ж) и всех t > 0. Пусть D линейное подпространство Dom(F '(0)) такое, что сужение оператора F '(0) на это подпространство замыкаемо, (L, Dom(L)) —соответствующее замыкание. Пусть (L, Dom(L)) является генератором сильно непрерывной полугруппы (Tt)t>0. Тогда для любого t0 > 0 последовательность (F(t/n))n)n£N сходится к (Tt)t>0 при n ^ ж в сильной операторной топологии, равномерно по t £ [0,t0], т.е. Tt = limn^^F(t/n))n. Здесь производная функции F : [0, e) ^ L(q), е > 0, это линейное отображение F'(0) : Dom (F'(0)) ^ X такое, что
F'(0)g := lim F(t)g — F(0)g,
t^0 t
где Dom (F'(0)) — векторное пространство элементов д G X, для которых этот предел существует.
Семейство операторов (F(t))t>о называется эквивалентным по Чернову полугруппе (Tt)t>0, если семейство удовлетворяет всем требованиям теоремы Чернова по отношению к полугруппе (Tt)t>0; т.е., если по теореме Чернова выполняется равенство
Tt = lim (F(t/n))n. (1)
В большинстве случаев операторы F (t) — это интегральные операторы и тогда в правой части равенства (1) стоит предел последовательности n-кратных интегралов. В этом случае равенство (1) называется формулой Фейнмана.
Определение 1. Формула Фейнмана — это представление решения начальной (или начально-краевой) задачи для эволюционного уравнения (или, что то же, представление соответствующей эволюционной полугруппы) в виде предела n-кратных интегралов при n ^ то.
Мы будем называть равенство Tt = limn^TO(F(t/n))n лагранжевой формулой Фейнмана, если F (t), t > 0, это интегральные операторы, ядра которых выражаются через элементарные функции; если F (t) — псевдо-диффе-ренциальные операторы (см. определение 3), то назовем это равенство га-мильтоновой формулой Фейнмана. Такая терминология связана с тем, что лагранжевы формулы Фейнмана аппроксимируют функциональне интегралы по множеству траекторий в конфигурационном пространстве системы (эволюция которой описывается полугруппой (Tt)t>0), в то время как гамильтоновы формулы Фейнмана дают аппроксимации функциональных интегралов по множеству траекторий в фазовом пространстве некоторой системы.
3. Псевдо-дифференциальные операторы, их символы и т-квантование
Пусть H : Rd х Rd ^ C измеримая функция и т G [0,1]. Мы определим псевдо- дифференциальный оператор HT( ,D) с т-символом H(q,p) на некотором банаховом пространстве (X, || • ||х) функций на Rd как
Ht(-,D)0(q) = (2n)-d^ J ^H(Tq + (1 - T)qbp)0(qi)dqidp, (2)
где область определения Оот(Нт(•, Л)) — это множество всех ф £ X таких, что правая часть формулы (2) корректно определена как элемент пространства (X, У • ||х). Далее, мы будем рассматривать псевдо-дифференциальные операторы на пространстве := {ф £ С(К^) : Нш|^0 ф(д) =0} с нормой ||ф||° := |ф(д)|. Мы всегда будем считать что множество пробных функций С°°(^) принадлежит области определения оператора Нт(•, Л).
Отображение • т : Н ^ Нт(•, Л) из пространства функций на ^ х ^ в пространство линейных операторов на (X, || • ||х) называется т-квантованием, оператор Нт(•, Л) называется т-квантованием функции Н. Заметим, что если Н — сумма функций, зависящих только от одной переменной д или р, то операторы Нт(-,Л) совпадают для всех т £ [0,1]. Если Н(д,р) = др = рд, д,р £ К1, то
И д д Нт(•, Л)ф(д) = -^тдддф(д) - ¿(1 - т)дд(дф(д)).
Таким образом, различные т соответствуют различному порядку применения некоммутирующих операторов. При этом, '^"-квантование соответствует случаю т = 1, '^"-квантование — случаю т = 0 и вейлевское квантование — случаю т = 1/2. Функция Н(д,р) обычно рассматривается как функция Гамильтона некоторой классической системы. Оператор Нт(•, Л) называется гамильтонианом квантовой системы, полученной т-квантованием классической системы с функцией Гамильтона Н.
4. Связь между т-квантованиями квадратичной функции Гамильтона
В дальнейшем р, д £ К1, С^ (К) — пространство к раз непрерывно дифференцируемых функций, ограниченных вместе со всеми своими производными. Мы будем рассматривать квадратичную функцию Гамильтона
Н (р, д) = а(д)(^р)2 + 6(д)гр + с(д),
где а( ) — непрерывная на К строго положительная ограниченная функция, причем а(д) не стремится к нулю на бесконечности, Ь(д) и с(д) — непрерывные и ограниченные на К функции.
Лемма 1. Пусть а(-) G Cb2(R), &(•) G C^R), c(^) G Cb(R). Тогда при каждом т G [0,1], оператор HT (,D) продолжается на множество C2 (R) по формуле
д2 д
H(q, D)0(q) = a(q) фМ + 2(1 - тK(q)^(q) +
д
+ (1 - т) V(q^(q) + b(q)^(q) + (1 - т)b/(q)0 + c(q)0(q). Пусть выполнены условия леммы. Тогда для ф G C(?0(R) имеем:
+то + то
Ht(q, D^(q) = 2п / / н(тq + (1 - тqiV(ql)dqldp =
—то -то +то +то
= 2П / / a^ + (1 - т)qi)(ip)2eip(<?-<?1V(qi)dqidp +
—то -то +то +то
+ 2л/ / ^ +(1 - т)qi)ipeip(q-qiV(ql)dqldp +
-то -то
+то +то
+ 2п// C(тq + (1 - тЫ^'-qiV(ql)dqldp.
-то -то
Пусть сначала a, b, c G Cто(К). Рассмотрим каждое слагаемое по отдельности:
+то +то
J J a(rq + (1 - т)ql)(ip)2eip(q-qi^(ql)dqldp =
-то -то
r
= — lim a^q + (1 - т)ql)(ip)2eip(q-'91 ^(q^dq^p =
2п Г^то J J -r sup ф
r
= — lim a^q + (1 - т^Оф^) (ip^e^'-qi)dpdqi_ = 2n r уто J J
supp ф -r
r
1 С d'2 С
= — lim a(rq + (1 - т/ eip(q-qi)dpdql =
2п г^то J dql2 J
supp ф - r
d2
5(q - ql), a(тq + (1 - тЫфЫ) =
= (6 (д- д^ а(тд +(1 - т МЖдО) =
= (1 - т)2о"(д)ф(д) + 2(1 - т)а,(д)ф,(д) + а(д)ф"(д). (3)
Здесь (/, д) обозначает действие обобщенной функции / £ (К) на основную функцию д £ Б (К), 6 — дельта-функция Дирака.
Аналогичным образом, для второго и третьего слагаемых получим:
+ТО +ТО
^ [ [ Ь(тд + (1 - т^¿ре^^ф^^д^р =
—то -то
г
11ш [ [ (тд + (1 - т^(¿рУ^-^фЫ^р = 2п г ^то у у
-г вирф
г
= 11ш [ Ь(тд + (1 - тЫф(дО / ¿регр(м)^д1 = 2п г ^то у у
вирр ф -г
г
= -— 11ш [ Ь(тд + (1 - т)д1 )ф(д1 / егр("^д1
2п г^то ] йд1 ]
вирр ф -г
_6(д - д^ ь(тд +(1 - т МЖдО) =
= -( 6(д - д^ +(1 - т МЖдО) =
= (1 - т )Ь,(д )ф(д) + Ь(д)ф,(д);
+ТО +ТО
71 / / с(тд + (1 - т)д1)егр(9-91)ф(д1)^д1^р
-° -°
г
= 11ш [ [ с(тд + (1 - т)д1)егр(м)ф(д1)^д1ф = 2п г ^то у у
-г вирф
г
= 11ш [ с(тд + (1 - тЫф(дО / е^^^р^ = 2п г^то у у
вирр ф -г
= (6(д - д1)с(тд +(1 - т Ыф(дО) = с(д)ф(д).
Объединяя слагаемые, получим:
д2 д
Нт(д, Я)ф(д) = а(д)^ф(д) + 2(1 - т)а%)д^Ф(^) +
д
+ (1 - т)2а//(д)ф(д) + Ь(д)-ф(д) + (1 - т)Ь/(д)ф + ф)ф(д).
Если теперь а(^) Е С(К), &(•) Е С (К), с(^) Е Сь(К), то существуют последовательности бесконечно дифференцируемых функций {ап(0}пЕм, {&п(0}пЕм, {сп(0}пЕм, локально равномерно на К сходящихся к а(-), &(•), с(^) вместе со всеми своими производными до второго, первого и нулевого порядка соответственно.
Тогда обобщенные функции
б2 б а(тд + (1 - т^^% - дО, Ь(тд + (1 - т- яО
и
с(тд + (1 - т)д1 )% - дО
являются слабыми пределами обобщенных функций
б2 б ап(тд + (1 - т)д1)^^% - яО, &п(тд + (1 - т- яО
и
сп(тд + (1 - т)д0% - дО соответственно, а для последних справедливы выкладки, приведенные выше. ►
Рассмотрим вид оператора Нт при некоторых т. При т = 1, т. е. при др-квантовании:
д2 д д? Ф(д) + Ь(д) ^ При т = 0, т.е. при рд-квантовании:
д2 д Яо(д,£)ф(д) = ^^^ [а(д)ф(д)] + ^ [%Жя)] + с(д)ф(д).
При т = 1, т.е. при квантовании Вейля:
д2 д 1 Нг (д, Я)ф(д) = ^Ф(я) + ^Ф(я) + 4«"(яЖя) +
д1
+ Ь(д) —Ф(д) + 2 Ь/(д)Ф + с(д)ф(д).
5. Задача Коши-Дирихле для семейства параболических уравнений второго порядка на отрезке
Пусть G = (a,ß) С R, G = [a,ß], т G [0,1]. Рассмотрим начально-краевую задачу Коши-Дирихле для параболического уравнения второго порядка в области G:
f (t,q) = HT(q,D)f(t,q), t > 0, q G G, f (0, q) = fo(q), q G G, (4)
f (t, q) = 0, t > 0, q G dG.
Здесь HT (•, D) — оператор, полученный т-квантованием функции
H(p, q) = a(q)(ip)2 + b(q)ip + c(q),
где a( ) — непрерывная на R строго положительная ограниченная функция, причем a(q) не стремится к нулю на бесконечности, b( ) и c( ) — непрерывные и ограниченные на R функции. Пусть Да — дифференциальный оператор на множестве дважды дифференцируемых функций на R, определяемый равенством: (Даф)^) = a(q)0(q). Будем также использовать обозначение V0(q) ДОя dq0(q) .
Мы предполагаем, что f0 G X и f (t, •) G X для всех t > 0, где X = C0 (G) — банахово пространство непрерывных на G функций, обращающихся в ноль на границе. Норма в пространстве X = C0(G) определена стандартным образом: ||f ||X = SUPqGG |f (q)|.
Пусть для каждого т G [0,1], семейство операторов (TtT )t>0 — это полугруппа на C0(G), разрешающая задачу Коши-Дирихле (4). Мы предполагаем, что коэффициенты рассматриваемого уравнения удовлетворяют достаточным условиям для того, чтобы полугруппа (TtT)t>0 была сильно непрерывной полугруппой на пространстве C0(G). Обозначим через (LT, Dom(LT)) генератор полугруппы (TtT)t>0. Пусть C0°(G) := CTO(G) П C0(G). Мы предполагаем, что множество D := {ф G C0°(G)| Ьф G C0(G)} — существенная область определения (LT, Dom (LT)).
Семейства операторов, эквивалентные по Чернову полугруппе, разрешающей задачу Коши-Дирихе при т = 1
Рассмотрим задачу Коши-Дирихле (4) при т = 1. Пространство C0(G) можно рассматривать как подпространство Cc(R) непрерывных в R функций с компактным носителем. Определим для любой ф G C0(G) непрерывное продолжение до функции ф G Cc(R), положив:
1.Для любых Ф1, Ф2, фз G Co(G) и Л, ß G R выполняется импликация: если Ф1 = Лф2 + дф3 на G, то <ф1 = Лф2 + дф3 на R.
2. 1 |ф||Со(G) = SUPqGM |ф(Я)|.
3. Если ф G C0°°(G), то ф G Cc°°(R).
Векторное пространство всех продолжений ф функций ф G C°(G) назовем D с Cc(R).
Определим на S(R) семейство операторов
ptc(q) г г ( )t 2
[F (¿)ф](д) = — е- 2 eitp6(q)eip(q-qoV(qo)dqodp.
-° -°
Лемма 2. Для любого ф из S(R) функция F(¿)ф принадлежит C°(R), где C°(R) — банахово пространство непрерывных в R функций таких, что
f (q) ^ 0 при |q| ^
М При t = 0 оператор принимает вид
(F(t^)(q) = / eipq[Ffl(p)dp = ф(д).
-°
Следовательно, F(0) = Id на S(R). Здесь и далее ^ф обозначает преобразование Фурье функции ф. При t > 0
(F(t^)(q) = etc(q)-L J е-^eip(q+t6(q))[F^](p)dp.
-°
Покажем, что функция F(^ф непрерывна и убывает к нулю на бесконечности.
2
lim q2(F(t^)(q) = lim etc(q)—= q2e-^eip(q+t6(q))[F^](p)dp =
— 00
+ ТО
= 11ш е*с(?)
\/2П
дУм е-^ еград[^ф](р)«р =
— 00
1
+ ТО
= 11ш . 7 „
Ы\° л/2^ У «р2
е-^ е^^ ф](р)«р.
Интегрированием по частям далее получим
Нш -е^/
Ы\° л/2п У «р2
— 00
е-^ [^ ф](р)
«р,
где
«р2
е-^ е*^ ф](р)
2
е 2 е
гр£Ь(д)
[^ф],,(р) - 2[^ф],(р)(р^а(д) - ¿*6(д)) +
+ [^ф](р)[р2^а(д) + ¿¿2р2Ь(д)а(д) - ¿а(д) - ¿¿2Ь(д)а(д) - ¿2Ь2(д)]
22
<
<
[Т-ф]"(р) - 2[^ф],(р)(р^а(д) - ¿*6(д)) +
+ [^ф](р)[р2£а(д) + ¿¿2р2Ь(д)а(д) - ¿а(д) - ¿¿2Ь(д)а(д) - £2Ь2(д)]
22
Таким образом,
11ш д2(^(*)ф)(д)
<
< евирдек с(зО
\/2П
[^ф]» - 2[^ф],(р)(р^а(д) - ¿*6(д)) +
+ [^ф](р)[р2^а(д) + ¿¿2р2Ь(д)а(д) - ¿а(д) - ¿¿2Ь(д)а(д) - ¿2Ь2(д)]
22
«р.
В виду того, что функции а(д), Ь(д),с(д) ограничены в К, ф, ^ф £ Б (К), последний интеграл существует и равен константе. Его значение обозначим С. В итоге
Нш д2(^(*)ф)(д)
< евир«ек с(9)С = const.
1
1
Непрерывность по д функции ^(£)ф очевидна. Таким образом, доказано что
(^(г)ф)(д) = е*с(?)I е-^е^+ад)[Т-ф](р)«р
л/2п 3
отображает Б (К) в С°(К). ►
Лемма 3. На множестве Б (К) сильная производная ^ (£) в нуле совпадает с Н71, то есть для любого ф £ Б (К) равномерно по д £ К
1
НшГ1 (^(г)ф(д) - ф(д)) = 2^аф(д) + Ь(д)Уф(д) + с(д)ф(д).
1 (^ (^)ф - ф)(д) =
= Нш -
¿\о г
1
е-^ Х^+ВД)!
\/2П
е""4^ I е 2 е*Р(Ч 1 "Ь(Ч))[^ф](р)«р-
\/2П
[^ ф](р)«р
= Нш 1
¿\о г
= Нш 1
¿\о г
л/2п
,-^^ +гр*Ь(д)+*с(д) _
1
-°
_\/2П
е*(-^ +фЬ(<7)+с(д))
1
[^ ф](р)«р
\Т ф](р)«р
Нш -
\/2П 3 *\о г
¿(-^+гр6(д)+с(д)) _
1
[^ ф](р)«р.
— 00
Перестановка предельных переходов имеет место в силу теоремы Лебега о мажорированной сходимости. Далее получаем
1
а(д )р2
Т"
+ ¿рЬ(д) + с(д)
е^ [^ ф](р)«р =
— 00
а(д)(гр)2
+ ¿рЬ(д) + с(д)
е^ [^ ф](р)«р =
1
1
1
1
2
+ТО +ТО
= / + -Л™ / ?[^ф]<* +
-то -то
+то +то
+ ) | е«= а|)^-2=1 е«[7-ф](р)ф =
-то -то
= 2 а(я)Ф"(я) + Ь(д)ф/(д) + с(д)ф(д).
Что и требовалось доказать. ►
Лемма 4. Для любой ф Е 5(К) и любого I > 0 норма
^(^)ф||сим) < еА||Ф||CUR),
где А = яир^ЕМ ф).
= вир
+то
лДП
Ш I е-^е*»(^%))
= вир
+то +то
1 еш
2пе
г»у
— 00 —00
Преобразуем выражение в степени экспоненты:
а(д)£р2
+ гр(д + £&(#) - у) = -
а(д)^р2 - 2гр(д + £&(д) - у)
2
р
2 2ф(д+гб(д)-у)
= - а(д)£
= -а(д)£
Гг12 _ 2ф(д+г6(д)-у) + ( ¿(д+г6(д)-у) )2 _ ( г(д+£Ь(д)-у) )2 р а(д)* + ( а(о)/ ) ( а(о)/ )
= - а(д)£
(р _ ¿(д+^ь(д)-у) )2 + (у-(д+^ь(д)) )2 (р ) + ( )
(у - (д + ¿%)))2 (Р
2а(д )£
¿(д+^Ь(д))-у) )2 а(д)г )
1
2
2
2
2
2
Таким образом,
1
2П
+ТО +ТО
е ш«рф(у)«у =
— 00 —00
+оо 2 (р г(д+*-ь(ч)-у) <2
1 I I 2а(д)г
2п
е
— 00 —00
+ ТО
1 [ (у-(я+мя)))2
= - е 2а(д)4
2п I
(1)г «рф(у)«у =
е а(<^ «рф(у)«у =
+то +то
1 (у-(ч+ЬЬ(д)))2 Г ■
— е 2а(ч)г е
2п У У
-то -то
1
а(ч)*'
аШ «
р
% + гЬ(д) - у)
а(д)г
фЫ«У =
2п у а(д)г
+ТО
(у-(я+£Ь(д)))
2
у/2па(д)г
е 2а(^ ф(у)«у.
В результате получим
(г)ф|
а
< вир
1
+то
(у-(я+£Ь(д)))
< вир
у/2па(д)г 1
е^с(«М е м^Г" ф(у)«у
<
-то
+ТО
/ е |ф(у)|«у <
у72па(д)г
+то
— 00
1 +с(я) [ (у-(я+*ь(я)))2 7
< вир , ==е*с(9) / е «
у72па(д)г
у - (д + ¿ь(д))
|ф|
^ (К) <
< ^ вирдек с(?)||ф|
Лемма доказана. ►
Лемма 5. Для любого г > 0 существует единственное продолжение оператора ^(г) до ограниченного линейного отображения из Сто(К) в Сто(К) с сохранением нормы.
а
1
М Согласно теореме об ограниченном линейном отображении (см. [15]), для ограниченного оператора, определенного на подпространстве банахова пространства, такое продолжение существует и единственно. Оператор Я(£), для любого I Е [0, +то), ограничен как отображение из Б (К) в СТО (К). Кроме того, Б (К) плотно в (К). Следовательно, существует единственное продолжение Я(£) из Сто (К) в себя с сохранением нормы. Это продолжение мы по-прежнему будем обозначать Я(£). ►
Лемма 6. Семейство Я(£) сильно непрерывно на СТО(К). М Пусть сначала ф Е Б (К).
11ш ||Я(*)ф - Я(¿с)ф||сто(м) =
= Нш вир
+ТО
1 Г
еЬс(9М е-^е^+^^ф^ф -
= Нш вир
1 ^ / е-егр(?+Ьо6(^))[^ф](р)^р
— ТО
+ТО
—то
+ТО
1 егм [Т"ф](р)[еад е-^ егрЬЬ(?) -
\/2П
еЬос(?)е-
2
Предел
Нш вир
22 е™[.ТфМе^е-^е^ - е^е-^е^
существует на каждом компакте [р*,р**] х [д*, д**] С К х К и равен нулю, в силу непрерывности функций под знаком предела.
Таким образом, выполнены условия теоремы о предельном переходе в несобственном интеграле и имеет место перестановка предельных переходов:
вир
+ТО
—= [ 11ш егм[Т"ф](р)[еаде-^егрЬЬ(?) - еЬос(?)е-^^егрЬоЬ(9)]ф V 2п ]
2
Интеграл равен нулю, так как подынтегральная функция равна нулю. Таким образом, lim ||F(£)ф — F(¿о)Ф||Соо(R) = 0 для любого ф G S(R).
Так как S(R) плотно в CTO(R), то для любого е > 0 и любой ф G CTO(R) существует функция фе G S(R), такая что ||ф — фб||сте(М) < зМ, где M наибольшая из норм ||F(to) || и ||F(t) ||. Тогда
l|F (^)ф — F Ш|сто(м) =
= ||F(^)ф — F(^)фе + F(to)фе — FВДфб + FВДфб — FВДф||Сто< < ||F(^)ф — F(to^e||cTO(R) + ||F(to)фе — Fйфб||Сто(к) +
+ ||Fффб — F(t)^||cTO(R) < < ||F(to)|Н|ф—фе||сто(м) + ||F(^)фе—F(¿)фе||сет(к) + ||Fй|Н|фе—ф||cto(r) <
< M • — + M • — + M • — = е.
- 3M 3M 3M
Из произвольности е следует что (F(t))t>o сильно непрерывно. ►
Пусть теперь e(t) — гладкая функция монотонно убывающая к нулю при t > 0, t — 0, такая что e(t) = o(t). Например, e(t) = c arctg(t) для некоторого c : 0 < c < 1 diam(G). Определим множество Ge(t) следующим образом: Ge(t) = {q G G : dist(q,dG) > e(t)}. Пусть R(t)(q)}t>o : R — [0,1] — семейство бесконечно дифференцируемых на G функций, таких что фф) (q) = 1 при q G Ge(t), фе(й) (q) = 0 при q G R \ G и фе(^) G [0,1] при q G G \ Ge(t). При t — 0 семейство фф) аппроксимирует индикатор области G в смысле поточечной сходимости. Будем считать, что отображение ^(t,q) = непрерывно по совокупности переменных t и q.
Рассмотрим композицию операторов: F(t) = Fe(t)F(t), где
Fe(t^(q)= фе(4)(9 )ф(q) для ф G Co(G) С CTO(R). Тогда
F(t) : Co(G) — Co(G). При этом, для любого ф G D справедливо представление
t ( )
[FF(t^](q) = фе(*)(^)в*С g I I ^eitP6(q)eiP(q—qo)F(qo)dqodp.
— OO —00
Теорема 2. Пусть а(-), &(•), с(-) Е СЬ(К). Пусть существует сильно непрерывная полугруппа (Т/)Ь>о на пространстве Со (С), разрешающая задачу Коши-Дирихле (4) при т =1. Пусть Л := {ф Е С0ТО(С)| ^ф Е Со (С)} — существенная область определения генератора (Ь1, Оот(Ь1)) этой полугруппы. Тогда семейство операторов Я(£)ь>0 на пространстве С0(С) эквивалентно по Чернову полугруппе, разрешающей задачу Коши-Дирихле (4) при т =1.
М В первую очередь, докажем, что сильная производная в нуле Я(£) совпадает с генератором полугруппы на существенной области определения I.
Покажем, что [Я(£)ф](д) = ф(д) + £^ф(д) + о(£) равномерно по д Е С при £ ^ 0. Используя результат леммы 3 получим:
Я(£)ф(д) = Яе(£)Я (£)ф(д) =
= (£)[ф(д) + 2 А«ф(д) + £&(д^ф(д) + £с(д)ф(д) + о(£)] =
= фе(ь)(ф(д) + 2 Ааф(д) + ¿Кд)Уф(д) + £с(д)ф(д) + о(£))
для любого ф Е I) равномерно по д Е К.
При д Е Се(ь) ^е(ь)(д) = 1 и
[Я(£)ф](д) = ф(д) + 1 Ааф(д) + Ь(д) Уф(д) + с(д)ф(дх) + о(£) равномерно по д Е Се(Ь) при £ ^ 0. При д Е К \ С, фе(Ь)(д) =0 и
[Я(£)ф](д) = 0 = ф(д) + 2Ааф(д) + Ь(д)Уф(д) + с(д)ф(д) + о(£)
равномерно по д Е К \ С при £ ^ 0. Если д Е С \ Се(Ь), то
|Я(£)ф(д) - ф(д) - £Ь1ф(д)| = |(1 - фф)(д))(ф(д) + £Ь1ф(д)) + о(£)|
при £ ^ 0. Функцию ф Е I на множестве С \ Се(Ь) можно разложить по формуле Тейлора в окрестности точки д Е С. Разложение примет вид:
ф(д) = ф(д*) + ф(1)(#д + (1 - 0д*))(д - д*),
где д* Е дС такое что ||д - д*|| < о(£),
|ф(д)| = |ф(1)(^д + (1 - 0д*))(д - д*)| < |ф(1)(^д + (1 - 0д*))|О(£) < о(£)
равномерно по д Е С \ Се(Ь).
Так как ^ф е Сд0 (С), то ^ф(д) ^ 0 равномерно по д е С \ Сф) при £ ^ 0. Следовательно, |^Р(£)ф(д) — ф(д) — ¿Ь1ф(д)| < о(£), равномерно по д Е С \ Сф). Таким образом, для ф е Л равномерно по д е С при £ ^ 0 выполняется равенство ^Р(£)ф(д) = ф(д) + £Ь1ф(д) + о(£). Этим доказано, что сильная производная в нуле на множестве Л, совпадает с генератором разрешающей полугруппы. Для композиции операторов
||^(£)|| = цед*1 (£)|| < ||ВД|| • (£)|| < еь8ир^кСЧ
так как ||*е(£)|| = 1. Таким образом, норма ||^Р(£)|| удовлетворяет условию теоремы Чернова.
В силу условий, наложенных на (фб(г)}, семейство *е(£) сильно непрерывно на (К). Тогда ^Р(£) сильно непрерывно как композиция сильно непрерывных семейств операторов. Так как Со (С) С СТО(К), то ^Р(£) сильно непрерывно на
С0 (С).
Таким образом, выполняются все условия теоремы Чернова и доказана эквивалентность по Чернову семейства (^Р(£))>0 для полугруппы (Т/)^>0, т.е.
Т/ = s- lim [F(t/n)]n.
Тем самым, доказательство теоремы 2 закончено. ►
Из доказанной теоремы следует, что для решения задачи 4 при т =1 справедлива гамильтонова формула Фейнмана:
f (t,qo) =
+TO +TO n n n n a( „.
1 r f А , , ч n E c(qj-i) in E H^-Op -n E ^^^^
1 / / "ГТ- , , n ^ "n ^ "W-u^o n ^ 2
= ьттойп I... 1n)(qj-i)e e e x
(2n )• . ,
— no
Преобразуем вид формулы Фейнмана. Прежде всего заметим, что при п =1
etc(q) / I e-^eip(q+6(q))e-ipy0(y)dydp
— 00 —00
+ТО
etc(q) / e (5)
— 00
y/2na(q)i
Распространяя преобразование на произвольное n, получим
„ „ n ^ V^ ^ b(qj-1)(qj-qj-1) f (t,qo) = Hm .. n^£(t/n)(qj-i)e "j=1 j ej=1 -(qj-1) x
R R j=1
pa(t/n, qo, qi )...pa(t/n, qn-i, qn)fo(qn)dqi...dqn, (6)
tb(q_i)2 1 (qj-1-qj)2
где V (qj-i) = c(qj-i) - ^aj^V и Pa(t,qj-bfc) = , ( , ,0 e- 2ta (qj-1) .
В работе [5] показано, что ни скорость сходимости e(t) ^ 0 (если она o(t)), ни выбор семейства {^6(i)j, аппроксимирующего область при t ^ 0, не влияют на предел в последней формуле, а поэтому
г Г n IV V(q -Л V b(qj-1)(q3 -qj-1)
„/ Л ... -i-T , / , n V (qj-1) a(q4_ 1)
f(t,qo) = hmn / ... J_J_^e(i/n)(qj-i)e j=1 ej=1 j x
R R j=i
pa(t/n, qo, qi)...pa(t/n, qn-i, qn)fo(qn)dqi...dqn =
r r - V V(0- 1) V b(qj-1)(qj-qj-1)
= lim / ... e j=1 ej=1 j x n^^J J
G G
x pa(t/n, qo, qi)...pa(t/n, qn-i, qn)fo(qn)dqi...dqn. (7) Выполняя обратное преобразование при G = (а, ß) С R, получим
+ ТО в + ТО в
n n n a (q • - )p2
^ , 1 f f + f n E c(qj-1) in E b(qj-1)pj -n E-^
f (t,qo) = iim^y J ... J J e j=1 e j=1 e j=1 x
n
-то а -то а
2n
¿Epj (qj-1 -qj)
x e j=1 fo(qn)dqndpndqn-idpn-i...dqidpi. (8)
6. Формула Фейнмана для полугруппы, разрешающей задачу Коши-Дирихле при произвольном т £ [0,1]
Теперь рассмотрим задачу Коши-Дирихле (4) при произвольном т Е [0,1].
Теорема 3. Пусть а(-) Е С62(К), &(•) Е СЬ1(К), с(^) е СЬ(К) и для любого т Е [0,1] существует сильно непрерывная полугруппа (Тт)^>0 на пространстве С0(С), разрешающая задачу Коши-Дирихле (4). Пусть := (ф Е С0°°(С)| ф Е С0(С)} — существенная область определения генератора (Ьт, Бош(Ьт)) полугруппы (Тт)^>0. Тогда решение задачи Коши-Дирихле (4) представляется в виде гамильтоновой формулы Фейнмана:
+<° в +<° в 1 Г Г Г Г |Е [(1—т)2а"(д^--1)+(1—т)Ь'(5^--1)+е(5^--1)]
/т<^> = -^-(2^7 У"УУе" Х
П—т>00
а а
2п
4-V-"
1)Р2
"1 Е [Ь(д^-1)+2(1—т)а'(Ф-1)]р,- — | Е ^^ х е э=1 е ■5-1 х
п
х е"йЛ'",/0(д„^д„^„^д„—
для любой /0 Е С0(С).
М Выше было доказано, что решение задачи Коши-Дирихле (4) при т =1 п представляется формулой Фейнмана
в в П П П а(п . . )р2
„ , , ° Г ° Г 1 й • 1 й «й-Л — 1 й ^^
/(£,®) = ,1'т (2П)П 4 е е 1-1 е х
а а
2п
"Т,Рз (Чо-1—Чо)
е ^ /з(дп)^дп^Рп^дп—1^Рп—1...^д1^Р1.
Заметим, что по лемме 1
д2 д Нт(д, Л) = а(д)^ф(д) + [Ь(д) + 2(1 — т)а'(д)]—ф(д) +
+ [(1 — т )2о"(д ) + (1 — т )Ь'(д)+ с(д )]ф(д) =
д2 д = а(д) ф(д) + ьт (д) ддф(д) + ст (д)ф(д^
где Ьт(д) = Ь(д) + 2(1 — т)а'(д), с(д)т = (1 — т)2а''(д) + (1 — т)Ь'(д) + с(д).
Таким образом, сужение оператора Нт(д, Л) на множество С°°(К) совпадает с 1-квантованием символа Нт(р, д) = а(д)(2р)2 + Ьт(д)ф + ст(д). Так
как а(-) е С£(К), &(•} е С1 (К), с(^) е СЬ(К), то а(д)А(д),Ст(д) удовлетворяют условиям теоремы 2. Тогда из теоремы 2, теоремы 3 и леммы 1 следует гамильтонова формула Фейнмана:
/ (^,дс) = ТТ /с Ы =
+ ТО в +ТО в П
1 Г Г Г Г П Е [(1-Т)2а"(Ф_1)+(1—т)6/(9з._1)+с(9з-_1)] = Пш .^ , ... е 5=1 х
п^то (2п)п } } } } —то а —то а
-v
2n
n n „(„. )p2 n
iП E [b(qj-i)+2(1-T)a'(qj_i)]pj -П E „j1^ i E Pj(qj-i-qj) x e j=1 e j=1 e j=1 x
x /o(qn)dqndpndqn-idpn-i...dqidpi, (9)
где fo G Co(G). ►
Замечание 1. Пользуясь конструкцией гамильтонова интеграла Фейнмана, введенной в работах [23], [7], [1], можно показать, что правая часть гамильтоновой формулы Фейнмана (9) совпадает с интегралом Фейнмана по множеству Eq0 траекторий (q(s), p(s)) в фазовом пространстве [а, в] x R таких, что q(0) = qo:
f (t,qo) = fo (qo) =
t
f /[(1-T )2a"(q(s)) + (1-т )6'(q(s))+c(q(s))]ds
= e0
E?0
t t 2
i/[b(q(s))+2(1-т)a'(q(s))]p(s)ds - / a(q(s)2)p (s)ds 1
x e 0 e 0 /o(q(i)№0(dqdp),
7. Заключение
В настоящей работе рассмотрено семейство параболических уравнений второго порядка, порожденных различными видами квантования квадратичной функции Гамильтона некоторой классической системы. Решение задачи Коши-Дирихле для рассмотренного семейства уравнений на отрезке представлено в виде гамильтоновой формулы Фейнмана, то есть в виде предела конечно-кратных интегралов от элементарных функций при стремлении кратности к бесконечности.
Тем самым, в работе получена новая формула, пригодная для непосредственных вычислений решения поставленной задачи и компьютерного моделирования соответствующей динамики. В работе также обсуждается связь между дифференциальными операторами, соответствующими различным типам квантования квадратичной функции Гамильтона, и связь полученной гамильтоно-вой формулы Фейнмана с интегралами Фейнмана по траекториям в фазовом пространстве.
Авторы благодарят профессора О.Г. Смолянова за полезные дискуссии и внимание. Работа выполнена при финансовой поддержке гранта Президента Российской Федерации MK-943.2010.1. и гранта РФФИ 10-01-00724-a.
Список литературы
1. Boettcher B., Butko Ya.A., Schilling R.L., Smolyanov O. G. Feynman formulae and path integrals for some evolution semigroups related to tau-quantization // Rus. J. Math. Phys., Vol. 18 N 4 (2011), Pp. 381-399.
2. Бутко Я.А. Формулы Фейнмана и функциональные интегралы для диффузии со сносом в области многообразия // Мат. Заметки. 2008 Т.83. N 3. С.333-349.
3. 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. V.151. N. 1.P. 2629-2638.
4. Бутко Я.А., Гротхаус М., Смолянов О.Г. Формула Фейнмана для класса параболических уравнения второго порядка в ограниченной области // Доклады РАН. 2008. Т. 421. N 6. C. 727-732.
5. Butko Ya.A., Grothaus M., Smolyanov O.G. Lagrangian Feynman Formulae for Second Order Parabolic Equations in Bounded and Unbounded Domains //IDAQP 2010. V.13. N. 3. P. 377-392.
6. Бутко Я.А.,Смолянов О.Г., Шиллинг Р.Л. Формулы Фейнмана для фелле-ровских полугрупп //Доклады РАН. 2010.Т.434. №1.С.7-11.
7. 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. V. 50. P. 2009-2018.
8. 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.
9. Chernoff P. Product formulas, nonlinear semigroups and addition of unbounded operators //Mem. Am. Math. Soc. 1974.140.
10. Feynman R.P. Space-time approach to nonrelativistic quantum mechanics //Rev. Mod. Phys. 1948. 20. P. 367-387.
11. Feynman R.P. An Operation Calculus Having Application in Quantum Electrodynamics//Phys. Rev. 1951.84.P.108-128.
12. Gadella M., Smolyanov O.G. Feynman Formulas for Particles with Position-Dependent Mass// Doklady Math. 2007. 77. N.1. P.120-123.
13. Jacob N. Pseudo-differential operators and Markov processes // Imperial College Press. 2001.Vol. 1-3.
14. Obrezkov O., Smolyanov O.G., Truman A., The Generalized Chernoff Theorem and Randomized Feynman Formula //Doklady Math. 2005. 71.N. 1. P.105-110.
15. Рид М., Саймон Б. Методы современной математической физики. Издательство Мир. Москва.1977. т.1
16. Рид М., Саймон Б. Методы современной математической физики. Издательство Мир. Москва.1977. т.2
17. Sakbaev V.G., Smolyanov O.G. Dynamics of a Quantum Particle with Discontinuous Position-Dependent Mass //Dokl. Math. 2010. 82. N. 1. P.630-634.
18. Smolyanov O.G., Feynman type formulae for quantum evolution and diffusion on manifolds and graphs//Quant. Bio-Informatics, World Sc. 2010.3. P.337-347.
19. Smolyanov O.G., Shamarov N.N. Feynman and Feynman-Kac formulae for evolution equations with Vladimirov operator // Doklady Math. 2008. 77. P.345-349.
20. Smolyanov O.G., Shamarov N. N. Hamiltonian Feynman Integrals for Equations with the Vladimirov Operator //Dokl. Math. 2010. 81. N. 2. P.209-214.
21. Смолянов О.Г., Шавгулидзе Е.Т. Континуальные интегралы. Изд-во МГУ. Москва. 1990.
22. Смолянов О.Г., Шавгулидзе Е.Т. Носитель симилектической меры Фейнмана и принцип неопределенности//Доклады А.Н. 1992. T.323.N6. С.1038-1042.
23. Smolyanov O.G., Tokarev A.G., Truman A. Hamiltonian Feynman path integrals via the Chernoff formula// J. Math. Phys. 2002. 43. N. 10. P. 5161-5171.
24. Smolyanov O.G., Weizsacker H.v. and Wittich O., Diffusion on compact Riemannian manifolds, and surface measures //Doklady Math.2000. V. 61. P.230-234.
25. 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. 29.P.589-602.
26. 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. Proc. 7th Intnl. Conf. Evolution Eqs and Appl., Levico Terme, Italy, Oct./Nov. 2000. Birkhauser, Prog. Nonlinear Differ. Eq. Appl. 2003. V. 55. P.349-358.
27. Smolyanov O.G., Weizsacker H.v., Wittich O. Surface Measures and Initial Boundary Value Problems Generated by Diffusions with Drift //Doklady Math. 2007. V. 76. N. 1.P.606-610.
28. Smolyanov O.G., Weizsacker H. v., Wittich O. Chernoff's Theorem and Discrete Time Approximations of Brownian Motion on Manifolds // 2007. Potent. Anal. V. 26. N. 1.P. 1-29.
29. Telyatnikov I.V., Smolyanov O.G., Weizsacker H.v. surface measures generated by diffusions on the set of trajectories in Riemannian manifolds //IDAQP 2008. V. 11. N .1. P.21-31.
electronic scientific and t echnical periodical
SCIENCE and EDUCATION
El № FS77 - 30569. №0421100025. ISSN 1994-0408
Feynman formulae for a family of parabolic equations corresponding to the tau-quantization of a quadratic Hamilton function
77-30569/251251
# 10, October 2011 Ya. Butko, A. Duryagin
Bauman Moscow State Technical University [email protected], [email protected]
In the present article a set of second order parabolic equations, generated by different types of quantization of a quadratic Hamilton function, is considered. The solution of the Cauchy-Dirichlet problem for this set of equations on a segment is represented by a Hamiltonian Feynman formula, that is by a limit of the n-fold iterated integrals of elementary functions, when n tends to infinity. Therefore, a new formula, suitable for the numerical solution of the considered problem and for the computer modeling of the corresponding dynamics, is obtained. The connection between the differential operators, corresponding to different types of quantization of a quadratic Hamilton function, and the connection between the obtained Hamiltonian Feynman formula and phase space Feynman path integrals are also discussed in this article.
References
1. Boettcher B., Butko Ya. A, Schilling R.L., Smolyanov O. G. Feynman formulae and path integrals for some evolution semigroups related to tau-quantization // Rus. J. Math. Phys., Vol. 18 N 4 (2011), Pp. 381-399.
2. Butko Ya.A. Feynman formulas and functional integrals for diffusion with drift in a domain on a manifold //Math. Notes. 2008. V. 83. N. 3-4. P. 301-316.
3. 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. V.151. N. 1.P. 2629-2638.
4. 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. V. 78. N. 1.P. 590-595.
5. Butko Ya.A, Grothaus M., Smolyanov O.G. Lagrangian Feynman Formulae for Second Order Parabolic Equations in Bounded and Unbounded Domains //IDAQP 2010. V.13.N. 3.P. 377-392.
6. Butko Ya.A, Schilling R.L., Smolyanov O.G. Feynman formulae for Feller semigroups //Doklady Math. 2010. V.82. N.2. P. 679-683.
7. 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, 2009-2018. V.50.
8. 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.
9. Chernoff P. Product formulas, nonlinear semigroups and addition of unbounded operators //Mem. Am. Math. Soc. 1974.V.140.
10. Feynman R.P. Space-time approach to nonrelativistic quantum mechanics //Rev. Mod. Phys. 1948. V.20. P. 367387.
11. Feynman R.P. An Operation Calculus Having Application in Quantum Electrodynamics//Phys. Rev. 1951.V.84.P.108-128.
12. Gadella M., Smolyanov O.G. Feynman Formulas for Particles with Position-Dependent Mass// Doklady Math. 2007. V.77. N.1.P.120-123.
13. Jacob N. Pseudo-differential operators and Markov processes // Imperial College Press. 2001.V.1-3.
14. Obrezkov O., Smolyanov O.G., Truman A., The Generalized Chernoff Theorem and Randomized Feynman Formula//Doklady Math. 2005. 71.N. V.1. P.105-110.
15. Sakbaev V.G., Smolyanov O. G. Dynamics of a Quantum Particle with Discontinuous Position-Dependent Mass //Dokl. Math. 2010. V.82. N. 1. P.630-634.
16. Smolyanov O.G., Feynman type formulae for quantum evolution and diffusion on manifolds and graphs //Quant. Bio-Informatics, World Sc. 2010. V.3. P.337-347.
17. Smolyanov O.G., Shamarov N.N. Feynman and Feynman-Kac formulae for evolution equations with Vladimirov operator // Doklady Math. 2008. V.77. P.345349.
18. Smolyanov O.G., Shamarov N. N. Hamiltonian Feynman Integrals for Equations with the Vladimirov Operator //Dokl. Math. 2010.V. 81. N. 2. P.209214.
19. Smolyanov O.G., Shavgulidze E.T. Kontinualnye integraly .in Russian. Moscow State University Press. 1990.
20. Smolyanov O.G., Shavgulidze E.T.. Th support of symplectic Feynman measure and uncertainty principle. //in Russian.Dokl. Acad. Nauk USSR. 1992. V.323. N.
6. P.1038-1042.
21. Smolyanov O.G., Tokarev A.G., Truman A. Hamiltonian Feynman path integrals via the Chernoff formula// J. Math. Phys. 2002. V.43. N. 10. P. 5161-5171.
22. Smolyanov O.G., Weizsacker H.v. and Wittich O., Diffusion on compact Riemannian manifolds, and surface measures //Doklady Math.2000.V.61. P.230-234.
23. 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. V.29.P.589-602.
24. 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. Proc. 7th Intnl. Conf. Evolution Eqs and Appl., Levico Terme, Italy, Oct./Nov. 2000. Birkhauser, Prog. Nonlinear Differ. Eq. Appl. 2003.V.55. P.349-358.
25. Smolyanov O.G., Weizsacker H.v., Wittich O. Surface Measures and Initial Boundary Value Problems Generated by Diffusions with Drift //Doklady Math. 2007. V.76.N. 1.P.606-610.
26. Smolyanov O.G., Weizsacker H. v., Wittich O. Chernoff's Theorem and Discrete Time Approximations of Brownian Motion on Manifolds // 2007. Potent. Anal. V.26.N. 1.P. 1-29.
27. Telyatnikov I.V., Smolyanov O.G., Weizsacker H.v. surface measures generated by diffusions on the set of trajectories in Riemannian manifolds //IDAQP 2008. V.11. N .1. P.21-31.
28. Reed M., Simon B. Methods of modern mathematical phisycs. Academic press New York London. 1972. V.1
29. Reed M., Simon B. Methods of modern mathematical phisycs. Academic press New York London. 1975. V.2