CC BY
9
IJDC 51-73+537.2 Вестник СПбГУ. Прикладная математика. Информатика... 2021. Т. 17. Вып. 2 MSC 35J05, 78А30
Mathematical modeling of triode system on the basis of field emitter with ellipsoid shape*
N. V. Egorov, E. M. Vinogradova
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Egorov N. V., Vinogradova E. M. Mathematical modeling of triode system on the basis of field emitter with ellipsoid shape. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2021, vol. 17, iss. 2, pp. 13H36. https://doi.org/10.21638/11701/spbul0.2021.203
In this paper the mathematical modeling of the triode emission axially symmetric system on the basis of field emitter is considered. Emitter is an ellipsoid of revolution, anode is a confocal ellipsoidal surface of revolution. Modulator is a part of the ellipsoidal surface of revolution, confocal with the cathode and anode surfaces. The boundary-value problem for the Laplace's equation in the prolate spheroidal coordinates with the boundary conditions of the first kind is solved. The variable separation method is applied to calculate the axisymmetrical electrostatic potential distribution. The potential distribution is represented as the Legendre functions expansion. The expansion coefficients are the solution of the system of linear equations. All geometrical dimensions of the system are the parameters of the problem.
Keywords: micro- and nanoelectronics, field emitter, field emission, mathematical modeling, electrostatic potential, boundary-value problem, Legendre functions.
1. Introduction. Field emitters are extensively used as electron sources of high current density for vacuum micro- and nanoelectronic devices [1-5]. This work is devoted to the mathematical modeling of an axially symmetric triode system on the basis of field emitter with ellipsoidal shape. The cathode's and anode's shapes are the confocal ellipsoidis of revolution. The modulator's shape is the part of the confocal ellipsoid of revolution (Figure). To calculate the axially symmetric electrostatic potential distribution the variable separation method in the prolate spheroidal coordinates (а, в) is used. In accordance with the theory of mathematical physics the boundary value problems with the mixed boundary-conditions for the bounded domain can be reduced to the linear algebraic equations system with constant coefficients [6, 7].
2. Mathematical model. The problem parameters: а = ai (в € [0, я]) — cathode surface; a = a3 (в € [0, я]) — anode surface; a = a2 (в € [во, я]) — modulator surface; U(1; ) = 0 — cathode boundary condition; U (а2,в) = f (в), (в € [в0, я]) — modulator boundary condition; U(a3,в) = 0 — anode boundary condition.
The electrostatic potential distribution U(а,, в) is a solution of the boundary-value problem for the Laplace's equation:
1
c2(sin в + sinh a)
d2U , dU d2U 3U
—+cotha—+ —+ cot/3—
(1)
* This work was supported by the Russian Foundation for Basic Research (project N 20-07-01086). © St. Petersburg State University, 2021
а3
а2
(Сш Л ///Л 1 ш
uZS / 1 >< J z
Figure. Schematic representation of the triode field emission system and boundary conditions:
U (ai ,в) = 0, U (as ,в) = 0, U (a2 ,в) = f (в),
в e [0, n], в e [0, n], в e [во, П]
(2)
where c — focal radius.
3. Solution of the problem. To solve the boundary-value problem (1), (2) the internal area of the system can be divide into following subdomains:
(1)-(a e [ai,a2], 3 e [0, n]);
(2) - (a e [a2, a3], /3 e [0, n])
(3) — (a e [ai, as], 3 e [0, 3o]).
Let
U(a, в) =
Ui(a, в), a e [ai, a2], в e [0, n] , U2(a, в), a e [a2, as], в e [0, n] , Us (a, в), a e [ai, as ], в e [0, во ] •
(3)
Making use of the variable separation method for the (i)-th subdomain, the potential distribution (3), Ui (a,/) (i =1, 2, 3) accordingly, can be represented as [8-11]:
n=o Wn (ai ,a2)
тт ( r\ sr л тэ r Wn(a,a3)
n=o Wn (a2 ,as)
(4)
Us (a, /3) =52 Bmw -i/2+irm (ai ,a)
m=i
-P— 1 /2-\-iTm (COS
l/2+гтm (COS (3o) '
where [12-14]
Wn(x,y) = P„(coshx)Q„(coshy) - P„(coshy)Q„(coshx), W-1/2+iTm (x,y) = P-i/2+iTm (cosh x)Q-l/2+iTm (coshy)-- P-i/2+iTm (cosh y)Q-i/2+iTm (cosh x);
Pn(cos¡3) — Legendre polynomials; P-1/2+iTm(cosx) — conical functions; Pn(coshx), P-i/2+irm (cosh x) — Legendre functions of the fi rst kinds; Qn (cosh x), Q-1/2+iTm (cosh x) — Legendre functions of the second kinds; rm — roots of a linear combination of Legendre functions of the first and second kind:
W-1/2+iTm (ab a3) = 0.
(6)
The potential distribution U(a, ¡3) (3) written in the form of series in eigenfunctions (4) satisfies the homogeneous boundary conditions on the cathode and anode surfaces U(ai,3)=0, U(as,3)=0 (3 e [0,
In order to fulfil the potential distribution continuity condition and to found the coefficients A^d Bm in the expantions (4)-(6) the additional boundary conditions can be used:
Ui(a2 ,в) = и2(а2,в)
U3(a2,0), 0 < в<во, f (в), в<во < п,
(7)
U3(a, во) =
Ui(a, во), ai ^ a ^ а2, U2(a, во), «2 < а < аз.
(8)
Conditions (7) and (8) provide not only potential distribution continuity, but also the continuity of its first derivative due to the fact that the boundary points of subdomains (1), (2) at a = a2 are internal points of subdomain (3) for 0 < ¡3 < ¡30, and that the boundary points of subdomain (3) at 3 = (a1 < a < a3) are internal points of subdomain (1) for a1 < a < a2 and subdomain (2) for a2 < a < a3.
Making use of formulas (1), (4), (7), (8) it turns out that
v^ I V BmW_1/2+iT (aua2)^ V2+ÎT (cos/3)
Y.AnPnicos/3) = 4 m l/2+lTmK 2)P_1/2+iTJcœe0y (9)
m=1
f (в),
¡3 <во < п,
BmW_1/2+iTm{al,a) = <
m=1
Ели/ a \ Wna)
AnPn cos/3o —-г..
Wn(a1, a2)
An Pn (cos во )
Wn(a, a3)
Wn(a2 ,a3)''
ai ^ a ^ a2,
a2 ^ a ^ a3.
(10)
The orthogonality of the eigenfunctions systems in the variables ¡3 (9) and a (10) leads to the linear algebraic equations system
An~(n+1/2)
P-i/2+iTm (cos ¡o)
Io
^y sin(3)P„ (cos 3)P-!/2+iTm (cos 3) = 0
n
= (n +1/2) J sin(3)P„ (cos 3)/(3) d3, (11)
lo
1 OO
Bm - J-f- AnPn (COS 100) X
n=0
«2
——-- / sinh(a) /2+jT (ai, a)Wn(o;i, a) da+
Wn(ai,a2) J
~ 0, (12)
+ T„ , 1-г [ sinh(a)W_1/2+i (a1,a)Wn(a, a3) da
Wn(a2, аз) J
where Nm = j sinha (W_1/2+iTm(ai, a))2 da — eigenfunction normalization constant.
4. Conclusion. This paper presents the mathematical model of the rotationally symmetrical field emission triode system on the base of field emitter with ellipsoid surface. The anode surface is an ellipsoid of revolution, the modulator surface is a part of an ellipsoid of revolution. To calculate the potential distribution, the boundary problem for the Laplace equation (1), (2) was solved by the separation variables method. The internal region of the emission system was divided into three overlapping subdomains. In each of the subdomains, the electrostatic potential distribution was presented in the form of series by the Legendre functions (3)-(6). Additional conditions (7), (8) ensuring the potential continuity and its first derivative made it possible to reduce the initial boundary problem to a system of linear equations (11), (12) relative to the unknown coefficients in the potential expansion (4). Thus, the potential distribution is found in the entire area of the field emission triode system.
References
1. Bachmann M., Dusberg F., Langer C., Herdl F., Bergbreiter L., Dams F., Miyakawa N., Eggert T., Pahlke A., Edler S., Prommesberger C., Lawrowski R., Hausladen M., Schreiner R. Vacuum-sealed field emission electron gun. Journal of Vacuum Science and Technology B. Nanotechnology and Microelectronics, 2020, vol. 38, iss. 2, no. 023203.
2. Tripathi P., Gupta B. K., Bankar P. K., More M. A., Late D. J., Srivastava O. N. Graphene nanosheets assisted carbon hollow cylinder for high-performance field emission applications. Materials Research Express, 2019, vol. 6, iss. 9, no. 095066.
3. Jones W. M., Lukin D., Scherer A. Practical nanoscale field emission devices for integrated circuits. Applied Physics Letters, 2017, vol. 110, iss. 26, no. 263101.
4. Sominskii G. G., Tumareva T. A., Taradaev E. P., Rukavitsyna A. A., Givargizov M. E., Stepa-nova A. N. Annular multi-tip field emitters with metal-fullerene protective coatings. Technical Physics, 2019, vol. 64, iss. 2, pp. 270-273.
5. Nguyen H. D., Kang J. S., Li M., Hu Y. High-performance field emission based on nanostructured tin selenide for nanoscale vacuum transistors. Nanoscalc, 2019, vol. 11, iss. 7, pp. 3129-3137.
6. Vinogradova E. M., Egorov N. V., Televnyi D. S. Calculation of a triode field-emission system with a modulator. Technical Physics, 2014, vol. 59, iss. 2, pp. 291-296.
7. Egorov N. V., Vinogradova E. M. Mathematical modeling of the electron beam formatting systems on the basis of field emission cathodes with various shapes. User Modeling and User-Adapted Interaction, 2003, vol. 72, iss. 2, pp. 103-111.
8. Majic M. R. A., Gray F., Auguié B., Ru E. C. L. Electrostatic limit of the T-matrix for electromagnetic scattering: Exact results for spheroidal particles. Journal of Quantitative Spectroscopy and Radiative Transfer, 2017, vol. 200, pp. 50-58.
9. Majic M. R. A., Gray F., Auguié B., Ru E. C. L. Spheroidal harmonic expansions for the solution of Laplace's equation for a point source near a sphere. Physical Review E, 2017, vol. 95, iss. 3, no. 033307.
10. Techaumnat B., Huynh V., Hidaka K. Three-dimensional lectromechanical analysis of a conducting prolate spheroid on a grounded plane. IEEE Transactions on Dielectrics and Electrical Insulation, 2014, vol. 21, iss. 1, pp. 80-87, no. 6740728.
11. Xue C., Edmiston R., Deng S. Image theory for Neumann functions in the prolate spheroidal geometry. Advances in Mathematical Physics, 2018, vol. 2018, no. 7683929.
12. Abramowitz M., Stegun I. A. Handbook of mathematical functions with formulas, graphs, and mathematical tables. Washington, DC, NBS Press, 1964, 1046 p. (Applied Mathematics. Series 55).
13. Bateman H., Erdelyi A. Higher transcendental functions. In 2 vol. New York, Toronto, London, McGraw-Hill Book Company Inc., 1953, vol. 1, 302 p.
14. Gradshteyn I. S., Ryzhik I. M. Table of integrals, series, and products. 7th ed. Amsterdam, Boston, Heidelberg, London, New York, Oxford, Paris, San Diego, San Francisco, Singapore, Sydney, Tokyo, Academic Press, 2007, 1171 p.
Received: December 18, 2020.
Accepted: April 05, 2021.
Authors' information:
Nickolay V. Egorov — Dr. Sci. in Physics and Mathematics, Professor; [email protected]
Ekaterina M. Vinogradova — Dr. Sci. in Physics and Mathematics, Professor; [email protected]
Математическое моделирование триодной системы на основе полевого катода эллипсоидальной формы*
Н. В. Егоров, E. М. Виноградова
Санкт-Петербургский государственный университет, Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7-9
Для цитирования: Egorov N. V., Vinogradova E. M. Mathematical modeling of triode system on the basis of field emitter with ellipsoid shape // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2021. Т. 17. Вып. 2. С. 13И36. https://doi.org/10.21638/11701/spbul0.2021.203
Данная статья посвящена математическому моделированию осесимметричной триодной эмиссионной системы на основе полевого эмиттера. Эмиттер представляет собой эллипсоид вращения, анод — конфокальную эллипсоидальную поверхность вращения, модулятор — часть эллипсоидальной поверхности вращения, конфокальной с поверхностями катода и анода. Граничная задача решена для уравнения Лапласа в вытяну-
* Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (грант № 20-07-01086).
тых сфероидальных координатах с граничными условиями первого рода. Для расчета распределения электростатического потенциала был использован метод разделения переменных. Распределение потенциала представлено в виде разложений по функциям Лежандра. Коэффициенты разложений являются решением системы линейных алгебраических уравнений. Все геометрические размеры системы являются параметрами задачи.
Ключевые слова: микро- и наноэлектроника, полевой эмиттер, полевая эмиссия, математическое моделирование, электростатический потенциал, граничная задача, функции Лежандра.
Контактная информация:
Егоров Николай Васильевич — д-р физ.-мат. наук, проф.; [email protected]
Виноградова Екатерина Михайловна — д-р физ.-мат. наук, проф.; [email protected]
