Собственные волны плоского симметричного анизотропного волновода Текст научной статьи по специальности «Физика»

#

Вестник РУДН. Серия МИФ

ИиБК Лоигпа1 of М1РЬ

http://journals.rudn.ru/miph

2018 У01. 26 N0. 2 119-128

Математическое моделирование

УДК 535.012.2

БГО: 10.22363/2312-9735-2018-26-2-119-128 Собственные волны плоского симметричного анизотропного

волновода

О. Н. Бикеев*, К. П. Ловецкий^, А. Л. Севастьянов^"

* Институт физических исследований и технологий Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198 ^ Кафедра прикладной информатики и теории вероятностей Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198

Получены точные дисперсионные уравнения для плоского симметричного диэлектрического анизотропного волновода, в котором изотропным является волноведущий слой, а обрамляющие среды предполагаются анизотропными одноосными средами. Тензоры диэлектрической проницаемости обрамляющих сред при этом не предполагаются диагональными, а именно — у одной из них этот тензор образован путём поворота диагонального тензора на некоторый угол между оптической осью анизотропной среды и направлением распространения электромагнитной волны. Тензор диэлектрической проницаемости другой анизотропной среды повёрнут на такой же угол, но в противоположном направлении, при этом оптические оси обеих обрамляющих сред лежат в плоскости, параллельной границам волноведущей структуры. Таким образом, в обрамляющих средах поддерживается существование шестикомпонентных электромагнитных волн.

В дисперсионных свойствах такого волновода наблюдаются некоторые особенности, по сравнению со случаем, когда обрамляющие среды предполагаются изотропными.

Обнаружено, что первая симметричная мода такого волновода имеет конечное замедление при нулевой толщине изотропного слоя, что говорит о возможности возникновения поверхностных электромагнитных волн (так называемых волн Дьяконова) на границах этого изотропного слоя. Отмечено, что переход антисимметричной моды в Дьяконовскую волну происходит при конечной толщине волноведущего слоя. Приведены зависимости величины замедления элементарной (симметричной) моды от угла поворота оптической оси анизотропных сред относительно направления распространения направляемой волны волноводной структуры.

Кроме дисперсионных свойств такого волновода исследованы поперечные распределения полей направляемой электромагнитной волны, приведены точные аналитические выражения для амплитуд полей этой волны.

Ключевые слова: поверхностная электромагнитная волна, анизотропная среда, тензор диэлектрической проницаемости, волноводные моды, поперечное распределение амплитуд, направляемые волны

Интерес исследователей к свойствам волн Дьяконова является постоянным на протяжении нескольких десятков лет. За это время изучены многие особенности таких волн, возникающих на границах анизотропных сред в различных геометриях. Наиболее важной особенностью волн Дьяконова является (в отличие от плазмонов) отсутствие потерь энергии по мере их распространения вдоль границы анизотропной среды.

Основные результаты исследования свойств таких волн достаточно подробно изложены в обзорах [1,2]. Следует отметить, что за столь долгий период исследования

Статья поступила в редакцию 9 января 2018 г.

Публикация подготовлена при поддержке Программы РУДН «5-100» и при частичной поддержке гранта РФФИ 18-07-00567.

1. Введение

волн Дьяконова экспериментальных работ проведено очень мало [3], что обусловлено, может быть, специфическими условиями существования этих волн.

Многие авторы считают, что для их наблюдения можно использовать свойства тонкоплёночных оптических волноводов, волноведущий слой которых либо является анизотропным, либо окружающие среды являются анизотропными. Показано [4,5], что при некоторых условиях удаётся выделить случай, когда перенос энергии по волноведущей структуре осуществляется только за счёт волн Дьяконова.

Настоящая работа ставит своей целью выбор такой анизотропной волноводной структуры, которая бы допускала распространение в ней волн Дьяконова и, с другой стороны, позволяла бы провести экспериментальную проверку полученных результатов.

2. Выбор волноведущей структуры

В качестве материала для волноведущего слоя возьмём изотропную среду, диэлектрическая проницаемость которой равна е8, а обрамляющие среды будем считать анизотропными с диэлектрической проницаемостью £1,2. Используемая система координат приведена на рис. 1. Предположим, что электромагнитные волны распространяются вдоль оси Z, а по оси У структура не ограничена. Все три среды являются немагнитными и непоглощающими.

Рис. 1. Система координат и исследуемая волноводная структура

Тензоры диэлектрической проницаемости обрамляющих сред определяются следующим образом. Исходный тензор в диагональном виде выглядит как:

£ =

£о 0 0

0

£о 0

0 0

£е

Оптическая ось направлена вдоль OZ. При повороте по часовой стрелке на угол ф вокруг оси ОХ компоненты тензора преобразуются в соответствии с матрицами поворота С 1(ф) и С2(ф) [6]

ад) =

1 0 0

0 СОВ ф — 8Ш ф , С2 (ф) = С1( —

0 8Ш ф СОВ ф

Результирующие тензоры диэлектрической проницаемости £1 и £2 обрамляющих сред вычисляются по формулам

£1

CieCf

£o 0 0 '

0 Е22(Ф) £2З(Ф) .0 £23(Ф) £эз(0).

£2(Ф) = £l(-

Таким образом, для среды при X > d/2 элементы тензора вычисляются по формулам

£22 (ф) = £o cos2 Ф + £e sin2 ф, £33 (ф) = £о sin2 ф + £e COS2 ф, £23(Ф) = £32 (Ф) = - Sin ф COS ф(Ее - во),

а для среды при X < -d/2 по формулам с учётом изменения знака угла поворота на обратный

£22 (ф) = £о COS2 Ф + £е sin2 ф, £33 (ф) = £о sin2 ф + £е COS2 ф, £2з(Ф) = £32(ф) = sin фCOS ф(Ее - £о).

В сформированной таким образом анизотропной структуре будем искать решение уравнений Максвелла для направляемых (волноводных) электромагнитных волн, распространяющихся вдоль оси OZ, которая является биссектрисой угла между оптическими осями верхней и нижней анизотропных сред. Пусть зависимость полей электромагнитных волн от времени и координат имеет вид

exp {-i [wt - (кхх + куу + kzz)]}, где к = (кх, ку, kz)т — волновой вектор.

Для каждой из этих волн должно выполняться волновое уравнение (в гауссовой системе), полученное из уравнений Максвелла в векторной форме:

кхН

= -koD,

кхЕ

= коН, (к • я) =0, ( к • я) =0.

Это уравнение имеет вид:

к х

кхЕ

+ к0 • D = 0, где ко = ^, D = е(ф) • Е.

Не теряя общности, можно предположить, что решением этого волнового уравнения является набор плоских волн с фазовыми фронтами параллельными оси ОУ, или д/ду = 0, а с ним и ку = 0.

3. Дисперсионные уравнения для направляемых мод

Имея в виду эти предварительные замечания, представим поля в волноведущем слое (ев) в виде суперпозиции ТЕ- и ТМ- волн, поскольку в обрамляющих анизотропных средах предполагается наличие шестикомпонентных волн. Пусть ТЕ-волна с положительным поперечным волновым числом кх имеет амплитуду А, В — амплитуда ТМ-волны с положительным кх. Для ТЕ- и ТМ-волн, имеющих отрицательный кх, обозначим амплитуды через С и Б соответственно. Тогда суммарные поля в изотропном слое можно записать следующим образом:

Еа(х) =

На(х)

- • В • eik'xX + • D • е

Кх &х

1 • Л • е

n—ikxx

- _i_ • с • е-—

^ . ^ikxx + ф ^ ikxx

,ikxx i X . Г1 . z~¡—ikx x к

ikxx \ £s ^ ^—ikxx

- i •A • pJkxX + i •c • e~ - Is • В • eikxX + ¿ • D • e~—

Kx Kx

ikxx + Q . p—ikx x

A

Здесь введены следующие обозначения: 7 — продольное волновое число, нормированное на ко, кх — поперечное волновое число волноведущего слоя, также нормированное на ко- А, В, С и Б — неизвестные амплитуды, х — поперечная координата, х — продольная координата. Множитель е-шЬ в дальнейшем будем подразумевать, не указывая в явном виде.

Для записи амплитуд напряжённостей полей в окружающих средах необходимо решить задачу на собственные значения и собственные волны выведенного в предыдущем разделе волнового уравнения, что подробно сделано в работах [7,8]. Здесь же приведём окончательные выражения.

В верхней среде (см. рис. 1) эти поля состоят из линейной комбинации обыкновенной и необыкновенной волн с различными амплитудами (Р и Е):

Et(x) =

%7

Е h e-T*(x-d) + FL- e-T°(x-d) 2 ,

' T 2

Ht(x) =

7 '

is,

—i ■

_go tgy Ep-Te(x-d)__— Fe-T°(x-d)

T2 ts v

Ee-4X-d) + Fe-T° (x-d) tsv Fp-Te(x-i) + -L- Fp-T° (x-i)

TQ2 tg v

T2 Ee-Te(x-d) + 1 Fe-T°(x-2)

• e

г'yz

T 2

T

Те ^Tg* Ee-Te(x-d) + Tot^ Fe-T°(x-d)

>tgv

o4Z

(2)

Обозначения здесь следующие: To — мнимая часть поперечного волнового числа обыкновенной волны, нормированная на к0, To = л/72 — eo; Те — мнимая часть поперечного волнового числа необыкновенной волны (нормировано на ко),

Те = \j72 (l + IT cos2 — £е, = £е — £o, Ф — угол поворота осей симметрии

верхней анизотропной среды (положительный); Е и F — неизвестные амплитуды необыкновенной и обыкновенной волн в верхней анизотропной среде.

В нижней анизотропной среде поля записываются аналогичным образом, но поперечные волновые числа должны быть отрицательными (что приведёт к экспоненциальному убыванию всех амплитуд при удалении от волноводного слоя), отрицательным должен быть и угол поворота

Ed(x) =

-г7

o t

T 2

G T2 eTe(x+d)+ H1 eTo(x+d)

№ GeT*(x+d) + _!_ HeTo(x+d У

I 2 1 tg v

Ge4 x+d) + HeTo( x+d)

0iyz

Hd(x) =

GeT*(x+d) + HeTo(x+d)

—i ■

7' To2 ^ - ' 1 tgv —,£o ■ GeT*(x+d) + ±HeTo(x+d)

Te GeT*(x+d) + To-1- HeTo(x+d)

>tg v

(3)

• e

г yyz

С и Н — амплитуды необыкновенной и обыкновенной волн в нижней анизотропной среде.

Сшивая тангенциальные компоненты полей на границах волноведущего слоя ±й/2, получим уравнение вида:

M(A,B,...,H)T = 0,

(4)

где М[е0,ее, ±ё/2, ±ф] — квадратная матрица размером 8 х 8, а (А,В,...,Н)т — вектор неизвестных амплитуд.

Требуя равенства нулю определителя матрицы М, получим дисперсионные уравнения. Поскольку исследуемая волноведущая структура симметрична, то, как и следовало ожидать, дисперсионных уравнений будет два, одно из них для симметричных мод структуры, другое — для антисимметричных. Каждое из этих уравнений является квадратным, поэтому даёт два корня. Итак, для симметричных мод:

tg ^1,2 =

м(То£s + Те£о) tVм2(To£s + Тее0)2 - 4NLT2esec

2krNen

(5)

где для краткости введены обозначения:

ТР

М = (TG2 - £о tg ^2), N = (ТеТ0 - £о tg ^2), L = (Т2 - ^£о tg ф

Для антисимметричных мод:

м(To£s + Те£о) tVМ2 (T0£s + Те£о)2 - 4NL£s£c

^ Фз,.

2krNec

, Ф = (6)

Графически зависимости замедлений этих мод от толщины изотропного слоя выглядят следующим образом:

Рис. 2. Зависимость замедления от толщины изотропного слоя: 71,2 — симметричные моды, 73,4 — антисимметричные

Уровень 7тт = -\/£е/(1 + (Де/£о) сов ) получается из условия Те = 0, при 7 < 7тт исследуемая структура переходит в режим излучения в окружающие анизотропные среды, так как Те становится мнимой величиной.

Отметим одну интересную особенность на рис. 2, а именно — простейшая симметричная мода волноводной структуры 71 при нулевой толщине изотропного слоя имеет замедление, превышающее 7т;п. Как было показано в [8-10], вдоль границы раздела анизотропных сред в этом случае может распространяться поверхностная (Дьяконов-ская [11]) электромагнитная волна, замедление которой зависит от угла поворота

Но и в исследуемой нами волноводной структуре простейшая мода 71 может трансформироваться при определённых условиях в поверхностную, Дьяконовскую, моду. Для этого необходимо лишь подобрать соответствующий угол поворота, что иллюстрирует рис. 3.

4.2

ч н

1 Ф-- 2.8 1 £□ = 4 гэ = 16 £е = 36

% ____ 4 11ШП

1 1 1

Толщина слоя, отн. ед

Рис. 3. Совокупность поверхностных (7т, 7-Оз) и волноводных мод (72, 73, 74)

исследуемой структуры

Естественно, что при вычислении 701,3 по дисперсионным уравнениям (5) и (6)

использована подстановка кх = гд/72 — е0 = гКх, ф = гФ, Ф = Кхй/2, позволяющая вычислять Дьяконовские моды. Как видно из рис. 3, в поверхностные моды переходят первая симметричная и первая антисимметричная волноводные моды. Отметим, что диапазон углов поворота, при которых существуют поверхностные моды довольно узкий, в чём можно убедиться, построив зависимости замедлений по (5) и (6) от угла поворота при фиксированной толщине изотропного слоя <1 (рис. 4).

ч) а а;

'С

со

4.5

3.5

1 (1:= 1.5 / ::: ::: ::: ::: :: /::::::::::::::::::::::

л/ёЁ УЫ/ Область :: излучателъных мод ;;

7 /: / : у.-... |. /А'/;;

1 1.5

Ф Угол поворота, рад.

Рис. 4. Зависимость замедления от угла поворота для простейших симметричной и антисимметричной мод. Переход в поверхностные моды происходит, когда замедление превышает у^ёТ. Пунктирная линия соответствует условию Те = 0

4. Собственные поля направляемых мод

Выяснив в основном дисперсионные свойства исследуемой структуры, перейдём теперь к вычислению амплитуд полей собственных электромагнитных волн для каждой из четырёх мод, описываемых вектором А, В,..., Н)т. Поскольку уже известны все собственные числа для матрицы М[е0,ее, ±й/2, ±ф\ в (4), то не составляет труда вычислить и неизвестные амплитуды из уравнения (4). В результате имеем: Для симметричных мод (5):

Л, =

T^Sg cos фЬ — Tekxeo sin tyM Т%кх tg(р

, Bs = —%£0 sin 'ф-

Te- Тп

Тп

Es = 2 sin фТо£° COS — к*Ео Sin , Fs = —2 sin — Tekx£o sin 4

k:,

T h

-L o™X

cs = As , Ds = —Bs , Gs = —Es , Hs = —Fs .

(7)

Для антисимметричных мод (6):

T0es sin ipL + TekxSo cos tyM Te — To Aas = г-^-;-;- , Bas = — £o COs ty-

T0kx tg , To£s sin ty + kxeo cos ty

Тп

Eas = 2 costy o~°"—r ' r , pas = —2costy o s

T^Ss sin ф + TekxSo cos ф (8)

kx

С as A-as , Das Das , Gas Das , Has Ec

Tokx

Lo^b X

Las x as .

Здесь, напомним, ф = кхй/2.

Теперь, воспользовавшись выражениями (1), (2) и (3), можно вычислить компоненты всех полей для любой из волноводных мод рассматриваемой структуры. В качестве примера на рис. 5 приведены графики поперечного распределения компонент напряжённости электрического поля простейшей симметричной моды при следующих параметрах волноведущей структуры: = 0,203 рад, е8 = 16, ее = 36, е0 = 4.

Рис. 5. Поперечные распределения компонент электрического поля для простейшей симметричной моды. Все величины выражены в относительных единицах, d — толщина изотропного слоя (в единицах длины волны)

Аналогичные зависимости для режима Дьяконовских мод можно получить с учётом указанной выше замены: кх = — £0 = гКх, ф = гФ. Эти зависимости

количественно будут подобны указанным на рис. 5, с той лишь разницей, что все тригонометрические функции станут гиперболическими.

5. Заключение

Таким образом, нам удалось показать, что в исследуемой волноведущей структуре возможна реализация двух типов распространения направляемых волн — обычное волноводное распространение и смешанный тип, когда наряду с обычными волно-водными присутствуют и поверхностные волны Дьяконовского типа. Переключение между указанными режимами можно обеспечить путём подбора угла между оптическими осями обрамляющих сред и направлением фазовой скорости волноводной волны. Наряду с этим получены точные аналитические выражения как для фазовых скоростей, так и для компонент напряжённостей полей направляемых волн различных типов. Отмечено, что даже при нулевой толщине изотропного слоя структура сохраняет способность поддерживать направляемую волну, которая в этом случае является поверхностной волной. Рассмотренная конструкция этой структуры, по мнению авторов, может быть реализована в эксперименте, причём в качестве изотропного слоя можно использовать жидкость с соответствующей диэлектрической проницаемостью. Последнее позволит плавно менять углы ^ оптических осей анизотропных сред относительно оси OZ.

Литература

1. Dyakonov Surface Waves: a Review / O. Takayma, L.-C. Crasovan, S. Johansen et al. // Electromagnetics. — 2008. — Vol. 28. — Pp. 126-145.

2. Polo Jr. J. A., Lakhtakia A. A Surface Electromagnetic Waves: a Review // Laser & Photonics Reviews. — 2011. — Vol. 5. — Pp. 234-246.

3. Observation of Dyakonov Surface Waves / O. Takayama, L. Crasovan, D. Artigas, L. Torner // Physical Review Letters. — 2009. — Vol. 102. — P. 043903.

4. New Waveguide Modes in Anisotropic Structures / L. Torner, C. Santos, J. P. Torres, D. Mihalache // Fiber and Integrated Optics. — 1993. — Vol. 13. — Pp. 271-280.

5. Boroujeni M. A, Shahabadi M. Modal Analysis of Multilayer Planar Lossy Anisotropic Optical Waveguides // Journal of Optics A: Pure and Applied Optics. — 2006. — Vol. 8. — Pp. 856-863.

6. Kweon G., Hwang-bo S., Kim C. Eigenmode of Anisotropic Planar Waveguide // Journal of the Optical Society of Korea. — 2004. — Vol. 8, No 3. — Pp. 137-146.

7. Lekner J. Reflection and Refraction by Uniaxial Crystals // Journal of Physics: Condensed Matter. — 1991. — Vol. 3. — Pp. 6121-6133.

8. Бикеев О. Н., Севастьянов Л. А. Поверхностные электромагнитные волны на границе анизотропных сред // Вестник РУДН. Серия: Математика. Информатика. Физика. — 2017. — Т. 25, № 2. — С. 141-148. — DOI: 10.22363/2312-9735-2017-252-141-148.

9. Nelatury S. R., Polo J. A., Lakhtakia Jr. A. Surface Waves with Simple Exponential Transverse Decay at a Biaxial Bicrystalline Interface: Errata // Journal of the Optical Society of America A. — 2007. — Vol. 24. — Pp. 2102-2102.

10. Walker D. B., Glytsis E. N., Gaylord T. K. Surface Mode at Isotropic-Uniaxial and Isotropic-Biaxial Interfaces // Journal of the Optical Society of America A. — 1998. — Vol. 15, No 1. — Pp. 248-260.

11. Аверкиев Н. С., Дьяконов М. И. Электромагнитные волны, локализованные на границе раздела прозрачных анизотропных сред // Оптика и спектроскопия. — 1990. — Т. 68, № 5. — С. 1118-1121.

UDC 535.012.2

DOI: 10.22363/2312-9735-2018-26-2-119-128

Eigen Waves of a Plane Symmetric Anisotropic Waveguide O. N. Bikeev*, K. P. Lovetskiy1", A. L. Sevastianov1"

* Institute of Physical Researches and Technologies Peoples' Friendship University of Russia (RUDN university) 6, Miklukho-Maklaya str., Moscow, 117198, Russian Federation ^ Department of Applied Probability and Informatics Peoples' Friendship University of Russia (RUDN university) 6, Miklukho-Maklaya str., Moscow, 117198, Russian Federation

Precise dispersion equations for a plane symmetric dielectric anisotropic waveguide are obtained, in which the waveguide layer is isotropic, and the framing media are assumed to be anisotropic uniaxial media. The tensors of the dielectric permittivity of the framing media are not assumed to be diagonal, namely, in one of them this tensor is formed by rotating the diagonal tensor by some angle between the optical axis of the anisotropic medium and the direction of propagation of the electromagnetic wave. The tensor of dielectric permittivity of another anisotropic medium is rotated by the same angle, but in the opposite direction, with the optical axes of both framing media lying in a plane parallel to the boundaries of the waveguiding structure. Thus, in framing media, the existence of six-component electromagnetic waves is maintained.

In the dispersion properties of such a waveguide, certain features are observed in comparison with the case when the framing media are assumed to be isotropic.

It is found that the first symmetric mode of such a waveguide has a finite deceleration with zero thickness of the isotropic layer, which indicates the possibility of the appearance of surface electromagnetic waves (the so-called Dyakonov waves) at the boundaries of this isotropic layer. It is noted that the transition of the antisymmetric mode to the Dyakonov wave occurs with a finite thickness of the waveguiding layer. Dependencies of the deceleration of the elementary (symmetric) mode on the angle of rotation of the optical axis of anisotropic media relative to the direction of propagation of the guided wave of the waveguide structure are given.

Key words and phrases: surface electromagnetic wave, anisotropic medium, permittivity tensor, waveguide modes, transverse amplitude distribution, guided waves

References

1. O. Takayma, L.-C. Crasovan, S. Johansen, D. Mihalache, D. Artigas, L. Torner, Dyakonov Surface Waves: A Review, Electromagnetics 28 (2008) 126-145.

2. J. A. Polo Jr., L. A., A Surface Electromagnetic Waves: a Review, Laser & Photonics Reviews 5 (2011) 234-246.

3. O. Takayama, L. Crasovan, D. Artigas, L. Torner, Observation of Dyakonov Surface Waves, Physical Review Letters 102 (2009) 043903.

4. L. Torner, C. Santos, J. P. Torres, D. Mihalache, New Waveguide Modes in Anisotropic Structures, Fiber and Integrated Optics 13 (1993) 271-280.

5. M. A. Boroujeni, M. Shahabadi, Modal Analysis of Multilayer Planar Lossy Anisotropic Optical Waveguides, Journal of Optics A: Pure and Applied Optics 8 (2006) 856-863.

6. G. Kweon, S. Hwang-bo, C. Kim, Eigenmode of Anisotropic Planar Waveguide, Journal of the Optical Society of Korea 8 (3) (2004) 137-146.

7. J. Lekner, Reflection and Refraction by Uniaxial Crystals, Journal of Physics: Condensed Matter 3 (1991) 6121-6133.

8. O. N. Bikeev, L. A. Sevastianov, Surface Electromagnetic Waves at the Interface of Two Anisotropic Media, RUDN Journal of Mathematics, Informational Sciences and Physics 25 (2) (2017) 141-148, in Russian. doi:10.22363/2312-9735-2017-25-2-141-148.

9. S. R. Nelatury, J. A. Polo, A. Lakhtakia Jr., Surface Waves with Simple Exponential Transverse Decay at a Biaxial Bicrystalline Interface: Errata, Journal of the Optical Society of America A 24 (2007) 2102-2102.

10. D. B. Walker, E. N. Glytsis, T. K. Gaylord, Surface Mode at Isotropic-Uniaxial and Isotropic-Biaxial Interfaces, Journal of the Optical Society of America A 15 (1) (1998) 248-260.

11. N. S. Averkiev, M. I. Dyakonov, Electromagnetic Waves Localized at the Boundary of Transparent Anisotropic Media, Optics and Spectroscopy 68 (5) (1990) 1118-1121.

© Бикеев О.Н., Ловецкий К. П., Севастьянов А. Л., 2018

Для цитирования:

Бикеев О. Н, Ловецкий К. П., Севастьянов А. Л. Собственные волны плоского симметричного анизотропного волновода // Вестник Российского университета дружбы народов. Серия: Математика. Информатика. Физика. — 2018. — Т. 26, № 2. — С. 119-128. — DOI: 10.22363/2312-9735-2018-26-2-119-128.

For citation:

Bikeev O.N., Lovetskiy K.P., Sevastianov A.L. Eigen Waves of a Plane Symmetric Anisotropic Waveguide, RUDN Journal of Mathematics, Information Sciences and Physics 26 (2) (2018) 119-128. DOI: 10.22363/2312-9735-2018-26-2-119-128. In Russian.

Сведения об авторах:

Бикеев Олег Николаевич — заведующий лабораторией Института физических исследований и технологий РУДН (e-mail: [email protected], тел.: +7 (905) 7458125) Ловецкий Константин Петрович — доцент, кандидат физико-математических наук, доцент кафедры прикладной информатики и теории вероятностей РУДН (e-mail: [email protected], тел.: +7(495)9522572)

Севастьянов Антон Леонидович — кандидат физико-математических наук, доцент кафедры прикладной информатики и теории вероятностей РУДН (e-mail: [email protected], тел.: +7(495)9522572)

Information about the authors:

Bikeev O. N. — Head of Laboratory of Institute of Physical Researches and Technologies of Peoples' Friendship University of Russia (RUDN university) (e-mail: [email protected], phone: +7 (905) 7458125)

Lovetskiy K. P. — Candidate of Physical and Mathematical Sciences, Associate Professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: [email protected], phone: +7(495)9522572)

Sevastianov A. L. — Candidate of Physical and Mathematical Sciences, assistant professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: [email protected], phone: +7(495)9522572)