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

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА

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

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

электронный научно-технический журнал

Формулы Фейнмана для параболического уравнения

с бигармоническим дифференциальным оператором

на конфигурационном пространстве

# 08, август 2012

Б01:10.7463/0812.0445534

Бузинов М. С., Бутко Я. А.

УДК 517.987.4

Россия, МГТУ им. Н.Э. Баумана Maxim.cad@gmail.com

Введение

Объектом исследования является параболическое уравнение с бигармоническим дифференциальным оператором Д(Д^) и аддитивным возмущением (дополнительное слагаемое, или "потенциал" в некоторых разделах физики) по переменной конфигурационного пространства. Подобные уравнения находят свое применение в разных областях физики, химии, биологии и компьютерных наук. Так, например, они используются в описании диффузии на поверхности тел [1], описании процессов образования снега [2] и движения веществ в легких [3], описании процессов изменения головного мозга [4], построении поверхностей в компьютерной геометрии [5], восстановлении изображений [6].

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

Термин "формула Фейнмана" в данном контексте был предложен в работе [7]; также в работах Смолянова и его соавторов [7, 8, 9, 10, 11] был предложен метод получения формул Фейнмана для широкого класса эволюционных уравнений. В настоящее время этот метод активно используется для описания классической, квантовой и стохастической динамики на различных геометрических объектах (см., например, [12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23]). Начало данному методу положила работа Ричарда Фейнмана 1948 г. [24], в которой рассматривалось уравнение Шредингера с потенциалом и (на эвристическом

уровне строгости) было получено решение задачи Коши для такого уравнения в виде предела конечнократных интегралов по декартовым степеням конфигурационного пространства соответствующей классической системы при стремлении кратности к бесконечности. Строгое математическое обоснование результатов Фейнмана было приведено Нельсоном [25] на основе формулы Троттера. Предложенный Смоляновым и его соавторами подход опирается на теорему Чернова [26], существенно обобщающую формулу Троттера.

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

Рассмотрим задачу Коши следующего вида:

öw. . т . . — х) = Lw(t x),

w(0, х) = w0(x).

Здесь L — некоторый линейный оператор, действующий на функцию w, зависящую от времени t > 0 и переменной х Е Q, где Q — некоторое множество, которое будем называть конфигурационным пространством эволюционной системы соответствующей задаче Коши (1).

Пусть L(X) — пространство всех линейных, непрерывных операторов T: X ^ X, где X — банахово пространство, а Dom(T) = X — область определения оператора T.

Определение 1. Семейство ограниченных линейных операторов (Tt)t>0 в X называется сильно непрерывной операторной полугруппой, если T0 = Id, Tt+s = TtTs для t, s > 0, x Е X, а также семейство сильно непрерывно, т.е. имеет место равенство lim ||Ttx — x||X = 0.

Оператор L := lim "X—~ с областью определения

Ttx — x'

Dom(L) := {x Е X: 3 Lx := lim }

называется генератором или производящим оператором полугруппы (Т4)4>0.

Можно показать [30], что для корректно поставленной в банаховом пространстве X задачи Коши (1) с начальным условием ^(0,я) = ^0(я), € Бош(Ь) решение представляется в виде х) = Т^0(х). Таким образом решение задачи Коши (1) равносильно построению сильно непрерывной полугруппы операторов (Т4)4>0 с заданным генератором Ь.

Явный вид полугруппы операторов (Т4)4>0, разрешающей задачу Коши (1), часто невозможно найти. В этом случае используются методы аппроксимации полугрупп. Будем использовать метод приближения, основанный на теореме Чернова.

Определение 2. Производная в нуле функции Р: [0, е) ^ £(Х), е > 0, — это линейное отображение Р'(0): Бош(Р'(0)) ^ X, определяемое следующим образом:

V Р(*) X - Р(0) X Р'(0)ж := 11т —^-, (2)

где Dom(F'(0)) — векторное пространство элементов из X, для которых данный предел существует (как сильный предел).

Теорема 1 (Чернов [12]). Пусть X — банахово пространство, F: [0; то) ^ L(X) — (сильно) непрерывное отображение, такое, что F(0) = Id и ||F(t)|| < eMt для некоторой константы M G R и t > 0. Пусть C — такое линейное подпространство Dom(F'(0)), что сужение оператора F'(0) на C замыкаемо. Пусть (L; Dom(L)) = F'(0) \ C — соответствующее замыкание. Если (L; Dom(L)) является генератором сильно непрерывной полугруппы (Tt)t>0, то для всех t0 > 0 последовательность операторов (F(t/n))n|neN сходится к (Tt)t>0 при n ^ то в сильной операторной топологии равномерно по t G [0; t0], т.е.

Tt = lim (F(*)Г (3)

Vn//

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

Определение 3. Равенство (3) будем называть формулой Фейнмана. Формула Фейн-мана называется гамилътоновой, если для всех t > 0 оператор F(t) является псевдодифференциальным оператором. Формула Фейнмана называется лагранжевой, если для всех t > 0 оператор F (t) является интегральным оператором, ядро которого выражается через элементарные функции.

Следующая теорема позволяет получать формулы Фейнмана для суммы (в смысле обычного сложения в L(X)) операторов, действующих на функцию.

Теорема 2 (Формула Фейнмана для аддитивных возмущений [35]). Рассмотрим банахово пространство (X, || ■ ||X) с соответствующей нормой. Пусть {Ttk}t>0, k = 1, ..., n, —

TA TT Tfc — - w

сильно непрерывные полугруппы на X с генераторами Lkи = lim-, и своими обла-

стями определения D(Lk) = G X: 3 lim T ^—— G . Предположим, что опера-

п

тор L = Li + L2 +... + Ln с областью определения D = П D(Lk) замыкаем, и что замыкание

fc=i

этого оператора является генератором сильно непрерывной полугруппы операторов (Tt)t>0 на X. Пусть Fk (t), k = 1, ..., m, — семейство операторов в X, эквивалентных по Чернову полугруппе (Ttk)t>0, т.е. для каждого k = 1, ..., m имеет место Fk(0) = Id, ||Fk(t)| < eMkk для некоторого Mk > 0. Пусть также Dk = core(Lk) С D(Lk), где core(Lk) — существен-

r т т-л т Fk (t)<p — ш т

ная область оператора Lk, и для всех ^ G Dk выполнено lim ----Lk^

m

Предположим, что существует D = core(L) и вместе с тем D СП Dk. Тогда семей-

fc=i

ство F(t) = F1 (t) о F2 (t) о... о Fn(t) эквивалентно по Чернову полугруппе (Tt)t>0 и для t > 0 имеет место формула Фейнмана

Tt = lim (f(-

= 0.

x

V, n

в сильной операторной топологии.

2. Задача Коши для параболического дифференциального уравнения с бигармоническим оператором на конфигурационном пространстве и разрешающая ее полугруппа

Перейдем к следующей задачи Коши иа конфигурационном пространстве Кп:

п п

-Д(Д^, х) = - £ £ ^^^ х).

¿=1 7=1 ® 3

Поставленная задача Коши (4) будет решаться в банаховом пространстве X = Ь2(Мп) квадратично интегрируемых функций со стандартной нормой

На конфигурационном пространстве Кп рассмотрим норму ||х|| = . Мы предпола-

гаем, что ^о, ■) € X для всех Ь > 0.

Сперва мы найдем семейство операторов, разрешающее задачу Коши (4) и обладающее свойствами сильно непрерывной полугруппы, которое действует в пространстве Б (Кп), всюду плотном в Ь2(Мп), и продолжим это семейство операторов с необходимыми свойствами на объемлющее пространство Ь2(Кп).

Для быстро убывающих функций определены прямое преобразование Фурье

К"

Будем считать, что , ■) € Б (Еп). Применим преобразование Фурье к уравнению задачи Коши (4). Получим

г,3

с соответствующими начальными условиями о;(0, -)(у) = ^о(у). Решив это ОДУ относительно Ь, получим его общее решение в области переменной у € Кп:

и обратное преобразование Фурье

(у) = ехр(- £ у2у2) С = е-«4 С.

Используя обратное преобразование Фурье, с учетом начальных условий, получаем явное решение задачи Коши (4):

и(£, х) = —/ е-'"У"4¿у. (5)

К"

Полученное решение (??) задает зависящее от £ семейство псевдо-дифференциальных операторов с символом ^(у) = ехр(—£ ^ уг2у2), определенных для всех и € Б(Ега) и £ > 0.

г.

Обозначим это семейство [И^ таким образом:

Я4[и](х) = (2^У е-4"у"4ег(х,у)и(у) ¿у.

К"

Теорема 3. Семейство операторов [И^ 0 задает сильно непрерывную полугруппу операторов действующую в пространстве Ь2(Кга).

Пусть Ц и Р — псевдо-дифференциальные операторы с символами р(у) и д(у). Тогда для любой и € Б(Кп)

Р о Ц[и](х) = —^т /р(уМу)ег(х,у)и(у) ¿у.

К"

В самом деле,

Р о Ц[и](х) = Р[Ц[и]](®) =

¿у ¿в =

1- [ р^в^ / е-^у)

(2п)3™/2

д(г)вг(у'г)и(г) ¿г 1

(2п)'

р(8)д(8)е'(ж'в)и(8) ¿в.

Заметим также, что Л,5(у)Л,т(у) = ехр( —(^ + т) ^ у2у|) = Л,5+т(у), з,т € К. Тем самым

г.

показано, что для всех и € Б(Кп) и ^ > 0 выполнено И о Ит = И5+т; также очевидно, что И0 = И.

0

Функции и, ехр(—£ принадлежат Б(Кп) и их произведение тоже принадлежит

г,.

этому пространству. Между тем (И4)4>0 можно представить так:

(Я4)М = ^ -1[е-«4 и(у)].

Прямое и обратное преобразование Фурье отображает Б(Кп) в Б(Кп), следовательно семейство операторов (И^ 0 тоже отображает Б (Кп) в Б (Кп). Мы показали, что (И^ 0 задает полугруппу операторов действующую в пространстве Б (Кп)

Покажем теперь, что (И4)4>0 продолжается до полугруппы операторов на банаховом

пространстве Ь2(Кга). Для всех £ > 0, п € N символ Л.4(у) = ехр(—¿^уг2у2) семейства

г,.

операторов (И4)4>0 принимает значения в [0; 1], поэтому для любой и Е Б(Кп), согласно теореме Планшереля, имеет место следующая оценка:

HHiMlb = ||F [ВДЦ2 = < ||и||2 = 1Mb.

Следовательно, для любой последовательности Коши {un} G S(Rn), при t G [0; и

||un — Um||2 < £ выполнено

||Н[и„] — Ht[wm]|2 = ||Ht[Wn — Wm]|2 <

t>0

Учитывая равенство ¿2(Ега) = Б(Ега), можно сказать, что семейство операторов (И) продолжается до ограниченного семейства операторов на Ь2(Кга). При этом (И4) 0 сохраняет свойство полугруппы операторов. Действительно, Б (Кп) всюду плотно в Ь2(Мга), следовательно при фиксированных т, 5 > 0 для любых и Е Ь2(Кга) и е > 0 можно указать такое ие Е Б(Мга), что ||и — и£|| < ^м, где М = шах(||И8+т||, ||И8 о Ит||). Поэтому при е > 0:

11И о Ит и — И5+т и 11 = 11И о Ит и — И о Ит и£ + И о Ит и£ — И5+т и£ + И8+т и£ — И8+т и 11 ¿2 < < ||И8 о Ит || ■ ||и — и£|| + ||И о Ити£ — И5+ти£|| + ||И8+т || ■ ||и£ — и|| < е.

Таким образом, (И4) 0 задает полугруппу операторов на Ь2(Кга), и для доказательства ее сильной непрерывности достаточно проверить это свойство в ¿0 = 0. Для всех и Е Б (К) имеем:

lim

H (t)M — и

lim

F -1[(h и)] — F -1[и ]

lim

F -1[(h — 1)и ]

lim t0

ei(x'y)(e-i"y"4 — 1)u(y) dy = limt||F-1^w] ||2 = 0.

t0

В самом деле, пользуясь разложением аналитической функции в ряд Тэйлора запишем:

exp

—11УЧ2) — 1 = (— У2У? + 2Еу»2У?)2exJ—s J]y2y2))t = ф(y) t

»j

■V

■ 2

где s G (0,t). Функция ф принадлежит S(Rn). Следовательно, фи G S(Rn) и, таким образом, ||F-1[фи]||2 = const.

Теперь, аналогично тому, как было доказано полугрупповое свойство, можно показать,

что limllHM — и||2 = 0 при и G L2(Rn). ►

Предложение 1. Генератором полугруппы (Ht) 0 служит замыкание (L, Dom(L)) опе-

д2 д 2 > ратора L = — ^ ——2, заданного на множестве S (Rn).

dXi dxj

Л Множество S (Rn) инвариантно относительно действия операторов полугруппы, а значит, согласно [28, Теорема X.49], является существенной областью определения оператора L.

2

2

2

Поэтому достаточно показать, что равенство lim

H«И-Ьш 2

0 выполнено при

и е Б(Ега).

Воспользуемся свойствами преобразования Фурье и разложением аналитической функции в ряд Тэйлора с остаточным членом в форме Лагранжа. Тогда при и е Б(Ега)

lim

t-> о

- ш

— Ьш

lim

t-> о

F-1[(htш)] — F-1[ш]

t

-Г

-1

ш

»J

Здесь

Поскольку £ G

lim

t-> о

1

я

1

ht — 1 ! sr^ 2 2

z^yjш

— I —

t

»j

lim t

t-> 0

F -1[£Ш ]

= 1 (J2 y2y|) exp ( — s £ y2y2 ) , s G (0, t)

ij

'),тои G S(Rn). Таким образом, ||F-1[£(y)w]||2 < то. ►

Тем самым показано, что для любой функции и0 е Бош(Ь), функция

ш(t, x)

1

(2п)n

e-t||y"4ei(x,y) шо(у) dy,

(6)

является решением задачи Коши (4). При этом

() = core(L), т.е.

() является существенной областью определения оператора (L,Dom(L)), где (L, Dom(L)) замыкание опера-

д 2 д2

тора — ^ —2 —2 S(Rn). Кроме того, для любой функции ш0 G L2(Rn) функция ш, определенная формулой (5), также является элементом L2 (Rn). Иногда ее называют решением "mild solution" [29] задачи Коши (4).

3. Гамильтонова и лагранжева Формулы Фейнмана для параболического дифференциального уравнения с бигармоническим оператором

Построим для задачи Коши (4) гамильтонову Формулу Фейнмана. Пусть Д(А, Ь) — резольвента оператора Ь задачи Коши (4). Тогда для Я(А, Ь) верно следующее утверждение.

Теорема 4 ([30]). Для любой сильно непрерывной полугруппы (Т4)4>0 на банаховом пространстве X с генератором (А, Бош(А)), имеет место равенство:

Ttx = lim (Пя(П, a)

t n—^ I t V t /

x = lim (Id--A

n—V n

x,

(7)

Обозначим Е(¿) = (Id — ¿Ь) 1 и представим Е(¿) в виде интегрального оператора. С помощью преобразования Фурье [31] получим следующий результат:

E (t)M(®) = (Id — tL)->o](x)

1

ei(x,y)

ш0(у)

(2n)n/2J 1+ t||y||4

dy,

т.е. E(t) представляет собой семейство псевдо-дифференциальных операторов с символом et(y) = ТГТ^Щ.

»J

2

2

t

2

0

n

n

Предложение 2. Семейство операторов E(t) продолжается до ограниченного по норме семейства операторов из L2(Rn) в L2 (Rn).

Л Заметим, что sup |et(y)| < 1. Поэтому норма операторов Et оценивается таким же

ye Rn

образом, как и при доказательстве ограниченности семейства операторов Ht в теореме 3. ►

Теорема 5. Решение u(t, x) задачи Коши (4) для любого и0 G L2(Rn) представимо в виде гамильтоновой формулы Фейнмана:

u(t, x) = lim (e(-))n[u0](x)

га^те V \И//

r. r. m

■■■ п

lim -—-т^те (2ffV

fc=1

1 + -L ||y

m

14

i Y1 uk •(vfc— Vfc_l)

e fc=1 u0(vm) dv1 du1... dvm dum, (8)

где х = -и0 Е

Утверждение непосредственно следует из предложения 2, и теоремы 4. ►

Для использования теоремы 2 в дальнейшем нам понадобится следующее утверждение. Предложение 3. Отображение Е(¿): [0; то) ^ £(Х), для любого и Е еоге(Ь) = Б(Кп)

удовлетворяет равенству lim

E(t)[w] — w

— Lu

0.

Л Согласно свойствам преобразования Фурье, для всех u G

) имеем:

lim

t-> 0

E (t)[u] — и

— lu

lim

t0

lim

t0

F -1[(etU)] — F -1[и ]

— F -1[(—E(y))U ]

Г

-1

—S(y) 1 + tE(y)

+ Е(у) и

lim t

t0

f -1[£и ]

Здесь E(y) = Y1 Уг2^! и на Rn выполнены соотношения £(y) = --

i,j 1 + t^(y)

вытекает, что £u G S(Rn) и, следовательно, ||F-1 [£(y)u]||2 < то. ►

< E2(y). Отсюда

Теперь мы построим формулу Фейнмана второго типа — лагранжеву формулу. Мы построим ее для задачи Коши (4) на прямой К, и аналогичным образом на пространстве К3. Представим оператор Е (¿) = (И — ¿¿)-1 в виде виде оператора свертки с функцией

1 Г е»(х,у)

Et(x)

(2n)n

1 + t| y|

dy

и рассмотрим задачу Коши (4) на прямой К. Функция Е будет выглядеть так:

1

Et(x) = —

2п J 1 + ty

-dy.

Вычислим интеграл, стоящий в представлении функции ЕДж). Подынтегральная функция имеет четыре полюса:

V2

V2.

У1,2 = -Г"(1 ± ^ Уз,4 = -Г" ( —1 ± i).

2t 4

2t 4

1

t

2

2

t

t

2

2

0

По теореме Коши о вычетах и лемме Жордана [32], верны равенства

1

2п J 1 + ty

-dy = lim —

2п ,/ 1 + tz4

dz = i ^^ Выч

L1 + ty4

где 7 — контур, составленный из отрезка [—г, г] и дуги окружности |г| = г, в области б1§п(х) > 0, а V — область ограниченная этим контуром. При х > 0 внутри контура 7

находятся полюсы у13, при х < 0 — полюсы у24. Следовательно, при х > 0

1

2п / 1 + tz4

dz =

¡pixy

ie

Чу — J4 (1 —i)) (y — (—1 + i)) (y — J4 (—1 — i))

+

y= 2^/4 (1+i)

ixy

+

ie

t(y — (1 + i)) (y — J4(1 — i)) (y — 2t1/4(—1 — i))

y= 2^74 (-1+l)

e 2t 4

2^2t 4

x

сое

V2

2t 4

+ в1п

x

V2

2t -4

Аналогично при х < 0 с учетом направления обхода контура получим:

1

xV2 1

2п У 1 + ¿г4 В результате для любого х е К

dz

e 2t4

2^2t 4

x

сое

V2

x

2t 4

— Э1П

V2

2t 4

Et(x) = e-C|x|(cos C|x| + sin C|x|),

где

с = 4

2* 4

Введем такое семейство интегральных операторов С1(*), что для любой и0 е Ь2

и t 0

G1 (t)[ш](x) = (Et * шо)(x) = (шо * Et)(x) =

C

с e-C|y|(cos C|y| + sin C|у|)шо(x — y) dy

= 2 I e N(cos |z| + sin |z|) шо(x — Cj dz.

e

y

e

e

Согласно теореме 4, С1(^) ~ (И4)4>0; здесь (И4)4>0 — найденная нами выше полугруппа, разрешающая исходную задачу Коши (4) на числовой прямой. Итак, справедливо следующее утверждение.

Предложение 4. В одномерном случае решение задачи Коши (4) при и0 Е Ь2 (К) может быть предствалено с помощью лагранжевой формулы Фейнмана:

w(t,x) = lim (üJ[wo](x)

га^те V \n//

n n

1 í* f —' | | t 1/4

Иш - / .. ek=1 k ^(cos|ufc| +sin|ufc|)wo(un+i^ufc)dui...du„, (9)

m m fc=i - fc=i

где х = ига+1.

Теперь рассмотрим задачу Коши (4) на пространстве К3 и получим для нее лагранжеву формулу Фейнмана. Функция Е в этом случае будет такой:

1 Г 6»(х,У)

£'<Х) = (2Л)3У ТГ«=

1 Г Г Г ег(х1У1+Х2У2+хзУз)

dyi dy2 dy3.

(2п)3 J J J 1 + t(y4 + y4 + y4 + y?y22 + y22y32 + y2y2)

R R R

Воспользуемся сферическими координатами:

yi = r sin $ cos $2, y2 = r sin $ sin $2, y3 = r cos $2

(при этом r2 = y2 + y2 + y2), а также четностью функции et и запишем интегральное выражение функции Et(x) в новых координатах:

п 2п

2 i' i' i' r2 sin $

Et(x) = 7—— expfirfe sin $ cos $2+x2 sin $ sin $2 +x3 cos -— drd^i d$2.

(2n)3 J J J 1 + tr2

0 0 R

Обозначим t = xi sin cos $2+x2 sin sin $2+x3 cos и найдем интеграл по переменной r.

Q2 , ч

Используя свойство F-i[y2ra/(y)](x) = (—(^F-i[/(y)](x)J, с учетом предыдущего результата для прямой R получаем:

r.2n írT ^2 Г ^irr

J (T ) = ~r——rdr = — т^ / —e-rdr

w I 1 + tr2 дт 2 / 1 + tr2

C д 2

C^(e-C|T|(cos C|t| + sin C|t|)) = C3e-C|T|(cos C|t| — sin C|t|), 2 дт2

где

C ^

2t1/4"

Вернемся к вычислению исходной функции Et(x) в сферических координатах. Используя равенство [33]

п 2п 1

У J f (т) sin tfid^A = 2* J f (Rs)ds, 0 0 -1

где R ^ \/СГ+С2Г+Сз, получаем:

п 2п

Et(x) = (2C33 J J ^e-C|T '(cos C It I - sin C |т |)) sin ^ d^ d^ = 00

i

2C3 Г

= ^^ e-C|Rs| (cos C|Rs| - sin C|Rs|) ds = (2n)2 J -1

CR

C2 [' 2C2

e-s (cos s — sin s) ds = e-CR sin(CR) =

Rn^^ ' R(2n)

0

t-1/2exp(—-(/ЧР

'V 5'2+Й+a

Введем такое семейство интегральнык операторов G3(t), что при wo G L2(R) и t > 0 G3(t)[w](x) = (Et * wo)(x) = (wo * Et)(x) =

_ Г Г Г t-1/2exp(—y^^+P) -(^^P)

^л/ёГ+ёГ+ё:

Wo(x — ё) d6 d6 d&

= ^ гехр(—нею81п(И^1)ио(х — ^2*1/40

п2 У IIе I

к3

где е = (е1,е2,ез) е к3.

Снова по теореме 4 имеем Сз(£) ~ (Н4)4>0, где (Н4)4>0 найденная выше полугруппа, разрешающая исходную задачу Коши (4) на пространстве К3. Таким образом, доказано следующее утверждение.

Предложение 5. Решение задачи Коши (4) для и0 е Ь2(К3) может быть получено с помощью следующей лагранжевой формулы Фейнмана:

/ * \ \ п -тз / 1 ^ г

w(t, x) = lim (Wn[wo](x) = Vn//

J™ ^ Í ... í exp(— Ы)П ^ Wo(un+1 — £ ufc)du1 ...dun, (10)

га^те n2i J J \ ¿—' / J-J- ||Ufc ||

где x = ui+1, x = (x1, x2, x3) G R3.

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

Рассмотрим схожую с (4) задачу Коши вида

ddt (t, x) = —A(Aw)(t, x) + V (xMt, x), (11)

w(o, x) = ^o(x),

где V(x) — ограниченная и непрерывная в Rn функция со значениями на сегменте [—а, а], а > 0. Здесь оператор умножения на функцию V можно рассматривать как ограниченное аддитивное возмущение [27], оператора L = — A(A). Пусть оператор L = L + V, с областью определения Dom(L) = Dom(L). Тогда, очевидно, оператор (L,Dom(L)) снова является генератором некоторой сильно непрерывной полугруппы (Tt)t>0 на L2(Rn).

Предложение 6. Если символ s(y) действующего в пространстве L2(Rn) псевдо-диффе-ренциального оператора S зависит от одной переменной и sup |s(y)| G [0; 1], то ||S|| < 1.

y€Rn

M По теореме Планшереля для пространства L2 (Rn)

||S|| = sup ||F-1(sw)||2 = sup ||sw|2 < |МЬ = 1.

(IM|2 = 1) (||wo||2 = 1)

Замечание 1. Из предыдущего предложения сразу следует, что нормы семейств операторов Ht, G1(t), G3(t), E(t) не превышают единицы.

Из теоремы 2 об аддитивных возмущениях вытекает следующая теорема, основная в этой работе.

Теорема 6. Семейства операторов (etV о (Ht)t>o), (etV оE(t)), (etV о G^t)), (etV оG3CO) эквивалентны по Чернову полугруппе (Tt) 0, разрешающей задачу Коши (??). Здесь etV — это оператор умножения на функцию etV(x).

Следовательно, теперь, мы можем записать выражения формул Фейнмана для задачи Коши (??).

Гамильтоновы формулы Фейнмана для w0 G L2(Rn): > V о (Я il^ I IWol(x)

(t \n enV о (Яt). J M(x)

n n - /

lim ^ I. . . I exp

n

in Rn

u + V(vfc))! eiuk(vk-vk-l)^0(Vn)du1 dv1... dundvn, (12)

Ln • fc=1

2n

n

w(t, x)= lim fenV о E f-^ U](x) = n^+те у \nJ J

1 Г Г n 1 Г t n 1

li^ ——n /.../Пmr^exp -Vv(vfc) eiuk(vk-Vk-l)W0(Vn)du1 d«1 ...dundvn, (13) (2n)nJ J H-1 + -uk Ln^f J

In Rn k=1

2n

где x = V0.

t

Лагранжева формула Фейнмана задачи Коши (??) на числовой прямой R при wo G L2

w(—,x) = lim (enV о Gif[wo](x) = га^+теV Vu//

n у- n n

1 f f E -K |+-t E V(vfc) n n ;

lim 2- / ■■ ek=1 n fc=1 Д (cos |vfc | + sin |vfc |) x

га^те 2n

TtD TtD k=i

— i n

x ^o(v„+i - V2(uu)4 £ vfc)dvi... dv-, (14)

и

fc=i

где x = Vn+i.

И, наконец, лагранжева формула Фейнмана задачи Коши (??) на пространстве R3 при

wo G L2(R3):

w(—, x) = lim fenV о üJ[wo](x) га^+те V Vn//

Iii - ™ ь || || -

= ralim -2- ... exp(-£ ||ufc || + V(ufc)) Д Sin Uf wo(u„+i - £ ufc) dui ...du-, (15)

ttd4 TTD4 k—i k—i k—i

где x = un+i, x = (xi, x2, x3) G R3.

Замечание 2. Равенства в полученных формулах Фейнмана следует понимать в смысле пространства L2(Rn), т.е. в смысле равенства

lim ||w(—,x) - Ф„(—,ж)||2 = 0,

га^те

где Фга(—, x) обозначает выражение под знаком предела в правой части формулы Фейнмана.

Замечание 3. Как следует из результатов работы [34], гамильтонову формулу Фейнмана (7) можно интерпретировать как интеграл Фейнмана по траекториям в фазовом пространстве. Следовательно, справедливо следующее представление полугруппы (Tt)t>o с помощью интеграла Фейнмана

= T мм = / e-'«' ^V0 „o(q(—))*№*), (16)

Ef

где пространство Etx и псевдомера Фейнмана Ф^ определены в работе [34].

Интеграл Фейнмана в формуле (??) может быть вычислен для w G L2(R) с помощью лагранжевой формулы Фейнмана (??), а для w G L2 (R3) — с помощью лагранжевой формулы Фейнмана (??).

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

В статье рассмотрено эволюционное параболическое уравнение с бигармоническим дифференциальным оператором по переменной конфигурационного пространства и соответствующая этому уравнению задача Коши на всем конфигурационном пространстве. Доказано существование сильно непрерывной полугруппы операторов на пространстве L2(Rn),

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

Надо отметить, что полученные лагранжевы формулы Фейнмана благодаря своему элементарному виду позволяют проводить численное моделирование динамики эволюционной системы, и, тем самым, являются новым инструментом численного моделирования эволюционных уравнений. Новизна такого метода располагает к разработке оптимальных вычислительных алгоритмов, реализующих формулы Фейнмана на ЭВМ, в том числе для параллельных вычислений. В то же время гамильтоновы Формулы Фейнмана связаны с интегралами Фейнмана по траекториям в фазовом пространстве квантовой системы. Такие интегралы являются важными объектами квантовой механики.

Исследования второго автора поддержаны федеральной целевой программой "Научные и научно-педагогические кадры инновационной России" на 2009-2013 гг. (соглашение № 14.B37.21.0370), а также грантом Президента Российской Федерации MK-4255.2012.1.

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

1. Tritscher P. An integrable fourth-order nonlinear evolution equation applied to surface redistribution due to capillarity // J. Austral. Math. Soc. Ser. B. 1997. V. 38, No4. P. 518-541.

2. Myers T.G., Charpin J.P.F. A mathematical model for atmospheric ice accretion and water flow on a cold surface // Int. J. Heat Mass Transf. 2004. V. 47. № 25.

3. Halpern D., Jensen O.E., Grotberg J.B. A theoretical study of surfactant and liquid delivery into the lung // J. Appl. Physiol. 1998. V 85. P. 333-352.

4. Toga A. Brain Warping. N.-Y.: Academic Press, 1998.

5. Monterde J., Ugail H. A general 4th-order PDE method to generate Bezier surfaces from the boundary // Comp. Aid. Geom. Des. 2006. V. 23. P. 208-225.

6. Kim S., Lim H. Fourth-order partial differential equations for effective image denoising // El. J. of Differential Equations. V. 17. P. 107-121. URL: http://ejde.math.txstate.edu (30.07.2012).

7. Smolyanov O.G., Tokarev A.G., Truman A. Hamiltonian Feynman path integrals via the Chernoff formula // J. Math. Phys. 2002. V. 43, № 10. P. 5161-5171.

8. 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.

9. Smolyanov O.G., Weizsäcker 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.

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

11. 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, № 1. P. 606-610.

12. Богачев В.И., Смолянов О.Г. Действительный и функциональный анализ. Университетский курс. М.-Ижевск: НИЦ «Регулярная и хаотическая динамика», Институт компьютерных исследований, 2009. 724 с.

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

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

15. 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, № 1. P. 2629-2638.

16. 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, № 3. P. 377392.

17. Gadella M., Smolyanov O.G. Feynman Formulas for Particles with Position-Dependent Mass // Doklady Math. 2007. V. 77, № 1. P. 120-123.

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

19. Sakbaev V.G., Smolyanov O.G. Dynamics of a Quantum Particle with Discontinuous Position-Dependent Mass // Doklady Math. 2010. V. 82, № 1. P. 630-634.

20. 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.

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

22. Smolyanov O.G., Shamarov N.N. Hamiltonian Feynman Integrals for Equations with the Vladimirov Operator//Doklady Math. 2010. V. 81, №2. P. 209-214.

23. 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.

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

25. Nelson E. Feynman integrals and the Schrodinger equation // J. Math. Phys. 1964. V. 3. P. 332-343.

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

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

28. Рид М., Саймон Б. Методы современной математической физики. Т. 2. М.: Мир, 1977. 393 с.

29. Pazy A. Semigroups of linear operators and Applications to partial differential equation. Springer-Verlag, 1983. 279 p.

30. Engel K.-J., Nagel R. One-Parameter Semigroups for Linear Evolution Equations. Springer, 2000. 609 p.

31. Зорич В.А. Математический анализ. Ч. II. М.: МЦНМО. 2002. 789 с.

32. Свешников А.Г., Тихонов А.Н. Теория функций комплексной переменной. М.: ФИЗМА-ТЛИТ, 2005. 335 с.

33. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: ФИЗМАТЛИТ, 1963. 1097 с.

34. 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. 2011. V. 18, №4. P. 381-399.

35. Butko Ya.A., Schilling R.L., Smolyanov O.G. Lagrangian and Hamiltonian Feynman formulae for some Feller semigroups and their perturbations // IDAQP. 2012. P. 1-19.

SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU

SCIENCE and EDUCATION

EL № FS77 - 48211. №0421200025. ISSN 1994-0408

electronic scientific and technical journal

Feynman formulae for a parabolic equation

with biharmonic differential operator

on a configuration space

# 08, August 2012

DOI: 10.7463/0812.0445534

Buzinov M. S., Butko Ya. A.

Russia, Bauman Moscow State Technical University

Maxim.cad@gmail.com

The Cauchy problem for a parabolic partial differential equation with biharmonic operator and additive perturbation is considered in this note. Such equations are used in different domains of physics, chemistry, biology, and computer sciences. The solution of the considered problem is represented by Feynman formulae, i.e. by limits of iterated integrals of elementary functions when multiplicity if integrals tends to infinity. The main part of these formulae is proved with the help of Chernoff's theorem; some formulae are obtained on the base of the Yosida approximations. Different types of Feynman formulae are presented in this work: Lagrangin and Hamiltonian. Lagrangian Feynman formulae are suitable for computer modeling of the considered dynamics. Hamiltonian Feynman formulae are related to some phase space Feynman path integrals; such integrals are important objects in quantum physics.

References

1. Tritscher P. An integrable fourth-order nonlinear evolution equation applied to surface redistribution due to capillarity. J. Austral. Math. Soc. Ser. B, 1997, vol. 38, No4. P. 518-541.

2. Myers T.G., Charpin J.P.F. A mathematical model for atmospheric ice accretion and water flow on a cold surface. Int. J. Heat Mass Transfer, 2004, vol. 47. no. 25, pp. 5483-5503.

3. Halpern D., Jensen O.E., Grotberg J.B. A theoretical study of surfactant and liquid delivery into the lung. J. Appl. Physio., 1998, vol. 85, pp. 333-352.

4. Toga A. Brain Warping. N.-Y., Academic Press, 1998.

5. Monterde J., Ugail H. A general 4th-order PDE method to generate Bezier surfaces from the boundary. Comp. Aid. Geom. Des., 2006. vol. 23, pp. 208-225.

6. Kim S., Lim H. Fourth-order partial differential equations for effective image denoising. El. J. of Differential Equations. 2009, conf. 17, pp. 107-121. Avaliableat: http://ejde.math.txstate.edu, accesed 18.08.2012.

7. 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.

8. 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. Vol/ 2. A Volume in Honor of Sergio Albeverio. Canadian Math. Society. Providence, AMS, 2000, pp. 589-602. (CMS Conference Proceedings, vol. 29).

9. Smolyanov O.G., Weizsacker H.v., Wittich O. Chernoff's theorem and the construction of semigroups. Proc. of Evolution Equations: Applications to Physics, Industry, Life Sciences and Economics: EVEQ2000 Conference, Levico Terme, Italy, Oct./Nov. 2000. Birkhauser Verlag, 2003. P. 349-358. (Progress in Nonlinear Differential Equations and Their Applicftions, vol. 55).

10. Smolyanov O.G., Weizsacker H.v., Wittich O. Diffusion on compact Riemannian manifolds, and surface measures. Doklady Math., 2000, vol. 61, pp. 230-234.

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

12. Bogachev V.I., Smoljanov O.G. Dejstvitel'nyj i funkcional'nyj analiz: universitetskij kurs. Moscow, Izhevsk, NIC "Reguljarnaya i haoticheskaya dinamika" [Scientific Research Center "Regular and Chaotic Dynamics"], Institut komp'juternyh issledovanij [Institute of Computer Science], 2009. 724 p.

13. Butko Ya.A. Feynman formulas and functional integrals for diffusion with drift in a domain of a Riemannian manifold. Math. Notes, 2008, vol. 83, no. 3, pp. 333-349.

14. Butko Ya.A., Grothaus M., Smolyanov O.G. Lagrangian Feynman formulae for second order parabolic equations in bounded and unbounded domains. The Reports of RAS, 2008. vol. 421. no. 6. pp. 727-732.

15. 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.

16. Butko Ya.A., Grothaus M., Smolyanov O.G. Lagrangian Feynman Formulae for Second Order Parabolic Equations in Bounded and Unbounded Domains. IDAQP, 2010, vol. 13, no. 3, pp. 377-392.

17. Gadella M., Smolyanov O.G. Feynman Formulas for Particles with Position-Dependent Mass. Doklady Math., 2007, vol. 77, no. 1, pp. 120-123.

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

19. Sakbaev V.G., Smolyanov O.G. Dynamics of a Quantum Particle with Discontinuous Position-Dependent Mass. Doklady Math., 2010, vol. 82, no. 1, pp. 630-634.

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

21. Smolyanov O.G., Shamarov N.N. Feynman and Feynman-Kac formulae for evolution equations with Vladimirov operator. Doklady Math., 2008, vol. 77, pp. 345-349.

22. Smolyanov O.G., Shamarov N.N. Hamiltonian Feynman Integrals for Equations with the Vladimirov Operator. Doklady Math., 2010, vol. 81, no. 2, pp. 209-214.

23. 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.

24. Feynman R.P. Space-time approach to nonrelativistic quantum mechanics. Rev. Mod. Phys., 1948, vol. 20, pp. 367-387.

25. Nelson E. Feynman integrals and the Schrodinger equation. J. Math. Phys., 1964, vol. 3, pp.332-343.

26. Chernoff P. Product formulas, nonlinear semigroups and addition of unbounded operators. Mem. Am. Math. Soc., 1974, vol. 140.

27. Butko Ya.A., Schilling R.L., Smolyanov O.G. Lagrangian and Hamiltonian Feynman formulae for some Feller semigroups and their perturbations. IDAQP, 2012, pp. 1-19.

28. Reed M., Simon B. Methods of modern mathematical phisycs. Vol. 2. Fourier Analysis, Self-Adjointness. Academic Press, 1975. 361 p. (Russ. ed.: RidM., SaimonB. Metody sovremennoi matematicheskoifiziki. T.2. Garmonicheskii analiz. Samosopriazhennost'. Moscow, Mir, 1977. 393 p.).

29. Pazy A. Semigroups of linear operators and Applications to partial differential equation. Springer-Verlag, 1983. 279 p.

30. Engel K.-J., Nagel R. One-Parameter Semigroups for Linear Evolution Equations. Springer, 2000. 609 p.

31. Zorich V.A. Matematichesky analiz. Ch. 2 [Mathematical analysis. Pt. 2]. Moscow, MCCME, 2002. 789 p.

32. Sveshikov A.G., Tikhonov A.N. Teoriya funkcy compleksnogoperemennogo [Theory of functions of a complex variable]. Moscow, PHYSMATHLIT, 2005. 335 p.

33. Gradshtein I.S., Ruezhik I.M. Tablitsy integralov, summ, riadov iproizvedenii [Integrals, series and sums tables]. Moscow, PHYSMATHLIT, 1963. 1097 p.

34. Butko Ya.A., Smolyanov O.G., Schilling R.L. Feynman formulae for Feller semigroups. Doklady Math., 2010, vol. 434, no. 1, pp. 7-11.

35. 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., 2011, vol. 18, no. 4, pp. 381-399.