Thomas-Fermi equation revisited: A variation on a theme by Majorana

ملخص تقني لورقة "إعادة النظر في معادلة توماس-فيرمي: تنويع على موضوع مايورانا" لبرتولد-جورج إنغلرت.

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

تتناول الورقة معادلة توماس-فيرمي (TF)، وهي معادلة تفاضلية غير خطية من الدرجة الثانية، تُعد مركزية في النموذج الإحصائي للذرات متعددة الإلكترونات. وبينما تم حل هذه المعادلة عددياً في الثمانينيات (تحديداً في مؤلف إنغلرت وشوينجر، ١٩٨٨)، إلا أن تلك الحسابات كانت مرهقة حاسوبياً واعتمدت على متسلسلات قوى ذات فترات تقارب محدودة وغير متداخلة.

يسعى المؤلف إلى:
١. إعادة زيارة وتفصيل طريقة اكتشفها إيتوري مايورانا في الثلاثينيات (نُشرت بعد وفاته) تستغل خصائص القياس (scaling properties) لمعادلة توماس-فيرمي لتحويل المعادلة التفاضلية من الدرجة الثانية إلى معادلة تفاضلية مكافئة من الدرجة الأولى.
٢. تطبيق هذه الطريقة على حالتين فيزيائيتين محددتين:

• الحالة (i): الذرات المتعادلة ($f(0)=1, f(\infty)=0$).
• الحالة (ii): الذرات المتأينة جزئياً (ضعيفة التأين) ($f(1)=0, f'(1)=-\Lambda^2$).
  ٣. إعادة حساب الثوابت والتكاملات الفيزيائية الرئيسية بدقة عالية للتحقق من النتائج العددية السابقة وتحسينها.

٢. المنهجية

أ. تحويل القياس والدوال الثابتة قياسياً

تُعطى معادلة توماس-فيرمي بالصيغة $f''(x) = x^{-1/2}f(x)^{3/2}$. تمتلك المعادلة تناظراً في القياس: إذا كان $f(x)$ حلاً، فإن $f_\lambda(x) = \lambda^3 f(\lambda x)$ هو أيضاً حل.

تتمثل رؤية مايورانا في بناء تركيبات ثابتة القياس للحل $f(x)$ ومشتقاته. يعرّف إنغلرت ثلاث دوال من هذا النوع:

• $P(x) = x^{3/2}f(x)^{1/2}$
• $Q(x) = -x^4 f'(x)$
• $R(x) = -f(x)^{-4/3}f'(x)$

تخضع هذه الدوال لعلاقات تفاضلية محددة. ومن خلال اختيار الزوج المناسب بناءً على الشروط الحدية، يتم اختزال المسألة من الدرجة الثانية إلى معادلة تفاضلية عادية (ODE) من الدرجة الأولى.

ب. الحالة (i): الذرات المتعادلة

بالنسبة للذرات المتعادلة، الشروط الحدية هي $f(0)=1$ و $f(\infty)=0$.

• التحويل: يسعى المؤلف لإيجاد $R$ كدالة في $P$. ومن خلال استبدال المتغيرات ($P = 12t^3$ و $R = (3/16)^{1/3}u(t)$)، تتقلص المنظومة إلى معادلة مايورانا:
  $$\frac{du}{dt} = -\frac{81 - t u(t)^2}{1 - t^2 u(t)}$$
  مع الشروط الحدية $u(1)=1$ و $u(0) = (16/3)^{1/3}B$، حيث $B$ هو ميل دالة توماس-فيرمي عند نقطة الأصل.
• إعادة البناء: يتم استعادة دالة توماس-فيرمي الأصلية $F(x)$ ومشتقتها عبر تكاملات تتضمن $u(t)$. والأهم من ذلك، أن المتغير $t$ يتراوح من $0$إلى$1$، مما يغطي النطاق الفيزيائي الكامل$x \in [0, \infty)$.

ج. الحالة (ii): الذرات المتأينة جزئياً

بالنسبة للذرات المتأينة جزئياً، يتلاشى الحل عند نصف قطر محدود $x=1$.

• التحويل: يسعى المؤلف لإيجاد $Q$ كدالة في $P$. وباستخدام الاستبدالات ($P = 12s$ و $Q = 432v(s)$)، يتم اشتقاق معادلة من الدرجة الأولى مشابهة:
  $$\frac{dv}{ds} = -\frac{8}{3} \frac{sv(s) - s^4}{v(s) - s^2}$$
  مع الشروط الحدية $v(1)=1$ و $v(0) = 1/(432\Lambda^2)$.

د. التنفيذ العددي عبر متسلسلة القوى

الميزة الرئيسية لهذه الطريقة هي خصائص التقارب للحلول $u(t)$ و $v(s)$.

• توسيع المتسلسلة: يقوم المؤلف بتوسيع $u(t)$ و $v(s)$ في قوى لـ $(1-t)$ و $(1-s)$ على التوالي.
• التقارب: على عكس متسلسلة القوى المباشرة لـ $F(x)$، التي لها فترات تقارب محدودة وغير متداخلة (تتطلب المطابقة عند نقاط متوسطة)، فإن متسلسلات $u(t)$ و $v(s)$ تتقارب عالمياً للنطاق الكامل $0 \le t, s \le 1$.
• علاقات التكرار: يتم تحديد معاملات هذه المتسلسلات عبر علاقات تكرارية مشتقة من المعادلات التفاضلية. وهذا يسمح بحساب قيم الثوابت الفيزيائية بدقة عالية.

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

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

• الميل $B$ للذرات المتعادلة: $1.58807102261...$
• المعلمة $\Lambda$ للذرات المتأينة: $32.729416116173...$
• المعلمات التقاربية $\beta$ و $\alpha$.
  ٤. تقييم التكاملات: تسهل الطريقة حساب مختلف التكاملات التي تتضمن قوى $F(x)$ و $x$، وهي ضرورية لصيغ الطاقة. وتشمل هذه التكاملات الأشكال $\int x^k F(x)^l dx$.

٤. النتائج

• التحقق من القيم السابقة: تطابق القيم التي تم الحصول عليها لـ $B$ و $\Lambda$ و $\beta$ و $\alpha$ مع الأرقام الموجودة في أدبيات الثمانينيات، مما يؤكد دقة الطرق القديمة الأكثر كلفة حاسوبياً.
• دقة جديدة: يقدم المؤلف توسعات عشرية ممتدة لهذه الثوابت والكميات المشتقة (مثل معاملات طاقة الربط).
• قيم التكاملات: يسرد الجدول ٤ في الورقة التكاملات المعاد حسابها (مثل $\int F(x) dx$ و $\int F(x)^2 dx$). وتؤكد النتائج القيم السابقة ولكن بدقة أعلى ودون عدم الاستقرار العددي المرتبط بالطرق القديمة.
• صيغ الطاقة: استُخدمت التكاملات المعاد حسابها لتنقيح المعاملات لـ:
  • طاقة الربط: $E \approx -c_7 Z^{7/3} + \frac{1}{2}Z^2 - c_5 Z^{5/3}$.
  • طاقة التأين: تم حساب الحد الرئيسي للذرات المتأينة جزئياً كـ $I \approx 3.151526$ إلكترون فولت.
  • التصحيحات النسبية: تم تحديث معاملات التصحيحات النسبية ($c_{rel}$) باستخدام قيم التكاملات الجديدة.

٥. الأهمية

• الأناقة الرياضية: تسلط الورقة الضوء على الرؤية العميقة لإيتوري مايورانا، مبينةً كيف يمكن لتماثلات القياس أن تبسط بشكل جذري المعادلات التفاضلية غير الخطية في الفيزياء.
• الكفاءة الحسابية: توفر الطريقة وسيلة "بسيطة" لحساب الأرقام الفيزيائية التي كانت تتطلب سابقاً تكاملاً عددياً معقداً ومطابقة للمتسلسلات. فالتقارب العالمي لمتسلسلة مايورانا يجعل الخوارزمية قوية وسهلة التنفيذ.
• القيمة التعليمية: تعمل الورقة كدليل شامل لتطبيق تقنية مايورانا، وتوسيعها ليشمل نظاماً فيزيائياً جديداً (الذرات المتأينة جزئياً) وتوفير علاقات تكرارية صريحة للحساب العملي.
• التطبيقات المستقبلية: يشير المؤلف إلى إمكانية تطبيق تقنية التحويل هذه على معادلات أخرى من الشكل $f''(x) = x^a f(x)^b$، مثل معادلة ليند-إيمدن لكرات الغاز البولتروبية، وعلى شروط حدود التأين الأكثر تعقيداً.

باختصار، تعيد ورقة إنغلرت إحياء تقنية رياضية تاريخية، مظهرةً فائدتها الفائقة للحسابات الحديثة عالية الدقة في الفيزياء الذرية، وتحديداً لنموذج توماس-فيرمي لكل من الذرات المتعادلة والمتأينة جزئياً.