كثافة الولايات

في فيزياء المادة المكثفة ، تصف كثافة الحالات ( DOS ) لنظام ما عدد الأنماط أو الحالات المسموح بها لكل وحدة نطاق طاقة. وتُعرَّف كثافة الحالات على النحو التالي:د(هـ)=شمال(هـ)/V{\displaystyle D(E)=N(E)/V}، أينشمال(هـ)دلتاهـ{\displaystyle N(E)\delta E}يمثل عدد الحالات في نظام الحجمV{\displaystyle V}والتي تتراوح طاقاتها بينهـ{\displaystyle E}لهـ+دلتاهـ{\displaystyle E+\delta E}يُعبَّر عنها رياضيًا كتوزيع بواسطة دالة كثافة احتمالية ، وهي عمومًا متوسطٌ على نطاقي المكان والزمان لمختلف الحالات التي يشغلها النظام. ترتبط كثافة الحالات ارتباطًا مباشرًا بعلاقات التشتت لخصائص النظام. تشير كثافة الحالات العالية عند مستوى طاقة محدد إلى توفر العديد من الحالات للإشغال.

بشكل عام، تكون كثافة حالات المادة متصلة. أما في الأنظمة المعزولة ، مثل الذرات أو الجزيئات في الحالة الغازية، فإن توزيع الكثافة يكون منفصلاً ، مثل الكثافة الطيفية . وتُعرف التغيرات الموضعية، التي غالباً ما تنتج عن تشوهات في النظام الأصلي، باسم كثافات الحالات الموضعية (LDOSs).

مقدمة

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

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

يُحدد هذا ما إذا كانت المادة عازلة أم معدنية في بُعد الانتشار. كما أن عدد الحالات في نطاق الطاقة مفيدٌ للتنبؤ بخصائص التوصيل. على سبيل المثال، في بنية بلورية أحادية البعد، يؤدي وجود عدد فردي من الإلكترونات لكل ذرة إلى امتلاء النطاق العلوي بنصف طاقته؛ وتوجد إلكترونات حرة عند مستوى فيرمي، مما ينتج عنه مادة معدنية. من ناحية أخرى، يؤدي وجود عدد زوجي من الإلكترونات إلى ملء عدد صحيح من نطاقات الطاقة، تاركًا الباقي فارغًا. إذا كان مستوى فيرمي يقع في فجوة نطاق طاقة مشغولة بين أعلى حالة مشغولة وأدنى حالة فارغة، فإن المادة ستكون عازلة أو شبه موصلة .

بحسب النظام الكمومي، يمكن حساب كثافة الحالات للإلكترونات أو الفوتونات أو الفونونات ، ويمكن التعبير عنها كدالة للطاقة أو متجه الموجة k . وللتحويل بين كثافة الحالات كدالة للطاقة وكثافة الحالات كدالة لمتجه الموجة، يجب معرفة علاقة تشتت الطاقة الخاصة بالنظام بين E و k .

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

تعريف

يتم تعريف كثافة الحالات المتعلقة بالحجم V ومستويات الطاقة N القابلة للعد على النحو التالي:د(هـ)=1Vأنا=1شمالدلتا(هـ-هـ(كأنا)).{\displaystyle D(E)={\frac {1}{V}}\,\sum _{i=1}^{N}\delta (EE({\mathbf {k} }_{i})).} لأن أصغر تغيير مسموح به في الزخمك{\displaystyle k}بالنسبة لجزيء في صندوق ذي أبعادد{\displaystyle d}والطولل{\displaystyle L}يكون(Δك)د=(2π/ل)د{\displaystyle (\Delta k)^{d}=({2\pi }/{L})^{d}}، يتم الحصول على كثافة الحالات المتعلقة بالحجم لمستويات الطاقة المستمرة في الحدل{\displaystyle L\to \infty }مثل د(هـ):=Rدددك(2π)ددلتا(هـ-هـ(ك))،{\displaystyle D(E):=\int _{\mathbb {R} ^{d}}{\frac {\mathrm {d} ^{d}k}{(2\pi )^{d}}}\cdot \delta (EE(\mathbf {k} )),} هنا،د{\displaystyle d}يمثل البعد المكاني للنظام المدروس وك{\displaystyle \mathbf {k} }متجه الموجة.

بالنسبة للأنظمة أحادية البعد المتناحية ذات تشتت الطاقة المكافئ، فإن كثافة الحالات هيد1د(هـ)=12π(2مهـ)1/2{\textstyle D_{1D}(E)={\tfrac {1}{2\pi \hbar }}({\tfrac {2m}{E}})^{1/2}}في بعدين، تكون كثافة الحالات ثابتة .د2د=م2π2{\displaystyle D_{2D}={\tfrac {m}{2\pi \hbar ^{2}}}}بينما في ثلاثة أبعاد يصبحد3د(هـ)=م2π23(2مهـ)1/2{\displaystyle D_{3D}(E)={\tfrac {m}{2\pi ^{2}\hbar ^{3}}}(2mE)^{1/2}}.

وبالمثل، يمكن فهم كثافة الحالات أيضًا على أنها مشتقة دالة التقسيم الميكروكانونيةZم(هـ){\displaystyle Z_{m}(E)}(أي العدد الإجمالي للحالات ذات الطاقة الأقل منهـ{\displaystyle E}) فيما يتعلق بالطاقة: د(هـ)=1VدZم(هـ)دهـ.{\displaystyle D(E)={\frac {1}{V}}\cdot {\frac {\mathrm {d} Z_{m}(E)}{\mathrm {d} E}}.}

عدد الحالات التي تحتوي على طاقةهـ{\displaystyle E'}(درجة الانحطاط) تُعطى بالصيغة التالية: ز(هـ)=ليمΔهـ0هـهـ+Δهـد(هـ)دهـ=ليمΔهـ0د(هـ)Δهـ،{\displaystyle g\left(E'\right)=\lim _{\Delta E\to 0}\int _{E'}^{E'+\Delta E}D(E)\,\mathrm {d} E=\lim _{\Delta E\to 0}D\left(E'\right)\Delta E,} حيث لا تنطبق المساواة الأخيرة إلا عندما تكون نظرية القيمة المتوسطة للتكاملات صالحة.

التناظر

المنطقة الأولى من بريلوين لشبكة FCC ، وهي عبارة عن ثماني السطوح مقطوع ، تُظهر علامات التناظر لخطوط ونقاط التناظر العالي

هناك مجموعة كبيرة ومتنوعة من الأنظمة وأنواع الحالات التي يمكن إجراء حسابات كثافة الحالات لها.

تتمتع بعض أنظمة المادة المكثفة بتناظر بنيوي على المستوى المجهري، ما يُسهّل حساب كثافة حالاتها. في الأنظمة ذات التناظر الكروي، تكون تكاملات الدوال أحادية البعد، لأن جميع المتغيرات في الحساب تعتمد فقط على المعامل القطري لعلاقة التشتت. تُعدّ السوائل والزجاج والمواد الصلبة غير المتبلورة أمثلة على نظام متناظر، حيث تتمتع علاقات التشتت فيه بتناظر دوراني.

المجسم الثماني.

تتطلب القياسات على المساحيق أو العينات متعددة البلورات تقييم وحساب الدوال والتكاملات على كامل نطاق علاقات التشتت للنظام المدروس ، وغالبًا ما يكون هذا النطاق هو نطاق بريلوين . في بعض الأحيان، يكون تناظر النظام عاليًا، مما يؤدي إلى ظهور شكل الدوال التي تصف علاقات التشتت عدة مرات على كامل نطاق علاقة التشتت. في مثل هذه الحالات، يمكن تقليل الجهد المبذول لحساب كثافة الحالات بشكل كبير عند حصر الحساب في نطاق مُصغّر أو نطاق أساسي . [ 1 ] يتميز نطاق بريلوين للشبكة المكعبة ذات المراكز الوجهية (FCC) في الشكل على اليمين بتناظر 48 طية للمجموعة النقطية O h مع تناظر ثماني السطوح كامل . هذا التكوين يعني أنه يمكن اختزال التكامل على كامل نطاق بريلوين إلى جزء من 48 من نطاق بريلوين بأكمله. كما يُظهر الجدول الدوري للبنية البلورية ، توجد العديد من العناصر ذات البنية البلورية المكعبة المتمركزة على الوجوه (FCC)، مثل الماس والسيليكون والبلاتين ، وتتميز مناطق بريلوين وعلاقات التشتت الخاصة بها بتناظر 48. ومن بين البنى البلورية الشائعة الأخرى، الشبكة المكعبة المتمركزة على الجسم (BCC) والبنية السداسية المتراصة (HCP)، حيث تتميز الأولى بشبكة مكعبة والثانية بشبكة سداسية. تتمتع بنية BCC بتناظر بيريتوهيدرالي 24 لمجموعة التناظر T h . أما بنية HCP، فتتمتع بتناظر ثنائي السطوح موشوري 12 لمجموعة التناظر D 3h . ويمكن الاطلاع على قائمة كاملة بخصائص التناظر لمجموعة التناظر في جداول خصائص مجموعات التناظر .

بشكل عام، يسهل حساب كثافة الحالات عندما يكون تناظر النظام أعلى وعدد الأبعاد الطوبولوجية لعلاقة التشتت أقل. غالبًا ما يمكن حساب كثافة الحالات لعلاقات التشتت ذات التناظر الدوراني تحليليًا. هذه النتيجة مفيدة، نظرًا لأن العديد من المواد ذات الأهمية العملية، مثل الفولاذ والسيليكون، تتمتع بتناظر عالٍ.

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

طوبولوجيا الفضاء k

الشكل 1: سطح كروي في فضاء k للإلكترونات في ثلاثة أبعاد.

تعتمد كثافة الحالات على الأبعاد الكلية للجسم. في نظام ثلاثي الأبعاد (ثلاثة أبعاد) تُقاس كثافة الحالات بوحدة الطاقة/ الحجم، وفي نظام ثنائي الأبعاد تُقاس بوحدة الطاقة/المساحة، وفي نظام أحادي الأبعاد تُقاس بوحدة الطاقة/الطول. يُشير الحجم المذكور إلى حجم فضاء k، وهو الفضاء المحصور بسطح الطاقة الثابتة للنظام، والمُستنتج من علاقة التشتت التي تربط الطاقة (E) بـ k . يوضح الشكل 1 مثالًا على فضاء k ثلاثي الأبعاد . يتضح من الشكل أن أبعاد النظام تُقيد زخم الجسيمات داخله.

كثافة حالات متجه الموجة (كرة)

تبدأ عملية حساب كثافة الحالات (DOS) بحساب عدد الحالات المسموح بها N عند قيمة k معينة ، والتي تقع ضمن النطاق [ k , k + dk ] داخل حجم النظام. تتم هذه العملية عن طريق اشتقاق حجم فضاء k بأكمله.Ωن،ك{\displaystyle \Omega _{n,k}}في فضاءات كروية ذات n بُعد عند قيمة k عشوائية ، بالنسبة إلى k . يُعبَّر عن الحجم أو المساحة أو الطول في فضاءات k الكروية ثلاثية الأبعاد أو ثنائية الأبعاد أو أحادية الأبعاد بالصيغة التالية:Ωن(ك)=جنكن{\displaystyle \Omega _{n}(k)=c_{n}k^{n}}

بالنسبة لفضاء k ذي n بُعد مع الثوابت المحددة طوبولوجيًا ج1=2، ج2=π، ج3=4π3{\displaystyle c_{1}=2,\ c_{2}=\pi ,\ c_{3}={\frac {4\pi }{3}}}بالنسبة للدوال الخطية والقرصية والكروية المتناظرة الشكل في فضاءات k الإقليدية أحادية وثنائية وثلاثية الأبعاد على التوالي.

وفقًا لهذا المخطط، فإن كثافة حالات متجه الموجة N هي، من خلال التفاضلΩن،ك{\displaystyle \Omega _{n,k}}بالنسبة إلى k ، معبرًا عنها بواسطة شمالن(ك)=دΩن(ك)دك=نجنكن-1{\displaystyle N_{n}(k)={\frac {\mathrm {d} \Omega _{n}(k)}{\mathrm {d} k}}=n\;c_{n}\;k^{n-1}}

تُكتب كثافة حالات متجه الموجة أحادية وثنائية وثلاثية الأبعاد لخط أو قرص أو كرة بشكل صريح على النحو التالي شمال1(ك)=2شمال2(ك)=2πكشمال3(ك)=4πك2{\displaystyle {\begin{aligned}N_{1}(k)&=2\\N_{2}(k)&=2\pi k\\N_{3}(k)&=4\pi k^{2}\end{aligned}}}

تكون الحالة الواحدة كبيرة بما يكفي لاحتواء جسيمات ذات طول موجي λ. ويرتبط الطول الموجي بـ k من خلال العلاقة التالية. ك=2πλ{\displaystyle k={\frac {2\pi }{\lambda }}}

في النظام الكمومي، يعتمد طول λ على المسافة المميزة للنظام L التي تحصر الجسيمات. وأخيرًا، تُضرب كثافة الحالات N بمعامل.s/Vك{\displaystyle s/V_{k}}حيث يمثل s عامل انحلال ثابتًا يأخذ في الاعتبار درجات الحرية الداخلية الناتجة عن ظواهر فيزيائية مثل الدوران أو الاستقطاب. إذا لم تكن هذه الظاهرة موجودة، فإنs=1{\displaystyle s=1}. V k هو الحجم في فضاء k الذي تكون متجهات الموجة الخاصة به أصغر من أصغر متجهات الموجة الممكنة التي تحددها المسافة المميزة للنظام.

كثافة حالات الطاقة

لإكمال حساب كثافة الحالات، أوجد عدد الحالات لكل وحدة حجم عينة عند طاقة معينة.هـ{\displaystyle E}داخل فترة زمنية[هـ،هـ+دهـ]{\displaystyle [E,E+\mathrm {d} E]}يُعطى الشكل العام لتوزيع كثافة الحالات (DOS) لنظام ما على النحو التالي: دن(هـ)=دΩن(هـ)دهـ{\displaystyle D_{n}\left(E\right)={\frac {\mathrm {d} \Omega _{n}(E)}{\mathrm {d} E}}} لا ينطبق المخطط الموضح حتى الآن إلا على علاقات التشتت المتزايدة بشكل رتيب والمتناظرة كرويًا . وبشكل عام، علاقة التشتتهـ(ك){\displaystyle E(k)}ليست متناظرة كرويًا، وفي كثير من الحالات لا ترتفع باستمرار. للتعبير عن D كدالة لـ فإن معكوس علاقة التشتتهـ(ك){\displaystyle E(k)}يجب استبدالها في تعبيرΩن(ك){\displaystyle \Omega _{n}(k)}كدالة لـ k للحصول على تعبيرΩن(هـ){\displaystyle \Omega _{n}(E)}كدالة للطاقة. إذا لم تكن علاقة التشتت متناظرة كرويًا أو متزايدة باستمرار، ولم يكن من السهل عكسها، ففي معظم الحالات يجب حساب كثافة الحالات عدديًا. تتوفر اشتقاقات أكثر تفصيلًا. [ 2 ] [ 3 ]

علاقات التشتت

يتم تحديد علاقة التشتت للإلكترونات في مادة صلبة من خلال بنية النطاق الإلكتروني .

تعتمد الطاقة الحركية للجسيم على مقدار واتجاه متجه الموجة k ، وخصائص الجسيم، والبيئة التي يتحرك فيها. على سبيل المثال، تُعطى الطاقة الحركية للإلكترون في غاز فيرمي بالعلاقة التالية:هـ=هـ0+(ك)22م ،{\displaystyle E=E_{0}+{\frac {\left(\hbar k\right)^{2}}{2m}}\ ,}

حيث m هي كتلة الإلكترون . علاقة التشتت عبارة عن قطع مكافئ متناظر كرويًا وهي تتزايد باستمرار، لذا يمكن حساب كثافة الحالات بسهولة.

الشكل 2: علاقة تشتت الفونون في السلسلة أحادية الذرة

بالنسبة للفونونات الطولية في سلسلة من الذرات، تُعطى علاقة التشتت للطاقة الحركية في فضاء k أحادي البعد ، كما هو موضح في الشكل 2، بالصيغة التالية: هـ=2ω0|الخطيئة(كأ2)|{\displaystyle E=2\hbar \omega _{0}\left|\sin \left({\frac {ka}{2}}\right)\right|} أينω0=كF/م{\textstyle \omega _{0}={\sqrt {k_{\text{F}}/m}}}تردد المذبذب،م{\displaystyle m}كتلة الذرات،كF{\displaystyle k_{\text{F}}}ثابت القوة بين الذرات وأ{\displaystyle a}المسافة بين الذرات. بالنسبة للقيم الصغيرة لـكπ/أ{\displaystyle k\ll \pi /a}العلاقة التشتتية خطية: هـ=ω0كأ{\displaystyle E=\hbar \omega _{0}ka}

متىكπ/أ{\displaystyle k\approx \pi /a}الطاقة هي هـ=2ω0|كوس(π-كأ2)|{\displaystyle E=2\hbar \omega _{0}\left|\cos \left({\frac {\pi -ka}{2}}\right)\right|}

مع التحولq=ك-π/أ{\displaystyle q=k-\pi /a}وصغيرةq{\displaystyle q}يمكن تحويل هذه العلاقة إلى هـ2ω0[1-(qأ2)2]{\displaystyle E\approx 2\hbar \omega _{0}\left[1-\left({\frac {qa}{2}}\right)^{2}\right]}

علاقات التشتت المتساوية الخواص

يمكن التعبير عن المثالين المذكورين هنا على النحو التالي: هـ=هـ0+جككص{\displaystyle E=E_{0}+c_{k}k^{p}}

هذا التعبير هو نوع من علاقات التشتت لأنه يربط بين خاصيتين موجيتين، وهو متناحٍ لأن طول متجه الموجة فقط، وليس اتجاهه، يظهر في التعبير. يرتبط مقدار متجه الموجة بالطاقة كما يلي: ك=(هـ-هـ0جك)1/ص،{\displaystyle k=\left({\frac {E-E_{0}}{c_{k}}}\right)^{{1}/{p}},}

وبناءً على ذلك، فإن حجم فضاء k ذي الأبعاد n الذي يحتوي على متجهات موجية أصغر من k هو: Ωن(ك)=جنكن{\displaystyle \Omega _{n}(k)=c_{n}k^{n}}

يُعطي استبدال علاقة الطاقة المتناحية حجم الحالات المشغولة Ωن(هـ)=جنجكنص(هـ-هـ0)ن/ص،{\displaystyle \Omega _{n}(E)={\frac {c_{n}}{{c_{k}}^{\frac {n}{p}}}}\left(E-E_{0}\right)^{{n}/{p}},}

يُعطي تفاضل هذا الحجم بالنسبة للطاقة تعبيرًا عن كثافة الحالات لعلاقة التشتت المتناحية دن(هـ)=ددهـΩن(هـ)=نجنصجكنص(هـ-هـ0)نص-1{\displaystyle D_{n}\left(E\right)={\frac {\mathrm {d} }{\mathrm {d} E}}\Omega _{n}(E)={\frac {nc_{n}}{p{c_{k}}^{\frac {n}{p}}}}\left(E-E_{0}\right)^{{\frac {n}{p}}-1}}

التشتت المكافئ

الشكل 3: كثافة حالات الإلكترون الحر في فضاء k ثلاثي الأبعاد

في حالة علاقة التشتت المكافئة ( p = 2)، كما هو الحال بالنسبة للإلكترونات الحرة في غاز فيرمي، فإن كثافة الحالات الناتجة،دن(هـ){\displaystyle D_{n}\left(E\right)}، بالنسبة للإلكترونات في الأنظمة ذات الأبعاد n هو

د1(هـ)=1جك(هـ-هـ0)د2(هـ)=πجكد3(هـ)=2πهـ-هـ0جك3.{\displaystyle {\begin{aligned}D_{1}\left(E\right)&={\frac {1}{\sqrt {c_{k}\left(E-E_{0}\right)}}}\\[1ex]D_{2}\left(E\right)&={\frac {\pi }{c_{k}}}\\[1ex]D_{3}\left(E\right)&=2\pi {\sqrt {\frac {E-E_{0}}{c_{k}^{3}}}}\,.\end{aligned}}}

لهـ>هـ0{\displaystyle E>E_{0}}، معد(هـ)=0{\displaystyle D(E)=0}لهـ<هـ0{\displaystyle E<E_{0}}.

في الأنظمة أحادية البعد، تتباعد كثافة الحالات عند أسفل النطاق كماهـ{\displaystyle E}انخفاض إلىهـ0{\displaystyle E_{0}}في الأنظمة ثنائية الأبعاد، يتبين أن كثافة الحالات مستقلة عنهـ{\displaystyle E}وأخيرًا، بالنسبة للأنظمة ثلاثية الأبعاد، ترتفع كثافة الحالات كجذر تربيعي للطاقة. [ 4 ]

بما في ذلك العامل المسبقs/Vك{\displaystyle s/V_{k}}، التعبير عن كثافة الحالات ثلاثية الأبعاد هو شمال(هـ)=V2π2(2م2)32هـ-هـ0،{\displaystyle N(E)={\frac {V}{2\pi ^{2}}}\left({\frac {2m}{\hbar ^{2}}}\right)^{\frac {3}{2}}{\sqrt {E-E_{0}}},}

أينV{\displaystyle V}هو الحجم الإجمالي، وشمال(هـ-هـ0){\displaystyle N(E-E_{0})}يتضمن ذلك انحلال الدوران الثنائي.

التشتت الخطي

في حالة العلاقة الخطية ( p = 1)، كما هو الحال بالنسبة للفوتونات أو الفونونات الصوتية أو بعض الأنواع الخاصة من النطاقات الإلكترونية في المواد الصلبة، فإن كثافة الحالات في الأنظمة أحادية وثنائية وثلاثية الأبعاد ترتبط بالطاقة على النحو التالي:

د1(هـ)=2πجكد2(هـ)=2πهـ-هـ0جك2د3(هـ)=4π(هـ-هـ0)2جك3{\displaystyle {\begin{aligned}D_{1}\left(E\right)&={\frac {2\pi }{c_{k}}}\\[1ex]D_{2}\left(E\right)&=2\pi {\frac {E-E_{0}}{c_{k}^{2}}}\\[1ex]D_{3}\left(E\right)&=4\pi {\frac {\left(E-E_{0}\right)^{2}}{c_{k}^{3}}}\end{aligned}}}

دوال التوزيع

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

الشكل 4: توزيع احتمالية فيرمي-ديراك، وكثافة الحالات، وحاصل ضربهما لأشباه الموصلات. يمثل الفص الأخضر السفلي طاقة الثقوب ، وبالتالي يستخدم      و(-x){\displaystyle f(-x)}كدالة توزيع.

إحصاءات فيرمي-ديراك : تُستخدم دالة توزيع احتمالية فيرمي-ديراك، الموضحة في الشكل 4، لإيجاد احتمالية وجود فرميون في حالة كمومية محددة ضمن نظام في حالة توازن حراري. الفرميونات هي جسيمات تخضع لمبدأ استبعاد باولي (مثل الإلكترونات والبروتونات والنيوترونات). يمكن كتابة دالة التوزيع على النحو التالي: وFد(هـ)=1خبرة(هـ-μكبتي)+1.{\displaystyle f_{\mathrm {FD} }(E)={\frac {1}{\exp \left({\frac {E-\mu }{k_{\mathrm {B} }T}}\right)+1}}.}

μ{\displaystyle \mu }هو الجهد الكيميائي (يرمز له أيضًا بـ E F ويسمى مستوى فيرمي عندما T = 0)،كب{\displaystyle k_{\mathrm {B} }}هو ثابت بولتزمان، وتي{\displaystyle T}درجة الحرارة. يوضح الشكل 4 كيف يمكن لناتج دالة توزيع فيرمي-ديراك وكثافة الحالات ثلاثية الأبعاد لأشباه الموصلات أن يعطي نظرة ثاقبة على الخصائص الفيزيائية مثل تركيز حاملات الشحنة وفجوات نطاق الطاقة.

إحصاءات بوز-أينشتاين : دالة توزيع بوز-أينشتاينوبهـ(هـ){\displaystyle f_{\mathrm {BE} }(E)}هو متوسط ​​عدد البوزونات التي تشغل أي حالة كمومية ذات طاقةهـ{\displaystyle E}في نظام في حالة توازن حراري. البوزونات هي جسيمات لا تخضع لمبدأ استبعاد باولي (مثل الفونونات والفوتونات)، لذلكوبهـ(هـ){\displaystyle f_{\mathrm {BE} }(E)}يمكن أن تكون أكبر من 1. يمكن كتابة دالة التوزيع على النحو التالي: وبهـ(هـ)=1خبرة(هـ-μكبتي)-1.{\displaystyle f_{\mathrm {BE} }(E)={\frac {1}{\exp \left({\frac {E-\mu }{k_{\text{B}}T}}\right)-1}}.}

من خلال هذين التوزيعين، يمكن حساب خصائص مثل الطاقة الداخلية لكل وحدة حجمu{\displaystyle u}عدد الجسيماتشمال{\displaystyle N}السعة الحرارية النوعيةج{\displaystyle c}، والتوصيل الحراريك{\displaystyle k}العلاقات بين هذه الخصائص وحاصل ضرب كثافة الحالات وتوزيع الاحتمال، مع الإشارة إلى كثافة الحالات بـز(هـ){\displaystyle g(E)}بدلاً مند(هـ){\displaystyle D(E)}يتم تقديمها بواسطة u=هـو(هـ)ز(هـ)دهـشمال=Vو(هـ)ز(هـ)دهـج=تيهـو(هـ)ز(هـ)دهـك=1دتيهـو(هـ)ز(هـ)ν(هـ)Λ(هـ)دهـ{\displaystyle {\begin{aligned}u&=\int E\,f(E)\,g(E)\,\mathrm {d} E\\[1ex]N&=V\int f(E)\,g(E)\,\mathrm {d} E\\[1ex]c&={\frac {\partial }{\partial T}}\int E\,f(E)\,g(E)\,\mathrm {d} E\\[1ex]k&={\frac {1}{d}}{\frac {\partial }{\partial T}}\int Ef(E)\,g(E)\,\nu (E)\,\Lambda (E)\,\mathrm {d} E\end{aligned}}}

د{\displaystyle d}الأبعاد،ν{\displaystyle \nu }هي سرعة الصوت وΛ{\displaystyle \Lambda }يعني المسار الحر .

التطبيقات

تظهر كثافة الحالات في العديد من مجالات الفيزياء، وتساعد في تفسير عدد من الظواهر الميكانيكية الكمومية.

التكميم

يُظهر حساب كثافة الحالات للهياكل الصغيرة أن توزيع الإلكترونات يتغير مع انخفاض الأبعاد. ففي الأسلاك الكمومية ، تصبح كثافة الحالات عند طاقات معينة أعلى من كثافة الحالات في أشباه الموصلات الصلبة، أما في النقاط الكمومية، فتُكمّم الإلكترونات عند طاقات معينة.

البلورات الضوئية

يمكن التحكم في كثافة حالات الفوتون باستخدام هياكل دورية بأطوال تقارب طول موجة الضوء. بعض هذه الهياكل قادرة على منع انتشار الضوء ذي ألوان (طاقات) معينة تمامًا، مما يُنشئ فجوة نطاق فوتوني: حيث تكون كثافة الحالات صفرًا لتلك الطاقات. بينما تُعيق هياكل أخرى انتشار الضوء في اتجاهات محددة فقط، لتُشكل مرايا، وموجهات موجية، وتجاويف. تُعرف هذه الهياكل الدورية بالبلورات الفوتونية . [ 5 ] [ 6 ] [ 7 ] [ 8 ] في الأوساط النانوية، غالبًا ما يكون مفهوم كثافة الحالات المحلية (LDOS) أكثر أهمية من مفهوم كثافة الحالات العامة (DOS)، نظرًا لاختلاف كثافة الحالات العامة اختلافًا كبيرًا من نقطة إلى أخرى.

الحساب الحسابي

تتسم الأنظمة المثيرة للاهتمام عمومًا بالتعقيد، مثل المركبات والجزيئات الحيوية والبوليمرات، وما إلى ذلك. ونظرًا لتعقيد هذه الأنظمة، فإن الحساب التحليلي لكثافة الحالات يكون مستحيلاً في معظم الحالات. توفر المحاكاة الحاسوبية مجموعة من الخوارزميات لتقييم كثافة الحالات بدقة عالية. إحدى هذه الخوارزميات تُعرف بخوارزمية وانغ ولاندو . [ 9 ]

في مخطط وانغ ولاندو، لا يُشترط معرفة مسبقة بكثافة الحالات. ويتم ذلك على النحو التالي: تُجزأ دالة التكلفة (مثل الطاقة) للنظام. وفي كل مرة يتم الوصول إلى الخانة i ، يتم تحديث الرسم البياني لكثافة الحالات.ز(أنا){\displaystyle g(i)}، بواسطة ز(أنا)ز(أنا)+و{\displaystyle g(i)\rightarrow g(i)+f} حيث يُطلق على f اسم عامل التعديل. بمجرد زيارة كل خانة في الرسم البياني عددًا معينًا من المرات (10-15)، يتم تقليل عامل التعديل وفقًا لمعيار معين، على سبيل المثال، ون+112ون{\displaystyle f_{n+1}\rightarrow {\frac {1}{2}}f_{n}} حيث يرمز n إلى خطوة التحديث رقم n . تنتهي المحاكاة عندما يكون عامل التعديل أقل من عتبة معينة، على سبيل المثالون<10-8{\displaystyle f_{n}<10^{-8}}.

تتميز خوارزمية وانغ ولاندو ببعض المزايا مقارنةً بالخوارزميات الشائعة الأخرى، مثل المحاكاة متعددة القنوات والتسخين المتوازي . فعلى سبيل المثال، تُحسب كثافة الحالات كناتج رئيسي للمحاكاة. إضافةً إلى ذلك، فإن محاكاة وانغ ولاندو مستقلة تمامًا عن درجة الحرارة. وتتيح هذه الميزة حساب كثافة الحالات لأنظمة ذات سطح طاقة غير منتظم للغاية، مثل البروتينات. [ 10 ]

رياضياً، تتم صياغة كثافة الحالات من حيث برج من الخرائط التغطية. [ 11 ]

الكثافة المحلية للحالات

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

مفهوم

يصف مصطلح كثافة الحالات المحلية (LDOS) كثافة الحالات المُحَلَّلة مكانيًا. في علم المواد، على سبيل المثال، يُعد هذا المصطلح مفيدًا عند تفسير البيانات المُستقاة من مجهر المسح النفقي (STM)، نظرًا لقدرة هذه الطريقة على تصوير كثافات حالات الإلكترون بدقة ذرية. وبناءً على بنية البلورة، يُمكن التنبؤ بهذه الكمية باستخدام الطرق الحسابية، مثل نظرية الكثافة الوظيفية .

تعريف عام

في الكثافة المحلية للحالات، يتم ترجيح مساهمة كل حالة بكثافة دالة الموجة الخاصة بها عند النقطة.شمال(هـ){\displaystyle N(E)}يصبحن(هـ،x){\displaystyle n(E,x)}

ن(هـ،x)=ج|ϕج(x)|2دلتا(هـ-εج){\displaystyle n(E,x)=\sum _{j}|\phi _{j}(x)|^{2}\delta (E-\varepsilon _{j})}

عامل|ϕج(x)|2{\displaystyle |\phi _{j}(x)|^{2}}وهذا يعني أن كل ولاية تساهم بشكل أكبر في المناطق ذات الكثافة السكانية العالية. متوسط ​​علىx{\displaystyle x}سيؤدي هذا التعبير إلى استعادة الصيغة المعتادة لكثافة الحالات. تُعد كثافة الحالات المحلية مفيدة في الأنظمة غير المتجانسة، حيثن(هـ،x){\displaystyle n(E,x)}يحتوي على معلومات أكثر منن(هـ){\displaystyle n(E)}وحيد.

بالنسبة لنظام أحادي البعد ذي جدار، تعطي الموجات الجيبية

ن1د(هـ،x)=2π2مهـالخطيئة2كx{\displaystyle n_{1D}(E,x)={\frac {2}{\pi \hbar }}{\sqrt {\frac {2m}{E}}}\sin ^{2}{kx}}

أينك=2مهـ/{\textstyle k={\sqrt {2mE}}/\hbar }.

في نظام ثلاثي الأبعاد معx>0{\displaystyle x>0}التعبير هو

ن3د(هـ،x)=(1-الخطيئة2كx2كx)ن3د(هـ){\displaystyle n_{3D}(E,x)=\left(1-{\frac {\sin {2kx}}{2kx}}\right)n_{3D}(E)}

في الواقع، يمكننا تعميم الكثافة المحلية للحالات بشكل أكبر إلى

ن(هـ،x،x)=جϕج(x)ϕج*(x)دلتا(هـ-εج){\displaystyle n(E,x,x')=\sum _{j}\phi _{j}(x)\phi _{j}^{*}(x')\delta (E-\varepsilon _{j})}

تُسمى هذه الدالة بالدالة الطيفية ، وهي دالة تحتوي على كل دالة موجية على حدة في متغيرها الخاص. في النظريات المتقدمة، ترتبط هذه الدالة بدوال غرين، وتوفر تمثيلاً موجزاً لبعض النتائج، مثل الامتصاص الضوئي .

كثافة الحالات المحلية المُحَدَّدة مكانيًا. سلسلة من الصور مع اختلاف انحياز البوابة في ترانزستور MOSFET نانوي السلك عند انحياز المصرف {{{1}}} لاحظ مستويات الطاقة المحصورة أثناء تحركها مع زيادة انحياز البوابة.

أجهزة الحالة الصلبة

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

البصريات والفوتونيات

في مجال البصريات والفوتونيات ، يشير مفهوم الكثافة الموضعية للحالات إلى الحالات التي يمكن أن يشغلها الفوتون. بالنسبة للضوء، تُقاس عادةً باستخدام طرق التألق، أو طرق المسح في المجال القريب، أو تقنيات التألق الكاثودي. تختلف سلوكيات الكثافة الموضعية للحالات باختلاف البنى الفوتونية، مما يؤثر بشكل متباين على الانبعاث التلقائي. في البلورات الفوتونية، يُتوقع أن تكون الكثافة الموضعية للحالات قريبة من الصفر، مما يُثبط الانبعاث التلقائي. [ 12 ] كما يُتوقع تعزيز مماثل للكثافة الموضعية للحالات في تجويف البلازمون. [ 13 ] مع ذلك، في البنى النانوية الفوتونية غير المنتظمة، تتصرف الكثافة الموضعية للحالات بشكل مختلف. فهي تتذبذب مكانيًا وفقًا لإحصائياتها، وتتناسب طرديًا مع قوة تشتت البنى. [ 14 ] إضافةً إلى ذلك، فإن العلاقة مع متوسط ​​المسار الحر للتشتت بسيطة، إذ يمكن أن تتأثر الكثافة الموضعية للحالات بشدة بالتفاصيل الدقيقة للاضطرابات القوية في صورة تعزيز قوي لانبعاث بورسل. [ 15 ] وأخيرًا، بالنسبة للاضطراب البلازموني، يكون هذا التأثير أقوى بكثير بالنسبة لتقلبات كثافة الحالات المحلية، حيث يمكن ملاحظته على أنه توضع قوي في المجال القريب. [ 16 ]

انظر أيضاً

مراجع

  1. والتر آشلي هاريسون (1989). البنية الإلكترونية وخواص المواد الصلبة . منشورات دوفر. ISBN 978-0-486-66021-9.
  2. حساب كثافة الحالات للعينة
  3. "حساب آخر لكثافة الحالات" . مؤرشف من الأصل بتاريخ 16-05-2008 . تم الاطلاع عليه بتاريخ 22-05-2008 .
  4. تشارلز كيتل (1996). مقدمة في فيزياء الحالة الصلبة ( الطبعة السابعة). وايلي. المعادلة (37)، ص 216. ISBN  978-0-471-11181-8.
  5. يابلونوفيتش، إي. (1987). "الانبعاث التلقائي المثبط في فيزياء الحالة الصلبة والإلكترونيات" . مجلة Physical Review Letters ، 58 (20): 2059-2062 . Bibcode : 1987PhRvL..58.2059Y . doi : 10.1103/PhysRevLett.58.2059 . PMID 10034639 . 
  6. جون، ساجيف؛ وانغ، جيان (1990). "الديناميكا الكهربائية الكمومية بالقرب من فجوة نطاق الفوتون: حالات الفوتون المقيدة والذرة الملبسة". مجلة Physical Review Letters ، 64 (20): 2418-2421 . Bibcode : 1990PhRvL..64.2418J . doi : 10.1103/PhysRevLett.64.2418 . PMID: 10041707 . 
  7. لودال، ب.؛ فان دريل، أ.ف.؛ نيكولايف، إ. (2004). "التحكم في ديناميكيات الانبعاث التلقائي من النقاط الكمومية بواسطة البلورات الضوئية". مجلة نيتشر . 430 (1): 654-657 . Bibcode : 2004Natur.430..654L . doi : 10.1038/nature02772 . hdl : 1874/16698 . PMID 15295594. S2CID 4334567 .  
  8. فوجيتا، ماسايوكي؛ تاكاهاشي، شيجيكي؛ تاناكا، يوشينوري؛ أسانو، تاكاشي؛ نودا، سوسومو (2005). "التثبيط وإعادة توزيع انبعاث الضوء التلقائي في البلورات الفوتونية". مجلة ساينس . 308 (5726): 1296-1298 . Bibcode : 2005Sci...308.1296F . doi : 10.1126/science.1110417 . PMID 15919989. S2CID 30116866 .  
  9. وانغ، فوغاو؛ لاندو، د.ب. (2001). "خوارزمية فعّالة للمشي العشوائي متعدد النطاقات لحساب كثافة الحالات". مجلة Physical Review Letters ، 86 (10): 2050-2053 . arXiv : cond-mat/0011174 . Bibcode : 2001PhRvL..86.2050W . doi : 10.1103/PhysRevLett.86.2050 . PMID : 11289852. S2CID : 2941153 .  
  10. أوجيدا، ب.؛ غارسيا، م. (2010). "التأثير الكهربائي على بنية البروتين الأصلية من نوع بيتا-شيت وتكوين بنية حلزونية" . مجلة الفيزياء الحيوية . 99 (2): 595-599 . Bibcode : 2010BpJ....99..595O . doi : 10.1016/j.bpj.2010.04.040 . PMC 2905109. PMID 20643079 .  
  11. أداتشي، ت. وسونادا، ت. (1993). "كثافة الحالات في الهندسة الطيفية". تعليقات الرياضيات الهيلفيتية 68 : 480-493 . doi : 10.1007/BF02565831 . S2CID 120828817 . 
  12. سبريك، ر.؛ فان تيغيلين، ب.أ.؛ لاغينديك، أ. (1996). "كثافة الحالات في الهندسة الطيفية" . رسائل الفيزياء الأوروبية. 35 ( 4): 265-270 . doi : 10.1209/epl/i1996-00564-y . S2CID 250854036 . 
  13. فرحاني، ج. ن.؛ بول، د. و.؛ إيسلر، هـ.-ج.؛ هيشت، ب. (2005). "نقطة كمومية مفردة مقترنة بهوائي بصري ماسح: باعث فائق قابل للضبط". مجلة Physical Review Letters، 95 ( 1 ) 017402. رمز Bibcode : 2005PhRvL..95a7402F . doi : 10.1103/PhysRevLett.95.017402 . PMID 16090656 . 
  14. بيرووسوتو، م.؛ سكيبتروف، س.؛ فوس، و.؛ موسك، أ. (2010). "ملاحظة التقلبات المكانية للكثافة المحلية للحالات في أوساط ضوئية عشوائية". مجلة Physical Review Letters ، 105 ( 1) 013904. arXiv : 1002.3186 . Bibcode : 2010PhRvL.105a3904B . doi : 10.1103/PhysRevLett.105.013904 . PMID 20867448. S2CID 25044558 .  
  15. سابينزا، ر.؛ بونداريف، ب.؛ بييرات، ر.؛ هابرت، ب.؛ كارميناتي، ر.؛ فان هولست، ن. ف. (2011). "إحصائيات الذيل الطويل لعامل بورسيل في الأوساط غير المنتظمة المدفوعة بتفاعلات المجال القريب". مجلة Physical Review Letters ، 106 (16)، 163902. Bibcode : 2011PhRvL.106p3902S . doi : 10.1103/PhysRevLett.106.163902 . PMID: 21599367 . 
  16. كراخمالنيكوف، ف.؛ كاستاني، إ.؛ دي وايلد، ي.؛ كارميناتي، ر. (2010). "إحصائيات الذيل الطويل لعامل بورسيل في الأوساط غير المنتظمة المدفوعة بتفاعلات المجال القريب". مجلة Physical Review Letters ، 105 (18) 183901. arXiv : 1007.3691 . doi : 10.1103/PhysRevLett.105.183901 . PMID 21231105. S2CID 15590513 .  

للمزيد من القراءة

  • تشين، غانغ. نقل الطاقة وتحويلها على المستوى النانوي. نيويورك: أكسفورد، 2005
  • ستريتمان، بن جي. وسانجاي بانيرجي. أجهزة إلكترونية صلبة الحالة. أبر سادل ريفر، نيوجيرسي: برنتيس هول، 2000.
  • مولر، ريتشارد إس. وثيودور آي. كامينز. إلكترونيات الأجهزة للدوائر المتكاملة. نيويورك: جون وايلي وأولاده، 2003.
  • كيتل، تشارلز وهربرت كرومر. الفيزياء الحرارية. نيويورك: دبليو إتش فريمان وشركاه، 1980
  • سزي، سيمون م. فيزياء أجهزة أشباه الموصلات. نيويورك: جون وايلي وأولاده، 1981