نظرية بلوخ

السطح المتساوي لمعامل التربيع لحالة بلوخ في شبكة السيليكون
الخط المتصل: رسم تخطيطي للجزء الحقيقي من حالة بلوخ النموذجية في بُعد واحد. الخط المنقط ناتج عن العامل e i k · r . الدوائر الفاتحة تمثل الذرات.

في فيزياء المادة المكثفة ، تنص نظرية بلوخ على أن حلول معادلة شرودنغر في جهد دوري يمكن التعبير عنها كموجات مستوية مُعدَّلة بدوال دورية . سُميت النظرية نسبةً إلى الفيزيائي السويسري فيليكس بلوخ ، الذي اكتشفها عام ١٩٢٩. [ ١ ] رياضياً، تُكتب [ ٢ ]

دالة بلوخ

ψ(ر)=هـأناكرu(ر){\displaystyle \psi (\mathbf {r} )=e^{i\mathbf {k} \cdot \mathbf {r} }u(\mathbf {r} )}

أينر{\displaystyle \mathbf {r} }هو الموقع،ψ{\displaystyle \psi }هي دالة الموجة ،u{\displaystyle u}هي دالة دورية لها نفس دورية البلورة، متجه الموجةك{\displaystyle \mathbf {k} }هو متجه زخم البلورة ،هـ{\displaystyle e}هو عدد أويلر ، وأنا{\displaystyle i}هي الوحدة التخيلية .

تُعرف الدوال من هذا الشكل باسم دوال بلوخ أو حالات بلوخ ، وهي بمثابة أساس مناسب للدوال الموجية أو حالات الإلكترونات في المواد الصلبة البلورية .

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

تُكتب هذه الحالات الذاتية باستخدام الرموز السفلية كما يلي:ψنك{\displaystyle \psi _{n\mathbf {k} }}، أينن{\displaystyle n}هو مؤشر منفصل، يُسمى مؤشر النطاق ، وهو موجود لأن هناك العديد من الدوال الموجية المختلفة التي لها نفسك{\displaystyle \mathbf {k} }(لكل منها مكون دوري مختلف)u{\displaystyle u}). ضمن نطاق (أي، لفترة ثابتةن{\displaystyle n})ψنك{\displaystyle \psi _{n\mathbf {k} }}يتغير باستمرار معك{\displaystyle \mathbf {k} }وكذلك طاقتها. أيضاً،ψنك{\displaystyle \psi _{n\mathbf {k} }}يكون فريدًا فقط حتى متجه الشبكة المقلوبة الثابتك{\displaystyle \mathbf {K} }، أو،ψنك=ψن(ك+ك){\displaystyle \psi _{n\mathbf {k} }=\psi _{n(\mathbf {k+K} )}}لذلك، فإن متجه الموجةك{\displaystyle \mathbf {k} }يمكن حصرها في منطقة بريلوين الأولى للشبكة المتبادلة دون فقدان العمومية .

التطبيقات والنتائج

قابلية التطبيق

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

متجه الموجة

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

لنفترض أن الإلكترون في حالة بلوخ ψ(ر)=هـأناكرu(ر)،{\displaystyle \psi (\mathbf {r} )=e^{i\mathbf {k} \cdot \mathbf {r} }u(\mathbf {r} ),} حيث تكون u دورية بنفس دورية الشبكة البلورية. وتُحدد الحالة الكمومية الفعلية للإلكترون بشكل كامل بواسطةψ{\displaystyle \psi }ليس k أو u بشكل مباشر. هذا مهم لأن k و u ليسا فريدين . تحديدًا، إذاψ{\displaystyle \psi }يمكن كتابة المتجه كما سبق باستخدام k ، ويمكن كتابته أيضًا باستخدام ( k + K ) ، حيث K هو أي متجه شبكي مقلوب (انظر الشكل على اليمين). لذلك، فإن المتجهات الموجية التي تختلف بمتجه شبكي مقلوب متكافئة، بمعنى أنها تميز نفس مجموعة حالات بلوخ.

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

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

مثال مفصل

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

إفادة

نظرية بلوخ - بالنسبة للإلكترونات في بلورة مثالية، يوجد أساس من الدوال الموجية بالخاصيتين التاليتين:

  • كل دالة من هذه الدوال الموجية هي حالة طاقة ذاتية،
  • كل دالة موجية من هذه الدوال هي حالة بلوخ، مما يعني أن هذه الدالة الموجيةψ{\displaystyle \psi }يمكن كتابتها بالشكل التاليψ(ر)=هـأناكرu(ر)،{\displaystyle \;\psi (\mathbf {r} )=e^{i\mathbf {k} \cdot \mathbf {r} }u(\mathbf {r} ),}أينu(ر){\displaystyle u(\mathbf {r} )}لها نفس الدورية التي تتميز بها البنية الذرية للبلورة، بحيثuك(x)=uك(x+نأ).{\displaystyle u_{\mathbf {k} }(\mathbf {x} )=u_{\mathbf {k} }(\mathbf {x} +\mathbf {n} \cdot \mathbf {a} ).}

وهناك طريقة ثانية ومكافئة لصياغة النظرية وهي التالية [ 3 ]

نظرية بلوخ لأي دالة موجية تحقق معادلة شرودنغر ولإزاحة متجه الشبكةأ{\displaystyle \mathbf {a} }يوجد متجه واحد على الأقلك{\displaystyle \mathbf {k} }بحيث: ψك(x+أ)=هـأناكأψك(x).{\displaystyle \psi _{\mathbf {k}}(\mathbf {x} +\mathbf {a} )=e^{i\mathbf {k} \cdot \mathbf {a} }\psi _{\mathbf {k} }(\mathbf {x} ).}

دليل

استخدام دورية الشبكة

إن نظرية بلوخ، كونها بيانًا حول دورية الشبكة، فإن جميع التناظرات في هذا البرهان مشفرة كتناظرات انتقالية للدالة الموجية نفسها.

البرهان باستخدام دورية الشبكة

المصدر: [ 4 ]

مقدمة: تناظرات البلورات، والشبكة البلورية، والشبكة المقلوبة

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

تحتوي البلورة ثلاثية الأبعاد على ثلاثة متجهات شبكية أولية a1 و a2 و a3 . إذا تم إزاحة البلورة بأي من هذه المتجهات الثلاثة، أو بمجموعة منها على النحو التالي :ن1أ1+ن2أ2+ن3أ3،{\displaystyle n_{1}\mathbf {a} _{1}+n_{2}\mathbf {a} _{2}+n_{3}\mathbf {a} _{3},} حيث n i هي ثلاثة أعداد صحيحة، فإن الذرات ستنتهي في نفس مجموعة المواقع التي بدأت بها.

من العناصر المفيدة الأخرى في البرهان متجهات الشبكة المقلوبة . وهي ثلاثة متجهات b₁ و b₂ و b₃ ( بوحدات معكوس الطول ) ، وتتميز بأن aᵢ · bᵢ = 2π ، بينما aᵢ · bⱼ = 0 عندما i j . ( للاطلاع على صيغة bᵢ ، انظر متجه الشبكة المقلوبة ) .

معضلة حول عوامل التحويل

يتركتي^ن1،ن2،ن3{\displaystyle {\hat {T}}_{n_{1},n_{2},n_{3}}}لنرمز إلى عامل إزاحة يُزيح كل دالة موجية بمقدار n₁a₁ + n₂a₂ + n₃a₃ (كما سبق، n₁ و n₂ أعداد صحيحة ) . الحقيقة التالية مفيدة لإثبات نظرية بلوخ :

اللمة - إذا كانت دالة الموجة ψ حالة ذاتية لجميع عوامل الإزاحة (في وقت واحد)، فإن ψ هي حالة بلوخ.

برهان اللمة

لنفترض أن لدينا دالة موجية ψ تمثل حالة ذاتية لجميع مؤثرات الإزاحة. وكحالة خاصة من ذلك، ψ(ر+أج)=ججψ(ر){\displaystyle \psi (\mathbf {r} +\mathbf {a} _{j})=C_{j}\psi (\mathbf {r} )} لـ j = 1، 2، 3 ، حيث C<sub> j</sub> هي ثلاثة أعداد ( القيم الذاتية ) لا تعتمد على r . من المفيد كتابة الأعداد C<sub> j</sub> بصيغة مختلفة، باختيار ثلاثة أعداد θ <sub>1</sub> ، θ <sub>2</sub> ، θ <sub>3 </sub> بحيث يكون e <sup> 2 </sup>πiθ<sub> j </sub> = C<sub> j </sub>. ψ(ر+أج)=هـ2πأناθجψ(ر){\displaystyle \psi (\mathbf {r} +\mathbf {a} _{j})=e^{2\pi i\theta _{j}}\psi (\mathbf {r} )}مرة أخرى ، θj هي ثلاثة أعداد لا تعتمد على r . عرّف k = θ1b1 + θ2b2 + θ3b3 ، حيث bj هي متجهات الشبكة المقلوبة ( انظر أعلاه ) . وأخيرًا ، عرّفu(ر)=هـ-أناكرψ(ر).{\displaystyle u(\mathbf {r} )=e^{-i\mathbf {k} \cdot \mathbf {r} }\psi (\mathbf {r} )\,.} ثم u(ر+أج)=هـ-أناك(ر+أج)ψ(ر+أج)=(هـ-أناكرهـ-أناكأج)(هـ2πأناθجψ(ر))=هـ-أناكرهـ-2πأناθجهـ2πأناθجψ(ر)=u(ر).// (}e^{-i\mathbf {k} \cdot \mathbf {r} }e^{-i\mathbf {k} \cdot \mathbf {a} _{j}}{\big )}{\big (}e^{2\pi i\theta _{j}}\psi (\mathbf {r} ){\big )}\\&=e^{-i\mathbf {ك} \cdot \mathbf {r} }e^{-2\pi i\theta _{j}}e^{2\pi i\theta _{j}}\psi (\mathbf {r} )\\&=u(\mathbf {r} ).\end{aligned}}} يثبت هذا أن u لها دورية الشبكة. بما أنψ(ر)=هـأناكرu(ر)،{\displaystyle \psi (\mathbf {r} )=e^{i\mathbf {k} \cdot \mathbf {r} }u(\mathbf {r} ),}وهذا يثبت أن الدولة هي دولة بلوخ.

وأخيراً، نحن جاهزون للبرهان الرئيسي لنظرية بلوخ والذي هو كما يلي.

كما سبق، دعتي^ن1،ن2،ن3{\displaystyle {\hat {T}}_{n_{1},n_{2},n_{3}}}لنرمز إلى مؤثر إزاحة يُزيح كل دالة موجية بمقدار n₁a₁ + n₂a₂ + n₃a₃ ، حيث nᵢ أعداد صحيحة . ولأن البلورة تتمتع بتناظر إزاحي، فإن هذا المؤثر يتبادل مع مؤثر هاميلتون . علاوة على ذلك، يتبادل كل مؤثر إزاحة من هذا النوع مع كل مؤثر آخر. لذلك، توجد قاعدة ذاتية متزامنة لمؤثر هاميلتون وكل دالة موجية ممكنة .تي^ن1،ن2،ن3{\displaystyle {\hat {T}}_{n_{1},n_{2},n_{3}}\!}المؤثر. هذا الأساس هو ما نبحث عنه. الدوال الموجية في هذا الأساس هي حالات طاقة ذاتية (لأنها حالات ذاتية للهاميلتوني)، وهي أيضًا حالات بلوخ (لأنها حالات ذاتية لمؤثرات الإزاحة؛ انظر اللمة أعلاه).

استخدام عوامل التشغيل

في هذا البرهان، يتم ترميز جميع التناظرات كخصائص تبادلية لمؤثرات الإزاحة.

البرهان باستخدام المعاملات

المصدر: [ 5 ]

نُعرّف عامل الإزاحة تي^نψ(ر)=ψ(ر+تين)=ψ(ر+ن1أ1+ن2أ2+ن3أ3)=ψ(ر+أن)// {a} _{1}+n_{2}\mathbf {a} _{2}+n_{3}\mathbf {a} _{3})\\&=\psi (\mathbf {r} +\mathbf {A} \mathbf {n} )\end{محاذاة}}} مع أ=[أ1أ2أ3]،ن=(ن1ن2ن3){\displaystyle \mathbf {A} ={\begin{bmatrix}\mathbf {a} _{1}&\mathbf {a} _{2}&\mathbf {a} _{3}\end{bmatrix}},\quad \mathbf {n} ={\begin{pmatrix}n_{1}\\n_{2}\\n_{3}\end{pmatrix}}} نستخدم فرضية الجهد الدوري المتوسط يو(x+تين)=يو(x){\displaystyle U(\mathbf {x} +\mathbf {T} _{\mathbf {n} })=U(\mathbf {x} )} وتقريب الإلكترون المستقل مع هاميلتوني ح^=ص^22م+يو(x){\displaystyle {\hat {H}}={\frac {{\hat {\mathbf {p} }}^{2}}{2m}}+U(\mathbf {x} )} بما أن الهاميلتوني ثابت بالنسبة للإزاحات، فإنه سيتبادل مع عامل الإزاحة. [ح^،تي^ن]=0{\displaystyle [{\hat {H}},{\hat {\mathbf {T} }}_{\mathbf {n} }]=0} ويجب أن يكون للمؤثرين مجموعة مشتركة من الدوال الذاتية. لذلك، نبدأ بالنظر إلى الدوال الذاتية لمؤثر الإزاحة: تي^نψ(x)=λنψ(x){\displaystyle {\hat {\mathbf {T} }}_{\mathbf {n} }\psi (\mathbf {x} )=\lambda _{\mathbf {n} }\psi (\mathbf {x} )} منحتي^ن{\displaystyle {\hat {\mathbf {T} }}_{\mathbf {n} }}هو عامل جمعي تي^ن1تي^ن2ψ(x)=ψ(x+أن1+أن2)=تي^ن1+ن2ψ(x){\displaystyle {\hat {\mathbf {T} }}_{\mathbf {n} _{1}}{\hat {\mathbf {T} }}_{\mathbf {n} _{2}}\psi (\mathbf {x} )=\psi (\mathbf {x} +\mathbf {A} \mathbf {n} _{1}+\mathbf {A} \mathbf {n} _{2})={\hat {\mathbf {T} }}_{\mathbf {n} _{1}+\mathbf {n} _{2}}\psi (\mathbf {x} )} إذا استبدلنا هنا معادلة القيم الذاتية وقسمنا كلا الطرفين علىψ(x){\displaystyle \psi (\mathbf {x} )}لدينا λن1λن2=λن1+ن2{\displaystyle \lambda _{\mathbf {n} _{1}}\lambda _{\mathbf {n} _{2}}=\lambda _{\mathbf {n} _{1}+\mathbf {n} _{2}}}

هذا صحيح بالنسبة لـ λن=هـsنأ{\displaystyle \lambda _{\mathbf {n} }=e^{s\mathbf {n} \cdot \mathbf {a} }} أينsج{\displaystyle s\in \mathbb {C} }إذا استخدمنا شرط التطبيع على خلية أولية واحدة بحجم V 1=V|ψ(x)|2دx=V|تي^نψ(x)|2دx=|λن|2V|ψ(x)|2دx{\displaystyle 1=\int _{V}|\psi (\mathbf {x} )|^{2}d\mathbf {x} =\int _{V}\left|{\hat {\mathbf {T} }}_{\mathbf {n} }\psi (\mathbf {x} )\right|^{2}d\mathbf {x} =|\lambda _{\mathbf {n} }|^{2}\int _{V}|\psi (\mathbf {x} )|^{2}d\mathbf {x} } وبالتالي 1=|λن|2{\displaystyle 1=|\lambda _{\mathbf {n} }|^{2}}وs=أناك{\displaystyle s=ik}أينكR{\displaystyle k\in \mathbb {R} }. أخيراً، تي^نψ(x)=ψ(x+نأ)=هـأناكنأψ(x)،{\displaystyle \mathbf {{\hat {T}}_{n}} \psi (\mathbf {x} )=\psi (\mathbf {x} +\mathbf {n} \cdot \mathbf {a} )=e^{ik\mathbf {n} \cdot \mathbf {a} }\psi (\mathbf {x} ),} وهذا ينطبق على موجة بلوخ، أي علىψك(x)=هـأناكxuك(x){\displaystyle \psi _{\mathbf {k} }(\mathbf {x} )=e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )}معuك(x)=uك(x+أن){\displaystyle u_{\mathbf {k} }(\mathbf {x} )=u_{\mathbf {k} }(\mathbf {x} +\mathbf {A} \mathbf {n} )}

استخدام نظرية الزمر

بغض النظر عن الجوانب التقنية لنظرية الزمر، فإن هذا البرهان مثير للاهتمام لأنه يوضح كيفية تعميم نظرية بلوخ للزمر التي لا تقتصر على الإزاحات فقط. يُطبق هذا عادةً على الزمر الفراغية التي تجمع بين إزاحة وزمرة نقطية ، ويُستخدم لحساب بنية النطاق والطيف والحرارة النوعية للبلورات، مع الأخذ في الاعتبار تناظر زمرة بلورية محددة مثل FCC أو BCC، وربما أساس إضافي . [ 6 ] : 365-367 [ 7 ] في هذا البرهان، يمكن أيضًا ملاحظة أهمية أن تكون الزمرة النقطية الإضافية مدفوعة بتناظر في الجهد الفعال، ولكنها يجب أن تتبادل مع الهاميلتوني.

البرهان باستخدام نظرية الشخصية [ 6 ] : 345-348

جميع التحويلات وحدوية وأبيلية . يمكن كتابة التحويلات بدلالة متجهات الوحدة .τ=أنا=13نأناأأنا{\displaystyle {\boldsymbol {\tau }}=\sum _{i=1}^{3}n_{i}\mathbf {a} _{i}} يمكننا اعتبار هؤلاء بمثابة مشغلي تنقلات τ^=τ^1τ^2τ^3{\displaystyle {\hat {\boldsymbol {\tau }}}={\hat {\boldsymbol {\tau }}}_{1}{\hat {\boldsymbol {\tau }}}_{2}{\hat {\boldsymbol {\tau }}}_{3}}أينτ^أنا=نأناأ^أنا{\displaystyle {\hat {\boldsymbol {\tau }}}_{i}=n_{i}{\hat {\mathbf {a} }}_{i}}

خاصية التبادلية لـτ^أنا{\displaystyle {\hat {\boldsymbol {\tau }}}_{i}}تُعطي المؤثرات ثلاث مجموعات فرعية دورية تبادلية (بافتراض إمكانية توليدها بواسطة عنصر واحد فقط) وهي مجموعات لانهائية، أحادية البعد، وأبيلية. جميع التمثيلات غير القابلة للاختزال للمجموعات الأبيلية أحادية البعد. [ 8 ]

بما أن هذه المجموعات أحادية البعد، فإن تمثيلها المصفوفي وخاصيتها متطابقان . الخاصية هي تمثيل المجموعة على الأعداد المركبة، أو هي أثر هذا التمثيل ، وهو في هذه الحالة مصفوفة أحادية البعد. جميع هذه المجموعات الفرعية، كونها دورية، لها خصائص تمثل جذور الوحدة . في الواقع، لها مولد واحد.γ{\displaystyle \gamma }الذي يجب أن يطيعγن=1{\displaystyle \gamma ^{n}=1}وبالتالي الشخصيةχ(γ)ن=1{\displaystyle \chi (\gamma )^{n}=1}لاحظ أن هذا واضح في حالة المجموعة الدورية المنتهية، ولكن في حالة المجموعة الدورية اللانهائية القابلة للعد (أي مجموعة الإزاحة هنا)، توجد نهاية لـن{\displaystyle n\to \infty }حيث تبقى الشخصية محدودة.

بما أن الحرف هو جذر الوحدة، فإنه يمكن كتابة الحرف لكل مجموعة فرعية على النحو التالي χك1(τ^1(ن1،أ1))=هـأناك1ن1أ1{\displaystyle \chi _{k_{1}}({\hat {\boldsymbol {\tau }}}_{1}(n_{1},a_{1}))=e^{ik_{1}n_{1}a_{1}}}

إذا أدخلنا شرط حدود بورن-فون كارمان على الجهد: V(ر+أناشمالأناأأنا)=V(ر+ل)=V(ر){\displaystyle V\left(\mathbf {r} +\sum _{i}N_{i}\mathbf {a} _{i}\right)=V(\mathbf {r} +\mathbf {L} )=V(\mathbf {r} )} حيث L هي دورية ماكروسكوبية في الاتجاهأ{\displaystyle \mathbf {a} }ويمكن اعتبار ذلك أيضًا مضاعفًا لـأأنا{\displaystyle a_{i}}أينل=أناشمالأناأأنا{\textstyle \mathbf {L} =\sum _{i}N_{i}\mathbf {a} _{i}}

هذا استبدال معادلة شرودنغر المستقلة عن الزمن بهاملتوني فعال بسيط ح^=-22م2+V(ر){\displaystyle {\hat {H}}=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V(\mathbf {r} )} يُحدث دورية مع الدالة الموجية: ψ(ر+أناشمالأناأأنا)=ψ(ر){\displaystyle \psi \left(\mathbf {r} +\sum _{i}N_{i}\mathbf {a} _{i}\right)=\psi (\mathbf {r} )}

ولكل بُعد، يوجد عامل إزاحة بفترة LP^ε|τأنا+لأنا=P^ε|τأنا{\displaystyle {\hat {P}}_{\varepsilon |\tau _{i}+L_{i}}={\hat {P}}_{\varepsilon |\tau _{i}}}

ومن هنا يمكننا أن نرى أن الحرف سيظل ثابتًا أيضًا عن طريق ترجمةلأنا{\displaystyle L_{i}}: هـأناك1ن1أ1=هـأناك1(ن1أ1+ل1){\displaystyle e^{ik_{1}n_{1}a_{1}}=e^{ik_{1}(n_{1}a_{1}+L_{1})}} ومن المعادلة الأخيرة نحصل على شرط دوري لكل بُعد: ك1ن1أ1=ك1(ن1أ1+ل1)-2πم1{\displaystyle k_{1}n_{1}a_{1}=k_{1}(n_{1}a_{1}+L_{1})-2\pi m_{1}} أينم1Z{\displaystyle m_{1}\in \mathbb {Z} }هو عدد صحيح وك1=2πم1ل1{\displaystyle k_{1}={\frac {2\pi m_{1}}{L_{1}}}}

متجه الموجةك1{\displaystyle k_{1}}تحديد التمثيل غير القابل للاختزال بنفس الطريقة التيم1{\displaystyle m_{1}}، ول1{\displaystyle L_{1}}يمثل طولًا دوريًا ماكروسكوبيًا للبلورة في اتجاهأ1{\displaystyle a_{1}}في هذا السياق، يعمل متجه الموجة كرقم كمي لمؤثر الإزاحة.

يمكننا تعميم هذا على ثلاثة أبعاد χك1(ن1،أ1)χك2(ن2،أ2)χك3(ن3،أ3)=هـأناكτ{\displaystyle \chi _{k_{1}}(n_{1},a_{1})\chi _{k_{2}}(n_{2},a_{2})\chi _{k_{3}}(n_{3},a_{3})=e^{i\mathbf {k} \cdot {\boldsymbol {\tau }}}} وتصبح الصيغة العامة للدالة الموجية كما يلي: P^Rψج=αψαχαج(R){\displaystyle {\hat {P}}_{R}\psi _{j}=\sum _{\alpha }\psi _{\alpha }\chi _{\alpha j}(R)} أي تخصيصها للترجمة P^ε|τψ(ر)=ψ(ر)هـأناكτ=ψ(ر+τ){\displaystyle {\hat {P}}_{\varepsilon |{\boldsymbol {\tau }}}\psi (\mathbf {r} )=\psi (\mathbf {r} )e^{i\mathbf {k} \cdot {\boldsymbol {\tau }}}=\psi (\mathbf {r} +{\boldsymbol {\tau }})} وقد أثبتنا نظرية بلوخ.

في النسخة المعممة من نظرية بلوخ، يتم تعميم تحويل فورييه، أي توسيع دالة الموجة، من تحويل فورييه منفصل ينطبق فقط على المجموعات الدورية، وبالتالي على الإزاحات، إلى توسيع مميز لدالة الموجة حيث يتم إعطاء المميزات من مجموعة النقاط المحدودة المحددة .

كما يمكن هنا أن نرى كيف يمكن التعامل مع الأحرف (باعتبارها ثوابت التمثيلات غير القابلة للاختزال) على أنها اللبنات الأساسية بدلاً من التمثيلات غير القابلة للاختزال نفسها. [ 9 ]

السرعة والكتلة الفعالة

إذا طبقنا معادلة شرودنغر غير المعتمدة على الزمن على دالة موجة بلوخ، فسنحصل على ح^كuك(ر)=[22م(-أنا+ك)2+يو(ر)]uك(ر)=εكuك(ر){\displaystyle {\hat {H}}_{\mathbf {k} }u_{\mathbf {k} }(\mathbf {r} )=\left[{\frac {\hbar ^{2}}{2m}}\left(-i\nabla +\mathbf {k} \right)^{2}+U(\mathbf {r} )\right]u_{\mathbf {k} }(\mathbf {r} )=\varepsilon _{\mathbf {k} }u_{\mathbf {k} }(\mathbf {r} )} مع الشروط الحدية uك(ر)=uك(ر+R){\displaystyle u_{\mathbf {k} }(\mathbf {r} )=u_{\mathbf {k} }(\mathbf {r} +\mathbf {R} )} بما أن هذا مُعرَّف في حجم محدود، فإننا نتوقع عائلة لانهائية من القيم الذاتية؛ هناك{\displaystyle {\mathbf {k} }}هو أحد معلمات الهاميلتوني، وبالتالي نصل إلى "عائلة متصلة" من القيم الذاتيةεن(ك){\displaystyle \varepsilon _{n}(\mathbf {k} )}يعتمد على المعامل المستمرك{\displaystyle {\mathbf {k} }}وبالتالي، فإن المفهوم الأساسي لبنية النطاق الإلكتروني هو...

البرهان [ 10 ]

هـك(هـأناكxuك(x))=[-22م2+يو(x)](هـأناكxuك(x)){\displaystyle E_{\mathbf {k} }\left(e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )\right)=\left[{\frac {-\hbar ^{2}}{2m}}\nabla ^{2}+U(\mathbf {x} )\right]\left(e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )\right)}

سنبقى مع هـكهـأناكxuك(x)=-22م(أناكهـأناكxuك(x)+هـأناكxuك(x))+يو(x)هـأناكxuك(x)هـكهـأناكxuك(x)=-22م(أناك(أناكهـأناكxuك(x)+هـأناكxuك(x))+أناكهـأناكxuك(x)+هـأناكx2uك(x))+يو(x)هـأناكxuك(x)هـكهـأناكxuك(x)=22م(ك2هـأناكxuك(x)-2أناكهـأناكxuك(x)-هـأناكx2uك(x))+يو(x)هـأناكxuك(x)هـكuك(x)=22م(-أنا+ك)2uك(x)+يو(x)uك(x){\displaystyle {\begin{aligned}E_{\mathbf {k} }e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )&={\frac {-\hbar ^{2}}{2m}}\nabla \cdot \left(i\mathbf {k} e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )+e^{i\mathbf {k} \cdot \mathbf {x} }\nabla u_{\mathbf {k} }(\mathbf {x} )\right)+U(\mathbf {x} )e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )\\[1.2ex]E_{\mathbf {k} }e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )&={\frac {-\hbar ^{2}}{2m}}\left(i\mathbf {k} \cdot \left(i\mathbf {k} e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )+e^{i\mathbf {k} \cdot \mathbf {x} }\nabla u_{\mathbf {k} }(\mathbf {x} )\right)+i\mathbf {k} \cdot e^{i\mathbf {k} \cdot \mathbf {x} }\nabla u_{\mathbf {k} }(\mathbf {x} )+e^{i\mathbf {k} \cdot \mathbf {x} }\nabla ^{2}u_{\mathbf {k} }(\mathbf {x} )\right)+U(\mathbf {x} )e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )\\[1.2ex]E_{\mathbf {k} }e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )&={\frac {\hbar ^{2}}{2m}}\left(\mathbf {k} ^{2}e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )-2i\mathbf {k} \cdot e^{i\mathbf {k} \cdot \mathbf {x} }\nabla u_{\mathbf {k} }(\mathbf {x} )-e^{i\mathbf {k} \cdot \mathbf {x} }\nabla ^{2}u_{\mathbf {k} }(\mathbf {x} )\right)+U(\mathbf {x} )e^{i\mathbf {k} \cdot \mathbf {x} }u_{\mathbf {k} }(\mathbf {x} )\\[1.2ex]E_{\mathbf {k} }u_{\mathbf {k} }(\mathbf {x} )&={\frac {\hbar ^{2}}{2m}}\left(-i\nabla +\mathbf {k} \right)^{2}u_{\mathbf {k} }(\mathbf {x} )+U(\mathbf {x} )u_{\mathbf {k} }(\mathbf {x} )\end{aligned}}}

يوضح هذا كيف يمكن اعتبار الزخم الفعال مكونًا من جزأين، ص^فعال=-أنا+ك،{\displaystyle {\hat {\mathbf {p} }}_{\text{eff}}=-i\hbar \nabla +\hbar \mathbf {k} ,} الزخم القياسي-أنا{\displaystyle -i\hbar \nabla }وزخم بلوريك{\displaystyle \hbar \mathbf {k} }وبشكل أدق، فإن الزخم البلوري ليس زخمًا ولكنه يمثل الزخم بنفس طريقة الزخم الكهرومغناطيسي في الحد الأدنى من الاقتران ، وكجزء من تحويل قانوني للزخم.

يمكننا اشتقاق السرعة الفعالة

متوسط ​​سرعة إلكترون بلوخ

εنك=2مدرψنك*(-أنا)ψنك=مص^=v^{\displaystyle {\frac {\partial \varepsilon _{n}}{\partial \mathbf {k} }}={\frac {\hbar ^{2}}{m}}\int d\mathbf {r} \,\psi _{n\mathbf {k} }^{*}(-i\nabla )\psi _{n\mathbf {k} }={\frac {\hbar }{m}}\langle {\hat {\mathbf {p} }}\rangle =\hbar \langle {\hat {\mathbf {v} }}\rangle }

البرهان [ 11 ]

نقوم بتقييم المشتقاتεنك{\displaystyle {\frac {\partial \varepsilon _{n}}{\partial \mathbf {k} }}}و2εن(ك)كأناكج{\displaystyle {\frac {\partial ^{2}\varepsilon _{n}(\mathbf {k} )}{\partial k_{i}\partial k_{j}}}} بافتراض أنها معاملات التوسع التالي في q حيث تعتبر q صغيرة بالنسبة إلى kεن(ك+q)=εن(ك)+أناεنكأناqأنا+12أناج2εنكأناكجqأناqج+يا(q3){\displaystyle \varepsilon _{n}(\mathbf {k} +\mathbf {q} )=\varepsilon _{n}(\mathbf {k} )+\sum _{i}{\frac {\partial \varepsilon _{n}}{\partial k_{i}}}q_{i}+{\frac {1}{2}}\sum _{ij}{\frac {\partial ^{2}\varepsilon _{n}}{\partial k_{i}\partial k_{j}}}q_{i}q_{j}+O(q^{3})} منحεن(ك+q){\displaystyle \varepsilon _{n}(\mathbf {k} +\mathbf {q} )}هي القيم الذاتية لـح^ك+q{\displaystyle {\hat {H}}_{\mathbf {k} +\mathbf {q} }} يمكننا النظر في مسألة الاضطراب التالية في q: ح^ك+q=ح^ك+2مq(-أنا+ك)+22مq2{\displaystyle {\hat {H}}_{\mathbf {k} +\mathbf {q} }={\hat {H}}_{\mathbf {k} }+{\frac {\hbar ^{2}}{m}}\mathbf {q} \cdot (-i\nabla +\mathbf {k} )+{\frac {\hbar ^{2}}{2m}}q^{2}} تنص نظرية الاضطراب من الدرجة الثانية على أن هـن=هـن0+درψن*V^ψن+نن|درψن*V^ψن|2هـن0-هـن0+...{\displaystyle E_{n}=E_{n}^{0}+\int d\mathbf {r} \,\psi _{n}^{*}{\hat {V}}\psi _{n}+\sum _{n'\neq n}{\frac {|\int d\mathbf {r} \,\psi _{n}^{*}{\hat {V}}\psi _{n}|^{2}}{E_{n}^{0}-E_{n'}^{0}}}+...} لحساب الرتبة الخطية في qأناεنكأناqأنا=أنادرuنك*2م(-أنا+ك)أناqأناuنك{\displaystyle \sum _{i}{\frac {\partial \varepsilon _{n}}{\partial k_{i}}}q_{i}=\sum _{i}\int d\mathbf {r} \,u_{n\mathbf {k} }^{*}{\frac {\hbar ^{2}}{m}}(-i\nabla +\mathbf {k} )_{i}q_{i}u_{n\mathbf {k} }} حيث تكون عمليات التكامل على خلية أولية أو على البلورة بأكملها، معطى إذا كان التكامل درuنك*uنك{\displaystyle \int d\mathbf {r} \,u_{n\mathbf {k} }^{*}u_{n\mathbf {k} }} يتم توحيدها عبر الخلية أو البلورة.

يمكننا التبسيط بالنسبة لـ q للحصول على εنك=2مدرuنك*(-أنا+ك)uنك{\displaystyle {\frac {\partial \varepsilon _{n}}{\partial \mathbf {k} }}={\frac {\hbar ^{2}}{m}}\int d\mathbf {r} \,u_{n\mathbf {k} }^{*}(-i\nabla +\mathbf {k} )u_{n\mathbf {k} }} ويمكننا إعادة إدخال الدوال الموجية الكاملة εنك=2مدرψنك*(-أنا)ψنك{\displaystyle {\frac {\partial \varepsilon _{n}}{\partial \mathbf {k} }}={\frac {\hbar ^{2}}{m}}\int d\mathbf {r} \,\psi _{n\mathbf {k} }^{*}(-i\nabla )\psi _{n\mathbf {k} }}

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

نظرية الكتلة الفعالة

2εن(ك)كأناكج=2مدلتاأناج+(2م)2نننك|-أناأنا|نكنك|-أناج|نك+نك|-أناج|نكنك|-أناأنا|نكεن(ك)-εن(ك){\displaystyle {\frac {\partial ^{2}\varepsilon _{n}(\mathbf {k} )}{\partial k_{i}\partial k_{j}}}={\frac {\hbar ^{2}}{m}}\delta _{ij}+\left({\frac {\hbar ^{2}}{m}}\right)^{2}\sum _{n'\neq n}{\frac {\langle n\mathbf {k} |-i\nabla _{i}|n'\mathbf {k} \rangle \langle n'\mathbf {k} |-i\nabla _{j}|n\mathbf {k} \rangle +\langle n\mathbf {k} |-i\nabla _{j}|n'\mathbf {k} \rangle \langle n'\mathbf {k} |-i\nabla _{i}|n\mathbf {k} \rangle }{\varepsilon _{n}(\mathbf {k} )-\varepsilon _{n'}(\mathbf {k} )}}}

البرهان [ 11 ]

الحد من الدرجة الثانية 12أناج2εنكأناكجqأناqج=22مq2+نن|درuنك*2مq(-أنا+ك)uنك|2εنك-εنك{\displaystyle {\frac {1}{2}}\sum _{ij}{\frac {\partial ^{2}\varepsilon _{n}}{\partial k_{i}\partial k_{j}}}q_{i}q_{j}={\frac {\hbar ^{2}}{2m}}q^{2}+\sum _{n'\neq n}{\frac {|\int d\mathbf {r} \,u_{n\mathbf {k} }^{*}{\frac {\hbar ^{2}}{m}}\mathbf {q} \cdot (-i\nabla +\mathbf {k} )u_{n'\mathbf {k} }|^{2}}{\varepsilon _{n\mathbf {k} }-\varepsilon _{n'\mathbf {k} }}}} مرة أخرى معψنك=|نك=هـأناكxuنك{\displaystyle \psi _{n\mathbf {k} }=|n\mathbf {k} \rangle =e^{i\mathbf {k} \mathbf {x} }u_{n\mathbf {k} }}12أناج2εنكأناكجqأناqج=22مq2+نن|نك|2مq(-أنا)|نك|2εنك-εنك{\displaystyle {\frac {1}{2}}\sum _{ij}{\frac {\partial ^{2}\varepsilon _{n}}{\partial k_{i}\partial k_{j}}}q_{i}q_{j}={\frac {\hbar ^{2}}{2m}}q^{2}+\sum _{n'\neq n}{\frac {|\langle n\mathbf {k} |{\frac {\hbar ^{2}}{m}}\mathbf {q} \cdot (-i\nabla )|n'\mathbf {k} \rangle |^{2}}{\varepsilon _{n\mathbf {k} }-\varepsilon _{n'\mathbf {k} }}}} الإزالةqأنا{\displaystyle q_{i}}وqج{\displaystyle q_{j}}لدينا النظرية 2εن(ك)كأناكج=2مدلتاأناج+(2م)2نننك|-أناأنا|نكنك|-أناج|نك+نك|-أناج|نكنك|-أناأنا|نكεن(ك)-εن(ك){\displaystyle {\frac {\partial ^{2}\varepsilon _{n}(\mathbf {k} )}{\partial k_{i}\partial k_{j}}}={\frac {\hbar ^{2}}{m}}\delta _{ij}+\left({\frac {\hbar ^{2}}{m}}\right)^{2}\sum _{n'\neq n}{\frac {\langle n\mathbf {k} |-i\nabla _{i}|n'\mathbf {k} \rangle \langle n'\mathbf {k} |-i\nabla _{j}|n\mathbf {k} \rangle +\langle n\mathbf {k} |-i\nabla _{j}|n'\mathbf {k} \rangle \langle n'\mathbf {k} |-i\nabla _{i}|n\mathbf {k} \rangle }{\varepsilon _{n}(\mathbf {k} )-\varepsilon _{n'}(\mathbf {k} )}}}

الكمية الموجودة على اليمين مضروبة في عامل12{\displaystyle {\frac {1}{\hbar ^{2}}}}يُطلق عليه اسم موتر الكتلة الفعالةم(ك){\displaystyle \mathbf {M} (\mathbf {k} )}[ 12 ] ويمكننا استخدامها لكتابة معادلة شبه كلاسيكية لحاملالشحنةفي نطاق [ 13 ]

معادلة الحركة شبه الكلاسيكية من الدرجة الثانية لحاملة الشحنة في نطاق

م(ك)أ=هـ(هـ+v(ك)×ب){\displaystyle \mathbf {M} (\mathbf {k} )\mathbf {a} =\mp e\left(\mathbf {E} +\mathbf {v} (\mathbf {k} )\times \mathbf {B} \right)}

أينأ{\displaystyle \mathbf {a} }يمثل تسارعًا . هذه المعادلة مماثلة لتقريب موجة دي برولي [ 14 ]

معادلة الحركة شبه الكلاسيكية من الدرجة الأولى للإلكترون في نطاق طاقة

ك˙=-هـ(هـ+v×ب){\displaystyle \hbar {\dot {k}}=-e\left(\mathbf {E} +\mathbf {v} \times \mathbf {B} \right)}

كتفسير بديهي، فإن المعادلتين السابقتين تشبهان من الناحية الشكلية وتوجدان في تشابه شبه كلاسيكي مع قانون نيوتن الثاني للإلكترون في قوة لورنتز خارجية .

تحذير رياضي

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

يمكن إجراء عمليات تحليل الدوال المستمرة الدورية - على غرار بلوخ - في فضاءات الدوال المحدودة أو الدوال المستمرة المحدودة، ولكن ليس في فضاءات الدوال القابلة للتكامل التربيعي على فضاء x الكامل، وهو ما يمثل إعداد فضاء هيلبرت المطلوب لميكانيكا الكم.

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

طُوِّر مفهوم حالة بلوخ على يد فيليكس بلوخ عام 1928 [ 17 ] لوصف توصيل الإلكترونات في المواد الصلبة البلورية. مع ذلك، اكتُشفت الرياضيات الأساسية نفسها بشكل مستقل عدة مرات: على يد جورج ويليام هيل (1877) [ 18 ] ، وجاستون فلوكي (1883) [ 19 ] ، وألكسندر ليابونوف (1892) [ 20 ] . ونتيجةً لذلك، شاع استخدام تسميات متنوعة: عند تطبيقها على المعادلات التفاضلية العادية ، تُسمى نظرية فلوكي (أو أحيانًا نظرية ليابونوف-فلوكي ). الصيغة العامة لمعادلة الجهد الدوري أحادي البعد هي معادلة هيل : [ 21 ]د2yدت2+و(ت)y=0،{\displaystyle {\frac {d^{2}y}{dt^{2}}}+f(t)y=0,} حيث f ( t ) هو جهد دوري. تشمل المعادلات الدورية أحادية البعد نموذج كرونيغ-بيني ومعادلة ماثيو .

رياضياً، تُفسَّر العديد من النظريات المشابهة لنظرية بلوخ، على سبيل المثال، من حيث الخصائص الوحدوية لمجموعة الشبكة، وتُطبَّق على الهندسة الطيفية . [ 22 ] [ 23 ] [ 24 ]

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

انظر أيضاً

مراجع

  1. ^ بلوخ، ف. (1929). Über die quantenmechanik der electronen in kristallgittern. Zeitschrift für physik, 52(7)، 555-600.
  2. كيتل، تشارلز (1996). مقدمة في فيزياء الحالة الصلبة . نيويورك: وايلي. ISBN 0-471-14286-7.
  3. زيمان، ج.م. (1972). مبادئ نظرية المواد الصلبة ( الطبعة الثانية). مطبعة جامعة كامبريدج. الصفحات 17-20 . ISBN   0521297338.
  4. ^ أشكروفت وميرمين 1976 ، ص. 134
  5. ^ أشكروفت وميرمين 1976 ، ص. 137
  6. 1 2 دريسلهاوس، إم إس (2002). "تطبيقات نظرية الزمر على فيزياء المواد الصلبة" (ملف PDF) . معهد ماساتشوستس للتكنولوجيا . مؤرشف (ملف PDF) من الأصل في 1 نوفمبر 2019. تم الاطلاع عليه في 12 سبتمبر 2020 .
  7. الطيف الاهتزازي والحرارة النوعية لبلورة مكعبة مركزية الوجوه، روبرت ب. لايتون
  8. روي، ريكي (2 مايو 2010). "نظرية التمثيل" (ملف PDF) . جامعة بيوجت ساوند.
  9. تمثيلات المجموعات والتحليل التوافقي من أويلر إلى لانجلاندز، الجزء الثاني
  10. ^ أشكروفت وميرمين 1976 ، ص. 140
  11. 1 2 أشكروفت وميرمين 1976 ، ص. 765 الملحق هـ
  12. ^ أشكروفت وميرمين 1976 ، ص. 228 
  13. ^ أشكروفت وميرمين 1976 ، ص. 229 
  14. ^ أشكروفت وميرمين 1976 ، ص. 227 
  15. ريد، مايكل؛ سيمون، باري (1978). أساليب الفيزياء الرياضية الحديثة - الجزء الرابع: تحليل المؤثرات . إلسيفير، أكاديميك برس. ص 279 وما بعدها. ISBN  978-0-08-057045-7.
  16. أفرون، ج.؛ جروسمان، أ.؛ رودريغيز، ر. (1974). "الهاميلتونيان في نظرية الإلكترون الواحد للمواد الصلبة. الجزء الأول". تقارير في الفيزياء الرياضية . 5 (1): 113-120 . doi : 10.1016/0034-4877(74)90020-2 .
  17. ^ فيليكس بلوخ (1928). "Über die Quantenmechanik der Elektronen in Kristallgittern". Zeitschrift für Physik (باللغة الألمانية). 52 ( 7– 8): 555–600 . بيب كود : 1929ZPhy...52..555B . دوى : 10.1007/BF01339455 . S2CID 120668259 . 
  18. جورج ويليام هيل (1886). "حول جزء حركة الحضيض القمري الذي هو دالة للحركات المتوسطة للشمس والقمر" . Acta Math . 8 : 1–36 . doi : 10.1007/BF02417081 .نُشر هذا العمل في البداية وتم توزيعه بشكل خاص في عام 1877.
  19. ^ غاستون فلوكيت (1883). "حول المعادلات التفاضلية الخطية للمعاملات الدورية" . الحوليات العلمية للمدرسة العليا للأساتذة . 12 : 47 – 88. دوى : 10.24033/asens.220 .
  20. ألكسندر ميهايلوفيتش ليابونوف (1992). المشكلة العامة لاستقرار الحركة . لندن: تايلور وفرانسيس.ترجمة أ. ت. فولر من الترجمة الفرنسية لإدوارد دافو (1907) للأطروحة الروسية الأصلية (1892).
  21. ماغنوس، دبليو ؛ وينكلر، إس (2004). معادلة هيل . كوريير دوفر. ص 11. ISBN  0-486-49565-5.
  22. كوشمنت، ب. (1982)، نظرية فلوكي للمعادلات التفاضلية الجزئية ، مجلة الرياضيات الروسية، 37، 1-60
  23. كاتسودا، أ.؛ سونادا، ت. (1987). "علم التماثل والجيوديسيات المغلقة في سطح ريمان مضغوط". المجلة الأمريكية للرياضيات 110 ( 1): 145-156 . doi : 10.2307/2374542 . JSTOR 2374542 . 
  24. كوتاني م؛ سونادا ت. (2000). "خرائط ألبانيز وسلوك تقاربي طويل الأمد غير قطري لنواة الحرارة". مجلة الاتصالات في الفيزياء الرياضية . 209 (3): 633-670 . رمز Bibcode : 2000CMaPh.209..633K . doi : 10.1007/s002200050033 . S2CID 121065949 . 

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