Исследование одной математической модели распределения потенциалов в кристаллическом полупроводнике Текст научной статьи по специальности «Математика»

УДК 517.9

DOI: 10.14529/mmp190213

ИССЛЕДОВАНИЕ ОДНОЙ МАТЕМАТИЧЕСКОЙ МОДЕЛИ РАСПРЕДЕЛЕНИЯ ПОТЕНЦИАЛОВ В КРИСТАЛЛИЧЕСКОМ ПОЛУПРОВОДНИКЕ

Н.А. Манакова1, К.В. Васючкова1

1 Южно-Уральский государственный университет, г. Челябинск, Российская Федерация

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

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

К 60-летию со дня рождения выдающегося математика Яцека Банасяка.

Введение

На современном этапе развития полупроводниковой промышленности не теряет своей важности изучение физики полупроводников. Главная особенность полупроводников заключается в том, что их электропроводность резко увеличивается с повышением температуры. Математическое моделирование играет важную роль в исследованиях данных процессов, так как проведение реальных натурных экспериментов требует больших финансовых, трудовых и иных затрат (в некоторых случаях проведение натурных экспериментов невозможно, так как нет возможности контролировать отдельные параметры) [1-3]. Многие процессы, протекающие в веществах полупроводникового типа, можно описать дифференциальными уравнениями, неразрешенными относительно старшей производной по времени. Пусть П С - ограниченная область с границей дП класса C^. Математическая модель распространения потенциалов в кристаллическом полупроводнике может быть описана задачей Коши - Дирихле для неклассического уравнения в частных производных:

• фазовая переменная x £ Ck(0, T; M) является решением нелинейного уравнения соболевского типа

(А - A)xt - a1 Ax - a2div(|Vx|2Vx) + а3|x|p-2x = f, p > 2, (1)

где функция x = x(s,t) описывает потенциал электрического поля, заданная функция f = f(s, t) характеризует внешнее воздействие, например, внешнее электрическое поле, параметры А £ R, а1 а2 £ R,a3 £{0}U R+;

• удовлетворяет начальному условию Коши

х(в, 0) = х0(в), 5 € П; (2)

• удовлетворяет краевому условию Дирихле

х(М) = 0, (в,*) € дП х (3)

Уравнение (1) является многомерным обобщением модельных одномерных уравнений, описывающих волны малой частоты и длинноволнового характера, уравнениями данного типа описываются также «дрейфовые> волны в плазме, возникновение которых вызвано различными факторами. Так как в плазме и в полупроводниках возникают взаимодействия заряженных частиц посредством электромагнитного поля, то данные взаимодействия схожи по своей природе и поэтому явления, наблюдаемые при этих процессах, описываются уравнениями одного типа. Математическая модель физического процесса распределения потенциалов в полукристаллическом полупроводнике представлена в работе [4]. Причем, в случае а\ = 1, а2 = 1, аз = 0 было показано, что при наличии источников тока свободных зарядов или отрицательности дифференциальной проводимости существует пробой полупроводников (вследствие того, что потенциальная энергия системы превосходит кинетическую).

Уравнение (1) относится к классу полулинейных уравнений соболевского типа. Одним из основных методов изучения уравнений данного типа является метод фазового пространства, который впервые был предложен Г.А. Свиридюком в работе [5]. В дальнейшем такой подход к исследованию полулинейных уравнений соболевского типа был применен при исследованиях различных моделей математической физики [6-8]. Согласно методу фазового пространства сначала строится множество М всех допустимых начальных значений х0 (т.е. множество, состоящее из векторов, для которых существует локальное (единственное) решение задачи Коши). Фазовое пространство определяется как замыкание множества М. Затем находятся такие условия, при которых фазовое пространство является простым банаховым многообразием, а исходное уравнение соболевского типа в данной окрестности редуцируется к регулярному уравнению вида й = Г(и). Существование однозначного локального решения задачи Коши для получившегося регулярного уравнения является следствием классической теоремы, описанной в [9].

Целью нашей работы является построение условий, при которых фазовое пространство М уравнения (1) является простым банаховым многообразием. Напомним, что банахово многообразие называется простым, если его атлас эквивалентен атласу, содержащему одну карту. Ограничимся рассмотрением решений задачи (1) — (3), определенных только на М. Такие решения будем называть квазистационарными полутраекториями [10] (т.е. траекториями, которые проходят через точку х0 и лежат поточечно в фазовом пространстве М). В работе получены условия существования квазистационарной полутраектории уравнения (1) и продолжимости локального решения по времени.

1. Построение математической модели

о

Рассмотрим функциональные пространства N = Ж^П), В = ЬР(П), Н =Ь2(П), определенные в области П. Через В* и N обозначим сопряженные пространства к В

и N относительно скалярного произведения < •, • > в Н соответственно. Если п > 2 4

и 2 < р < 2 Н--в силу теоремы Соболева имеют место плотные и непрерывные

п — 2

вложения

N ^ В ^ Н ^ В* ^ Ш*. (4)

Определим в Н скалярное произведение формулой < х,у >= / ху^з, Ух, у € Н. Операторы Ь, М, N1, N определим следующим образом:

п

< Lx, y >= J (Axy + Vx ■ VyVx, y G N, п

< Mx, y >= — ai J Vx ■ Vyds, Vx, y G N, п

<Ni(x),y>= —.a2 J |Vx|2Vx ■ Vyds, Vx,y G N, п

<N2(x),y>= — a3 J |x|p-2xyds, Vx,y G B. п

По построению операторов и с учетом выбора функциональных пространств задача (1) — (3) редуцируется к задаче Коши

x(0) = xo (5)

для абстрактного полулинейного уравнения соболевского типа

Lx = Mx + Ni(x) + N2(x) + f, ker L = {0}. (6)

Обозначим через } последовательность собственных функций однородной задачи Дирихле для оператора (—А) в области П, а через {A&} - соответствующую им последовательность собственных значений, занумерованную по неубыванию с учетом кратности. Система собственных векторов } образует базис в пространстве

о

W^n) и в пространствах N, H в силу вложений (4). Тогда ker L = span{^i,..., <}, где {^i,..., <} - собственные векторы, соответствующие Ai, в случае A = Ai.

Замечание 1. В дальнейшем мы будем рассматривать диссипативные и р-коэрцитивные операторы. Оператор N2 называется диссипативным, если < N2 (x) — N2(y),x — y >< 0, Vx,y G B. Причем, диссипативность оператора N2 эквивалентна монотонности оператора (—N2). Заметим также, что из s-монотонности (т.е. < —N2xy, y >> 0 Vx,y G B) оператора (—N2) вытекает его строгая монотонность [11] и, как следствие, диссипативность.

Оператор N2 называется р-коэрцитивным [11], если 3 CN2, CN2 и 3 р > 2 такие, что < N2(x),x >> CN2 || x ||p и || N2(x) ||*< CN2 || x ||p-i Vx G B. Следует отметить, что из р-коэрцитивности оператора вытекает его сильная коэрцитивность.

Лемма 1. (i) При всех A > —Ai оператор L G L (N, N*) самосопряженный, фред-гольмов и неотрицательно определенный.

(ii) При всех ai G R+ оператор M G L (N, N*) диссипативный, а оператор (—M) 2-коэрцитивный.

(т) При всех а2,а3 Е оператор N Е С^ (ЭТ, диссипативный, а оператор (—N) 4-коэрцитивный. Оператор N Е Сте(В, В*) диссипативный, а оператор (—N) р-коэрцитивный.

Доказательство. Утверждение (г) является классическим результатом, так как

|<Ьу,у>|= / (Лу2 + Уу-Уу) <| Л Н1 У |Ц2(П) + II У Н\ < С II у Vy Е N. п

В утверждении (гг) диссипативность оператора М Е £ (ЭТ, ЭТ*) является следствием 5-монотонности оператора (—М):

< — Мху, у >= а^ У У у -Уу^> 0, Vx,y Е N х,у = 0, п

кроме того, в работе [11] была получена 2-коэрцитивность оператора (—М). (ггг) Диссипативность операторов N1 и N2 является следствием 5-монотонности операторов (—Ж1), (—соответственно:

< — у >= а Л [2(Ух ■ Уу )2 + |Ух|2Уу ■ Уу] ^ > 0, V», у Е N х, у = 0,

< —N2^у,у >= а3(р — 1)у |х|р 2у2 ^ > 0, Vx,y Е В, х,у = 0.

п

Кроме того, оператор (—N1) является 4-коэрцитивным, так как

< N(х), х >= а2 J |Ух|2Ух ■ Ух= а2 || х V» Е N. п

В [11] показана р-коэрцитивность оператора (—N2). □

2. Исследование фазового пространства

В силу самосопряженности и фредгольмовости оператора Ь отождествим N кегЬ = сокегЬ С ЭТ*, тогда N = сокегЬ 0 1шЬ и N = N° 0 где N° = кегЬ, N1 = со1ш Ь. Следовательно,

\УГ, если Л>-Ль , —-

1ш Ь = < к =1,1,

1{х Е N*| < х, >= 0}, если Л = —Л1,

| У1, если Л > —Л1, —-

со1ш Ь = < к =1,1.

1{х Е N < х, ^ >= 0}, если Л = —Л1,

Построим проектор Р в пространстве N вдоль N° на N и проектор Q в пространстве N вдоль сокег Ь на 1ш Ь:

Р = д = (/, если Л> —ЛЬ

I I — < >= 0}, если Л = —Л1.

Вестник !Ю"УрГ"У. Серия «Математическое моделирование 1 53

и программирование» (Вестник ЮУрГУ ММП). 2019. Т. 12, № 2. С. 150-157 53

В дальнейшем будем считать, что

(I — Q)/ не зависит от t G (0,T]. (7)

Построим множество

{ж G N : /(—aАж — a2div(|Vx|2Vx) + a3|ж|р-2ж)^ds = / /^kds}, Л = —Ai, M = п n

N, A > —Ai.

Определение 1. Вектор-функцию ж G C1((0,T]; N), k G N U {то} назовем решением уравнения (1), если она при некотором T G R+ удовлетворяет ему. Решение ж = x(t) уравнения (1) назовем решением задачи Коши (1), (2), если оно также удовлетворяет начальному условию (2).

Замечание 2. Немного отходя от стандарта [10], решение ж = x(t) уравнения (1) будем называть квазистационарной полутраекторией уравнения (1), если x(t) G M при всех t G (0,T]. Если квазистационарная полутраектория уравнения (1) удовлетворяет условию (3), то она называется квазистационарной полутраекторией уравнения (1), проходящей через точку жо.

Определение 2. [10] Пусть точка ж0 G M, положим ж1 = Рж0 G N1. Множество M в точке ж0 является банаховым Ck-многообразием, если существуют окрестности OM С M и O0 С N1 точек ж0 и ж1 соответственно и Ck-диффеоморфизм 6 : O0 ^ OM такой, что 6-1 равен сужению проектора Р на OM. Множество M называется банаховым Ck-многообразием, моделируемым подпространством N1, если оно является банаховым Ck-многообразием в каждой своей точке. Связное банахово Ck-многообразие называется простым, если любой его атлас эквивалентен атласу, содержащему единственную карту.

4

Теорема 1. Пусть п>2,2<р<2-\--, Л > —Ai, ai,a2,a3 G R+, и выполнено

n — 2

условие (7), тогда

(i) множество M является простым банаховым C1 -многообразием, моделируемым подпространством coimL;

(ii) существует единственная квазистационарная полутраектория ж G C1 ((0, M) уравнения (2), проходящая через точку ж0.

Доказательство. Рассмотрим доказательство утверждения (i). Так как уравнение (1) редуцируется к уравнению соболевского типа (6) и операторы M, N1 и N2 удовлетворяют условиям [11, Лемма 1.2], то построенное ранее множество M является простым банаховым C^многообразием всюду, за исключением, может быть, точки нуль. Идея доказательства основывается на свойствах диссипативности и р-коэрцитивности операторов M, N1 и N2, теореме Вишика - Минти - Браудера [12], а также теореме

0 неявной функции.

Рассмотрим доказательство утверждения (ii). В силу представления пространства N = N0 ф N1 и существования проекторов P и Q любой вектор можно представить в виде ж = ж0 + ж1, где (I — P)ж = ж0 G N0, а Рж = ж1 G N1. Так как оператор L линеен, то

Lж = L^1 + ж0) = L^1,

где - сужение оператора Ь на подпространство Ш1. Пусть О^ и О1 являются окрестностями точек хо € М и х0 € Рх0 соответственно. Так как множество М в точке х0 является простым банаховым С1 -многообразием, то существует диффеоморфизм 5 € С 1(О0; 0М). Следовательно, при каждой фиксированной точке х1 € 00 оператор 5х1 € ^(Ш1; ТхМ), где 5хх - производная Фреше оператора 5 в точке х1 € О,1,, а ТхМ - касательное пространство в точке х = 5(х1).

В окрестности 00 уравнение (1) редуцируется к регулярному уравнению х = Р(х), где Р = 5'х х Ь-1 ^(М+N1 + N2) : х ^ Тх М, тогда Р € С 1(0М; Т)М), где Т)М -сужение касательного расслоения ТМ на О^. В силу классической теоремы существования единственного решение задачи Коши [9] получим однозначную локальную разрешимость задачи (1) — (3).

Далее необходимо показать, что полученное решение является квазистационарной полутраекторией. Поскольку х = 5(х1), то х(£) € М при каждом £ € (0,Т). Вектор-функция х1(^) = Рх(£) удовлетворяет уравнению (1), так как 5'хх € ^(Ш1; ТхМ) -топлинейный изоморфизм. То есть, если х = х(£) - решение задачи (1) — (3), то х = х(£) - квазистационарная полутраектория уравнения (1), проходящая через точку х0. Рассмотрим стационарное уравнение

Му + ^(у) + N2 (у) + / = 0. (8)

В силу теоремы Вишика - Минти - Браудера [12] У/ € Ш* существует единственное

н-

решение уравнения (8). Введем в coim L норму |u|2 = (Lu, u) . В силу принципа Кура: та эта норма эквивалента норме, индуцированной из надпространства H. Применив свойство диссипативности операторов, получим

1 d 1 d

2 Jt^ = 2 Jt^X~ У|2 - {МХ + Nl{x) + N2{X) ~МУ + NlM + У) < 0,

где y = y (s) - решение уравнения (8), в свою очередь, являющееся стационарным решением уравнения (1), x = x(s, t) - решение уравнения (1). Отсюда, в силу единственности решения задачи Коши (1) — (3), получим продолжимость решения на интервал (0, +œ). □

Статья выполнена при поддержке Правительства РФ (Постановление № 211 от 16.03.2013 г.), соглашение № 02.А03.21.0011.

Литература

1. Al'shin, A.B. Blow-Up in Nonlinear Sobolev Type Equations / A.B. Al'shin, M.O. Korpusov, A.G. Sveshnikov. - Berlin: Walter de Gruyter, 2011.

2. Banasiak, J. Asynchronous Exponential Growth of a General Structured Population Model / J. Banasiak, K. Pichor, R. Rudnicki // Acta Applicandae Mathematicae. - 2012. - V. 119, № 1. - P. 149-166.

3. Banasiak, J. Asymptotic State Lumping in Transport and Diffusion Problems on Networks with Applications to Population Problems / J. Banasiak, A. Falkiewicz, P. Namayanja // Mathematical Models and Methods in Applied Sciences. - 2016. - V. 26, № 2. - P. 215-247.

4. Корпусов, М.О. Трехмерные нелинейные эволюционные уравнения псевдопараболического типа в задачах математической физики / М.О. Корпусов, А.Г. Свешников // Журнал вычислительной математики и математической физики. - 2003. - T. 43, № 12. -С. 1835-1869.

5. Свиридюк, Г.А. Многообразие решений одного сингулярного псевдопараболического уравнения / Г.А. Свиридюк // ДАН СССР. - 1986. - Т. 289, № 6. - С. 1-31.

6. Манакова, Н.А. Неклассические уравнения математической физики. Фазовые пространства полулинейных уравнений соболевского типа / Н.А. Манакова, Г.А. Свиридюк // Вестник ЮУрГУ. Серия: Математика, Механика, Физика. - 2016. - Т. 8, № 3. -С. 31-51.

7. Kondyukov, A.O. Phase Space of the Initial-Boundary Value Problem for the Oskolkov System of Nonzero Order / T.G. Sukacheva, A.O. Kondyukov // Computational Mathematics and Mathematical Physics. - 2015. - V. 55, № 5. - P. 823-828.

8. Свиридюк, Г.А. О складке фазового пространства одного неклассического уравнения / Г.А. Свиридюк, А.Ф. Карамова // Дифференциальные уравнения. - 2005. - Т. 41, № 10. - С. 1400-1405.

9. Leng, S. Introduction to Differentiable Manifolds / S. Leng. - New York: Springer, 2002.

10. Свиридюк, Г.А. Квазистационарные траектории полулинейных динамических уравнений типа Соболева / Г.А. Свиридюк // Известия РАН. Серия математическая. - 1994. -Т. 42, № 3. - С. 601-614.

11. Свиридюк, Г.А. Фазовое пространство уравнений типа соболева с s-монотонными и p-коэрцитивными операторами / Г.А. Свиридюк, М.В. Климентьев // Известия высших учебных заведений. Математика. - 1994. - Т. 38, № 11. - С. 72-79.

12. Гаевский, Х. Нелинейные операторные уравнения и операторные дифференциальные уравнения / Х. Гаевский, К. Грегер, К. Захариас. - М.: Мир, 1978.

Наталья Александровна Манакова, доктор физико-математических наук, доцент, кафедра «Уравнения математической физики>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), manakovana@susu.ru.

Ксения Владимировна Васючкова, аспирант, кафедра «Уравнения математической физики>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), vasiuchkovakv@susu.ru.

Поступила в редакцию 27 декабря 2018 г.

MSC 60H30 DOI: 10.14529/mmp190213

RESEARCH OF ONE MATHEMATICAL MODEL OF THE DISTRIBUTION OF POTENTIALS IN A CRYSTALLINE SEMICONDUCTOR

N.A. Manakova1, K.V. Vasiuchkova1

1 South Ural State University, Chelyabinsk, Russian Federation, E-mails: manakovana@susu.ru, vasiuchkovakv@susu.ru

The article is devoted to the research of the Cauchy problem for a mathematical model of the distribution of potentials in a crystalline semiconductor. By a semiconductor we mean

a substance with finite electrical conductivity, which rapidly increases with increase in the temperature. The mathematical model of the distribution of potentials is based on the semilinear Sobolev type equation supplemented by the Dirichlet and Cauchy conditions. We use the phase space method to construct sufficient conditions for the existence of the solution to the model under study. The conditions for the continuability of the solution are given.

Keywords: Sobolev type equations; mathematical model of distribution of potentials in crystalline semiconductor; phase space method; quasi-stationary trajectories.

References

1. Al'shin A.B., Korpusov M.O., Sveshnikov A.G. Blow-Up in Nonlinear Sobolev Type Equations. Berlin, Walter de Gruyter, 2011. DOI:10.1134/S001226610603013X

2. Banasiak J., Pichor K., Rudnicki R. Asynchronous Exponential Growth of a General Structured Population Model. Acta Applicandae Mathematicae, 2012, vol. 119, no. 1, pp. 149-166. DOI: 10.1007/s10440-011-9666-y

3. Banasiak J., Falkiewicz A., Namayanya P. Asymptotic State Lumping in Transport and Diffusion Problems on Networks with Applications to Population Problems. Mathematical Models and Methods in Applied Sciences, 2016, vol. 26, no. 2, pp. 215-247. DOI: 10.1142/S0218202516400017

4. Korpusov M.O., Sveshnikov A.G. Three-Dimensional Nonlinear Evolution Equations of Pseudoparabolic Type in Problems of Mathematical Physics. II. Computational Mathematics and Mathematical Physics, 2003, vol. 43, no. 12, pp. 1835-1869.

5. Sviridyuk G.A. The Manifold of Solutions of an Operator Singular Pseudoparabolic Equation. DAN SSSR, 1986, vol. 289, no. 6, pp. 1-31. (in Russian)

6. Manakova N.A., Sviridyuk G.A. Non-Classical Equations of Mathematical Physics. Phase Spaces of Semilinear Sobolev Equations. Bulletin of the South Ural State University. Series: Mathematics. Mechanics. Physics, 2016, vol. 8, no. 3, pp. 31-51. DOI:10.14529/mmph160304 (in Russian)

7. Kondyukov A.O., Sukacheva T.G. Phase Space of the Initial-Boundary Value Problem for the Oskolkov System of Nonzero Order. Computational Mathematics and Mathematical Physics, 2015, vol. 55, no. 5, pp. 823-828. DOI:10.1134/S0965542515050127

8. Sviridyuk G.A., Karamova A.F. On the Phase Space Fold of a Nonclassical Equation. Differential Equation, 2005, vol. 41, no. 10, pp. 1476-1481. DOI: 10.1007/s10625-005-0300-5

9. Leng S. Introduction to Differentiable Manifolds. N.Y., Springer, 2002.

10. Sviridyuk G.A. Quasistationary Trajectories of Semilinear Dynamical Equations of Sobolev Type. Izvestiya RAN. Seriya Matematicheskaya, 1994, vol. 42, no. 3, pp. 601-614.

11. Sviridyuk G.A., Klimentev M.V. Phase Spaces of Sobolev-Type Equations with s-Monotone or Strongly Coercive Operators. Russian Mathematics, 1994, vol. 38, no 11, pp. 72-79.

12. Gaevskiy H., Zaharias K. Nichtlineare operatorgleichungen und operatordifferentialgleichungen. Berlin, Akademie Verlag, 1974. (in German)

Received December 27, 2018