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

यहाँ बर्थोल्ड-जॉर्ज एंग्लर्ट के शोध पत्र "थॉमस–फर्मी समीकरण का पुनरावलोकन: मेजरोना द्वारा एक विषय पर भिन्नता" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण

यह शोध पत्र थॉमस–फर्मी (TF) समीकरण को संबोधित करता है, जो कई-इलेक्ट्रॉन परमाणुओं के सांख्यिकीय मॉडल में केंद्रीय भूमिका निभाने वाला एक द्वितीय-क्रम का गैर-रेखीय अवकल समीकरण (second-order non-linear differential equation) है। हालाँकि 1980 के दशक में इस समीकरण को संख्यात्मक रूप से हल किया गया था (विशेष रूप से एंग्लर्ट और श्विंगर, 1988 के मोनोग्राफ में), वे गणनाएँ अत्यधिक श्रमसाध्य थीं और सीमित, गैर-अतिव्यापी अभिसरण अंतरालों (non-overlapping intervals of convergence) वाले पावर सीरीज़ पर निर्भर थीं।

लेखक का लक्ष्य है:

1. एटोर मेजरोना द्वारा 1930 के दशक में खोजे गए (मरणोपरांत प्रकाशित) एक तरीके का पुनरावलोकन और विस्तार करना, जो थॉमस–फर्मी समीकरण के स्केलिंग गुणों का उपयोग करके द्वितीय-क्रम के अवकल समीकरण को एक समतुल्य प्रथम-क्रम के अवकल समीकरण में बदलने का लाभ उठाता है।
2. दो विशिष्ट भौतिक मामलों पर इस विधि को लागू करना:
   • मामला (i): उदासीन परमाणु ($f(0)=1, f(\infty)=0$)।
   • मामला (ii): दुर्बल आयनित परमाणु ($f(1)=0, f'(1)=-\Lambda^2$)।
3. पिछले संख्यात्मक परिणामों को सत्यापित करने और सुधारने के लिए प्रमुख भौतिक स्थिरांकों और समाकलों (integrals) की उच्च सटीकता के साथ पुनर्गणना करना।

2. कार्यप्रणाली

A. स्केलिंग रूपांतरण और स्केल-इनवेरिएंट फलन

TF समीकरण $f''(x) = x^{-1/2}f(x)^{3/2}$ द्वारा दिया गया है। इस समीकरण में एक स्केलिंग समरूपता (scaling symmetry) है: यदि $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)$

ये फलन विशिष्ट अवकल संबंधों का पालन करते हैं। सीमा शर्तों (boundary conditions) के आधार पर उपयुक्त युग्म का चयन करके, द्वितीय-क्रम की समस्या को प्रथम-क्रम के साधारण अवकल समीकरण (ODE) में बदल दिया जाता है।

B. मामला (i): उदासीन परमाणु

उदासीन परमाणुओं के लिए, सीमा शर्तें $f(0)=1$ और $f(\infty)=0$ हैं।

• रूपांतरण: लेखक $P$ के फलन के रूप में $R$ की खोज करता है। चर प्रतिस्थापन ($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$ मूल बिंदु (origin) पर TF फलन का ढाल (slope) है।
• पुनर्गठन (Reconstruction): मूल TF फलन $F(x)$ और उसके डेरिवेटिव को $u(t)$ से जुड़े समाकलों के माध्यम से पुनर्प्राप्त किया जाता है। महत्वपूर्ण रूप से, चर $t$ का मान $0$से$1$तक होता है, जो संपूर्ण भौतिक डोमेन$x \in [0, \infty)$ को कवर करता है।

C. मामला (ii): दुर्बल आयनित परमाणु

दुर्बल आयनित परमाणुओं के लिए, समाधान एक परिमित त्रिज्या $x=1$ पर शून्य हो जाता है।

• रूपांतरण: लेखक $P$ के फलन के रूप में $Q$ की खोज करता है। प्रतिस्थापन ($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)$ हैं।

D. पावर सीरीज़ के माध्यम से संख्यात्मक कार्यान्वयन

इस विधि का मुख्य लाभ $u(t)$ और $v(s)$ के अभिसरण गुण (convergence properties) हैं।

• सीरीज विस्तार: लेखक $(1-t)$ और $(1-s)$ की घातों के रूप में $u(t)$ और $v(s)$ का विस्तार करता है।
• अभिसरण: $F(x)$ के प्रत्यक्ष पावर सीरीज़ के विपरीत, जिनके अभिसरण अंतराल सीमित और गैर-अतिव्यापी होते हैं (जिसके लिए मध्यवर्ती बिंदुओं पर मिलान की आवश्यकता होती है), $u(t)$ और $v(s)$ के लिए सीरीज़ पूरे भौतिक डोमेन $0 \le t, s \le 1$ के लिए वैश्विक (globally) रूप से अभिसरित होती है।
• पुनरावृत्ति संबंध (Recurrence Relations): इन सीरीज़ के गुणांकों को अवकल समीकरणों से प्राप्त पुनरावृत्ति संबंधों द्वारा निर्धारित किया जाता है। यह भौतिक स्थिरांकों के उच्च-सटीकता वाले मानों की गणना करने की अनुमति देता है।

3. प्रमुख योगदान

1. आयनित मामले का व्युत्पन्न: जबकि मेजरोना की विधि उदासीन परमाणुओं के लिए ज्ञात थी, एंग्लर्ट स्पष्ट रूप से दुर्बल आयनित परमाणुओं (मामला ii) के लिए समतुल्य प्रथम-क्रम समीकरण का विश्लेषण और व्युत्पन्न करते हैं, जो सीमित प्रणालियों के लिए प्रासंगिक समाधान है।
2. वैश्विक अभिसरण: शोध पत्र प्रदर्शित करता है कि मेजरोना रूपांतरण समस्या को एक ऐसी समस्या में बदल देता है जिसे पूरे भौतिक डोमेन के लिए अभिसरित होने वाली पावर सीरीज़ द्वारा हल किया जा सकता है, जिससे पिछले संख्यात्मक अध्ययनों में उपयोग की जाने वाली "Tedious" (थकाऊ) मिलान प्रक्रियाओं की आवश्यकता समाप्त हो जाती है।
3. स्थिरांकों की पुनर्गणना: लेखक इस नई सीरीज़ विधि का उपयोग करके मौलिक स्थिरांकों की उच्च सटीकता के साथ पुनर्गणना करता है:
   • उदासीन परमाणुओं के लिए ढाल (Slope) $B$: $1.58807102261...$
   • आयनित परमाणुओं के लिए पैरामीटर $\Lambda$: $32.729416116173...$
   • एसिम्प्टोटिक पैरामीटर $\beta$ और $\alpha$।
4. समाकलों का मूल्यांकन: यह विधि $F(x)$ और $x$ की घातों वाले विभिन्न समाकलों की गणना को सुगम बनाती है, जो ऊर्जा सूत्रों के लिए आवश्यक हैं। इनमें $\int x^k F(x)^l dx$ के रूप के समाकल शामिल हैं।

4. परिणाम

• पिछले मूल्यों का सत्यापन: $B, \Lambda, \beta,$ और $\alpha$ के लिए प्राप्त मान 1980 के दशक के साहित्य में पाए गए अंकों से मेल खाते हैं, जो पुराने, अधिक गणनात्मक रूप से महंगे तरीकों की सटीकता की पुष्टि करते हैं।
• नई सटीकता: लेखक इन स्थिरांकों और व्युत्पन्न मात्राओं (जैसे बाइंडिंग एनर्जी गुणांक) के लिए विस्तारित दशमलव विस्तार प्रदान करता है।
• समाकल मान: शोध पत्र में तालिका 4 में पुनर्गणना किए गए समाकल (जैसे $\int F(x) dx, \int F(x)^2 dx$) सूचीबद्ध हैं। परिणाम पिछले मूल्यों की पुष्टि करते हैं लेकिन उच्च सटीकता के साथ और बिना पुरानी विधियों से जुड़ी संख्यात्मक अस्थिरता के।
• ऊर्जा सूत्र: पुनर्गणना किए गए समाकलों का उपयोग निम्नलिखित के लिए गुणांकों को परिष्कृत करने के लिए किया जाता है:
  • बाइंडिंग एनर्जी (Binding Energy): $E \approx -c_7 Z^{7/3} + \frac{1}{2}Z^2 - c_5 Z^{5/3}$।
  • आयनीकरण ऊर्जा (Ionization Energy): दुर्बल आयनित परमाणुओं के लिए अग्रणी पद की गणना $I \approx 3.151526$ eV के रूप में की गई है।
  • सापेक्षिक सुधार (Relativistic Corrections): सापेक्षिक सुधारों ($c_{rel}$) के लिए गुणांकों को नए समाकल मानों का उपयोग करके अपडेट किया गया है।

5. महत्व

• गणितीय सुंदरता: यह शोध पत्र एटोर मेजरोना की गहन अंतर्दृष्टि को उजागर करता है, जो दिखाता है कि कैसे स्केलिंग समरूपता (scaling symmetries) भौतिकी में जटिल गैर-रेखीय अवकल समीकरणों को काफी सरल बना सकती है।
• गणना दक्षता: यह विधि एक ऐसा "सरल तरीका" प्रदान करती है जिससे उन भौतिक संख्याओं की गणना की जा सकती है जिनके लिए पहले जटिल संख्यात्मक एकीकरण और सीरीज़ मिलान की आवश्यकता होती थी। मेजरोना सीरीज़ का वैश्विक अभिसरण एल्गोरिदम को मजबूत और लागू करने में आसान बनाता है।
• शैक्षणिक मूल्य: यह शोध पत्र मेजरोना की तकनीक को लागू करने के लिए एक व्यापक मार्गदर्शिका के रूप में कार्य करता है, इसे एक नए भौतिक क्षेत्र (आयनित परमाणु) तक विस्तारित करता है और व्यावहारिक गणना के लिए स्पष्ट पुनरावृत्ति संबंध प्रदान करता है।
• भविष्य के अनुप्रयोग: लेखक सुझाव देते हैं कि इस रूपांतरण तकनीक को अन्य गैर-रेखीय समीकरणों, जैसे कि पॉलीट्रोपिक गैस बॉल्स के लिए लेन-एम्डेन समीकरण, और अधिक जटिल आयनीकरण सीमा स्थितियों पर लागू किया जा सकता है।

संक्षेप में, एंग्लर्ट का शोध पत्र एक ऐतिहासिक गणितीय तकनीक को पुनर्जीवित करता है, जो उदासीन और दुर्बल आयनित दोनों परमाणुओं के लिए थॉमस–फर्मी मॉडल हेतु आधुनिक उच्च-सटीकता गणनाओं के लिए इसकी श्रेष्ठ उपयोगिता को प्रदर्शित करता है।