Министерство науки и высшего образования Российской Федерации
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ ИТМО
ITMO University
ВЫПУСКНАЯ КВАЛИФИКАЦИОННАЯ РАБОТА/GRADUATION THESIS
Динамика полевых наблюдаемых в квантовых каналах с учетом термализации
Автор/ Author
Медведева Светлана Сергеевна
Направленность (профиль) образовательной программы/Major
Оптические и квантовые технологии в коммуникациях 2017
Квалификация/ Degree level
Бакалавр
Руководитель ВКР/ Thesis supervisor
Гайдаш Андрей Алексеевич, кандидат физико-математических наук, Университет ИТМО,
лаборатория квантовых процессов и измерений, ведущий научный сотрудник
Группа/Group
V3400
Факультет/институт/кластер/ Faculty/Institute/Cluster
факультет фотоники
Направление подготовки/ Subject area
12.03.03 Фотоника и оптоинформатика
Обучающийся/Student
Документ
подписан
Медведева
Светлана
Сергеевна
14.05.2021
(эл. подпись/ signature)
Руководитель ВКР/ Head
of Graduate Project
Документ
подписан
Гайдаш Андрей
Алексеевич
14.05.2021
(эл. подпись/ signature)
Медведева
Светлана
Сергеевна
(Фамилия И.О./ name
and surname)
Гайдаш Андрей
Алексеевич
(Фамилия И.О./ name
and surname)
Министерство науки и высшего образования Российской Федерации
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ ИТМО
ITMO University
ЗАДАНИЕ НА ВЫПУСКНУЮ КВАЛИФИКАЦИОННУЮ РАБОТУ /
OBJECTIVES FOR A GRADUATION THESIS
Обучающийся / Student Медведева Светлана Сергеевна
Группа/Group V3400
Факультет/институт/кластер/
Faculty/Institute/Cluster
факультет
фотоники
и
оптоинформатики
Квалификация/ Degree level Бакалавр
Направление подготовки/ Subject area 12.03.03 Фотоника и оптоинформатика
Направленность (профиль) образовательной программы/Major Оптические и
квантовые технологии в коммуникациях 2017
Специализация/ Specialization
Тема ВКР/ Thesis topic Динамика полевых наблюдаемых в квантовых каналах с учетом
термализации
Руководитель ВКР/ Thesis supervisor Гайдаш Андрей Алексеевич, кандидат физикоматематических наук, Университет ИТМО, лаборатория квантовых процессов и измерений,
ведущий научный сотрудник
Срок сдачи студентом законченной работы до / Deadline for submission of complete
thesis 31.05.2021
Техническое задание и исходные данные к работе/ Requirements and premise for the
thesis
Целью выполнения работы является исследование свойств оптоволоконных квантовых
каналов, в частности динамики полевых наблюдаемых с учетом термализации. Задачи,
решаемые при выполнении работы: 1. Вывод и анализ выражений, описывающих
динамику среднего числа фотонов и его дисперсии для сжатых вакуумных состояний. 2.
Вывод и анализ выражений, описывающих динамику среднего числа фотонов и его
дисперсии для сжатых когерентных состояний. 3. Вывод и анализ выражений,
описывающих динамику среднего числа фотонов и его дисперсии для суперпозиции
когерентных состояний. 4. Вывод и анализ выражений, описывающих динамику средних
значений операторов квадратур в случае сжатых вакуумных состояний. 5. Вывод и анализ
выражений, описывающих динамику средних значений операторов квадратур в случае
сжатых когерентных состояний. 6. Вывод и анализ выражений, описывающих динамику
средних значений операторов квадратур в случае суперпозиции когерентных состояний.
Содержание выпускной квалификационной работы (перечень подлежащих
разработке вопросов)/ Content of the thesis (list of key issues)
1. Исследование динамики открытых квантовых систем.
2. Исследование динамики среднего числа фотонов и его дисперсии в оптическом волокне.
3. Исследование динамики средних значений операторов квадратур в оптическом волокне.
Перечень графического материала (с указанием обязательного материала) / List of
graphic materials (with a list of required material)
Исходные материалы и пособия / Source materials and publications
Pirandola S. et al. Advances in quantum cryptography //Advances in Optics and Photonics. –
2020. – Т. 12. – №. 4. – С. 1012-1236.
Rivas A., Huelga S. F. Open quantum systems. – Berlin : Springer, 2012. – Т. 10. – С. 978-3.
Schoeller H. Dynamics of open quantum systems, Computing Solids: Models, Ab-initio Methods
and Supercomputing, Lecture Notes of the 45th Spring School 2014. – 2014.
Zhang W. M. Exact master equation and general non-Markovian dynamics in open quantum
systems //The European Physical Journal Special Topics. – 2019. – Т. 227. – №. 15. – С. 18491867.
Megier N., Smirne A., Vacchini B. Evolution equations for quantum semi-Markov dynamics
//Entropy. – 2020. – Т. 22. – №. 7. – С. 796.
Дата выдачи задания/ Objectives issued on 09.03.2021
СОГЛАСОВАНО / AGREED:
Руководитель ВКР/
Thesis supervisor
Документ
подписан
Гайдаш Андрей
Алексеевич
09.03.2021
Гайдаш Андрей
Алексеевич
(эл. подпись)
Задание принял к
исполнению/ Objectives
assumed by
Документ
подписан
Медведева
Светлана
Сергеевна
19.03.2021
Медведева
Светлана
Сергеевна
(эл. подпись)
Руководитель ОП/ Head
of educational program
Документ
подписан
Козлов Сергей
Аркадьевич
29.04.2021
(эл. подпись)
Козлов Сергей
Аркадьевич
Министерство науки и высшего образования Российской Федерации
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ ИТМО
ITMO University
АННОТАЦИЯ
ВЫПУСКНОЙ КВАЛИФИКАЦИОННОЙ РАБОТЫ /
SUMMARY OF A GRADUATION THESIS
Обучающийся/ Student
Медведева Светлана Сергеевна
Наименование темы ВКР / Title of the thesis
Динамика полевых наблюдаемых в квантовых каналах с учетом термализации
Наименование организации, где выполнена ВКР/ Name of organization
Университет ИТМО
ХАРАКТЕРИСТИКА ВЫПУСКНОЙ КВАЛИФИКАЦИОННОЙ РАБОТЫ/
DESCRIPTION OF THE GRADUATION THESIS
1. Цель исследования / Research objective
Исследование свойств оптоволоконных квантовых каналов, в частности динамики полевых
наблюдаемых с учетом термализации.
2. Задачи, решаемые в ВКР / Research tasks
1. Вывод и анализ выражений, описывающих динамику среднего числа фотонов и его
дисперсии для сжатых вакуумных состояний. 2. Вывод и анализ выражений, описывающих
динамику среднего числа фотонов и его дисперсии для сжатых когерентных состояний. 3.
Вывод и анализ выражений, описывающих динамику среднего числа фотонов и его
дисперсии для суперпозиции когерентных состояний. 4. Вывод и анализ выражений,
описывающих динамику средних значений операторов квадратур в случае сжатых
вакуумных состояний. 5.Вывод и анализ выражений, описывающих динамику средних
значений операторов квадратур в случае сжатых когерентных состояний. 6. Вывод и анализ
выражений, описывающих динамику средних значений операторов квадратур в случае
суперпозиции когерентных состояний.
3. Краткая характеристика полученных результатов / Short summary of
results/conclusions
1. Были выведены и проанализированы выражения, описывающие динамику дисперсии
среднего числа фотонов сжатого вакуумного состояния, сжатого когерентного состояния и
суперпозиции когерентных состояний при прохождении через оптоволокно. 2. Были
определены расстояния распространения, на которых происходит потеря квантовых
свойств состояния, и расстояния, на которых можно считать произошедшей термализацию
системы. 3. Были выведены и проанализированы выражения, описывающие динамику
средних значений наблюдаемых операторов квадратур для сжатого когерентного состояния
и суперпозиции когерентных состояний при прохождении через оптоволокно. 4. Была
получена зависимость уменьшения амплитуды средних значений наблюдаемых от
расстояния распространения.
4. Наличие публикаций по теме выпускной работы/ Have you produced any
publications on the topic of the thesis
1
2
Gaidash A., Kozubov A., Medvedeva S., Miroshnichenko G. The Influence of Signal
Polarization on Quantum Bit Error Rate for Subcarrier Wave Quantum Key Distribution
Protocol//Entropy, 2020, Vol. 22, No. 12, pp. 1393 (Статья; Scopus, Web of Science)
Медведева С.С., Гайдаш А.А., Мирошниченко Г.П. Анализ динамики дисперсии
среднего числа фотонов произвольного квантового состояния в оптическом волокне//
Сборник трудов XII Международной конференции «Фундаментальные проблемы
оптики – 2020» - 2020. - С. 95 (Тезисы)
5. Наличие выступлений на конференциях по теме выпускной работы/ Have you
produced any conference reports on the topic of the thesis
1
2
КМУ 2021, 14.04.2021 - 17.04.2021 (Конференция, статус - всероссийский)
XII Международная конференция "Фундаментальные проблемы оптики" (ФПО 2020),
19.10.2020 - 23.10.2020 (Конференция, статус - международный)
6. Полученные гранты, при выполнении работы/ Grants received while working on the
thesis
7. Дополнительные сведения/ Additional information
Обучающийся/Student
Документ
подписан
Медведева
Светлана
Сергеевна
14.05.2021
(эл. подпись/ signature)
Руководитель ВКР/ Head
of Graduate Project
Документ
подписан
Гайдаш Андрей
Алексеевич
14.05.2021
(эл. подпись/ signature)
Медведева
Светлана
Сергеевна
(Фамилия И.О./ name
and surname)
Гайдаш Андрей
Алексеевич
(Фамилия И.О./ name
and surname)
СОДЕРЖАНИЕ
ВВЕДЕНИЕ .............................................................................................................. 7
1
2
ДИНАМИКА ОТКРЫТЫХ КВАНТОВЫХ СИСТЕМ ................................. 9
1.1
Открытая квантовая система ..................................................................... 9
1.2
Методы описания динамики открытых квантовых систем .................. 10
1.3
Уравнение Лиувилля – фон Неймана ..................................................... 11
1.4
Уравнение Накадзимы – Цванцига ......................................................... 12
1.5
Уравнение Горини – Коссаковского – Сударшана – Линдблада ........ 13
1.6
Динамика одномодового состояния в оптоволокне .............................. 14
ДИНАМИКА ДИСПЕРСИИ СРЕДНЕГО ЧИСЛА ФОТОНОВ ................ 17
2.1
Уравнение динамики дисперсии среднего числа фотонов состояния 17
2.2 Динамика дисперсии среднего числа фотонов сжатого вакуумного
состояния............................................................................................................. 18
2.3 Динамика дисперсии среднего числа фотонов сжатого когерентного
состояния............................................................................................................. 22
2.4 Динамика дисперсии среднего числа фотонов суперпозиции
когерентных состояний ..................................................................................... 29
2.5 Динамика дисперсии среднего числа фотонов состояния: основные
выводы................................................................................................................. 33
3
ДИНАМИКА НАБЛЮДАЕМЫХ ОПЕРАТОРОВ КВАДРАТУР ............. 35
3.1 Уравнение динамики средних значений наблюдаемых операторов
квадратур состояния .......................................................................................... 35
3.2 Динамика средних значений наблюдаемых операторов квадратур
сжатого когерентного состояния ...................................................................... 36
3.3 Динамика средних значений наблюдаемых операторов квадратур
суперпозиции когерентных состояний ............................................................ 39
3.4 Динамика среднего значения наблюдаемых операторов квадратур
состояния: основные выводы ............................................................................ 42
ЗАКЛЮЧЕНИЕ ..................................................................................................... 43
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ ........................................... 45
6
ВВЕДЕНИЕ
В настоящее время квантовые технологии являются одним из наиболее
активно развивающихся направлений исследований [1]. Можно выделить
четыре области исследований: квантовые коммуникации, опирающиеся на
постулаты квантовой физики для защищённой передачи данных [2]; квантовая
симуляция, где доступная для контроля квантовая система используется для
моделирования системы, управление характеристиками которой не доступно
в лаборатории [3]; квантовые вычисления, использующие уникальные
квантовые эффекты для значительного ускорения вычислений, например,
факторизации чисел [4]; квантовая метрология, где высокочувствительная к
внешнему воздействию квантовая система помогает повысить точность
производимых измерений физических величин [5].
Такой
широкий
спектр
применений
обусловлен
возможностью
использования не имеющих аналогов в классической физике квантовых
эффектов. Явления суперпозиции и запутанности выступают уникальным
преимуществом в решении ряда задач. Однако квантовые состояния
подвержены декогеренции – большинство квантовых неклассических
особенностей могут быть разрушены воздействием окружающей среды на
квантовую систему [6]. При описании изменения во времени квантовой
системы необходимо учитывать явление декогеренции, что можно корректно
сделать в рамках теории открытых квантовых систем. Поэтому особый
интерес не только в сфере фундаментальных исследований, но и в сфере
непосредственного
применения
квантовых
технологий
представляет
рассмотрение открытых квантовых систем.
Так в данной работе с помощью теории открытых квантовых систем
рассматриваются особенности распространения квантового состояния света в
оптоволоконном
канале.
А
именно
производится
вывод
и
анализ
аналитических выражений, описывающих динамику полевых наблюдаемых:
дисперсии среднего числа фотонов и средних значений операторов квадратур.
7
Проводится исследование важной энергетической характеристики квантового
состояния – дисперсии среднего числа фотонов в моде, так как наблюдение
данной величины позволяет дать оценку времени термализации квантового
состояния в канале. Кроме того, в работе проводится исследование изменения
во времени состояния с помощью анализа динамики средних значений
наблюдаемых
операторов
квадратур.
Представлены
результаты
моделирования найденных динамик для нескольких конкретных примеров
квантовых
состояний,
наиболее
активно
используемых
как
в
фундаментальных исследованиях, так и в технических решениях: сжатого
вакуумного, сжатого когерентного и суперпозиции когерентных состояний.
8
1 ДИНАМИКА ОТКРЫТЫХ КВАНТОВЫХ СИСТЕМ
1.1 Открытая квантовая система
Для
описания
квантовомеханических
систем
с
учетом
их
взаимодействия с внешней средой используется теория открытых квантовых
систем. Открытые квантовые системы – квантовые системы, которые
обмениваются энергией или частицами с окружающей средой. Открытые
квантовые системы характеризуются неунитарной динамикой, тогда как
изолированные квантовые системы характеризуются унитарной динамикой
[7].
Состояние квантовой системы S в Гильбертовом пространстве
S
может быть описано с помощью матрицы плотности ( t ) . Тогда изменение
энергии этой системы может быть представлено временной эволюцией
. Окружающая среда
наблюдаемой гамильтониана H ( t ) = Hˆ ( t ) = tr Hˆ ( t )
B (от heat bath, термостат: система с бесконечным числом степеней своды, так
что частоты её мод образуют континуум, находящаяся в состоянии
термического равновесия), с которой входит во взаимодействие открытая
система S , в своем Гильбертовом пространстве
временной эволюцией матрицы плотности
B
может быть описана
B ( t ) . Тогда наблюдаемая
гамильтониана H B ( t ) = Hˆ B ( t ) = tr Hˆ B B ( t ) описывает изменение энергии.
Можно перейти к рассмотрению полной изолированной системы S + B , в
расширенном Гильбертовом пространстве
S
B
описываемой матрицей
плотности SB ( t ) . В этом случае интересующая квантовая система может быть
представлена приведенной матрицей плотности ( t ) = trB SB ( t ) . При этом
изменение энергии полной системы описывает наблюдаемая гамильтониана
H SB ( t ) = Hˆ SB ( t ) = tr Hˆ SB SB ( t ) , где оператор учитывает взаимодействие
открытой системы с окружающей средой:
9
Hˆ SB = Hˆ IˆB + Iˆ Hˆ B + Hˆ I ,
(1)
где гамильтониан Hˆ I ( t ) отвечает за взаимодействие системы с окружающей
средой.
1.2 Методы описания динамики открытых квантовых систем
Динамика открытых квантовых систем может быть описана различными
способами. Отличия существующих подходов заключаются в ряде критериев.
Например, по представлению системы, рассматриваемой в задаче, подходы
можно разделить на представляющие квантовую систему непосредственно как
открытую и на производящие замену открытой системы рассмотрением
большей, но формально изолированной системы. Может различаться
математический
представления
корреляционные
объект,
матриц
динамика
плотности
функции.
которого
исследуется:
используются
Объединяет
различные
волновые
подходы
кроме
функции,
наличие
некоторых приближений, учитывающих важные для конкретной задачи
условия и отбрасывающих те, которыми можно пренебречь [8].
Одно из актуальных направлений описания динамики открытых
квантовых систем состоит в рассмотрении полной волновой функции
многочастичной изолированной системы. Для представления в явном виде
взаимодействия квантовой системы с окружающей средой в общем случае
требуется создание моделей динамики квантовой системы со значительным
числом степеней свободы, что требует высоких вычислительных мощностей.
Поэтому большинство этих методов не являются аналитическими, а
позволяют получить численное моделирование в рамках конкретных задач.
Примером могут служить: усовершенствованные подходы, основанные на
методе Хартри-Фока [9, 10], рассмотрение усеченного Гильбертово
пространства [11], анзац Давыдова [12]. Эти методы особенно популярны в
сфере исследований физической химии.
10
Другим
подходом
может
быть
рассмотрение
различных
корреляционных функций. Примером может выступать использование
функции Грина для частиц-фермионов в случае рассмотрения задач квантовой
термодинамики [13].
Методом, наиболее часто используемым в сфере квантовой оптики,
является представление приведенной матрицы плотности. В основе этого
подхода стоит исключение степеней свободы окружающей среды за счёт
перехода
от
рассмотрения
матрицы
плотности
полной
системы
к
исследованию динамики подсистемы (являющейся интересующей открытой
системой). Тогда временная эволюция открытой квантовой системы
описывается основным кинетическим уравнением [14].
1.3 Уравнение Лиувилля – фон Неймана
Полная система S + B , состоящая из системы
S , рассмотрение которой
интересует, и окружающей среды B являет собой изолированную квантовую
систему SB ( t ) , эволюция которой описывается уравнением Лиувилля – фон
Неймана (здесь и далее постоянная Планка
= 1):
SB ( t ) = −i Hˆ SB , SB ( t ) .
t
Стоит отметить, что найденная матрица плотности
(2)
SB ( t ) содержит
предсказание динамики всех возможных наблюдаемых. После нахождения
решения уравнения может быть произведен переход к системе
S
( t ) = trB SB ( t ). Однако в общем случае получение решения уравнения (2) с
учетом особенности гамильтониана (1) является трудоемкой задачей. И так как
целью выступает рассмотрение динамики именно системы S , удобнее
произвести переход к уравнению, в явном виде описывающему динамику
системы S . В таком уравнении – основном кинетическом – ( t ) выступает
величиной, эволюция которого описывается в явном виде, в то время как
11
характеристики окружающей среды B становятся параметрами. Переход к
основному кинетическому уравнению может быть произведен с помощью
ряда приближений, вводимых в рассмотрение уравнения Лиувилля – фон
Неймана.
1.4 Уравнение Накадзимы – Цванцига
В рамках статистической физики существует подход к разбиению
динамики сложной системы на релевантную и нерелевантную части с
помощью операторов проекции – формализм Цванцига – Мори [15]. Этот
метод в рамках механики текучих сред и физики конденсированного
состояния позволяет найти наиболее точные решения, описывающие
динамику релевантной части.
При рассмотрении квантовомеханических систем основные идеи
формализма
Цванцига-Мори
используется
для
выведения
уравнения
Накадзимы – Цванцига [16, 17]. Вводимый оператор проекции ˆ разбивает
(
)
матрицу плотности полной системы SB ( t ) = ˆ + ˆ SB ( t ) на релевантную
ˆ
SB
(t )
и
нерелевантную
ˆ
SB
(t )
части,
где
(
)
ˆ Iˆ − ˆ .
Далее
рассматривается уравнение Лиувилля – фон Неймана для части ˆ SB ( t ) .
Найденное точное интегро-дифференциальное уравнение, описывающее
эволюцию ˆ SB ( t ) , – уравнение Накадзимы – Цванцига – в представлении
взаимодействия [7] имеет следующий вид:
ˆ
SB ( t ) = ˆ ˆ ( t ) ˆ ( t0 ) ˆ SB ( t0 ) + ˆ ˆ ( t ) ˆ SB ( t ) +
t
t
+ ds ˆ ( t , s ) SB ( s ),
(3)
2
t0
где – параметр взаимодействия при рассмотрении гамильтониана полной
системы вида: Hˆ SB ( t ) = Hˆ SB0 ( t ) + Hˆ I ( t ) , где Hˆ SB0 ( t ) отвечает за несвязанные
12
систему и среду; ˆ ( t ) – лиувиллиан системы: ˆ ( t ) −i Hˆ I ( t ) ,
; ˆ (t, s ) –
t
введенный пропагатор: ˆ ( t , s ) T exp ds ' ˆ ˆ ( s ') , где T – оператор
s
хронологического порядка; ˆ ( t , s ) – ядро: ˆ ( t , s ) = ˆ ˆ ( t ) ˆ ( t , s ) ˆ ˆ ( s ) ˆ .
Однако поиск решения уравнения Накадзимы – Цванцига выступает
настолько же сложной задачей, как в случае рассмотрения уравнения
Лиувилля – фон Неймана для полной системы.
Поэтому часто вводятся
некоторые приближения, позволяющие упростить вид уравнения (3).
Например, состояние нерелевантной системы в начальный момент времени
принимается нулевым ˆ SB ( t0 ) = 0 , изменятся ядро ˆ ( t , s ) или производится
разложение по малому параметру [7]. В месте с этим исследуются и другие
возможности внесения приближений в уравнение Накадзимы – Цванцига [18].
1.5 Уравнение Горини – Коссаковского – Сударшана – Линдблада
Другой подход к поиску решения ( t ) состоит в непосредственном
применении последовательных приближений к уравнению Лиувилля – фон
Неймана для полной системы. Каждое последующее приближение порождает
подгруппу уравнений нового типа [19].
Первое из приближений, которые возможно применить к уравнению
Лиувилля – фон Неймана, состоит в предположении, что взаимозависимость
системы и окружающей среды имеет медленный темп развития. Это
рассмотрение
систем,
характеризующихся
малой
величиной
силы
взаимодействия со средой, носящее название приближения Борна, позволяет
отбросить члены по степени параметра взаимодействия старше квадратичных.
Следующее приближение – приближение Маркова – может быть
сформулировано как отсутствие памяти в рассматриваемой системе.
Формально происходит подынтегральная замена матрицы плотности,
зависящей от параметра интегрирования (интегрирование происходит по
13
времени) на матрицу плотности, зависящую только от времени в момент
наблюдения. Примером уравнения такого типа может служить уравнение
Редфилда [20, 21]. Система, удовлетворяющая этому приближению, носит
название марковской системы. Стоит отметить, что в настоящее время интерес
представляют также и не Марковские системы [22, 23].
Третье, секулярное приближение отбрасывает члены, колебания
которых происходят много быстрее времени эволюции интересующей
Марковской системы. Это приближение позволяет получить уравнение
Горини – Коссаковского – Сударшана – Линдблада (уравнение Линдблада)
[24, 25], которое может быть представлено в следующем виде [26]:
1
( t ) = −i Hˆ + Hˆ LS , ( t ) + Lˆi ( t ) Lˆi † − Lˆi † Lˆi , ( t )
t
2
i
,
(4)
где Hˆ LS – гамильтониан, введенный для перенормировки системы; Lˆi , Lˆi † –
операторы квантового скачка.
В данное время многие фундаментальные теоретические работы [27–29]
посвящены рассмотрению именно уравнения Линдблада. Однако прикладные
вопросы так же требуют анализа. Например, далее будет приведен вид
уравнения типа Линдблада, использующийся применительно к задаче
описания электромагнитного излучения в оптоволоконном канале.
1.6 Динамика одномодового состояния в оптоволокне
Конкретный вид гамильтонианов, описывающих рассматриваемую
систему, определяет и вид уравнения Линдблада. Тогда, как показано в [14],
для случая одномодового состояния монохроматического электромагнитного
поля с частотой в оптоволоконном канале возможно рассмотрение
следующего основного кинетического уравнения:
14
ˆ
ˆ
t ( t ) = −i H , ( t ) + Г ( t )
( t ) = 0
t =0
(5)
для гамильтониана системы, имеющего вид:
1
Hˆ = aˆ †aˆ +
2
(6)
где aˆ † , aˆ – оператор рождения и оператор уничтожения; с использованием
следующего оператора:
Г̂ ( t ) = −
†
†
†
( nT + 1) ( aˆ aˆ ( t ) + ( t ) aˆ aˆ − 2aˆ ( t ) aˆ )
2 + nT ( aa
ˆ ˆ † ( t ) + ( t ) aa
ˆ ˆ † − 2aˆ † ( t ) aˆ )
(7)
где – скорость термализации [30]:
e
где
− t
L
= 10 10 ,
(8)
– коэффициент, определяющий потери оптоволокна; L =
ct
–
расстояние, пройденное излучением в оптоволокне, где c – скорость света,
– константа распространения;
nT – среднее число тепловых фотонов:
nT =
1
e
kT
−1
,
(9)
где k – постоянная Больцмана; T – температура окружающей среды.
В представлении Лиувилля [31]:
Aˆ Aˆ ,
Aˆ Aˆ ,
система (5) может быть записана как
(10)
15
ˆ
t ( t ) = L ( t )
,
( t ) = 0
t =0
где
(
(11)
)
Lˆ = −i aˆ †aˆ − aˆ †aˆ + Гˆ ,
(
(
) (
(12)
))
ˆ ˆ † + nT aa
ˆ ˆ † + aa
ˆ ˆ † − 2aˆ † aˆ .
Гˆ = − ( nT + 1) aˆ †aˆ + aˆ †aˆ − 2aa
(13)
2
Задача поиска решения системы (9) нетривиальна и достаточно сложна:
возможны подходы как непосредственного поиска
( t ) [30], так и
подстановки в выражения операторов интересуемых физических величин
(например, оператор Â ) для последующего получения средних значений –
наблюдаемых [32]:
ˆ ˆ ( t ) → Aˆ ( t ) .
tr Aˆ ( t ) = tr AL
t
(14)
16
2
ДИНАМИКА ДИСПЕРСИИ СРЕДНЕГО ЧИСЛА ФОТОНОВ
2.1 Уравнение динамики дисперсии среднего числа фотонов состояния
Первым рассматриваемым наблюдаемым будет дисперсия среднего
числа фотонов некоторого квантового состояния:
n ( t ) = nˆ 2 ( t ) − nˆ ( t ) ,
2
(15)
где n̂ – оператор среднего числа фотонов [33]:
nˆ = aˆ † aˆ.
(16)
Именно эта физическая величина была выбрана для исследования, так как
анализ дисперсии среднего числа фотонов состояния позволяет оценить
изменение мощности оптического излучения.
Тогда в уравнения системы (11) можно произвести следующую подстановку:
nˆ ( t ) = x ( t ) ,
ˆ ˆ (t ) ,
x ( t ) = t r nL
t
(17)
в этом случае исследуемое наблюдаемое – среднее число фотонов состояния;
и
nˆ 2 ( t ) = y ( t ) ,
y ( t ) = t r nˆ 2 Lˆ ( t ) ,
t
(18)
здесь исследуемое наблюдаемое – квадрат среднего числа фотонов состояния.
Подстановка (17), (18) в (11) позволяет получить две системы
уравнений, состоящие из линейного дифференциального уравнения первого
порядка и начального условия для двух наблюдаемых – две задачи Коши.
Решение полученных систем имеет вид:
x ( t ) = nT (1 − e− t ) + x0e− t ,
y ( t ) = (1 + 2nT ) nT + (1 + 4nT ) ( x0 − nT ) e− t + ( y0 − (1 + 4nT ) x0 + 2nT 2 ) e −2 t ,
(19)
(20)
где x0 , y0 – константы, вид которых определяется для конкретного состояния
условиями:
17
x0 = x ( t ) t =0 = nˆ ( t )
t =0
y0 = y ( t ) t =0 = nˆ 2 ( t )
= t r nˆ 0 ,
t =0
(21)
= t r nˆ 2 0 ,
Возможно произвести проверку полученных выше выражений x ( t ) и y ( t ) ,
сравнив с предельными значениями выражений, составляющими:
x ( t ) t → = t r nˆ T ( t ) = nT ,
(22)
y ( t ) t → = t r nˆ 2 T ( t ) = nT (1 + 2nT ) ,
где T – матрица плотности теплового состояния:
1
T =
1 + nT
m
nT
m m.
1
+
n
m
T
(23)
2.2 Динамика дисперсии среднего числа фотонов сжатого вакуумного
состояния
Вакуумное
состояние
представляет
собой
состояние
моды
электромагнитного излучения с наименьшей энергией и точно определенным
числом фотонов в моде, равным 0 [33]. Матрица плотности вакуумного
состояния имеет вид:
0(вак) = 0 0 .
(24)
Интересен переход к рассмотрению сжатого вакуумного состояния, так как его
наблюдение обнаруживает ненулевое среднее число фотонов в моде, но все
еще нулевую амплитуду.
Для сжатия состояния на него нужно подействовать соответствующим
оператором:
1
2
aˆ
Sˆ ( ) = e 2
1
− aˆ †2
2
,
(25)
где – комплексный параметр сжатия: = r ei .
Действие оператора сжатия на вакуумное состояние приводит к появлению
сжатого вакуумного состояния, имеющего в формализме матриц плотности
следующий вид:
18
0(сж.вак) = Sˆ ( ) 0 0 Sˆ † ( ).
(26)
Тогда рассмотрение средних значений t r nˆ 0(сж .вак) и t r nˆ 2 0(сж.вак)
позволяет найти величину констант x 0 и y 0 :
x0 = sinh 2 r ,
y0 = 2cosh 2 r sinh 2 r + sinh 4 r.
(27)
Таким образом, динамика дисперсии сжатого вакуумного состояния зависит
только от модуля параметра сжатия r , выбор фазы параметра сжатия не
оказывает влияния на характер динамики.
Моделирование динамики системы здесь и далее производилось при
использовании следующих значений параметров. Скорость термализации
при значениях параметров
= −0.2 дБ/км, = 1.5, характеризующих
стандартное телекоммуникационное оптоволокно, например, SMF-28 [32],
согласно выражению (8) составило = 10 кГц . В соответствии с выражением
(9) для длины волны излучения 1550 нм и температуры окружающей среды
300 К среднее число тепловых фотонов составило nT = 10−13 . Для удобства
графики приведены с нормированной шкалой t по оси абсцисс.
Полученная динамика дисперсии среднего числа фотонов n ( t ) сжатого
вакуумного состояния (сж .вак) ( t ) для разных значений модуля параметра
сжатия r (рис. 1) приводится в линейном и логарифмическом масштабах (рис.
2, 3). Можно отметить, что при r = 0 , т.е. при отсутствии сжатия вакуумного
состояния, динамика дисперсии среднего числа фотонов соответствует
термализации системы: приходу некоторого числа тепловых фотонов при
изначальном нулевом значении среднего числа фотонов в системе (рис. 4).
Было определено, что на малом временном интервале (рис.2, 4)
динамика системы может быть описана выражением:
n ( t )
(y
0
− x02 ) e −2 t ,
(28)
19
тогда как на большем временном интервале (рис.3, 4) динамика соответствует:
n ( t )
Представляет
интерес
x0e − t .
рассмотрение
(29)
пересечения
кривых
этих
приближенных значений динамики. Так было установлено, что значение по
оси абсцисс, соответствующее пересечению (28) и (29), т.е. соответствующее
смене скорости затухания, смене динамик, зависит от начальных условий
следующим образом:
tq = − ln
x0
.
y0 − x02
(30)
Пересечение линий, соответствующих динамике (29) и значений среднего
числа тепловых фотонов происходит при:
tT = − ln
nT
.
x0
(31)
Рисунок 1 – Зависимость дисперсии среднего числа фотонов n0 сжатого
вакуумного состояния 0(сж .вак) от модуля параметра сжатия r в
логарифмическом масштабе
20
Рисунок 2 – Зависимость динамики дисперсии среднего числа фотонов n ( t )
сжатого вакуумного состояния (сж .вак) ( t ) от модуля параметра сжатия r в
линейном масштабе
Рисунок 3 – Зависимость динамики дисперсии среднего числа фотонов n ( t )
сжатого вакуумного состояния (сж .вак) ( t ) от модуля параметра сжатия r в
логарифмическом масштабе
21
Рисунок 4 –Динамика дисперсии среднего числа фотонов n ( t ) вакуумного
(вак) ( t ) (красный) и сжатого вакуумного состоянии (сж .вак) ( t ) при r = 1
(зеленый) в логарифмическом масштабе, приближение на малом временном
отрезке (синий пунктир), приближение на большом временном отрезке
(черный пунктир)
2.3 Динамика дисперсии среднего числа фотонов сжатого когерентного
состояния
Когерентное
состояние
является
наиболее
близким
описанием
когерентного излучения в формализме квантовой оптики:
2
0(когер) = , = e
−
2
n
n.
n
!
n =0
(32)
Однако сжатое когерентное состояние
0(сж.когер) = Sˆ ( ) Sˆ † ( )
(33)
не имеет классического аналога по своим уникальным особенностям.
22
Значения констант x 0 и y 0 , найденные с помощью средних значений
t r nˆ 0(сж.когер) и t r nˆ 2 0(сж.когер) , составили:
x0 =
2
( cosh
2
r + sinh 2 r ) − ( ) ei cosh r sinh r −
2
− 2e − i cosh r sinh r + sinh 2 r ,
(1 + ) cosh r + (1 + 3
( ) 2 + e
(
)
y0 =
2
2
4
2
+
4
) sinh r −
4
2
( ) 2
−
cosh 3 r sinh r +
2
2
e − i
+ ( ) + 2
2
+ 2 +
4 + 4 e 2i + 4 e −2i
cosh 2 r sinh 2 r −
+
2
2
2
4
+ 2 + + 2 + 4 +
( )2 4 + 2 2 ei +
cosh r sinh 3 r −
−
2
2
− i
4 + 2 e
2
i
2
(
(
(
(
(
− ( ) 2 +
2
)
) (
)
2
(34)
)
)
).
Параметр сжатия, модуль параметра сжатия r и фаза параметра сжатия
, оказывает значительное влияние на начальные условия рассмотрения
динамики дисперсии среднего числа фотонов n ( t ) , что приведено на рис.5 и
рис.6. Поэтому динамика дисперсии среднего числа фотонов n ( t ) сжатого
когерентного состояния
(сж .когер) ( t ) была рассмотрена для различных
значений модуля параметра сжатия r и фазы параметра сжатия в линейном
и логарифмическом масштабах (рис.7–10).
Сравнение динамики дисперсии среднего числа фотонов
n ( t )
когерентного (когер) ( t ) и сжатого когерентного (сж .когер) ( t ) состояний в
логарифмическом (рис.11) масштабах показывает, что сжатое когерентное
состояние, как и сжатое вакуумное состояние, рассмотренное выше,
23
характеризуется на малом временном интервале быстрым затуханием,
описываемым выражением (29). Динамика дисперсии среднего числа фотонов
когерентного состояния, не претерпевшего сжатия, описывается выражением
(28), как и динамика дисперсии среднего числа фотонов сжатого когерентного
состояния на большем временном отрезке.
Рисунок 5 – Зависимость дисперсии среднего числа фотонов n0 сжатого
когерентного состояния 0(сж .когер) от модуля параметра сжатия r при
различных значениях фазы параметра сжатия в логарифмическом масштабе
Рисунок 6 – Зависимость дисперсии среднего числа фотонов n0 сжатого
когерентного состояния 0(сж .когер) от фазы параметра сжатия при различных
значениях модуля параметра сжатия r в линейном масштабе
24
Рисунок 7 – Зависимость динамики дисперсии среднего числа фотонов n ( t )
сжатого когерентного состояния (сж .когер) ( t ) от модуля параметра сжатия r
при значении фазы параметра сжатия =
3
в линейном масштабе
4
25
Рисунок 8 – Зависимость динамики дисперсии среднего числа фотонов n ( t )
сжатого когерентного состояния (сж .когер) ( t ) от модуля параметра сжатия r
при значении фазы параметра сжатия =
3
в логарифмическом масштабе
4
26
Рисунок 9 – Зависимость динамики дисперсии среднего числа фотонов n ( t )
сжатого когерентного состояния (сж .когер) ( t ) от фазы параметра сжатия при
значении модуля параметра сжатия r = 1 в линейном масштабе
27
Рисунок 10 – Зависимость динамики дисперсии среднего числа фотонов
n ( t ) сжатого когерентного состояния (сж .когер) ( t ) от фазы параметра сжатия
при значении модуля параметра сжатия r = 1 в логарифмическом масштабе
28
Рисунок 11 –Динамика дисперсии среднего числа фотонов n ( t )
когерентного состояния (когер) ( t ) (красный) и сжатого когерентного
состояния (сж .когер) ( t ) при r = 1 и =
3
(зеленый) в логарифмическом
4
масштабе, приближение на малом временном отрезке (синий пунктир),
приближение на большом временном отрезке (черный пунктир)
2.4 Динамика
дисперсии
среднего
числа
фотонов
суперпозиции
когерентных состояний
Суперпозиция когерентных состояний
=
(
1
2 1 + cos ( ) e
−2
2
)
(
+ ei − ) ,
(35)
в представлении матрицы плотности имеет вид:
0(суперпоз.) =
1
(
+ e−i − + ei − + − − ) ,
(36)
29
где нормировочный коэффициент
составляет:
(
).
2
= 2 1+ e
−2
2
cos
(37)
Это состояние также являет собой пример уникального квантового эффекта,
заключающегося
в
нахождении
квантового
состояния
в
двух
взаимоисключающих состояниях одновременно.
= 2 n, n
При значении разности фаз
такую суперпозицию
называют четной (при разложении этого состояния по базису Фока сумма
содержит только слагаемые с четными значениями числа фотонов в моде):
even
=
1
(
2 1+ e
тогда как при = (1 + 2n ) , n
−2
2
)
(
+ − ) ,
(38)
нечетной (при разложении этого состояния
по базису Фока сумма содержит только слагаемые с нечетными значениями
числа фотонов в моде):
odd
=
1
(
2 1− e
−2
2
)
(
− − ) .
(39)
Рассмотрение средних значений t r nˆ 0(суперпоз.) и t r nˆ 2 0(суперпоз.)
позволило найти необходимые константы x 0 и y 0 :
x0 =
2
2
−2
2
2
2
(
2 1+ e
y0 =
(
2 1+ e
−2
2
cos
cos
)
)
2
2
(
1− e
−2
(1 +
2
2
)
cos ,
−e
−2
2
( cos −
2
))
(40)
cos .
30
Динамика дисперсии среднего числа фотонов n ( t ) суперпозиции
когерентных состояний (суперпоз.) ( t ) приведена для различных значений
разности фаз (рис. 12) в линейном (рис. 13) и логарифмическом (рис. 14)
масштабах.
Рассмотрение динамики дисперсии среднего числа фотонов
суперпозиции когерентных состояний в логарифмическом масштабе (рис. 15)
позволяет определить, что выражение (29) и в этом случае способно дать
описание динамики на большем временном интервале. Нетривиальным
результатом оказалось наблюдение увеличения значения дисперсии среднего
числа фотонов при нечётной суперпозиции когерентных состояний, что
показано на рис. 15.
Рисунок 12 – Зависимость дисперсии среднего числа фотонов n0
суперпозиции когерентных состояний 0(суперпоз.) от разности фаз в
линейном масштабе
31
Рисунок 13 – Зависимость динамики дисперсии среднего числа фотонов
n ( t ) суперпозиции когерентных состояний (суперпоз.) ( t ) от разности фаз в
линейном масштабе
Рисунок 14 – Зависимость динамики дисперсии среднего числа фотонов
n ( t ) суперпозиции когерентных состояний (суперпоз.) ( t ) от разности фаз в
логарифмическом масштабе
32
Рисунок 15 – Динамика дисперсии среднего числа фотонов n ( t )
суперпозиции когерентных состояний (суперпоз.) ( t ) при значении разности
фаз = (красный) в логарифмическом масштабе, приближение на большом
временном отрезке (черный пунктир)
2.5 Динамика дисперсии среднего числа фотонов состояния: основные
выводы
В данной части работы была рассмотрена динамика дисперсии среднего
числа фотонов квантового состояния. С помощью тории открытых квантовых
систем и формализма матриц плотности были выведены аналитические
выражения, способные дать описание этой динамики.
Было установлено, что динамика дисперсии среднего числа фотонов
характеризуется экспоненциальным затуханием. Скорость затухания различна
на двух временных отрезках. Были найдены аппроксимирующие динамик на
двух временных отрезках.
В качестве примеров были приведены численные модели динамики
дисперсии среднего числа фотонов сжатого вакуумного, сжатого когерентного
состояний и суперпозиции когерентных состояний. Было установлено влияние
33
значений параметра сжатия на начальные значения величины дисперсии: при
увеличении по модулю значения модуля параметра сжатия r происходит
увеличение значений дисперсии среднего числа фотонов.
С помощью выражений (30) и (31) были определены значения дистанции
Lq – расстояние, на котором происходит смена динамики затухания значения
дисперсии среднего числа фотонов, и LT – расстояние, на котором можно
говорить о термализации системы. Стоит отметить, что полная термализация
системы – присутствие только тепловых фотонов в моде – наступает
асимптотически при стремлении времени распространения к бесконечности.
Поэтому за фактический момент термализации далее будет принята
координата по оси абсцисс пересечения приближения динамики дисперсии
(29) и значения среднего числа тепловых фотонов. Было определено, что для
сжатого вакуумного состояния, характеризующегося r = 1 характерные
расстояния составили: Lq (сж .вак)
34
км,
LT (сж .вак)
690
когерентного состояния, характеризующегося r = 1 и =
км; для сжатого
3
: Lq (сж .когер)
4
36 км,
LT (сж.когер)
700 км; для суперпозиции когерентных состояний при =
LT (суперпоз.)
630 км.
3
:
4
34
3 ДИНАМИКА НАБЛЮДАЕМЫХ ОПЕРАТОРОВ КВАДРАТУР
3.1 Уравнение динамики средних значений наблюдаемых операторов
квадратур состояния
Представляет интерес рассмотрение средних значений наблюдаемых
операторов квадратур:
1 †
( aˆ + aˆ ) ,
2
i
pˆ =
aˆ † − aˆ ) ,
(
2
qˆ =
(41)
так как операторы квадратур q̂ и p̂ представляют вещественную и мнимую
части оператора комплексной амплитуды aˆ =
1
( qˆ + ipˆ ) .
2
Используя подход, приведенный в выражении (14), на основе уравнений
(11) – (13) были выведены уравнения динамики средних значений
наблюдаемых операторов квадратур:
q (t ) = e
− t
2
p (t ) = e
− t
2
( q0 cos t + p0 sin t ) ,
( p0 cos t − q0 sin t ) ,
(42)
где константы q 0 и p0 могут быть найдены для конкретного состояния из
начальных условий:
q0 = qˆ ( t )
p0 = pˆ ( t )
t =0
= t r qˆ 0 ,
t =0
= t r pˆ 0 .
(43)
Далее будут рассмотрены примеры анализа динамики средних значений
наблюдаемых операторов квадратур сжатого когерентного состояния и
суперпозиции когерентных состояний. Сжатое вакуумное состояние имеет
35
нулевую амплитуду = 0 , как было отмечено ранее, поэтому не представляет
интереса в рассмотрении.
3.2 Динамика средних значений наблюдаемых операторов квадратур
сжатого когерентного состояния
Приведенная ранее матрица плотности сжатого когерентного состояния
0(сж .когер) позволяет найти константы q 0 и p0 с помощью применения
операции взятия следа t r qˆ 0(сж.когер) и t r pˆ 0(сж.когер) :
1
+ ) cosh r − ( e−i + ei ) sinh r ,
(
2
i
p0 =
− ) cosh r − ( e−i − ei ) sinh r .
(
2
q0 =
(44)
Динамика средних значений наблюдаемых операторов квадратур q ( t ) и
p ( t ) сжатого когерентного состояния (сж .когер) ( t ) находится в зависимости от
модуля параметра сжатия r и фазы параметра сжатия . На рис. 16, 17
приведено численное моделирование зависимости динамки средних значений
наблюдаемых операторов квадратур q ( t ) и p ( t ) от значений параметра
сжатия при = 10 рад/с для демонстрации характера осцилляций. Значение ,
соответствующее рассматриваемой длине волны излучения 1550 нм, из-за
очень большой величины неудобно для численного расчета. Главный интерес
при анализе динамики представляет именно огибающая осцилляций. С
течением времени термализация системы приводит к затуханию осцилляций,
из-за чего средние значения наблюдаемых асимптотически стремятся к
нулевым значениям (рис. 18). Из анализа выражений (42) было установлено,
что динамика средних значений наблюдаемых операторов квадратур q ( t ) и
p ( t ) характеризуется экспоненциальным затуханием с коэффициентом,
зависящим от начальных условий следующим образом:
36
f огиб (сжат.когер) ( t ) = q0 2 + p0 2 e
− t
2
.
(45)
Анализ выражения огибающей осцилляций (45) позволил дать следующую
оценку: за отрезок безразмерного времени
t = 2
(46)
происходит уменьшение амплитуды осцилляций в e раз.
Рисунок 16 – Зависимость динамики средних значений наблюдаемых
операторов квадратур q ( t ) (синий) и p ( t ) (красный) сжатого когерентного
состояния (сж .когер) ( t ) от модуля параметра сжатия r при значении фазы
параметра сжатия =
3
в линейном масштабе
4
37
Рисунок 17 – Зависимость динамики средних значений наблюдаемых
операторов квадратур q ( t ) (синий) и p ( t ) (красный) сжатого когерентного
состояния (сж .когер) ( t ) от фазы параметра сжатия при значении модуля
параметра сжатия r = 1 в линейном масштабе
38
Рисунок 18 – Огибающая динамики среднего значения наблюдаемого
оператора квадратуры q ( t ) сжатого когерентного состояния (сж .когер) ( t ) при
значении модуля параметра сжатия r = 1 и фазы параметра сжатия =
3
4
осциллирующих значений f огиб (сжат.когер) ( t ) в линейном масштабе
3.3 Динамика средних значений наблюдаемых операторов квадратур
суперпозиции когерентных состояний
Значения констант q 0 и p0 для суперпозиции когерентных состояний
0(суперпоз.) , найденные с помощью t r qˆ 0(суперпоз.) и t r pˆ 0(суперпоз.) , составили:
39
q0 =
2ie
(
2 1+ e
p0 =
−2
−2
− 2e
(
2 1+ e
2
cos
−2
−2
2
2
)
2
2
cos
)
2
(
(
− ) sin ,
+ ) sin .
(47)
Приведённая динамика средних значений наблюдаемых операторов
квадратур q ( t ) и p ( t ) суперпозиции когерентных состояний (суперпоз.) ( t ) ,
находящаяся в зависимости от разности фаз (рис.19) также имеет характер
осцилляций, экспоненциально затухающих с течением времени (рис.20).
Аналитический вид огибающей экспоненциально затухающих колебаний
является аналогичным рассмотренному выше выражению (45) огибающей для
случая сжатых когерентных состояний. Так же, как и в случае сжатых
когерентных состояний, уменьшение амплитуды осцилляций в e раз по
сравнению с начальным значением происходит при координате по оси
абсцисс, описываемой выражением (46).
40
Рисунок 19 – Зависимость динамики средних значений наблюдаемых
операторов квадратур q ( t ) (синий) и p ( t ) (красный) суперпозиции
когерентных состояний (суперпоз.) ( t ) от разности фаз в линейном масштабе
41
Рисунок 20 – Огибающая динамики среднего значения наблюдаемого
оператора квадратуры q ( t ) суперпозиции когерентных состояний
(суперпоз.) ( t ) при значении разности фаз =
3
в линейном масштабе
2
3.4 Динамика среднего значения наблюдаемых операторов квадратур
состояния: основные выводы
В ходе исследования средних значений наблюдаемых операторов
квадратур
подтвердился
общий
предполагаемый
характер
динамики:
затухающие осцилляции. Удалось получить уточнение, показывающее, что
затухание осцилляций является экспоненциальным с показателем, равным
половине скорости термализации.
Было найдено аналитическое выражение, описывающее огибающую
затухающих осцилляций, что позволяет определить, например, расстояние Lе ,
за которое происходит уменьшение величины амплитуды осцилляций в e раз.
Как для случая сжатых когерентных состояний, так и для случая суперпозиции
когерентных состояний оно составило порядка 45 км.
42
ЗАКЛЮЧЕНИЕ
Описание динамики открытых квантовых систем в настоящее время
является комплексной и при этом высокомотивированной задачей. Несмотря
на то, что этот вопрос в исторической ретроспективе имеет значительную
историю развития, и сейчас остаются вопросы, требующие нахождения
решения.
Методы, предлагаемые исследователями на протяжении развития этой
области,
достаточно
разнообразны,
и
их
выбор
аргументирован
характеристиками конкретного случая. Кроме того, с помощью приближений
возможно осуществить переход от более общих уравнений, описывающих
эволюцию
открытых
квантовых
систем,
к
рассмотрению
частных
интересующих примеров.
Так
в
ходе
работы
был
произведен
анализ
особенностей
распространения неклассических состояний света в оптоволоконном канале, в
рамках приближений, вводимых решением уравнения Линдблада. Были
выведены и проанализированы выражения, описывающие динамику полевых
наблюдаемых: дисперсии среднего числа фотонов и средних значений
операторов квадратур.
Было установлено, что динамика дисперсии среднего числа фотонов
характеризуется экспоненциальным затуханием, скорость которого различна
на двух временных отрезках. Было определено, что динамика средних
значений операторов квадратур имеет характер затухающих осцилляций.
Аналитическое выражение огибающей осцилляций было найдено. Найденные
аппроксимирующие динамик позволили определить характерные расстояния
распространения неклассических состояний: за
Lе ~45 км происходит
уменьшение амплитуды средних значений операторов квадратур в е раз, на
расстоянии Lq ~30-40 км происходит смена динамики убывания значений
дисперсии среднего числа фотонов, а на LT ~600-700 км можно говорить о
43
полной термализации состояния. Приведены интервалы значений, так как на
величину расстояния распространения влияют амплитуда состояний, значения
параметра сжатия или разности фаз.
Полученные результаты будут важны для учета явлений диссипации, а
именно термализации, и декогеренции в открытых квантовых системах при
рассмотрении широкого круга применений, например, квантовых протоколов
шифрования,
алгоритмов
квантовых
вычислений,
гомодинном
детектировании и создании систем квантового имаджинга. Таким образом,
произведенная оценка характерных времён распространения квантовых
состояний через оптоволоконный канал представляет интерес как с точки
зрения фундаментальных исследований, так и прикладных разработок.
В дальнейшей работе планируется проведение эксперимента для
сопоставления экспериментально полученных данных с теоретическими.
Также представляет интерес анализ динамики дисперсии средних значений
операторов квадратур для рассмотрения принципа неопределённости
Гейзенберга и уточнения времени декогеренции. Кроме того, в силу
универсальности метода планируется рассмотрение других операторов и
возможности его применения для случая двумодовых и многомодовых
состояний.
44
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1. Acín A. et al. The quantum technologies roadmap: a European community
view //New Journal of Physics. – 2018. – Т. 20. – №. 8. – С. 080201.
2. Pirandola S. et al. Advances in quantum cryptography //Advances in Optics
and Photonics. – 2020. – Т. 12. – №. 4. – С. 1012-1236.
3. Georgescu I. M., Ashhab S., Nori F. Quantum simulation //Reviews of
Modern Physics. – 2014. – Т. 86. – №. 1. – С. 153.
4. Preskill J. Quantum Computing in the NISQ era and beyond //Quantum. –
2018. – Т. 2. – С. 79.
5. Giovannetti V., Lloyd S., Maccone L. Advances in quantum metrology
//Nature photonics. – 2011. – Т. 5. – №. 4. – С. 222.
6. Joos E. et al. Decoherence and the appearance of a classical world in quantum
theory. – Springer Science & Business Media, 2013.
7. Breuer H. P. et al. The theory of open quantum systems. – Oxford University
Press on Demand, 2002.
8. Berkelbach T. C., Thoss M. Special topic on dynamics of open quantum
systems. – 2020.
9. Wang H., Thoss M. Multilayer formulation of the multiconfiguration timedependent Hartree theory //The Journal of chemical physics. – 2003. – Т. 119.
– №. 3. – С. 1289-1299.
10.Picconi D., Burghardt I. Open system dynamics using Gaussian-based
multiconfigurational time-dependent Hartree wavefunctions: Application to
environment-modulated tunneling //The Journal of chemical physics. – 2019.
– Т. 150. – №. 22. – С. 224106.
11.Bouakline F., Fischer E. W., Saalfrank P. A quantum-mechanical tier model
for phonon-driven vibrational relaxation dynamics of adsorbates at surfaces
//The Journal of chemical physics. – 2019. – Т. 150. – №. 24. – С. 244105.
12.Werther M. et al. Davydov-Ansatz for Landau-Zener-Stueckelberg-Majorana
transitions in an environment: Tuning the survival probability via number
45
state excitation //The Journal of chemical physics. – 2019. – Т. 150. – №. 23.
– С. 234109.
13.Shastry A., Xu Y., Stafford C. A. The third law of thermodynamics in open
quantum systems //The Journal of Chemical Physics. – 2019. – Т. 151. – №.
6. – С. 064115.
14.Carmichael H. An open systems approach to quantum optics: lectures
presented at the Université Libre de Bruxelles, October 28 to November 4,
1991. – Springer Science & Business Media, 2009. – Т. 18.
15.Zwanzig R. Nonequilibrium statistical mechanics. – Oxford University Press,
2001.
16.Nakajima S. On quantum theory of transport phenomena: steady diffusion
//Progress of Theoretical Physics. – 1958. – Т. 20. – №. 6. – С. 948-959.
17.Zwanzig R. Ensemble method in the theory of irreversibility //The Journal of
Chemical Physics. – 1960. – Т. 33. – №. 5. – С. 1338-1341.
18.Xu M. et al. Convergence of high order memory kernels in the NakajimaZwanzig generalized master equation and rate constants: Case study of the
spin-boson model //The Journal of chemical physics. – 2018. – Т. 148. – №.
16. – С. 164101.
19.Albert V. V. et al. Lindbladians with multiple steady states: theory and
applications //Bulletin of the American Physical Society. – 2017. – Т. 62.
20.Redfield A. G. The theory of relaxation processes //Advances in Magnetic and
Optical Resonance. – Academic Press, 1965. – Т. 1. – С. 1-32.
21.Trushechkin A. S. Higher-order corrections to the Redfield equation with
respect to the system-bath coupling based on the hierarchical equations of
motion //Lobachevskii Journal of Mathematics. – 2019. – Т. 40. – №. 10. – С.
1606-1618.
22.Zhang W. M. Exact master equation and general non-Markovian dynamics in
open quantum systems //The European Physical Journal Special Topics. –
2019. – Т. 227. – №. 15. – С. 1849-1867.
46
23.Megier N., Smirne A., Vacchini B. Evolution equations for quantum semiMarkov dynamics //Entropy. – 2020. – Т. 22. – №. 7. – С. 796.
24.Gorini V., Kossakowski A., Sudarshan E. C. G. Completely positive
dynamical semigroups of N‐level systems //Journal of Mathematical Physics.
– 1976. – Т. 17. – №. 5. – С. 821-825.
25.Lindblad G. On the generators of quantum dynamical semigroups
//Communications in Mathematical Physics. – 1976. – Т. 48. – №. 2. – С. 119130.
26.Manzano D. A short introduction to the Lindblad master equation //AIP
Advances. – 2020. – Т. 10. – №. 2. – С. 025106.
27.Teretenkov A. E. Dynamics of moments for quadratic GKSL generators
//Mathematical Notes. – 2019. – Т. 106. – №. 1. – С. 151-155.
28.Teretenkov A. E. Irreversible quantum evolution with quadratic generator
//Infinite Dimensional Analysis, Quantum Probability and Related Topics. –
2019. – Т. 22. – №. 04. – С. 1930001.
29.Trushechkin A. S. Decoherence and Coherence Preservation in the Solutions
of the GKSL Equation in the Theory of Open Quantum Systems
//Mathematical Notes. – 2019. – Т. 106. – №. 5. – С. 986-993.
30.Gaidash A., Kozubov A., Miroshnichenko G. Dissipative dynamics of
quantum states in the fiber channel //Physical Review A. – 2020. – Т. 102. –
№. 2. – С. 023711.
31.Royer A. Aspects of open quantum dynamics //Irreversible Quantum
Dynamics. – Springer, Berlin, Heidelberg, 2003. – С. 47-63.
32.Gaidash A. et al. The Influence of Signal Polarization on Quantum Bit Error
Rate for Subcarrier Wave Quantum Key Distribution Protocol //Entropy. –
2020. – Т. 22. – №. 12. – С. 1393.
33.Scully M. O., Zubairy M. S. Quantum optics. – 1999.
47
Отзывы:
Авторизуйтесь, чтобы оставить отзыв