DOI: https://doi.org/10.1038/s41524-025-01855-4
تاريخ النشر: 2026-01-10
المؤلف: Max Großmann وآخرون
الموضوع الرئيسي: الهياكل البلورية والخصائص
نظرة عامة
في هذه الدراسة، يقوم المؤلفون بتقييم مختلف أساليب نظرية الاضطراب ذات الجسم المتعدد (MBPT) مقابل نظرية الوظائف الكثافة (DFT) للتنبؤ بفجوات النطاق للمواد الصلبة. يقومون بتقييم أربعة متغيرات GW بشكل منهجي: G0W0 باستخدام تقريب قطب البلازما Godby-Needs (G0W0-PPA)، وquasiparticle الكامل G0W0 (QPG0W0)، وquasiparticle الذاتي المتسق الكامل GW (QSGW)، وQSGW مع تصحيحات الرأس (QSGŴ). يتم مقارنة أداء هذه الطرق مع الوظائف المعروفة جيدًا DFT مثل mBJ وHSE06.
تشير النتائج إلى أن G0W0-PPA يوفر تحسينات هامشية فقط في الدقة مقارنة بأفضل طرق DFT، على الرغم من تكلفته الحسابية الأعلى. ومع ذلك، فإن استبدال PPA بالتكامل الكامل للتشويش الكهربائي يعزز بشكل كبير دقة التنبؤ، مما يتماشى تقريبًا مع النتائج من QSGŴ. بينما يقوم QSGW بتخفيف انحياز نقطة البداية بشكل فعال، فإنه يميل إلى المبالغة في تقدير فجوات النطاق التجريبية بحوالي 15%. إن إدخال تصحيحات الرأس في التفاعل الكولومبي المشوش من خلال QSGŴ يصحح هذه المبالغة، مما ينتج عنه فجوات نطاق ليست دقيقة فحسب، بل قادرة أيضًا على تحديد البيانات التجريبية المشكوك فيها.
مقدمة
تؤكد مقدمة هذه الورقة البحثية على أهمية التنبؤ بدقة بفجوة النطاق للمواد شبه الموصلة والعازلة للتطبيقات البصرية والإلكترونية الضوئية. بينما تُستخدم نظرية الوظائف الكثافة (DFT) على نطاق واسع لمثل هذه التنبؤات، فإنها غالبًا ما تقلل من تقدير فجوات النطاق بسبب تفسير فجوة كوهين-شام. أظهرت دراسة مرجعية أجراها بورليدو وآخرون أن الوظائف المتقدمة لـ DFT، وخاصة تقريب التدرج العام الميتا (meta-GGA) والوظائف الهجينة، يمكن أن تخفف من هذه المشكلة، على الرغم من أن هذه التحسينات غالبًا ما تستند إلى تعديلات تجريبية بدلاً من أسس نظرية صلبة. بالمقابل، تقدم نظرية الاضطراب ذات الجسم المتعدد (MBPT)، وبشكل خاص تقريب GW، نهجًا أكثر صرامة لحساب الهياكل الإلكترونية وفجوات النطاق، على الرغم من تكلفته الحسابية الأعلى.
يقترح المؤلفون إجراء تقييم منهجي يقارن MBPT مع أفضل الوظائف أداءً لـ DFT، مع التركيز بشكل خاص على تقريب GW. يستخدمون مجموعة بيانات واسعة من فجوات النطاق التجريبية لـ 472 مادة غير مغناطيسية، تم تجميعها في الأصل بواسطة بورليدو وآخرين، لتسهيل المقارنة المباشرة. تهدف الدراسة إلى معالجة الحاجة الملحة لمجموعات بيانات موثوقة في سياق تطبيقات التعلم الآلي في علوم المواد، حيث تعتمد التنبؤات الدقيقة على جودة بيانات التدريب. تسلط الأبحاث الضوء على ضرورة إنشاء معيار لـ MBPT مقابل DFT لإبلاغ اختيار الطرق الحسابية لتوليد مجموعات بيانات عالية الدقة، خاصة للتعلم الانتقالي في نماذج التعلم الآلي.
طرق
في هذا القسم، يقوم المؤلفون بتقييم مختلف الطرق الحسابية للتنبؤ بفجوات النطاق مقابل البيانات التجريبية، مستخدمين عدة مقاييس خطأ إحصائية: متوسط الخطأ (ME)، والانحراف المعياري للخطأ (σ)، ومتوسط الخطأ المطلق (MAE)، وجذر متوسط مربع الخطأ (RMSE)، ومتوسط الخطأ النسبي المطلق (MAPE). يتم تعريف الخطأ على أنه الفرق بين فجوات النطاق المحسوبة والتجريبية. يتم أيضًا فحص تأثير اقتران الدوران (SOC) على النتائج، مع تطبيق SOC كتصحيح لهياكل النطاق الكوانتي. تشير النتائج إلى أن طريقة QSGW تتفوق باستمرار على الآخرين، مثل QPG0W0 وmBJ وHSE06، التي تميل إلى تقليل تقدير فجوات النطاق. من الجدير بالذكر أن تضمين SOC يحسن دقة نتائج QPG0W0 وQSGW، خاصة للمواد ذات العناصر الثقيلة.
يكشف التحليل أنه بينما يبدو أن MAPE مواتية لطرق مثل mBJ وHSE06، إلا أنها مضللة بسبب حساسيتها للمواد ذات فجوات النطاق الصغيرة. يجادل المؤلفون بضرورة إعطاء الأولوية لـ MAE وRMSE كمقاييس أكثر موثوقية لتقييم الطرق الحسابية. كما يناقشون مزايا استخدام طريقة QPG0W0، التي، على الرغم من عدم وجود اتساق ذاتي، توفر نتائج قابلة للمقارنة مع طريقة QSGW الأكثر كثافة حسابيًا. يختتم القسم بالتأكيد على أهمية تضمين تصحيحات SOC للتنبؤات الدقيقة في المواد ذات العناصر الثقيلة، كما يتضح من التحسينات الكبيرة في حسابات فجوة النطاق للمواد مثل PbTe وSnTe.
نتائج
يقيم قسم النتائج أداء أربع طرق GW: G$_0$W$_0$-PPA وQPG$_0$W$_0$ وQSGW وQSGŴ، مع التركيز على دقتها في التنبؤ بفجوات النطاق. تحقق طريقة G$_0$W$_0$-PPA التقارب ضمن 25 ميلي فولت، بينما تتقارب الطرق الأكثر تعقيدًا (ii-iv) ضمن 100 ميلي فولت. يتم تقديم تحليل شامل للخطأ، يوضح مقاييس مثل متوسط الخطأ (ME) والانحراف المعياري (σ) ومتوسط الخطأ المطلق (MAE) وجذر متوسط مربع الخطأ (RMSE) ومتوسط الخطأ النسبي المطلق (MAPE) عبر مجموعات بيانات مختلفة.
تشير النتائج إلى أن طريقة QSGŴ، خاصة عند تضمين اقتران الدوران (SOC)، تظهر أقل MAE (0.30 eV) وRMSE (0.47 eV) بين الطرق التي تم تقييمها، مما يشير إلى دقة متفوقة. من الجدير بالذكر أن الفروق في تصنيفات فجوة النطاق بين تطبيقات DFT المختلفة تمت مناقشتها، والتي تُعزى بشكل أساسي إلى معالجة الحالات شبه الأساسية في الحسابات. يسلط هذا التحليل الضوء على أهمية الخيارات المنهجية في علوم المواد الحسابية وتأثيرها على تفسير الخصائص الإلكترونية.
مناقشة
في هذا القسم، يناقش المؤلفون تحديد وتحليل القيم الشاذة في حسابات فجوة النطاق باستخدام طرق QPG 0 W 0 وQSGW، التي أظهرت دقة عالية. يركزون على المواد التي تظهر انحرافات كبيرة عن فجوات النطاق التجريبية المبلغ عنها، وخاصة تلك التي يقلل فيها QSGW من تقدير فجوة النطاق، على عكس ميله النموذجي للمبالغة. تشمل القيم الشاذة الملحوظة AgF وCuLaO2، حيث أن فجوة النطاق لـ QPG 0 W 0 للأولى البالغة 0.91 eV أقل بكثير من القيمة التجريبية البالغة 2.8 eV، والتي تُعزى إلى حالتها الأرضية المعدنية المتوقعة بواسطة LDA. بالمثل، فإن قيمة QPG 0 W 0 لـ CuLaO2 البالغة 4.35 eV تفرط في تقدير فجوة النطاق التجريبية البالغة 2.34 eV، مع اقتراح المؤلفين أن الفروق قد تنشأ من تفاعلات الإلكترون-فونون وطبيعة فجوات النطاق البصرية مقابل الأساسية.
يستفيض المؤلفون في توضيح قيم شاذة أخرى، مثل Cu2S وCu2Se، التي وُجد أنها شبه معدنية في حسابات QSGW، مما يعقد تحديد فجوات النطاق التجريبية بسبب الاضطراب الهيكلي والارتباطات الإلكترونية القوية. كما يناقشون RbAu، حيث تنحرف فجوة النطاق المحسوبة لـ QSGW البالغة 1.61 eV عن القيمة التجريبية البالغة 2.75 eV، مما يشير إلى الحاجة لمزيد من التحقيق التجريبي. يختتم القسم بملخص لقيود طرق QSGW وQSGW، مؤكدًا أن الفروق غالبًا ما تنشأ من نقص تصحيحات الرأس ومعالجة تفاعلات الإلكترون-فونون، التي تعتبر حاسمة لتحقيق توقعات دقيقة لفجوة النطاق. يدعو المؤلفون إلى نهج أكثر شمولاً لدراسة هذه المواد، مع إمكانية دمج طرق الجسم المتعدد المتقدمة ونماذج هيكلية واقعية لمواءمة التوقعات النظرية مع النتائج التجريبية.
DOI: https://doi.org/10.1038/s41524-025-01855-4
Publication Date: 2026-01-10
Author(s): Max Großmann et al.
Primary Topic: Crystal Structures and Properties
Overview
In this study, the authors benchmark various many-body perturbation theory (MBPT) approaches against density functional theory (DFT) for predicting the band gaps of solids. They systematically evaluate four GW variants: G0W0 using the Godby-Needs plasmon-pole approximation (G0W0-PPA), full-frequency quasiparticle G0W0 (QPG0W0), full-frequency quasiparticle self-consistent GW (QSGW), and QSGW with vertex corrections (QSGŴ). The performance of these methods is compared to the well-regarded DFT functionals mBJ and HSE06.
The findings indicate that G0W0-PPA provides only marginal improvements in accuracy over the best DFT methods, albeit at a higher computational cost. However, substituting the PPA with full-frequency integration of dielectric screening significantly enhances prediction accuracy, nearly aligning with the results from QSGŴ. While QSGW effectively mitigates starting-point bias, it tends to overestimate experimental band gaps by approximately 15%. The introduction of vertex corrections in the screened Coulomb interaction through QSGŴ rectifies this overestimation, yielding band gaps that are not only accurate but also capable of identifying questionable experimental data.
Introduction
The introduction of this research paper emphasizes the significance of accurately predicting the band gap of semiconductors and insulators for optical and optoelectronic applications. While density functional theory (DFT) is widely used for such predictions, it often underestimates band gaps due to the interpretation of the Kohn-Sham gap. A benchmark study by Borlido et al. demonstrated that advanced DFT functionals, particularly meta-generalized gradient approximation (meta-GGA) and hybrid functionals, can mitigate this issue, although these improvements are often based on empirical adjustments rather than solid theoretical foundations. In contrast, many-body perturbation theory (MBPT), specifically the GW approximation, offers a more rigorous approach to calculating electronic structures and band gaps, albeit at a higher computational cost.
The authors propose a systematic benchmark comparing MBPT with the best-performing DFT functionals, specifically focusing on the GW approximation. They utilize an extensive dataset of experimental band gaps for 472 non-magnetic materials, originally curated by Borlido et al., to facilitate a direct comparison. The study aims to address the pressing need for reliable datasets in the context of machine learning applications in materials science, where accurate predictions depend on the quality of training data. The research highlights the necessity of establishing a benchmark for MBPT against DFT to inform the choice of computational methods for generating high-fidelity datasets, particularly for transfer learning in machine learning models.
Methods
In this section, the authors benchmark various computational methods for predicting band gaps against experimental data, employing several statistical error metrics: mean error (ME), standard deviation of the error (σ), mean absolute error (MAE), root mean squared error (RMSE), and mean absolute percentage error (MAPE). The error is defined as the difference between calculated and experimental band gaps. The impact of spin-orbit coupling (SOC) on the results is also examined, with SOC applied as a correction to the quasiparticle band structures. The findings indicate that the QSGW method consistently outperforms others, such as QPG0W0, mBJ, and HSE06, which tend to underestimate band gaps. Notably, the inclusion of SOC improves the accuracy of QPG0W0 and QSGW results, particularly for materials with heavy elements.
The analysis reveals that while MAPE appears favorable for methods like mBJ and HSE06, it is misleading due to its sensitivity to materials with small band gaps. The authors argue for prioritizing MAE and RMSE as more reliable metrics for evaluating computational methods. They also discuss the advantages of using the QPG0W0 method, which, despite lacking self-consistency, provides results comparable to the more computationally intensive QSGW method. The section concludes by emphasizing the importance of including SOC corrections for accurate predictions in materials with heavy elements, as demonstrated by significant improvements in band gap calculations for compounds like PbTe and SnTe.
Results
The results section evaluates the performance of four GW methods: G$_0$W$_0$-PPA, QPG$_0$W$_0$, QSGW, and QSGŴ, with a focus on their accuracy in predicting band gaps. The G$_0$W$_0$-PPA method achieves convergence within 25 meV, while the more complex methods (ii-iv) converge within 100 meV. A comprehensive error analysis is presented, detailing metrics such as mean error (ME), standard deviation (σ), mean absolute error (MAE), root mean square error (RMSE), and mean absolute percentage error (MAPE) across various datasets.
The findings indicate that the QSGŴ method, particularly when including spin-orbit coupling (SOC), exhibits the lowest MAE (0.30 eV) and RMSE (0.47 eV) among the methods evaluated, suggesting superior accuracy. Notably, discrepancies in band gap classifications between different DFT implementations are discussed, attributed primarily to the treatment of semicore states in the calculations. This analysis highlights the importance of methodological choices in computational materials science and their impact on the interpretation of electronic properties.
Discussion
In this section, the authors discuss the identification and analysis of outliers in band gap calculations using the QPG 0 W 0 and QSGW methods, which demonstrated high accuracy. They focus on materials that exhibit significant deviations from reported experimental band gaps, particularly those where QSGW underestimates the band gap, contrary to its typical tendency to overestimate. Notable outliers include AgF and CuLaO2, where the former’s QPG 0 W 0 band gap of 0.91 eV is significantly lower than the experimental value of 2.8 eV, attributed to its metallic ground state predicted by LDA. Similarly, CuLaO2’s QPG 0 W 0 value of 4.35 eV overestimates the experimental band gap of 2.34 eV, with the authors suggesting that the discrepancies may arise from electron-phonon interactions and the nature of optical versus fundamental band gaps.
The authors further elaborate on other outliers, such as Cu2S and Cu2Se, which are found to be semimetallic in QSGW calculations, complicating the determination of their experimental band gaps due to structural disorder and strong electronic correlations. They also discuss RbAu, where the calculated QSGW band gap of 1.61 eV deviates from the experimental value of 2.75 eV, suggesting a need for further experimental investigation. The section concludes with a summary of the limitations of the QSGW and QSGW methods, emphasizing that discrepancies often arise from missing vertex corrections and the treatment of electron-phonon interactions, which are critical for achieving accurate band gap predictions. The authors advocate for a more comprehensive approach to studying these materials, potentially incorporating advanced many-body methods and realistic structural models to align theoretical predictions with experimental results.