CC BY
15
2018, Т. 160, кн. 1 С. 7-16
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)
УДК 535.14
СТОХАСТИЧЕСКАЯ МОДЕЛЬ ДИНАМИКИ ЭЛЕКТРОМАГНИТНЫХ ПОЛЕЙ, СВЯЗАННЫХ ДРУГ С ДРУГОМ НА ЗЕРКАЛЕ
А.М. Башаров
НИЦ «Курчатовский институт», г. Москва, 123182, Россия Московский физико-технический институт, г. Долгопрудный , 141701, Россия
Аннотация
Сформулированы стохастические дифференциальные уравнения и на их основе выведены кинетические уравнения для фотонов микрорезонатора в случае накачки классическим или квантованным сжатыми полями. Показано, что оператор связи полей на зеркале микрорезонатора выводится из оператора нерезонансного взаимодействия полей с квантовой частицей методами алгебраической теории возмущений, при этом возможен учет различных интерференционных процессов.
Ключевые слова: марковское приближение, квантовое стохастическое дифференциальное уравнение, кинетическое уравнение, нерезонансное взаимодействие
Введение
Многие задачи оптической спектроскопии и квантовой оптики описывают процессы со случайными воздействиями окружающей среды на исследуемую систему. В качестве исследуемой системы обычно рассматривают атомы или молекулы, тогда описание их динамики в марковском приближении дается на основе кинетического уравнения в форме Линдблада [1—3]. Эти уравнения легко получаются, если оператор эволюции всей системы и ее окружения представлять стохастическим дифференциальным уравнением (СДУ). Заметим, что если марковское приближение формулировать сразу непосредственно в уравнении эволюции, то иного и быть не может, поскольку тогда уравнения эволюции представляют собой уравнения типа Ланжевена и их математический статус строго определяется переходом к СДУ. Другой путь состоит в переходе от уравнения с оператором эволюции к кинетическому уравнению для матрицы плотности, методов, реализующих такой переход, много (см., например, [4, 5]). Существуют также методы, позволяющие анализировать формальное интегральное представление оператора эволюции, использующее полугрупповое свойство оператора эволюции, и представления его в виде произведения большого числа сомножителей-операторов эволюции на малых подинтервалах разбиения исходного временного промежутка эволюции [6]. Здесь к марковскому приближению, как и в других методах, прибегают на различных этапах вычислений. Однако если при этом ввести квантовые случайные процессы, то их алгебра позволяет легко просуммировать формальный ряд интегрального представления и получить СДУ как обычного винеровского типа, так и обобщенного невинеровского типа [3].
Если в качестве исследуемой системы рассматривать электромагнитное поле в одномодовом резонаторе, то одним из механизмов воздействия окружающей среды на микрорезонатор является связь внешних полей с полем резонатора на
зеркале резонатора, при этом говорят о накачке резонатора. Но можно говорить и о взаимодействии электромагнитных полей на границе раздела различных пространств. В настоящей работе показано, как в случае одномодового резонатора, такое взаимодействие описывается СДУ, причем вывод СДУ представлен в одной и той же форме и для классического, и для квантового случайных полей с учетом сжатия полей. Показаны отличия классической и квантовой накачек и области, в которых квантовую накачку можно моделировать классической. Рассмотрено нерезонансное взаимодействие одномодового и широкополосного полей с квантовой частицей и показано, что последовательное применение алгебраической теории возмущений позволяет рассматривать связь полей на зеркале микрорезонатора как интерференцию нерезонансного взаимодействия полей с квантовой частицей. Это, в свою очередь, указывает на необходимость перенормировки параметров резонатора и дает технику написания слагаемых эффективного гамильтониана, описывающих в подобных задачах различные возможные интерференционные процессы, в том числе и определяющие новые каналы релаксации.
1. СДУ для накачки микрорезонатора
Пусть имеется одномодовый микрорезонатор, на зеркало которого подается классическое электромагнитное поле напряженности Е = Е(г) ехр(-М) + С.с., V - несущая частота, Е(г) - медленно меняющаяся амплитуда. Традиционно вся система описывается гамильтонианом в картине Дирака (представление взаимодействия):
НМ1ггог(г) = кЕ(г) ехр(-М)с+ ехр('шсг) + Н.с.
где с+ и с - операторы рождения и уничтожения фотона микрорезонатора частоты шс, [с, с+] = 1, к - константа связи полей. Буквами Н.с. и С.с. обозначены слагаемые, эрмитово и комплексно сопряженные к предыдущему слагаемому.
Гамильтониан НМ1ГГОГ(г) имеет стандартную простую интерпретацию: фотон из классического поля уничтожается, рождая фотон микрорезонатора, и обратно. Таким образом, идет накачка микрорезонатора.
Уравнением, подлежащим анализу, является уравнение Шредингера для оператора эволюции и(г) или для волнового вектора |Ф(г)> в картине Дирака (считаем все величины безразмерными)
ш (г) = -гНМ1ттот(г)и (г), и (0) = 1, №(г)> = и (г)|Ф(0)>. (1)
dt
Пусть теперь классическое поле - случайное. Для различных случаев естественным требованием, накладываемым на E(t), является стационарность этой величины, то есть зависимость коррелятора (E(t)E* (t')) только от разности времен t -1'.
Если разложить E(t) в интеграл Фурье, то стационарность воздействия влечет за собой некоррелированность (независимость) спектральных компонент, то есть
(Ъ(ш)Ъ*(ш')) = 5(ш - и')S(и), Ъ(ш) = — /exp(iwt)£(t) dt,
2п J
где S(t) - спектр мощности случайного воздействия, Ъ(ш) - фурье-компоненты.
Если принять естественное приближение S(w) = const для определенного спектрального диапазона и распространить его на всю область спектра (приближение Найквиста), то нетрудно заметить, что в соответствующей нормировке
г
{£(г)£*(г')) = пЗ(г — г'). Тогда величина У Е(г') ¿г' = Ш1(г) + Ш2(г), пред-
0
ставленная через действительную и мнимую части, непрерывна, но нигде не дифференцируема (^>о - постоянная фаза). Разными способами можно показать, что плотность функции распределения вероятности р(лу ,г) величин Шу удовлетворяет кинетическому уравнению (уравнению Фоккера-Планка)
др(лу ,г) В у д2р(ту ,г)
дг = ~2 Щ '
решение которого есть
рНМ = — го)Ву ехЧ — 2П,(г — го) )
с начальным условием р(лу,го\лоу,го) = 5(ту — шоу). Величину Ш(г) называют (нестандартным) комплексным винеровским процессом.
Таким образом, уравнение (1) оказывается математически неопределенным в приближении Найквиста. Но можно строго определить его интегральный вариант
и (г) — и (го) = т ехр [—гк(с+(Ш (г) — Ш (го)) + Н.С.)] (2)
при разумном определении интеграла общего вида по винеровскому процессу
I о(г')зш (г').
Если обычным образом определить интеграл У О(г') ¿Ш(г') как предел частич-
п
ных сумм Бп =
)(Ш (гк) — ш (г к —1)), то результат будет зависеть от кон-
к = 1
кретного выбора промежуточных точек гк—1 < Тк < гк на интервалах разбиения. В этом нетрудно убедиться на примере интеграла У Шу (г') ¿Шу (г'), если использовать свойства винеровского процесса Ш(г) - статистическую независимость его приращений Ш(гк) — Ш(гк—1) друг от друга и от Ш(гк—1). Для так называемых неупреждающих функций О(г), статистически независимых в момент времени г от будущего поведения винеровского процесса, удобен и самосогласован выбор Ито Тк = гк—1, при этом предел интегральных сумм понимается как среднеквадратичный предел. Тогда можно ввести согласованные дифференциалы Ито
¿Ш (г)ш *(г) = (В1 + В2^г = пйг,
¿ш (гум *(г) = (В1 + ехр'и^вун = таг, (3)
¿ш(г)ё,г = аtаt = о, \т\ < п, {¿ш(г)) = о.
Здесь дифференциальные соотношения следует понимать как выполнение интегральных равенств типа (для любых неупреждающих функций 0(г')):
г г
J о(г') ¿Ш(г') ¿Ш(г') = ^ о(г') ¿г'.
г0 г0
Через интегральное уравнение (2) определяется СДУ и его решение. Для современной математики характерно сначала определить интеграл по мере, потом производные.
Вывод СДУ для оператора эволюции состоит в определении дифференциала Ито
¿и (г) = и (г + ¿г) - и (г) = (ехр (-гк(с+¿ш (г) + саш *(г))) - \)и (г) = = (-гкс+вш(г) - гксШ*(г))и(г)--7^((сс+ + с+с)п + с2т* + с+2т)и(г). (4)
При выводе использована алгебра (3). В уравнении (4) отсутствует какая-либо отстройка, поскольку в предположении широкополосности спектра внешнего электрического поля его центральная частота будет всегда в резонансе.
Уравнение (4) представляет собой классическое СДУ винеровского типа. Кинетическое уравнение для матрицы плотности получается обычным образом:
вр(г) = р(г + ¿г) -р(г), р(г + ¿г) = |Ф(г + ¿г)){^(г + ¿г)\, вр(г) = ви (г)\Щ0)){Щ0)\и +(г) + и (г)\ъ(0)){ъ(0)\си +(г)+ (5)
¿и (г)\^(0)){Ф(0)\ви +(г).
В результате простых вычислений, после усреднения ¿Ш(г) получаем кинетическое уравнение для матрицы плотности р(г), описывающей фотоны в микрорезонаторе. Форма этого уравнения имеет обобщенную форму Линдблада:
¿¿^ = К-(п(2ср(г)с+ - с+ср(г) - р(г)с+с) + т(2с+р(г)с+ - с+2р(г) - р(г)с+2))+ К
+ —(п(2с+р(г)с - сс+р(г) - р(г)сс+) + т*(2ср(г)с - с2р(г) - р(г)с2)). (6)
По аналогии с классическим винеровским процессом Ш(г) введем квантовый «винеровский» процесс В(г).
Пусть на зеркало подается квантованное электромагнитное поле напряженности Е(г) = J Г(ш)Ьш ехр (-гшг) ¿ш + Н.с. Состоянием поляризации пренебрегаем и
стандартным образом осуществляем переход в суммировании от волнового вектора к частоте [2]. Коммутационные соотношения имеют вид [Ьш ,Ь+, ] = 5(ш - ш'),
К, Ьс] = К ,Ь+] = 0 .
Обычно взаимодействие на зеркале рассматривается сразу в резонансном и марковском приближениях
НЦМотог(г) = к I Ьшс+ ехр(-г(ш - шс)г) ¿ш + Н.с. (7)
Введем квантовые рождающий В+(г) и уничтожающий В(г) процессы
Ь(г) = _ / Ьш ехр(-г(ш - шс)г) ¿ш,
\/2п }
} (8) В(г) = Ь(г') ¿г', ¿В (г) = В(г + ¿г) - В(г),
0
аналогичные тем, которые предшествовали введению классического винеровского процесса:
г
Е(г) = ! Ь(ш)в-*шг ¿ш, ш(г) = ! Е(г') ¿г', ¿ш(г) = ш(г + ¿г) - ш(г).
Интеграл в смысле Ито аналогичен классическому, то есть
N
N'
N
/ I(г) ¿в(г) = Ит VI(гк—1)(В(гк) — в(гк—1).
/ N^00 *—'
^ 1__-I
к=1
Алгебра Ито операторов ¿В(г) и ¿В+(г) дается соотношениями
¿в(г)¿в+(г) = (п + шг, ¿в+(г)¿в (г) = тЛ, в(гыв(г) = m¿г,
(9)
¿в+(г)¿г = ¿в(г)йг = ¿г¿г = о, {¿в(г)) = о.
Важным отличием от классического случая является условие \т\ < \/п(п+1) . Алгебру (9) принято также называть алгеброй Гардинера - Коллет [1].
Квантовое СДУ для дифференциала Ито ¿и(г) следует из формального решения и (г) = Т ехр(—^ НЯМггог(г') ¿г'). Тогда
к2 ¿г
¿и(г) = —iк(c+¿B(г) + сЛв+(г))--— (с+с(п +1) + сс+п + с2т + с+2т*)и(г). (10)
При этом нетрудно убедиться, что выполнено правило дифференцирования Ито (как и в классическом случае)
¿(и (г)и+(г)) = (¿и (г))и +(г) + и (г^и +(г) + (¿и (г))^и + (г)).
Кинетическое уравнение для матрицы плотности получается в результате цепочки преобразований (5) в виде
¿¿^ = у((п+1)(2ср(г)с+ —с+ср(г)—р(г)с+с)+т(2с+р(г)с+ —с+2р(г)—р(г)с+2))+ к2
+ —(п(2с+р(г)с — сс+р(г) — р(г)сс+) + т*(2ср(г)с — с2р(г) — р(г)с2)). (11)
Сравнение случаев накачки классическим шумовым (4), (6) и квантовым (10), (11) полями показывает, что различие возникает только в двух случаях:
1) в случае малого числа фотонов квантованного поля;
2) в случае большого числа фотонов и максимального сжатия \т\ = . Для классического поля второй случай вообще невозможен - чисто формально
он приводит к отрицательности диагональных элементов матрицы плотности, чего быть не должно. Однако с точки зрения моделирования процессов динамики фотонов вдали от областей значения параметров, отвечающих максимальному сжатию и большому числу фотонов, разницы практически нет, так что в этом случае удобно моделировать воздействие квантованного шума классическим винеровским процессом.
2. Модель взаимодействия полей при нерезонансном воздействии на атом
Пусть теперь квантованное широкополосное электромагнитное поле и фотоны одномодового поля нерезонансно взаимодействуют с атомами, локализованными в области пространства с размерами, много меньшими длины волны. Атомы заселяют только нижний энергетический уровень. Исходный оператор взаимодействия
в картине Дирака есть оператор электродипольного взаимодействия, который складывается из оператора взаимодействия Нл-с(г) атомов с локализованной модой и оператора Нл-Р (г) взаимодействия атомов с широкополосным полем:
Нл-с(г) = 9с(с+е^г + св-^г)^2 в^ г\Ек Е, = Ек - Е;, Нл-Р(г) = I ¿шГ(ш) (Ьшв-^г + Ь+в™)^ ¿кзв™»г\Ек)(1){Ез\(€).
г,3к
Собственные вектора г -го атома удовлетворяют обычным условиям
£\Еи )(){Еи\() = 1« {Ек\() Ез)« = 5кз, к
где квантовые невырожденные состояния \Ек)(г) энергии Ек характеризуют г-й атом; через ¿кз обозначены матричные элементы оператора дипольного момента атома. Верхний индекс у векторов состояний отмечает пространство состояний г-го атома, а суммирование по г выполняется по всем атомам обсуждаемого ансамбля.
Следуя унитарной симметрии квантовой теории, совершим унитарное преобразование исходной полной волновой функции [7]
^ (г)) =ехр(^(г))\щг)), Я+(г) = Я(г).
Уравнение Шредингера для преобразованной волновой функции открытой системы и окружения
^ = Н т,«»
определяется преобразованным гамильтонианом
Н (г)=ехр(^(г))(Нл-с(г)+ Нл-Р (г)) ехр^(г)) - г ехр(^(г)) ' ехр^(г)).
Далее в соответствии с понятием открытой квантовой системы оператор Q(г) и преобразованный гамильтониан Н(г) разлагаются в ряды теории возмущений по константам взаимодействия подсистем открытой квантовой системы с внешними полями и между собой, порядок слагаемых которых отмечается верхними индексами [8]. Для определенности считаем, что первый (левый) индекс в группе верхних индексов слагаемых указывает на порядок по параметру связи атомов с широкополосным полем окружения, второй - на порядок по параметру связи атома с резонаторной модой:
Q(г) = Q(0'1) + Q(1'0) + ...,
(12)
н (г) = н (0'1) (г) + н(1'0)(г) + н (1'1)(г) + н (2'0)(г) + ■■■.
Критерий [7, 8] отличается от известных, но формально приводит к тому же виду слагаемых первого порядка по константе связи с окружением. Будем считать, что в слагаемых (12) отсутствуют быстропеременные временные множители, или, что то же самое, все слагаемые (12) медленно меняются во времени. В случае широкополосного квантованного (или классического) поля в широкополосном поле всегда найдутся частоты, которые будут резонансными некоторым атомным переходам с нижнего энергетического уровня на возбужденные. Если в квантованном широкополосном поле плотность числа фотонов нулевая, то это не приведет к реальным
переходам в атоме. Но такие переходы разбивают широкополосное поле на совокупность независимых шумовых источников. Чтобы не включать в рассмотрение побочные процессы, будем рассматривать только один из шумовых источников и в дальнейшем о нем будем говорить как о широкополосном электромагнитном поле. Этот источник будет нерезонансен всем атомным переходам, но, чтобы он проявился в задаче нерезонансного взаимодействия, его центральная частота должна совпадать с частотой локализованной фотонной моды. Требование отсутствия быстропеременных временных множителей приводит к уравнениям
i,kj
н (10)(t) = Hi (01 = o
0kj + —c 0kj — —c
Q(0,i'4t) = )<«j'( c+j± + ¿j),
/—J V 0kj + —c 0kj — —c J
e„ (Qkj +u)t e-i(Q.kj-u)t
d— dkj \Ek)() Ej \(i)(b+ 1Г~— + К 0 „
f-r V 0lkj + — Okj — —
i kj
где знак штрих у суммы по атомным энергетическим уровням означает отсутствие слагаемых с резонансными знаменателями с 0kj ± — ~ 0.
По формулам
н(2,0)(t) = — 2 [Q(1 , 0)(t),HA-F (t)]' — 2 [Q(1 , 0)(t),H(1,0)(t)]', H(0,2)(t) = — 2 [Q(0 •1)(t),HA-c(t)]' — 2 [Q(0 '1)(t),H(0 •1)(t)]',
Н (1'1)(t) = — 2 [Q(1'0)(t),HA-c(t)]> — 2 [Q(0'1)(t),HA-F (t)]'—
— 2 [Q(1'0)(t),H (0'1)(t)]' — 2 [Q(0'1)(t),H (10)(t)]',
где штрих отмечает отсутствие быстро меняющихся во времени слагаемых, нетрудно найти слагаемые эффективного гамильтониана второго порядка по константам взаимодействия с электромагнитными полями. Именно требование медленного изменения во времени слагаемых
H(2'0)(t) и H(1'1)(t) выделяет из спектра широкополосного электромагнитного поля учтенный независимый шумовой источник с центральной частотой —c. В отсутствие каких-либо комбинационных резонансов полученные формулы содержат поправки лэмбовского типа к частотам атомных переходов и к частоте локализованной фотонной моды, а также следующий оператор (Щ (—) - параметр теории [3, 7])
H(11\t) = J'd— V(—)gcbu c+e-i(u-uc)^ nk—)J2 \Ek)(i) Ek\(i) + H.c..
ki
Если в полученном выражении сделать марковское приближение, то есть потребовать выполнения нижеприведенных условий, то получим уравнения (10), (11), в которых коэффициент связи полей будет выражен через параметры нерезонансных атомов среды зеркала. Марковские условия для уравнения эволюции с оператором Н(11)(t) = HNonresonant имеют вид
Г— = const, (Ьш b+,) = (n + 1)£(— — —'), (b+ Ьш) = nS(— — —'), (bu Ьш>) = mS(2—c — — — —')■
Заключение
Оператор связи полей на зеркале обычно записывается сразу в резонансном и марковском приближениях. Представленный нами подход позволяет рассмотреть оператор связи на зеркале без подобных приближений. Тогда при помощи алгебраической теории возмущений определяются различные интерференционные слагаемые, описанные в работах [7, 8], в том числе и поправки к сдвигам частот. В настоящей статье было показано, что именно таким интерференционным слагаемым в случае нерезонансного взаимодействия с атомным ансамблем и является общепринятая связь полей на зеркале.
Полученные результаты необходимы для формулировки теории для описания динамики пространственно распределенных систем, теория которых на основе СДУ невинеровского типа строго сформулирована ранее лишь в случае, когда открытая система локализована в области пространства с размерами, много меньшими длины волны [3, 7]. Если заменить операторы рождения и уничтожения фотонов одномодового поля на операторы рождения и уничтожения второго квантованного широкополосного поля или компоненты Фурье классического случайного поля (с учетом аналогий разд. 1), то получим взаимосвязь широкополосных полей. Модели, в которых одно из таких полей представляет собой классический белый шум, а другое поле - квантовый шум, нам неизвестны. Некоторые построения обобщенной алгебры для случая двух различных широкополосных квантованных полей с нулевой плотностью фотонов рассмотрены в [9].
Благодарности. Автор выражает благодарность профессору В.В. Самарцеву за проявленный интерес и внимание к настоящей работе.
Работа выполнена при частичной финансовой поддержке РФФИ (проект № 16-02-00453а).
Литература
1. Gardiner C.W., Collet M.J. Input and output in damped quantum systems: Quantum stochastic differential equations and the master equation // Phys. Rev. A. - 1985. - V. 31, No 6. - P. 3761-3774. - doi: 10.1103/PhysRevA.31.3761.
2. Gardiner C.W., Zoller P. Quantum noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics. - Berlin: SpringerVerlag, 2000. - 438 p.
3. Basharov A.M. Stark interaction of identical particles with the vacuum electromagnetic field as quantum Poisson process suppressing collective spontaneous emission // Phys. Rev. A. - 2011. - V. 84, No 1. - Art. 013801, P. 013801-1-013801-14. - doi: 10.1103/PhysRevA.84.013801.
4. Bonitz M. Quantum Kinetic Theory. - Stuttgart: B.G.Teubner, 1989.
5. Кузнецов Д.В., Рерих Вл.К., Гладуш М.Г. Локальное поле и скорость радиационной релаксации в диэлектрической среде // Журн. эксперим. и теор. физики. - 2011. -Т. 140, Вып. 4. - С. 742-754.
6. Pichler H., Zoller P. Photonic circuits with time delays and quantum feedback // Phys. Rev. Lett. - 2016. - V. 116, No 9. - Art. 093601, P. 093601-1-093601-6. - doi: 10.1103/PhysRevLett.116.093601.
7. Maimistov A.I., Basharov A.M. Nonlinear optical waves. - Netherlands: Springer, 1999.-XIV, 656 p. - doi: 10.1007/978-94-017-2448-7.
8. Башаров А.М. Квантовая теория открытых систем на основе стохастических дифференциальных уравнений обобщенного ланжевеновского (невинеровского) типа // Журн. эксперим. и теор. физики. - 2012. - Т. 142, Вып. 3 - С. 419-441.
9. Башаров А.М. Теория открытых систем на основе стохастических дифференциальных уравнений // Оптика и спектроскопия. - 2014. - Т. 116, № 4. - С. 532-540. - doi: 10.7868/S0030403414040059.
Поступила в редакцию 04.12.17
Башаров Асхат Масхудович, кандидат физико-математических наук, старший научный сотрудник; доцент кафедры математики и математических методов физики НИЦ «Курчатовский институт»
пл. Академика Курчатова, д. 1, г. Москва, 123182, Россия Московский физико-технический институт
Институтский пер., д. 9, г. Долгопрудный, 141701, Россия E-mail: basharov@gmail.com
ISSN 2541-7746 (Print)
ISSN 2500-2198 (Online)
UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI
(Proceedings of Kazan University. Physics and Mathematics Series)
2018, vol. 160, no. 1, pp. 7-16
A Stochastic Model for the Dynamics of Electromagnetic Waves Connected to Each Other on the Mirror
A.M. Basharov
aMoscow Aviation Institute (National Research University), Moscow, 125993 Russia bInstitute of Applied Mechanics, Russian Academy of Sciences, Moscow, 125040 Russia
E-mail: basharov@gmail.com
Received December 4, 2017 Abstract
The paper develops stochastic differential equations serving as the basis for derivation of master equations for microcavity photons pumped with either classic or quantized contracted fields. The interaction operator on the microcavity mirror is derived from the non-resonant field operator interacting with the quantum particle by means of the algebraic perturbation theory, with allowance for various interference processes.
Keywords: Markovian approximation, quantum stochastic differential equation, master equation, nonresonant interaction
Acknowledgments. We thank Professor V.V. Samartsev for his interest in and support for this research.
The study was supported in part by the Russian Foundation for Basic Research (project no. 16-02-00453a).
References
1. Gardiner C.W., Collet M.J.Input and output in damped quantum systems: Quantum stochastic differential equations and the master equation. Phys. Rev. A., 1985, vol. 31, no. 6, pp. 3761-3774. doi: 10.1103/PhysRevA.31.3761.
2. Gardiner C.W., Zoller P. Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics. Berlin, SpringerVerlag, 2000. 438 p.
3. Basharov A.M. Stark interaction of identical particles with the vacuum electromagnetic field as quantum Poisson process suppressing collective spontaneous emission. Phys. Rev. A., 2011, vol. 84, no. 1, art. 013801, pp. 013801-1-013801-14. doi: 10.1103/Phys-RevA.84.013801.
4. Bonitz M. Quantum Kinetic Theory. Stuttgart, B.G.Teubner, 1989.
5. Kuznetsov D.V., Roerich Vl.K., Gladush M.G. Local field and radiative relaxation rate in a dielectric medium. J. Exp. Theor. Phys., 2011, vol. 113, no. 4, pp. 647-658. doi: 10.1134/S1063776111100050.
6. Pichler H., Zoller P. Photonic circuits with time delays and quantum feedback. Phys. Rev. Lett., 2016, vol. 116, no. 9, art. 093601, pp. 093601-1-093601-6. doi: 10.1103/Phys-RevLett.116.093601.
7. Maimistov A.I., Basharov A.M. Nonlinear Optical Waves. Netherlands, Springer, 1999. XIV, 656 p. doi: 10.1007/978-94-017-2448-7.
8. Basharov A.M. Quantum theory of open systems based on stochastic differential equations of generalized Langevin (non-Wiener) type. J. Exp. Theor. Phys., 2012, vol. 115, no. 3, pp. 371-391. doi: 10.1134/S1063776112070035.
9. Basharov A.M. A theory of open systems based on stochastic differential equations. Opt. Spectrosc., 2014, vol. 116, no. 4, pp. 495-503. doi: 10.1134/S0030400X14040055.
I Для цитирования: Башаров А.М. Стохастическая модель динамики электромаг-( нитных полей, связанных друг с другом на зеркале // Учен. зап. Казан. ун-та. Сер. \ Физ.-матем. науки. - 2018. - Т. 160, кн. 1. - С. 7-16.
For citation : Basharov A.M. A stochastic model for the dynamics of electromagnetic waves connected to each other on the mirror. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2017, vol. 160, no. 1, pp. 7-16. (In Russian)
