نظرية الاضطراب k·p

في فيزياء الحالة الصلبة ، تُعدّ نظرية الاضطراب k·p منهجًا شبه تجريبي تقريبيًا لحساب بنية النطاقات (وخاصة الكتلة الفعّالة ) والخصائص البصرية للمواد الصلبة البلورية. [ 1 ] [ 2 ] [ 3 ] تُنطق "k dot p"، وتُعرف أيضًا باسم طريقة k·p . وقد طُبّقت هذه النظرية تحديدًا في إطار نموذج لوتينجر-كون (نسبةً إلى خواكين مازداك لوتينجر ووالتر كون )، ونموذج كين (نسبةً إلى إيفان أو. كين ).

الخلفية والاشتقاق

نظرية بلوخ ومتجهات الموجة

وفقًا لميكانيكا الكم (في تقريب الإلكترون الواحد )، تتميز الإلكترونات شبه الحرة في أي مادة صلبة بدوال موجية تمثل حالات ذاتية لمعادلة شرودنغر الثابتة التالية :

(ص22م+V)ψ=هـψ{\displaystyle \left({\frac {p^{2}}{2m}}+V\right)\psi =E\psi }

حيث p هو مؤثر الزخم الكمومي ، وV هو الجهد ، و m هي كتلة الإلكترون في الفراغ. (تتجاهل هذه المعادلة تأثير الدوران المداري ؛ انظر أدناه.)

في المواد الصلبة البلورية ، تكون V دالة دورية ، ولها نفس دورية الشبكة البلورية . تثبت نظرية بلوخ أن حلول هذه المعادلة التفاضلية يمكن كتابتها على النحو التالي:

ψن،ك(x)=هـأناكxuن،ك(x){\displaystyle \psi _{n,\mathbf {k} }(\mathbf {x} )=e^{i\mathbf {k} \cdot \mathbf {x} }u_{n,\mathbf {k} }(\mathbf {x} )}

حيث k هو متجه (يسمى متجه الموجةو n هو مؤشر منفصل (يسمى مؤشر النطاق )، و u n ، k هي دالة لها نفس الدورية مثل الشبكة البلورية.

لكل قيمة معطاة لـ n ، تُسمى الحالات المرتبطة بها نطاقًا . في كل نطاق، توجد علاقة بين متجه الموجة k وطاقة الحالة E <sub>n , k </sub>، تُسمى تشتت النطاق . يُعد حساب هذا التشتت أحد التطبيقات الأساسية لنظرية الاضطراب k · p .

نظرية الاضطراب

الدالة الدورية u n , k تحقق معادلة شرودنغر التالية (ببساطة، توسيع مباشر لمعادلة شرودنغر بدالة موجية من نوع بلوخ): [ 1 ]

حكuن،ك=هـن،كuن،ك{\displaystyle H_{\mathbf {k}}u_{n,\mathbf {k} }=E_{n,\mathbf {k} }u_{n,\mathbf {k} }}

حيث يكون الهاميلتوني

حك=ص22م+كصم+2ك22م+V{\displaystyle H_{\mathbf {k} }={\frac {p^{2}}{2m}}+{\frac {\hbar \mathbf {k} \cdot \mathbf {p} }{m}}+{\frac {\hbar ^{2}k^{2}}{2m}}+V}

لاحظ أن k متجه يتكون من ثلاثة أعداد حقيقية بأبعاد تساوي طولها العكسي ، بينما p متجه من المؤثرات؛ وللتوضيح،

كص=كx(-أناx)+كy(-أناy)+كz(-أنا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}})}

على أي حال، نكتب هذه الدالة الهاميلتونية كمجموع حدين:

حك=ح0+حك،ح0=ص22م+V،حك=2ك22م+كصم// {\hbar \mathbf {k} \cdot \mathbf {p} }{m}}}

هذا التعبير هو أساس نظرية الاضطراب . الهاميلتوني "غير المضطرب" هو H₀ ، والذي يساوي في الواقع الهاميلتوني الدقيق عند k  =  0 (أي عند النقطة غاما ). "الاضطراب" هو الحدحك{\displaystyle H_{\mathbf {k} }'}يُطلق على التحليل الناتج اسم نظرية الاضطراب k·p ، نسبةً إلى الحد المتناسب مع k·p . وتتمثل نتيجة هذا التحليل في تعبير عن E <sub>n , k</sub> و u <sub>n , k </sub> بدلالة الطاقات والدوال الموجية عند k  =  0.

لاحظ أن مصطلح "الاضطراب"حك{\displaystyle H_{\mathbf {k} }'}يتناقص حجم k تدريجيًا كلما اقتربت قيمة k من الصفر. لذا، تكون نظرية الاضطراب k·p أكثر دقةً عند القيم الصغيرة لـ k . مع ذلك، إذا أُضيفت حدود كافية إلى متسلسلة الاضطراب ، فإن النظرية قد تكون دقيقةً إلى حد معقول لأي قيمة لـ k في منطقة بريلوين بأكملها .

تعبير عن نطاق غير متدهور

بالنسبة لحزمة غير متدهورة (أي حزمة لها طاقة مختلفة عند k  =  0 عن أي حزمة أخرى)، مع قيمة قصوى عند k  =  0، وبدون اقتران الدوران المداري ، فإن نتيجة نظرية الاضطراب k · p هي (إلى أدنى رتبة غير تافهة ): [ 1 ]

uن،ك=uن،0+مننuن،0|كص|uن،0هـن،0-هـن،0uن،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}}
هـن،ك=هـن،0+2ك22م+2م2نن|uن،0|كص|uن،0|2هـن،0-هـن،0// {ص} |u_{n',0}\rangle |^{2}}{E_{n,0}-E_{n',0}}}}

بما أن k عبارة عن متجه من الأعداد الحقيقية (بدلاً من متجه من المؤثرات الخطية الأكثر تعقيدًا)، يمكن إعادة كتابة عنصر المصفوفة في هذه التعبيرات على النحو التالي:

uن،0|كص|uن،0=كuن،0|ص|uن،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 }

لذلك، يمكن حساب الطاقة عند أي قيمة k باستخدام عدد قليل من المعاملات المجهولة، وهي E n ,0 وuن،0|ص|uن،0{\displaystyle \langle u_{n,0}|\mathbf {p} |u_{n',0}\rangle }تُسمى هذه الأخيرة "عناصر المصفوفة البصرية"، وهي مرتبطة ارتباطًا وثيقًا بعزوم ثنائي القطب الانتقالي . وعادةً ما يتم استنتاج هذه المعلمات من البيانات التجريبية.

عمليًا، غالبًا ما يشمل المجموع على n نطاقًا واحدًا أو نطاقين فقط، نظرًا لأهميتهما (بسبب المقام). مع ذلك، لتحسين الدقة، خاصةً عند قيم k الكبيرة ، يجب تضمين نطاقات أكثر، بالإضافة إلى حدود أكثر في التوسع الاضطرابي من تلك المذكورة أعلاه.

الكتلة الفعالة

باستخدام الصيغة المذكورة أعلاه لعلاقة تشتت الطاقة، يمكن إيجاد صيغة مبسطة للكتلة الفعالة في نطاق التوصيل لأشباه الموصلات. [ 3 ] لتقريب علاقة التشتت في حالة نطاق التوصيل، نعتبر الطاقة E <sub>n0 </sub> هي الحد الأدنى لطاقة نطاق التوصيل E <sub> c0 </sub>، ونُضمّن في المجموع فقط الحدود ذات الطاقات القريبة من الحد الأقصى لنطاق التكافؤ، حيث يكون فرق الطاقة في المقام أصغر ما يمكن. (هذه الحدود هي المساهمات الأكبر في المجموع). ثم يُقرّب هذا المقام بفجوة النطاق E<sub> g</sub> ، مما يؤدي إلى صيغة الطاقة التالية:

هـج(ك)هـج0+(ك)22م+2هـزم2ن|uج،0|كص|uن،0|2{\displaystyle E_{c}({\boldsymbol {k}})\approx 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م=12م2هـج(ك)ككم1م+2هـزم2م، نuج،0|ص|uن،0uن،0|صم|uج،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}}}\approx {\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 }}

بغض النظر عن تفاصيل عناصر المصفوفة، فإن النتائج الرئيسية هي أن الكتلة الفعالة تتغير بتغير أصغر فجوة نطاق وتؤول إلى الصفر عندما تؤول الفجوة إلى الصفر. [ 3 ] ويُعدّ التقريب التالي مفيدًا لعناصر المصفوفة في أشباه الموصلات ذات الفجوة المباشرة : [ 4 ]

2هـزم2م، ن|uج،0|ص|uن،0||uج،0|صم|uن،0|20هـV1مهـز ،{\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 |}\approx 20\mathrm {eV} {\frac {1}{mE_{g}}}\ ,}

وينطبق ذلك بنسبة 15% أو أفضل على معظم أشباه الموصلات من المجموعات الرابعة والثالثة-الخامسة والثانية-السادسة. [ 5 ]

على النقيض من هذا التقريب البسيط، في حالة طاقة نطاق التكافؤ، يجب إدخال تفاعل الدوران المداري (انظر أدناه)، كما يجب مراعاة العديد من النطاقات الأخرى بشكل فردي. يُقدَّم الحساب في بحث يو وكاردونا [ 6 ] . في نطاق التكافؤ، تكون حاملات الشحنة المتحركة عبارة عن فجوات . يوجد نوعان من الفجوات، تُسمى الفجوات الثقيلة والفجوات الخفيفة ، ذات كتل متباينة الخواص.

نموذج k·p مع تفاعل الدوران المداري

وبإضافة تفاعل الدوران المداري ، تكون معادلة شرودنغر لـ u هي: [ 2 ]

حكuن،ك=هـن،كuن،ك{\displaystyle H_{\mathbf {k} }u_{n,\mathbf {k} }=E_{n,\mathbf {k} }u_{n,\mathbf {k} }}

حيث [ 7 ]

حك=ص22م+مكص+2ك22م+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 }}}

أينσ=(σx،σy،σz){\displaystyle {\vec {\sigma }}=(\sigma _{x},\sigma _{y},\sigma _{z})}هو متجه يتكون من مصفوفات باولي الثلاث . يمكن إخضاع هذا الهاميلتوني لنفس نوع تحليل نظرية الاضطراب المذكور أعلاه.

الحساب في الحالة المنحلة

بالنسبة للنطاقات المتدهورة أو شبه المتدهورة، ولا سيما نطاقات التكافؤ في بعض المواد مثل زرنيخيد الغاليوم ، يمكن تحليل المعادلات باستخدام طرق نظرية الاضطراب المتدهور . [ 1 ] [ 2 ] تشمل نماذج هذا النوع " نموذج لوتينجر-كون " (المعروف أيضًا باسم "نموذج كون-لوتينجر")، [ 8 ] و" نموذج كين ". [ 7 ]

بشكل عام، هاميلتوني فعالحهـوو{\displaystyle H^{\rm {eff}}}يتم تقديمها، ويمكن التعبير عن عناصر مصفوفتها، من الدرجة الأولى، على النحو التالي:

حك،منهـوو=uم،0|ح0|uن،0+كuم،0|كحك|uن،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 }

بعد حلها، يتم الحصول على الدوال الموجية ونطاقات الطاقة.

انظر أيضاً

ملاحظات ومراجع

  1. 1 2 3 4 ب. يو، م. كاردونا (2005). أساسيات أشباه الموصلات: الفيزياء وخواص المواد ( الطبعة الثالثة). سبرينغر . القسم 2.6، الصفحات 68 وما بعدها . ISBN  3-540-25470-6.
  2. 1 2 3 سي. كيتل (1987). نظرية الكم للمواد الصلبة (الطبعة الثانية المنقحة ). نيويورك: وايلي . الصفحات 186-190 . ISBN   0-471-62412-8.
  3. 1 2 3 دبليو. بي. هاريسون (1989) [1980]. البنية الإلكترونية وخواص المواد الصلبة (طبعة مُعاد طباعتها ). منشورات دوفر . ص 158 وما بعدها . ISBN   0-486-66021-4.
  4. أشباه الموصلات ذات الفجوة المباشرة هي تلك التي يحدث فيها الحد الأقصى لنطاق التكافؤ والحد الأدنى لنطاق التوصيل في نفس الموضع في فضاء k ، وعادة ما تكون ما يسمى بنقطة Γ حيث k = 0.
  5. انظر الجدول 2.22 في Yu & Cardona، المرجع السابق.
  6. ^ انظر يو وكاردونا، مرجع سابق. سيتي. ص 75-82
  7. 1 2 إيفان أو. كين (1957). "بنية نطاقات إنديوم أنتيمونيد". مجلة فيزياء وكيمياء المواد الصلبة . 1 (4): 249-261 . Bibcode : 1957JPCS....1..249K . doi : 10.1016/0022-3697(57)90013-6 .
  8. ج. م. لوتينجر، و. كون (1955). "حركة الإلكترونات والفجوات في الحقول الدورية المضطربة". مجلة Physical Review . 97 (4): 869–883 . Bibcode : 1955PhRv...97..869L . doi : 10.1103/PhysRev.97.869 .