نهج DFT-1/2 دقيق لحالات الشوائب الضحلّة: حساب كفء لطاقات ربط المانحين في السليكون

لماذا تهم التعديلات الطفيفة في السيليكون

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

الذرات التي تمنح إلكترونات إضافية

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

لماذا تفشل الحسابات القياسية

نماذج الحاسوب المبنية على نظرية الدوال الكثافة (DFT) هي حجر الأساس في تصميم المواد، لكنها تميل إلى التقليل من مدى تمركز الإلكترونات وإلى وضع حواف نطاقات الطاقة في أشباه الموصلات في مواضع خاطئة. بالنسبة للمانحين الضحلين، الذين تمتد سحب إلكتروناتهم عبر عشرات الذرات، يعني هذا أن DFT عادةً ما تتنبأ بطاقة ربط أقل بكثير من الواقع. يمكن للطرائق المتقدمة، مثل الدوال الهجينة وحسابات GW، أن تصلح هذه المشاكل لكنها باهظة التكلفة الحسابية، لا سيما عندما تكون صناديق المحاكاة كبيرة وتضم آلاف الذرات لاحتواء حالة المانح الممتدة. سبق أن اضطرّت نهج «التناوب» إلى مزج مستويات نظرية مختلفة في أحجام خلايا مختلفة ثم تجميع النتائج، مما جعل سير العمل معقَّداً ومعتمداً على النظام المدروس.

تصحيح بسيط بعائد كبير

يبني المؤلفون على تقنية تُدعى DFT-1/2، التي تضيف تصحيحاً تقريبيّاً للطاقة الذاتية مباشرة إلى DFT القياسي. عملياً، يقومون بتعديل طفيف في الكمون الفعّال لذرات محددة عبر إزالة مفهومية لنصف إلكترون من مدارات ذرية معينة. أولاً، يطبقون هذا التصحيح على السيليكون البلّوري بحيث يطابق الفجوة الطاقية المحسوبة القياسات التجريبية بشكل أفضل ويأمن مرجعاً موثوقاً لحالة النطاق التوصيلي. بعد ذلك، يفحصون الطابع الإلكتروني لحالة المانح ويجدون أنها، لجميع مشوّهات المجموعة V، تُهيمن عليها المدار s للذرة المانحة. يطبقون تصحيح نصف الإلكترون المخصص لذلك المدار ويضبطون بدقّة نصف قطر قطع وحيد لتعظيم الفاصل بين مستوى المانح وأقرب حالة موصِلة فارغة. والأهم أن هذا التصحيح المُحسن يبقى صالحاً عند تكبير صندوق المحاكاة، فيمكن إعادة استخدامه في فوق خلايا تحوي آلاف الذرات.

مدى فعالية الطريقة

بهذا التصحيح ذي الخطوتين — أولاً للمادة المضيفة ثم للمانح نفسه — تنتج الطريقة طاقات ربط للمانحين تطابق قيماً تجريبية عن كثب. بالنسبة للزرنيخ في السيليكون، يختلف الطاقة المتوقعة عن التجربة بمقدار 0.3 ميلي إلكترون فولت فقط، اتفاق شبه مثالي ويقارن بنتائج دوال هجينة أغلى بكثير. بالنسبة للأنتبوم والفسفور، تكون الأخطاء نحو 5 و8 ميلي إلكترون فولت على التوالي، تحسّن كبير على DFT غير المصححة. بالنسبة للبزموت، وهو مشوٍّ ثقيل جداً، يدرج المؤلفون أيضاً ارتباط الدوران-المدار (spin-orbit coupling)، وهو أثر نسبي يعيد تشكيل مستويات الطاقة قليلاً. هذا يخفض طاقة الربط المحسوبة لتقع ضمن نحو 5 ميلي إلكترون فولت من القيم التجريبية ويبرز فيزيائياً عناصر أغفلتها طرق سابقة وأكثر تكلفة. ولإظهار أن النهج لا يقتصر على السيليكون، يطبقون سير العمل نفسه بنجاح على مانح هيدروجين في أكسيد الزنك، مع إعادة إنتاج طاقات ربط مقاسة ضمن بضعة ميلي إلكترون فولت.

أداة عملية لتصميم شرائح المستقبل

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

الاستشهاد: Claes, J., Partoens, B., Lamoen, D. et al. An accurate DFT-1/2 approach for shallow defect states: efficient calculation of donor binding energies in silicon. npj Comput Mater 12, 153 (2026). https://doi.org/10.1038/s41524-026-02003-2

الكلمات المفتاحية: المانحون الضحلون, مشوّات السيليكون, نظرية الدوال الكثافة, المواد الكمومية, عيوب أشباه الموصلات