Диамагнетизм двухэлектронных примесных центров в полупроводниковых квантовых точках Текст научной статьи по специальности «Физика»

ФИЗИКА

УДК 535.8; 537.9; 539.33

В. Д. Кревчик, А. В. Разумов, П. С. Будянский

ДИАМАГНЕТИЗМ ДВУХЭЛЕКТРОННЫХ ПРИМЕСНЫХ ЦЕНТРОВ В ПОЛУПРОВОДНИКОВЫХ КВАНТОВЫХ ТОЧКАХ

Аннотация.

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

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

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

Выводы. Наличие квантового размерного эффекта приводит к уменьшению диамагнитной восприимчивости как одноэлектронных, так и двухэлектронных

110

University proceedings. Volga region

№ 1 (33), 2015

Физико-математические науки. Физика

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

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

V. D. Krevchik, A. V. Razumov, P. S. Budyanskiy

DIAMAGNETISM OF DIELECTRONIC IMPURITY CENTERS IN SEMICONDUCTOR QUANTUM DOTS

Abstract.

Background. A large number of impurity centers in semiconductors may possess not one, but two electrons in the bound state. The magnetic moment of such impurity centers, causing paramagnetism, equals to zero. In this connection the dielectronic impurity centers may possess only diamagnetic properties, researching of which allows to obtain valuable information on the pattern of semiconductor’s structural changes depending on changes in the content of the dielectronic impurity therein.

The interest to invrstigation of magnetic properties of quantum dots with dielectronic impurity centers is caused, first of all, by a new physical situation, associated with quantum size effect. From the practical point of view, such systems may be used in development of qubits and lazer structures on impurity transitions. The aim of the work is to theoretically investigate the features of diamagnetic susceptibility of a neutral dielectronic impurity center in a quantum dot, associated with size quanti-zetion. The authors compare the obtained results with a case of diamagnetic susceptibility of the one-electron D- -center in a quantum dot, as well as with a case of a neutral impurity center in a bulk semiconductor.

Materials and methods. Theoretical claculations of diamagnetic susceptibility were carried out for a case of the InSb quantum dot. On the basis of the two-electron model of the impurity center, the potential of which appears to be a zero radius potential, the authors conducted a variation calculation of the energy of the impurity’s main condition in a spherical quantum dot. The second potential of ionization, determining the energy of the main condition of such an impurity center in a quantum dot, was considered as an empirical parameter. Diamagnetic susceptibility of the two-electron impurity center was calculated by the Langevin-Pauli formula.

Results. It is shown that the diamagnetic susceptibility of a neutral impurity center in a quantum dot is lower the corresponding value in a bulk semiconductor. Herewith, with the decrease of the second potential of ionization their relation slowly increases, and with the increase of the quantum dot radius it tends to unity. It is revealed that the value of diamagnetic susceptibility of the two-electron impurity center in a quantum dot is vereal times greater than the value of diamagnetic susceptibility

of the one-electron D- -center, and with the increase of the quantum dot radius this difference becomes stronger. It is found that the diamagnetic susceptibility quite rapidly fades away at transition to deeper impurity centers.

Conclusions. The presence of the quantum size effect leads to the decrease of diamagnetic susceptibility of both one-electron and two-electron impurity centers at the transition “bulk semiconductor-quantum dot”. In quantum dots, the typical size of which exceeds the exciton’s Bohr radius, the contribution of two-electron impurity centers into diamagnetic susceptibility may be significant.

Key words: two-electron impurity centers, ionization potential, diamagnetic susceptibility, quantum dot.

Physics and mathematics sciences. Physics

111

Известия высших учебных заведений. Поволжский регион

Введение

Диамагнитная восприимчивость (ДВ) для двухэлектронных заряженных глубоких примесных центров (ПЦ) была рассчитана в работе [1] для случая объемного полупроводника в рамках модели короткодействующего потенциала. Была выявлена высокая чувствительность ДВ к величине второго

потенциала ионизации [1]. В случае квантовых ям ДВ одноэлектронных D -центров теоретически исследовалась в работе [2], где было показано, что за

счет более жесткого критерия перекрытия электронных орбит ДВ D- -центров в квантовой яме может значительно превосходить соответствующую величину в объемном полупроводнике. Интерес к исследованию ДВ двухэлектронных ПЦ в квантовых точках (КТ) обусловлен прежде всего новой физической ситуацией, связанной с квантовым размерным эффектом.

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

Расчет диамагнитной восприимчивости двухэлектронного примесного центра в полупроводниковой квантовой точке

Удерживающий потенциал КТ описывался в рамках модели «жестких стенок». На основе двухэлектронной модели ПЦ, потенциал которого представляет собой потенциал нулевого радиуса, нами ранее [3] проведен вариационный расчет энергии основного состояния нейтральной донорной примеси в сферической КТ. Второй потенциал ионизации £2 , входящий в параметр

f = yj!'-d(Ed - эффективная боровская энергия) и определяющий энергию основного состояния такого ПЦ после удаления одного электрона, принимался в качестве эмпирического параметра, заменяющего недостающую информацию о потенциале нулевого радиуса [4]. В работе [3] было получено выражение для энергии основного состояния двухэлектронного ПЦ в КТ, которое может быть записано через экстремальное значение величины f в виде

х

sh4 (R0f 1)(th (R0f 1)- R0f 1csch (R0f 1))

х -4f| 1R0 -8ch3(f 1R0)sh(f 1R0)Chi(2f 1R0) + 2ln(2exp(1 + C)f| 1R0)x xsh (2f-1R) + (Chi (4f-1R) - ln (2)) sh (4ff-1R0) - Shi (2f-1R) x

112

University proceedings. Volga region

№ 1 (33), 2015

Физико-математические науки. Физика

(l - 2ch (2n -R) - ch (4n-R)) + ch (4fl-R) Shi (4n-1Rq* )]},

(1)

здесь C - постоянная Эйлера; R* = Rq/aд ; Rq - радиус КТ; Chi(x) и Shi (x) - гиперболические интегральные косинус и синус соответственно; ад - эффективный боровский радиус.

Диамагнитная восприимчивость двухэлектронного ПЦ в полупроводниковой КТ может быть вычислена по формуле Ланжевена - Паули [5]:

x = Ndelad /_*2 . _*2

xq2 =----Г^Т r1 + r2

6m c '

(2)

где e - заряд электрона; m - эффективная масса электрона; с - скорость света в вакууме; Nd = 1/V, V - объем КТ; ц = ч/ад , г - координаты электронов (i = 1,2).

Задача вычисления ДВ сводится к расчету следующих интегралов:

Г*2 + Г? = jj¥*(r*,Ц)(r” + Ц*2)*(п.Ц2)dr*dr2 . (3)

(* * \

, Г2 ) имеет вид [3]:

¥(ц*, Ц* ) = ¥(ц* )^(r*). (4)

Для расчета ДВ используются одноэлектронные волновые функции в (4), записанные через экстремальное значение вариационного параметра n :

¥

(п )=4

Ron -pdn ) sin (xori sh (R*fi-1) sin (орд)

sh (R0n-1 - r*f\-1) *

- ri — Pd,

ri ^Pd,

(5)

здесь

*

B

sh (R0n 1

= i2n7 (th (R0n-1)- R0n-1csch ( -1)) ; X о = ^

о -n

-2

V* = V*/Ed ; csch(z) - гиперболический косеканс.

В результате (4) с учетом (1), (3) и (5) запишется в виде

¥

(* *\ r1 , r2 ) =

13 2

* * r1 r2

sh2 (1 -PdЛ 1) sin (x0r1* ) sin (x0Л, * *

r1 , r2 ^Pd,

sh

:(R*n-1

sin2

( oPd)

sh (RoTl 1 - ЛЛ 1) sh (R0n 1 - r2*n 1) * *

2 / * 1 \ , r1 , r2 — Pd.

sh2(Ron-1)

(6)

Physics and mathematics sciences. Physics

113

Известия высших учебных заведений. Поволжский регион

Тогда для (1) имеем

2 2 g 4 п 2п п 2п

Xq2 =- ^*02** ,4/ „ *д j d 02 sin 02 j d Ф2 j d 01 sin 01 j d Ф1

6т*о2 sh4 (Я°г\ 1)'

X

О О

X

sh4 № 1 _Pdfl 1 )pd *Pd. *

S1„4 (X„Pd ) j dr j dn(Г + '2 )Sin 0Г )Sin 0,2 )

+

* *

+ j dr1* j d,2* (,1*2 + r*2) sh2 (R0 ff1 _ ,1* ff1) sh2 (R„fi-1 _ ,2л-1)

pd pd

Для интегралов в (7) можно записать

pd pd

j d,1* j d,2 (r1*2 + r2*2)sin2 (5C0r1*) sm2 (X0r2* ) = 2_53_1%„4 x

(7)

0 0

X

-3 +

6x2pd +16X0p4. _ 24XM cos(2X0pd) + (3 _ 6x2p2 )cos (4X0pd)

+

+12X 0pdsin ( 0pd) _ 32X 0pdsin ( 0pd) + 6X 0pdsin (4X0pd)

(8)

;]< ;]< R «0

j dr* j d, (+ ?)sh2 («0л-1 - r*n-1)sh2 (R*n-1 _ ^-1) = 2-53_1n4 X

pd pd

X

3 + 24«02f|_2 + 16«04fi_4 _ 24R*fi_2pd _ 16«03fi_4pd _ 6fi_2pd

_16«0л_4p3 + 16ff4pd _ 24n_2 («0* _ pd)pdch(2Л-1 («0 _ pd)) + (3 + 6n_2pd)X Xch(4Л-1 («0 _ pd)) _ 24Л_1R°sh(2Л-1 («0* _ pd)) _ 8Л_3«°3sh(2fl-1 («0 _ pd +12л-1pd sh (2л-1 («0 _ pd)) _ 24nr3«0pd sh (2л-1 («0 _ pd

+32i 3pdsh (2П 1 («0* _pd))+6л 1pdsh (4л 1 («0* _pd

(9)

Тогда выражение для ДВ будет иметь вид

Nde2 a2d

Xq2 ^9 * 9

3 т о _г.4/п*~_1

2—2 п п

-X

sh

(«0п 1)(th («0п 1)_ «0fi 1csch («0тП 1

0

114

University proceedings. Volga region

№ 1 (33), 2015

Физико-математические науки. Физика

х <

sh4 (1 -pd fi 1 )х X() sin4 ( oPd)

3 + 6XoP2 + 16XM - 24XoPd cos (2%Qpd) +

+ (3 - 6X ciPd) cos (oPd) +12X oPd sin ( Qpd)- 32X ifcd sin ( Qpd) +

+6XoPd sin (Qpd )] + Л4 -3 + 24rq2ti 2 +16Rq4ti -24RqT1 2pd -

16Ro3fi 4Pd-6n 2pd- 16Ron 4pd + 16л 4pd-24fi2(R*-Pd)x

xPdch(2Л-1 (RQ- Pd)) + (3 + 6n"2Pd)ch(4Л-1 (R - Pd )) -

-24Л-1RQsh (2Ti-1 (R - Pd)) - 8f|-3RQ3sh (2Ti-1 (R - Pd)) +

+12л-1Pd sh (2л-1 ( rQ - Pd))- 24л-3R*Pd sh (2л-1 ( rQ - Pd))+ +32n-3Pdsh(2Л-1 (rQ - Pd)) + 6^-1Pdsh (4Л-1 (R*- Pd )) } •

(1Q)

Выполняя в (1Q) предельный переход Pd ^ 0 (Pd - радиус потенциальной ямы, в которой аппроксимировался короткодействующий потенциал [3]) для ДВ окончательно получим

Хд2 Хо

п

-2

sh4 (Ron 1)(th (Ron 1)- Ron 1csch (R*fi 1

-3 + 24 RQ2f|-2 +

+16RQ4^-4 + 3ch(4-n-1Rq*) - 8тП-1R* (3 + ff2R*2)sh (2f|-1R0) ,

(11)

здесь Xo =-Nd^/32m*cl •

Для сравнения выполним расчет ДВ D -центра в полупроводниковой КТ. Волновая функция электрона, локализованного на короткодействующем потенциале в центре КТ, имеет следующий вид [6]:

\

th (ni 1rQ )

2Щ1___________

-n11RQcsch (n-1R*)

х

х

sh (пД*) ch (n- 1Rq ) r* sh (n-1RQ) ’

(12)

где П1 =4\Ex\/Ed ; Ex - энергия связи D -состояния, отсчитываемая от дна КТ.

Physics and mathematics sciences. Physics

115

Известия высших учебных заведений. Поволжский регион

ДВ определяется следующим выражением:

J У* (г*) r*2^/)

*

dr .

(13)

Вычисление интеграла в (13) приводит к следующему выражению для Xq1:

На рис. 1 представлена зависимость отношения %^2/%v2 ДВ двухэлектронного ПЦ в КТ Xq2 и объемном полупроводнике %v2 от параметра п для

*

различных значений радиуса КТ Rq . Величина Xq2 была рассчитана с использованием формулы (11), а формула для xv2 была взята из работы [1]. Из рисунка видно, что ДВ нейтрального двухэлектронного ПЦ в КТ меньше соответствующей величины в объемном полупроводнике. При этом с уменьшением второго потенциала ионизации их отношение медленно возрастает, а с ростом радиуса КТ стремится к единице (ср. кривые 1 и 3). На рис. 2 приведена зависимость отношения %q2 /Xq1 ДВ двухэлектронного Xq2 и одно*

электронного %q1 ПЦ в КТ от радиуса КТ Rq для различных значений параметров п и п, характеризующего мощность потенциала нулевого радиуса [6]. Величина %q1 рассчитывалась с использованием формулы (14). Из рис. 2 видно, что величина ДВ нейтрального двухэлектронного ПЦ в КТ в несколько раз больше, чем величина ДВ одноэлектронного D -центра и с ростом радиуса КТ это различие усиливается. Можно также видеть, что величина ДВ достаточно быстро убывает при переходе к более глубоким примесным центрам (ср. кривые 1, 2).

В рамках модели потенциала нулевого радиуса для двухэлектронного примесного центра в КТ получена зависимость ДВ такого примесного центра от второго потенциала ионизации и радиуса КТ.

Показано, что наличие квантового размерного эффекта приводит к уменьшению ДВ как одноэлектронных, так и двухэлектронных ПЦ при переходе объемный полупроводник - квантовая точка. В КТ, характерный размер которых превышает боровский радиус экситона, вклад двухэлектронных ПЦ в диамагнитную восприимчивость может быть значительным.

%q1 4 th (m 1RQ) - П11RQcsch (m 1Ro*)

х(-4RQ3 - 6п2RQ + 3n3sh (2n-1RQ)).

m1 (csch (m 1r0 ))

(14)

Зависимость диамагнитной восприимчивости от второго потенциала ионизации и радиуса квантовой точки

Заключение

116

University proceedings. Volga region

№ 1 (33), 2015

Физико-математические науки. Физика

*

для различных значений радиуса КТ Щ ; кривые: 1 - 3; 2 - 2; 3 - 1

для различных значений параметров п и П; : кривые: 1 - п = 1, П; = 6; 2 - п = 1, П; = 8; 3 - п = 0,5, П; = 6

Physics and mathematics sciences. Physics

117

Известия высших учебных заведений. Поволжский регион

Список литературы

1. Белорусец, Е. Д. Диамагнетизм двухэлектронных глубоких примесных центров в полупроводниках / Е. Д. Белорусец, Э. З. Имамов, Б. Р. Маматкулов // Физика и техника полупроводников. - 1982. - Т. 16, № 1. - С. 89-92.

2. Imamov, E. Z. On Diamagnetism of Deep Impurities in Quantized Semiconductor Film / E. Z. Imamov, V. D. Krevchik // Phys. St. Sol. (b). - 1982. - Vol. 114. - P. 201207.

3. Кревчик, В. Д. Двойная фотоионизация двухэлектронных примесных центров в квазинульмерных наноструктурах / В. Д. Кревчик, А. В. Разумов, П. С. Будянский // Известия высших учебных заведений. Поволжский регион. Физикоматематические науки. - 2014. - № 3 (31). - С. 228-240.

4. Гринберг, А. А. Об энергетическом спектре многозарядных примесных центров в полупроводниках / А. А. Гринберг, Е. Д. Белорусец // Физика и техника полупроводников. - 1978. - Т. 12, № 10. - С. 1970-1978.

5. Флюгге, З. Задачи по квантовой механике / З. Флюгге. - М. : Мир, 1974. - Т. 1.

6. Кревчик, В. Д. Энергетический спектр и оптические свойства комплекса квантовая точка - примесный центр / В. Д. Кревчик, А. В. Левашов // Физика и техника полупроводников. - 2002. - Т. 36, № 2. - С. 216-220.

References

1. Belorusets E. D., Imamov E. Z., Mamatkulov B. R. Fizika i tekhnika poluprovodnikov [Physics and technology of semiconductors]. 1982, vol. 16, no. 1, pp. 89-92.

2. Imamov E. Z., Krevchik V. D. Phys. St. Sol. (b). 1982, vol. 114, pp. 201-207.

3. Krevchik V. D., Razumov A. V., Budyanskiy P. S. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2014, no. 3 (31), pp. 228-240.

4. Grinberg A. A., Belorusets E. D. Fizika i tekhnika poluprovodnikov [Physics and technology of semiconductors]. 1978, vol. 12, no. 10, pp. 1970-1978.

5. Flyugge Z. Zadachi po kvantovoy mekhanike [Problems on quantum mechanics]. Moscow: Mir, 1974, vol. 1.

6. Krevchik V. D., Levashov A. V. Fizika i tekhnika poluprovodnikov [Physics and technology of semiconductors]. 2002, vol. 36, no. 2, pp. 216-220.

Кревчик Владимир Дмитриевич

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

E-mail: physics@pnzgu.ru

Разумов Алексей Викторович

кандидат физико-математических наук, доцент, кафедра общей физики и методики обучения физике, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)

E-mail: physics@pnzgu.ru

Krevchik Vladimir Dmitrievich Doctor of physical and mathematical sciences, professor, dean of the faculty of instrument engineering, information technology and electronics, Penza State University (40 Krasnaya street, Penza, Russia)

Razumov Aleksey Viktorovich Candidate of physical and mathematical sciences, associate professor, sub-department of general physics and physics teaching technique, Penza State University (40 Krasnaya street, Penza, Russia)

118

University proceedings. Volga region

№ 1 (33), 2015

Физико-математические науки. Физика

Будянский Павел Сергеевич Budyanskiy Pavel Sergeevich

аспирант, Пензенский государственный Postgraduate student, Penza State

университет (Россия, г. Пенза, University (40 Krasnaya street,

ул. Красная, 40) Penza, Russia)

E-mail: physics@pnzgu.ru

УДК 535.8; 537.9; 539.33 Кревчик, В. Д.

Диамагнетизм двухэлектронных примесных центров в полупроводниковых квантовых точках / В. Д. Кревчик, А. В. Разумов, П. С. Будянский // Известия высших учебных заведений. Поволжский регион. Физикоматематические науки. - 2015. - № 1 (33). - С. 110-119.

Physics and mathematics sciences. Physics

119