Explicit Two-Sided Eigenvalue Bounds for Schrödinger Operators with Singular Potentials via Finite Element Method

ملخص تقني: حدود القيم الذاتية الصريحة من الجانبين لمؤثرات شرودنغر ذات الجهود المنفردة عبر طريقة العناصر المحدودة

بيان المشكلة
تتناول الورقة البحثية حساب حدود صريحة وقابلة للحساب من الجانبين للقيم الذاتية $\lambda_k$ لمؤثر شرودنغر $H = -\Delta + V$ على نطاقات غير محدودة $\mathbb{R}^N$ ($N=2, 3$). يكمن التحدي الرئيسي في التعامل مع الجهود الجاذبة المنفردة من نوع كولوم، وتحديداً $V(x) = -Z/|x|$ (الذي يمثل ذرة الهيدروجين)، والجهود ثنائية المركز $V(x) = -Z_1/|x-a_1| - Z_2/|x-a_2|$ (الذي يمثل أيون الهيدروجين الموجب $H_2^+$).

توفر الطرق العددية الحالية عادةً حدوداً عليا باستخدام طريقة غاليركين المتوافقة (رايلي-ريتز). أما الحصول على حدود دنيا صارمة فهو أكثر صعوبة بكثير، لا سيما في النطاقات غير المحدودة والجهود المنفردة حيث يكون الجهد غير محدود من الأسفل. تفشل تقنيات الحدود الدنيا الكلاسيكية (مثل ليهمان-مايلي-غوريش، وكاتو-تمبل) غالباً لأنها تتطلب طرقاً غير متوافقة أو شروطاً محددة تفشل عندما يكون الجهد $V$ سالباً ومنفرداً. علاوة على ذلك، تعتمد طرق العناصر المحدودة القياسية للحدود الدنيا غالباً على إزاحة الجهد بمقدار معيار $L^\infty$ لجزءه السالب ($\gamma_h = \|c_h^-\|_{L^\infty}$). وبالنسبة لانفرادات كولوم، تتباعد $\gamma_h$ مع تنعيم الشبكة ($O(h^{-1})$)، مما يؤدي إلى فشل الحد الأدنى في التقارب.

المنهجية
تجمع الخوارزمية المقترحة بين قطع النطاق وتوسيع لطريقة العناصر المحدودة من نوع "كروزي-ريفيار المركبة والمُثرية" (CECR). وتتكون هذه المنهجية من ثلاثة مكونات رئيسية:

1. قطع النطاق: يتم قطع النطاق غير المحدود $\mathbb{R}^N$ إلى نطاق محدود $D(R)$.

   • الحدود العليا: يتم الحصول عليها عبر شروط حدودية ديريكليه ($u=0$ على $\partial D(R)$) باستخدام عناصر $P_1$ غاليركين المتوافقة.
   • الحدود الدنيا: يتم الحصول عليها عبر شروط حدودية نيومان ($\partial u/\partial n = 0$). تثبت لمة (Lemma) أن إذا تجاوز أدنى قيمة للجهد الخارجي $\sigma_{ext}(D(R))$ القيمة الذاتية $k$-الرابعة، فإن قيمة نيومان الذاتية على $D(R)$ توفر حداً أدنى صالحاً للقيمة الذاتية في الفضاء الكامل. ويتم التحقق من هذا الشرط لاحقاً (a posteriori) باستخدام الحد العلوي المحسوب.

2. إطار عمل CECR للجهود متغيرة الإشارة:

   • تستخدم الطريقة إطار عمل فضاء الزوج (pair-space framework) حيث تُرفع الدوال إلى أزواج $(u, \Pi_{0,h}u)$، مما يفصل طاقة الحركة (التدرج) عن طاقة التفاعل (الجهد).
   • للتعامل مع الانفراد دون إزاحات متباعدة، قدم المؤلفون ثابت تفاعل لكل عنصر $\Gamma_h := \max_K c_{h|K}^- \cdot h_K^2 / \pi^2$. وبخلاف الإزاحة العالمية $\gamma_h$، فإن $\Gamma_h$ يتلاشى ($O(h^2)$ أو $O(h)$ اعتماداً على تدرج الشبكة) مع تنعيم الشبكة، حتى بالنسبة لانفرادات كولوم.
   • يتم تطبيق إزاحة ثابتة $\sigma > 0$ (مستقلة عن حجم الشبكة $h$) على الجهد لضمان تماسك الصيغة الثنائية الخطية المزاحة. وهذا يتجنب الحاجة إلى إزاحات تعتمد على الشبكة وتتباعد.

3. تحليل الاضطراب:

   • يتم تقريب الجهد المستمر $V$ بمتوسط عناصره الثابت لكل قطعة $c_h$.
   • تحدد لمة الاضطراب الخطأ الناتج عن استبدال $V$ بـ $c_h$. وقد ثبت أن نسبة الاضطراب $\epsilon_h$ تتقارب إلى الصفر ($O(h^2)$ في الشبكات المتدرجة) تحت افتراضات محددة لتدرج الشبكة بالقرب من الانفرادات.

المساهمات الرئيسية
تدعي الورقة البحثية المساهمات المحددة التالية:

• أول حدود دنيا صريحة ومتقاربة: تقدم أول خوارزمية عددية توفر حدوداً من الجانبين، صريحة ومتقاربة، لمؤثرات شرودنغر ذات جهود منفردة غير محدودة على نطاقات غير محدودة.
• مبرهنة الفضاء الزوجي المجردة: طوروا إطار عمل فضاء الزوج لرفع تحليل CECR إلى فضاء مكون من عنصرين، مما يسمح بإثبات ذاتي لـ "الحد الأقصى-الأدنى" (max-min) لمبرهنة الحد الأدنى للجهود ذات الإشارة غير المحددة.
• ثابت التفاعل لكل عنصر ($\Gamma_h$): قدموا $\Gamma_h$ ليحل محل الإزاحة المتباعدة لـ $L^\infty$. وهذا يسمح لفجوة CECR بالانغلاق بمعدل $O(h^2)$ على الشبكات المتدرجة، متجاوزين الحاجة إلى ثوابت معيار الشكل غير القابلة للحساب.
• الحد الثنائي الجانب ذو الإزاحة الثابتة: يجمع العمل حصرًا كاملاً قابلاً للحساب $L_k \le \lambda_k \le U_k$ لا يتطلب التحكم في $L^\infty$ للجزء السالب من الجهد.
• التطبيق على $H_2^+$: خصصوا النظرية لجهد كولوم ثنائي المركز، موفرين صيغاً صريحة للتحويلات ومعاملات الاضطراب في كل من 2D و 3D.

النتائج العددية
تم التحقق من صحة الخوارزمية عبر أربع حالات اختبار:

1. ذرة الهيدروجين في 2D ($V = -1/|x|$): نجحت الطريقة في حصر الحالة الأرضية الدقيقة $\lambda_1 = -1$. وتضيق فترة الحد الأدنى المعتمد بشكل ملحوظ مع تنعيم الشبكة، حيث حققت حدود الإزاحة المثلى عرضاً قدره $0.002$ عند أدنى مستوى.
2. أيون $H_2^+$ في 2D: تظهر نتائج جهد كولوم ثنائي المركز التقارب، حيث تضيق الفجوة المعتمدة للحالة الأرضية من $[-4.43, -1.317]$ إلى $[-2.36, -1.317]$ مع تنعيم الشبكة.
3. ذرة الهيدروجين في 3D: على نطاق كروي، تحصر الطريقة $\lambda_1 = -0.25$. ويتقارب عرض حصر نيومان المعتمد بمعدل $O(h^2)$.
4. أيون $H_2^+$ في 3D: بالنسبة للأيون الجزيئي على نطاق صندوقي، توفر الطريقة فجوة معتمدة $[-1.045, -0.544]$ (عرض $0.501$) للحالة الأرضية، والتي تحتوي على القيمة المرجعية$\approx -0.5513$.

في جميع الحالات، لوحظ أن فجوة CECR الدنيا $(\mu_{k,h}^\sigma - \sigma) - L_{k}^{\mu,\sigma}$ تنغلق بالمعدل المتوقع $O(h^2)$ على الشبكات المتدرجة. وتشير الورقة إلى أنه بينما تكون الفجوات في 3D أوسع منها في 2D بسبب التلاشي الأبطأ للثابت لكل عنصر في شبكات رباعي الأوجه الخشنة، إلا أن سلوك التقارب متسق مع النظرية.

الأهمية والادعاءات
تؤكد الورقة أن هذا العمل يشكل أول إطار عمل للعناصر المحددة قادر على توفير حدود دنيا صريحة ومتقاربة للقيم الذاتية المرتبطة بجهود منفردة غير محدودة. تكمن الأهمية في التغلب على "فجوة التماسك" الناتجة عن الجهود الجاذبة المنفردة، والتي تسببت سابقاً في تباعد تقنيات الحدود الدنيا القياسية.

يؤكد المؤلفون أن الطريقة صريحة: فجميع الثوابت (ثوابت الاستيفاء، حدود الاضطراب، معاملات الإزاحة) مشتقة في صيغة مغلقة من بيانات الشبكة والجهد. وبينما يستخدم التنفيذ الحالي الحسابات بدقة مزدوجة (double-precision)، يذكر المؤلفون صراحة أن التحويل إلى حصر معتمد بالكامل (آمن وفق معيار IEEE-754 وحساب الفترات) هو خطوة تحقق منفصلة تركت للعمل المستقبلي. كما تشير الورقة إلى أنه بينما تعتمد هذه الحدود على الإسقاط، إلا أنها تعمل كنقطة انطلاق صارمة ($\rho$) لتقنيات المعالجة اللاحقة الأكثر حدة مثل حدود ليهمان-غوريش، والتي يمكن أن تضيق الفترات بشكل أكبر.

لا تدعي الورقة حل مشكلة متعددة الإلكترونات مباشرة، بل تحدد امتداد النظرية إلى الهاملتونيان المغناطيسي والأنظمة متعددة الإلكترونات مع قيود التماثل كاتجاهات مستقبلية طبيعية.