'

Квантовые траектории в диссипативной динамике одиночного кубита

Понравилась презентация – покажи это...





Слайд 0

Квантовые траектории в диссипативной динамике одиночного кубита Гельман Александр Иосифович, нс 170 отд, ИПФ РАН, г. Н. Новгород


Слайд 1

Предмет исследования и актуальность темы Квантовая оптика Квантовая информация Нанотехнологии Квантовые эффекты становятся определяющими в работе современных высокотехнологичных устройств: Полупроводниковая техника Компьютеры (процессоры) Сверхпроводящие объекты Вся наноэлектроника Закон Мура Сверхсовременное направление – квантовая информация ROADMAP №1, Europe, www.qurope.net, апрель 2010 ROADMAP №2, USA, www.ostp.gov, 2009, 2004 1. Квантовая криптография Принципиально не дешифруемая передача данных, доведено до уровня готовых приборов; www.idquantique.com, www.magiqtech.com 2. Квантовый компьютер Квантовые сверхскоростные вычисления, взлом существующих систем кодирования Основной рабочий элемент – квантовый бит, кубит


Слайд 2

Квантовая информация


Слайд 3

Среднесрочные планы


Слайд 4

Постановка задачи Кубит является основным рабочим элементом квантового компьютера, физически реализуется в виде двухуровневой квантовой системы, например фотона, атома, донора в полупроводнике, ядерного спина, серхпроводящего контура Важнейшей характеристикой кубита является время декогерентности Наиболее перспективные устройства для реализации квантового компьютера - кубиты на основе джозефсоновских переходов M. Nakahara and T. Ohmi Quantum computing: from linear algebra to physical realizations. – London, 2008 Проблемы: - измерение скоростей релаксации кубита - управление динамикой кубитов в сильном переменном поле - измерения состояния кубита в условиях шума


Слайд 5

J.E. Mooij, et.al,Science 285,1036 (1999) Yu. Makhlin, et.al., Rev. Mod.Phys. 73, 357 (2001) 3JJ qubit: гамильтониан Искусственный атом W.D.Oliver,et.al.,Quant inf Process 8,261(2009) Данные состояния могут быть измерены, соответствуют току в кубите по- и против часовой стрелки


Слайд 6

Шум в системе M.Sillanpaa, et al.,PRL 96,187002(2006),W.D.Oliver,et.al., Science 310,1653(2005), D. M. Berns, et al., Nature 455,51(2008) Механизмы релаксации в кубите: Флуктуация заряда на джозефсоновских контактах Квазичастицы на островках сверхпроводимости (конечное сопротивление) Ядерные спины в подложке (флуктуация магнитного поля) Радиационное затухание, связь с управляющим полем Флуктуация потока Метод Монте-Карло (квантовых траекторий) C. W. Gardiner, P. Zoller. Quantum noise. – Berlin: Springer, 2000. M. B. Plenio, P. L. Knight. Rev. Mod. Phys. 1998. V. 70, №1. P. 101–143. H.-P. Breuer, F. Petruccione. The theory of open quantum systems. – Cambridge: Oxford University Press, 2002. А.И. Гельман, В.А. Миронов. Численное моделирование квантовой релаксации в многоуровневых атомных системах методом Монте-Карло: препринт №773; ИПФ РАН. – Н. Новгород, 2008. – 42 с


Слайд 7

Квантовая теория релаксации: метод Монте-Карло Общий вид уравнения для оператора плотности в марковском приближении - форма Линдблада неэрмитов эффективный Гамильтониан оператор квантового скачка число скоростей релакс. оператор системы, взаимодействующий с резервуаром


Слайд 8

1. Вероятность излучения фотона нет излучения излучился фотон 3. повтор шагов 1,2 для расчета эволюции за требуемое время – квантовая траектория 4. повтор шагов 1-3 n раз для получения статистического ансамбля квантовых траекторий – усредненная динамика системы, соответствует уравнению для оператора плотности Метод МК: алгоритм численного моделирования Расчет ведется для волновой функции системы Доказательство эквивалентности Усредняем по n реализациям 2. генерация случайного числа r, равномерно распределенного на отрезке [0,1]


Слайд 9

Метод МК: преимущества Квантовая теорема регрессии K. Molmer,Y. Castin, J. Dalibard, J. Opt. Soc. Am. B, 10, 524 (1993) Расчет среднего Расчет двухвременных корреляционных функций Расчет единичных реализаций процессов в квантовых системах Расчет диссипативной динамики одиночных квантовых систем Размерность системы ДУ для волновой функции N вместо N2 для матрицы плотности Простота и эффективность распараллеливания алгоритма Расчет систем в немарковском приближении W. T. Strunz et.al, Phys. Rev. Lett, 82, 1999 Точность Статистически независимые реализации


Слайд 10

Результаты тестирования 1. Двухуровневый атом под действием резонансного поля, в вакууме 0 1 2 3 0 1 2 Гt Гt 10000 реализаций Динамика одного атома (одна из реализаций метода) Усредненная динамика системы (ансамбля атомов) 1000 реализаций 0 3


Слайд 11

3JJ qubit: квантовые траектории Эффективный гамильтониан в резонансном приближении С увеличением Г динамика существенно меняется. Кубит может возбуждаться на верхний уровень. Нет пленения на ~Т/2 (между пересечениями уровней). Даже при сильном поле влияние шума существенно КПТ Rabi+LZ КПТ Rabi+LZ


Слайд 12

3JJ qubit: усредненная динамика После усреднения по 3000 реализация метода МК (соответствует условиям эксперимента) наблюдается классическая Раби - динамика Видно отличие релаксационной динамики в одной реализации от усредненной динамики, когда наблюдается насыщение и выход населенности на стационарное значение A. Gelman, A.M. Satanin, JETP lett. 91, 535-540 (2010) Совпадение с моделью и экспериментом: D. M. Berns atal. PRL 97,150502(2006) КПТ КПТ


Слайд 13

3JJ qubit: усредненная динамика При большом шуме эффекты когерентности (КПТ) исчезают Перекрытие резонансов Kayanuma, Y. Phys. Rev B., 47, 9940 (1992)


Слайд 14

Приложение к амплитудной спектроскопии N=3000 realizations Подгонка параметров шума при прямом численном моделировании под результаты эксперимента позволит восстановить параметры образца с хорошей точностью


Слайд 15

В хорошем соответствии с N=3000 в предыдущем рассмотрении и эксперименте, где обычно N=3000-10000. A. I. Gelman, A.M. Satanin, JETP lett. 91, 535-540 (2010). Приложение к амплитудной спектроскопии: измерение состояния кубита Зависимость интерференционной картины от числа реализаций метода (числа измерений в эксперименте).


Слайд 16

Зависимость интерференционной картины от флуктуации начальной фазы импульса P(t) P(t) Равномерное распределение, в каждой реализации - случайное В условиях шума при усреднении по 2000 реализаций населенность по окончании импульса не зависит от начальной фазы импульса. Усреднение по мелкомасштабным осцилляциям Ландау-Зинера. Глобальная динамика (насыщение и выход населенности на стационарное значение) сохраняется D. M. Berns, et.al. Phys. Rev. Lett. 97, 150502 (2006).


Слайд 17

Дана интерпретация переходам ЛЗ в условиях шума и методу АС на примере сверхпроводящих кубитов с точки зрения единичных реализаций, а также установлена связь с усредненной динамикой системы, наблюдающейся в экспериментах. Рассмотрено влияние различного уровня шума на населенности кубитов и зависимость резкости интерференционной картины метода АС от числа измерений состояния кубитов. Показана возможность контрастного формирования интерференционной картины уже при 100 реализациях, что может быть существенным при проведении эксперимента. Промоделирован процесс измерения кубита, включая в рассмотрение классический шум, вызванный флуктуациями начальной фазы возбуждающего импульса в методе АС. Показано что такой шум не влияет на усредненную по реализациям интерференционную картину. Программный комплекс протестирован на известных задачах квантовой оптики, использовались вычислительные кластерные системы ИПФ РАН и ННГУ Проект получил финансовую поддержку на конкурсе «У.М.Н.И.К.», 2010 г. Элементы методики расчета использованы при чтении спецкурса «Квантовая оптика» автором на физическом факультете ННГУ, для студентов 4-5 курса в весеннем семестре 2010 г. Диплом 15-ой Нижегородской сессии молодых ученых (естественные науки), 2010 Основные результаты


Слайд 18

Список публикаций А. И. Гельман, А. М. Сатанин. Релаксационная динамика сверхпроводящих джозефсоновских кубитов в сильном переменном поле // ФТТ. 2010.Т.52.С. 2094 А. И. Гельман, А. М. Сатанин. Квантовые скачки при переходах Ландау-Зинера в диссипативной динамике сверхпроводящего кубита // Письма в ЖЭТФ. 2010. Т. 91. С. 584. А. И. Гельман, А. М. Сатанин. Квантовые скачки при спонтанной релаксации сверхпроводящего кубита под действием сильного ВЧ поля // Вестник ННГУ. 2010. А. И. Гельман, М.В. Денисенко, А. М. Сатанин. Динамический контроль квантовых состояний джозефсоновских кубитов // Вестник ННГУ. 2010. А. И. Гельман, В. А. Миронов. Подавление шума в атомной системе под действием поля в сжатом когерентном состоянии // ЖЭТФ. 2010. Т. 137. С. 627. A. Gelman, V. Mironov. Noise suppression in three-level atomic system driven by quantized field // Proc. SPIE. 2009. V. 7521. P. 75210F. A. I. Gelman, M. V. Denisenko, A. M. Satanin, and F. Nori Amplitude and phase effects in Josephson qubits driven by a biharmonic electromagnetic field (in preparation) A. Shvetsov, A. M. Satanin, A Gelman, A. Zagoskin, S. Savel'ev, and F. Nori. Quantum photonic crystals with tuneable electromagnetic properties (in preparation) A. Shvetsov,A. M. Satanin, A. Gelman, A. Zagoskin, S. Savel'ev, and F. Nori. Control of superconducting metamaterials (in preparation) А. И. Гельман, Диссертация «Диссипативная динамика и контролируемая релаксация в одиночных квантовых системах», ИПФ РАН, 2010


Слайд 19

Квантовая теория релаксации: методы исследования Природа диссипации - взаимодействие системы с резервуаром с большим числом степеней свободы система резервуар взаимодействие система-резервуар Оператор плотности «система+резервуар» Оператор плотности системы C.W.Gardiner, P.Zoller, Quantum noise, Springer, 2000 Скалли М. О., Зубайри М. С., Квантовая оптика, М., Физматлит, 2003


Слайд 20

Решение уравнений для элементов оператора плотности, n2 штук Метод оператора плотности Метод Гейзенберга-Ланжевена, квантовый шум в представлении взаимодействия Решение стохастических дифференциальных уравнений для операторов Соответствует флуктуационно-диссипационной теореме Квантовая теория релаксации: методы исследования Оба метода изначально предполагают описание ансамбля систем


Слайд 21

Резонансная флуоресценция одиночного иона 138Ba+ на переходе |0>-|1>. Оптические поля действуют на оба перехода. Если ион переходит в состояние |2> (метастабильное), флуоресценция исчезает. Через промежуток времени, обычно равный времени жизни уровня 2 (32 сек. в данном эксперименте) атом возвращается в состояние 0, и флуоресценция на переходе |0>-|1> восстанавливается. Одиночный V-атом: флуоресценция t, с Для описания такого поведения потребовалась новая теория (у ансамбля атомов скачков не будет) Bergquist, J. C. et. al. Observation of quantum jumps in a single atom // PRL. 1986. V. 57. P. 1699 Nagourney, W. et.al Shelved optical electron amplifier: Observation of quantum jumps //PRL.1986,56,P.2797 Sauter, T. et. al. Observation of quantum jumps Phys. Rev. Lett. 1986. V. 57. P. 1696. Отсчет, 1/с


Слайд 22

оператор системы, взаимодействующий с резервуаром Метод МК: интерпретация Общий вид уравнения для оператора плотности системы с учетом релаксации в марковском приближении - форма Линдблада эффективный Гамильтониан оператор квантового скачка число скоростей релакс.


Слайд 23

Метод МК: преимущества Квантовая теорема регрессии K. Molmer,Y. Castin, J. Dalibard, J. Opt. Soc. Am. B, 10, 524 (1993) Расчет среднего Расчет двухвременных корреляционных функций Расчет единичных реализаций процессов в квантовых системах Расчет диссипативной динамики одиночных квантовых систем Размерность системы ДУ для волновой функции N вместо N2 для матрицы плотности Простота и эффективность распараллеливания алгоритма Расчет систем в немарковском приближении W. T. Strunz et.al, Phys. Rev. Lett, 82, 1999 Точность Статистически независимые реализации


Слайд 24

Результаты тестирования 1. Двухуровневый атом под действием резонансного поля, в вакууме 1 2 3 1 2 Гt Гt 10000 реализаций Динамика одного атома (одна из реализаций метода) Динамика ансамбля атомов (усредненная динамика) 1000 реализаций 3


Слайд 25

1. Вероятность излучения фотона нет излучения излучился фотон 3. повтор шагов 1,2 для расчета эволюции за требуемое время – квантовая траектория 4. повтор шагов 1-3 n раз для получения статистического ансамбля квантовых траекторий – усредненная динамика системы, соответствует уравнению для оператора плотности Метод МК: алгоритм численного моделирования Расчет ведется для волновой функции системы Доказательство эквивалентности Усредняем по n реализациям 2. генерация случайного числа r, равномерно распределенного на отрезке [0,1]


Слайд 26

Метод МК – простейший вывод Расcмотрим Аналогично В среднем динамика унитарна, существует строгая теория КСДУ, теория измерения


Слайд 27

Известные приложения метода Прочитан спецкурс на Физическом факультете ННГУ, весенний семестр 2010 г. (4-5 курсы) Победа в конкурсе УМНИК: численное моделирование динамики квантовых систем с использованием суперкомпьютерных технологий Лазерное охлаждение (учет импульса) Dalibard, J. Wave-function approach to dissipative processes in quantum optics // Phys. Rev. Lett. 1992. V. 68. P. 580. 2. Взаимодействие квантованного света с атомными системами Nakano, M. // Phys. Rev. A. 2004. V. 70. P. 033407. Гельман, А. И. // ЖЭТФ. 2010. Т. 137. С. 627-636. 3. Устойчивость квантовых протоколов в зашумленной среде (одиночные системы) Jun J. // Phys. Rev. A. 2006. V. 73. P. 064301. Goto, H. et.al.// Phys. Rev. A. 2005. V. 72. P. 054301. Carlo, G. et.al. // Phys. Rev. Lett. 2003. V. 91. P. 257903. Witthout, D. et.al.// Phys. Rev. A 2009. V. 79. P. 033621 4. Динамика Бозе-конденсата в зашумленной среде А.И. Гельман, В.А. Миронов. Численное моделирование квантовой релаксации в многоуровневых атомных системах методом Монте-Карло:препринт №773 ИПФ РАН.2008.(42 с) 5. Динамика сверхпроводящих кубитов в условиях шума A. Gelman, A.M. Satanin, JETP lett. 91, 535-540 (2010)


Слайд 28

Результаты тестирования 2. Трехуровневый атом с ?-конфигурацией электронных уровней под действием резонансных полей, в вакууме Эффект когерентного пленения населенностей N=5000, N=10000, результат неотличим от точного решения Агапьев Б.Д и др., УФН 163, 9 (1993)


Слайд 29

Джозефсоновский переход В зависимости от того, какая энергия больше, выделяют зарядовый и фазовый (потоковый) кубиты (определяется геометрией) Кулоновская энергия Джозефсоновская энергия Число куперовских пар Переход к квантовому описанию, соотношение неопределенностей Сверхпроводящий электрод 1 Сверхпроводящий электрод 2 Изолятор


Слайд 30

Сверхпроводящий потоковый кубит (3JJ qubit) J.E. Mooij, et.al, Science 285, 1036 (1999). Yu. Makhlin, G. Schon, and A. Shnirman, Rev. Mod.Phys. 73, 357 (2001)


Слайд 31

Частица с анизотропной массой в 2d джозефсоновском потенциале самоиндукция 3JJ qubit: гамильтониан


Слайд 32

Движение возможно только в одном направлении – эффективный двухъямный потенциал J.E. Mooij, et.al, Science 285, 1036 (1999). Сверхпроводящий ток в противоположных направлениях, одинаков через каждый переход 3JJ qubit: гамильтониан


Слайд 33

Данные состояния могут быть измерены, соответствуют току в кубите по- и против часовой стрелки Искусственный атом W.D. Oliver, et.al., Quant inf Process 8, 261 (2009) Учет нижнего уровня в каждой яме 3JJ qubit: гамильтониан


Слайд 34

Гамильтониан справедлив для описания динамики всех типов сверхпроводящих кубитов (не только потокового) Различие заключается в способе управления внешними параметрами В случае 3JJ кубита – путем изменения амплитуды внешних полей fdc и fac W.D. Oliver, S.O. Valenzuela Quant Inf Process 8, 261 (2009) 3JJ qubit: управление динамикой Воздействие на кубит внешними полями: постоянным магнитным полем fdc переменным ВЧ электромагнитным полем fac


Слайд 35

3JJ qubit: резонансное приближение Аналогичный гамильтониан широко известен в оптике: резонансное взаимодействие 2-хуровневого атома с частотно-модулированным лазерным полем Средняя по времени населенность верхнего уровня кубит. Ф. Бесселя приводит к интерференционной картине в зависимости населенности от амплитуды постоянного и переменного полей (для определения параметров кубита) G. S. Agarwal et.al., Phys. Rev. A 50, R4465 (1994) S. Ashhab, et.al., Phys. Rev. A 75, 063414 (2007) отличие от оптики


Слайд 36

Вероятность перехода в пределе бесконечно большого времени Скачки населенности происходят при каждом пересечении уровней. Периодическое пересечение приводит к интерференционной картине Основная цель работы – изучение влияния шума на динамику кубита и переходы ЛЗ L. D. Landau, Phys. Z. Sowjetunion 2, 46 (1932) C. Zener. Proc. R. Soc. A137, 696 (1932) W. D. Oliver, Y. Yu, J. C. Lee, et.al., Science 310, 1653 (2005) 3JJ qubit: переходы Ландау-Зинера Впервые исследованы при рассмотрении пересечения уровней при столкновении атомов


Слайд 37

3JJ qubit: квантовый шум в системе Кинетическое уравнения в марковском приближении M. Sillanpaa, et al., Phys. Rev.Lett. 96, 187002 (2006). W. D. Oliver, et al., Science 310, 1653 (2005). D. M. Berns, et al., Phys. Rev. Lett. 97, 150502 (2006). D. M. Berns, et al., Nature 455, 51 (2008).


Слайд 38

3JJ qubit: шум в системе Leggett, A., et.al. (1987) Dynamics of the dissipative two-state system, Rev. Mod. Phys. 59, 1. Schon, G., and A. D. Zaikin, 1990, Phys. Rep. 198, 237; Weiss, U. (1999) Quantum Dissipative Systems


Слайд 39

3JJ qubit: квантовые траектории Эффективный гамильтониан в резонансном приближении С увеличением Г динамика существенно меняется. Кубит может возбуждаться на верхний уровень. Нет пленения населенностей на временах ~Т/2 (между пересечениями уровней) КПТ Rabi+LZ Кубит приготовлен в состоянии |0> J. You, et.al, Phys. Rev. Lett. 100, 047001(2008)


Слайд 40

3JJ qubit: квантовые траектории КПТ Rabi+LZ Кубит приготовлен в состоянии |0> J. You, et.al, Phys. Rev. Lett. 100, 047001(2008) Даже при сильном управляющем поле влияние шума существенно С увеличением Г динамика существенно меняется. Кубит может возбуждаться на верхний уровень. Нет пленения населенностей на временах ~Т/2 (между пересечениями уровней)


Слайд 41

3JJ qubit: усредненная динамика Согласно подходу скоростных уравнений В случае сильного внешнего воздействия: D. M. Berns at.al.PRL 97, 150502 (2006) (кубит на ниж- нем уровне)


Слайд 42

3JJ qubit: усредненная динамика После усреднения по 3000 реализация метода МК (что соответствует условиям эксперимента) наблюдается классическая Раби - динамика Видно отличие релаксационной динамики в одной реализации – стохастическое случайное движение от усредненной динамики, когда наблюдается насыщение и выход населенности на стационарное значение A. Gelman, A.M. Satanin, JETP lett. 91, 535-540 (2010) Совпадение с моделью и экспериментом: D. M. Berns at.al. PRL 97, 150502 (2006) КПТ (динамическая локализация)


Слайд 43

PRL 97, 150502 (2006) D. M. Berns at.al. Параметры кубита могут быть получены из анализа результатов эксперимента (частота Раби, ширина резонанса, населенности уровней) Even for CDT regime После усреднения по 3000 реализация метода МК (что соответствует условиям эксперимента) наблюдается классическая Раби - динамика Видно отличие релаксационной динамики в одной реализации – стохастическое случайное движение от усредненной динамики, когда наблюдается насыщение и выход населенности на стационарное значение A. Gelman, A.M. Satanin, JETP lett. 91, 535-540 (2010) 3JJ qubit: усредненная динамика


Слайд 44

3JJ qubit: усредненная динамика При большом шуме эффекты когерентности (КПТ) исчезают Перекрытие резонансов Kayanuma, Y. Phys. Rev B., 47, 9940 (1992)


Слайд 45

Приложение к амплитудной спектроскопии N=3000 realizations


Слайд 46

Более контрастная картина наблюдается для резонансов высокого порядка. Хорошее совпадение с экспериментом (D.Berns et.al.,PRL 97, 150502 (2006)) Приложение к амплитудной спектроскопии Подгонка параметров шума при прямом численном моделировании под результаты эксперимента позволит восстановить параметры образца с хорошей точностью


Слайд 47

In good correspondence with N=3000 in previous consideration and experiments. In experiments usually N=3000-10000. A. I. Gelman, A.M. Satanin, JETP lett. 91, 535-540 (2010). Приложение к амплитудной спектроскопии: измерение состояния кубита Зависимость интерференционной картины от числа реализаций метода (числа измерений в эксперименте). Расчет без усреднения по времени.


Слайд 48

Приложение к амплитудной спектроскопии: фазовый контроль населенности M. V. Denisenko, A. I. Gelman, A. M. Satanin (in preparation).


×

HTML:





Ссылка: