DOI: https://doi.org/10.1038/s41524-024-01503-3
PMID: https://pubmed.ncbi.nlm.nih.gov/39872025
تاريخ النشر: 2025-01-24
المؤلف: Sahar Pakdel وآخرون
الموضوع الرئيسي: المواد ثنائية الأبعاد والتطبيقات
طرق
في هذه الدراسة، يستكشف المؤلفون الخصائص الإلكترونية والمغناطيسية للمواد ثنائية الأبعاد (2D) التي تحتوي على معادن الانتقال 3d (V، Cr، Mn، Fe، Co، Ni، Cu) باستخدام نهج حسابي عالي الإنتاجية. تم اختيار المواد من قاعدة بيانات C2DB، مما يضمن أن كل طبقة أحادية لها طاقة فوق الحدود المحدبة أقل من 0.2 eV/ذرة، كما تحددها دالة Perdew-Burke-Ernzerhof (PBE). يضمن هذا المعيار أن المواد مستقرة كيميائيًا وأن النتائج ليست مشوهة بهياكل غير واقعية. بالإضافة إلى ذلك، تم حساب الخصائص فقط للطبقات الأحادية حيث أسفرت طريقة PBE+U عن فجوة نطاق نهائية، حيث أن تصحيحات Hubbard U تنطبق فقط على الأنظمة غير المعدنية.
استخدمت المنهجية الحسابية إطار عمل وصفات المحاكاة الذرية (ASR) بلغة بايثون واستفادت من كود هيكل GPAW الإلكتروني، مع قطع موجي مسطح يبلغ 800 eV وكثافة نقاط k تبلغ 12 Å. تضمنت حسابات DFT+U معلمة U فعالة، تعرف على أنها \( U_{\text{eff}} = U – J \)، لأخذ تفاعلات كولومب وتبادل الموقع في الاعتبار. يتم إعطاء تعبير الطاقة الكلية لطريقة DFT+U بواسطة \( E_{\text{DFT+U}} = E_{\text{DFT}} + \sum_{a,i} \frac{1}{2} U_{\text{eff},i} \text{Tr}(\rho_{a,i} – \rho_{a,i} \rho_{a,i}) \)، حيث يكون المجموع على الذرات \( a \) وأوربتالاتها 3d المقابلة. قام المؤلفون بتعيين \( U_{\text{eff}} \) إلى 4 eV لهذه الأوربتالات وتأكدوا من أن جميع الهياكل تم استرخاؤها إلى قوة قصوى على الذرات أقل من 0.01 eV/Å.
نتائج
في قسم النتائج، يستكشف المؤلفون أولاً اعتماد خصائص مختلفة على معلمة Hubbard U لمجموعة مختارة من المواد. يحددون قيمة تصحيح U تبلغ 4 eV للحسابات اللاحقة، مسترشدين بمبدأ شفرة أوكام. ومع ذلك، يجدون أن تضمين هذه القيمة U يؤدي إلى توافق أقل مع البيانات التجريبية فيما يتعلق بثوابت الشبكة في الطبقات الأحادية. نتيجة لذلك، يتم إجراء جميع حسابات الخصائص اللاحقة باستخدام الطبقات الأحادية المسترخية PBE.
بالإضافة إلى ذلك، يفحص المؤلفون تأثير Hubbard U على فجوة النطاق والخصائص المغناطيسية عبر جميع الطبقات الأحادية المدروسة. تهدف هذه التحليل الشامل إلى توضيح دور U في التأثير على الخصائص الإلكترونية والمغناطيسية الرئيسية، مما يوفر رؤى حول خصائص المواد قيد النظر.
نقاش
في هذه الدراسة، يستكشف باكدل وأولسن وثيجسن بشكل منهجي تأثير تصحيحات Hubbard U على الهيكل الإلكتروني لـ 638 مادة ثنائية الأبعاد (2D) تحتوي على معادن الانتقال 3d، باستخدام نظرية الوظائف الكثافة (DFT) مع كل من تقريبات PBE و PBE+U. تشير النتائج إلى أنه بينما يؤدي تضمين تصحيح U عمومًا إلى زيادة فجوات النطاق—مما يؤدي إلى انتقال من المعدن إلى العازل لـ 134 مادة (21%)—فإنه يؤثر سلبًا على دقة ثوابت الشبكة، مما يؤدي إلى تفضيل هياكل PBE في التقييمات اللاحقة. من الجدير بالذكر أن العزوم المغناطيسية تظهر اعتمادًا ضئيلًا على U، بينما يتم تقليل كل من اقتران التبادل المغناطيسي ومعلمات التباين المغناطيسي بشكل منهجي، مما يشير إلى أن تصحيح Hubbard U قد يخفض درجات حرارة كوري المتوقعة في هذه المواد.
يكشف التحليل أن تصحيح U يؤثر بشكل كبير على الخصائص الإلكترونية، حيث تزداد فجوات النطاق عادةً بمقدار 0-0.5 eV، على الرغم من أن بعض المواد تظهر سلوكيات أكثر تعقيدًا. تؤكد الدراسة على أن الخصائص المغناطيسية تظل غير متغيرة إلى حد كبير بالنسبة لـ U، حيث تظل العزوم المغناطيسية قيمًا صحيحة للحالات العازلة. علاوة على ذلك، تسلط الأبحاث الضوء على التحديات المتمثلة في وصف المواد ثنائية الأبعاد بدقة مع الإلكترونات المحلية، حيث تكافح طرق DFT التقليدية مع الأخطاء الناتجة عن التفاعل الذاتي والتفكيك. توفر النتائج، التي تم تجميعها في قاعدة بيانات المواد ثنائية الأبعاد الحاسوبية (C2DB)، موردًا شاملاً لفهم الخصائص الإلكترونية والمغناطيسية للمواد ثنائية الأبعاد تحت تصحيحات Hubbard U، بينما تبرز أيضًا قيود طرق DFT+U في التقاط تعقيدات الأنظمة المرتبطة بشدة.
DOI: https://doi.org/10.1038/s41524-024-01503-3
PMID: https://pubmed.ncbi.nlm.nih.gov/39872025
Publication Date: 2025-01-24
Author(s): Sahar Pakdel et al.
Primary Topic: 2D Materials and Applications
Methods
In this study, the authors investigate the electronic and magnetic properties of two-dimensional (2D) materials containing 3d transition metals (V, Cr, Mn, Fe, Co, Ni, Cu) using a high-throughput computational approach. The materials were selected from the C2DB database, ensuring that each monolayer has an energy above the convex hull of less than 0.2 eV/atom, as determined by the Perdew-Burke-Ernzerhof (PBE) functional. This criterion guarantees that the materials are chemically stable and that the findings are not skewed by unrealistic structures. Additionally, properties were calculated only for monolayers where the PBE+U method resulted in a finite band gap, as the Hubbard U corrections are applicable solely to non-metallic systems.
The computational methodology employed the Atomic Simulation Recipes (ASR) Python framework and utilized the GPAW electronic structure code, with a plane wave cutoff of 800 eV and a k-point density of 12 Å. The DFT+U calculations incorporated an effective U-parameter, defined as \( U_{\text{eff}} = U – J \), to account for on-site Coulomb and exchange interactions. The total energy expression for the DFT+U method is given by \( E_{\text{DFT+U}} = E_{\text{DFT}} + \sum_{a,i} \frac{1}{2} U_{\text{eff},i} \text{Tr}(\rho_{a,i} – \rho_{a,i} \rho_{a,i}) \), where the sum is over atoms \( a \) and their corresponding 3d orbitals. The authors set \( U_{\text{eff}} \) to 4 eV for these orbitals and ensured that all structures were relaxed to a maximum force on the atoms of less than 0.01 eV/Å.
Results
In the results section, the authors first investigate the dependency of various properties on the Hubbard U parameter for a select group of materials. They determine a U-correction value of 4 eV for subsequent calculations, guided by the principle of Occam’s razor. However, they find that incorporating this U value leads to poorer alignment with experimental data regarding the in-plane lattice constants of the monolayers. As a result, all further property calculations are conducted using the PBE relaxed monolayers.
Additionally, the authors examine the impact of the Hubbard U on the band gap and magnetic properties across all monolayers studied. This comprehensive analysis aims to clarify the role of U in influencing key electronic and magnetic characteristics, thereby providing insights into the material properties under consideration.
Discussion
In this study, Pakdel, Olsen, and Thygesen systematically explore the impact of Hubbard U corrections on the electronic structure of 638 two-dimensional (2D) materials containing 3d transition metals, utilizing density functional theory (DFT) with both PBE and PBE+U approximations. The findings indicate that while the inclusion of the U correction generally increases band gaps—resulting in a metal-to-insulator transition for 134 materials (21%)—it adversely affects the accuracy of lattice constants, leading to a preference for PBE structures in subsequent evaluations. Notably, the magnetic moments exhibit minimal dependence on U, while both magnetic exchange coupling and magnetic anisotropy parameters are systematically reduced, suggesting that the Hubbard U correction may lower predicted Curie temperatures in these materials.
The analysis reveals that the U correction significantly influences the electronic properties, with band gaps typically increasing by 0-0.5 eV, although some materials display more complex behaviors. The study emphasizes that the magnetic properties remain largely unchanged with respect to U, as the magnetic moments remain integer values for insulating states. Furthermore, the research highlights the challenges of accurately describing 2D materials with localized electrons, as traditional DFT methods struggle with self-interaction and delocalization errors. The results, compiled in the Computational 2D Materials Database (C2DB), provide a comprehensive resource for understanding the electronic and magnetic characteristics of 2D materials under Hubbard U corrections, while also underscoring the limitations of DFT+U methods in capturing the complexities of strongly correlated systems.