ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ В КОГЕРЕНТНО ВОЗБУЖДАЕМЫХ НАНОСТРУКТУРАХ И МЕТАМАТЕРИАЛАХ
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
ОБЗОРНАЯ СТАТЬЯ
УДК 538.9
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ В КОГЕРЕНТНО ВОЗБУЖДАЕМЫХ НАНОСТРУКТУРАХ
И МЕТАМАТЕРИАЛАХ
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов, Р.Е. Носков, Ю.С. Кившарь
Розанов Николай Николаевич, членкорреспондент Российской академии наук, доктор физико-математических наук, профессор, начальник отдела ФГУП «НПК «ГОИ им. С.И. Вавилова», заведующий кафедрой Оптики лазеров СПбНИУ ИТМО. Лауреат премии им. Д.С. Рождественского РАН (2007). Область научных интересов – физическая и нелинейная оптика, лазерная физика, физика солитонов.
Высотина Нина Вениаминовна, научный сотрудник ФГУП «НПК «ГОИ им. С.И. Вавилова». Об-
ласть научных интересов – вычислительная математика, численное моделирование.
Шацев Анатолий Натанович, кандидат физико-математических наук, старший научный сотрудник ФГУП «НПК «ГОИ им. С.И. Вавилова». Область научных интересов – нелинейная оптика, лазерная физика, физика солитонов.
Десятников Антон Сергеевич, кандидат физико-математических
наук, старший научный сотрудник Австралийского национального университета.
Шадривов Илья Владимирович, кандидат физико-математических наук, старший научный сотрудник Австралийского национального университета, стипендиат Королевы Елизаветы II.
Кившарь Юрий Семенович, академик Австралийской академии наук, профессор, руководитель отдела, ведущий ученый, научный руководитель лаборатории Австралийского национального университета, обладатель мегагранта Министерства образования и науки Российской Федерации.
Носков Роман Евгеньевич, кандидат физико-математических наук, старший научный сотрудник Санкт-Петербургского национального исследовательского университета информационных технологий, механики и оптики. Область научных интересов – оптика наноструктурированных композитных сред на основе металлических наночастиц, левые среды, электромагнитная невидимость.
Проведен сравнительный анализ локализованных структур – волн переключения и диссипативных солитонов – в трех дискретных нелинейных системах, резонансно возбуждаемых когерентным поддерживающим излучением: цепочке молекул, цепочке металлических наношаров и решетке разомкнутых кольцевых резонаторов. Ключевые слова: диссипативные солитоны, дискретные солитоны, волны переключения.
Введение
Диссипативные солитоны, представляющие локализованные за счет баланса притока и оттока энергии структуры волн, к настоящему времени более полно изучены в таких широкоапертурных оптических системах, как нелинейные интерферометры с когерентным поддерживающим излучением и лазеры с насыщающимся поглощением [1]. В этих континуальных системах продемонстрированы одномер-
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
1
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
ные и двумерные пространственные солитоны и предложены их приложения для обработки информации. Однако размеры этих солитонов составляют не менее 10 мкм, что не делает такие системы конкурентноспособными по сравнению с традиционной электроникой. Кроме того, в подобных системах сложно формировать трехмерные, в особенности топологические солитоны.
Указанные ограничения могут быть преодолены при использовании периодических наноструктур и метаматериалов, которые представляют собой набор слабо связанных нелинейных элементов, возбуждаемых внешним когерентным излучением [2–4] (рис. 1). В этой статье мы рассмотрим три родственные системы, в которых могут быть сформированы дискретные диссипативные структуры.
Модель 1 – это цепочка молекул с трехуровневой схемой электронных переходов. Внешнее излучение резонансно возбуждает молекулы, под его действием молекулы излучают, подобно диполям, и это излучение частично перехватывается соседними молекулами (диполь-дипольное взаимодействие молекул). К настоящему времени в расчетах продемонстрированы солитоны с размерами в несколько нанометров и возможность управлять их положением в цепочке и движением с помощью наклонного падения внешнего излучения [1, 5]. Перспективны для реализации таких систем J-агрегаты цианиновых красителей, обладающие коллективным (экситонным) механизмом их резонансного возбуждения лазерным излучением. В то же время еще предстоит разработка более полной теории этих принципиально квантовых объектов с выходом за рамки одночастичного приближения и учетом межмолекулярных корреляций.
у
z
E0
k
x
Рис. 1. Одномерная дискретная система. Цепочка нелинейных элементов, изображенных кружками, расположена вдоль оси x. E0 – вектор напряженности электрического поля поддерживающего излучения;
волновой вектор k падающей волны лежит в плоскости xz; φ – угол между волновым вектором k и осью z;
θ – угол поворота вектора поляризации поля относительно оси y
m–1,n
n–1 n n+1
m,n–1 m,n m,n+1
m+1,n
Рис. 2. Цепочки (слева) и двумерные решетки (справа) разомкнутых кольцевых резонаторов с различной ориентацией
Модель 2 близка к предыдущей, но в ней молекулы заменяются металлическими наношарами, погруженными в диэлектрик [6, 7]. В отличие от молекулярных J-агрегатов с фиксированным направлением дипольного момента, наношары обладают изотропным откликом, так что здесь добавляется поляризационная степень свободы, и солитоны становятся в общем случае векторными. К настоящему времени в литературе для моделей 1 и 2 изучены только одномерные системы, и рассмотрение хотя бы двумерных
2 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
схем представляет важную задачу дальнейших исследований. В то же время в модели 3 – цепочке или решетке возбуждаемых когерентным излучением разомкнутых кольцевых резонаторов (рис. 2) с нелинейными элементами, помещенными в разъемы [8–16] – возможна реализация дискретных диссипативных солитонов любой геометрической размерности. В моделях 1 и 2 аналогичен вид связи элементов (диполь-дипольное взаимодействие), но различается вид нелинейности. Последний совпадает (кубическая нелинейность) для моделей 2 и 3. Ниже мы проанализируем общие и специфические свойства солитонов в трех указанных моделях, уделяя большее внимание менее изученным в литературе вариантам – модели 2 для одномерной геометрии и модели 3 для двумерной и трехмерной геометрии.
Управляющие уравнения
Модель 1. J-агрегат моделируется линейной цепочкой N трехуровневых молекул, взаимодейст-
вующих между собой через электромагнитное поле. Два нижних уровня формируют оптический переход,
частота которого близка к лазерной частоте. Третий молекулярный уровень служит для описания анни-
гиляции двух возбуждений на соседних молекулах: одно из возбуждений дезактивируется, а другое пе-
реходит на третий уровень с последующей релаксацией либо обратно на второй, либо на первый уровень.
В результате соответственно либо одно, либо оба возбуждения исчезают из системы. Частота перехода из
основного состояния (1) в третье (3) предполагается приблизительно равной удвоенной частоте перехода
1 2 . Расстояние между молекулами d много меньше длины волны излучения. Частота лазерного излу-
чения ω0, поляризация излучения линейная. В подходе одночастичной матрицы плотности управляющие уравнения для элементов матрицы
плотности ρik, i, k = 1, 2, …, N, в приближении медленно меняющихся амплитуд имеют вид (см. [1, 5] и приведенные там ссылки)
1(1k )
1 4
N l 1
(lk
ilk )Rl Rk*
i
4
Rk*
c.c.
s
(k 22
)
[2(2k
1)
(k 22
1)
]
313(3k
)
(k
21 22
)
(k 22
)
,
2(2k )
1 4
N l 1
(lk
ilk )Rl Rk*
i
4
Rk*
c.c.
2s
(k 22
)
[2(2k
1)
(k 22
1)
]
32
(k 33
)
(k
21 22
)
(k 22
)
,
3(3k)
s
(k 22
)
[2(2k
1)
(k 22
1)
]
31 32
(k ) 33
,
Rk
( i0 )Rk
N
(lk
ilk )Rl [2(2k )
(k 11
)
]
i [2(2k )
(k 11
)
]
s
Rk
[2(2k
1)
(k 22
1)
]
.
l 1
Точка в левых частях уравнений означает временную производную; с.с. – знак комплексного со-
пряжения; 0 – значение частотной расстройки для молекул в цепочке; Rk 22(1k) exp(i0t) . Матрицы lk
и lk ( l k ) описывают запаздывающее диполь-дипольное взаимодействие молекул в агрегате:
lk
2 d 3
k0
d
cos
k0d l k lk 2
sin
k0d l k lk 3
1 3cos2
k0d 2 sin
k0d l k lk
sin2 , l k,
lk
2 d 3
cos
k0d l k l k 3
sin k0d
k0d l k lk 2
1 3cos2
k0d 2 cos
k0d l k lk
sin2 , l k,
(1)
ll 0, ll 0.
Волновое число k0 0 c , где c – скорость света в вакууме. Величина ( / )Ei – частота Раби, пропорциональная комплексной амплитуде поля, так что вещественная напряженность электрическо-
го поля в падающей волне
E Re[Ei (t) exp(i0t)] . Другими параметрами задачи являются матричный элемент дипольного момента перехода μ, угол
между дипольным моментом и осью цепочки θ, релаксационные постоянные γ и Г, параметр экситон-
экситонной аннигиляции s . Принятый подход оправдан, если уровень возбуждения молекул слабый:
1
ρ(k 11
)
1
(диагональные
элементы
матрицы
плотности
пропорциональны
населенностям
соответст-
вующих уровней).
В линейных агрегатах взаимодействие с излучением молекул у краев цепочки существенно отли-
чается от случая молекул в середине цепочки, главным образом, из-за эффекта двухэкситонной анниги-
ляции (параметр s ). В типичных условиях эффекты краев захватывают несколько десятков молекул, поэтому для средней части достаточно длинной цепочки от краевых эффектов можно отвлечься. Чтобы
исключить их в расчетах, можно представить, что цепочка замкнута в кольцо (периодические граничные
условия). Хотя молекулярные агрегаты с кольцевой геометрией и известны, молекулы в них неодинаково
взаимодействуют с излучением с линейной поляризацией (частота Раби пропорциональна косинусу угла
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
3
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
между направлениями дипольного момента молекулы и поляризации излучения). Поэтому для обеспечения одинаковых условий для различных молекул пришлось бы перейти к их взаимодействию с излучением, обладающим, например, круговой поляризацией.
В модели 2 управляющие уравнения имеют вид [6, 7]
i
dPn d
i Pn 2
Pn
G n,m
Pm
En ,
mn
i
dPn d
i Pn 2
Pn Gn,m Pm En ,
mn
(2)
где Pn 2 Pn 2 Pn 2 , верхние индексы указывают составляющие векторов отклика шаров (Р) и под-
держивающего излучения (Е) по отношению к оси цепочки (рис. 1), параметры γ, Ω, k0 и d имеют тот же смысл, что и в модели 1, а коэффициенты связи элементов аналогичны (1):
G n,m
2
k0d
2
ik0d nm
n
1 m
2
exp
ik0d n m nm
,
G n,m
ik0d nm
n
1 m
2
exp
ik0d n m nm
.
Параметр η определяется отношением диаметра шаров и шага решетки, а также диэлектрической
проницаемостью шаров и матрицы-диэлектрика [6, 7].
Управляющие уравнения для модели 3 формулируются для медленно меняющихся амплитуд электрического тока или намагниченности отдельных резонаторов [8]:
i
dn dt
(2 i | n
|2 )n
Cn
.
Здесь время t нормировано на период собственных колебаний изолированного резонатора, Ω – от-
стройка собственной частоты колебаний от частоты когерентного поддерживающего излучения, норми-
рованная на собственную частоту, γ – затухание резонатора с той же нормировкой, 1 – знак коэф-
фициента керровской нелинейности. Коэффициент κ определяет силу связи ближайших резонаторов, и
Cn 2[(nx , ny , nz 1) (nx , ny , nz 1)] [(nx 1, ny , nz ) (nx 1, ny , nz )]
[(nx , ny 1, nz ) (nx , ny 1, nz )], n (nx , ny , nz ).
Целочисленные индексы nx , ny , nz определяют положение резонатора в решетке. Управляющие
уравнения для каждой модели дополняются начальными условиями, роль которых особо важна в услови-
ях бистабильности или тристабильности (см. ниже).
Одномерные локализованные структуры (модель 2)
Квазиоднородные распределения. Для идеальной цепочки бесконечной длины с одинаковым стационарным состоянием всех элементов (пространственно-однородные стационарные состояния) из (2)
при Pn P0 и нулевых производных получаем алгебраическое уравнение 5-й степени:
P0
2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
E0
2
cos
2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
sin 2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
,
где
Aj j3 ik0d j2 k0d 2 j1 exp ik0d j ,
Aj 2 j3 ik0d j2 exp ik0d j .
При 0, / 2 степень полинома снижается до третьей, тогда одному и тому же значению интен-
сивности поддерживающего излучения могут отвечать три значения отклика. С учетом заведомой неустойчивости промежуточного значения это означает бистабильность (рис. 3, а). Этот вариант характерен
4 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
также для скалярных моделей 1 и 3, тогда как тристабильность становится возможной для векторной модели 2 (рис. 3, б).
0,2 0,3
0,16 0,12 0,2
|P0,|| |2 |P0 |2
0,08 0,1 0,04
00
0 0,0001 0,0002 0,0003 0,0004
0 0,0002 0,0004 0,0006 0,0008
|Е0,|| |2 а
|Е0 |2 б
Рис. 3. Зависимость поперечной и продольной поляризации P0 2 , P0 2 от интенсивности внешнего поля
I E0 2 , E0 2 для бесконечной цепочки; перпендикулярная поляризация падающей волны (θ = 0) –
сплошная кривая, параллельная поляризация (θ = /2) – пунктирная кривая, 0,1 (а).
Зависимость P0 2 от E0 2 при 0, 2 , 0, 2 для бесконечной цепочки (б). Устойчивость для
конечной цепочки проверялась численно. Верхняя ветвь устойчива целиком, нижняя ветвь устойчива при
интенсивностях падающей волны от 0 до ~0,0004, остальные ветви неустойчивы (расчет при N 100 )
|Pn|2
0,16 0,01 0,008
0,12 0,006 0,08 0,004
|Pn|||2
0,04 0,002
|Pn|2
00
0
0,1 0,08
20 40 60 80 100 n
а
0
0,1 0,08
20 40 60 80 100 n
б
0,06 0,06
|Pn|||2
0,04 0,04
0,02 0,02
00
0 20 40 60 80 100 n
0 20 40 60 80 100 n
вг
Рис. 4. Интенсивность перпендикулярной (а, в) и параллельной (б, г) составляющих отклика
при E0 2 0, 000088 и 0,12 . Для (в, г) начальные условия нулевые, для (а, б) ненулевые
0,2 0,2
0,16 0,16 0,12 0,12
|P|2
|P|2
0,08 0,04 b a
0,08 0,04 b a
00
0 40 80 120 160 200
0 40 80 120 160 200
nn
аб Рис. 5. Профили неподвижной (a) и движущейся (b) волны переключения. I = 0,55×10–4 , скорость V < 0
(область 2 рис. 7) (а); I = 0,88×10–4,V > 0 (область 5 рис. 7) (б).
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
5
|P|2 w
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
Для цепочки конечной длины вследствие краевых эффектов однородные распределения заменяются на неоднородные. При этом следует говорить уже о пространственной бистабильности (бистабильности пространственных структур, устанавливающихся при различных начальных условиях) [1]. Еще одно важное условие реализации режима состоит в требовании его устойчивости относительно малых возмущений. Для бесконечной цепочки линейный анализ устойчивости стационарных однородных распределений выполняется стандартным образом [6, 7], в других случаях обычно необходимы численные расчеты. В результате развития неустойчивости в цепочке формируются и устанавливаются квазипериодические структуры. Интересной особенностью, известной и для континуальных систем [1], является возможность установления при одних и тех же параметрах системы, но при различных начальных условиях, либо квазиоднородной (рис. 4, а, б), либо квазипериодической (рис. 4, в, г) структуры.
Волны переключения в скалярной системе. Как и в континуальных системах, внутри интервала бистабильности квазиоднородных распределений существуют волны переключения (ВП) (другие названия – доменные стенки и кинки). При удалении от фронта волны переключения в одну или другую сторону распределение асимптотически приближается к одному из возможных квазиоднородных распределений, различающихся по разные стороны от фронта. Фронт волны может быть неподвижным (неподвижная волна переключения, рис. 5, кривая а, или движущимся (движущаяся волна переключения, рис. 5, кривая б). Если для континуальных систем скорость волны переключения обращается в нуль при одном («максвелловском») значении интенсивности поддерживающего излучения [1], то в дискретных системах это происходит на целом интервале интенсивности. Общей особенностью дискретных волн переключения служит то, что форма фронта и скорость движущейся волны меняются со временем квазипериодически (рис. 6). В связи с этим имеет смысл говорить о средней скорости таких волн.
0,3 12
0,2 8
0,1 4
00 500 600 700 800 900 1000
t Рис. 6. Временное изменение максимальной интенсивности отклика на фронте (штриховая кривая, левая шкала) и ширины фронта w (сплошная кривая, правая шкала, в единицах расстояния между элементами)
движущейся волны переключения в условиях рис. 5, б (кривая b, I = 0,88 10–4). Ширина фронта w
определялась как длина промежутка, на котором интенсивность волны спадает от максимального значения не менее чем в 100 раз
При изменении интенсивности поддерживающего излучения скорость волны переключения меняется, в том числе по знаку. При этом мы считаем скорость V положительной, если со временем при передвижении фронта увеличивается число элементов с большей амплитудой отклика (рис. 5, б, в отличие от рис. 5, а, где V < 0). Зависимость скорости волны переключения V от интенсивности падающего излучения I приведена на рис. 7. Для удобства диапазон изменения I разделен на 6 областей. Хотя в области 1 (0,0381×10–4 < I < 0,46×10–4) и имеется бистабильность однородных состояний бесконечной цепочки, повидимому, область притяжения верхнего квазиоднородного состояния столь узка, что практически оно неустойчиво (из-за неизбежных в расчете возмущений). При этом исходное ступенчатое распределение быстро сваливается в нижнее квазиоднородное состояние, так что волны переключения отсутствуют. В области 2 (0,46×10–4 < I < 0,56×10–4) существует как неподвижная, так и отрицательная (V < 0) волны переключения (рис. 5, а). В области 3 (0,56×10–4 < I < 0,836×10–4) существует только неподвижная волна переключения. В области 4 (0,836×10–4 < I < 0,9216×10–4) существует как неподвижная, так и движущаяся с положительной скоростью (V > 0) волна переключения (рис. 5, б). При I > 0,9216×10–4 происходит скачкообразный переход из режима с V = 0 в режим с V > 0. В области 5 (0,9216×10–4 < I < 2,045×10–4) существуют только волны переключения с положительной скоростью. Заметим, что приблизительно при I > 1,4×10–4 уже трудно получить одиночную волну переключения, так как из-за краевых эффектов неизбежно порождается встречная волна.
Все эти области 1–5 находятся внутри области пространственной бистабильности квазиоднородных состояний. Однако и вне этой области, но внутри области бистабильности для бесконечной цепочки, имеются метастабильные волны переключения. Дело в том, что дестабилизирующее действие краевых
6 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
эффектов происходит не мгновенно, и если время их развития превышает время распространения фронта по конечной цепочке, то неустойчивость будет существенно подавлена. Поэтому при 2,045×10–4 < I < < 2,9×10–4 (область 6) имеются метастабильные волны переключения с положительной скоростью. Наконец, при I > 2,9×10–4 даже метастабильных (со значительным временем жизни) волн переключения нет, переход на исходной ступеньке из нижнего состояния в верхнее происходит практически сразу, скачком.
0,08 1 2 3 4
5
6
0,08 1 2
3 45
0,06 0,06
0,04 0,04
0,02 0,02
V/c |Pn|2
V/c |Pn|||2
00
–0,02 –0,04 –0,06
12 E|E00|21100–44
3 –0,02 –0,04 –0,06
0,6 0,8 |EE00|22101–044
1
–0,08
–0,08
–0,1 –0,1
–0,12
–0,12
–0,14
–0,14
аб
Рис. 7. Зависимость скорости волны переключения V c от интенсивности поддерживающего излучения
I = |E0|2 (а). То же, укрупненно (б)
Векторные волны переключения ( 0 ). Если внешнее излучение содержит обе поляризационные компоненты, то, как указывалось выше, возникают следующие особенности. Во-первых, в соответствии с рис. 3, б, может иметь место мультистабильность (тристабильность) квазиоднородных состояний, из-за чего возможны новые типы волны переключения и, как результат их столкновения, новые типы солитонов. Во-вторых, развитие модуляционной неустойчивости некоторых из квазиоднородных состояний может при нелинейном ограничении приводить к бистабильности между квазиоднородным и пространственно промодулированным состоянием и соответственно к волнам переключения между такими состояниями. В этом разделе мы ограничимся изучением только второй возможности, выбирая параметры так, чтобы они отвечали бистабильности квазиоднородных распределений.
При 0,12 и E0 2 0, 0001 квазиоднородное состояние на верхней ветви устойчиво, а на ниж-
ней – модуляционно неустойчиво, со временем оно превращается в стационарное пространственно промодулированное состояние. Как показывает рис. 8, исходное ступенчатое распределение возбуждает волну переключения, фронт которой со временем надвигается на промодулированное состояние. Для цепочки конечной длины фронт останавливается на некотором расстоянии от края цепочки.
0,2 0,08
0,16 3
0,12 0,08
12
0,04
0,06 0,04 0,02 1 2
3
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 8. Интенсивность перпендикулярной (a) и параллельной (б) составляющей поля при 0,12 и
E0 2 0, 0001 . Исходное распределение – ступенчатое. Кривые 1, 2, 3 – соответственно распределения полей в моменты времени t 500 , 2000 и 3500 (установившаяся структура)
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
7
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
При меньшем значении интенсивности E0 2 0, 000088 и 0,12 фронт между квазиоднород-
ным и промодулированным состояниями неподвижен даже вдали от краев цепочки. Устанавливается стационарное распределение (рис. 9), т.е. ступенька практически сразу останавливается.
Наконец, при E0 2 0, 0001 и 0,1 и нижнее, и верхнее квазиоднородные состояния устойчи-
вы. В этом случае волны переключения – того же типа, как в скалярном случае ( 0 ). Движение фронта ВП не останавливается краевыми эффектами (рис. 10).
0,16 0,04
0,12 0,03
|Pn|||2
|Pn|2
0,08 0,02
0,04 0,01
0 0 20 40 60 80 100 n
0 0 20 40 60 80 n
аб
Рис. 9. Интенсивность перпендикулярной (a) и параллельной (б) составляющей отклика
при E0 2 0, 000088 и 0,12
100
0,2 0,03
|Pn|2
0,16
0,12
0,08
12 3
4
0,04
|Pn|||2
0,02
0,01
12 3
4
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 10. Движущаяся волна переключения в условиях классической бистабильности
для перпендикулярной (a) и параллельной (б) составляющей отклика. Кривые 1, 2, 3, 4 – соответственно
распределения полей в моменты времени t 500, 2500, 4500, 7000 ; E0 2 0, 0001 , 0,1
Скалярные дискретные диссипативные солитоны (модель 2, θ = 0)
Важным сценарием формирования диссипативных солитонов служит столкновение волн пере-
ключения. Поскольку в рассматриваемой системе возможны такие волны различных типов, то число
сценариев увеличивается.
Столкновение волн переключения с нулевой и положительной скоростью ( I 0,88 104 , область 4
рис. 7). На рис. 11, а, справа движется волна навстречу неподвижной волне в левой половине цепочки.
0,25 0,2
0,2 0,16 5
|Pn|2 |Pn|2
0,15 0,1 4 3 2 1
0,05
0,12 0,08 0,04
1234 4321
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 11. Динамика столкновения волн переключения. Кривые 1, 2, 3, 4 – соответственно распределения
интенсивности отклика в моменты времени t 1000 , 2000, 3000 и 4000, I 0,88 104 (а),
кривые 1, 2, 3, 4, 5 – соответственно распределения интенсивности отклика в моменты времени
t 1000 , 1500, 2000, 2500 и 5000, I 0,95 104 (б)
Движущаяся справа волна останавливается с образованием темного солитона, включающего 15
элементов. Левая «неподвижная» волна в самом конце столкновения сдвигается и останавливается. Ана-
8 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
|Pn|2 |Pn|2
логичные, но более узкие темные солитоны формируются при бóльших интенсивностях I, вплоть до
I 0,94 104 , когда устанавливается солитон с шириной, отвечающей 6 элементам. Но уже при
I 0,95 104 сталкивающиеся волны переключения не останавливаются, а аннигилируют, результи-
рующее поле стационарно и имеет квазиоднородное распределение (кривая 5 на рис. 11, б). Столкновение волн переключения с нулевой и отрицательной скоростью (область 2 рис. 7). Ана-
логично предыдущему сценарию, в результате взаимодействия двух волн переключения образуются солитоны, но теперь они светлые. На рис. 12, а, приведены профили результирующих стационарных распределений поля для различных значений интенсивности поддерживающего излучения I. Самый узкий
солитон, соответствующий I 0,5510 104 , имеет ширину 2 элемента. Динамика формирования таких
солитонов показана на рис. 11, б. Аннигиляцию волн переключения с нулевой и отрицательной скоростью иллюстрирует рис. 13.
0,16 0,16
0,12 0,12
0,08 0,08
0,04
0,04
43
21
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 12. Установившиеся солитоны при различных значениях интенсивности: пунктирная, сплошная, пунк-
тирная жирная, сплошная жирная кривые соответствуют значениям
E0 2 104 0,555; 0,552; 0,5513; 0,5510 (а); динамика образования солитона с шириной w 2 при столк-
новении волны переключения, движущейся налево (кривые 1–4 отвечают моментам времени t = 0, 500,
1500, 3000) с неподвижной волной переключения; I 0, 0000551 (б).
0,16
0,12
|Pn|2
0,08 3 2 1 0,04 4
0 0 40 80 120 160 200 n
Рис. 13. Левый фронт практически неподвижен, а правый фронт движется налево (кривые 1–4 соответствуют моментам t = 200, 1600, 2800, 4000), где сталкивается с неподвижным фронтом, и фронты
схлопываются; I = |E0|2 =0,5476×10–4
0,2
0,16
0,12
|Pn|2
0,08
0,04
0
0 20 40 60 80 100
n Рис. 14. Одностороннее расширение области возбуждения цепочки между фронтами неподвижной ВП (слева) и ВП, движущейся с положительной скоростью. Моменты (слева направо) t = 500, 1000, 1500, 2000;
I = |E0|2 =0,8491×10–4 (область 3 рис. 7)
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
9
|Pn|2 |Pn|2 |Pn|2 |Pn|2
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
Если фронт движущейся волны переключения не приближается к фронту неподвижной волны, а удаляется от него, и начальное расстояние между фронтами не слишком мало, то возникает расширяющаяся со временем область повышенного возбуждения (рис. 14). При значительном расстоянии между фронтами эти две волны переключения практически не взаимодействуют друг с другом (расширение в одну сторону).
Диссипативные солитоны могут формироваться и по другим сценариям, например, при подаче на цепочку дополнительного «записывающего» узкого пучка-импульса или нескольких пучков-импульсов. Примеры установившихся структур приведены на рис. 15. При выходе за границы области устойчивости, например, для модели 1, характерно появление пульсирующих солитонов. Рис. 16 демонстрирует, также для модели 1, возможность управляемого передвижения устойчивых солитонов по цепочке.
0,12
0,16 0,12
0,12
0,08
0,12 0,08
0,08
0,08
0,04
0,04 0,04
0,04
0
00
0
0 20 40 60 80 100 n
0 20 40 60 80 100 0 20 40 60 80 100 nn
0 20 40n 60 80 100
Рис. 15. Примеры узких светлых солитонов; I = |E0|2 =0,5329×10–4 (область 1 рис. 7)
11 11 11 11 0,98 0,98 0,98 0,98
0,94 0,94 0,94 0,94
0,9 0
0,9
100 200 300 0 n а
0,9 100 200 300 0
n
б
100 200 n
0,9 300 0
в
100 200 300 n г
Рис. 16 (модель 1). Управление движением солитона при наклонном падении поддерживающего излучения. Профили населенности основного состояния молекул в цепочке в различные моменты времени t: t = 0, угол падения θ = 0, солитон неподвижен (а); t = 50, угол падения θ = 0,1, солитон движется направо (б); t = 300, угол падения θ = 0,1, солитон останавливается у правого края цепочки (в); t = 700,
угол падения θ = – 0,1, солитон двигается налево и останавливается у левого края цепочки (г)
Двумерные локализованные структуры (модель 3)
Более удобна для анализа двумерных решеток схема разомкнутых кольцевых резонаторов (модель 3, рис. 2) В этом случае схема существенно анизотропна, так как взаимодействие резонаторов зависит от их относительной ориентации. При фиксированной ориентации устойчивыми оказываются только сравнительно узкие двумерные солитоны (рис. 17). Широкие же исходные структуры обычно распадаются на набор узких солитонов вследствие модуляционной неустойчивости.
0,54 0,55
0,34 0,35
0,14 0,15
Рис. 17. Два примера двумерных солитонов в модели 3; шкала для намагниченности | ψ |2
указана справа от рисунков, положение резонаторов задается целыми индексами nx и nz
10 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
Трехмерные локализованные структуры (модель 3)
Трехмерные структуры также проще анализировать в решетке разомкнутых кольцевых резонаторов (модель 3). Ограничимся здесь демонстрацией, по-видимому, наиболее яркого примера топологических – узловых – дискретных диссипативных солитонов (рис. 18). Условия отвечают бистабильности, когда каждый из слабо взаимодействующих элементов может находиться в основном или возбужденном состоянии. Темными шариками обозначены резонаторы в возбужденном состоянии. Такие шарики соединяются линией, если расстояние между ними не превышает шага решетки. Линия будет узловой, если она замкнута, не содержит самопересечений и не превращается в окружность при плавных деформациях. Представленные на рис. 18 примеры свидетельствуют о чрезвычайном разнообразии дискретных диссипативных солитонов, которое можно наблюдать в нелинейных магнитных метаматериалах.
аб Рис. 18. Устойчивые узловые солитоны в модели 3: «трилистник» (а) и «звезда» (б). На врезках показаны
проекции узловой линии на плоскость
Заключение
Три рассмотренных модели дискретных нелинейных диссипативных систем, резонансно раскачиваемых когерентным поддерживающим излучением, обладают качественно аналогичными свойствами, хотя и при резко различающихся масштабах (единицы нанометров для модели 1 и доли миллиметра для модели 3). Совместное рассмотрение этих моделей позволяет получить более полную картину формирования разнообразных структур. Если в модели 1 уже продемонстрирована возможность управления положением солитонов с помощью наклонного падения поддерживающего излучения, то в модели 2 привлекает появление дополнительной – поляризационной – степени свободы, а в модели 3 – сравнительная легкость формирования топологических трехмерных солитонов. Представляется, что экспериментально более просто реализовать формирование солитонов в модели 3, что послужит проверкой правильности общей картины и будет способствовать разработке новых подходов к обработке информации.
Работа поддержана грантами Министерства образования и науки Российской Федерации и НИУ ИТМО №411402.
Литература
1. Розанов Н.Н. Диссипативные оптические солитоны. От микро- к нано- и атто-. – М.: Наука, 2011. – 536 с.
2. Denz C., Flach S., Kivshar Yu.S., eds. Nonlinearities in Periodic Structures and Metamaterials. – Berlin: Springer-Verlag, 2010. – 202 p.
3. Носков Р.Е. Метаматериалы. Oт левосторонних сред к нелинейной плазмонике. – Lambert Academic Publishing, 2011. – 116 с.
4. Litchinitser N.M., Gabitov I.R., Maimistov A.I., Shalaev V.M. Linear and Nonlinear Metamaterials and Transformation Optics, in book «Tutorials in Metamaterials», Eds. M.A. Noginov and V.A. Podolskiy; series in Nano-Optics and Nanophotonics, Series Eds. S. Kawata and V.M. Shalaev. – CRC Press, Taylor & Francis Group, New York, NY, 2012. – Р. 1–27.
5. Rosanov N.N., Fedorov S.V., Shatsev A.N., Vyssotina N.V. Dissipative molecular solitons // Eur. Phys. J. D. – 2010. – V. 59. – P. 3–12.
6. Noskov R.E., Belov P.A., and Kivshar Yu. S. Subwavelength modulational instability and plasmon oscillons in nanoparticle arrays // Phys. Rev. Lett. – 2012. – V. 108. – P. 093901.
7. Noskov R.E., Belov P.A., and Kivshar Yu.S. Subwavelength plasmonic kinks in arrays of metallic nanoparticles // Optics Express. – 2011. – V. 20. – P. 2733–2739.
8. Shadrivov I.V., Zharov A.A., Zharova N.A., Kivshar Yu.S. Nonlinear magnetoinductive waves and domain walls in composite metamaterials // Photonics and nanostructures – Fundamentals and Applications. – 2006. – V. 4. – P. 69–74.
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
11
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
9. Lazarides N., Eleftheriou M., Tsironis G.P. Discrete breathers in nonlinear magnetic metamaterials // Phys. Rev. Lett. – 2006. – V. 97. – P. 157406.
10. Eleftheriou M., Lazarides N., Tsironis G.P. Magnetoinductive breathers in metamaterials // Phys. Rev. E. – 2008. – V. 77. – P. 036608.
11. Lazarides N., Tsironis G.P., Kivshar Yu.S. Surface breathers in discrete magnetic metamaterials // Phys. Rev. E. – 2008. – V. 77. – P. 065601.
12. Molina M., Lazarides N., Tsironis G.P. Bulk and surface magnetoinductive breathers in binary metamaterials // Phys. Rev. E. – 2009. – V. 80. – P. 046605.
13. Cui W., Zhu Y., Li H., Liu S. Soliton excitation in a one-dimensional nonlinear diatomic chain of split-ring resonators // Phys. Rev. E. – 2010. – V. 81. – P. 016604.
14. Розанов Н.Н., Высотина Н. В., Шацев А.Н., Шадривов И.В., Кившарь Ю.С. Гистерезис волн переключения и диссипативные солитоны в нелинейных магнитных метаматериалах // Письма в ЖЭТФ. – 2011. – Т. 93. – С. 838–841.
15. Rosanov N.N., Vysotina N.V., Shatsev A.N., Shadrivov I.V., Powell D.A., Kivshar Yu.S. Discrete dissipative localized modes in nonlinear magnetic metamaterials // Optics Express. – 2011. – V. 19. – P. 26500– 26505.
16. Rosanov N.N., Vysotina N.V., Shatsev A.N., Desyatnikov A.S., Kivshar Yu.S. Knotted solitons in nonlinear magnetic metamaterials // Phys. Rev. Lett. – 2012. – V. 108. – P. 133902.
Розанов Николай Николаевич
Высотина Нина Вениаминовна Шацев Анатолий Натанович Десятников Антон Сергеевич Шадривов Илья Владимирович Носков Роман Евгеньевич
Кившарь Юрий Семенович
– НПК «ГОИ им. С.И. Вавилова», член-корр. РАН, доктор физ.-мат. наук, начальник отдела; НИУ ИТМО, зав. кафедрой, nrosanov@yahoo.com
– НПК «ГОИ им. С.И. Вавилова», научный сотрудник, nvvyss@yahoo.com
– НПК «ГОИ им. С.И. Вавилова», кандидат физ.-мат. наук, старший научный сотрудник, anshat@yahoo.com
– Австралийский национальный университет, кандидат физ.-мат. наук, старший научный сотрудник, asd124@physics.anu.edu.au
– Австралийский национальный университет, кандидат физ.-мат. наук, старший научный сотрудник, ivs124@physics.anu.edu.au
– Санкт-Петербургский Национальный исследовательский университет информационных технологий, механики и оптики, кандидат физ.-мат. наук, старший научный сотрудник, nanometa@gmail.com
– Австралийский национальный университет, академик Австралийской академии наук, профессор, ysk@physics.anu.edu.au
12 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
ОБЗОРНАЯ СТАТЬЯ
УДК 538.9
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ В КОГЕРЕНТНО ВОЗБУЖДАЕМЫХ НАНОСТРУКТУРАХ
И МЕТАМАТЕРИАЛАХ
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов, Р.Е. Носков, Ю.С. Кившарь
Розанов Николай Николаевич, членкорреспондент Российской академии наук, доктор физико-математических наук, профессор, начальник отдела ФГУП «НПК «ГОИ им. С.И. Вавилова», заведующий кафедрой Оптики лазеров СПбНИУ ИТМО. Лауреат премии им. Д.С. Рождественского РАН (2007). Область научных интересов – физическая и нелинейная оптика, лазерная физика, физика солитонов.
Высотина Нина Вениаминовна, научный сотрудник ФГУП «НПК «ГОИ им. С.И. Вавилова». Об-
ласть научных интересов – вычислительная математика, численное моделирование.
Шацев Анатолий Натанович, кандидат физико-математических наук, старший научный сотрудник ФГУП «НПК «ГОИ им. С.И. Вавилова». Область научных интересов – нелинейная оптика, лазерная физика, физика солитонов.
Десятников Антон Сергеевич, кандидат физико-математических
наук, старший научный сотрудник Австралийского национального университета.
Шадривов Илья Владимирович, кандидат физико-математических наук, старший научный сотрудник Австралийского национального университета, стипендиат Королевы Елизаветы II.
Кившарь Юрий Семенович, академик Австралийской академии наук, профессор, руководитель отдела, ведущий ученый, научный руководитель лаборатории Австралийского национального университета, обладатель мегагранта Министерства образования и науки Российской Федерации.
Носков Роман Евгеньевич, кандидат физико-математических наук, старший научный сотрудник Санкт-Петербургского национального исследовательского университета информационных технологий, механики и оптики. Область научных интересов – оптика наноструктурированных композитных сред на основе металлических наночастиц, левые среды, электромагнитная невидимость.
Проведен сравнительный анализ локализованных структур – волн переключения и диссипативных солитонов – в трех дискретных нелинейных системах, резонансно возбуждаемых когерентным поддерживающим излучением: цепочке молекул, цепочке металлических наношаров и решетке разомкнутых кольцевых резонаторов. Ключевые слова: диссипативные солитоны, дискретные солитоны, волны переключения.
Введение
Диссипативные солитоны, представляющие локализованные за счет баланса притока и оттока энергии структуры волн, к настоящему времени более полно изучены в таких широкоапертурных оптических системах, как нелинейные интерферометры с когерентным поддерживающим излучением и лазеры с насыщающимся поглощением [1]. В этих континуальных системах продемонстрированы одномер-
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
1
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
ные и двумерные пространственные солитоны и предложены их приложения для обработки информации. Однако размеры этих солитонов составляют не менее 10 мкм, что не делает такие системы конкурентноспособными по сравнению с традиционной электроникой. Кроме того, в подобных системах сложно формировать трехмерные, в особенности топологические солитоны.
Указанные ограничения могут быть преодолены при использовании периодических наноструктур и метаматериалов, которые представляют собой набор слабо связанных нелинейных элементов, возбуждаемых внешним когерентным излучением [2–4] (рис. 1). В этой статье мы рассмотрим три родственные системы, в которых могут быть сформированы дискретные диссипативные структуры.
Модель 1 – это цепочка молекул с трехуровневой схемой электронных переходов. Внешнее излучение резонансно возбуждает молекулы, под его действием молекулы излучают, подобно диполям, и это излучение частично перехватывается соседними молекулами (диполь-дипольное взаимодействие молекул). К настоящему времени в расчетах продемонстрированы солитоны с размерами в несколько нанометров и возможность управлять их положением в цепочке и движением с помощью наклонного падения внешнего излучения [1, 5]. Перспективны для реализации таких систем J-агрегаты цианиновых красителей, обладающие коллективным (экситонным) механизмом их резонансного возбуждения лазерным излучением. В то же время еще предстоит разработка более полной теории этих принципиально квантовых объектов с выходом за рамки одночастичного приближения и учетом межмолекулярных корреляций.
у
z
E0
k
x
Рис. 1. Одномерная дискретная система. Цепочка нелинейных элементов, изображенных кружками, расположена вдоль оси x. E0 – вектор напряженности электрического поля поддерживающего излучения;
волновой вектор k падающей волны лежит в плоскости xz; φ – угол между волновым вектором k и осью z;
θ – угол поворота вектора поляризации поля относительно оси y
m–1,n
n–1 n n+1
m,n–1 m,n m,n+1
m+1,n
Рис. 2. Цепочки (слева) и двумерные решетки (справа) разомкнутых кольцевых резонаторов с различной ориентацией
Модель 2 близка к предыдущей, но в ней молекулы заменяются металлическими наношарами, погруженными в диэлектрик [6, 7]. В отличие от молекулярных J-агрегатов с фиксированным направлением дипольного момента, наношары обладают изотропным откликом, так что здесь добавляется поляризационная степень свободы, и солитоны становятся в общем случае векторными. К настоящему времени в литературе для моделей 1 и 2 изучены только одномерные системы, и рассмотрение хотя бы двумерных
2 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
схем представляет важную задачу дальнейших исследований. В то же время в модели 3 – цепочке или решетке возбуждаемых когерентным излучением разомкнутых кольцевых резонаторов (рис. 2) с нелинейными элементами, помещенными в разъемы [8–16] – возможна реализация дискретных диссипативных солитонов любой геометрической размерности. В моделях 1 и 2 аналогичен вид связи элементов (диполь-дипольное взаимодействие), но различается вид нелинейности. Последний совпадает (кубическая нелинейность) для моделей 2 и 3. Ниже мы проанализируем общие и специфические свойства солитонов в трех указанных моделях, уделяя большее внимание менее изученным в литературе вариантам – модели 2 для одномерной геометрии и модели 3 для двумерной и трехмерной геометрии.
Управляющие уравнения
Модель 1. J-агрегат моделируется линейной цепочкой N трехуровневых молекул, взаимодейст-
вующих между собой через электромагнитное поле. Два нижних уровня формируют оптический переход,
частота которого близка к лазерной частоте. Третий молекулярный уровень служит для описания анни-
гиляции двух возбуждений на соседних молекулах: одно из возбуждений дезактивируется, а другое пе-
реходит на третий уровень с последующей релаксацией либо обратно на второй, либо на первый уровень.
В результате соответственно либо одно, либо оба возбуждения исчезают из системы. Частота перехода из
основного состояния (1) в третье (3) предполагается приблизительно равной удвоенной частоте перехода
1 2 . Расстояние между молекулами d много меньше длины волны излучения. Частота лазерного излу-
чения ω0, поляризация излучения линейная. В подходе одночастичной матрицы плотности управляющие уравнения для элементов матрицы
плотности ρik, i, k = 1, 2, …, N, в приближении медленно меняющихся амплитуд имеют вид (см. [1, 5] и приведенные там ссылки)
1(1k )
1 4
N l 1
(lk
ilk )Rl Rk*
i
4
Rk*
c.c.
s
(k 22
)
[2(2k
1)
(k 22
1)
]
313(3k
)
(k
21 22
)
(k 22
)
,
2(2k )
1 4
N l 1
(lk
ilk )Rl Rk*
i
4
Rk*
c.c.
2s
(k 22
)
[2(2k
1)
(k 22
1)
]
32
(k 33
)
(k
21 22
)
(k 22
)
,
3(3k)
s
(k 22
)
[2(2k
1)
(k 22
1)
]
31 32
(k ) 33
,
Rk
( i0 )Rk
N
(lk
ilk )Rl [2(2k )
(k 11
)
]
i [2(2k )
(k 11
)
]
s
Rk
[2(2k
1)
(k 22
1)
]
.
l 1
Точка в левых частях уравнений означает временную производную; с.с. – знак комплексного со-
пряжения; 0 – значение частотной расстройки для молекул в цепочке; Rk 22(1k) exp(i0t) . Матрицы lk
и lk ( l k ) описывают запаздывающее диполь-дипольное взаимодействие молекул в агрегате:
lk
2 d 3
k0
d
cos
k0d l k lk 2
sin
k0d l k lk 3
1 3cos2
k0d 2 sin
k0d l k lk
sin2 , l k,
lk
2 d 3
cos
k0d l k l k 3
sin k0d
k0d l k lk 2
1 3cos2
k0d 2 cos
k0d l k lk
sin2 , l k,
(1)
ll 0, ll 0.
Волновое число k0 0 c , где c – скорость света в вакууме. Величина ( / )Ei – частота Раби, пропорциональная комплексной амплитуде поля, так что вещественная напряженность электрическо-
го поля в падающей волне
E Re[Ei (t) exp(i0t)] . Другими параметрами задачи являются матричный элемент дипольного момента перехода μ, угол
между дипольным моментом и осью цепочки θ, релаксационные постоянные γ и Г, параметр экситон-
экситонной аннигиляции s . Принятый подход оправдан, если уровень возбуждения молекул слабый:
1
ρ(k 11
)
1
(диагональные
элементы
матрицы
плотности
пропорциональны
населенностям
соответст-
вующих уровней).
В линейных агрегатах взаимодействие с излучением молекул у краев цепочки существенно отли-
чается от случая молекул в середине цепочки, главным образом, из-за эффекта двухэкситонной анниги-
ляции (параметр s ). В типичных условиях эффекты краев захватывают несколько десятков молекул, поэтому для средней части достаточно длинной цепочки от краевых эффектов можно отвлечься. Чтобы
исключить их в расчетах, можно представить, что цепочка замкнута в кольцо (периодические граничные
условия). Хотя молекулярные агрегаты с кольцевой геометрией и известны, молекулы в них неодинаково
взаимодействуют с излучением с линейной поляризацией (частота Раби пропорциональна косинусу угла
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
3
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
между направлениями дипольного момента молекулы и поляризации излучения). Поэтому для обеспечения одинаковых условий для различных молекул пришлось бы перейти к их взаимодействию с излучением, обладающим, например, круговой поляризацией.
В модели 2 управляющие уравнения имеют вид [6, 7]
i
dPn d
i Pn 2
Pn
G n,m
Pm
En ,
mn
i
dPn d
i Pn 2
Pn Gn,m Pm En ,
mn
(2)
где Pn 2 Pn 2 Pn 2 , верхние индексы указывают составляющие векторов отклика шаров (Р) и под-
держивающего излучения (Е) по отношению к оси цепочки (рис. 1), параметры γ, Ω, k0 и d имеют тот же смысл, что и в модели 1, а коэффициенты связи элементов аналогичны (1):
G n,m
2
k0d
2
ik0d nm
n
1 m
2
exp
ik0d n m nm
,
G n,m
ik0d nm
n
1 m
2
exp
ik0d n m nm
.
Параметр η определяется отношением диаметра шаров и шага решетки, а также диэлектрической
проницаемостью шаров и матрицы-диэлектрика [6, 7].
Управляющие уравнения для модели 3 формулируются для медленно меняющихся амплитуд электрического тока или намагниченности отдельных резонаторов [8]:
i
dn dt
(2 i | n
|2 )n
Cn
.
Здесь время t нормировано на период собственных колебаний изолированного резонатора, Ω – от-
стройка собственной частоты колебаний от частоты когерентного поддерживающего излучения, норми-
рованная на собственную частоту, γ – затухание резонатора с той же нормировкой, 1 – знак коэф-
фициента керровской нелинейности. Коэффициент κ определяет силу связи ближайших резонаторов, и
Cn 2[(nx , ny , nz 1) (nx , ny , nz 1)] [(nx 1, ny , nz ) (nx 1, ny , nz )]
[(nx , ny 1, nz ) (nx , ny 1, nz )], n (nx , ny , nz ).
Целочисленные индексы nx , ny , nz определяют положение резонатора в решетке. Управляющие
уравнения для каждой модели дополняются начальными условиями, роль которых особо важна в услови-
ях бистабильности или тристабильности (см. ниже).
Одномерные локализованные структуры (модель 2)
Квазиоднородные распределения. Для идеальной цепочки бесконечной длины с одинаковым стационарным состоянием всех элементов (пространственно-однородные стационарные состояния) из (2)
при Pn P0 и нулевых производных получаем алгебраическое уравнение 5-й степени:
P0
2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
E0
2
cos
2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
sin 2
P0
2
Re
i 1
Ai
2
Im
i 1
Ai
2
,
где
Aj j3 ik0d j2 k0d 2 j1 exp ik0d j ,
Aj 2 j3 ik0d j2 exp ik0d j .
При 0, / 2 степень полинома снижается до третьей, тогда одному и тому же значению интен-
сивности поддерживающего излучения могут отвечать три значения отклика. С учетом заведомой неустойчивости промежуточного значения это означает бистабильность (рис. 3, а). Этот вариант характерен
4 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
также для скалярных моделей 1 и 3, тогда как тристабильность становится возможной для векторной модели 2 (рис. 3, б).
0,2 0,3
0,16 0,12 0,2
|P0,|| |2 |P0 |2
0,08 0,1 0,04
00
0 0,0001 0,0002 0,0003 0,0004
0 0,0002 0,0004 0,0006 0,0008
|Е0,|| |2 а
|Е0 |2 б
Рис. 3. Зависимость поперечной и продольной поляризации P0 2 , P0 2 от интенсивности внешнего поля
I E0 2 , E0 2 для бесконечной цепочки; перпендикулярная поляризация падающей волны (θ = 0) –
сплошная кривая, параллельная поляризация (θ = /2) – пунктирная кривая, 0,1 (а).
Зависимость P0 2 от E0 2 при 0, 2 , 0, 2 для бесконечной цепочки (б). Устойчивость для
конечной цепочки проверялась численно. Верхняя ветвь устойчива целиком, нижняя ветвь устойчива при
интенсивностях падающей волны от 0 до ~0,0004, остальные ветви неустойчивы (расчет при N 100 )
|Pn|2
0,16 0,01 0,008
0,12 0,006 0,08 0,004
|Pn|||2
0,04 0,002
|Pn|2
00
0
0,1 0,08
20 40 60 80 100 n
а
0
0,1 0,08
20 40 60 80 100 n
б
0,06 0,06
|Pn|||2
0,04 0,04
0,02 0,02
00
0 20 40 60 80 100 n
0 20 40 60 80 100 n
вг
Рис. 4. Интенсивность перпендикулярной (а, в) и параллельной (б, г) составляющих отклика
при E0 2 0, 000088 и 0,12 . Для (в, г) начальные условия нулевые, для (а, б) ненулевые
0,2 0,2
0,16 0,16 0,12 0,12
|P|2
|P|2
0,08 0,04 b a
0,08 0,04 b a
00
0 40 80 120 160 200
0 40 80 120 160 200
nn
аб Рис. 5. Профили неподвижной (a) и движущейся (b) волны переключения. I = 0,55×10–4 , скорость V < 0
(область 2 рис. 7) (а); I = 0,88×10–4,V > 0 (область 5 рис. 7) (б).
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
5
|P|2 w
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
Для цепочки конечной длины вследствие краевых эффектов однородные распределения заменяются на неоднородные. При этом следует говорить уже о пространственной бистабильности (бистабильности пространственных структур, устанавливающихся при различных начальных условиях) [1]. Еще одно важное условие реализации режима состоит в требовании его устойчивости относительно малых возмущений. Для бесконечной цепочки линейный анализ устойчивости стационарных однородных распределений выполняется стандартным образом [6, 7], в других случаях обычно необходимы численные расчеты. В результате развития неустойчивости в цепочке формируются и устанавливаются квазипериодические структуры. Интересной особенностью, известной и для континуальных систем [1], является возможность установления при одних и тех же параметрах системы, но при различных начальных условиях, либо квазиоднородной (рис. 4, а, б), либо квазипериодической (рис. 4, в, г) структуры.
Волны переключения в скалярной системе. Как и в континуальных системах, внутри интервала бистабильности квазиоднородных распределений существуют волны переключения (ВП) (другие названия – доменные стенки и кинки). При удалении от фронта волны переключения в одну или другую сторону распределение асимптотически приближается к одному из возможных квазиоднородных распределений, различающихся по разные стороны от фронта. Фронт волны может быть неподвижным (неподвижная волна переключения, рис. 5, кривая а, или движущимся (движущаяся волна переключения, рис. 5, кривая б). Если для континуальных систем скорость волны переключения обращается в нуль при одном («максвелловском») значении интенсивности поддерживающего излучения [1], то в дискретных системах это происходит на целом интервале интенсивности. Общей особенностью дискретных волн переключения служит то, что форма фронта и скорость движущейся волны меняются со временем квазипериодически (рис. 6). В связи с этим имеет смысл говорить о средней скорости таких волн.
0,3 12
0,2 8
0,1 4
00 500 600 700 800 900 1000
t Рис. 6. Временное изменение максимальной интенсивности отклика на фронте (штриховая кривая, левая шкала) и ширины фронта w (сплошная кривая, правая шкала, в единицах расстояния между элементами)
движущейся волны переключения в условиях рис. 5, б (кривая b, I = 0,88 10–4). Ширина фронта w
определялась как длина промежутка, на котором интенсивность волны спадает от максимального значения не менее чем в 100 раз
При изменении интенсивности поддерживающего излучения скорость волны переключения меняется, в том числе по знаку. При этом мы считаем скорость V положительной, если со временем при передвижении фронта увеличивается число элементов с большей амплитудой отклика (рис. 5, б, в отличие от рис. 5, а, где V < 0). Зависимость скорости волны переключения V от интенсивности падающего излучения I приведена на рис. 7. Для удобства диапазон изменения I разделен на 6 областей. Хотя в области 1 (0,0381×10–4 < I < 0,46×10–4) и имеется бистабильность однородных состояний бесконечной цепочки, повидимому, область притяжения верхнего квазиоднородного состояния столь узка, что практически оно неустойчиво (из-за неизбежных в расчете возмущений). При этом исходное ступенчатое распределение быстро сваливается в нижнее квазиоднородное состояние, так что волны переключения отсутствуют. В области 2 (0,46×10–4 < I < 0,56×10–4) существует как неподвижная, так и отрицательная (V < 0) волны переключения (рис. 5, а). В области 3 (0,56×10–4 < I < 0,836×10–4) существует только неподвижная волна переключения. В области 4 (0,836×10–4 < I < 0,9216×10–4) существует как неподвижная, так и движущаяся с положительной скоростью (V > 0) волна переключения (рис. 5, б). При I > 0,9216×10–4 происходит скачкообразный переход из режима с V = 0 в режим с V > 0. В области 5 (0,9216×10–4 < I < 2,045×10–4) существуют только волны переключения с положительной скоростью. Заметим, что приблизительно при I > 1,4×10–4 уже трудно получить одиночную волну переключения, так как из-за краевых эффектов неизбежно порождается встречная волна.
Все эти области 1–5 находятся внутри области пространственной бистабильности квазиоднородных состояний. Однако и вне этой области, но внутри области бистабильности для бесконечной цепочки, имеются метастабильные волны переключения. Дело в том, что дестабилизирующее действие краевых
6 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
эффектов происходит не мгновенно, и если время их развития превышает время распространения фронта по конечной цепочке, то неустойчивость будет существенно подавлена. Поэтому при 2,045×10–4 < I < < 2,9×10–4 (область 6) имеются метастабильные волны переключения с положительной скоростью. Наконец, при I > 2,9×10–4 даже метастабильных (со значительным временем жизни) волн переключения нет, переход на исходной ступеньке из нижнего состояния в верхнее происходит практически сразу, скачком.
0,08 1 2 3 4
5
6
0,08 1 2
3 45
0,06 0,06
0,04 0,04
0,02 0,02
V/c |Pn|2
V/c |Pn|||2
00
–0,02 –0,04 –0,06
12 E|E00|21100–44
3 –0,02 –0,04 –0,06
0,6 0,8 |EE00|22101–044
1
–0,08
–0,08
–0,1 –0,1
–0,12
–0,12
–0,14
–0,14
аб
Рис. 7. Зависимость скорости волны переключения V c от интенсивности поддерживающего излучения
I = |E0|2 (а). То же, укрупненно (б)
Векторные волны переключения ( 0 ). Если внешнее излучение содержит обе поляризационные компоненты, то, как указывалось выше, возникают следующие особенности. Во-первых, в соответствии с рис. 3, б, может иметь место мультистабильность (тристабильность) квазиоднородных состояний, из-за чего возможны новые типы волны переключения и, как результат их столкновения, новые типы солитонов. Во-вторых, развитие модуляционной неустойчивости некоторых из квазиоднородных состояний может при нелинейном ограничении приводить к бистабильности между квазиоднородным и пространственно промодулированным состоянием и соответственно к волнам переключения между такими состояниями. В этом разделе мы ограничимся изучением только второй возможности, выбирая параметры так, чтобы они отвечали бистабильности квазиоднородных распределений.
При 0,12 и E0 2 0, 0001 квазиоднородное состояние на верхней ветви устойчиво, а на ниж-
ней – модуляционно неустойчиво, со временем оно превращается в стационарное пространственно промодулированное состояние. Как показывает рис. 8, исходное ступенчатое распределение возбуждает волну переключения, фронт которой со временем надвигается на промодулированное состояние. Для цепочки конечной длины фронт останавливается на некотором расстоянии от края цепочки.
0,2 0,08
0,16 3
0,12 0,08
12
0,04
0,06 0,04 0,02 1 2
3
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 8. Интенсивность перпендикулярной (a) и параллельной (б) составляющей поля при 0,12 и
E0 2 0, 0001 . Исходное распределение – ступенчатое. Кривые 1, 2, 3 – соответственно распределения полей в моменты времени t 500 , 2000 и 3500 (установившаяся структура)
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
7
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
При меньшем значении интенсивности E0 2 0, 000088 и 0,12 фронт между квазиоднород-
ным и промодулированным состояниями неподвижен даже вдали от краев цепочки. Устанавливается стационарное распределение (рис. 9), т.е. ступенька практически сразу останавливается.
Наконец, при E0 2 0, 0001 и 0,1 и нижнее, и верхнее квазиоднородные состояния устойчи-
вы. В этом случае волны переключения – того же типа, как в скалярном случае ( 0 ). Движение фронта ВП не останавливается краевыми эффектами (рис. 10).
0,16 0,04
0,12 0,03
|Pn|||2
|Pn|2
0,08 0,02
0,04 0,01
0 0 20 40 60 80 100 n
0 0 20 40 60 80 n
аб
Рис. 9. Интенсивность перпендикулярной (a) и параллельной (б) составляющей отклика
при E0 2 0, 000088 и 0,12
100
0,2 0,03
|Pn|2
0,16
0,12
0,08
12 3
4
0,04
|Pn|||2
0,02
0,01
12 3
4
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 10. Движущаяся волна переключения в условиях классической бистабильности
для перпендикулярной (a) и параллельной (б) составляющей отклика. Кривые 1, 2, 3, 4 – соответственно
распределения полей в моменты времени t 500, 2500, 4500, 7000 ; E0 2 0, 0001 , 0,1
Скалярные дискретные диссипативные солитоны (модель 2, θ = 0)
Важным сценарием формирования диссипативных солитонов служит столкновение волн пере-
ключения. Поскольку в рассматриваемой системе возможны такие волны различных типов, то число
сценариев увеличивается.
Столкновение волн переключения с нулевой и положительной скоростью ( I 0,88 104 , область 4
рис. 7). На рис. 11, а, справа движется волна навстречу неподвижной волне в левой половине цепочки.
0,25 0,2
0,2 0,16 5
|Pn|2 |Pn|2
0,15 0,1 4 3 2 1
0,05
0,12 0,08 0,04
1234 4321
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 11. Динамика столкновения волн переключения. Кривые 1, 2, 3, 4 – соответственно распределения
интенсивности отклика в моменты времени t 1000 , 2000, 3000 и 4000, I 0,88 104 (а),
кривые 1, 2, 3, 4, 5 – соответственно распределения интенсивности отклика в моменты времени
t 1000 , 1500, 2000, 2500 и 5000, I 0,95 104 (б)
Движущаяся справа волна останавливается с образованием темного солитона, включающего 15
элементов. Левая «неподвижная» волна в самом конце столкновения сдвигается и останавливается. Ана-
8 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
|Pn|2 |Pn|2
логичные, но более узкие темные солитоны формируются при бóльших интенсивностях I, вплоть до
I 0,94 104 , когда устанавливается солитон с шириной, отвечающей 6 элементам. Но уже при
I 0,95 104 сталкивающиеся волны переключения не останавливаются, а аннигилируют, результи-
рующее поле стационарно и имеет квазиоднородное распределение (кривая 5 на рис. 11, б). Столкновение волн переключения с нулевой и отрицательной скоростью (область 2 рис. 7). Ана-
логично предыдущему сценарию, в результате взаимодействия двух волн переключения образуются солитоны, но теперь они светлые. На рис. 12, а, приведены профили результирующих стационарных распределений поля для различных значений интенсивности поддерживающего излучения I. Самый узкий
солитон, соответствующий I 0,5510 104 , имеет ширину 2 элемента. Динамика формирования таких
солитонов показана на рис. 11, б. Аннигиляцию волн переключения с нулевой и отрицательной скоростью иллюстрирует рис. 13.
0,16 0,16
0,12 0,12
0,08 0,08
0,04
0,04
43
21
0 0 20 40 60 80 100 n
0 0 20 40 60 80 100 n
аб
Рис. 12. Установившиеся солитоны при различных значениях интенсивности: пунктирная, сплошная, пунк-
тирная жирная, сплошная жирная кривые соответствуют значениям
E0 2 104 0,555; 0,552; 0,5513; 0,5510 (а); динамика образования солитона с шириной w 2 при столк-
новении волны переключения, движущейся налево (кривые 1–4 отвечают моментам времени t = 0, 500,
1500, 3000) с неподвижной волной переключения; I 0, 0000551 (б).
0,16
0,12
|Pn|2
0,08 3 2 1 0,04 4
0 0 40 80 120 160 200 n
Рис. 13. Левый фронт практически неподвижен, а правый фронт движется налево (кривые 1–4 соответствуют моментам t = 200, 1600, 2800, 4000), где сталкивается с неподвижным фронтом, и фронты
схлопываются; I = |E0|2 =0,5476×10–4
0,2
0,16
0,12
|Pn|2
0,08
0,04
0
0 20 40 60 80 100
n Рис. 14. Одностороннее расширение области возбуждения цепочки между фронтами неподвижной ВП (слева) и ВП, движущейся с положительной скоростью. Моменты (слева направо) t = 500, 1000, 1500, 2000;
I = |E0|2 =0,8491×10–4 (область 3 рис. 7)
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
9
|Pn|2 |Pn|2 |Pn|2 |Pn|2
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
Если фронт движущейся волны переключения не приближается к фронту неподвижной волны, а удаляется от него, и начальное расстояние между фронтами не слишком мало, то возникает расширяющаяся со временем область повышенного возбуждения (рис. 14). При значительном расстоянии между фронтами эти две волны переключения практически не взаимодействуют друг с другом (расширение в одну сторону).
Диссипативные солитоны могут формироваться и по другим сценариям, например, при подаче на цепочку дополнительного «записывающего» узкого пучка-импульса или нескольких пучков-импульсов. Примеры установившихся структур приведены на рис. 15. При выходе за границы области устойчивости, например, для модели 1, характерно появление пульсирующих солитонов. Рис. 16 демонстрирует, также для модели 1, возможность управляемого передвижения устойчивых солитонов по цепочке.
0,12
0,16 0,12
0,12
0,08
0,12 0,08
0,08
0,08
0,04
0,04 0,04
0,04
0
00
0
0 20 40 60 80 100 n
0 20 40 60 80 100 0 20 40 60 80 100 nn
0 20 40n 60 80 100
Рис. 15. Примеры узких светлых солитонов; I = |E0|2 =0,5329×10–4 (область 1 рис. 7)
11 11 11 11 0,98 0,98 0,98 0,98
0,94 0,94 0,94 0,94
0,9 0
0,9
100 200 300 0 n а
0,9 100 200 300 0
n
б
100 200 n
0,9 300 0
в
100 200 300 n г
Рис. 16 (модель 1). Управление движением солитона при наклонном падении поддерживающего излучения. Профили населенности основного состояния молекул в цепочке в различные моменты времени t: t = 0, угол падения θ = 0, солитон неподвижен (а); t = 50, угол падения θ = 0,1, солитон движется направо (б); t = 300, угол падения θ = 0,1, солитон останавливается у правого края цепочки (в); t = 700,
угол падения θ = – 0,1, солитон двигается налево и останавливается у левого края цепочки (г)
Двумерные локализованные структуры (модель 3)
Более удобна для анализа двумерных решеток схема разомкнутых кольцевых резонаторов (модель 3, рис. 2) В этом случае схема существенно анизотропна, так как взаимодействие резонаторов зависит от их относительной ориентации. При фиксированной ориентации устойчивыми оказываются только сравнительно узкие двумерные солитоны (рис. 17). Широкие же исходные структуры обычно распадаются на набор узких солитонов вследствие модуляционной неустойчивости.
0,54 0,55
0,34 0,35
0,14 0,15
Рис. 17. Два примера двумерных солитонов в модели 3; шкала для намагниченности | ψ |2
указана справа от рисунков, положение резонаторов задается целыми индексами nx и nz
10 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)
Н.Н. Розанов, Н.В. Высотина, А.Н. Шацев, А.С. Десятников, И.В. Шадривов…
Трехмерные локализованные структуры (модель 3)
Трехмерные структуры также проще анализировать в решетке разомкнутых кольцевых резонаторов (модель 3). Ограничимся здесь демонстрацией, по-видимому, наиболее яркого примера топологических – узловых – дискретных диссипативных солитонов (рис. 18). Условия отвечают бистабильности, когда каждый из слабо взаимодействующих элементов может находиться в основном или возбужденном состоянии. Темными шариками обозначены резонаторы в возбужденном состоянии. Такие шарики соединяются линией, если расстояние между ними не превышает шага решетки. Линия будет узловой, если она замкнута, не содержит самопересечений и не превращается в окружность при плавных деформациях. Представленные на рис. 18 примеры свидетельствуют о чрезвычайном разнообразии дискретных диссипативных солитонов, которое можно наблюдать в нелинейных магнитных метаматериалах.
аб Рис. 18. Устойчивые узловые солитоны в модели 3: «трилистник» (а) и «звезда» (б). На врезках показаны
проекции узловой линии на плоскость
Заключение
Три рассмотренных модели дискретных нелинейных диссипативных систем, резонансно раскачиваемых когерентным поддерживающим излучением, обладают качественно аналогичными свойствами, хотя и при резко различающихся масштабах (единицы нанометров для модели 1 и доли миллиметра для модели 3). Совместное рассмотрение этих моделей позволяет получить более полную картину формирования разнообразных структур. Если в модели 1 уже продемонстрирована возможность управления положением солитонов с помощью наклонного падения поддерживающего излучения, то в модели 2 привлекает появление дополнительной – поляризационной – степени свободы, а в модели 3 – сравнительная легкость формирования топологических трехмерных солитонов. Представляется, что экспериментально более просто реализовать формирование солитонов в модели 3, что послужит проверкой правильности общей картины и будет способствовать разработке новых подходов к обработке информации.
Работа поддержана грантами Министерства образования и науки Российской Федерации и НИУ ИТМО №411402.
Литература
1. Розанов Н.Н. Диссипативные оптические солитоны. От микро- к нано- и атто-. – М.: Наука, 2011. – 536 с.
2. Denz C., Flach S., Kivshar Yu.S., eds. Nonlinearities in Periodic Structures and Metamaterials. – Berlin: Springer-Verlag, 2010. – 202 p.
3. Носков Р.Е. Метаматериалы. Oт левосторонних сред к нелинейной плазмонике. – Lambert Academic Publishing, 2011. – 116 с.
4. Litchinitser N.M., Gabitov I.R., Maimistov A.I., Shalaev V.M. Linear and Nonlinear Metamaterials and Transformation Optics, in book «Tutorials in Metamaterials», Eds. M.A. Noginov and V.A. Podolskiy; series in Nano-Optics and Nanophotonics, Series Eds. S. Kawata and V.M. Shalaev. – CRC Press, Taylor & Francis Group, New York, NY, 2012. – Р. 1–27.
5. Rosanov N.N., Fedorov S.V., Shatsev A.N., Vyssotina N.V. Dissipative molecular solitons // Eur. Phys. J. D. – 2010. – V. 59. – P. 3–12.
6. Noskov R.E., Belov P.A., and Kivshar Yu. S. Subwavelength modulational instability and plasmon oscillons in nanoparticle arrays // Phys. Rev. Lett. – 2012. – V. 108. – P. 093901.
7. Noskov R.E., Belov P.A., and Kivshar Yu.S. Subwavelength plasmonic kinks in arrays of metallic nanoparticles // Optics Express. – 2011. – V. 20. – P. 2733–2739.
8. Shadrivov I.V., Zharov A.A., Zharova N.A., Kivshar Yu.S. Nonlinear magnetoinductive waves and domain walls in composite metamaterials // Photonics and nanostructures – Fundamentals and Applications. – 2006. – V. 4. – P. 69–74.
Научно-технический вестник информационных технологий, механики и оптики, 2012, № 4 (80)
11
ДИСКРЕТНЫЕ ВОЛНЫ ПЕРЕКЛЮЧЕНИЯ И ДИССИПАТИВНЫЕ СОЛИТОНЫ …
9. Lazarides N., Eleftheriou M., Tsironis G.P. Discrete breathers in nonlinear magnetic metamaterials // Phys. Rev. Lett. – 2006. – V. 97. – P. 157406.
10. Eleftheriou M., Lazarides N., Tsironis G.P. Magnetoinductive breathers in metamaterials // Phys. Rev. E. – 2008. – V. 77. – P. 036608.
11. Lazarides N., Tsironis G.P., Kivshar Yu.S. Surface breathers in discrete magnetic metamaterials // Phys. Rev. E. – 2008. – V. 77. – P. 065601.
12. Molina M., Lazarides N., Tsironis G.P. Bulk and surface magnetoinductive breathers in binary metamaterials // Phys. Rev. E. – 2009. – V. 80. – P. 046605.
13. Cui W., Zhu Y., Li H., Liu S. Soliton excitation in a one-dimensional nonlinear diatomic chain of split-ring resonators // Phys. Rev. E. – 2010. – V. 81. – P. 016604.
14. Розанов Н.Н., Высотина Н. В., Шацев А.Н., Шадривов И.В., Кившарь Ю.С. Гистерезис волн переключения и диссипативные солитоны в нелинейных магнитных метаматериалах // Письма в ЖЭТФ. – 2011. – Т. 93. – С. 838–841.
15. Rosanov N.N., Vysotina N.V., Shatsev A.N., Shadrivov I.V., Powell D.A., Kivshar Yu.S. Discrete dissipative localized modes in nonlinear magnetic metamaterials // Optics Express. – 2011. – V. 19. – P. 26500– 26505.
16. Rosanov N.N., Vysotina N.V., Shatsev A.N., Desyatnikov A.S., Kivshar Yu.S. Knotted solitons in nonlinear magnetic metamaterials // Phys. Rev. Lett. – 2012. – V. 108. – P. 133902.
Розанов Николай Николаевич
Высотина Нина Вениаминовна Шацев Анатолий Натанович Десятников Антон Сергеевич Шадривов Илья Владимирович Носков Роман Евгеньевич
Кившарь Юрий Семенович
– НПК «ГОИ им. С.И. Вавилова», член-корр. РАН, доктор физ.-мат. наук, начальник отдела; НИУ ИТМО, зав. кафедрой, nrosanov@yahoo.com
– НПК «ГОИ им. С.И. Вавилова», научный сотрудник, nvvyss@yahoo.com
– НПК «ГОИ им. С.И. Вавилова», кандидат физ.-мат. наук, старший научный сотрудник, anshat@yahoo.com
– Австралийский национальный университет, кандидат физ.-мат. наук, старший научный сотрудник, asd124@physics.anu.edu.au
– Австралийский национальный университет, кандидат физ.-мат. наук, старший научный сотрудник, ivs124@physics.anu.edu.au
– Санкт-Петербургский Национальный исследовательский университет информационных технологий, механики и оптики, кандидат физ.-мат. наук, старший научный сотрудник, nanometa@gmail.com
– Австралийский национальный университет, академик Австралийской академии наук, профессор, ysk@physics.anu.edu.au
12 Научно-технический вестник информационных технологий, механики и оптики,
2012, № 4 (80)