Beyond the conventional Emery model: crucial role of long-range hopping for cuprate superconductivity

ملخص تقني: ما وراء نموذج إيمري التقليدي

بيان المشكلة
يُعد نموذج إيمري، وهو نموذج ثلاثي النطاقات يتكون من مدارات $d_{x^2-y^2}$ للنحاس ومدارات $p_x, p_y$ للأكسجين في المستوى، الإطار المعياري لوصف الموصلات الفائقة من نوع الكوبرات. ومع ذلك، فإن التطبيقات التقليدية لهذا النموذج عادة ما تقتصر على معاملات القفز (hopping parameters) بين أقرب جيران النحاس والأكسجين ($t_{pd}$)، وبين أقرب وأبعد جيران الأكسجين ($t_{pp}, t'_{pp}$). وقد واجهت الدراسات السابقة التي استخدمت هذه المجموعة "التقليدية" من المعاملات صعوبة في إعادة إنتاج مخطط الطور للموصلية الفائقة للكوبرات كمياً، وغالباً ما اعتمدت على مجموعات عنقودية صغيرة تعاني من آثار الحجم المحدود، أو على نماذج هوبارد أحادية النطاق قد تفتقر إلى الواقعية المادية اللازمة. وهناك سؤال مفتوح حاسم: هل معاملات القفز التقليدية كافية لوصف قبة الموصلية الفائقة من نوع $d$ وتماثل معامل الترتيب الصحيح، لا سيما عند مستويات تشويب الثقوب العالية؟

المنهجية
يستخدم المؤلفون تقريب الرأس الديناميكي (D$\Gamma$A) لدراسة نموذج إيمري. وهذا الامتداد المخططي لنظرية المجال المتوسط الديناميكي (DMFT) يسمح بمعالجة الارتباطات غير الموضعية على شبكات زخم دقيقة (تصل إلى $130 \times 130$)، مما يحاكي نظاماً لانهائياً بفعالية ويتجنب مشكلات الإشارة المرتبطة بمحاكاة مونت كارلو الكمية على العناقيد الكبيرة.

تستخدم الدراسة تجسيداً مادياً محدداً وهو $\text{CaCuO}_2$. تم اشتقاق الهاملتوني غير المتفاعل من حسابات نظرية الدالة الكثافية (DFT) باستخدام تقريب التدرج المعمم (GGA) والإسقاط على دالات وانييه الموضعية كحد أقصى. كما تم حساب قوة التفاعل $U = 6.5$ إلكترون فولت عبر تقريب ران-كابتر-رين (cRPA). ويقارن المؤلفون بين ثلاثة نماذج متميزة:

1. نموذج إيمري التقليدي: يتضمن فقط $t_{pd}$ و $t_{pp}$ و $t'_{pp}$.
2. نموذج إيمري التساهمي: مجموعة معاملات تقليدية مشتقة لمركب $\text{Bi}_2\text{Sr}_2\text{CaCu}_2\text{O}_8$ (Bi-2212) من الأدبيات السابقة.
3. نموذج d-p الكامل: يتضمن جميع عناصر القفز المشتقة من إسقاط DFT-Wannier، بما يتجاوز الجيران الأقرب والأبعد التاليين.

يتم تحديد عدم استقرار الموصلية الفائقة من خلال حل معادلة إيلياشبيرغ الخطية ضمن صياغة السلم لـ D$\Gamma$A، واستخراج القيمة الذاتية للموصلية الفائقة $\lambda_{SC}$. وتشير قيمة $\lambda_{SC} \to 1$ إلى درجة حرارة الانتقال $T_c$.

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

• عدم كفاية القفز التقليدي: يوضح المؤلفون أن نموذج إيمري التقليدي يقلل بشكل كبير من معاملات القفز الفعالة عند اختزالها إلى نموذج هوبارد أحادي النطاق فعال. وتحديداً، يكون القفز بين الجيران غير الأقرب $t'$ في النموذج المختزل صغيراً جداً (على سبيل المثال، $|t'/t| \approx 0.07$ لـ $\text{CaCuO}_2$) مقارنة بالقيم الواقعية ($|t'/t| \approx 0.15 - 0.25$). وينشأ هذا القصور لأن النموذج التقليدي يهمل مسارات القفز طويلة المدى الموجودة فعلياً في هاملتوني d-p الكامل.
• تصحيحات بنية النطاق: يؤدي تضمين معاملات القفز طويلة المدى في نموذج d-p الكامل إلى تصحيح إزاحة تفرد فان هوف (van Hove singularity) نحو طاقات أقل، مما يوافق نتائج بنية النطاق لـ DFT-Wannier. في المقابل، يضع النموذج التقليدي تفرد فان هوف قريباً جداً من طاقة فيرمي.
• مخطط طور الموصلية الفائقة:
  • النموذج التقليدي: يفشل في إظهار عدم استقرار الموصلية الفائقة عند تشويب الثقوب $\delta \geq 0.15$ (نطاق التشويب الأمثل). وهو يعطي قبة موصلية فائقة ضيقة جداً (9%–14% تشويب) ويتنبأ بمعامل ترتيب من نوع $d$ معدل بقوة مع وزن طيفي منخفض عند نقاط الـ antinodes. ويُعزى هذا التعديل إلى تطور فجوة شبه (pseudogap) قبل أوانها بسبب وضع تفرد فان هوف في مكان خاطئ.
  • نموذج d-p الكامل: ينجح في إعادة إنتاج قبة موصلية فائقة عريضة تتوافق مع ظواهر الكوبرات، حيث تمتد تقريباً من 7% إلى 22% تشويب ثقوب. وعند $\delta = 0.15$، يُلاحظ بوضوح عدم استقرار الموصلية الفائقة. ويؤدي معامل الترتيب من نوع $d$ الناتج إلى شكل تقليدي، دون التعديل الزائف المشاهد في النماذج المقيدة.
• دور القفز طويل المدى: تثبت الدراسة أن معاملات القفز طويلة المدى (التي تتجاوز المعاملات الثلاثة التقليدية) ضرورية كمياً للحصول على مخطط الطور الصحيح ومعامل الترتيب. وهذه المعاملات تقع في نفس نطاق حد $t'$ في نموذج هوبارد، وهو أمر معروف بأهميته الحاسمة للموصلية الفائقة.

الأهمية
يزعم البحث أنه بينما يعد نموذج إيمري التقليدي مفيداً نوعياً، فإنه غير كافٍ كمياً لوصف الموصلية الفائقة في الكوبرات عند مستويات التشويب العالية. ويجادل المؤلفون بأن قصر الحسابات على معاملات القفز الثلاثة التقليدية يؤدي إلى تنبؤات خاطئة فيما يتعلق بمدى قبة الموصلية الفائقة وطبيعة معامل الترتيب. ومن خلال استخدام نموذج d-p الكامل المشتق من المبادئ الأولى وطريقة D$\Gamma$A، يقدم العمل وصفاً أكثر واقعية للمواد يتضمن طبيعياً مدارات الأكسجين وتأثيرات القفز طويلة المدى. يحل هذا النهج التناقضات المتعلقة بموقع تفرد فان هوف وسلوك الفجوة شبه الناتجة، مما يمهد الطريق نحو فهم نظري أكثر دقة للموصلية الفائقة عالية الحرارة دون اللجوء إلى ملاءمة المعاملات التجريبية.