Ground station efficiency for quantum communications
- Authors: Barbyshev K.A.1,2, Duplinsky A.V.1,3, Khmelev A.V.1, Kurochkin V.L.1,4
-
Affiliations:
- QSpace Technologies LLC
- Moscow Power Engineering Institute
- HSE University
- MISIS National University of Science and Technology
- Issue: Vol 88, No 6 (2024)
- Pages: 967-974
- Section: Luminescence and Laser Physics
- URL: https://medbiosci.ru/0367-6765/article/view/276205
- DOI: https://doi.org/10.31857/S0367676524060193
- EDN: https://elibrary.ru/PFOTOB
- ID: 276205
Cite item
Full Text
Abstract
We investigated the possibility of practical implementation of a quantum communication link for encryption key distribution between a Micius satellite and a mobile ground optical station. Thanks to theoretical estimation numerical values of the main parameters of such a communication link are obtained: total loss, key generation rate and secret key length, as well as quantum bit error rate.
Full Text
ВВЕДЕНИЕ
Для обеспечения безопасной передачи и хранения данных на сегодняшний день предложено использовать технологию квантового распределения ключей (КРК). Она направлена на устранение угрозы информационной безопасности, в том числе со стороны активно разрабатываемых в данный момент времени квантовых компьютеров. Технология КРК при обмене сгенерированными истинно случайными битовыми последовательностями между удаленными пользователями опирается на фундаментальные принципы квантовой механики, а не на сложность решения ряда математических задач. Поэтому такой метод обеспечения защиты информации не зависит от вычислительных мощностей злоумышленника и может обеспечить абсолютную секретность передаваемых данных [1, 2].
Главным ограничением современных квантово-криптографических систем является расстояние передачи квантового сигнала, у которого существует предел из-за поглощения фотонов в среде оптического волокна [3]. На данный момент максимальный продемонстрированный диапазон составляет несколько сотен километров [4]. В то же время использование низкоорбитальных космических аппаратов позволяет обойти главное ограничение волоконной квантовой криптографии и соединить абсолютно защищенным каналом связи любые точки планеты [5].
Работы, в которых описано применение научного спутника Micius в качестве промежуточного доверенного узла, демонстрируют возможность создания глобальной сети КРК [6—9]. Однако практическая реализация подобных спутниковых систем по-прежнему сталкивается с рядом инженерно-технических проблем, такими как большие размеры и вес наземных оптических станций, необходимых для приема квантовых сигналов. В данном исследовании мы сообщаем о результатах моделирования процедуры КРК по линии «космос‒земля» на основе сеанса квантовой связи между спутником Micius и портативной приемной наземной станцией. Вес станции составляет 150 кг, что позволяет разместить ее в черте мегаполиса, например, на крыше зданий государственных органов и финансовых учреждений.
МОДЕЛЬ СПУТНИКОВОГО КВАНТОВОГО КАНАЛА СВЯЗИ
Выбор сценария сеанса распределения ключей
Для моделирования выберем из существующих вариантов [10] наиболее оптимальный сценарий КРК с нисходящей линией связи, т. е. когда передатчик, установленный на спутнике, посылает информационные сигналы на Землю. В этом случае основной причиной потерь в канале связи является дифракционная расходимость излучения, которая зависит от длины волны лазерного источника, а также размера и конструкции передающей оптической системы. Отклонение лазерного луча от первоначальной траектории, вызванное влиянием атмосферной турбулентности, происходит в самом конце пути передачи (в приземном слое атмосферы толщиной 20 км), где размер пучка из-за дифракции обычно намного больше, чем его отклонение. Следовательно, нисходящая линия характеризуется меньшим расширением луча по сравнению с восходящей и, таким образом, имеет более высокую эффективность связи.
Характеристики информационного канала также зависят от типа орбиты космического аппарата, которые делятся по орбитальной высоте на три основных класса: низкая околоземная орбита (500—2 000 км), средняя околоземная (2 000—36 000 км) и геостационарная орбиты (> 36 000 км). Для нашего случая оптимальным будет расположение спутника на низкой орбите (на высоте порядка 500 км), преимущества которой заключаются в меньших оптических потерях и меньшей стоимости достижения, а также успешной эксплуатации по сравнению с более высокими орбитами. Для снижения уровня фонового шума канал КРК работает только в ночное время. Поэтому была выбрана круговая солнечно-синхронная полуночная орбита, пересекающая экватор примерно в полночь по местному времени, тем самым оптимизируя время, проведенное в тени Земли.
Оценка длины канала
Представленный в работе метод анализа характеристик канала квантовой связи отличается от взятой за основу модели [11] тем, что учитывает реальные значения эффективности оптической системы и детекторов одиночных фотонов, которые были определены экспериментальным путем. Для того чтобы при помощи математической модели [11] определить эффективность нашего спутникового канала КРК, необходимо изначально рассчитать зависящее от времени расстояние между передатчиком и приемником, а также угол возвышения спутника над горизонтом. Рассмотрим наилучший вариант прохождения спутника — через зенитную точку над станцией наблюдения. Здесь и далее все вычисления выполнены с учетом того, что в качестве передающего узла выступает первый в мире космический аппарат для КРК — китайский спутник Micius.
Угловая скорость спутника на низкой околоземной орбите определяется выражением
, (1)
где g — ускорение свободного падения, R — расстояние от центра Земли до космического аппарата, ωE — угловая скорость Земли, β — угол наклона орбиты космического аппарата.
Численные значения параметров орбиты спутника приведены в табл. 1.
Таблица 1. Параметры орбиты спутника Micius
h, км | RE, км | R, км | β, ° | g, м·с-2 | ωE, рад·с-1 | ω, рад·с-1 |
500 | 6371 | 6871 | 97.3 | 9.8 | 7.3∙10–5 | 120.4∙10–5 |
Зависимость угла возвышения спутника над горизонтом от времени рассчитана, исходя из простых геометрических соображений:
, (2)
где RE — радиус Земли.
Затем была определена зависимость расстояния между спутником и наземной станцией от времени:
, (3)
где h — высота орбиты космического аппарата.
При моделировании необходимо учесть, что сеанс распределения ключей возможен не на всей траектории прохождения спутника над пунктом наблюдения, т. к. наземной станции необходимо некоторое время, чтобы навестись на объект слежки.
На рис. 1 показаны зависимости угла места спутника и расстояния до него от времени, полученные для зенитного прохода. Результаты представлены для того времени, когда спутник находился в поле зрения наземной станции, т. е. когда угол возвышения над горизонтом составлял свыше 20°.
Рис. 1. Угол возвышения спутника (1, сплошная кривая) и расстояние до него (2, штриховая кривая)
Из графиков, представленных на рис. 1, видно, что расстояние между спутником и приемником уменьшается с увеличением угла места. Строго говоря, лишь проход спутника через зенитную точку обеспечивает максимальную продолжительность квантовой связи и наименьшую длину канала. В нашем случае длительность сеанса составила примерно 270 с, а расстояние «спутник–Земля» варьировалось от 500 до 1200 км соответственно. Для более наглядного представления полученные зависимости симметрично выровнены относительно наибольшего угла восхождения (зенитного положения), что соответствует нулевому времени.
Оценка потерь в канале связи
Потери в канале квантовой связи между космическим аппаратом и наземной приемной станцией динамически изменяются в ходе сеанса КРК, что существенно отличает этот метод от распределения ключей по оптическому волокну. При оценке уровня потерь в линии спутниковой квантовой связи необходимо учитывать следующие факторы:
1) дифракцию луча;
2) атмосферное поглощение и рассеяние;
3) эффективность оптического тракта приемного модуля;
4) квантовую эффективность детекторов одиночных фотонов.
Как правило, наибольший вклад в суммарное ослабление сигнала вносят дифракционные потери, которые зависят от изменяющейся протяженности канала. Кроме того, меняющийся со временем угол места приводит к различной эффективной толщине атмосферы, что в свою очередь отражается на коэффициенте пропускания света.
Потери из-за дифракции были рассчитаны по формуле
, (4)
где e — коэффициент обструкции приемного телескопа, Dt — диаметр приемного телескопа, γ — угловая расходимость лазерного источника, d — длина канала связи.
Следующее немаловажное влияние на ослабление оптического сигнала оказывают потери на прохождение излучения сквозь толщу атмосферы. Для точной оценки необходимо учесть явления поглощения и рассеяния.
Рассеяние лазерного излучения на осадках различных типов таких, как дождь, снег, туман, пыль, дымка, является одной из основных причин ослабления сигнала. Присутствие таких частиц в воздухе приводит к угловому перераспределению светового потока и уменьшению дальности его распространения. Каждый тип рассеяния определяется физическим размером частиц по отношению к длине волны лазерного излучения.
Рэлеевское рассеяние — это упругое рассеяние излучения, связанное с молекулами воздуха, которые значительно меньше по размеру, чем длина волны. Коэффициент молекулярного рассеяния рассчитывается согласно формуле
, (5)
где n — показатель преломления среды, N — число молекул в единице объема, λ — длина волны излучения, δ — фактор деполяризации рассеянного излучения.
По формуле (5) видно, что искомый коэффициент обратно пропорционально зависит от четвертой степени длины волны. Это означает, что при использовании информационного сигнала с λ = 850 нм можно пренебречь энергетическими потерями за счет рэлеевского рассеяния.
Еще один тип рассеяния происходит на частицах аэрозоля, которые больше, чем молекулы, но достаточно малы, чтобы оставаться в воздухе на протяжении длительного периода времени — это рассеяние Ми. Размеры аэрозольных частиц в атмосфере колеблются в диапазоне от 2 нм до 100 мкм. Для оценки затухания сигнала, вызванного аэрозолями, необходимо учитывать их состав, концентрацию и распределение. Из-за значительных экспериментальных трудностей, связанных с определением всех этих параметров, были созданы модели, описывающие аэрозольные условия в зависимости от метеорологических или локальных параметров окружающей среды. Одна из таких моделей, предложенная авторами работы [12], описывает затухание сигнала в атмосфере выражением
, (6)
где V — метеорологическая оптическая дальность, q — параметр, зависящий от метеорологической оптической дальности следующим образом:
. (7)
Чтобы учесть тот факт, что свет проходит большее расстояние через атмосферу у горизонта, чем в зените, воспользуемся соотношением Янга‒Ирвина, которое авторы вывели в работе [13] для связи воздушной массы с углом места. Под воздушной массой в данном случае понимают длину пути через атмосферу. Уравнение
(8)
предполагает использование истинного зенитного угла, т. е. угла, отсчитываемого от горизонта к наблюдаемому объекту в отсутствие влияния рефракции атмосферы.
С учетом всего вышеизложенного итоговая эффективность канала была определена с помощью выражения
(9)
где ηopt — оптическая эффективность приемного модуля, ηdet — квантовая эффективность детекторов одиночных фотонов.
В табл. 2 представлены все необходимые параметры для оценки суммарных потерь при использовании приемной оптической станции с апертурой 300 мм.
Таблица 2. Параметры для оценки суммарных потерь в линии связи
Апертура, Dt, м | Обструкция, e | Расходимость, γ, рад | Коэф. экстинкции атмосферы, α | Эффективность оптики, ηopt | Эффективность детекторов, ηdet |
0.3 | 0.81 | 10–5 | 0.1—1 | 0.5 | 0.62 |
Для более наглядного представления информации выразим потери, рассчитанные по формуле (9), в дБ:
. (10)
По графикам, представленным на рис. 2 можно определить, что разница потерь для положения спутника на краю траектории (θel = 20°, t = 135 с) и в зените (θel = 90°, t = 0 с) составляет 9 дБ для наилучших погодных условий и 15 дБ для наихудших.
Рис. 2. Суммарные потери в линии связи в случае отличных погодных условий (1, сплошная кривая), хороших (2, штриховая кривая), средних (3, штрихпунктирная кривая), плохих (4, пунктирная кривая)
Передача однофотонного сигнала по линии связи «спутник–Земля»
Для оценки ключевых параметров сеанса КРК была использована модель, представленная в работе [14]. Согласно этой модели, передача i-фотонного состояния по каналу связи, описывается выражением
, (11)
где i = 0, 1, 2…
Согласно теории, число фотонов в каждом слабом когерентном импульсе лазерного источника имеет распределение Пуассона [6]. В нашем случае при реализации протокола BB84 Decoy State сигнальные импульсы имеют интенсивность сигнала μ = 0.8, состояния-приманки — ν = 0.1 и состояния вакуума — λ = 0, которые испускаются случайным образом с вероятностями ps = 0.5, pd = 0.25 и pν = 0.25 соответственно.
Величина фонового шума Y0, который представляет собой темновые срабатывания детектора, была определена экспериментальным путем при сопровождении спутника с выключенным передатчиком КРК. Уровень шума был постоянным на протяжении всего сеанса связи и составил 250 кликов/с для наземной станции с апертурой 300 мм.
Кроме того, ошибка детектирования (вероятность попадания фотона в ошибочный детектор) edet, вызванная точностью измерений приемника, не зависела от расстояния и составила порядка 1 %. Значения, использованные для моделирования сеанса КРК, перечислены в табл. 3.
Таблица 3. Параметры для моделирования сеанса распределения ключей
Тип импульса | Интенсивность, фотон·имп-1 | Вероятность отправки, p | Частота импульсов, f, МГц | Фоновый шум, Y0, клик·с-1 | Ошибка детектирования, edet |
Сигнал, μ | 0.8 | 0.50 | 100 | 250 | 0.01 |
Приманка, ν | 0.1 | 0.25 | |||
Вакуум, λ | 0 | 0.25 |
Скорость генерации просеянного ключа Rsift и коэффициент квантовых битовых ошибок QBER были рассчитаны с учетом параметров оптической линии связи, а также технических характеристик устройств Алисы и Боба по оценочным выражениям, предложенным в работе [14]:
, (12)
, (13)
где f — частота следования импульсов, e0 — ошибка, равная 0.5 для протокола BB84 Decoy State, η1(t) — эффективность передачи однофотонного состояния.
На рис. 3 представлены расчетные значения Rsift и QBER для зенитного пролета спутника в зависимости от погодных условий. Видно, что скорость распределения просеянного ключа и коэффициент ошибок имеют сильную зависимость от угла места и изменяются неравномерно по ходу проведения сеанса КРК.
Рис. 3. Скорость генерации просеянного ключа (а) и коэффициент квантовых битовых ошибок (б) в случае отличных погодных условий (1, сплошная кривая), хороших (2, штриховая кривая), средних (3, штрихпунктирная кривая), плохих (4, пунктирная кривая)
Оценка длины ключа
Длину просеянного ключа lsift после процедуры сверки базисов можно определить, исходя из информации о скорости его распределения между абонентами, а также длительности сеанса связи:
, (14)
где t1 и t2 — время начала и окончания сеанса КРК, соответственно.
Финальный же ключ пользователи получают из просеянного после выполнения ряда процедур и математических операций по коррекции ошибок и усилению секретности. Нижний предел длины секретного ключа lsec оценивается выражением
, (15)
где N — общее количество импульсов, излучаемых источником, ps — вероятность испускания сигнального импульса, q — ошибка выбора базиса (для протокола BB84 Decoy State q = 0.5), Q1 — усиление однофотонных состояний, Qμ — усиление состояний сигнала с интенсивностью μ, H2(p) — функция бинарной энтропии Шеннона, e1 — коэффициент ошибок однофотонных состояний, Eμ — среднее значение квантовых битовых ошибок за сеанс, ε — коэффициент эффективности коррекции ошибок, равный 1.44 [15].
Для того, чтобы определить длину ключа, необходимо знать, от чего зависят и как рассчитываются величины, входящие в состав выражения (15). Например, под бинарной энтропией Шеннона понимают меру неопределенности случайной величины, принимающей значения 0 и 1. Она показывает количество информации, содержащейся в сообщении, и определяется выражением
, (16)
где p — вероятность того, что случайная бинарная величина равна 1.
Если значение p близко к 0 или 1, то энтропия Шеннона будет близка к 0. Это означает, что сообщение содержит мало информации и практически известно. Если же вероятность p равна 0.5, то энтропия будет максимальна и равна 1.
Важным параметром является вероятность обнаружения на стороне Боба i-фотонного состояния, посланного Алисой, Yi (здесь и далее i = 1). Эта величина зависит от уровней фонового шума Y0 и полезного квантового сигнала ηi:
. (17)
Усиление i-фотонного состояния Qi представляет собой произведение вероятности того, что Алиса отправит i-состояние в соответствии с распределением Пуассона и вероятности события обнаружения этого состояния Бобом:
. (18)
Частота ошибок i-фотонного состояния ei определяется выражением
, (19)
где e0 — частота ошибок фона, равная 0.5, т. к. фон случаен.
Усиление состояний сигнала Qμ с интенсивностью μ вычисляется по формуле
. (20)
Среднее значение квантовых битовых ошибок за сеанс Eμ:
. (21)
В табл. 4 представлены итоговые результаты моделирования сеанса КРК для случая пролета спутника над наземной станцией через зенит (сеанс начинается, когда спутник достигает угловой высоты, равной 20°). Значения Rsift и QBER из данной таблицы согласовываются с графиками на рис. 3 для временного промежутка от –135 с до 135 с, что соответствует длительности сеанса КРК 270 с.
Таблица 4. Результаты моделирования сеанса распределения ключей для случая пролета спутника над наземной станцией через зенит
Параметр | Погодные условия | |||
Отличные | Хорошие | Средние | Плохие | |
Коэффициент экстинкции атмосферы, α | 0.1—0.3 | 0.31—0.50 | 0.51—0.70 | 0.71—1.00 |
Потери в канале, ηlink, дБ | 31—40 | 32—43 | 33—45 | 34—49 |
Скорость генерации ключа, Rsift, кбит·с-1 | 1.9—15.0 | 1.2—12.6 | 0.7—10.4 | 0.3—7.6 |
Длина просеянного ключа, lsift, кбит | 2000 | 1582 | 1245 | 830 |
Длина секретного ключа, lsec, кбит | 456 | 310 | 192 | 36 |
QBER, % | 1.20—2.60 | 1.24—3.64 | 1.29—5.34 | 1.4—10.7 |
Описание оптического приемного модуля системы КРК
Оптическая схема наземного приемного модуля для регистрации однофотонных сигналов со спутника показана на рис. 4а. Данный модуль объединяет в себе следующие функциональные узлы:
1) телескоп Шмидта‒Кассегрена (SCT) с фокусным расстоянием 3 м и апертурой 0.3 м для наведения и слежения за космическими аппаратами;
2) модуль анализа и обработки оптических сигналов (APS), где разделяются квантовый и опорный сигналы от спутника;
3) поляризационный анализатор (PA), представляющий собой блок декодера по протоколу BB84;
4) блок электроники с детекторами одиночных фотонов (SPD) для регистрации сигнального излучения;
5) телескоп-гид (GT) для первичного грубого наведения на спутник;
6) блок приема опорного сигнала (SCM) для синхронизации времени отправки и приема квантовых состояний.
Рис. 4. Наземная оптическая станция [11]: оптическая схема (а) и внешний вид (б)
Принцип работы приемного модуля заключается в следующем. Когда спутник выходит из тени Земли, в систему автоматического наведения и сопровождения наземной станции передаются TLE-координаты космического аппарата для начала процедуры отслеживания целевого объекта. TLE (Two-Line Element set) — это формат описания орбитальных параметров и прогнозирования положений и скоростей искусственных спутников Земли в режиме реального времени. Опорно-поворотное устройство станции осуществляет поворот оптико-механического блока по заданному TLE-координатами направлению. Начинается процедура грубого слежения за спутником, в ходе которой происходит удержание целевой точки в поле зрения камеры телескопа-гида за счет лишь прецизионной механики, позволяющей добиться плавного и равномерного ведения трубы телескопа. Далее включается наземный лазерный маяк с длиной волны 671 нм, работающий в непрерывном режиме. После обнаружения этого сигнала системой наведения спутника, космический аппарат активирует свой лазер-маяк с длиной волны 532 нм и частотой следования импульсов 10.7 кГц. Это излучение также используется приемным терминалом для синхронизации отправленных и принятых квантовых состояний. Передача квантового сигнала на длине волны 850 нм начинается не ранее, чем спутник поднимется на 20° над горизонтом. Собранное телескопом Шмидта‒Кассегрена (SCT) излучение разделяется по спектру дихроичным зеркалом (DM). Сигнал маяка проходит сквозь данное зеркало без отражения и попадает на матрицу камеры (CAM1), обеспечивающей обратную связь для контура точного наведения. В данном случае специальный алгоритм вычисляет положение центра лазерного пятна относительно центра матрицы и на основе полученных значений отклонений корректирует скорость поворота оптико-механического блока. Квантовый же сигнал, отражаясь от дихроичного зеркала, проходит через набор моторизованных волновых пластин (CP(m)), предназначенных для компенсации поворота плоскости поляризации фотонов, узкополосный интерференционный фильтр (IF1) и попадает в блок дешифратора по протоколу BB84 (PA). Там фотоны случайным образом проходят в одно из плеч светоделительного куба (BS), который делит падающий на него поток излучения в соотношении 50:50, реализуя тем самым пассивный выбор измерительного базиса. В каждом из базисов расположен поляризационный светоделительный куб (PBS), который разделяет ортогональные поляризационные состояния на два канала. Для увеличения коэффициента экстинкции поляризации и, соответственно, уменьшения уровня квантовых битовых ошибок, в отражающем плече такого куба дополнительно установлен поляризатор (POL). Таким образом, попадая в один из четырех каналов, сигнал при помощи коллиматора заводится в оптическое волокно и по нему направляется на приемную площадку детектора одиночных фотонов (SPD). При регистрации фотона детектор посылает электрический импульс на время-цифровой преобразователь (TDC), который фиксирует время срабатывания детектора.
ЗАКЛЮЧЕНИЕ
Описана мобильная наземная оптическая станция для приема квантовых сигналов со спутника, разработанная с учетом требований к миниатюризации. Вес станции составляет 150 кг, что позволяет разместить ее и протестировать на крышах зданий в черте мегаполиса.
Осуществлено математическое моделирование процедуры спутникового квантового распределения ключей, выполнена оценка длины секретного ключа и уровня ошибок в ходе сеанса квантовой связи между спутником Micius и нашим наземным приемным терминалом. Полученные результаты показывают перспективу использования оптических станций малого формата для решения задач спутниковой квантовой криптографии. Применение таких портативных станций позволит в будущем развернуть глобальную квантовую сеть. Кроме того, наш оптический модуль может быть использован в качестве приемного терминала для высокоскоростной космической лазерной связи.
Работа выполнена при поддержке Министерства образования и науки Российской Федерации в рамках Программы стратегического академического лидерства «Приоритет 2030» (Стратегический проект «Квантовый Интернет»).
About the authors
K. A. Barbyshev
QSpace Technologies LLC; Moscow Power Engineering Institute
Author for correspondence.
Email: k.barbyshev@goqrate.com
Russian Federation, Moscow; Moscow
A. V. Duplinsky
QSpace Technologies LLC; HSE University
Email: k.barbyshev@goqrate.com
Russian Federation, Moscow; Moscow
A. V. Khmelev
QSpace Technologies LLC
Email: k.barbyshev@goqrate.com
Russian Federation, Moscow
V. L. Kurochkin
QSpace Technologies LLC; MISIS National University of Science and Technology
Email: k.barbyshev@goqrate.com
Russian Federation, Moscow; Moscow
References
- Трушечкин А.С., Киктенко Е.О., Кронберг Д.А., Федоров А.К. // УФН. 2021. Т. 191. № 1. С. 93; Trushechkin A.S., Kiktenko E.O., Kronberg D.A., Fedorov A.K. // Phys. Usp. 2021. V. 64. No. 1. P. 88.
- Балыгин К.А., Кулик С.П., Молотков С.Н. // Письма в ЖЭТФ. 2022. Т. 116. № 2. С. 128; Balygin K.A., Kulik S.P., Molotkov S.N. // JETP Lett. 2022. V. 116. No. 2. P. 128.
- Курочкин В.Л., Неизвестный И.Г. // Изв. РАН. Сер. физ. 2015. Т. 79. № 2. С. 195; Kurochkin V.L., Neizvestnyj I.G. // Bull. Russ. Acad. Sci. Phys. 2015. V. 79. No. 2. P. 173.
- Wang S., Yin Z.Q., He D.Y. et al. // Nature Photonics. 2022. V. 16. No. 2. P. 154.
- Хмелев А.В., Дуплинский А.В., Майборода В.Ф. и др. // Письма в ЖТФ. 2021. Т. 47. № 17. С. 46; Khmelev A.V., Duplinsky A.V., Mayboroda V.F. et al. // Tech. Phys. Lett. 2021. V. 47. No. 12. P. 858.
- Liao S.K., Cai W.Q., Liu W.Y. et al. // Nature. 2017. V. 549. No. 7670. P. 43.
- Yin J., Li Y.H., Liao S.K. et al. // Nature. 2020. V. 582. No. 7813. P. 501.
- Liao S.K., Cai W.Q., Handsteiner J. et al. // Phys. Rev. Lett. 2018. V. 120. No. 3. Art. No. 030501.
- Liao S.K., Yong H.L., Liu C. et al. // Nature Photonics. 2017. V. 11. No. 8. P. 509.
- Bedington R., Arrazola J.M., Ling A. // NPJ Quantum Inf. 2017. V. 3. Art. No. 30.
- Khmelev A.V., Ivchenko E.I., Miller A.V. et al. // Entropy. 2023. V. 25. No. 4. Art. No. 670.
- Kim I.I., McArthur B., Korevaar E.J. // Proc. SPIE. 2001. V. 4214. P. 26.
- Young A.T., Irvine W.M. // Astron. J. 1967. V. 72. No. 8. P. 945.
- Ma X., Qi B., Zhao Y., Lo H.K. // Phys. Rev. A. 2005. V. 72. No. 1. Art. No. 012326.
- Kiktenko E.O., Trushechkin A.S., Lim C.C.W. et al. // Phys. Rev. Appl. 2017. V. 8. No. 4. Art. No. 044017.
Supplementary files
