Научная статья на тему 'МОДЕЛИРОВАНИЕ ОРБИТАЛЕЙ В ТЕТРАЭДРЕ КРИСТАЛЛА КУБИЧЕСКОГО УГЛЕРОДА'

МОДЕЛИРОВАНИЕ ОРБИТАЛЕЙ В ТЕТРАЭДРЕ КРИСТАЛЛА КУБИЧЕСКОГО УГЛЕРОДА Текст научной статьи по специальности «Физика»

CC BY
19
3
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГЕКСАГОНАЛЬНАЯ ПРИЗМА / ЭЛЛИПТИЧЕСКИЕ ОРБИТАЛИ СВЯЗЕЙ

Аннотация научной статьи по физике, автор научной работы — Неустроев Степан Архипович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MODELING OF ORBITALS IN C-C CRYSTAL TETRAHEDRON

The practical application of the elliptical orbital model for calculation of the potential distribution between adjacent atoms will permit to determine the energy of the particles moving to the surface from space while depositing c-C films. The visualization of c-C cell in hexagonal amorphous has been presented. The visualization has revealed the tetrahedrons as well as the orbitals location and form. The orbitals look as flat, convex, closed curves, which envelope each pair of adjacent atoms. The parameters of these curves have been determined. The calculated values of orbitals parameters of c-C tetrahedrons are consistent with the ellipse mathematical definition. The process of producing diamond films depends on the state of the gas phase particles and on the deposition surface topology, as well as on position and the form of orbitals in cubic carbon crystal tetrahedron.

Текст научной работы на тему «МОДЕЛИРОВАНИЕ ОРБИТАЛЕЙ В ТЕТРАЭДРЕ КРИСТАЛЛА КУБИЧЕСКОГО УГЛЕРОДА»

МАТЕРИАЛЫ ЭЛЕКТРОНИКИ ELECTRONICS MATERIALS

УДК 548.3: 621.382 DOI: 10.24151/1561-5405-2017-22-2-120-127

Моделирование орбиталей в тетраэдре кристалла кубического углерода

С.А. Неустроев

Национальный исследовательский университет «МИЭТ», г. Москва, Россия

doloto@list.ru

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

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

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

Ключевые слова: кубический углерод; гексагональная призма; тетраэдры; эллиптические орбитали связей; алмазные пленки.

Для цитирования: Неустроев С.А. Моделирование орбиталей в тетраэдре кристалла кубического углерода // Изв. вузов. Электроника. - 2017. - Т.22. -№2. - С. 120-127. БО!: 10.214151/1561-5405-2017-22-2-120-127

© С.А. Неустроев, 2017

Modeling of Orbitals in c-C Crystal Tetrahedron

S.A.Neoustroev

National Research University of Electronic Technology, Moscow, Russia doloto@list.ru

The practical application of the elliptical orbital model for calculation of the potential distribution between adjacent atoms will permit to determine the energy of the particles moving to the surface from space while depositing c-C films.

The visualization of c-C cell in hexagonal amorphous has been presented. The visualization has revealed the tetrahedrons as well as the orbitals location and form. The orbitals look as flat, convex, closed curves, which envelope each pair of adjacent atoms. The parameters of these curves have been determined. The calculated values of orbitals parameters of c-C tetrahedrons are consistent with the ellipse mathematical definition.

The process of producing diamond films depends on the state of the gas phase particles and on the deposition surface topology, as well as on position and the form of orbitals in cubic carbon crystal tetrahedron.

Keywords: c-C; hexahedron prism; tetrahedrons; ellipse orbital; bonds.

For citation: Neoustroev S.A. Modeling of Orbitals in c-C Crystal Tetrahedron // Proc. of universities. Electronics. - 2017. - Vol.22. - №2. - P. 120-127. DOI: 10.214151/1561-5405-2017-22-2-120-127

Введение. Получение алмазных пленок и покрытий, формируемых из газовой фазы, - перспективное направление электроники [1]. Применение таких пленок улучшает конструктивные решения и характеристики выпускаемых и создаваемых приборов, а также обеспечивает многократное удешевление производства и уменьшение количества вредных отходов. Процесс получения пленок алмаза зависит от энергетического состояния частиц газовой фазы и топологии поверхности осаждения. Энергия атомов углерода, вступающих в химическую связь с поверхностными атомами кристалла, должна содержать и кинетическую составляющую. В работе [2] показано, что значение энергии движения атомов, приближающихся к месту вхождения в кристалл, должно быть достаточным для преодоления потенциального барьера над ним. В местах с большей энергией такие атомы отталкиваются. Можно говорить, что кинетическая энергия таких частиц ограничивается нижним и верхним пределом.

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

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

Ячейка кубического углерода. Ячейка кубического углерода с-С - разновидность алмаза - относится к примитивной сингонии с плотнейшей упаковкой, как и некоторые полупроводники и металлы. Это куб, в углах и середине граней которого находятся

атомы углерода, его пространственная постоянная а = 2,826195716-10-10 м [3] и определена по плотности. В гексагональной анаморфозе такой куб выглядит трехслойной структурой ...АВСА..., ячейку которой можно представить призмой с двумя базовыми плоскостями А; и Аг+1 в виде шестиугольников, а также с двумя промежуточными плоскостями В; и С;. Атомы углерода располагаются в центре и углах гексагонов базовых плоскостей, атомы промежуточных плоскостей занимают лунки, образованные сферами атомов базовых плоскостей. Лунки и атомы промежуточных плоскостей смещены относительно расположения базовых плоскостей, но их гексагональная конфигурация расположения на плоскостях В; и подобна расположению атомов на базовых плоскостях. В работе [2] подробно приведено построение призмы.

Пространственные постоянные призмы аь и еь определяются по значению постоянной куба: сторона гексагона аь= а /Л =1,998422156-10 м [3], она также равна радиусу его описанной окружности; высота призмы еь определяется исходя из объема куба Укуб с учетом количества атомов в ячейках (в кубе их четыре, в призме - девять). Объем гексагональной призмы Упр и постоянная пространственной ячейки еь вычисляются по формуле Упр = Укуб 9/4 = где ^ - площадь гексагона базовой плоскоскости, состоящей из шести равносторонних треугольников со стороной аь.

В результате вычислений получены следующие значения: ^ = 10,37591387-10-20 м2; Упр = 50,79128721 10-30 м3; а еь = 4,895114575 • 10-10 м.

Гексагон БЕЕОШ (рис.1) окружают подобные гексагоны, им присвоены номера от 1 до 17. Нумерация атомов, находящихся в окружающих гексагонах, соответствует номерам гексагона. Расстояние между плоскостями в гексагональной призме равно 1/3 части ее высоты (еь/3 = 1,631704858-10-10 м).

Из таких призм состоит совершенный кристалл с-С, и они повторяются как вдоль главной оси, так и в перпендикулярных к этой оси направлениях неограниченно.

Рис.1. Расположение тетраэдров на базовой плоскости А; в призме А; В; С; А+

В призме можно выявить гипотетические тетраэдры [4]. Построение тетраэдра сводится к соединению отрезками прямых центров четырех сомкнутых сфер атомов углерода. Центры трех сомкнутых сфер на плоскости являются углами равностороннего треугольника - это основание тетраэдра, а центр четвертого атома, находящегося в лунке между атомами основания, оказывается вершиной тетраэдра. Все грани этого тетраэдра равны, а длина ребер соответствует постоянной ан призмы. Атомы, расположенные в вершинах тетраэдра, образуют по три внутренние связи с соседними атомами. На их образование каждый атом тетраэдра затрачивает по три валентности. При этом четвертая валентность каждого атома насыщается взаимодействием с атомами окружения - это могут быть атомы или тетраэдры.

Расположение тетраэдров на базовой плоскости Аг кристалла с-С показано на рис. 1. Правильный шестиугольник Б15Е15Е15О15Н15З15 вписан в окружность радиуса ан. Атомы углерода занимают центр и углы гексагона. Большие диагонали образуют шесть равносторонних треугольников, через их ортоцентры проведена окружность радиусом ^ортоц. Положение ортоцентров треугольников совпадает с центрами лунок, образованными атомами плоскости Аг. Атомы углерода, находящиеся в углах гексагона Б15, Е15 и Н15, являются вершинами тетраэдров, а их основания расположены на плоскости Вг выше плоскости Аг и обозначены штрихом. Вершинами тетраэдров гексагона Б^Е^Е^О^Н^З^ являются атомы углов Е9, О9 и З9, их основания, обозначенные двумя штрихами, расположены на плоскости Сг-1 ниже плоскости Аг. Гексагон БЕЕОНЗ - совмещенное изображение гексагонов 15 и 9. Углы гексагона являются вершинами тетраэдров, а их основания находятся выше и ниже плоскости Аг. Расположив одну над другой две такие сборки, можно получить структуру С1-1А1В1С1А1+1..., которая является рассматриваемой призмой. В ее объеме содержатся три тетраэдра с основаниями на плоскости Вг и три тетраэдра с основаниями на плоскости Сг. Вершины этих тетраэдров находятся на плоскостях Аг и Аг-+1 соответственно. В границах гексагона БЕЕОНЗ, расположенного на плоскости Аг, размещается только по одному углу каждого тетраэдра, а другие два угла каждого тетраэдра размещаются в прилегающих гексагонах.

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

Таблица 1

Параметры кристалла кубического углерода в гексагональной анаморфозе (пространственная постоянная ячейки а = 2,826195716 •Ю-10 м)

Фигура Параметр Значение

Тетраэдр Радиус описанной окружности основания Яоп, м Высота Н, м Длина ребра тетраэдра ан, м Орбиталь тетраэдра: - эксцентриситет еорб - фокусное расстояние сорб, м - малая полуось Ьорб, м - большая полуось /орб, м 1,153789571 • 10-10 1,631704856-10-10 1,998422156-10-10 0,816496581 0,999211078 10-10 0,706548929-10-10 1,223778643 10-10

Призма Радиус описанной окружности гексагона ан, м Высота сн, м 1,998422156-10-10 4,895114575 • 10-10

На рис.2 на фоне гексагона БЕЕОНЗ приведены внешние связи тетраэдров, обеспечивающие насыщение валентностей атомов, входящих в призму:

- связи между центральным атомом гексагона О и вершинами тетраэдров Н и Е, а также аналогичные связи в смежных гексагонах 1, 2, 4 и 5. Эти связи между атомами расположены на плоскости Аг и обозначены штрихпунктиром. Длина связей равна ан;

Рис.2. Расположение связей между тетраэдрами, тетраэдрами и центральным атомом

гексагона, а также между призмами:---- - связи между центральным атомом гексагона

О и атомами его большой диагонали H и E;---□ - связи между центральными атомами

базовых плоскостей призм А1Ч А,, А, - А j+i и т.д.; — • — проекция связи между атомами оснований тетраэдров, расположенных на плоскости В,, и атомами оснований тетраэдров, расположенных на плоскости C,

- связи между центральным атомом О гексагона DEFGHJ, расположенного на плоскости А,, с центральными атомами гексагонов поскостей Аг-1 и А,+1, которые находятся ниже и выше плоскости А,. Эти связи обозначены линиями с квадратом. Длина связей равна ch;

- связи между атомами оснований тетраэдров, находящихся на плоскости В,, c атомами оснований тетраэдров, которые находятся на плоскости C,. Эти связи обозначены линиями O'-D", F'-J"2, G'-O" и F'5-J". Длина связей равна геометрической сумме радиуса описанной окружности основания тетраэдра (Яоп = ^ортоц, см. рис.1) и расстояний между основаниями тетраэдров, расположенных на рассматриваемых плоскостях:

д/яоЛ + H2 = 1,998422156-10-10 м.

При построении всех связей учтен запрет на их взаимное пересечение.

Приведенные построения обеспечивают полное насыщение валентностей атомов углерода в тетраэдрах, гексагонах и призме, а также в кристалле с-С. Все связи проявляются как простые химические. Ординарность, направленность и обособление, а также отсутствие свободных зарядов позволяют сравнивать эти связи с простыми связями двухатомных молекул. Электронное облако таких молекул охватывает два положительных заряда, а его ось (прямая между ними) является осью симметрии бесконечного порядка с бесконечным числом плоскостей симметрии [5]. Взаимодействие электронов связи между собой (расталкивание) обусловливает легкость перехода из одной плоскости в другую. По каждой из этих плоскостей могут проходить силовые линии, совокупность которых имеет вид облака.

В отличие от электронов связей свободных молекул пересечения осей связей между атомами в с-С отсутствуют, а возможность бокового направления движения электронов связи в с-С ограниченна. Следствием этого является образование плоских траекторий их движения. Из этого следует, что связь в с-С (твердом диэлектрике с плотнейшей упаковкой) возможна только в виде плоской выпуклой замкнутой кривой, которая охватывает пару смежных атомов тетраэдра. Ребро тетраэдра является осью, а середина оси - центром. Орбиталь находится в плоскости, проходящей через ребро и центр тетраэдра. Подобное стеснение орбиталей наблюдается в алканах, например в октане. Боковые орбитали связей углерода с водородом ... - СН2 - СН2 - располагаются в плоскости, перпендикулярной главной оси.

Распределение потенциалей между двумя равными и одноименными зарядами. B работе [6] в координатах х, у приведена зависимость потенциала от расстояния между двумя равными и одноименными зарядами. Это кривая второго порядка - гипербола, ее фокусное расстояние с. Асимптотами являются прямые, проходящие через центр под углом в к оси абсцисс ^Р = 42 ). Вращение асимптот вокруг оси гиперболы ведет к образованию конуса с двухполостным гиперболоидом; вращение асимптот вокруг мнимой оси гиперболы выявляет конус с однополостным гиперболоидом.

Конические сечения - кривые второго порядка - в законченной форме отражены в канонических уравнениях гиперболы и эллипса. Их уравнения отличаются только знаком у второго члена (у эллипса - «плюс», у гиперболы - «минус»). Это и отражает их ортогональность и сопряженность - они расположены во взаимно ортогональных плоскостях.

Сечение конуса с однополостным гиперболоидом плоскостью, перпендикулярной его оси, выявляет окружность с радиусом, равным расстоянию от центра до заряда. Эта окружность является предельным случаем эллипса с нулевым эксцентриситетом. Другие сечения с углами в диапазоне 0° < а < 90° выявляют набор эллипсов с соответствующими эксцентриситетами. Сечение конуса плоскостью, перпендикулярной его образующей, выявляет эллипс, сопряженный с гиперболой с фокусным расстоянием сорб. Угол наклона а этой плоскости есть дополнительный к углу р. Его тангенс равен отношению малой полуоси орбитали к ее фокусному расстоянию: tgа = Ъорб/сорб = И42.

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

На рис.3 изображена искомая кривая и ее параметры. Это центральная кривая с центром Р, ее ось - ребро 0'0\ тетраэдра, вершиной которого является угол Б гексагона БЕЕОШ (см. рис.1, 2). Допустимо считать, что отрезок прямой от центра ребра Р до каждого атома О' и О'1 является фокусным расстоянием сорб, а вместе с продолже-

нием каждого отрезка до пересечения с ор-биталью в точках К и Ь - полуосью длиной /орб. Перпендикулярно к ней располагается малая ось МЫ и малые полуоси РМ и РЫ длиной ¿орб.

Искомая зависимость представлена в виде плоской замкнутой выпуклой кривой, охватывающей два заряда атомов ребра тетраэдра, расположенных от срединной точки ребра тетраэдра Р на расстоянии ^/2 =сорб; она симметрична относительно осей координат Р-х и Р-у. Можно представить, что пересечение кривой координатных осей происходит по нормали - это точки К и Ь на оси абсцисс и точки М и N на оси ординат.

Расстояние от срединной точки Р до точки М, а также точки N равно: Ьорб= сор бЛ/2 =

-10 г г

= 0,706548929-10 м, расстояние от срединной точки Р до точки К и Ь равно: /орб = 1,223778643 10-1° м. Эти расстояния соответствуют понятиям малой и большой полуоси кривой. Вследствие того, что орбиталь охватывает два заряда, она характеризуется соответствующим эксцентриситетом: еорб = сорб//орб = 0,816496581.

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

Следует отметить, что существование четырех точек искомой кривой не исключает изменений монотонности ее хода в диапазоне от центра Р до вершины Ь в промежутках. Для выявления таких «провалов» в кривой от фокусов ^ и Р2 к текущей точке Т могут быть проведены векторы, а их сумма |г1| + |г2| должна быть равна большой оси.

В табл.2 приведены рассчитанные значения для шести точек кривой, включая точку, соответствующую фокусному расстоянию. Все определяемые точки расположены в одном из квадрантов орбитали (см. рис.3), их положение является суммой расстояний до соответствующих фокусов и сравнивается с удвоенной длиной большой полуоси. Эти расстояния вычисляются по известным соотношениям:

г 1 = = /орб - еорб х; г 2 = Г^г = /орб + еорб х,

где х - координата (расстояние от центра орбитали / вдоль большой полуоси к вершине Ь).

Таким образом, параметры орбитали не только совпадают с параметрами эллипса, они удовлетворяют также его математическому определению.

Таблица 2

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

Рис.3. Схема для расчета параметров эллиптической орбитали в тетраэдре с-С

х ^рб х гь 10-10, м Г2, 10-10, м (г1 + Г2)'10-1°, м

0,1 0,081649658 1,142128985 1,305428301 2,447557286

0,2 0,163299316 1,060479327 1,387077959 2,447557286

0,4 0,326598632 0,89718001 1,550377275 2,447557285

0,6 0,489897948 0,733880694 1,673404812 2,447557285

0,8 0,653197264 0,570581378 1,876975908 2,447557285

0,999211078 0,815852428 0,407926214 2,039631071 2,447557285

Примечание. Исходные данные:

/орб = 1,223778643

10-10 м; в,

орб = 0,816496581; тх +г2 =

= 2/орб = 2,447557286-10-

м.

Заключение. Пространственное ограничение движения электронов связи в тетраэдрах кристалла с-С в направлении, перпендикулярном плоскости связи, в полной мере относится к орбиталям, охватывающим центральные атомы гексагонов A;-i - Ai, Ai -Ai+1 и т.д. Это связи между призмами, их положение соответствует главной оси кристалла. Ось орбитали совпадает с главной осью призмы, а ее плоскость перпендикулярна базовой плоскости, большому диаметру гексагона (см. рис.2, НЕ) и проходит перпендикулярно его малому диаметру.

Сопротивление кристалла с-С внешним нагрузкам зависит от количества связей и их энергии. Например, при растяжении призмы вдоль главной оси необходимо учитывать энергию связей между центральными атомами гексагонов базовых плоскостей и энергию связей между основаниями тетраэдров, расположенных на плоскостях Bi и Ci. К таким связям в рассматриваемой призме можно отнести три связи: две расположены в границах гексагона O'-D" и G'-O" и две F'-J'2 и F'5-J" в гексагонах 2 и 5. Энергия этих связей вместе с энергией связи между атомами базовых плоскостей меньше энергии связей других сечений призмы. Это, возможно, и выражается в меньшей прочности связей кристалла кубического углерода в направлении его главной оси.

Автор выражает благодарность А.Н. Неустроеву за полезную дискуссию и помощь в оформлении статьи.

Литература

1. Твердотельный автоэмиссионный диод / В.А. Беспалов. Э.А. Ильичев, А.Е. Кулешов и др. // Письма в ЖТФ. - 2013. - Т. 39. - Вып. 4. - С. 46-52.

2. Неустроев С.А. Влияние упаковки на воспроизводство сингонии кристалла кубического углерода // Изв. вузов. Электроника. - 2016. - Т. 21. - № 6. - С. 515-519.

3. Неустроев С.А. Уточнение параметров ячейки кристалла кубического углерода // Изв. вузов. Электроника. - 2014. - № 6 (110). - С. 86-87.

4. Неустроев С.А. Тетраэдрические связи кубического углерода // Изв. вузов. Электроника. -2010. - № 4 (84). - С. 86-88.

5. КрасновК.С. Молекулы и химическая связь. - М.: Высшая школа, 1984. - 295 с.

6. Говорков В.А. Электрические и магнитные поля. - M.: Энергия, 1968. - 487 с.

Поступила 17.08.2016 г.; принята к публикации 24.01.2017 г.

Неустроев Степан Архипович - доктор технических наук, профессор кафедры материалов функциональной электроники МИЭТ (Россия, 124498, г. Москва, г. Зеленоград, пл. Шокина, д. 1), doloto@list.ru

References

1. Bespalov V.A., Il'ichev Е.А., Kuleshov A.E., et al. Tverdotel'nyj avtoehmissionnyj diod [Solidsteit avtoemissive diod]. Pis'ma v ZHTF- Technical Physics Letters, 2013, vol. 39, iss. 4, pp. 46-52. (In Russian).

2. Neoustroev S.A. Vliyanie upakovki na vosproizvodstvo singonii kristalla kubicheskogo ugleroda [Influence of Packing on Reproduction of Singonia of c-C Crystal]. Izv. vuzov. Elektronika - Proceedings of universities. Electronics, 2016, vol. 21, no 6, pp. 515-519. (In Russian). (In Russian).

3. Neoustroev S.A. Utochnenie parametrov yachejki kristalla kubicheskogo ugleroda [More Precise Definition of Paramters in c-C Crystal]. Izv. vuzov. Elektronika - Proceedings of universities. Electronics, 2014, no. 6 (110), pp. 86-87. (In Russian).

4. Neoustroev S.A. Tetraehdricheskie svyazi kubicheskogo ugleroda [Cubic Carbon Tetrahedral Bonds] Izv. vuzov. Elektronika - Proceedings of universities. Electronics, 2010, no. 4 (84), pp. 86-88. (In Russian).

5. Krasnov K.S. Molekuly i khimicheskaya svyaz' [Moleculs and chemical bonds]. Moscow, Vysshaya shkola, 1984, 295 p. (In Russian).

6. Govorkov V.A. Elektricheskie i magnitnye polya [Electrical and magnetic fields]. Moscow, Energiya, 1968, 487 p. (In Russian).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Submitted 17.08.2016; accepted 24.01.2017

i Надоели баннеры? Вы всегда можете отключить рекламу.