Full Minimal Coupling All-Electron Real-Time TDDFT for X-Ray-Matter Interactions

ملخص تقني: الاقتران الأدنى الكامل (Full Minimal Coupling) لنموذج "الكل-إلكترون" في نظرية الكثافة الوظيفية المعتمدة على الزمن في الوقت الحقيقي (RT-TDDFT) لتفاعلات الأشعة السينية مع المادة

بيان المشكلة
تعتمد نظرية الكثافة الوظيفية المعتمدة على الزمن في الوقت الحقيقي (RT-TDDFT) القياسية عادةً على اقتران المجالات الكهرومغناطيسية بالمادة عبر جهد متجه (dipole) منتظم مكانياً، $\mathbf{A}(t)$. يفترض هذا التقريب أن الطول الموجي للفوتون، $\lambda$، أكبر بكثير من الامتداد المكاني للنظام. وبينما يكون هذا التقريب صالحاً في أطياف الضوء المرئي وفوق البنفسجي (UV)، إلا أنه ينهار في أطياف الأشعة السينية والأشعة فوق البنفسجية المتطرفة (XUV) حيث يقترب $\lambda$ من المقاييس بين الذرية (على سبيل المثال، $\lambda \approx 2.5$ Å عند 5 keV). وفي هذه الحالة، لا يمكن إهمال الاعتماد المكاني للجهد المتجه $\mathbf{A}(\mathbf{r}, t)$. إن إهمال هذا الاعتماد يؤدي إلى فقدان زخم الفوتون المحدود ($k_\gamma$)، مما يجعل الطريقة غير قادرة على وصف عمليات انتقال الزخم، وعدم التماثل غير ثنائي القطب (nondipole asymmetries)، والظواهر غير الخطية مثل توليد التوافق الثاني للأشعة السينية (SHG) أو خلط الموجات التي تحدث عند متجهات موجات غير صفرية في الفضاء المتبادل. علاوة على ذلك، تعتمد التنفيذات القياسية غالباً على الجهود الكاذبة (pseudopotentials) أو الطرق القائمة على الشبكات (grid-based) التي تواجه صعوبة في التعامل مع التكلفة الحسابية لمعالجة حالات اللب (core states) شديدة الترابط أو عدم استقرار القواعد المرتبط بالجهود الكاذبة شبه-اللب المطلوبة لمطيافية الأشعة السينية.

المنهجية
يقدم المؤلفون تنفيذاً للاقتران الأدنى الكامل ضمن كود FHI-aims، باستخدام مجموعة قواعد مدارات ذرية ممركزة عددية (NAO). هذا النهج هو نهج "الكل-إلكترون" (all-electron)، مما يتجنب تعقيدات الجهود الكاذبة لحالات اللب. تتضمن المنهجية أربعة مكونات تقنية رئيسية:

1. هاملتوني الاقتران الأدنى الكامل: يتضمن انتشار معادلة كوهن-شارك المعتمدة على الزمن الجهد المتجه المتغير مكانياً $\mathbf{A}(\mathbf{r}, t)$ في كل من حدي التفاعل البارامغناطيسي ($\mathbf{A} \cdot \hat{\mathbf{p}}$) والديامغناطيسي ($\mathbf{A}^2$). يتم التعامل مع التعديل المكاني عبر عناصر مصفوفة التدرج المعدلة ($M^{\cos}$ و $M^{\sin}$) في قاعدة المدارات الذرية، مع الاحتفاظ بعوامل الطور $e^{\pm i \mathbf{k}_\gamma \cdot \mathbf{r}}$.
2. اقتران النطاق الجانبي لزخم الفوتون: لاستعادة استجابة الكثافة عند متجه موجة الفوتون $\mathbf{k}_\gamma$ (التي تُفقد في تقريب $k$-diagonal القياسي)، قدم المؤلفون توسعاً اضطرابياً لمصفوفة الكثافة. يقومون بنشر الكتل خارج القطر لـ $k$، وتحديداً النطاق الجانبي $\ell = +1$ المتمثل في $\rho(\mathbf{k} + \mathbf{k}_\gamma, \mathbf{k}; t)$، مدفوعاً بالكتلة $\ell=0$. يسمح هذا بحساب عوامل شكل الكثافة عند متجهات الشبكة المتبادلة غير المتوافقة.
3. تشخيص التيار في الفضاء المتبادل: يقوم التنفيذ بحساب كثافة التيار الفيزيائي $\mathbf{J}(\mathbf{r}, t)$ المفككة إلى مكونات بارامغناطيسية وديامغناطيسية. ومن الأهمية بمكان أنه يحلل هذه التيارات في الفضاء المتبادل عند متجهات $\mathbf{G}_n = n\mathbf{k}_\gamma$ محددة (تحديداً $\mathbf{G}_1 = \mathbf{k}_\gamma$ و $\mathbf{G}_2 = 2\mathbf{k}_\gamma$). وهذا يتيح مراقبة هياكل متعددة الأقطاب والاستجابات التوافقية غير الخطية التي تكون غير مرئية للتيار العياني ($G=0$).
4. قاعدة المدارات العددية الممركزة (NAO) للكل-إلكترون: من خلال استخدام NAOs، تصف الطريقة حالات اللب شديدة الترابط بكفاءة دون ظهور الحالات الشبحية (ghost states) أو عدم استقرار القواعد المرتبط بالجهود الكاذبة شبه-اللب، مما يسمح بمعالجة الجزيئات المعزولة، والصفائح، والمواد الصلبة (bulk) على قدم المساوايد.

النتائج الرئيسية والاختبارات المرجعية
تم التحقق من صحة التنفيذ عبر أربعة اختبارات مرجعية تمتد من 100 إلكترون فولت إلى 5 كيلو إلكترون فولت:

• البنزين عند حافة كربون K (263–303 إلكترون فولت): تم إجراء تحقق مباشر عبر الكود مقابل مرجع Maxwell–TDDFT ذاتي الاتساق (Bonafé et al.). أظهر ثنائي القطب المستحث الناتج عن الاقتران الأدنى الكامل (FMC) تصحيحاً يتجاوز ثنائي القطب بنسبة 1–3.5% بالنسبة لتقريب ثنائي القطب (ED)، وهو ما يتوافق مع التقديرات الهندسية لـ $(k_\gamma r)^2/2$. أكد هذا التنفيذ الصحيح للاقتران $\mathbf{A} \cdot \hat{\mathbf{p}}$ المتغير مكانياً.
• الثيوفين عبر حواف الكبريت والكربون (100 إلكترون فولت – 5 كيلو إلكترون فولت): رسمت الدراسة خريطة للتنافس بين التيارات البارامغناطيسية والديامغناطيسية. وقد أعادت إنتاج العبور غير الاضطرابي حيث يهيمن التيار البارامغناطيسي الرنيني على الخلفية الديامغناطيسية (Thomson) عند حواف اللب (S L$_{2,3}$ و C K)، بينما يهيمن القناة الديامغناطيسية عند المناطق خارج الرنين. ظلت القناة الديامغناطيسية خلفية مسطحة ($\chi_{dia} = N_e$)، مما أثبت صحة تنفيذ حد $\mathbf{A}^2$.
• توليد التوافق الثاني للماس (1 كيلو إلكترون فولت): في بلورة الماس متمركزة التناظر، يُمنع توليد التوافق الثاني الأمامي ($G=0$) تماماً في حد ثنائي القطب. نجح تنفيذ FMC في توليد إشارة $2\omega_0$ متماسكة تتناسب مع $A_0^2$، والتي كانت غائبة في مرجع ED. كانت الإشارة مهيمنة بالديامغناطيسية (آلية Thomson)، وهو ما يتوافق مع نظرية البلازما الباردة. كما نجحت الطريقة في فصل مساهمات الكتلة (bulk) والسطح عبر مقارنة الصفائح المتناظرة وغير المتناظرة. وبشكل محدد، تجاوزت إشارة FMC مرجع حد ثنائي القطب بعامل قدره ~490 في السعة وللصفيحة المتناظرة، وحوالي ~2200 للمادة الصلبة (bulk).
• CO/O على سطح Ru(0001) عند حافة أكسجين K (~530 إلكترون فولت): على سطح حفاز معدني، قامت الطريقة بفصل التيار التفاعلي المسموح به ثنائياً (والذي طابق أشكال الخطوط التجريبية) عن استجابة الطاقة الامتصاصية. كشف FMC عن تصحيح امتصاصي بحت يتجاوز ثنائي القطب (~14%) في الطاقة الممتصة، ويُعزى ذلك إلى انتقال الرباعي الأقطاب (quadrupole) المسموح به O 1s $\to$ d. هذا المسار ممنوع في حد ثنائي القطب ويصبح متاحاً بسبب الاحتفاظ بالطور المكاني للمجال.

الأهمية والادعاءات
يزعم البحث أنه يضع حجر الأساس لمعالجة ab initio عامة لتفاعلات الضوء والمادة في الأشعة السينية عبر الجزيئات، والأسطح، والمواد الصلبة. من خلال الاحتفاظ بزخم الفوتون على أساس "الكل-إلكترون"، توفر الطريقة إمكانية الوصول إلى:

• تصحيحات "ما وراء ثنائي القطب" (beyond-dipole) القابلة للقياس لثنائيات القطب والامتصاص.
• فتح القنوات الممنوعة في حد ثنائي القطب (مثل الانتقالات رباعية الأقطاب، والتوافق الثاني الممنوع تناظرياً).
• تفكيك ميكانيكي للاستجابات غير الخطية إلى مساهمات بارامغناطيسية وديامغناطيسية محلولة في الفضاء المتبادل.

يشير المؤلفون إلى قيود وتقريبات محددة: يتم التعامل مع اقتران النطاق الجانبي بشكل اضطرابي (فقط $\ell = \pm 1$)، ويقتصر مجموع الشبكة على حد الموقع الواحد ($R=0$)، ولا يتم حل التغذية الراجعة للمجال المستحث بشكل ذاتي الاتساق (على عكس Maxwell–TDDFT المقترن). بناءً على ذلك، فإن الطريقة هي الأنسب للمجالات الضعيفة إلى المتوسطة حيث تكون تأثيرات "ما وراء ثنائي القطب" من الدرجة الأولى هي المهيمنة، ويتم تفسير تيارات $G$ المحدودة عبر نافذة زمنية كأدوات تشخيص ميكانيكية وليس كملاحظات عيانية كاملة للهارمونيات العليا. ورغم هذه التقريبات، تظهر النتائج أن الط फेज المكاني لمجال الأشعة السينية له تبعات فيزيائية ملموسة عبر طيف الأشعة السينية اللينة إلى الصلبة.