CC BY
17
Труды БГТУ, 2016, № 6, с. 101-105
101
УДК 537.633.2
В. Р. Мадьяров
Белорусский государственный технологический университет
ИЗМЕРЕНИЕ ПАРАМЕТРОВ ЭЛЕКТРОННОГО ПЕРЕНОСА В ПОЛУПРОВОДНИКАХ С ПОМОЩЬЮ ЭФФЕКТА ФАРАДЕЯ В МИЛЛИМЕТРОВОМ ДИАПАЗОНЕ
Применение магнитооптического эффекта для исследования полупроводников позволяет определить такие параметры, как эффективная масса, подвижность и время релаксации. В диапазоне частот 50-75 ГГц поглощение электромагнитных волн в объеме полупроводника практически полностью обусловлено свободными носителями, и в спектре поглощения отсутствуют области с аномальной дисперсией. Поэтому распространение плоско-поляризованной электромагнитной волны описывается системой линейных дифференциальных уравнений для немагнитной проводящей среды. Для продольного магнитного поля (геометрии Фарадея) получена система уравнений для угла поворота плоскости поляризации зондирующей волны, решение которой дает эффективное значение подвижности и время релаксации. Измерения угла поворота плоскости поляризации на пластинках n-Si и n-Ge показывают, что его зависимость от магнитной индукции близка к линейной, а постоянная Верде возрастает с повышением частоты. В исследуемом диапазоне частот ее величина изменяется в пределах 0,02-0,06 град/(мТл ■ мм) для n-Si и 0,02-0,20 град/(мТл ■ мм) для n-Ge. Значение эффективной массы для n-Si в направлении оси [100] хорошо согласуется с расчетным значением, полученным на основе известных данных для продольной и поперечной масс.
Ключевые слова: волны миллиметрового диапазона, эффект Фарадея, подвижность, удельная электропроводность.
V. R. Madiarov
Belarusian State Technological University
MEASUREMENT OF ELECTRONIC TRANSFER PARAMETERS IN SEMICONDUCTORS USING FARADAY EFFECT IN MILLIMETER-WAVE BAND
The application of magnetooptic effect for semiconductors studies allows to determine such parameters as effective mass, mobility and relaxation time. In the frequencies range of 50-75 GHz the absorption of electromagnetic waves in semiconductor volume is almost completely due to free carriers and there are no abnormal dispersion areas in absorption spectrum. Propagation of plane polarized electromagnetic wave is described by a system of the linear differential equations for nonmagnetic medium. For a longitudinal applied magnetic field (Faraday geometry) the equation system for rotation angle of the probing wave polarization plane is obtained the solving of which gives effective mobility and relaxation time values. The measurements of rotation angle on n-Si and n-Ge wafers show that its dependence on magnetic induction is very close to the linear, Verdet constant values increase with frequency. In the studied frequency range it varies in the region of 0.020.06 grad/(mTl ■ mm) for n-Si and 0.02-0.20 grad/(mTl ■ mm) for n-Ge. Effective mass value obtained from the measurements for n-Si in [100] axis direction is in good agreement with that predicted on the basis of the known values of longitudional and transverse masses.
Key words: millimeter waves, Faraday effect, mobility, specific conductivity.
Введение. Для исследования полупроводниковых материалов широко применяются бесконтактные методы с использованием излучения видимого диапазона, включая лазерно-интерференционный метод определения диффузионной длины [1], а также измерение скорости поверхностной рекомбинации и объемного времени жизни по кинетике избыточного теплового излучения [2]. В СВЧ-области применяются ре-зонаторные и волноводные методы [3]. Из всех методов с использованием зондирующих излу-
чений СВЧ-диапазона наиболее привлекательны те, которые основаны на зависимости амплитуды и фазы отраженной или прошедшей зондирующей волны от параметров полупроводникового образца: проводимости, подвижности, коэффициентов диффузии, объемной и поверхностной рекомбинации [4, 5]. Применение магнитооптического эффекта позволяет определить такие параметры, как эффективная масса и подвижность носителей заряда. В области частот ~1010-1011 Гц отсутствуют полосы
поглощения зондирующего излучения [6]. Поэтому распространение плоско-поляризованных зондирующих волн можно описать системой линейных дифференциальных уравнений Максвелла для проводящих сред, что упрощает интерпретацию данных измерений. В данной работе представлены результаты измерений подвижности и времени релаксации для германия и кремния с электронной проводимостью, полученные на основе измерения угла поворота плоскости поляризации в миллиметровом диапазоне длин волн.
Основная часть. В измерениях использовалось зондирующее излучение в диапазоне частот 50-75 ГГц. В этом диапазоне поглощение электромагнитных волн в объеме полупроводника практически полностью обусловлено свободными носителями, и в спектре поглощения отсутствуют области с аномальной дисперсией. Поэтому распространение плоско-поляризованной электромагнитной волны с частотой ю в направлении оси г описывается дифференциальными уравнениями, которые для немагнитной проводящей среды с удельной электропроводностью с имеют вид [7]:
Е = Е0 ехр 7(га / - кг), Н = Н0 ехр 7(га / - кг), юХЕ = -д(мН) = -/гац0 Н,
д/
(1)
ГОН = 88,
е
"а
- сЕ = 7Ю88„Е + сЕ,
где к - волновое число; е = е(ю) - диэлектрическая проницаемость среды, зависящая от частоты.
Исключив из уравнений (1) напряженность магнитного поля Н, для геометрии Фарадея (внешнее магнитное поле с индукцией В ориентировано в направлении распространения,^ а волновой вектор к перпендикулярен Е и Н волны) получим:
( г<2 \
,2 ®
к -8— + 7ц0юс С
Е = (к-Е)к = 0, (2)
где С
= иф
80М 0
скорость электромагнитной
волны в вакууме.
Удельная электропроводность с в уравнениях (1) и (2) определяется матрицей [7]:
(
с =
с -с
ху
\
уу У
где схх - проводимость, обусловленная дрейфом носителей под действием электрического поля зондирующей волны; сху - холловская проводимость полупроводникового образца, вызванная взаимодействием тока проводимости с постоянным магнитным полем с индукцией В.
Компоненты схх и сху в приближении сферических энергетических поверхностей определяются соотношениями
-7Ю
с = с =ст
хх уу 0
с = -с =ст
ху ух 0
(т 1 + 7га) +Ю2
га
(т 1 + /га)2 +га2 '
(3)
где с0 - удельная электропроводность на постоянном токе; т - время релаксации; юс = еВ / т -циклотронная частота; т - эффективная масса носителей заряда.
Для проекций электрического поля Ех и Еу на оси х и у получим:
С
к2 -8^
V С2 у
Ех + ^0га(сххЕх +схуЕу ) = 0,
(
га
2 Л
(4)
к2-8—
Еу + '^0га (-схуЕх + сххЕу ) =
В системе уравнений (4) с^ и сху выражаются формулами (3). Система (4) имеет решение, если определитель равен нулю:
(
к2 -8
га
V
-7^0гас х
+ ('М^гас ху )2= 0. (5)
Из (5) следует решение для к:
га
к+ =8~Г - 7М0гасхх ±М0гасху .
(6)
Введя полную проводимость с_ = с хх + /с ху, получим для правой и левой круговой поляризации:
2
к 2 = — К = С2
1
8„ - 7-
га 1 + /(га + гаС )т
(7)
где е, - диэлектрическая проницаемость в постоянном электрическом поле.
Угол поворота а плоскости поляризации после прохождения образца с толщиной й выражается с помощью соотношения (7) через углы правого и левого вращения а+ и а_:
а =1 (а+-а-), (8)
где
а± = Яе(к± )й =
2
гай
(82±+82± )2
"1±
с42
81± =8, -82± (га + гас)т = 8, -82± (гат + цВ)т, (9)
1
0га [1 + (га + гас)2 т2 ]
8пга
2
В формулах (9) частота юс выражена через эффективную подвижность носителей заряда ц = ет / т , в которую входят дрейфовая подвижность ц^ и холловская подвижность цн = гц^, где г - холл-фактор. Для частот ю ~ 3 • 1011 с-1 (~50 ГГц) квадратичное слагаемое в знаменателе (ю ± юс)2т2 << 1 и им можно пренебречь.
Используя формулы (8) и (9) и значения а, измеренные при двух значениях индукции магнитного поля В, получим систему уравнений для определения подвижности и времени релаксации:
а1 = 2[а+ (В1) -а-С^ХЬ
а2 = -[а+ (В2) -а-(В2)].
(10)
В качестве источника зондирующего излучения применялся генератор Г4-142 (диапазон частот 50-78 ГГц). Волновод имел стандартное прямоугольное сечение и был рассчитан на прохождение моды ТЕ11. Использовались пластинки образцов и-кремния и и-германия толщиной 0,5 мм. Образцы помещались в неподвижную цилиндрическую переходную латунную секцию, состоящую из двух коаксиальных цилиндров. Один из цилиндров мог поворачиваться и был соединен со второй прямоугольной секцией. На подвижную цилиндрическую секцию одевался соленоид, который создавал продольное магнитное поле в месте расположения образца. Плоскость образца кремния перпендикулярна оси [100], для которой известны продольная и поперечная эффективные массы. Направление магнитного поля совпадало с указанной осью. Угол поворота цилиндрической секции измерялся по смещению светового пятна лазерного луча, отраженного от зеркальца, жестко соединенного с вращающейся секцией. Мощность прошедшей волны измерялась детекторным блоком, соединенным со второй прямоугольной секцией. Предварительно проводилась проверка прохождения волны без образца. Регистрируемая мощность изменялась пропорционально соэ2ф (ф - угол между главными плоскостями подвижной и неподвижной секций волноводов). Такая же зависимость получалась и со вставленным образцом без приложенного магнитного поля. При смене направления магнитного поля путем переключения полярности источника питания соленоида направление поворота плоскости поляризации волны менялось на противоположное. Воспроизводимость измеренного угла поворота была не хуже 5%.
На рис. 1 и 2 представлены зависимости угла поворота плоскости поляризации а от магнитной индукции, измеренные при комнат-
ной температуре для образцов с различными удельными электропроводностями в диапазоне частот 53-63 ГГц.
а, град
В, мТл
а, град
0 10
0,6 0,4 0,21 0
20 а
30
40
1
2
3
В, мТл
0
10
20 б
30
40
Рис. 1. Зависимости угла магнитного вращения от магнитного поля для п-германия с удельной электропроводностью с0 = 86 См/м (а) и 12 См/м (б) на частотах: 1 - 63 ГГц; 2 - 56 ГГц; 3 - 53 ГГц
1,0 0,8 0,6 -I 0,4 0.2
0
а, град
В, мТл
0
10
20 30
40
0,25 0,20 0,15 0,10 0,05 0
б
Рис. 2. Зависимости угла магнитного вращения от магнитного поля для п-кремния с удельной электропроводностью с0 = 96 См/м (а) и 25 См/м (б) на частотах: 1 - 63 ГГц; 2 - 56 ГГц; 3 - 53 ГГц
а
Наблюдаемые зависимости практически являются линейными. Видно, что угол поворота плоскости поляризации значительно больше для образцов с более высокой удельной электропроводностью. Для германия измеренный угол а значительно больше, чем для кремния при одинаковой толщине. Видно также, что с ростом зондирующей частоты наклон этих зависимостей заметно возрастает.
На рис. 3 и 4 представлены зависимости постоянной Верде V = (1 / й)(Да / ДВ) от частоты в диапазоне 50-70 ГГц для образцов германия и кремния с различными удельными электро-проводностями. Значения V значительно больше для образцов с высокой электропроводностью. Для германия с более высокой удельной электропроводностью с0 значения V примерно в 4 раза превышают значения, полученные для кремния с близким значением с0. Для образцов с меньшей электропроводностью (с0 = 12 и 25 См/м) это различие заметно уменьшается, однако для «-Ое по-прежнему наблюдается большее значение V.
0,25 0,20
0,15 0,10 0,05
0
К град/(мТл • мм)
2
V, ГГц
45
50
55
60
65
70
Рис. 3. Зависимости постоянной Верде от частоты для п-германия: с0 = 86 См/м (1) и 12 См/м (2)
0,06 0,05 0,04.
0,03 0,02 0,01 0
К град/(мТл • мм)
1
2
V, ГГц
45
50
55
60
65
70
Рис. 4. Зависимости постоянной Верде от частоты для «-кремния: с0 = 96 См/м (1) и 25 См/м (2)
В таблице приведены значения эффективной подвижности цэф и времени релаксации т,
которые определялись с помощью системы уравнений (10) по значениям а, измеренным в магнитном поле с индукцией 15 и 30 мТл. В расчетах использовались значения е,, равные 12,4 и 16 для «-81 и «-Ое соответственно. Значения параметров ют и цВ не превышали 0,12 и 0,1 соответственно.
Значения эффективной подвижности времени релаксации т и эффективной массы т на частотах 53 и 63 ГГц
Данные на частотах 53/63 ГГц Дзс^ м2/(В • с) т, пс * т * т расч
ве 0,37/0,36 0,35/0,33 0,157/0,150т0 -
81 0,13/0,12 0,21/0,19 0,276/0,268т0 0,262т0
Значение эффективной массы т = ет / ц для кремния хорошо согласуется со значением т расч, полученным путем расчета по формуле
*
т =
V т1
_1
т.
где т1 = 1,59т0 и т2 = 0,19т0 - известные значения эффективной массы в направлении, параллельном и перпендикулярном оси [100] кремния соответственно (т0 - масса покоя электрона, т0 = 1,6 •Ю-19 кг).
Заключение. Показана возможность оценки подвижности и времени релаксации носителей заряда в «-81 и «-Ое с помощью магнитного вращения плоскости поляризации зондирующей волны в области частот 5070 ГГц. Доказано, что для образцов с толщиной около 0,5 мм угол поворота плоскости поляризации в «-Ое и «-81 в этом диапазоне частот пропорционален индукции внешнего магнитного поля. Поэтому для измерений на данной частоте можно применять калибровочный образец и использовать слабые магнитные поля.
Получены значения постоянной Верде, которые в исследуемом диапазоне составили 0,02-0,06 град/(мТл • мм) для кремния и 0,020,20 град/(мТл • мм) для германия. Значения эффективной массы т , полученные для кремния с помощью метода магнитного вращения, хорошо согласуются с расчетными для оси [100].
Для интерпретации результатов измерений с достаточной степенью точности можно использовать приближение сферических энергетических поверхностей.
-1
1
Литература
1. Манухов В. В., Федорцов А. Б., Иванов А. С. Лазерно-интерференционный метод определения длины диффузии носителей заряда в полупроводниках // Физика и техника полупроводников. 2015. Т. 49, вып. 9. С. 1153-1158.
2. Богатыренко В. В., Зиновчук А. В. Измерение скорости поверхностной рекомбинации и объемного времени жизни в пластинах Si по кинетике избыточного теплового излучения // Физика и техника полупроводников. 2011. Т. 45, вып. 1. С. 62-66.
3. Кучис Е. В. Гальваномагнитные эффекты и методы их исследования. М.: Советское радио, 1990. 263 с.
4. Сверхбыстродействующие фотоприемники на основе эффектов взаимодействия микроволнового электромагнитного излучения с фотовозбужденной плазмой в полупроводниках / В. В. Антонов [и др.] // Журнал технической физики. 1998. Т. 68, № 11. С. 94-99.
5. Кац Л. И., Альтшуллер Е. Ю., Чупис В. Н. Отражение электромагнитной волны от тонкой полупроводниковой пластины с управляемой ударной ионизацией концентрационной неоднородностью свободных носителей // Радиотехника и электроника. 1992. Т. 37, № 3. С. 560-566.
6. Уханов Ю. И. Оптические свойства полупроводников. М.: Наука, 1977. 368 с.
7. Ансельм А. И. Введение в теорию полупроводников. М.: Наука, 1978. 616 c.
References
1. Manukhov V. V., Fedortsov A. B., Ivanov A. S. Laser-interferometry method of determination of charge carrier diffusion length in semiconductors. Fizika i tekhnika poluprovodnikov [Semiconductor physics and engineering], 2015, vol. 49, issue 9, pp. 1153-1158 (In Russian).
2. Bogatyrenko V. V., Zinovchuk A. V. Measurement of surface recombination velocity and bulk lifetime in Si wafers by transient behavior of excess thermal emission. Fizika i tekhnika poluprovodnikov [Semiconductor physics and engineering], 2011, vol. 45, issue 1, pp. 62-66 (In Russian).
3. Kuchis E. V. Gal'vanomagnitnyye effekty i metody ikh issledovaniya [Galvanomagnetic effects and methods to study them]. Moscow, Sovetskoe radio Publ., 1990. 263 p.
4. Antonov V. V., Ivanov S. V., Tsarev V. P., Chupis V. N. Ultrahigh-speed photodetectors on the basis of interaction effects of microwave electromagnetic radiation with photoexcited plasma in semiconductors. Zhurnal tekhnicheskoy fiziki [Journal of technical physics], 1998, vol. 68, no. 11, pp. 9499 (In Russian).
5. Kats L. I., Al'tshuller Е. Yu., Chupis V. N. Reflection of an electromagnetic wave from thin semiconductor wafer with free carrier concentration nonuniformity controlled by impact ionization. Radiotekhnika i elektronika [Radio engineering and electronics], 1992, vol. 37, no. 3, pp. 560-566 (In Russian).
6. Ukhanov Ju. I. Opticheskiye svoystva poluprovodnikov [Optical properties of semiconductors]. Moscow, Nauka Publ., 1977. 368 p.
7. Ansel'm A. I. Vvedenie v teoriyu poluprovodnikov [Introduction to the theory of semiconductors]. Moscow, Nauka Publ., 1978. 616 p.
Информация об авторе
Мадьяров Владимир Рафкатович - кандидат физико-математических наук, доцент, доцент кафедры физики. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: [email protected]
Information about the author
Madiarov Vladimir Rafkatovich - PhD (Physics and Mathematics), Assistant Professor, Assistant Professor, the Department of Physics. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Поступила 02.03.2016
