First-Principles Study of Fe Adsorption and Its Effects on the Mechanical and Electrical Properties of Monolayer and Bilayer Biphenylene Networks

ملخص تقني: دراسة من المبادئ الأولى لامتزاز الحديد (Fe) وتأثيراته على الخصائص الميكانيكية والكهربائية لشبكات البيفينيلين (Biphenylene) أحادية وطبقة ثنائية

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

المنهجية
تعتمد هذه الدراسة على حسابات المبادئ الأولى القائمة على نظرية وظيفية الكثافة (DFT) باستخدام حزمة محاكاة "فيينا" للبدايات من الأول (VASP). يتم التعامل مع طاقة الارتباط التبادلي ضمن تقريب التدرج المعمم (PBE functional)، مع تصحيح تفاعلات فان دير فالس بعيدة المدى عبر طريقة DFT-D4. كما تُجرى حسابات مستقطبة مغناطيسياً (Spin-polarized) مع حد قطع للطاقة الحركية للموجة المستوية يبلغ 500 إلكترون فولت.

يركز البحث على خلايا فائقة (supercells) بحجم 2×2 لشبكات BPN أحادية وطبقة ثنائية. ولتحديد الاستقرار الديناميكي الحراري، أُجريت تحليلات "الغلاف المحدب" (convex hull) باستخدام متوسط طاقات الامتزاز ($E_{ad}$) لأحجام مختلفة من عناقيد الحديد ($n$). كما أُجريت عمليات تحسين هيكلية حتى انخفضت القوى إلى أقل من $10^{-3}$ إلكترون فولت/أنجستروم. ويتم تقييم الخصائص الميكانيكية باستخدام طريقة "الطاقة-الإجهاد" لحساب الثوابت المرنة ($C_{ij}$)، والتي تُشتق منها معاملات يونغ ومعاملات بواسون المعتمدة على الاتجاه. ويتم حساب التوصيل الكهربائي عند 300 كلفن باستخدام نظرية بولتزمان للنقل ضمن تقريب زمن الاسترخاء الثابت، بافتراض زمن استرخاء حامل الشحنة ($\tau$) قدره $10^{-15}$ ثانية.

المساهمات الرئيسية والنتائج

• الديناميكا الحرارية للامتزاز وتفضيل الموقع:

  • BPN أحادية الطبقة: يتناقص متوسط طاقة الامتزاز (يصبح أكثر سالبية) مع زيادة تغطية الحديد، مما يشير إلى تقوية التفاعلات بين الحديد والركيزة عند مستويات التغطية المتوسطة. تم تحديد التكوين الأكثر استقراراً عند نسبة حديد/كربون تبلغ 50% (12 ذرة حديد لكل خلية فائقة 2×2). وتعتمد تفضيلات المواقع على التغطية: فعند التغطية المنخفضة، يفضل الحديد مركز الحلقات سداسية الأعضاء، بينما يصبح مركز الحلقة رباعية الأعضاء هو المفضل عند التغطية المنخفضة جداً (مثل خلية فائقة 1/66).
  • BPN ثنائية الطبقة: يُفضل الامتزاز بين الطبقات طاقياً على الامتزاز السطحي. الموقع الأكثر استقراراً لذرة حديد واحدة هو مركز الحلقة رباعية الأعضاء الموجودة بين الطبقات. ويظل أقصى تحميل مستقر هو 12 ذرة حديد لكل خلية فائقة 2×2. ومن الجدير بالذكر أن امتزاز الحديد في الفراغ بين الطبقات يلغي الإزاحة التراكمية الأصلية للطبقة الثنائية، مما يؤدي إلى استقرار تكوين أكثر تماثلاً.

• الخصائص الميكانيكية:

  • الاستقرار داخل المستوي: تُظهر شبكات BPN الأصلية (أحادية وثنائية الطبقة) معاملات يونغ وقص عالية داخل المستوي، مما يؤكد الاستقرار الميكانيكي الممتاز. ويُظهر النظام ثنائي الطبقة معامل يونغ يعادل ضعف معامل الطبقة الأحادية تقريباً في اتجاه x، و1.8 مرة في اتجاه y. كما تُظهر كلا الأنظمة تباينًا ملحوظًا في الخصائص داخل المستوي، حيث يكون اتجاه y أكثر صلابة.
  • تأثير الحديد على ميكانيكا داخل المستوي: لامتزاز الحديد تأثير طفيف على الخصائص الميكانيكية داخل المستوي لكل من BPN أحادية وثنائية الطبقة عند التغطية المنخفضة والمتوسطة، مما يشير إلى أن الصلابة داخل المستوي محكومة بإطار الكربون الأصلي. وفقط عند تركيزات الحديد العالية (أكثر من 7 ذرات في خلية 2×2) يزداد معامل يونغ داخل المستوي بشكل ملحوظ.
  • الميكانيكا خارج المستوي: حُسب الثابت المرن خارج المستوي ($C_{33}$) لشبكة BPN ثنائية الطبقة الأصلية بـ 24.59 جيجا باسكال، مما يشير إلى تفاعلات ضعيفة بين الطبقات. ومع ذلك، فإن امتزاز الحديد بين الطبقات يعزز هذه الخاصية بشكل كبير؛ فعند امتزاز 12 ذرة حديد (نسبة حديد/كربون 25%)، يرتفع $C_{33}$ إلى 515.63 جيجا باسكال، مما يثبت أن الحديد الممتز بين الطبقات يمكنه بفعالية تعزيز الاقتران بين الطبقات والصلابة خارج المستوي.

• الخصائص الكهربائية:

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

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