Photogalvanic currents from first-principles real-time density-matrix dynamics

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

يولد التأثير الكهروضوئي الفوتوني (PGE)، أو تأثير الفولتوضوئية الحجمي، تيارًا مستمرًا (DC) في المواد غير المركزية التماثل تحت الإضاءة دون وجود مجالات خارجية. وبينما ربطت الأبحاث الحديثة التأثير الكهروضوئي الفوتوني بالخصائص الهندسية الكمومية (انحناء بيري، والمقياس الكمومي) في المواد الطوبولوجية، تواجه الدراسات الحالية القائمة على المبادئ الأولى (ab initio) قيودًا كبيرة:

• التركيز على الاستثارة الضوئية: ركزت النظريات السابقة بشكل كبير على عملية الاستثارة الضوئية الأولية فقط، متجاهلة أحداث التشتت اللاحقة.
• نماذج تشتت غير كافية: تتضمن المواد الواقعية عمليات تشتت معقدة تتوسطها البوزونات (الفوتونات والفونونات)، بما في ذلك الاسترخاء داخل النطاق وبين النطاقات وإعادة اتحاد الإلكترون-فجوة. وتفشل التقريبات الشائعة، مثل زمن الاسترخاء الوحيد المستقل عن الحالة ($\tau_0$)، في وصف هذه الآليات بشكل متسق ذاتيًا.
• الديناميكيات العابرة: هناك نقص في الأطر النظرية الصارمة لتفسير تيارات الفوتو-جلفانية العابرة (التي تُلاحظ عبر انبعاثات THz)، لا سيما أصل استجابات التيار "ثنائية القطب" (التي تعكس اتجاهها بمرور الوقت).
• مفهوم خاطئ حول تيار الإزاحة (Shift Current): يُفترض عمومًا أن تيار الإزاحة (تحت الضوء المستقطب خطيًا) هو تأثير ناتج عن الاستثارة الضوئية فقط، مع تجاهل المساهمة المحتملة للتشتت بوساطة الفونونات.

2. المنهجية: ديناميكيات مصفوفة الكثافة في الوقت الحقيقي من المبادئ الأولى (FPDMD)

طوّر المؤلفون إطارًا نظريًا شاملاً يعتمد على معادلات الماستر الكمومية (QME) لمحاكاة ديناميكيات الإلكترون-بوزون في الوقت الحقيقي.

• الهاملتوني (Hamiltonian): يتم نمذجة النظام باستخدام هاملتوني كلي $H_T = H_e^0 + H_b^0 + H'$، حيث $H_e^0$ هو الهاملتوني الإلكتروني غير المضطرب (المستمد من نظرية الكثافة الوظيفية لـ Kohn-Sham)، و $H_b^0$ هو هاملتوني البوزون (الفونون/الفوتون)، و $H'$ يمثل التفاعل بين الإلكترون والبوزون.
• حدود التفاعل:
  • الإلكترون-الفوتون: يتم التعامل معه في مقياس السرعة (velocity gauge) باستخدام تقريب ثنائي القطب.
  • الإلكترون-الفونون: يتم حسابه باستخدام عناصر المصفوفة من الدرجة الأولى المستمدة من نظرية الاضطراب الوظيفي للكثافة (DFPT).
• الديناميكيات: يتم التحكم في التطور الزمني لمصفوفة الكثافة المختزلة ($\rho$) بواسطة معادلة الماستر الكمومية (QME):
  $$\frac{d\rho}{dt} = -\frac{i}{\hbar}[H_e^0, \rho] + \left.\frac{d\rho}{dt}\right|_{exc} + \left.\frac{d\rho}{dt}\right|_{ph} + \left.\frac{d\rho}{dt}\right|_{rec}$$
  وهذا يتضمن صراحةً حدود الاستثارة الضوئية، وتشتت الإلكترون-فونون، وإعادة اتحاد الإلكترون-فجوة.
• حساب التيار: يتم حساب كثافة التيار الشحني $J(t)$ كأثر للمصفوفة لكثافة المصفوفة مضروبة في مؤثر السرعة. ويتم تفكيك التيار إلى مكونات قطرية للنطاق (تيار الحقن) ومكونات خارج قطر النطاق (تيار الإزاحة).
• الاتساق الذاتي: على عكس النهج الاضطرابي، تحل هذه الطريقة لتوزيعات الإلكترونات في حالة الاستقرار غير المتوازن وأزمنة الاسترخاء المعتمدة على الحالة ($\tau_{ks}$) بشكل متسق ذاتيًا.

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

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

4. النتائج الرئيسية (دراسة حالة: BaTiO$_3$)

طبق المؤلفون صياغة FPDMD على تيتانات الباريوم الرباعية (BaTiO$_3$)، وهي مادة كهرضغطية نموذجية.

• التأثير الكهروضوئي الخطي (تيار الإزاحة):

  • مساهمة الفونون: خلافًا للاعتقاد بأن تيار الإزاحة هو تأثير استثارة بحت، وجدت الدراسة أن التشتت بوساطة الفونونات يساهم بشكل كبير في تيار الإزاحة، خاصة عند ترددات الضوء الأعلى.
  • الآلية: ينشأ تيار إزاحة الفونون من تشتت الإلكترونات من سطح الاستثارة إلى حافة النطاق. وتكون magnitude هذه المساهمة مقاربة لتيار إزاحة الاستثارة، بل وتتفوق عليه عند الترددات العالية.
  • الهندسة الكمومية: يرتبط متجه إزاحة الفونون بـ انحناء بيري ($\Omega_k$). ويرتبط دوران انحناء بيري في BaTiO$_3$ بكهرباء الفيروكهربية (ferroelectricity) الخاصة به؛ حيث يؤدي عكس الدوران إلى عكس اتجاه التيار.
  • التوافق مع التجربة: إن تضمين مساهمة الفونون يجعل النظرية الإجمالية لتيار الفوتو-جلفاني أكثر توافقًا مع البيانات التجريبية عبر نطاق طاقة أوسع فوق حافة النطاق.

• التأثير الكهروضوئي الدائري (تيار الحقن):

  • تم تطوير نظرية متسقة ذاتيًا لتيار الحقن المستقر باستخدام أزمنة الاسترخاء المعتمدة على الحالة ($\tau_{ks}$) المستمدة من توزيع FPDMD.
  • أظهر أن موصلية تيار الحقن هي توليفة من خصائص هندسية كمومية: انحناء بيري (على سطح الاستثارة الضوئية) والمقياس الكمومي (المتعلق بتشتت الإلكترون-فونون داخل النطاق).

• التيار الكهروضوئي الفوتوني العابر:

  • ثنائي القطب مقابل أحادي القطب: تفسر المحاكاة التيارات العابرة ثنائية القطب ($J(t)$) الملاحظة تجريبيًا كعملية انتقال ديناميكي بين آليتين مختلفتين في التأخر الزمني:
    1. تيار إزاحة الاستثارة ($J_{exc}$): يستجيب بسرعة ($\tau_{eph} \approx 0.75$ فمتوثانية).
    2. تيار إزاحة الفونون ($J_{ph}$): يستجيب بتأخر يتوافق مع استرخاء الطاقة ($\tau_{er} \approx 25$ فمتوثانية).
  • العبور (Crossover): إذا كان لـ $J_{exc}$ و $J_{ph}$ إشارات متعاكسة، يكون التيار الإجمالي ثنائي القطب. وإذا كان لهما نفس الإشارة، يكون التيار أحادي القطب. وهذا يفسر السلوك المعتمد على التردد الملاحظ في مطيافية انبعاث THz.

5. الأهمية والأثر

• حل النزاعات: يعمل هذا العمل على حل الارتباك المتعلق بآليات التشتت التي تساهم في PGE، مثبتًا أن تشتت الفونونات ليس مجرد تصحيح طفيف، بل هو عامل مهيمن في المواد الواقعية مثل المواد الكهرضغطية.
• القدرة التنبؤية: نهج FPDMD تنبؤي للمواد الواقعية حيث لا يمكن تجاهل اقتران الإلكترون-فونون، مما يمهد الطريق لتصميم كواشف ضوئية وخلايا شمسية أفضل (خاصة للفوتونات منخفضة الطاقة).
• ما وراء تيارات الشحنة: الإطار قابل للتعميم على مراقبات أخرى، مثل التيارات الضوئية المغزلية (spin photocurrents) والتيارات الضوئية المدارية (orbital photocurrents)، والتي تعد بالغة الأهمية في المجالات الناشئة لتقنيات السبين-أوبترونيكس والأوربيت-رونيكس.
• البصريات غير الخطية: توفر النظرية طريقًا مباشرًا لحساب الاستجابات البصرية من الرتب العليا (مثل التوليد التوافقي الثاني) وتوصيف الخصائص الطوبولوجية في المواد الكمومية الجديدة.

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