Вакуумная наноэлектроника (ВНЭ) - новая область микро и наноэлектроники, получившая развитие в последние десятилетия. ВНЭ основана на использовании электронов в вакууме при размере активных элементов в десятые и сотые доли микрона. Практическим активным элементом в системах ВНЭ является автоэлектронный эмиттер. ВНЭ включает в себя микроскопические устройства и активные компоненты систем, соизмеримые с нанометровым масштабом. ВНЭ является своего рода альтернативной электроникой по отношению к твердотельной и использует управляемое баллистическое движение электронов в вакууме, в отличие от переноса носителей (электронов и дырок) в полупроводниках. Малость автоэмиттера позволяет миниатюризировать активные элементы вплоть до 10-8...10-9 м. В настоящее время есть экспериментальные результаты, показывающие, что для особых материалов (в частности, нанотрубок углерода и т.д.) размер активных элементов может составлять десятки атомных размеров. Малые размеры автокатода позволяют реализовать плотность элементов до 108 ...1010 см-2, а с использованием самоорганизующихся систем (фуллеренов, нанотрубок т. д.) до 1012 см 2.
Полевая эмиссия характеризуется следующими преимуществами [ ]:
1) отсутствие рассеяния энергии при переносе электронов в вакууме, которое обеспечивает преимущество при создании целого класса мощных, в том числе СВЧ приборов;
2) безынерционность автоэмиссионного процесса и короткий период переноса электронов в вакуумном зазоре, позволяющие создавать быстро-действующие высокочастотные устройства;
3) резкая нелинейность вольт-амперных характеристик автоэмиссии и возможность создания за счет этого преобразователей и умножителей частоты;
4) отсутствие энергетических затрат па эмиссионный акт вследствие квантовой природы туннелирования;
5) радиационная и термическая стойкость эмиттера.
Однако, полевая эмиссия также обладает и своими недостатками:
1) требуется высокая глубина вакуума, так как при сорбировании остаточных газов па поверхности автокатода увеличивается работа выхода материала эмиттера, что может значительно уменьшать эмиссионный ток;
2) сложность изготовления папоразмерных эмиттеров, а также их деформация в ходе работы, которая влечет приход всего устройства в негодность;
3) область высоких напряжений, где электроны приобретают наибольшую энергию, лежит в зоне острия эмиттера. В результате электрон получает вектор скорости, далеко не всегда совпадающий с вектором напряженности поля в приборе. Электронный пучок расширяется, по поскольку используется мпогоострийная или многолезвийная система, то общий поток электронов получается пересекающимся, из-за чего требуется прибегать к различным технологическим приемам для сужения и выпрямления пучка;
В настоящее время перспективными источниками электронов в вакуум- пых приборах па основе автоэлектронной эмиссии являются эмисионные катоды в виде мпогоэмиттерных массивов [ ]. В данной работе рассматривается ячейка тонкопленочного мпогоэмиттерного автокатода с лезвийной структурой эмиттеров вертикального типа [ ] в триодной копигурации (Рис. 1).
В данной работе описана математическая модель триодной ячейки массива автоэмиссионного катода, используя методы частицы в сетке и расщепления задачи на физические процессы. Модель включает в себя как формализацию, так и описание численных схем, необходимых для решения полученных задач.
Рассмотрен метод численного решения уравнений в частных производных эллиптического типа с помощью конечно-разностной аппроксимации высокого порядка точности на неравномерной сетке, также проанализированы методы решения полученной сеточной системы линейных алгебраических уравнений, в частности подробно изучен метод последовательной верхней релаксации, который является одним из наиболее эффективными при решении сеточных систем такого типа.
Рассмотрена схема численного интегрирования системы обыкновенных дифференциальных уравнений с перешагиванием. Схема устойчива на длинных траекториях, в силу того, что ошибки аппроксимации законов сохранения энергии и импульсов ограничены. Исходя из задачи интегрирования траекторий макрочастиц (39), данная схема подходит под специфику этой задачи, так как требуется рассчитывать относительно длинные траектории. Также в сравнении с другими схемами, метод с перешагиванием экономичен в вычислительном смысле, для обеспечения второго порядка точности по шагу необходимо только одно вычисление правой части системы обыкновенных дифференциальных уравнений.
В вычислительный пакет DAISI добавлены следующие модификации:
1) Расчёт плотности тока согласно модели Фаулера-Нордгейма, а также последующая инжекция частиц в расчётную область на основе эффекта полевой эмиссии.
2) Возможность задания начального распределения скоростей частиц.
3) Пакет адаптирован для запуска серий вычислительных экспериментов в высокопроизводительной вычислительной системе с общей памятью.
В ходе модификации пакета DAISI изучены способы организации объектно-ориентированного PIC-кода на примере уже существующего и работающего вычислительного пакета, написанного на C++.
Проведены следующие вычислительные эксперименты, используя модифицированный пакет DAISI:
1) Осуществлён обоснованный выбор вычислительной сетки.
2) Смоделировано начальное распределение скоростей, на основе экспериментальных спектров начальных энергий частиц.
3) Подобран коэффициент эффективного усиления поля в виде функции напряженности на основе экспериментальных данных.
4) Вычислены зависимость анодного тока от затворного напряжения и такие свойства ячейки, как потенциал в межэлектродном пространстве, напряженность электрического поля с учетом распределения заряда автоэлектронов в межэлектродном пространстве, плотность заряда и общий ток эмиссии (Рис. 17, 18, 19, 20, 21, 22)
Работа выполнена с использованием оборудования ресурсного центра Научного парка СПбГУ Вычислительный центр и Междисциплинарного ресурсного центра по направлению «Нанотехнологии» СПбГУ.
[1] Фурсей Г.Н. Автоэлектронная эмиссия СПб. Лань, 2012. 321 с.
[2] Егоров Н. В., Шешин Е. П. Автоэлектронная эмиссия. Принципы и приборы. М.: Интеллект, 2011. 704 с.
[3] Никифоров К. А., Егоров Н. В., Сайфуллин М. Ф. Математическое моделирование диодной системы с матричным автоэлектронным катодом // Журнал технической физики, 2015, Т. 85, С. 55-59.
[4] Алцыбеев В. В. Моделирование и оптимизация динамики интенсивных пучков заряженных частиц: дис. ... канд. физ.-мат. наук: 05.13.18. СПб., 2016. 130 с.
[5] Altsybeyev V. et al. Numerical simulation of a triode source of intense radial converging electron beam // Journal of Applied Physics, 2016. Vol. 120, No 14. P. 336-338.
[6] Feng Y., Verboncoeur J. P. A model for effective field enhancement for Fowler-Nordheim field emission // Phys. Plasmas 2005 Vol. 12, 103301.
[7] Fowler R.H., Dr. L. Nordheim Electron Emission in Intense Electric Fields // Proceedings of the Royal Society A., 1928, V. 119 (781), 173-181 c.
[8] Власов А. А. Теория многих частиц // M. : Наука, 1950, 350 с.
[9] Hockney R., Eastwood J. Computer Simulation Using Particles Francis, 1988. 540 p.
[10] Григорьев Ю.Н., Вшивков В.А., Федорук М.П. Численное моделирование методами частиц в ячейках Новосибирск:СО РАН, 2004, 360 с.
[11] Ponomarev, V., Altsybeyev, V. Development of 2D Poisson equation C++ finite-difference solver for particle-in-cell method //In "Stability and Control Processes"in Memory of V.I. Zubov (SCP), 2015, P. 195-197.
[12] В. В. Алцыбеев Оптимизационный алгоритм расчета плотности тока эмиссии // Вести. С.-Петербург, ун-та. Сер. 10. Прикл. матем. Информ. Проц, упр., 2015, А5 4, С. 56—71.
[13] Kevin L. Jensen Introduction to the Physics of Electron Emission John Wiley & Sons, 2017, 712 c.
[14] А. А. Самарский, П. H. Вабищевич, П. П. Матус Разностные схемы повышенного порядка точности на неравномерных сетках // Диффе- ренциальные уравнения, 1996. Т. 32, А5 2. С. 265-274.
[15] Вержбицкий В.М. Основы численных методов М. : Высш, шк, 2002. 840 с.
[16] Тихонов А. Н., Самарский А. А. Уравнения математической физики (5-еизд.). М. : Наука, 1977. - 742 с.
[17] Hong Q., et al. Why is Boris algorithm so good? // Phys. Plasmas, 2013. Vol. 20, no. 084503.
[18] G.L. Delzanno, E. Camporeale On particle movers in cylindrical geometry for Particle-In-Cell simulations // Journal of Computational Physics V. 253, 2013, P. 259-277.
[19] Поттер Д. Вычислительные методы в физике, М.: Мир, 1975. 392 с.
[20] Skeel R. D. Variable Step Size Destabilizes the Stomer/Leapfrog/Verlet Method // BIT Numerical Mathematics, 1993, Vol. 33, P. 172-175.
[21] R. G. Forbes, J. Deane Correction for Forbes and Deane, Reformulation of the standard theory of Fowler-Nordheim tunnelling and cold field electron emission // Mathematical, Physical and Engineering Sciences 2007, Vol. 463, No. 2087, 2907-2927 c.
[22] Зартдинов A. H., Никифоров К. А. Усиление электрического поля на наноструктурной эмиссионной поверхности
Труды 44-ой международной научной конференции аспирантов и студентов, С.-Петерб.: 2013, 215-219с
[23] Б. В. Гнеденко. Курс теории вероятностей. М.: ФИЗМАТГИЗ, 1961, 406 с
[24] Tsung-Wei Huang, Dian-Lun Lin, Yibo Lin, Chun-Xun Lin Cpp-Taskflow v2: A General-purpose Parallel and Heterogeneous Task Programming System at Scale arXiv, 2004, 19 c.
[25] A.B. Батраков, И.В. Пегелв, Д.И. Проскуровский Ограничение плотности тока автоэлектронной эмиссии пространственным зарядом эмитированных электронов // Писвма в ЖТФ, 1999, том 25, вып. 11, С. 78-82.