УДК 535.14
ФАЗОВЫЕ ПЕРЕХОДЫ И ХРАНЕНИЕ КВАНТОВОЙ ОПТИЧЕСКОЙ ИНФОРМАЦИИ НА ОСНОВЕ ПОЛЯРИТОНОВ В ПРОСТРАНСТВЕННО-ПЕРИОДИЧЕСКИХ
СТРУКТУРАХ
А. П. Алоджанц1, С. М. Аракелян1, С. Н. Багаев2
владимирский государственный университет, Владимир 2 Институт лазерной физики РАН, Новосибирск
[email protected], [email protected], [email protected] PACS 42.50.Nn, 05.30.Jp, 32.70.Jz
В работе обсуждаются различные аспекты фазовых переходов и бозе-эйнштейновской конденсации (БЭК) для двумерного газа поляритонов, формирующихся в полупроводниковых и атомных микро- и наноструктурах; рассмотрение ведется с точки зрения существующих экспериментальных возможностей и имеющихся технологических достижений. Предложен новый тип пространственно-периодической структуры - поляритонный кристалл (ПК), образованный двумерной решеткой удерживаемых двухуровневых атомов и взаимодействующих с электромагнитным полем в решетке туннельно-связанных микрорезонаторов. Основной акцент сделан на выяснении механизмов формирования когерентных поляритонов, определяемых связанными состояниями оптического поля и атомной среды для нижней ветви дисперсионной кривой поляритонов. Проанализированы случаи существенного изменения (замедления) групповой скорости оптических волновых пакетов в условиях квазиконденсации (или истинной конденсации) поляритонов. Показано, что структура ПК позволяет полностью локализовать поляритоны, что может быть использовано как для получения квазиконденсации поляритонов нижней дисперсионной ветви, так и для записи и хранения квантовой оптической информации распространяющегося светового излучения.
Ключевые слова: фазовые переходы, пространственно-периодические среды и фотонные кристаллы, поля-ритоны, квантовая информация.
1. Введение. Фазовые переходы в системе двумерного газа поляритонов
Эксперименты по бозе-эйнштейновской конденсации (БЭК) макроскопически большого числа атомов (см. например, [1]) оказали большое влияние на развитие различных направлений современной квантовой и лазерной физики, а также в приложениях для задач квантовой информации и фотоники. При реализации БЭК, когда в условиях температурного фазового перехода макроскопически большое число атомов находится на основном (нижнем) квантовом уровне, формируется новое когерентное состояние вещества. Это проявляется в том, что, например, в предельном случае температуры Т = 0 ансамбль атомов конденсата, как и каждый атом в отдельности, описываются общей волновой функцией Ф(г, ¿), соответствующей когерентному состоянию. В данном аспекте явление БЭК аналогично процессу лазерной генерации, когда, например, в резонаторных схемах происходит фазовая синхронизация мод светового излучения (ср. с [2]). Отличительной особенностью состояния БЭК атомов является их равновесное состояние, что для «стандартных» лазерных систем (с инверсией населенностей) является практически не выполнимым ввиду
не равновесного характера взаимодействия среды и квантованного поля излучения в резонаторе (ср. с [3]). Для источников когерентного излучения без инверсии населенностей, таких как лазеры без инверсии населенностей [4], низкопороговые поляритонные лазеры [5] и др. вопросы получения термодинамического равновесия связанных состояний среды и поля являются на сегодняшний день чрезвычайно актуальными и требуют своего отдельного рассмотрения в каждом конкретном случае, в т. ч. и с учетом временного масштаба протекающих физических процессов (ср. с [6]). Именно такие, безинверсные системы анализируются в данной работе.
Несмотря на впечатляющие успехи в плане поучения БЭК атомов, необходимость поддерживания сверхнизких температур (десятки нК) налагает принципиальные ограничения на возможность их практического использования. В этой связи чрезвычайно важным становится вопрос о получении макроскопически когерентных состояний вещества при относительно высоких (вплоть до комнатных) температурах.
Возможным решением данной проблемы может быть использование квазиконденсации двумерного слабовзаимодействующего бозе-газа поляритонов; такая квазиконденсация происходит в результате фазового перехода Костерлица—Таулеса (Kosterlitz—Thouless) к сверхтекучему состоянию двумерных бозе-систем, для которых истинная (true) бозе-эйнштейновская конденсация (в отсутствии удержания частиц газа в ловушке) невозможна. Хотя данные состояния нельзя считать конденсатом в строгом термодинамическом смысле из-за неравновесного состояния системы в целом, однако при определенных условиях, налагаемых на характер реализуемых взаимодействий, и в заданных временных интервалах термализации, поляритоны, действительно, могут конденсироваться; их распределение удовлетворяет функции распределения Бозе—Эйнштейна идеального газа бозонов (см. [7]).
В настоящее время определенный успех в получении квазиконденсации поляри-тонов достигнут для полупроводниковых микро- и наноструктур — квантовых ям, помещенных в полупроводниковые микрорезонаторы, которым посвящено большое количество публикаций — см, например, [8-11]. В частности, в работе [8] для этих целей был задействован полупроводниковый микрорезонатор на основе структуры CdTe/CdMgTe с помещенным внутрь образцом, состоящим из 16 квантовых ям. В эксперименте регистрировалось распределение фотолюминесценции как функции от угла в, определяемой выражением k\\ = (E/hc) sine, где к\\ - составляющая волнового вектора поляритонов с энергией Е, параллельно зеркалам резонатора. При достижении пороговой мощности Pth — 1, 67 кВт/см2 и критической температуре Тс = 19 K поляритоны начинают заполнять низшее энергетическое состояние, соответствующее их нижней дисперсионной ветви. При этом необходимым условием является так называемое условие сильной связи между экситонами и квантованным полем в резонаторе, когда выполняется неравенство (ср. с [2]):
где шс — кооперативная частота, определяющая коллективные эффекты взаимодействия эк-ситонов и поля, ^о — частота перехода, й — дипольный момент перехода, п- плотность экси-тонного газа, тсон - характерное время когерентности среды. При этом само поле является слабым (по числу фотонов).
В таблице 1 приведены численные данные сравнительных характеристик получения бозе-конденсации для двух систем — атомного газа и газа поляритонов. Принципиальным моментом здесь является то, что относительно высокие температуры фазового перехода для поляритонов фактически обусловлены их малой массой.
(1)
Системы Атомный газ Резонаторные поляри-
Типичные тоны в полупроводни-
характеристики ковых наноструктурах
Эффективная масса m/mat, 1 Ю-9 ~ Ю-10
Tnat — масса свободного атома
Боровский радиус ав 0,01 нм 10 нм
Расстояние между частицами n~1/D, D — размерность газа 0,1 мкм 1 мкм
Время жизни 1с 1 ~ 10 пс
Время термализации 1 мс 1 ~ 10 пс
Критическая температура Тс 1 нК ~ 1 мкК 1 К ~ 300 К
ТАБЛИЦА 1. Сравнительные характеристики БЭК для атомов и поляритонов, образующихся в полупроводниковых структурах, ср. с [1,7-14]
Обратимся теперь к основным статистическим свойствам ансамбля частиц в двумерном пространстве, к которым можно отнести резонаторные поляритоны при соблюдении условия сильной связи (1).
Переход к квантовому газу бозе-частиц, образующихся на двумерной поверхности, в этом случае реализуется при Т ^ Та, где
_ 2тгП2п2
- -;--(Я)
ткв
представляет температуру квантового вырождения газа, п2- двумерная плотность частиц в плоскости, перпендикулярной оси резонатора.
Однако, для идеального газа, даже при наличии вырождения, как это известно, не существует (отличной от нуля) температуры фазового перехода к конденсации. Вместе с тем истинная бозе-эйнштейновская конденсация всё же возможна при наличии внешнего потенциала удерживающего частицы в ловушке. В этом случае критическая температура перехода к истинной конденсации может быть определена из выражения (ср. с [12,13]):
I-,
Т= ! (3)
с k2Bmrlg2(2,0y К )
где д2(2, 0) ~ 1, 645, г0 — поперечный размер области удержания, N — число частиц, и0 — потенциал удержания в центре ловушки.
На практике обычно используют гармоническую ловушку с потенциалом удержания
частиц:
U (г) = С/о — , (4)
Г0
где г — поперечная координата. Например, авторы работы [9] для целей удержания поляритонов в условиях (4) использовали постоянное механическое напряжение, приложенное к полупроводниковому образцу. При атомно-оптических взаимодействиях удержание поляритонов может быть достигнуто с помощью специальных волноводов (см. [14]). Заметим, что в условиях имеющихся экспериментов по получению БЭК для связанных состояний среды и поля тепловая энергия квТ существенно больше энергии Ьш^ ~ л/Щ), определяющей удержание частиц в ловушке, ср. с [6,9].
В результате при температурах Т ^ Тс двумерный бозе-газ, удерживаемый гармонической ловушкой, находится в основном состоянии; число частиц в этом состоянии М0
зависит от температуры в соответствии с соотношением:
N0 « N (1 - (Т/Тс)2) .
(5)
С учетом взаимодействия частиц картина существенно меняется. Даже слабое их взаимодействие приводит к так называемой квазиконденсации двумерного бозе-газа, которую называют фазовым переходом Костерлица—Таулеса (см., например, [13]). Суть данного явления заключается в том, что на двумерной поверхности, ограничивающей неидеальный газ бозе-частиц при температурах ниже Ткт, образуются капельки БЭК, обладающие всеми свойствами сверхтекучести, но не являющиеся фазово коррелированными. При этом сама критическая температура Т^т может быть найдена из выражения:
где пв — плотность сверхтекучей жидкости на двухмерной поверхности. Как показывают расчёты (см., например, [13]), величина Ткт приблизительно в 4 раза меньше температуры вырождения для идеального бозе-газа.
Таким образом, в двумерном случае бозе-газ может претерпевать два различных типа фазовых перехода. При этом для атомных систем значения критических температур истинной конденсации (3), а также квазиконденсации (6) являются чрезвычайно низкими и по-прежнему лежат в областях десятков (сотен) нанокельвинов. С другой стороны, для рассматриваемых в работе бозе-частиц — поляритонов, фазовые переходы могут происходить при значительно более высоких температурах. Это связано в первую очередь с тем, что их эффективная масса существенно меньше массы атомов (см. таблицу 1). Именно данное обстоятельство делает их привлекательными для физических исследований и перспективными для возможных практических приложений.
Следует однако, иметь ввиду, что для полупроводниковых структур характерное время Тсоь достаточно мало и составляет всего 10_12-10_11 с (см. таблицу 1). Кроме того, эффекты ионизации экситонов, а также наличие примесей в образце, которые ограничивают возможности микроэлектронных систем в рамках существующих технологий, играют негативную роль при получении макроскопического конденсата поляритонов особенно при высоких температурах. Это делает проблематичным его использование в приложениях, например, для целей квантовой обработки информации. В этой связи актуальным представляется получение поляритонного конденсата в атомных микро- и наноструктурах, где время жизни ограничено временем жизни возбужденных атомных состояний, составляющих десятки наносекунд. В частности, в работе [15] впервые было предсказано время термализации связанных (одетых) атомно-оптических состояний и продемонстрировано формирование ими термодинамически равновесного состояния при взаимодействии атомов рубидия с частицами буферного газа при высоком давлении и высокой температуре (530 К), находящихся в условиях оптических столкновений. Явление ограничивалось только спонтанными переходами для возбужденных состояний этих атомов.
В настоящей работе будут рассмотрены квантовые статистические свойства газа поляритонов, формирующихся в атомных микроструктурах, температура которых определяется временем удержания атомов в решетке. Материал статьи распределен следующим образом.
В разделе 2 будет предложена решеточная модель для поляритонов, основанная на удержании двухуровневых атомов в оптической решетке и взаимодействующих с световым полем в системе резонаторов, что позволяет осуществить их (поляритонов) квазиконденсацию. Предельно возможная групповая скорость поляритонов с учетом их взаимодействия, а также неравновесных процессов релаксации обсуждается в разделе 3. Физические основы
(6)
записи информации на основе поляритонов нижней дисперсионной ветви, формирующихся в рассматриваемой пространственно-периодической структуре, будут предложены в разделе 4. В заключении сформулированы основные выводы работы.
2. Поляритоны в атомных микроструктурах
Пространственно-периодические структуры атомов, заключенных в оптические резонаторы, открывают новые перспективы для исследования критических явлений в связанных атомно-оптических системах. Создание и развитие методов удержания атомов в оптических решетках различной конфигурации способствовало широкому развитию исследований различных аспектов фазовых переходов для связанных атомно-оптических состояний [16]. В частности, речь идет о квантовых фазовых переходах, происходящих при абсолютном нуле температуры, о переходах типа металл-диэлектрик, о фазовых переходах к сверхизлучению, и др. (см., например, [17-21]).
С другой стороны, современный уровень развития нанотехнологий и нанофотоники делает возможным создание периодических структур на основе цепочек связанных микрорезонаторов — так называемые связанные волноводные оптические резонаторы (coupled resonator optical waveguides), в которых помещены двух- или трех- уровневые атомы [22]. При этом наиболее перспективным выглядит использование микрорезонаторов нового поколения (см. например, [23]) с высоким Q-фактором, для которых условие сильной связи между квантованным излучением и атомами выполняется вплоть до одиночного атома.
В этой связи недавно нами был предложен так называемый поляритонный кристалл (ПК), позволяющий полностью локализовать поляритоны в своей структуре, аналогично «запиранию» фотонов в фотонном кристалле [21].
На рис. 1 представлена модель двумерного поляритонного кристалла, рассматриваемого нами. Она состоит из решетки М туннельно-связанных микрорезонаторов в плоскости XY, образованных в структуре двумерного фотонного кристалла с дефектами. При этом каждый резонатор содержит небольшое (но макроскопическое) количество двухуровневых атомов, взаимодействующих с электромагнитным полем. Кроме того, система атомов снаружи помещена в двумерную оптическую решетку, минимумы которой совпадают с местоположением резонаторов, образованных дефектами фотонного кристалла.
Можно показать, что в приближении так называемой тесной связи (tight-binding approximation) —см., например, [24], полный Гамильтониан Н = Hat + Hph + Hint обсуждаемой атомно-оптической системы может быть представлен в виде Hat, — описывает атомную и оптическую части, соответственно, Hint — Гамильтониан взаимодействия):
(7б)
(7в)
(7а)
(б)
РИС. 1. (а) — схематическое представление геометрии решетки ПК с допустимыми связями между соседними резонаторами, (б) — практическая реализация двумерного поляритонного кристалла. Обозначения: целые числа пх, пу нумеруют ячейки решетки — микрорезонаторы с двухуровневыми атома-
где Мх и Му — число резонаторов по осям X и У
"У
ми; п.
х,у
1, 2,
, Nх,у,
соответственно, I — постоянная решетки
где ЫПхПу
(а,Ь)
'пхпу, аЬ
(^апхпу апхпу + Ьпхпу Ьпхпу^
среднее число атомов в узлах с номерами пх,пу;
ш,
— частота атомов в решетке, коэффициенты связи (ЗпхЬпу (&Пх,пу) характеризуют туннелирование атомов (фотонов) между соседними ячейками и определяются интегралами перекрытия их волновых функций и производных, соответственно; д характеризует атомно-оптическую связь. В (7а), (7б) аПхПуи ЬПхПу — бозе-операторы уничтожения атомов на основном и возбужденном состояниях в ячейке с индексом пх, пу, соответственно; фПхпу —
оператор уничтожения фотонов в той же ячейке. Далее, для простоты мы полагаем все резонаторы идентичными и содержащими одинаковое число атомов N = МПхПу. Исходя из полной симметрии задачи (см. рис. 1а), можно также положить д = дПхПу, ¡Зп^ — /Зпу'^ = [1
С помощью так называемых преобразований Хольштейна—Примакова (см. [25]) гамильтониан (7а)-(7в) можно привести к виду
Hat = Н ^ilny, at - , at + iß(o)) ф]ПхПу фПхпу
nx,ny (8а)
-ß(Ь) (ф1пу Фпх+1пу + ф]ПхПу ФпхПу+1 + Н.С.)
Hph = [иПхПу,рнф]ПхПуФпхпу - а (ф1хПуфПх+1пу + ф]ПхПуФпхпу+1 + Н.С(8б)
Hint = hg Y \Ф]п
хПу Ф^хПу + ФпхПу ФпхПу , (8в)
пх ,Пу
где фп, (п = [пх,пу}) — операторы возмущений в атомной решетке, удовлетворяющие коммутационным соотношениям:
[Ф = $тп, П,т ={пх,пу} . (9)
При выводе (8а)-(8в) учтено, что среднее число атомов на нижнем уровне значительно превосходит их количество в возбужденном состоянии, т. е.
{Ъ]пЪп) < (afnan) . (10)
Данное условие иногда связывают с приближением малой плотности для газа поляритонов (ср. с [21]). ^
Перейдем в (8а)-(8в) к импульсному ^-представлению. Учитывая периодичность структуры поляритонного кристалла, операторы фп = фПхпу и фп = фПхпу можно представить в виде:
где / — вектор решетки. В двумерном случае имеем: к I = кхпх£+купуI, пх = 1, 2,..., Ых, пу = 1, 2,..., Иу, М = Мх х Иу (далее полагаем Их = Иу) — общее число ячеек (резонаторов) в решетке, I — постоянная решетки.
Подставляя (11а,б) в (8а)-(8в) получаем:
Я
П Y1 Фк + Фк + Я (Ф1 Фк + Ф{Фк^ , (12)
где = ^ph(k) и uat = wat(k) определяют дисперсионные соотношения для фотонной и атомной подсистем структуры ПК соответственно, и могут быть записаны в виде
¿atк) = ^1пу, at - ^п1пу, at + 4(3^ - ^ (cOS ^ + COS %) , (13а)
¿Phk) = ш-хПу, Ph2a (cos К +cos Ю ■ (13б)
Выражение (12) в приближении тесной связи представляет многочастичный гамильтониан в импульсном представлении, описывающий периодическую структуру рис. 1, и может
и аПх ~ аПу = а.
быть проанализирован в терминах поляритонов верхней (Е1 —) и нижней (е2 —) дисперсионных ветвей. Он приводится к диагональному виду с помощью линейных преобразований Боголюбова:
Е1,- = ХиФи + Сфъ Е2,— = ХпФи - С^Фь (14)
где Хг и Сг — коэффициенты Хопфилда, определяемые как (см., например, [21])
Xg = -И 1 + / k = , Cf=-F 1--/ k 1 , (15а>6)
1/2 / \ 1/2
ll. 5г \ „II. 5t
Lfc ~~ "
V2 \ . /4п2 + С . , ,
к/ \ V ^ к,
= ш(нк) — к) — отстройка частоты, зависящая от квазиимпульса к.
Операторы , 52£ в (16а,б) характеризуют два типа элементарных возмущений с характерными частотами П1)2(к), которые определяют дисперсионные соотношения, а также зонную структуру поляритонного кристалла. Они имеют вид:
П1>2(к) = ^ [ujat{k) + ujph{k) ± +Vj . (16a,6)
В результате масса поляритонов верхней (1) и нижней (2) ветвей определяется выражением:
2mat mphJ К2 + 4g2 mh2 =-. -j, (17)
(mat + mph)y К2 + 4g2 ± (mat — mph) К
где К = К — 2 {а — [3—независящая от квазиимпульса отстройка, К = 8к|g_0 =
шпхпу, ph nlny, at - ыnlny, at + 4/3(<l)) • При этом выражения mat = Щр, трН = ^^ определяют эффективные массы атомов и фотонов в данной решетке, соответственно.
Физически особый интерес представляют дисперсионные свойства поляритонов нижней ветви. В квазиравновесном состоянии именно она является макроскопически заполненной поляритонами и поэтому — наблюдаемой в эксперименте [8]. Вопросы взаимодействия поляритонов, принадлежащих различным дисперсионным ветвям, которые могут происходить при достаточно высоких температурах и больших атомных нелинейностях, подробно были проанализированы в работе [26], и здесь мы их не обсуждаем. Таким образом, далее будут обсуждаться свойства поляритонов только нижней ветви (индекс «2» у соответствующих величин мы далее для краткости опускаем).
На рис. 2 представлена трехмерная дисперсионная зависимость характерной частоты ВД = ^2 к) для нижней поляритонной ветви в первой зоне Бриллюэна (ЗБ). Принципиальным моментом здесь является наличие минимума при к = 0.
При малых значений квазиимпульса, когда оказывается справедливым условие к£ ^ 1, на основе (12), (16а,б) получаем:
Ьк2
Ньв = П^(к)Е^, ВД«—. (18а,б)
к
Выражения (18а,б) описывают идеальный газ свободных частиц — поляритонов , находящихся на дне дисперсионной поверхности 0,(к), приведенной на рис. 2. Именно данное обстоятельство является принципиальным при обсуждении вопроса о конденсации
поляритонов в состояние с к как:
0. Химический потенциал ц такого газа определяется [13]
ц = кв Т ln(1 - е~П2Хт) = к в Т ln(1 - e-Td/T),
(19)
где Тл — температура квантового вырождения газа, определенная в (2), \т = (2пН2/тквТ)1/2 длина тепловой волны де Бройля. Температура в (18а,б) определяется условием, когда тепловая длина волны Ат имеет тот же порядок величины, что и среднее расстояние между частицами п2 1/2 (см. таблицу 1).
РИС. 2. Дисперсионная зависимость П (к) для поляритонов нижней ветви. Параметры системы: коэффициент атомно-оптической связи д/2п ~ 12.2 ГГц, число атомов в каждом резонаторе N = 100, число резонаторов М = 100, постоянная решетки I = 3 мкм и отстройка частоты А = 0
Приведем численные оценки для наблюдения рассматриваемых явлений на примере взаимодействия двухуровневых атомов рубидия с электромагнитным полем в структуре ПК, средняя частота перехода между уровнями D-линий которых составляет 382 ТГц. Например, в случае резонансного атомно-оптического взаимодействия при А = 0, масса поляритонов может быть найдена из выражения (17) и имеет вид:
2т,
m =--. (20)
1 + mph/mat
В эксперименте, когда выполняется условие mph/mat ^ 1 (а ^ ¡3), и масса поляритонов нижней ветви достаточно мала, а именно, m ~ 2mph — 5.6 х 10-36 кг при массе атомов рубидия mat — 1.44 х 10-25 кг. Поэтому температура вырождения двумерного газа поляритонов Т^ может быть достаточно высокой и составляет приблизительно 448 K при концентрации
двумерного газа поляритонов в решетке п2 ~ 0.5 х 108 см-2 и тепловой волне де Бройля Ат ~ 1.4 мкм (ср. с таблицей 1).
Из (19) следует, что, в условиях сильного вырождения газа поляритонов, наступающего при п2ХГ ^ 1 (Т ^ То), химический потенциал газа ц является аналитической функцией от температуры и возрастает при Т ^ 0, оставаясь при этом отрицательным. В результате конденсация двумерного газа поляритонов отсутствует при конечной (ненулевой) температуре. Вместе с тем, как это уже было отмечено выше, уже при температуре Ткт (см. (6)), в двумерном бозе-газе происходит фазовый переход Костерлица—Тоулесса к состоянию сверхтекучести. В этом случае, с учетом даже слабого взаимодействия между поляритонами (не идеальный газ) дополнительный параметр азсаг, который является длиной рассеяния и зависит от потенциала взаимодействия; он влияет на энергетический спектр над основным состоянием, который становится фононным, ср. с [13]. При этом эффективный процесс релаксации поляритонов на дно «дисперсионной ямы» может быть связан с весьма интенсивным процессом поляритон-поляритонного взаимодействия, обсуждаемого в целом ряде работ (в основном для задач с полупроводниковыми микрорезонаторами [27]).
Следует, однако, иметь в виду, что пространственно-периодическая структура ПК, приведенная на рис. 1, имеет большое время когерентности только при относительно низких температурах, когда удержание атомов в оптической решетке становится возможным. Именно большое время когерентности при формировании когерентного бозе-конденсата поляритонов в структуре ПК представляет интерес для задачи распределенной записи и хранения квантовой оптической информации, рассмотренной нами в [28].
3. Предельная групповая скорость поляритонов в двухуровневой атомной среде
В связи с задачами квантовой информации в квантовой и атомной оптике с недавних пор большое внимание уделяется динамике распространения оптических волновых пакетов, в том числе эффекту замедления («остановки») света, заключающемуся в существенном уменьшении наблюдаемой групповой скорости светового импульса, распространяющегося в атомной среде [29]. Особый интерес к таким динамическим задачам вызван в первую очередь возможностью создания квантовых устройств для копирования и хранения информации, предложенных, например, в работе [28].
Эффект «остановки» света был теоретически предсказан и экспериментально обнаружен для широкого круга задач, для которых обеспечивается макроскопическая когерентность двух- или трехуровневой атомной системы при ее взаимодействии с квантованным световым полем — см. например, [30]. В этом аспекте следует также отметить и среды типа фотонных кристаллов, а также близкие к ним системы связанных резонаторов [31,32]. При этом существенная редукция групповой скорости светового излучения в таких системах как раз и является физической основой для устройств хранения информации. Вместе с тем, остается открытым вопрос о нижнем, по сути квантовом, пределе минимальной скорости распространения светового импульса в такой среде. Именно в задачах копирования (записи) и хранения квантовой информации эта проблема является принципиальной в аспекте выяснения перспектив практического использования подобных систем. Само изменение величины групповой скорости распространяющегося в таких системах светового импульса непосредственно связано с формированием поляритонов, образованных линейной суперпозицией фотонов внешнего (пробного) светового поля и макроскопических (когерентных) возмущений среды. В таких связанных состояниях поле + среда происходит эффективная перекачка энергии между ними, и, при определенных условиях, становится возможным случай, когда вся энергия квантованного светового поля передается возмущениям среды.
Следует, однако, иметь в виду, что в условиях насыщения атомно-оптической системы ее поляритонная модель теряет физический смысл. В этой связи интересной представляется работа [33], авторы которой связывают экспериментальные результаты работы [34], касающиеся квантовой записи и хранения оптической информации, именно со свойством среды как насыщаемого поглотителя.
Для структуры ПК, изображенной на рис. 1, выполняется предел поляритонного газа малой плотности, так что атомно-оптическая система далека от насыщения. Именно этот факт позволяет использовать поляритоны нижней дисперсионной ветви для целей когерентного копирования (записи) квантовой информации, содержащейся в световом импульсе, распространяющимся в такой атомной среде, на основе управляемого изменения его групповой скорости.
Действительно, групповая скорость поляритонов в структуре ПК, определяется выражением
= Пу/ъш2{кхг) + $т2{куг) 21таЬтрН
+ «пь-пЛ
(Ч + 4»2) !
(21)
в котором мы пока не учитываем процессы поляритон-поляритонного рассеяния и декоге-рентности.
Зависимости приведенной (к скорости света в вакууме) групповой скорости поляри-тонов для первой ЗБ показаны на рис. 3. Значение скорости V = 0 достигается на границах ЗБ при кх,у£ = пж, п = 0, ±1,.... Это означает полную локализацию поляритонов внутри зоны, что в свою очередь создает возможность формирования «медленного» света на основе структуры ПК. В идеальном случае для малых значений квазиимпульсов (таких,
Нк
что кХуУ£ 1), зависимость групповой скорости поляритонов линейна, т.е. у (к) ~ — и
соответствует параболическому закону дисперсии (18а,б). Таким образом, в пределе к ^ 0 групповая скорость поляритонов v(k) также стремится к нулю.
Однако, на практике для достижения предельно низких значений групповой скорости поляритонов, которые и определяют наблюдаемую скорость распространения светового импульса в атомной среде, принципиальным является выполнение ряда условий, определяемых как поляритон-поляритонным рассеянием, так и неравновесными (релаксационными) процессами. Детальный анализ влияния последних можно найти в [35], где показано, что они сводятся к неравенствам
5П(к) < П(к). (22а)
Нк2
5П(к) <-, (226)
2т
в которых величина 60,(к) (мнимая часть частоты) характеризует релаксационные процессы в газе поляритонов нижней дисперсионной ветви.
Физически условие (22а) означает наличие узких спектральных линий в широком диапазоне значений отстройки Д. Именно в этом пределе можно полагать справедливым выражение (21) для групповой скорости поляритонов, полученное на основе обобщения на двумерный случай известного определения V = д&(к)/дк — ср. с [36]. Если условие (21) не выполняется, то групповая скорость оптического волнового пакета может быть определена через скорость переноса энергии электромагнитного излучения в среде; данная проблема была подробно исследована в работе [37]. Однако, в нашем случае при учете скорости спонтанных переходов для атомов рубидия, составляющей ~ 37 МГц, неравенство (22а)
РИС. 3. Приведенная групповая скорость поляритонов v/c как функция от нормированного квазиимпульса кх1 при куI = п (сплошная линия) и куI = 0 (пунктирная линия); с — скорость света в вакууме. Остальные параметры те же, что и на рис. 2. На вставке показана нормированная скорость V = v0 сверхтекучих поляритонов с учетом их взаимодействия, определяемого параметром связи и0/2кН ~ 61 МГц, взятого в единицах частоты, и среднего числа поляритонов из расчета на один резонатор про1 = 10
выполняется с большим запасом (характерная частота П(к) находится в оптическом диапазоне частот —см. рис. 2).
Соотношению (22б) можно придать несколько иной физический смысл, связав его с неопределенностью волновых векторов 5к поляритонов (ср. с [38]). Действительно, стандартная «неопределенность» частоты поляритонов 5П(к) приводит к неопределенности их волновых векторов 5к, которая и определяет мнимую часть волнового вектора. Она может быть найдена из выражения 8к = . Такая неопределенность волновых векторов приводит к нежелательному уширению функции распределения для значений импульсов поляритонов. В этой связи волновые векторы к можно полагать «пригодными» для квантового рассмотрения, только если выполняется условие 5к << к, что, по сути, аналогичное (22б), которое фиксирует приемлемые значения параметров релаксации атомной (уаь) и оптической (урь) подсистем для пространственно-периодической структуры на рис. 1, а также области значений волновых векторов к, для которых поляритоны сохраняют когерентные свойства и могут быть описаны квантовым образом, в частности, с помощью когерентных волновых пакетов.
Проиллюстрируем сказанное на примере преобразования групповой скорости по-ляритонов, соответствующих нижней ветви дисперсионной кривой вблизи ее минимума к = 0, изображенного на рис. 4.
Пусть в начальный момент времени значение отстройки Д соответствует полностью
М;2
фотоноподобным поляритонам нижней ветви с массой т ~ ш7,/, и частотой 0(7,- ) ~ -.
2три
Тогда в области высоких (оптических) частот П^ на рис. 4 для спектра входного светового импульса имеем 5 = Брь ~ {ш-г,к2/1т н)2+т • Если неравенства (22а), (226) выполняются, то такие поляритоны представляют квантовый (когерентный) волновой пакет со средним значением волнового вектора .
РИС. 4. Дисперсионные зависимости Q(k) фотоноподобных (Qph) и эксито-ноподобных (Qat ) поляритонов нижней ветви
Далее, изменяя Д, переходим к другому пределу — экситоноподобным поляритонам
hk2
(нижняя ветвь зависимости на рис. 4) с массой m ~ mat и частотой Q(k) ~ -. В
2mat
этом пределе спектр исходного оптического волнового пакета существенно сужается, преобразуясь в спектр экситонов (двухуровневых возмущений атомной среды), определяемый
выражением S = S(!/ ~ - -—-g- в области низких частот Поэтому, как
(ш - hk2/2mat) + lit
это видно из рис. 4, уширение спектра может оказаться существенным, и условие (22б) нарушается, что означает потерю когерентности (дефазировку (dephasing)) распространяющегося волнового пакета.
Разумеется, можно выбрать другую группу поляритонов с большим значением волнового вектора fc(2), при котором условие (22б) выполняется для обоих рассматриваемых
по величине Д случаев. Однако, при этом также растет и нижний уровень для значений групповой скорости поляритонов.
В частности, в предельном случае 5к — к неопределенности волновых векторов 5к
для атомоподобных поляритонов нижней ветви дисперсионной кривой, когда выполняются
2
"МрИ У
неравенства- -г— <<с; 1, для групповой скорости у поляритонов и для неопределен-
та1 Д2
ности их волновых векторов 8к получаем соотношения:
у = Пкд2 ^ = ^а1трНК2 ^
тр1г А2 Щ
Для атомов рубидия при значении квазиимпульса к — 8 х 106 1/м, предельное значение групповой скорости равно V2 — 4.6 м/с, что соответствует минимально возможной скорости атомоподобных поляритонов в среде с учетом сохранения их когерентности.
Обсудим теперь ограничения, возникающие из-за влияния поляритон-поляритонного рассеяния на групповую скорость поляритонов нижней дисперсионной ветви. Оно может быть учтено, например, в рамках подхода Боголюбова для слабо взаимодействующего газа поляритонов (ср. с [5,11]). В этом случае минимальная групповая скорость поляритонов ограничена снизу скоростью элементарных возмущений, а именно, скоростью звука, определяемого выражением
V = (24)
т
где щ = —( 'и ХЦ (X,, = Х;- (), С0 = С^=0) представляет нелинейный параметр, характери-
зующии взаимодеиствие поляритонов, про1 = --безразмерная плотность поляритонов
в рассматриваемой структуре.
На вставке к рис. 3 представлена численная зависимость нормированной групповой скорости поляритонов V для рассматриваемой задачи взаимодействия светового излучения с атомами рубидия при учете поляритон-поляритонного рассеяния. В частности, видно, что поведение величины V демонстрирует сверхтекучесть газа поляритонов в структуре ПК. При этом типичные значения V для элементарных возмущений в газе поляритонов существенно больше таковых для ансамбля атомов, находящихся в состоянии БЭК (ср. с [1]).
4. Физические принципы хранения квантовой оптической информации на основе поляритонов в структуре ПК
Рассмотрим проблему физического копирования (записи), хранения и считывания квантовой информации на основе структуры ПК. Под копированием квантового состояния оптической системы здесь подразумевается по возможности полное его отображение на состояние атомной среды при условии также сохранения формы распространяющегося волнового пакета. Физически процедура копирования, эквивалентна определенному измерению, проводимому над квантовой системой, в результате которого исходное (копируемое) состояние светового излучения перестает существовать.
В условиях сохранения когерентности поляритонов во всей области изменения приведенного параметра Д /д поляритоны нижней дисперсионной ветви с соответствующим законом дисперсии (18) представим как волновой пакет, описываемый при Ь = 0 волновой
функцией вида (здесь для простоты (см. рис. 5) мы рассматриваем одномерную задачу):
(25)
Ф (ж, i = 0) = exp
2
х р0х + г-
2 Р ' ' К
где f характеризует ширину оптического волнового пакета, т.е. определяется шириной светового пучка на входе в среду. На основе решения уравнения Шредингера получаем следующее выражение для плотности вероятности:
{х - Уро1£)2
2 1 1ф(М)1 =7rl/4/-/(t)exp
f 2(/т
(26)
где f(t) = \J 1 + t2/r2 — безразмерная ширина волнового пакета в момент времени /.,
тъ = mf2/К — эффективное время его расплывания и vo = —- = — — групповая скорость
m m
центра волнового пакета вдоль оси х.
Хранение квантовой оптической информации с помощью структуры ПК предполагает трехшаговый физический алгоритм реализации — см. рис. 5а.
Во-первых (стадия 1), для записи квантовой информации выбираем такую отрицательную отстройку Д, чтобы выполнялось условие Д ^ —2 \д\; на рис. 5а —это интервал времени 0 ^ Т ^ 0.25. В данном случае поляритоны нижней дисперсионной ветви являются фотоноподобными, т. е. m ~ mph, и распространяются благодаря туннелированию фотонов в решетке ПК. Волновой пакет при этом распространяется со скоростью v, определяемой выражением
hk
Y=J^ = 2af2k, (27)
mph
и смещается вдоль оси х — из точки х = 0 в точку х = vot - рис. 5б. В этом случае для характерного времени расплывания имеем оценку ть = mphf 2/h ~ 0.27 не при ширине входного пучка f ~ 10_4 м.
Во-вторых (стадия 2), для отображения оптической информации на когерентные атомные возмущения среды мы переключаем отстройку, сделав ее положительной, т. е. Д ^ 2 \д\. В этом пределе поляритоны нижней ветви становятся атомоподобными. При
mPh $ „ 92 ы ¡Т\
-= — их групповая скорость описывается выражением (23а,б).
mat а Д2
2
тт с 92 „ mph „ ,
При выполнении более строгого условия -г— <<с;-<<с; 1 групповая скорость поля-
Д2 mat
ритонов нижней ветви будет определяться выражением
hk
v = — = 2 f3f2k, (28)
mat
что соответствует скорости перемещения атомов в решетке. Фактически выражение (28) определяет нижний допустимый предел скорости оптического волнового пакета в структуре ПК. В этом пределе информация, содержащаяся в световом пучке, полностью записывается на атомные возмущения. Тогда для характерного времени расплывания пакета будет ть = mat f2/h ~ 13.7 е при той же ширине входного пучка. Данный временной масштаб (интервал 0.25 ^ Т ^ 0.75 на рис. 5) и определяет максимальное время хранения информации квантовым газом двухуровневых атомов рубидия. Точнее, необходимым условием такой записи информации является выполнение неравенства rrec,stor ^ Ч = matf 2/h, которое означает сохранение формы распространяющегося волнового пакета в течение всего времени отображения и хранения квантового состояния, где Tstor определяет время хранения информации атомной системой.
Наконец, в-третьих, для восстановления оптической информации на выходе из среды через интервал времени Tstor мы должны теперь сделать поляритоны фотоноподобными, переключая отстройку Д в обратном порядке. На рис.5а характерное время обратного переключения волнового пакета rretr определяется промежутком времени 0.75 ^ Т ^ 0.8. На выходе из среды (при Т =1) опять имеем оптический волновой пакет, но смещенный вдоль оси х в плоскости, перпендикулярной оси резонатора.
В данной работе мы не приводим оценку качества (надежности) сохраняемой информации. В [35] нами рассчитана оценка параметра надежности (fidelity) сохраняемой информации при учете только изменения формы волнового пакета. Однако, в общем случае этого недостаточно. Необходимым является также анализ преобразования квантового состояния светового поля на основе структуры ПК с учетом ее декогерентности. Анализ данной задачи представляет самостоятельный интерес и выходит за рамки настоящей работы. В этой связи лишь отметим, что в реальной экспериментальной ситуации интервал времени записи, хранения и считывания информации все же ограничен временем декоге-рентности структуры ПК. При использовании ультрахолодных атомов, система обладает достаточно большим временем макроскопической когерентности в целом, составляющем десятки микросекунд по экспериментальным данным [30].
5. Заключение
В работе с единых позиций рассмотрены различные аспекты фазовых переходов для низкоразмерного (двумерного) газа поляритонов формирующихся в полупроводниковых и атомных средах. В этой связи предложены решеточные модели для поляритонов в пространственно-периодической структуре — поляритонном кристалле, образованном двумерной решеткой резонаторов с помещенными в них макроскопически небольшими порциями двухуровневых атомов, эффективно взаимодействующих с электромагнитным полем в рамках приближения сильной (тесной) связи между соседними резонаторами. Показано, что структура ПК позволяет полностью локализовать поляритоны нижней дисперсионной ветви, что, во- первых, может быть использовано для получения бозе-конденсации по-ляритонов, и, во-вторых, приводит к существенному замедлению светового импульса в рассматриваемой среде. В частности, показано, что как квазиконденсация, так и истинная конденсация двумерного газа поляритонов может быть в принципе достигнута при относительно высоких температурах при взаимодействии двухуровневых атомов с квантованным световым излучением в условиях сильной связи между атомами и полем. Кроме того, проанализировано влияние однородного уширения, а также поляритон-поляритонного рассеяния на процесс конденсации поляритонов и на величину их групповой скорости. Когерентные свойства ансамбля поляритонов нами рассматриваются в аспекте реализации возможности пространственно-распределенной квантовой записи, хранения и считывания информации, ассоциированной с распространяющимся оптическим волновым пакетом. В этой связи предложены физические принципы хранения квантовой оптической информации на основе когерентных поляритонов в структуре ПК. Показано, что запись квантовой информации в среде может быть осуществлена адиабатическим образом на основе управления групповой скоростью когерентного волнового пакета (поляритонов) — путем замедления («остановки») наблюдаемого распространения светового импульса, энергия которого преобразуется в когерентные атомные возмущения. При этом расплывание волнового пакета происходит на временах, значительно превышающих характерное время когерентности конденсата атомов, что позволяет восстановить на выходе волновой пакет практически без изменения его формы.
(б)
РИС. 5. Распространение одномерного поляритонного волнового пакета в структуре ПК; (а) — зависимость параметра И = Д/2 \д\ от нормированного времени Т = Ш/трь,$2 и (б) — зависимость огибающей пакета (плотности вероятности) 5 = \Ф(ж, ¿)\2/\Ф(0, 0)\2 от Т и нормированной координаты X = х/! для комбинации параметров = 10. На рисунке обозначены стадии записи информации (1), ее хранения — (2) , и считывания (восстановления) оптической информации — (3).
Работа поддержана грантами РФФИ ННИО_а-09-02-91350 и РТ_оми-10-02-13300, а также проектами федеральных программ Минобрнауки РФ.
Литература
[1] Кеттерле B. Когда атомы ведут себя как волны. Бозе-эйнштейновская конденсация и атомный лазер. УФН, 2003, 173 (14), 1339-1358.
[2] Bagayev S.N., Vasil'ev V.V., Egorov V.S., Lebedev V.N., Mekhov I.B., Moroshkin P.V., Fedorov A.N., Chekhonin I.A. Coherent light sources under strong field-matter coupling in an optically dense resonant medium without population inversion. Laser Physics, 2005, 15(7), 975-982.
[3] Ораевский А.Н. Бозе-конденсаты с точки зрения лазерной физики. Квантовая электроника, 2001, 31, 1038-1057.
[4] Скалли М.О., Зубайри М.С., Квантовая оптика, Москва, Физматлит, 2003, 512с
[5] Kavokin A, Malpuech G., Laussy F.P. Polariton laser and polariton superfluidity in microcavities. Phys. Lett. A, 2003, 306, 187-199.
[6] Klaers, J., Vewinger, F., Weitz, M. Thermalization of a two-dimensional photonic gas in a 'white-wall' photon box. Nature Physics, 2010, 6, 512-515
[7] Deng H., Weihs G., Snoke D., Bloch J., Yamamoto Y. Polariton lasing vs. photon lasing in a semiconductor microcavity. PNAS, 2003, 100 (24), 15318-15323.
[8] Kasprzak J., Richard M., Kundermann S., Baas A., Jeambrun P., Keeling J.M.J., Marchetti F.M., Szymanska M.H., Andre R., Staehli J.L., Savona V., Littlewood P.B., Deveaud B., Dang Le Si. Bose-Einstein condensation of exciton polaritons. Nature, 2006, 443, 409-414.
[9] Balili R., Hartwell V., Snoke D., Pfeiffer L., West K. Bose-Einstein Condensation of microcavity polaritons in a trap. Science, 2007, 316, 1007-1010.
[10] Amo A., Lefrere J., Pigeon S., Adrados C., Ciuti C., Carusotto I., Houdre R., Giacobino E., Bramati A. Superfluidity of polaritons in semiconductor microcavities. Nature Physics, 2009, 5(11), 805-810.
[11] Utsunomiya S., Tian L., Roumpos G., Lai C.W., Kumada N., Fujisawa T., Kuwata-Gonokami M., Loffler A., Hofling S., Forchel A., Yamamoto Y. Observation of Bogoliubov excitations in exciton-polariton condensates. Nature Physics, 2008, 4(12), 700-705.
[12] Bagnato V., Kleppner D.K. Bose-Einstein condensation in low-dimensional traps. Phys. Rev. A, 1991, 44 (11), 7439-7441.
[13] Petrov D.S., Gangardt G.M., Shlyapnikov G.V. Low-dimensional trapped gases. J.Phys. IV France, 2004, 116, 3-45 .
[14] Аверченко В.А., Алоджанц A.n., Аракелян С.М., Багаев С.Н., Виноградов Е.А., Егоров В.С., Столяров А.И., Чехонин И.А. Высокотемпературная бозе-эйнштейновская конденсация поляритонов: реализация в условиях внутрирезонаторной лазерной накачки вещества. Квантовая Электроника, 2006, 36(6), 532-538.
[15] Chestnov I. Yu., Alodjants A. P., Arakelian S. M., Nipper J., Vogl U., Vewinger F., Weitz M. Thermalization of coupled atom-light states in the presence of optical collisions. Phys. Rev. A., 2010, 81(5), 053843-1-13.
[16] Bloch I. Exploring quantum matter with ultracold atoms in optical lattices. J. Phys. B: At. Mol. Opt. Phys., 2005, 38, S629-S643; Jaksch D., Zoller P. The cold atom Hubbard toolbox. Annals of Phys., 2005, 315(1), 52-79.
[17] Silver A.O., Hohenadler M., Bhaseen M. J., Simons B.D. Bose-Hubbard models coupled to cavity light fields. Phys. Rev. A, 2010, 81(2), 023617-1-15.
[18] Hartmann M.J., Brandao F., Plenio M.B. Strongly interacting polaritons in coupled arrays of cavities. Nature Physics, 2006, 2(14), 849-855; Greentree A.D., Tahan C., Cole J.H., Hollenberg L.C.L. Quantum phase transitions of light. Nature Physics, 2006, 2(14), 856-861.
[19] Soi-Chan Lei, Ray-Kuang Lee. Quantum phase transitions of light in the Dicke-Bose-Hubbard model. Phys. Rev. A, 2008, 77(3), 033827-1-8.
[20] Fleischhauer M., Otterbach J., Unanyan R.G. Bose-Einstein Condensation of Stationary-Light Polaritons. Phys. Rev. Lett., 2008, 101(16), 163601-01-04.
[21] Alodjants A.P., Barinov I.O., Arakelian S.M. Strongly localized polaritons in an array of trapped two-level atoms interacting with a light field. J. Phys. B: At. Mol. Opt. Phys., 2010, 43(12), 095502-1-10.
[22] Painter O., Vuckovic J., Scherer A. Defect modes of a two-dimensional photonic crystal in an optically thin dielectric slab. J. Opt. Soc. Am. B, 1999, 16(2), 275-285; Lee R.K., Painter O., Kitzke B., Scherer A., Yariv
A. Emission properties of a defect cavity in a two-dimensional photonic bandgap crystal slab. J. Opt. Soc. Am.
B, 2000, 17(4), 629-633.
[23] Braginsky V.B., Gorodetsky M.L., Ilchenko V.S. Quality factor and nonlinear properties of optical whispering-gallery modes. Phys. Lett. A, 137,1989,393-397; Vernooy D.W., Kimble H.J. Quantum structure and dynamics for atom galleries. Phys. Rev. A, 1997, 55(2), 1239-1261; Buck J.R., Kimble H.J. Optimal sizes of dielectric microspheres for cavity QED with strong coupling. Phys. Rev. A, 2003, 67(3), 033806-1-033806-11; Deych L.I., Roslyak O. Photonic band mixing in linear chains of optically coupled microspheres. Phys. Rev. E, 2006, 73(3), 036606-1-036606-12.
[24] Smerzi A., Trombettoni A., Kevrekidis P.G., Bishop A.R. Dynamical Superfluid-insulator transition in a chain of weakly coupled Bose-Einstein Condensates. Phys. Rev. Lett., 2002, 89(17), 170402-1-4.
[25] Holstein T., Primakoff H. Field Dependence of the intrinsic domain magnetization of a ferromagnet. Phys. Rev., 1940, 58(14), 1098-1113.
[26] Alodjants A.P., Arakelian S.M., Bagayev S.N., Egorov V.S., Leksin A.Yu. Josephson dynamics for coupled polariton modes under the atom-field interaction in the cavity. Applied Physics B, 2007, 89(1), 81-89.
[27] Karr J. Ph., Baas A., Giacobino E. Twin polaritons in semiconductor microcavities. Phys. Rev. A, 2004, 69(6), 063807-1-12.
[28] Alodjants A.P., Arakelian S.M., Bagayev S.N., Chekhonin I.A., Egorov V.S. Quantum cloning in coupled states of an optical field and an atomic ensemble by means of quasi-condensation of polaritons. J. of Rus. Laser Research, 2006, 27(5), 482-491.
[29] Fleischauer M., Lukin M.D. Quantum memory for photons: dark-state polaritons. Phys. Rev. A, 2002, 65(2), 022314-1-12 .
[30] Hau L.N., Harris S.E., Dutton Z., Behroozi C.H. Light speed reduction to 17 meters per second in an ultracold atomic gas. Lett. to Nature, 1999, 397, 594-596.
[31] Gersen H., Karle T.J., Engelen R.J.P., Bogaerts W., Korterik J.P., van Hulst N.F., Krauss T.F., Kuipers L. Real-space observation of ultraslow light in photonic crystal waveguides. Pys. Rev. Lett., 2005, 94(7), 0739031-073903-4.
[32] Yanik M.F., Fan S. Stopping and storing light coherently. Phys. Rev. A, 2005, 71(1), 013803-1-10.
[33] Александров Е.Б., Запасский В.С. В погоне за «медленным светом». УФН, 2004, 174, 1093-1102.
[34] Phillips D. F., Fleischhauer A., Mair A., Walsworth R. L., Lukin M. D. Storage of light in atomic vapor, Phys.Rev.Lett. 2001, 86 (5), 783-787
[35] Alodjants A.P., Arakelian S.M., Leksin A.Yu. Storage of quantum optical information based on the intracavity polaritons under the Bose-Einstein Condensation conditions. Laser Physics,2007, 17(11), 1-9.
[36] Zhang Z. M., Park K., On the group front and group velocity in a dispersive medium upon refraction from a nondispersive medium. J. of Heat Transfer, 2004, 126, 244-249
[37] Loudon R. The propagation of electromagnetic energy through an absorbing dielectric, J. of Phys. A, 1970, 233245.
[38] Agranovich V.M., Litinskaia M, Lidzey D.G. Cavity polaritons in microcavities containing disordered organic semiconductors. Phys. Rev. B., 2003, 67(8), 085311-1-10.