Extraction of the self energy and Eliashberg function from angle resolved photoemission spectroscopy using the xARPES code

1. بيان المشكلة

تعد مطيافية الانبعاث الضوئي بالزاوية المحللة (ARPES) أداة رئيسية لدراسة التفاعلات بين الإلكترون والبوزون، لا سيالما اقتران الإلكترون-فونون (EPC)، من خلال تحليل الدالة الطيفية للإلكترون. ومع ذلك، فإن استخراج المعلمات المتعددة للأجسام (تحديداً الطاقة الذاتية $\Sigma(E)$ ودالة إلياشبرج $\alpha^2F(\omega)$) من البيانات التجريبية يواجه عدة تحديات:

• التشتتات غير الخطية: غالباً ما تعتمد الطرق التقليدية تقريبات لتشتت النطاقات غير المتفاعلة ($\varepsilon_n(k)$) كدوال خطية. يؤدي هذا إلى انحياز في استخراج الطاقة الذاتية عندما تكون النطاقات منحنية (قطع مكافئ أو رتب أعلى)، حيث لا تكون الدالة الطيفية مجرد دالة لورنتزية بسيطة في فضاء الزخم.
• التعيين اليدوي للمعلمات: يعتمد تفكيك إجمالي الطاقة الذاتية إلى مساهمات الفونون ($\Sigma_{ph}$)، والإلكترون-إلكترون ($\Sigma_{el}$)، والشوائب ($\Sigma_{imp}$) عادةً على الفحص البصري اليدوي أو الافتراضات التعسفية، مما يؤدي إلى إدخال تحيز بشري.
• تأثيرات عناصر المصفوفة: إن تجاهل عناصر مصفوفة الانبعاث الضوئي (PMEs) يمكن أن يشوه الطاقة الذاتية المستخرجة بشكل كبير، وأحياناً بعامل يصل إلى ضعفين.
• عدم استقرار المسألة العكسية: يعد استخراج $\alpha^2F(\omega)$ من $\Sigma(E)$ مسألة عكسية سيئة التحديد. تتطلب طرق "الحد الأقصى للإنتروبيا" (MEM) القياسية ضبطاً دقيقاً للمعلمات الفائقة والمعرفة المسبقة، وغالباً ما تفتقر إلى إطار إحصائي صارم لتحسين المعلمات.

2. المنهجية

يقدم المؤلفون xARPES، وهو كود بلغة بايثون (مرخص بموجب رخصة GPL-v3) ينفذ سير عمل صارماً ومؤتمتاً لاستخراج الطاقات الذاتية ودوال إلياشبرج. تتكون المنهجية من ثلاث مراحل رئيسية:

أ. ملاءمة MDC متقدمة مع تشتتات منحنية

• بدلاً من افتراض نطاقات خطية، يقوم xARPES بملاءمة منحنيات توزيع الزخم (MDCs) باستخدام تشتتات غير متفاعلة متعددة الحدود (مثل القطع المكافئ).
• يقوم الكود برسم خرائط زوايا الكاشف إلى زخوم البلورة، مع مراعاة الهندسة التجريبية ودوران العينة.
• يتضمن تصحيحات عنصر مصفوفة الانبعاث الضوئي (PME)، مما يسمح للمستخدم بتحديد عناصر المصفوفة المعتمدة على الزاوية $|M_n(\eta)|^2$ لتصحيح التباينات في الوزن الطيفي الناتجة عن الخصائص المدارية واستقطاب الضوء.
• تستخرج عملية الملاءمة هذه مواقع القمم غير البعدية ($\tilde{r}_n$) وعروضها ($\tilde{r}_n$)، والتي يتم تحويلها بعد ذلك إلى الأجزاء الحقيقية والتخيلية للطاقة الذاتية، $\tilde{\Sigma}'_n(E)$ و $-\tilde{\Sigma}''_n(E)$.

ب. تفكيك الطاقة الذاتية

يتم تفكيك إجمالي الطاقة الذاتية وفقاً لقاعدة ماتيسين:
$$\Sigma_n(E) = \Sigma_{ph,n}(E) + \Sigma_{el,n}(E) + \Sigma_{imp,n}(E)$$

• $\Sigma_{imp}$: تُنمذج كمعامل اضمحلال تخيلي ثابت ($-i\Gamma_{imp}$).
• $\Sigma_{el}$: تُنمذج باستخدام دالة ظاهراتية متسقة مع علاقة كرامرز-كرونيغ تحقق سلوك سائل فيرمي عند الطاقات المنخفضة وتضمحل كـ $|E|^{-2}$ عند الطاقات العالية لتجنب التباعد فوق البنفسجي.
• $\Sigma_{ph}$: تهيمن عليها حد فان-ميغدال الديناميكي، المرتبط بدالة إلياشبرج عبر معادلة تكامل النواة.

ج. الاستدلال البايزي وطريقة الحد الأقصى للإنتروبيا (MEM)

• استخراج بضربة واحدة: يستخدم الكود طريقة MEM لعكس معادلة التكامل التي تربط $\Sigma_{ph}(E)$ بـ $\alpha^2F(\omega)$. يتضمن هذا معرفة مسبقة (مثل شبه الإيجابية لـ $\alpha^2F$) عبر حد إنتروبيا "شانون-جاينز" المعمم.
• التحسين ثنائي المستوى: المساهمة الجديدة هي دمج MEM ضمن حلقة استدلال بايزي.
  • الحلقة الداخلية تعمل على تحسين $\alpha^2F(\omega)$ لمجموعة ثابتة من معلمات النموذج.
  • الحلقة الخارجية تعمل على تحسين معلمات النموذج (كتلة غير التفاعل $m_b$، متجه موجة فيرمي $k_F$، قوة الشوائب $\Gamma_{imp}$، اقتران الإلكترون-إلكترون $\lambda_{el}$، وارتفاع دالة النموذج $h$) لتعظيم الاحتمالية البعدية.
  • هذا يلغي الضبط اليدوي للمعلمات، مما يوفر تحديداً قوياً إحصائياً للمعلمات الأكثر احتمالاً.

3. المساهمات الرئيسية

1. برمجية xARPES: إصدار حزمة بايثون شاملة مفتوحة المصدر لأتمتة استخراج معلمات الأجسام المتعددة من بيانات ARPES، مع دعم كل من التشتتات الخطية والمنحنية (القطع المكافئ).
2. التعامل مع التشتت المنحني: تقديم صياغة تتعامل بشكل صحيح مع تشتتات النطاقات غير الخطية، مما يزيل الانحياز المتأصل في ملاءمة لورنتز التقليدية للنطاقات المنحنية.
3. التحسين البايزي المؤتمت: إطار عمل جديد يحسن في آن واحد معلمات التشتت غير المتفاعل ومقدار قنوات التشتت (الإلكترون-إلكترون، الشوائب) جنباً إلى جنب مع دالة إلياشبرج، مما يزيل التحيز البشري.
4. تصحيح عنصر المصفوفة: التنفيذ الصريح لتصحيحات PME، مما يوضح أهميتها الحاسمة للتحليل الكمي للطاقة الذاتية.
5. مقياس المقارنة الكمية: تقديم مقياس المسافة الإقليدية للمقارنة الكمية بين دوال إلياشبرج المستخرجة من فروع أو مجموعات بيانات مختلفة.

4. النتائج والتحقق

أ. التحقق من البيانات الاصطناعية

• باستخدام بيانات اصطناعية ذات معلمات حقيقية معروفة، أظهر xARPES قدرة على استعادة الطاقة الذاتية الحقيقية ودالة إلياشبرج بدقة عالية.
• أظهرت الطريقة أن ملاءمة التشتت المنحني تعطي نتائج غير منحازة (تداخل بنسبة 95% في فاصل الثقة) حتى مع وجود دقة طاقة محدودة، بينما تؤدي ملاءمة لورنتز/الخطية إلى انحياز كبير عند طاقات الارتباط العالية.
• نجحت الحلقة البايزية في استعادة معلمات النموذج (الكتلة، $k_F$، قوى الاقتران) ضمن حوالي 5% من القيم الحقيقية لمستويات الضوضاء الواقعية.

ب. دراسة حالة 1: SrTiO$_3$ المنتهي بـ TiO$_2$

• النظام: سائل إلكتروني ثنائي الأبعاد (2DEL) على Nb-doped SrTiO$_3$.
• النتائج:
  • أظهرت دوال إلياشبرج المستخرجة من الفروع الداخلية اليسرى والداخلية اليمنى لنطاق $d_{xy}$ تشابهاً عالياً ($M \approx 0.01$)، مما يشير إلى تناظر كامن.
  • تأثير PME: أدى تجاهل تصحيحات عناصر المصفوفة إلى جعل الطاقات الذاتية أكبر بنحو ضعف القيم المصححة، مما يبرز ضرورة PME في التحليل الكمي.
  • تم تحديد أنماط الفونون (TO4, LO4) بما يتوافق مع الدراسات السابقة ولكن بدقة محسنة وآثار جانبية أقل.

ج. دراسة حالة 2: الجرافين المطعّم بالليثيوم

• النظام: جرافين شبه حر شبه مستقر مطعم بالليثيوم على Co(0001).
• النتائج:
  • تم تطبيق الطريقة على تشتتات مخروط ديراك الخطية.
  • أظهرت دوال إلياشبرج المستخرجة من اثنين من الكينات (kinks) المتناظرة بنيوياً (الفروع L و R) توافقاً غير مسبوق ($M \approx 0.0067$)، وهو أفضل بكثير من المقارنات التجريبية السابقة.
  • أكدت النتائج أن عدم التماثل التجريبي (مثل إزاحات حافة فيرمي) كان المصدر الرئيسي للاختلافات الطفيفة، مما يثبت قوة طريقة الاستخراج.
  • تم تحديد أنماط الفونون المهيمنة حول 166 ميلي إلكترون فولت (أنماط A'$_1$ و E$_{2g}$).

5. الأهمية

• الربط بين التجربة والنظرية: يوفر xARPES طريقة معيارية وغير منحازة لاستخراج دوال إلياشبرج، مما يتيح مقارنة مباشرة وعادلة مع حسابات المبادئ الأولى (مثل DFPT, GW, DMFT).
• قابلية التكرار: من خلال أتمتة اختيار المعلمات عبر الاستدلال البايزي، يزيل الكود الطبيعة "الصندوق الأسود" للملاءمة اليدوية، مما يضمن أن النتائج قابلة للتكرار وموضوعية.
• معيار للمجتمع العلمي: يضع إصدار xARPES معياراً جديداً لتحليل بيانات ARPES، خاصة للأنظمة ذات هياكل النطاقات المعقدة أو التفاعلات المتعددة للأجسام القوية.
• الاتجاهات المستقبلية: تسمح بنية الكود النمطية بتوسعات مستقبلية، مثل مصطلحات التشتت ذات الرتب الأعلى، وتضمين اقتران البلازمون، والمقارنة المباشرة مع تشتتات غير التفاعل النظرية.