НАУЧНО-ТЕХНИЧЕСКИИ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ май-июнь 2016 Том 16 № 3 ISSN 2226-1494 http://ntv.ifmo.ru/
SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS
May-June 2016
Vol. 16 No 3 ISSN 2226-1494
http://ntv.ifmo.ru/en
УДК 51-72, 538.915, 538.911 КВАНТОВО-МЕХАНИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОСТРАНСТВЕННОЙ И ЗОННОЙ СТРУКТУРЫ СЦИН ТИЛЛЯЦИОННОГО КРИСТАЛЛА Y3AI5O12
И.И. Врубель3, Р.Г. Полозков3, И.А. Шелыхь
a Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация ь Исландский университет, Институт Науки, Рейкьявик, IS-107, Исландия Адрес для переписки: [email protected] Информация о статье
Поступила в редакцию 01.04.16, принята к печати 18.04.16 doi: 10.17586/2226-1494-2016-16-3-409-415 Язык статьи - русский
Ссылка для цитирования: Врубель И.И., Полозков Р.Г., Шелых И.А. Квантово-механическое моделирование пространственной и зонной структуры сцинтилляционного кристалла Y3Al5Oi2 // Научно-технический вестник информационных технологий, механики и оптики. 2016. Т. 16. № 3. С. 409-415. doi: 10.17586/2226-1494-2016-16-3-409-415
Аннотация
Исследована электронная и пространственная структура элементарной ячейки иттрий-алюминиевого граната. Представлена квантово-механическая модель. Для расчета геометрии элементарной ячейки использованы полуэмпирические методы с параметризацией PM6 и PM7. Расчет зонной структуры выполнен в рамках теории функционала плотности с использованием обменно-корреляционного PBE-функционала. Построены гистограммы длин связей металл-кислород для оптимизированной геометрии элементарной ячейки кристалла. Проведен сравнительный анализ использованных методов, сформирована рекомендация об их применимости к решению задач подобного класса. Произведен расчет одноэлектронных волновых функций и энергии соответствующих состояний. Получена оценка ширины запрещенной зоны. Построены графики конфигурации валентной зоны. Показано, что предложенный алгоритм расчета корректно описывает пространственную и зонную структуру сцинтилляционного кристалла Y3Al5O12. Результаты работы могут быть использованы для улучшения характеристик сцинтилляционных кристаллов на основе гранатов. Ключевые слова
иттрий-алюминиевый гранат, сцинтилляционный кристалл, оптимизация геометрии, теория функционала плотности,
зонная структура
Благодарности
Работа выполнена при государственной финансовой поддержке ведущих университетов Российской Федерации (субсидия 074-U01). Авторы также выражают благодарность финансовой поддержке Министерства образования и науки РФ (Соглашение № 14.587.21.0020, уникальный идентификатор проекта RFMEFI58715X0020). И.А.Ш. выражает благодарность программе «ITMO Visiting Professorship».
QUANTUM-MECHANICAL MODELING OF SPATIAL AND BAND STRUCTURE OF Y3AL5O12 SCINTILLATION CRYSTAL I.I. Vrubel 3, R.G. Polozkov 3, I.A. Shelykh ь
3 ITMO University, Saint Petersburg, 197101, Russian Federation b University of Iceland, Science Institute, Reykjavik, IS-107, Iceland Corresponding author: [email protected] Article info
Received 01.04.16, accepted 18.04.16 doi: 10.17586/2226-1494-2016-16-3-409-415 Article in Russian
For cit3tion: Vrubel I.I., Polozkov R.G., Shelykh I.A. Quantum-mechanical modeling of spatial and band structure of Y3Al5O12 scintillation crystal. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2016, vol. 16, no. 3, pp. 409-415. doi: 10.17586/2226-1494-2016-16-3-409-415
Abstract
Spatial and electronic structures of a unit cell of yttrium-aluminum garnet have been studied. Quantum-mechanical model have been presented. Semi-empirical methods PM6 and PM7 have been used for geometry optimization of the crystal unit cell. Band structure has been calculated within density functional theory with the use of PBE exchange-correlation functional. Histograms of metal-oxygen distances for equilibrium geometry have been constructed. Comparison of the used methods has been carried out and recommendation about their applicability for such problems was given. The single-particle wave functions and energies have been calculated. The bandgap was estimated. The band structure was plotted. It was shown that
the method gives reliable results for spatial and band structure of Y3Al5012 scintillation crystal. The results of this work can
be used for improvement of characteristics of garnet scintillation crystals.
Keywords
yttrium-aluminum garnet, scintillation crystal, geometry optimization, density functional theory, band structure Acknowledgements
This work was financially supported by the Government of the Russian Federation (Grant 074-U01).The authors acknowledge support from Russia's Federal Program "Scientific and Educational Manpower for Innovative Russia" (grant No RFMEFI58715X0020). I.A.Sh. thanks "ITMO Visiting Professorship" scheme.
Введение
Гранаты - широко известный класс кристаллов, используемый в качестве люминофоров, твердотельных лазеров или сцинтилляторов [1]. Сцинтиллятор - это объект, преобразующий энергию ионизирующего излучения в фотоны оптического диапазона, которые могут быть легко зарегистрированы современными приемниками и обработаны электронной аппаратурой. Именно сцинтилляционные материалы в настоящее время являются ограничивающим элементом детектирующих систем, построенных на этом принципе. В случае регистрации гамма-излучения к сцинтилляционному кристаллу предъявляется множество требований [2] : большая плотность и эффективный атомный номер вещества, высокая конверсионная эффективность, максимальная скорость люминесценции, отсутствие послесвечения и другие. Улучшение любого из представленных параметров или поиск оптимального соотношения между ними позволяют добиться улучшения технических параметров приборов, таких как динамический диапазон, быстродействие или дозовая нагрузка ионизирующего излучения. Характеристики материалов, соответствующие предъявляемым требованиям, прямо или косвенно определяются химическим составом и электронной структурой материала, а значит, исследование этих параметров с помощью квантово-механического моделирования представляет собой одну из важных компонент многопараметрической задачи оптимизации сцинтилляторов.
Простейшей структурой элементарной ячейки граната является многогранник, в центре которого располагается атом металла, а в вершинах - атомы кислорода. Изучаемый класс гранатов содержит в ячейке комбинацию из 24 додекаэдров, 16 октаэдров и 24 тетраэдров [3], что составляет химическую формулу D0D3(0CT2TET3)012. При этом на позициях в центрах додекаэдров «D0D» могут находиться атомы редкоземельных элементов Lu, У и Gd, а на позициях в центрах октаэдров и тетраэдров «0CT» и «TET» - Al или Ga соответственно [4-6]. В зависимости от атома металла, образующего многогранник, могут получаться фигуры с различной длиной связи и количеством вершин. При согласовании объемов и количества вершин многогранников атомы кислорода способны обобществляться между атомами металла и образовывать элементарную ячейку кристалла. Целью настоящей работы является исследование пространственной и электронной структуры классического материала указанного семейства - иттрий-алюминиевого граната Y3Al5012 (YAG) с помощью методов квантово-механического моделирования.
Используемые методы и приближения
Как известно, задача о расчете пространственной и электронной структуры любой многоатомной системы сводится к численному решению молекулярного уравнения Шредингера. Существуют два подхода к решению этой задачи: первый составляют методы расчета из первых принципов (ab-initio), второй - полуэмпирические. Первые обладают высокой точностью, но, зачастую, очень большим временем расчетов, вторые проигрывают в точности, но существенно выигрывают в скорости вычислений. При этом принципиально с точки зрения подхода к решению уравнения Шредингера эти методы не отличаются. Разница состоит лишь в том, что в процессе расчета вводятся приближения, позволяющие сократить общее количество двухэлектронных интегралов, на вычисление которых расходуется основная часть машинного времени; кроме того, некоторые двухэлектронные и одноэлектронные интегралы, а также интегралы перекрытия не рассчитываются, а оцениваются, опираясь на экспериментальные данные.
Так как процесс оптимизации геометрии предполагает поиск на каждом шаге самосогласованного решения для текущих координат атомов электронной структуры, то использование подхода ab-initio становится крайне ресурсозатратным в смысле машинного времени. Исходя из этого, для оптимизации пространственной структуры элементарной ячейки YAG в настоящей работе предлагается использовать полуэмпирические методы расчета, а для вычисления электронной структуры в рамках уже оптимизированной равновесной геометрии - подход ab-initio.
Семейство полуэмпирических методов (AM1, PM3, PM6, PM7) квантово-механического моделирования [7] основывается на расчете молекулярных орбиталей соединений с помощью параметризованных наборов, описывающих различные межатомные взаимодействия. Исторически это направление развивается по пути увеличения числа параметров, их точности и расширению набора опорных (референсных) химических соединений, используемых для вычисления величин параметров. Среди достоинств метода стоит отметить высокую вычислительную эффективность и приемлемую достоверность результатов [8] для «характерных» для метода соединений, а именно, схожих с опорными соединениями, использован-
ными для параметризации. Однако следует понимать, что при использовании полуэмпирических методов для моделирования «нехарактерных» соединений вычисления могут давать сомнительные результаты.
Наиболее известной реализацией подхода ab-initio является теория функционала плотности. Основу данной теории составляют две теоремы. Первая теорема была сформулирована в 1964 г. Хоэнбергом и Коном [9], и ее суть заключается в том, что все свойства электронной структуры в основном состоянии определяются электронной плотностью п(г). Причем полная энергия системы, записанная как функционал этой плотности, имеет единственный минимум, равный энергии основного состояния:
5£[п(г)] = 5(Г[п(г)] + Uext[n(r)] + №«] + Ех<Мг)]) = 0, где Г[п(г)] - вклад кинетической энергии; Wext[n(r)] - потенциальная энергия во внешнем поле; UH[n(r)] - энергия межэлектронного кулоновского отталкивания; £,Хс[п(г)] - обменно-корреляционный функционал, вид которого в общем случае неизвестен.
Вторая - теорема Кона-Шэма [10] - имеет практическое значение и позволяет вывести с помощью вариационного принципа, примененного к полной энергии системы, одночастичные уравнения: 1
2А + 7eff фг(г) = егф.(г), Г n(r')dr'
VeffGO = Vext(r) + J | + ^xc(r);
7xc(r) = 5£-[n(r)]
5n(r)
|2
n(r) = 2Г=1|фг(г)Г
где ф¿ - волновые функции молекулярных орбиталей; - соответствующие им энергии; Vext(r) - потенциал внешнего поля; V^c(r) - обменно-корреляционный потенциал. Очевидно, что полная электронная плотность п(г) определяется видом молекулярных орбиталей ф^г).
Общий вид обменно-корреляционного функционала неизвестен в принципе, поэтому существует множество его аппроксимаций. Самой простой и исторически первой аппроксимацией являлось приближение локальной плотности (Local Density Approximation) [10], в рамках которого обменно-корреляционный функционал определяется из модели однородного электронного газа. Как следствие, разрабатываемые аппроксимации обменно-корреляционной энергии были полностью локальными, т.е. величина функционала полностью определялась величиной электронной плотности в заданной точке. Позже был разработан класс аппроксимаций GGA (generalized gradient approximation) [11], которые учитывали в функциональной зависимости также градиентную компоненту электронной плотности. В наших расчетах мы использовали именно GGA-аппроксимацию, сформулированную группой ученых J. Perdew, K. Burke и M. Ernzerhof, имеющую название PBE [11]. Преимуществами этого приближения являются отсутствие параметризации и хорошие результаты вычислений электронной структуры для широкого класса веществ.
Расчет равновесной пространственной структуры ячейки YAG
Для расчета равновесной геометрии всего класса композиций на базе гранатов была сформирована стартовая геометрия (рис. 1), в которой присутствуют 64 позиции для атомов металлов и 96 - для атомов кислорода, координаты которых масштабируются в соответствии с постоянной решетки структуры. Позиции всех атомов были подвергнуты возмущению (не более 10% от длины связи металл-кислород), что позволило получить больше информации о корректности работы используемых методов и создать стартовую напряженность в ячейке, требующую релаксации. Оптимизация геометрии искаженной ячейки производилась с помощью двух полуэмпирических методов PM6 и PM7 с использованием программного пакета MOPAC2012 [12, 13] с целью сравнения и формирования вывода об их применимости. Так как оба метода являются параметрическими, то априори неизвестно, какой из них даст лучший результат для конкретной системы. Из общих соображений можно лишь сказать, что метод PM7 более адаптирован к моделированию кристаллических структур, чем PM6, за счет включения в процедуру расчета большего количества параметров при учете парных взаимодействий. При моделировании были использованы некоторые общие ограничения на позиции атомов. Первое ограничение связано с размерами ячейки, а именно, шаг решетки был фиксирован при расчетах, а его величина взята из экспериментальных данных [3]. Второе ограничение заключается в том, что атомы, находящиеся на плоскостях ячейки, образованных векторами трансляции и началом координат, могут перемещаться только внутри них.
Для проверки результатов оптимизации геометрии были использованы экспериментальные данные о длине связи металл-кислород [14]. Для алюминиевых тетраэдрических позиций длина связи равна 0,178 нм, а для октаэдрических позиций - 0,191 нм. Додекаэдр, образованный атомами иттрия, является геометрически искаженным и имеет два характерных размера: 0,230 нм и 0,243 нм. Для полноценного сравнения результатов расчета с экспериментальными данными к элементарной ячейке было добавлено 26 копий, смещенных на различные комбинации векторов трансляции, в результате чего была образована
«макроячейка». В этой структуре для металлов, принадлежавших к исходной ячейке, были построены гистограммы распределений длин связей до ближайших атомов кислорода. Результаты для стартовой и оптимизированных с помощью РМ6 и РМ7 методов структур приведены рис. 2 и 3.
1,4
Х, мм
Рис.
1. ХТ-проекция структуры граната Y3Al5O12, где позиции «У» обозначены синими кругами, позиции «А1» - зелеными, а атомы кислорода - красными
100
н
3
<й «
о а о
4
о
5 И и о S
о g
о и н о
о
80
60
40
20
(1)
эксперимент
ш
ъ
I
I
(2)
0,165 0,170 0,175 0,180 0,185 0,190 0,195 0,200 0,205 0,210 Длина связи металл-кислород, нм
стартовая геометрия оптимизация РМ6 ^Я оптимизация РМ7
Рис. 2. Гистограмма длин связи Al-O в макроячейке YAG.
(1) - тетраэдрические позиции; (2) - октаэдрические позиции
На рис. 2 каждая структура имеет два распределения на гистограмме, отвечающие тетраэдриче-ским (1) и октаэдрическим (2) позициям атомов Al. Соответствующие экспериментальные значения [14] обозначены штрихованными областями. Число атомов кислорода в обоих пиках равно количеству вершин соответствующих образующих многогранников: 24^4=96 для тетраэдрических позиций и 16^6=96 для октаэдрических позиций. Необходимо отметить, что в результате оптимизации геометрии характер распределений изменяется. До оптимизации распределение имеет большую дисперсию и не может быть с высокой точностью характеризовано одним определенным значением длины связи. В результате оптимизации геометрии распределения принимают форму монолиний с точностью до дискрета гистограммиро-вания. Из рис. 2 видно, что результаты, полученные с помощью метода PM6, несколько лучше согласуются с экспериментальными данными, чем с помощью метода PM7.
0
120
н
В
CÖ
п о л о
4
о
к и и о
5 о
S3
о и н о
CD
100
80
60
40
о
20
0,224
0,228 0,232
0,236
0,240 0,244
0,248
Длина связи металл-кислород, нм
стартовая геометрия оптимизация РМ6 оптимизация РМ7
Рис. 3. Гистограмма длин связи Y-O в макроячейке YAG
Распределение на рис. 3 соответствует позициям атомов Y. Интеграл под всеми тремя распределениями равен 192, что согласуется с суммарным количеством вершин в иттриевых позициях (24^8=192). Видно, что после оптимизации с использованием алгоритма PM7 распределение длин связей Y-O разбивается на две группы, центры которых хорошо согласуются с экспериментальными значениями. При этом распределение, полученное методом PM6, не имеет подобного расщепления. Такое отсутствие расщепления длин связей металл-кислород для додекаэдрических позиций свидетельствует о наличии существенных недостатков в параметризации PM6.
Расчет электронной структуры ячейки YAG
Расчет электронной структуры элементарной ячейки граната был произведен в рамках теории функционала плотности в базисе плоских волн с использованием PBE-аппроксимации обменно-корреляционного функционала для каждого типа геометрии. Вычисления производились с помощью программ «pw» (planewave) и «dos» (density of states) из пакета квантово-механического моделирования QuantumEspresso [15], с использованием псевдопотенциалов атомов с конфигурациями Y(5s24d1), Al(3s23p1) и O(2s22p4). Полученные в результате расчетов числа заполнения молекулярных орбиталей и соответствующие им значения энергий использовались для построения распределения электронных состояний по энергии для всех типов геометрий.
к к о
* 5
CD Д
* £
(Ъ о
Й S
н о о
400
300
200
100 0
-20
400
300
200 100 0
-20
-10
0
10
400
300
200
100 0
-20
Энергия, отсчитываемая от энергии Ферми, эВ а б в
Рис. 4. Плотность энергетических состояний электронов в элементарной ячейке YAG, рассчитанная для трех типов геометрий: стартовая геометрия (а); геометрия, оптимизированная методом PM6 (б);
геометрия, оптимизированная методом PM7 (в)
На рис. 4 представлена плотность электронных состояний в зависимости от энергии для трех типов геометрии. Во всех случаях энергия Ферми систем приведена к нулю. Диапазону отрицательных
0
энергий соответствуют заполненные молекулярные орбитали, а положительному - незаполненные. Таким образом, на рис. 4 ширина запрещенной зоны соответствует левой границе состояний в области положительных энергий. В результате анализа рассчитанных плотностей электронных состояний были получены значения ширины запрещенной зоны в трех случаях: без оптимизации, с оптимизацией в рамках PM6 и PM7. Полученные результаты сравнивались с известным экспериментальным значением ширины запрещенной зоны для YAG, равным по разным данным [1, 3] порядка 7 эВ. Наиболее точный результат, 5,5 эВ, достигнут с использованием пространственной структуры, полученной в рамках метода PM7. Оптимизация геометрии в рамках метода PM6 и стартовая искаженная геометрия дают примерно равную и менее точную величину - 5 эВ. Стоит отметить, что при расчете с использованием теории функционала плотности с учетом только локальных корреляций, даже в рамках достаточно точного приближения PBE, значение энергии нижнего незаполненного состояния определяется с систематической ошибкой, обусловленной ограничением теории. Как следствие, в расчетных данных возникает систематическая ошибка ширины запрещенной зоны порядка 20%, т.е. в рассматриваемом случае - 7*0,2=1,4 эВ. А значит, ожидаемое расчетное значение должно быть равным примерно 5,6 эВ, что позволяет сделать вывод об адекватном описании ширины запрещенной зоны материала с помощью комбинации методов PM7 и PBE.
Заключение
Пространственная структура и плотность электронных состояний элементарной ячейки YAG были исследованы с помощью комплекса методов квантово-механического моделирования. Для определения равновесной геометрии был использован полуэмпирический метод расчета с параметрическими наборами PM6 и PM7. В результате оптимизации геометрии элементарной ячейки YAG установлено, что параметризация PM7 является более предпочтительной. Квантово-механическое моделирование электронной структуры из первых принципов позволили визуализировать электронный спектр элементарной ячейки и оценить ширину запрещенной зоны материала. Расчет показал актуальность выбора корректного алгоритма оптимизации геометрии, так как ошибка в 0,01 нм (по сравнению с длиной связи 0,2 нм), допущенная при оптимизации геометрии, приводит к значительным искажениям (0,5 эВ при рассчитанной величине запрещенной зоны 5,5 эВ) в зонной структуре элементарной ячейки.
В результате проведенного исследования был сформирован алгоритм квантово-механического моделирования, позволяющий в дальнейшем проводить более сложные многопараметрические исследования многокомпонентных гранатов. Представлен способ контроля результатов оптимизации геометрии, основанный на построении гистограммы длин связи металл-кислород всех позиций элементарной ячейки. Показана актуальность достоверности оптимизации пространственной структуры на примере ее влияния на зонную структуру.
Литература
1. Chen L., Chen X., Liu F. et al. Charge deformation and orbital hybridization: intrinsic mechanisms on tunable chromaticity of Y3Al5O12:Ce3+ luminescence by doping Gd3+ for warm white LEDs // Scientific Reports. 2015. V. 5. Art. 11514. doi:10.1038/srep11514
2. Ronda C., Wieczorek H., Khanin V., Rodnyi P. Review-scintillators for medical imaging: a tutorial overview // ECS Journal of Solid State Science and Technology. 2016. V. 5. N 1. P. 3121-3125. doi:10.1149/2.0131601jss
3. Huang Z., Zhang L., Pan W. Antisite defect in nonstoichiometric yttrium aluminum garnet: Experimental and first-principles calculation // Journal of the European Ceramic Society. 2014. V. 34. N 3. P. 783-790. doi:10.1016/j.jeurceramsoc.2013.09.007
4. Dorenbos P. Electronic structure and optical properties of the lanthanide activated RE3(Al1-xGax)5O12 (RE=Gd, Y, Lu) garnet compounds // Journal of Luminescence. 2013. V. 134. P. 310-318. doi: 10.1016/j.jlumin.2012.08.028
5. Fasoli M., Vedda A., Nikl M., Jiang C., Uberuaga B.P., Andersson D.A., McClellan K.J., Stanek C.R. Bandgap engineering for removing shallow traps in rare-earth Lu3Al5O12 garnet scintillators using Ga3+ doping // Physical Review B. 2011. V. 84. Art. 081102(R). doi: 10.1103/PhysRevB.84.081102
6. Munoz-Garcia A.B., Seijo L. Structural, electronic, and spectroscopic effects of Ga codoping on Ce-doped yttrium aluminum garnet: first-principles study // Physical Review B. 2010. V. 82. Art. 184118. doi: 10.1103/PhysRevB.82.184118
7. Stewart J.J.P. Optimization of parameters for semiempirical methods VI: more modifications to the NDDO approximations and re-optimization of parameters // Journal of Molecular Modeling. 2013. V. 19. N 1. P. 132. doi: 10.1007/s00894-012-1667-x
8. Puzyn T., Suzuki N., Haranczyk M., Rak J. Calculation of quantum-mechanical descriptors for QSPR at the DFT level: is it necessary? // Journal of Chemical Information and Modeling. 2008. V. 48. N 6. P. 11741180. doi: 10.1021/ci800021p
9. Hohenberg P., Kohn W. Inhomogeneous electron gas // Physical Review. 1964. V. 136. N 3B. P. B864-B871. doi: 10.1103/PhysRev.136.B864
10.Kohn W., Sham L.J. Self-consistent equations including exchange and correlation effects // Journal of Physical Review. 1965. V. 140. N 4A. P. 1133. doi: 10.1103/PhysRev.140.A1133
11.Perdew J.P., Burke K., Ernzerhof M. Generalized gradient approximation made simple // Physical Review Letters. 1996. V.77. N 18. P. 3865-3868. doi: 10.1103/PhysRevLett.77.3865
12.Stewart J.J.P. General description of M0PAC2012 [Электронный ресурс]. Режим доступа: http:// openmopac.net/manual/index.html, свободный. Яз. англ. (дата обращения 21.04.2016).
13.Maia J.D.C., Carvalho G.A.U., Mangueira Jr. C.P. et al. GPU linear algebra libraries and GPGPU programming for accelerating MOPAC semiempirical quantum chemistry calculations // Journal of Chemical Theory and Computation. 2012. V. 8. N 9. P. 3072-3081. doi: 10.1021/ct3004645
14.Furman J.D., Gundiah G., Page K., Pizarro N., Cheetham A.K. Local structure and time-resolved photoluminescence of emulsion prepared YAG nanoparticles // Chemical Physics Letters. 2008. V. 465. N 1-3. P. 6772. doi:10.1016/j.cplett.2008.09.045
15. Giannozzi P., Baroni S., Bonini N. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials // Journal of Physics: Condensed Matter. 2009. V. 21. N 39. Art. 395502. doi: 10.1088/0953-8984/21/39/395502
Врубель Иван Игоревич Полозков Роман Григорьевич
Шелых Иван Андреевич
аспирант, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
кандидат физико-математических наук, ведущий научный сотрудник, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected] кандидат физико-математических наук, профессор, Исландский университет, Институт Науки, Рейкьявик, 1Б-107, Исландия, [email protected]
Ivan I. Vrubel Roman G. Polozkov Ivan A. Shelykh
postgraduate, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
PhD, leading researcher, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
PhD, Professor, University of Iceland, Science Institute, Reykjavik, IS-107, Iceland, [email protected]