Объяснение дисперсионных свойств различных материалов является одним из важных фундаментальных разделов курса общей физики, преподаваемого в вузах. Это объяснение базируется на классической теории Максвелла. Дисперсия, как зависимость показателя преломления от длины волны (частоты) падающего света, связана с неодинаковой диэлектрической проницаемостью вещества для света различной длины волны. Традиционно в лекционном курсе излагается общая теория, объясняющая как нормальную, так и аномальную дисперсию вынужденными колебаниями заряженных гармонических осцилляторов под действием электрического поля электромагнитной волны. Однако методическая общность такого теоретического описания нарушается в лабораторном практикуме, который предлагает студентам только работу по изучению нормальной дисперсии. Создание лабораторной работы по изучению аномальной дисперсии позволило бы не только добиться методической целостности рассмотрения дисперсии в курсе общей физики, но и подчеркнуть особенности методики исследований, связанных с изучением аномальной дисперсии.
Поскольку аномальная дисперсия наблюдается вблизи полосы поглощения, для ее экспериментального изучения в курсе общей физики необходимо выбрать вещества, имеющие достаточно широкую и изолированную полосу поглощения. Многие полупроводники имеют такую полосу в инфракрасной области. Для обоснования методики экспериментального исследования рассмотрим одну из физических моделей, описывающую дисперсионные свойства диэлектрической проницаемости, а значит и показателя преломления, полупроводников в инфракрасной области - однофононный резонанс.
Физическая модель
В отличие от нормальной дисперсии, где основным механизмом изменения диэлектрической проницаемости вещества принято считать вынужденные колебания оптического электрона атома [2], одной из причин аномальной дисперсии в полупроводниках можно считать вынужденные колебания ионов кристаллической решетки [4].
Элементарная ячейка алмазоподобных полупроводников состоит из двух ионов массами mA и mB с зарядами +q и -q, соответственно. Задача об относительных колебаниях ионов традиционным образом сводится к колебаниям материальной точки массой m = mAmB/(mA + mB) (приведенная масса) и зарядом q относительно неподвижного центра инерции [1]. Благодаря строению кристаллической решетки, тепловые (собственные) относительные колебания ионов происходят с одной разрешенной частотой wT, соответствующей инфракрасному диапазону спектра (при нормальных условиях), и распространяются по кристаллу в виде поперечных волн - фононов.
В классическом приближении можно считать, что относительное смещение ионов, возникающее вследствие вынужденных колебаний под действием внешней гармонической силы (электрического поля
электромагнитной волны), описывается уравнением движения:
или
где - вынуждающая гармоническая сила, действующая со стороны электрического поля
на заряд q;
- суммарная квазиупругая сила, удерживающая ионы в узлах решетки и заставляющая их совершать собственные колебания с частотой wT. Колебания ионов распространяются по кристаллу в виде упругих волн, что приводит к рассеянию энергии, которое можно описать как работу силы «сопротивления»
, вызывающей затухание собственных колебаний ионов (b - коэффициент затухания).
Из теории дифференциальных уравнений известно, что решение полученного неоднородного уравнения представляет собой сумму общего решения однородного уравнения (уравнения затухающих колебаний) и частного решения неоднородного уравнения. Хорошо известно, что общее решение однородного уравнения (уравнения затухающих колебаний) не зависит от частоты вынуждающей силы (частоты света), поэтому не влияет на дисперсионные свойства вещества. Частное решение неоднородного уравнения ищут в виде, аналогичном виду вынуждающей силы, которое после подстановки в исходное дифференциальное уравнение, можно записать следующим образом
или
Относительное смещение ионов создает дипольный момент (поляризованность) элементарной ячейки кристаллической решетки. Поляризованность
полупроводника в целом равна сумме дипольных моментов
элементарных ячеек кристаллической решетки:
, где N - концентрация элементарных ячеек кристаллической решетки в единице объема. Тогда, учитывая связь диэлектрической проницаемости e полупроводника с дипольным моментом элементарной ячейки кристаллической решетки, получим зависимость диэлектрической проницаемости от частоты света:
Проанализировав полученную зависимость e(ω), можно заметить, что вдали от частоты собственных колебаний ωT ионов кристаллической решетки диэлектрическая проницаемость стремится к константам. Если диэлектрическую проницаемость при w → ∞ обозначить через ε∞, а при w → 0 через εs, то можно представить дисперсионную (частотную) зависимость диэлектрической проницаемости в виде:
а дисперсионную зависимость вещественной e1 и мнимой e2 составляющих диэлектрической проницаемости соответственно:
С другой стороны, диэлектрическая проницаемость связана с показателем преломления n и показателем поглощения k вещества:
e = (n - ik)2 = (n2 - k2) - i×2nk.
Из сравнения двух комплексных представлений диэлектрической проницаемости e (e1 - i e2 = n2 - k 2 - i×2nk) легко получить связь n и k с вещественной e1 и мнимой e2 составляющими диэлектрической проницаемости:
e1 = n2 - k2, e2 = 2nk.
Эта связь, с одной стороны, раскрывает характер зависимости показателя преломления и показателя поглощения от частоты падающего света n(w), k(w), а с другой, создает основу для экспериментальных исследований дисперсионных свойств полупроводников, так как n и k могут быть измерены экспериментально.
Методика изучения дисперсионных свойств полупроводника
Для определения параметров, входящих в аналитические выражения для e1(w) и e2(w), исследуем их на экстремум. Приравняв производную de1/dw к нулю, получим, что функция e1(w) имеет максимум и минимум соответственно при частотах
и
Аналогично можно показать, что мнимая составляющая e2 диэлектрической проницаемости имеет максимум при частоте
Максимум наблюдается при частоте близкой, но отличной от частоты собственной колебаний wT ионов. Это отличие тем больше, чем больше коэффициент затухания b (при b ® 0 w2 ® wT). Функции e1(w) и e2(w) немонотонны, значит связанные с ними функции n(w) и k(w) тоже немонотонны, и существует область вблизи частоты собственных колебаний wT ионов, где dn/dw < 0 - область аномальной дисперсии.
Вдали от частоты собственных колебаний ионов (|wT2 - w2| >> bw) мнимая составляющая e2 диэлектрической проницаемости, а значит и показатель поглощения k, обращается в ноль. Поэтому для получения параметров e¥, es, wT и b, определяющих зависимость диэлектрической проницаемости от длины волны (частоты) падающего света, необходимо экспериментально измерить показатель преломления света ns при w << wT, n¥ при w >> wT и определить, при каких частотах wmax и wmin наблюдаются максимальное и минимальное значения вещественной составляющей e1 диэлектрической проницаемости. Тогда получим:
Так как связь между парами величин n, k и e1, e2 нелинейная, максимумы и минимумы зависимостей n(w) и k(w) смещены относительно максимумов и минимумов зависимостей e1(w) и e2(w), и не совпадают ни с wT, ни с w2, wmax и wmin. Поэтому нельзя получить wT и b непосредственно из n(w), k(w), и приходится выполнять восстановление зависимостей e1(w), e2(w) как промежуточный этап исследования.
Используя параметры (e¥, es, wT) дисперсионной кривой, полученные по экспериментальным данным, можно рассчитать заряд q ионов кристаллической решетки полупроводника. Поскольку
тогда
где Ω = 1/N = a3/4 - объем элементарной ячейки, a - период кристаллической решетки.
Экспериментальное измерение показателя преломления и показателя поглощения
Для экспериментального определения показателя преломления и показателя поглощения рассмотрим прямоугольную призму, сечением которой является полукруг (полуцилиндр), изготовленную из изу- чаемого полупроводникового материала (рис. 1). Тогда показатель преломления n можно определить из закона преломления света n = sin a/sin g, а показатель поглощения - из закона Бугера-Ламберта-Бера
где коэффициент поглощения xλ связан с показателем поглощения k следующим соотношением: xλ = 4πk/λ.
Так как показатель поглощения существенно меняется при удалении от частоты поглощения, близкой к частоте собственных колебаний ионов кристаллической решетки, то относительная интенсивность I0/I света, прошедшего через призму (см. рис. 1), может меняться на несколько порядков. Поэтому для практических измерений более удобной является не относительная интенсивность I0/I, а оптическая плотность Dλ - безразмерная физическая величина, характеризующая степень непрозрачности слоя вещества толщиной R для световых лучей, численно равная десятичному логарифму относительной интенсивности света, прошедшего через поглощающую среду толщиной R:
Тогда измерив оптическую плотность Dl призмы при прохождении через нее света фиксированной длины волны l, можно определить показатель поглощения:
Рис. 1. Прохождение луча монохроматического света через призму
Таким образом, изложенный выше материал показывает, что для алмазоподобных полупроводников можно предложить достаточно простую, адекватную курсу общей физики физическую модель и методику экспериментального исследования аномальной дисперсии. Ее реализация в лабораторном практикуме позволит транслировать общность теоретического освещения данного вопроса на лабораторные исследования.
Компьютерная модель экспериментальной установки и результаты исследования
На основании изложенной физической модели и методики исследования авторами разработана компьютерная лабораторная работа «Аномальная дисперсия света в полупроводниках» (рис. 2).
Как показано выше, для определения параметров дисперсионной зависимости диэлектрической проницаемости необходимо выполнить измерения, как в области аномальной дисперсии, так и вдали от нее. Поэтому в работе используется источник инфракрасного монохроматического света, длину волны которого можно изменять в пределах 5-245 мкм. Вначале исследования необходимо определить, в каком диапазоне длин волн для изучаемого полупроводника наблюдается аномальная дисперсия. Для этого, установив неподвижно источник света так, чтобы луч падал под углом к плоской грани призмы, плавно увеличивают длину волны. В области нормальной дисперсии (dn/dλ < 0) угол преломления медленно увеличивается с увеличением длины волны падающего света. В области аномальной дисперсии (dn/dλ > 0), наоборот, - резко уменьшается. Таким образом, наблюдая за преломлением луча, можно определить, при какой длине волны угол преломления имеет минимальное и максимальное значения. Эти значения длин волн соответствуют диапазону, где наблюдается аномальная дисперсия для данного полупроводникового материала.
За счет нелинейного характера взаимосвязи диэлектрической проницаемости с показателем поглощения и показателем преломления максимум и минимум вещественной составляющей диэлектрической проницаемости e1 смещены относительно максимума и минимума показателя преломления n в сторону больших длин волн (рис. 3). Поэтому для корректного определения положения экстремумов зависимости вещественной составляющей диэлектрической проницаемости от длины волны необходимо выполнить ряд измерений величин n и k в области несколько более широкой, чем область аномальной дисперсии. Рассчитав по полученным для каждой длины волны данным значения вещественной и мнимой составляющих диэлектрической проницаемости, можно определить, при каких длинах волн (частотах wmax, wmin, w2) эти функции имеют максимум и минимум, на основании чего - рассчитать частоту wT собственных колебаний ионов и коэффициент затухания b этих колебаний в полупроводнике.
Для вычисления заряда ионов кроме wT необходимо знать показатель преломления ns и n¥ для длин волн, максимально удаленных от области аномальной дисперсии (для w << wT и w >> wT). Поэтому описанные выше измерения и расчеты дополняют измерениями показателя преломления и показателя поглощения для длин волн, равных 5 и 245 мкм. Рассчитанные на основании предложенной методики значения заряда ионов кристаллической решетки алмазоподобных полупроводников отличаются от значений, приведенных в научных публикациях [4], не более чем на 3-5%.
Рис. 2. Компьютерная лабораторная работа «Аномальная дисперсия света в полупроводниках»
Рис. 3. Экспериментальные дисперсионные зависимости
Компьютерная лабораторная работа «Аномальная дисперсия света в полупроводниках» позволяет выполнять исследования для следующих полупроводниковых материалов, обладающих различными дисперсионными свойствами: GaP, AlAs, InP, GaAs, GaSb, InAs, InSb. На рис. 3 приведены полученные при выполнении данной работы зависимости показателя преломления и показателя поглощения (рис. 3,а), а также вещественной и мнимой составляющих диэлектрической проницаемости (рис. 3,б) от длины волны падающего света для арсенида алюминия AlAs. По экспериментальным точкам, приведенным на рис. 3, видно, что для экстремумов вещественной и мнимой составляющих диэлектрической проницаемости наблюдается длинноволновое смещение по сравнению с экстремумами показателя преломления и показателя поглощения. Это общее свойство всех изучаемых в работе полупроводников.
Заключение
Теория, на которой базируется данная работа, иллюстрирует влияние взаимодействия электрического поля электромагнитной волны с ионами кристаллической решетки на дисперсионные свойства вещества, не рассматриваемое в курсе общей физики. Выполнение данной работы в рамках лабораторного практикума по общей физике позволяет студентам закрепить теоретическую общность классического объяснения нормальной и аномальной дисперсии света, расширить свои представления о важности различных этапов экспериментальных исследований, таких, например, как выделение в имеющемся диапазоне длин волн области аномальной дисперсии, или получение частотной зависимости непосредственно экспериментально неизмеряемой величины (диэлектрической проницаемости).
В методическом плане данную работу можно использовать как компьютерную модель - методическую основу для создания аналогичной натурной лабораторной работы. Разработка компьютерных лабораторных работ подобного типа позволяет не только выбрать физическую модель изучаемого явления, соответствующую курсу общей физики, но и апробировать в учебном процессе методику ее лабораторного исследования, доступную студентам в учебной лаборатории. Использование такого педагогического опыта существенно упрощает процесс постановки новых натурных работ лабораторного практикума. Это способствует переходу модернизации лабораторного практикума на качественно новый научно-методический и технический уровень [3].
Данная работа входит в комплекс лабораторных работ по изучению моделей физических процессов и явлений на компьютере [5], который разрабатывается на кафедре теоретической и экспериментальной физики Томского политехнического университета с 2002 г. В настоящее время комплекс включает 26 лабораторных работ.
Рецензенты:
-
Кистенев Ю.В., д.ф.-м.н., профессор кафедры теоретической и экспериментальной физики Томского государственного университета (ТГУ), г. Томск;
-
Кунашенко Ю.П., д.ф.-м.н., профессор кафедры теоретической физики Томского государственного педагогического университета (ТГПУ), г. Томск.
Работа поступила в редакцию 30.04.2012.