k · p теория возмущений - k·p perturbation theory

редактировать
Физическая модель твердого тела

В физике твердого тела, k · p теория возмущений представляет собой приближенный полуэмпирический подход для расчета зонной структуры (в частности, эффективной массы ) и оптических свойств кристаллических твердых тел. Он произносится как «k dot p», и его также называют «методом k · p ». Эта теория была применена специально в рамках модели Латтингера – Кона (после Хоакина Маздака Латтинджера и Уолтера Кона ) и модели Кейна модель (после Эван О. Кейн ).

Содержание

  • 1 Предпосылки и вывод
    • 1.1 Теорема Блоха и волновые векторы
    • 1.2 Теория возмущений
    • 1.3 Выражение для невырожденной полосы
      • 1.3.1 Эффективная масса
    • Модель 1,4 k · p со спин-орбитальным взаимодействием
    • 1.5 Расчет в вырожденном случае
  • 2 См. также
  • 3 Примечания и ссылки

Предпосылки и вывод

Теорема Блоха и волновые векторы

Согласно квантовая механикаодноэлектронном приближении ) квазисвободные электроны в любом твердом теле характеризуются волновыми функциями, которые являются собственными состояниями следующего стационарного уравнения Шредингера :

(p 2 2 m + V) ψ = E ψ {\ displaystyle \ left ({\ frac {p ^ {2}} {2m}} + V \ right) \ psi = E \ psi}\ left ({\ frac {p ^ {2}} {2m}} + V \ right) \ psi = E \ psi

где p - квантово-механический оператор импульса, V - потенциал, а m - вакуумная масса электрон. (Это уравнение не учитывает спин-орбитальный эффект ; см. Ниже.)

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

ψ n, k (x) = eik ⋅ xun, k (x) {\ displaystyle \ psi _ {n, \ mathbf {k}} (\ mathbf {x}) = e ^ {i \ mathbf {k} \ cdot \ mathbf {x}} u_ {n, \ mathbf {k}} (\ mathbf {x})}\ psi _ {{n, {\ mathbf {k}}}} ({\ mathbf {x}}) = e ^ {{я {\ mathbf {k}} \ cdot {\ mathbf {x}}}} u _ {{n, {\ mathbf {k}}}} ({\ mathbf {x}})

, где k - вектор (называемый волновым вектором), n - дискретный индекс (называемый индексом band ), а u n,k- функция с той же периодичностью, что и кристаллическая решетка.

Для любого заданного n связанные состояния называются полосой. В каждой полосе будет связь между волновым вектором k и энергией состояния E n,k, называемая полосной дисперсией . Вычисление этой дисперсии - одно из основных приложений теории возмущений k · p.

Теория возмущений

Периодическая функция u n,kудовлетворяет следующему уравнению типа Шредингера (просто прямому разложению уравнения Шредингера с волновой функцией типа Блоха):

ЧАС Кун, К = Е N, Кун, К {\ Displaystyle H _ {\ mathbf {k}} u_ {n, \ mathbf {k}} = E_ {n, \ mathbf {k}} u_ {n, \ mathbf { k}}}H _ {{{\ mathbf {k}}}} u _ {{n, {\ mathbf {k}}}} = E_ {{n, {\ mathbf {k}}}} u _ {{n, {\ mathbf {k}}}}

где гамильтониан равен

H k = p 2 2 m + ℏ k ⋅ pm + ℏ 2 k 2 2 m + V {\ displaystyle H _ {\ mathbf {k }} = {\ frac {p ^ {2}} {2m}} + {\ frac {\ hbar \ mathbf {k} \ cdot \ mathbf {p}} {m}} + {\ frac {\ hbar ^ { 2} k ^ {2}} {2m}} + V}H _ {{{\ mathbf {k} }}} = {\ frac {p ^ {2}} {2m}} + {\ frac {\ hbar {\ mathbf {k}} \ cdot {\ mathbf {p}}} {m}} + {\ frac {\ hbar ^ {2} k ^ {2}} {2m}} + V

Обратите внимание, что k - это вектор, состоящий из трех действительных чисел с размерностями обратной длины, а p - вектор операторов; чтобы быть явным,

k ⋅ p = kx (- i ℏ ∂ ∂ x) + ky (- i ℏ ∂ ∂ y) + kz (- i ℏ ∂ ∂ z) {\ displaystyle \ mathbf {k} \ cdot \ mathbf {p} = k_ {x} (- i \ hbar {\ frac {\ partial} {\ partial x}}) + k_ {y} (- i \ hbar {\ frac {\ partial} {\ partial y }}) + k_ {z} (- i \ hbar {\ frac {\ partial} {\ partial z}})}{\ mathbf {k}} \ cdot {\ mathbf {p}} = k_ {x} (- i \ hbar {\ frac {\ partial} {\ partial x}}) + k_ {y} (- i \ hbar {\ frac {\ partial} {\ partial y}}) + k_ {z} (- i \ hbar {\ frac {\ partial} {\ partial z}})

В любом случае мы запишем этот гамильтониан как сумму двух членов:

H Знак равно ЧАС 0 + ЧАС К ', ЧАС 0 знак равно п 2 2 м + В, ЧАС К' = ℏ 2 К 2 2 м + ℏ К ⋅ pm {\ Displaystyle Н = Н_ {0} + Н _ {\ mathbf {k} } ', \; \; H_ {0} = {\ frac {p ^ {2}} {2m}} + V, \; \; H _ {\ mathbf {k}}' = {\ frac {\ hbar ^ {2} k ^ {2}} {2m}} + {\ frac {\ hbar \ mathbf {k} \ cdot \ mathbf {p}} {m}}}H=H_{0}+H_{{{\mathbf {k}}}}',\;\;H_{0}={\frac {p^{2}}{2m}}+V,\;\;H_{{{\mathbf {k}}}}'={\frac {\hbar ^{2}k^{2}}{2m}}+{\frac {\hbar {\mathbf {k}}\cdot {\mathbf {p}}}{m}}

Это выражение является основой для теория возмущений. «Невозмущенный гамильтониан» равен H 0, что фактически равно точному гамильтониану при k = 0 (то есть в гамма-точке ). «Возмущение» - это термин H k ′ {\ displaystyle H _ {\ mathbf {k}} '}H_{{{\mathbf {k}}}}'. Полученный в результате анализ называется «k · p теорией возмущений» из-за члена, пропорционального k · p . Результатом этого анализа является выражение для E n,kи u n,kв терминах энергий и волновых функций при k = 0.

Обратите внимание, что термин «возмущение» H k ′ {\ displaystyle H _ {\ mathbf {k}} '}H_{{{\mathbf {k}}}}'становится все меньше по мере того, как k приближается к нулю. Следовательно, теория возмущений k · p является наиболее точной для малых значений k . Однако, если в пертурбативное разложение включено достаточно членов, то теория может быть достаточно точной для любого значения k во всей зоне Бриллюэна.

Выражение для невырожденной полосы

Для невырожденной полосы (т. е. полосы, которая имеет энергию, отличную от энергии при k = 0 от любой другой полосы), с экстремумом при k = 0, и без спин-орбитальной связи результат теории возмущений k · p (до низшего нетривиального порядка ):

ООН, К знак равно ООН, 0 + ℏ м ∑ N '≠ N ⟨UN, 0 | К ⋅ p | un', 0 ′ E n, 0 - E n ', 0 un', 0 {\ displaystyle u_ {n, \ mathbf {k}} = u_ {n, 0} + {\ frac {\ hbar} {m}} \ sum _ {n '\ neq n} {\ frac {\ langle u_ {n, 0} | \ mathbf {k} \ cdot \ mathbf {p} | u_ {n ', 0} \ rangle} {E_ {n, 0} -E_ {n', 0}}} u_ {n ', 0}}u_{{n,{\mathbf {k}}}}=u_{{n,0}}+{\frac {\hbar }{m}}\sum _{{n'\neq n}}{\frac {\langle u_{{n,0}}|{\mathbf {k}}\cdot {\mathbf {p}}|u_{{n',0}}\rangle }{E_{{n,0}}-E_{{n',0}}}}u_{{n',0}}
E n, k = E n, 0 + ℏ 2 k 2 2 m + ℏ 2 m 2 ∑ n ′ ≠ n | ⟨un, 0 | k ⋅ p | un ′, 0⟩ | 2 E n, 0 - E n ′, 0 {\ displaystyle E_ {n, \ mathbf {k}} = E_ {n, 0} + {\ frac {\ hbar ^ {2} k ^ {2}} {2m}} + {\ frac {\ hbar ^ {2}} {м ^ {2}} } \ sum _ {n '\ neq n} {\ frac {| \ langle u_ {n, 0} | \ mathbf {k} \ cdot \ mathbf {p} | u_ {n', 0} \ rangle | ^ { 2}} {E_ {n, 0} -E_ {n ', 0}}}}E_{{n,{\mathbf {k}}}}=E_{{n,0}}+{\frac {\hbar ^{2}k^{2}}{2m}}+{\frac {\hbar ^{2}}{m^{2}}}\sum _{{n'\neq n}}{\frac {|\langle u_{{n,0}}|{\mathbf {k}}\cdot {\mathbf {p}}|u_{{n',0}}\rangle |^{2}}{E_{{n,0}}-E_{{n',0}}}}

Поскольку k является вектором действительных чисел (а не вектором более сложных линейных операторов), матричный элемент в этих выражениях можно переписать как:

⟨un, 0 | k ⋅ p | u n ′, 0⟩ = k ⋅ ⟨u n, 0 | p | un ', 0⟩ {\ displaystyle \ langle u_ {n, 0} | \ mathbf {k} \ cdot \ mathbf {p} | u_ {n', 0} \ rangle = \ mathbf {k} \ cdot \ langle u_ {n, 0} | \ mathbf {p} | u_ {n ', 0} \ rangle}\langle u_{{n,0}}|{\mathbf {k}}\cdot {\mathbf {p}}|u_{{n',0}}\rangle ={\mathbf {k}}\cdot \langle u_{{n,0}}|{\mathbf {p}}|u_{{n',0}}\rangle

Следовательно, можно вычислить энергию при любом k, используя только несколько неизвестных параметров, а именно E n, 0 и ⟨un, 0 | p | u n ′, 0⟩ {\ displaystyle \ langle u_ {n, 0} | \ mathbf {p} | u_ {n ', 0} \ rangle}\langle u_{{n,0}}|{\mathbf {p}}|u_{{n',0}}\rangle . Последние называются «оптическими матричными элементами» и тесно связаны с дипольными моментами перехода . Эти параметры обычно выводятся из экспериментальных данных.

На практике сумма по n часто включает только ближайшую одну или две полосы, поскольку они, как правило, являются наиболее важными (из-за знаменателя). Однако для повышения точности, особенно при больших k, необходимо включить больше полос, а также больше членов в пертурбативное разложение, чем написано выше.

Эффективная масса

Используя приведенное выше выражение для зависимости дисперсии энергии, можно найти упрощенное выражение для эффективной массы в зоне проводимости полупроводника. Чтобы аппроксимировать дисперсионное соотношение в случае зоны проводимости, возьмите энергию E n0 как минимальную энергию зоны проводимости E c0 и включите в суммирование только члены с энергиями, близкими к валентной максимум полосы, где разница энергий в знаменателе наименьшая. (Эти члены вносят наибольший вклад в суммирование.) Затем этот знаменатель аппроксимируется как ширина запрещенной зоны E g, что приводит к выражению энергии:

E c (k) ≈ E c 0 + ( ℏ k) 2 2 m + ℏ 2 E gm 2 n | ⟨U c, 0 | k ⋅ p | u n, 0⟩ | 2 {\ displaystyle E_ {c} ({\ boldsymbol {k}}) \ приблизительно E_ {c0} + {\ frac {(\ hbar k) ^ {2}} {2m}} + {\ frac {\ hbar ^ {2}} {{E_ {g}} m ^ {2}}} \ sum _ {n} {| \ langle u_ {c, 0} | \ mathbf {k} \ cdot \ mathbf {p} | u_ { n, 0} \ rangle | ^ {2}}}E_ {c} ({\ boldsymbol k}) \ приблизительно E _ {{c0}} + {\ frac {(\ hbar k) ^ {2}} {2m}} + {\ frac {\ hbar ^ {2 }} {{E_ {g}} m ^ {2}}} \ sum _ {n} {| \ langle u _ {{c, 0}} | {\ mathbf {k}} \ cdot {\ mathbf {p} } | u _ {{n, 0}} \ rangle | ^ {2}}

Тогда эффективная масса в направлении ℓ равна:

1 m ℓ = 1 ℏ 2 ∑ m ⋅ ∂ 2 E c (k) ∂ k ℓ ∂ km ≈ 1 м + 2 E gm 2 ∑ m, n ⟨uc, 0 | p ℓ | u n, 0⟩ ⟨u n, 0 | п м | uc, 0⟩ {\ displaystyle {\ frac {1} {{m} _ {\ ell}}} = {{1} \ over {\ hbar ^ {2}}} \ sum _ {m} \ cdot {{ \ partial ^ {2} E_ {c} ({\ boldsymbol {k}})} \ over {\ partial k _ {\ ell} \ partial k_ {m}}} \ приблизительно {\ frac {1} {m}} + {\ frac {2} {E_ {g} m ^ {2}}} \ sum _ {m, \ n} {\ langle u_ {c, 0} | p _ {\ ell} | u_ {n, 0} \ rangle} {\ langle u_ {n, 0} | p_ {m} | u_ {c, 0} \ rangle}}{\ displaystyle {\ frac {1} {{m} _ {\ ell}}} = {{1} \ over {\ hbar ^ {2}}} \ sum _ {m} \ cdot {{\ partial ^ {2} E_ {c} ({\ boldsymbol {k}})} \ over {\ partial k _ {\ ell} \ partial k_ {m}}} \ приблизительно {\ frac {1} {m}} + {\ frac {2} {E_ {g} m ^ {2}}} \ sum _ {m, \ n} {\ langle u_ {c, 0} | p _ {\ ell } | u_ {n, 0} \ rangle} {\ langle u_ {n, 0} | p_ {m} | u_ {c, 0} \ rangle}}

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

2 E g m 2 ∑ m, n | ⟨U c, 0 | p ℓ | u n, 0⟩ | | ⟨U c, 0 | п м | u n, 0⟩ | ≈ 20 е В 1 м E g, {\ displaystyle {\ frac {2} {E_ {g} m ^ {2}}} \ sum _ {m, \ n} {| \ langle u_ {c, 0} | p _ {\ ell} | u_ {n, 0} \ rangle |} {| \ langle u_ {c, 0} | p_ {m} | u_ {n, 0} \ rangle |} \ приблизительно 20 \ mathrm {eV} {\ frac {1} {mE_ {g}}} \,}{\ frac {2} {E_ {g} m ^ {2} }} \ sum _ {{m, \ n}} {| \ langle u _ {{c, 0}} | p _ {{\ ell}} | u _ {{n, 0}} \ rangle |} {| \ langle u _ {{c, 0}} | p _ {{m}} | u _ {{n, 0}} \ rangle |} \ около 20 {\ mathrm {eV}} {\ frac {1} {mE _ {{g} }}} \,

, который применяется в пределах 15% или лучше к большинству полупроводников групп IV, III-V и II-VI.

В отличие от этого В этом простом приближении в случае энергии валентной зоны необходимо ввести спин-орбитальное взаимодействие (см. ниже), и нужно отдельно рассматривать еще много зон. Расчет представлен в Yu и Cardona. В валентной зоне мобильные носители - это дырки. Было обнаружено, что есть два типа дыр, названные тяжелыми и легкими, с анизотропными массами.

модель k · p со спин-орбитальным взаимодействием

Включая спин-орбитальное взаимодействие, уравнение Шредингера для u:

H kun, k = E п, кун, к {\ displaystyle H _ {\ mathbf {k}} u_ {n, \ mathbf {k}} = E_ {n, \ mathbf {k}} u_ {n, \ mathbf {k}}}H _ {{{\ mathbf {k}}}} u _ {{n, {\ mathbf {k}}}} = E_ {{n, {\ mathbf {k}}}} u _ {{n, {\ mathbf {k}}}}

где

ЧАС К знак равно п 2 2 м + ℏ мк ⋅ п + ℏ 2 к 2 2 м + V + ℏ 4 м 2 с 2 (∇ V × (п + ℏ к)) ⋅ σ → {\ displaystyle H _ {\ mathbf {k}} = {\ frac {p ^ {2}} {2m}} + {\ frac {\ hbar} {m}} \ mathbf {k} \ cdot \ mathbf {p} + {\ frac {\ hbar ^ {2} k ^ {2}} {2m}} + V + {\ frac {\ hbar} {4m ^ {2} c ^ {2}}} (\ nabla V \ times (\ mathbf { p} + \ hbar \ mathbf {k})) \ cdot {\ vec {\ sigma}}}H _ {{{\ mathbf {k}}}} = {\ frac {p ^ {2}} {2m}} + { \ frac {\ hbar} {m}} {\ mathbf {k}} \ cdot {\ mathbf {p}} + {\ frac {\ hbar ^ {2} k ^ {2}} {2m}} + V + {\ frac {\ hbar} {4m ^ {2} c ^ {2}}} (\ nabla V \ times ({\ mathbf {p}} + \ hbar {\ mathbf {k}})) \ cdot {\ vec \ sigma}

где σ → = (σ x, σ y, σ z) {\ displaystyle {\ vec { \ sigma}} = (\ sigma _ {x}, \ sigma _ {y}, \ sigma _ {z})}{\ vec \ sigma} = (\ sigma _ {x}, \ sigma _ {y}, \ sigma _ {z}) - вектор, состоящий из трех матриц Паули. Этот гамильтониан можно подвергнуть тому же анализу теории возмущений, что и выше.

Расчет в вырожденном случае

Для вырожденных или почти вырожденных полос, в частности, валентных зон в некоторых материалах, таких как арсенид галлия, уравнения могут быть проанализированы методами теории вырожденных возмущений. К моделям этого типа относятся «модель Латтинджера – Кона » (также известная как «модель Кона – Латтинджера») и «модель Кейна ".

Как правило, эффективный гамильтониан H eff { \ displaystyle H ^ {\ rm {eff}}}{\ displaystyle H ^ {\ rm {eff}}} , и для первого порядка его матричные элементы могут быть выражены как

H k, mneff = ⟨um, 0 | H 0 | un, 0⟩ + К ⋅ ⟨um, 0 | ∇ К ЧАС к '| un, 0⟩ {\ Displaystyle H _ {\ mathbf {k}, mn} ^ {\ rm {eff}} = \ langle u_ {m, 0} | H_ {0} | u_ {n, 0} \ rangle + \ mathbf {k} \ cdot \ langle u_ {m, 0} | \ nabla _ {\ mathbf {k}} H _ {\ mathbf {k} } '| u_ {n, 0} \ rangle}{\displaystyle H_{\mathbf {k},mn}^{\rm {eff}}=\langle u_{m,0}|H_{0}|u_{n,0}\rangle +\mathbf {k} \cdot \langle u_{m,0}|\nabla _{\mathbf {k} }H_{\mathbf {k} }'|u_{n,0}\rangle }

После решения получаются волновые функции и энергетические зоны.

См. также

Электронная зонная структура

Свойства зон

Волновые функции

Фундаментальная теория

N Примечания и ссылки

Последняя правка сделана 2021-05-26 04:30:17
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте