CC BY
24
ШБН 2079-3316 ПРОГРАММНЫЕ СИСТЕМЫ: ТЕОРИЯ И ПРИЛОЖЕНИЯ № 1(5), 2011, с. 51-61
УДК 517.977
О. В. Батурина, О. В. Моржин Оптимизация управления квантовой системой на модели Ландау—Зинера
Аннотация. Статья посвящена вычислительным экспериментам по применению метода глобального (нелокального) улучшения В.Ф. Кротова к исследованию одной задачи оптимального управления квантовой системой, описываемой уравнением Шредингера с гамильтонианом в форме Ландау— Зинера.
Ключевые слова и фразы: квантовые системы, глобальные улучшения.
Введение
В работах [1,2] для класса общих нелинейных задач оптимального управления дифференциальными системами со свободным правым концом траектории предложен метод глобального улучшения в рамках теории достаточных условий оптимальности [3]. В методе не осуществляется варьирование относительно улучшаемого управления с параметрической оптимизацией, в отличие, например, от методов условного градиента, игольчатого варьирования [4].
Решение задач управления квантовыми системами [5] в контексте современных исследований по нанотехнологиям, квантовым компьютерам [6] представляет актуальную, перспективную тематику. Публикация в [2, 7] результатов моделирования оптимального управления молекулярным состоянием вещества с применением глобального метода вызвала интерес физиков [8] к данному методу.
Цель данной статьи — исследование вопросов применения метода глобального улучшения [2,9-12] для проблемы оптимального управления спинами [8] на примере случая гамильтониана в форме Лан-дау-Зинера [8], причем реализация метода улучшения проводится с линейной по фазовой переменной разрешающей функцией (Кротова), следуя теоретическому аппарату, изложенному в статьях [9-11].
Работа выполнена при финансовой поддержке РФФИ (проект 10-01-90718-моб_ст).
© О. В. Батурина, О. В. Моржин, 2011 © Программные системы: теория и приложения, 2011
1. Постановка задачи
Рассматривается квантовая система [В, 9]
(1) іпФії) = H[u(t)Mt), Ф(о) = е,
где h — постоянная Планка, H[u(t)] —гамильтониан, t Є [О, T]. Функция u(t), t Є [0, T], кусочно-непрерывная, характеризует воздействие внешним полем. В рамках статьи функция u действительная и скалярная. Состояния квантовой системы изучаются на комплексной сфере единичного радиуса [5].
Вводится целевой критерий с вещественнозначным функционалом:
(2) 1 ^^Н1 -\(фG,ф(T ))|2 =
() 1 — (W(T),Wg)(Wg,W(T))^ inf, \фg\2 = 1,
где Wg —заданный вектор, \ • \2 —модуль комплексного числа.
Будем исследовать случай n = 2. Рассматривается гамильтониан в форме Ландау-Зинера, использованный также в работе [ ]: H =
( U Ш I . Система (1) в покомпонентной форме: ш —u I
(3) ihWl = пфі + шф2, ihW2 = шФі — пф2,
Фі(0)= Єї, ф2(0)= Є2.
Целевой критерий (2) в двумерном случае:
(4) I(ф,и) = 1 - |(фа)^i(T) + (Фо)2Ф2(T)|2 ^ inf .
Представляет интерес изучение вычислительных аспектов реализации глобального метода улучшения на примере задачи (3), (4).
2. Конструкции метода глобального улучшения
В [9,10] конструкции глобального метода записаны в терминах комплексных переменных для класса задач оптимального управления квантовыми системами. В контексте работы [11] применительно к изучаемому случаю п = 2 удобно перейти к задаче в терминах действительных переменных. Производится замена фі = г\ + *23, Ф2 = 22 + *24.
Система (3) представляется в виде
21 = 1 (М23 + ^4) , 22 = - (^3 - М^) ,
1 ^ I
(5) 23 = - (-М21 - ^22) , -24 = -(-^21 + «22 ) ,
2(0) = Хз = (Ив фі(0), Ив ф2(0), ІШ Фі(0), 1т ф2(0))Т .
Полагаем в (4) фа = І аі I + * | 61 |, где аі, 6і, а2, 62 —дей-
. . V 62 /
ствительные числа, удовлетворяющие условию а2 + 62 + а2 + 62 = -. Критерий принимает следующий вид:
I(2,и) = 1 - [(а2 + 62)22(Т) + 2(аіа2 + 6і62)21(7^(Т) + +2(аіб2 - 6^2)21 (Т)^4(Т) + (а2 + 62)22(Т)+
(6) +2(а26і - 62аі )22 (Т )2з(т) + (6і + а2)2:з(т) +
+ 2(6і62 + аіа2)23(Т )24 (Т) + (62 + а2)24(Т)) =
= 1 - (2^), Ьг(Т)} ^ М, где матрица коэффициентов
( (а2 + 62) (аіа2 + 6і&2) 0 (аіб2 - 6іа2) \
(аіа2 + 6162) (а2 + 62) (а2бі - б2аі) 0
0 (а26і - 62аі) (62 + а2) (6162 + аіа2)
V (аі62 - 6іа2) 0 (6162 + аіа2) (62 + а2) )
Функция Понтрягина для (5), (6):
Н(р, 2,и) = - [рі(и2з + ^24) + Р2(^^3 - «24) +
I
+Рз(-М2і - ^22) + Р4(-^2і + «23)] .
Сопряженная система
(7)
Рі = Т («Р3 + ^Р4) , Р2 = 7 ФР3 - «Р4) ,
\ \
р3 = -- (ирі + ^Р2) , Р4 = Т («Р2 - ^рі) ,
I I
рі(Т) = 2(а2 + 62)2і(Т) + 2(аіа2 + 6і62)22(7) + 2(аі62 - 6^2)24(7), Р2(Т) = 2(аіа2 + 6162)21(7) + 2(а2 + 62)22(Т) + 2(а26і - 62аі)^3(Т), Р3(Т) = 2(а26і - 62аі)22(Т) + 2(62 + а2)2:3(Т) + 2(6162 + аіа2)24(Т), Р4(Т) = 2(аі62 - 6іа2)2і(Т) + 2(6162 + аіа2)23(Т) + 2(62 + а^гфТ).
Для линейной по u задачи (5), (6) условно вводится ограничение
(8) u(t) G [—v, v], t G [0, T], v> 0,
в плане построения зависимости w(t, z) = v sign K(p(t),z), t G [0,T], где функция переключения K(p(t), z) = ff„(p(t), z, u).
Отображение w(t,z), как известно, записывается в форме
Г — v, K(p(t),z) < 0,
(9) u(t, z) = < любое из [—v, v], K(p(t),z)=0,
1 v, K(p(t),z) > 0.
В случае K(p(t),z) = 0 можно полагать м равным — v или v; при этом частые переключения управления с — v на v на особом режиме могут означать практическую нереализуемость расчетного управления. Поэтому вместо реализации таких переключений будем использовать специальную формулу для управления, предложенную в статьях [10,11].
С вычислительной точки зрения равенство нулю понимается в смысле принадлежности некоторой достаточно малой е-окрестности нуля.
Алгоритм. На k-й итерации глобального метода для улучшения управления u(k)(t), t G [0,T], k > 0, в задаче (5), (6), (8) выполняются следующие шаги:
1) интегрирование системы (5) при управлении u(k) для расчета траектории z(k);
2) вычисление решения сопряженной системы (7) на процессе (z(k),u(k));
3) интегрирование системы (5) с отображением (t, z), t G [0, T], в плане расчета траектории z(fc+1);
4) построение нового приближения u(fc+1)(t) = M(fc)(t,z(fc+1)(t)), t G [0,T].
В алгоритме изучается поведение знака производной функции переключения, следуя формулам из статьи [11].
Линейность задачи по u является положительным моментом для разрешения операции максимизации функции Понтрягина. С другой стороны, интегрирование системы (5) с зависимостью (9) имеет известные особенности [11]. В этой связи на практике известен подход
с регуляризацией целевого функционала, следуя А.Н. Тихонову [13]. Вводится в рассмотрение вспомогательный целевой критерий
т
2
(10) Je(z,u) = I(z, u) + в/u2(t)dt —— inf, в > 0.
0
Функция Понтрягина He (t, p, z, u) = H(t, p, z, u) — eu2. Дифференцируя эту функцию по u и рассматривая p(k)(t), приходим к зависимости
(11) ue(t,z) = ^K(t,p(k)(t),z), t G [0,T], в> 0.
2в
Глобальный метод применительно к вспомогательной задаче (10), (5), (6) не требует учета условно введенного ограничения (8). В отличие от стандартного глобального алгоритма, необходим подбор достаточно малого значения в > 0 в данной модификации. Вопросы реализации подхода, основанного на (10), (11), относительно билинейных систем в рамках глобального метода изучались также в работе [12].
3. Вычислительные эксперименты
Проведена программная реализация глобального метода (на языках Ма^аЬ, Еогігап). Расчеты проведены для различных условных значений фа, £. Представим сравнительные результаты расчетов с условным набором
^ = (1/^2,1/^2,0,0)Т, аі = 0.6, 6і = -0.3, а2 = 0.1, 62 = ^0.54, ш = I =1, Т =1.5.
Приближение и(0)(і), і Є [0,Т], представляет собой решение краевой задачи [8]
U(0) (t) = 4y (w2 + u2(t)),
1
Y =
4Tw
( u(0) (T )\ / u(0) (0) \
arctg f ------------j — arctg f --------------j
/(0) (0) = _,,(°)(
2TWarctg (, u(0)(0) = —u(0)(T).
Формирование начального приближения таким образом исходит из адиабатического приближения в квантовой механике [5]. Общее ре-
/ и(0)(і)\
шение агсіг -------- = 47ші + С, частное решение
V ш )
(12) и(0) (і) = шtg (27ш(2і — Т)), і Є [0, Т].
Функция и(0)(і), і Є [0, Т], формируется численно по формуле (12) в компьютерной программе.
Полагаем є = 0.001. Интегрирование систем методом Рунге-Кут-та 4-го порядка.
1. Глобальный метод с применением формулы для управления при особом реж.име [11]. Считаем и(0)(0) = —30 в формуле
(12) для 7. Задаем V = 30 в условии (8).
На начальном приближении и(0) (і) имеем значение I(0) = 0.7681. В таблице 1 приведены результаты последовательных улучшений по первым итерациям (ГлМ — глобальный метод, ГрМ — градиентный метод). Рис. 1 иллюстрирует высокую сравнительную эффективность глобального метода. В отличие от градиентного метода, глобальный метод не зависит от алгоритмического параметра. Глобальный метод обеспечивает быструю релаксацию: на 9-й итерации получено значение I = 0.000952, что составляет 0.124 % от I(0). Сравнение графиков управлений и(1), и(9), полученных за первую итерацию и в результате серии итераций соответственно, показывает, что участок [0,іі), на котором управление принимает значение —V, важен в плане реализации выхода на особый режим.
Таблица 1.
Номер к I(г(й),и(й)) по ГлМ I(г(й),и(й)) по ГрМ
0 0.7681 0.7681
1 0.1401 0.6911
2 0.0040 0.6107
3 0.0021 0.5421
4 0.0015 0.4913
На рис. 2 приведены графики управлений и(0) (і) и и (9)(і). В точке іі = 0.0667 происходит переключение с граничного управления и(9)(і) = —30 на внутреннее и(9)(і) = и (9)(і), где функция и(9) (і)
П Я
ГрМ
0.6 \ ^_
П А
\
П 9
\ 1ЛМ
Л - \
1 2 3 4 5 *
Рис. 1.
30 20 10 0 -10 -20 -30
0 , 0.25 0.5 0.75 1 1.25 1.5 ,
Рис. 2.
определяется согласно теории [ ]. Участку [0, ^ 1) соответствуют до-
статочно большие по модулю отрицательные значения функции переключения к(9)(г) = К(р(8)(г),г(9) (г)).
На рис. 3 показаны графики функций г(9)(г), ] = 1, 4, соответствующих управлению и(9)(г). Графики иллюстрируют достаточно быстрое изменение траектории на участке [0, г 1), где применяется управление, равное — V. Это необходимо для выхода на особый режим при заданном начальном состоянии.
Сравнительные эксперименты с привлечением градиентного метода [10] показывают, что глобальный метод надежнее, так как, во-первых, учитывает специфику особого оптимального управления и,
Рис. 3.
во-вторых, не требует настройки алгоритмического параметра (в не-регуляризованной, исходной формулировке).
2. Глобальный метод с регуляризацией по управлению. Полагаем значение 0) = —30 в формуле для 7 в (12). Рассматриваемая модификация глобального метода не требует условного ограничения (8); тем не менее, в контексте с с изложенными выше результатами расчетов полагаем V = 30 в ограничении (8).
Рассматриваем сначала /3 = 0.0055. В таблице 2 показаны последовательные улучшения по итерациям, а на рис. 4 — график управления г№(і).
Регуляризация требует настройки параметра, причем проигрывает стандартному глобальному методу и по шагам уменьшения функционала /, что обусловлено влиянием добавочного слагаемого в функционале Iі3. Положительной чертой модификации является непрерывность функции «(*>(*,£) по г.
Далее, рассматриваем меньшее значение /3 = 0.002. В таблице 3 показаны последовательные улучшения по итерациям. Рис. 5 иллюстрирует учет ограничения (8) при данном значении /3.
Рис. 4.
Рис. 5.
Рис. 4, 5 показывают наличие участка управления с достаточно быстрым изменением в некоторой полуокрестности [0,#і), что обусловлено учетом начального условия г(0) = гв-
Таблица 2.
Номер к Значение I(z(fc),u(fc)) Номер к Значение I(z(fc) ,u(k))
0 0.7681 5 0.0027
1 0.0614 6 0.0017
2 0.0242 7 0.0014
3 0.0119 8 0.0013
4 0.0057 9 0.0012
Таблица 3.
Номер к Значение I(z(fc),u(fc)) Номер к Значение I(z(fc) ,u(k))
0 0.7681 5 0.0027
1 0.0876 6 0.0019
2 0.0122 7 0.0015
3 0.0068 8 0.0013
4 0.0042 9 0.0012
Заключение
В статье для задачи оптимального управления квантовой системой на модели Ландау-Зинера проведены реализации стандартного метода глобального улучшения [2, 11], а также модификации с регуляризацией функционала по управлению. Основным результатом статьи является демонстрация высокой эффективности глобального метода применительно к модели, характерной для проблем оптимального управления системами спинов, описываемыми уравнением Шредингера.
Представляется перспективным распространение результатов работы на задачи об управлении цепочками спинов с рассмотрением различных прикладных аспектов. В связи с большими размерностями реальных прикладных задач представляется перспективным привлечение суперкомпьютерных технологий (MPI, T+—+ и т.п.).
Список литературы
[1] Кротов В. Ф., Фельдман Н. Н. Итерационный метод решения задач оптимального управления // Изв. АН СССР. Техн. киберн., 1983, № 2, с. 160—168. ?[]
[2] Krotov V. F. Global methods in optimal control theory. New York : Marcel Dekker, 1996.— 408 p. ?[], 3
[3] Кротов В. Ф., Гурман В. И. Методы и задачи оптимального управления. М. : Наука, 1973.— 448 с. ?[]
[4] Васильев О. В. Лекции по методам оптимизации. Иркутск : Изд-во Иркутск. ун-та, 1994.— 344 с. ?[]
[5] Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. Т. 3: Квантовая механика (нерелятивистская теория). 6-е изд. M., 2008.— 700 с. ?[], 1, 3
[6] Валиев К. А. Квантовые компьютеры и квантовые вычисления // Успехи физ. наук, 2005. 175, с. 3—39. ?[]
[7] Казаков В. А., Кротов В. Ф. Оптимальное управление взаимодействием света и вещества // Автоматика и телемеханика, 1987, № 4, с. 9—15. ?[]
[8] Caneva T., Murphy M., Calarco T., eta. Optimal control at the quantum speed limit // Physkal Review Lett., 2009. 103. ?[], 1, 1, 3
[9] Кротов В. Ф. Об оптимизации управления квантовыми системами // Доклады РАН, 2008. 423, № 3, с. 316-319. ?[], 1, 2
[10] Кротов В. Ф. Управление квантовыми системами и некоторые идеи теории оптимального управления // Автоматика и телемеханика, 2009, № 3, с. 15-23. ?2, 2, 3
[11] Батурина О. В., Булатов А. В., Кротов В. Ф. К оптимизации линейных систем с управляемыми коэффициентами // Автоматика и телемеханика, 2011 (в печати). ?[], 2, 2, 3, 3, 3
[12] Батурина О. В. Билинейные динамические системы: исследование итеративных методов оптимизации // Проблемы управления, 2010, № 5, с. 2227. ?[], 2
[13] Тихонов А. Н. О методах регуляризации задач оптимального управления // Докл. АН СССР, 1965. 162, № 4, с. 763-765. ?2
O. V. Baturina, O. V. Morzhin. Quantum system control optimization for the Landau-Ziner model.
Abstract. The article is devoted to computer experiments according to application of the global (nonlocal) improvement method by V.F. Krotov for investigation of an optimal control problem for quantum system with Schrodinger’s equation and hamiltonian in the form of Landau-Ziner.
Key Words and Phrases: quantum systems, global improvements.
Образец ссылки на статью:
О. В. Батурина, О. В. Моржин. Оптимизация управления квантовой системой на модели Ландау-Зинера // Программные системы: теория и приложения : электрон. научн. журн. 2011. № 1(5), с. 51-61. URL: http://psta.psiras .ru/read/psta2011_1_51-61.pdf
