डार्विन-फाउलर विधि

सांख्यिकीय यांत्रिकी में, माध्य संभाव्यता के साथ वितरण फ़ंक्शन (भौतिकी) प्राप्त करने के लिए डार्विन-फाउलर विधि का उपयोग किया जाता है। इसे 1922-1923 में चार्ल्स गैल्टन डार्विन और राल्फ एच. फाउलर द्वारा विकसित किया गया था।

वितरण कार्यों का उपयोग सांख्यिकीय भौतिकी में ऊर्जा स्तर पर रहने वाले कणों की औसत संख्या का अनुमान लगाने के लिए किया जाता है (इसलिए इसे व्यवसाय संख्या भी कहा जाता है)। ये वितरण अधिकतर उन संख्याओं के रूप में प्राप्त होते हैं जिनके लिए विचाराधीन प्रणाली अधिकतम संभावना की स्थिति में होती है। किंतु वास्तव में किसी को औसत संख्या की आवश्यकता होती है। ये औसत संख्याएं डार्विन-फाउलर विधि द्वारा प्राप्त की जा सकती हैं। अवश्य ही, सांख्यिकीय यांत्रिकी की तरह, थर्मोडायनामिक सीमा (कणों की बड़ी संख्या) में प्रणाली के लिए, परिणाम अधिकतमकरण के समान ही होते हैं।

डार्विन-फाउलर विधि
सांख्यिकीय यांत्रिकी पर अधिकांश ग्रंथों में सांख्यिकीय वितरण कार्य $$f$$ करता है मैक्सवेल-बोल्ट्ज़मैन आँकड़े, बोस-आइंस्टीन आँकड़े, फ़र्मी-डिराक आँकड़े) उन लोगों को निर्धारित करके प्राप्त किए जाते हैं जिनके लिए प्रणाली अधिकतम संभावना की स्थिति में है। किंतु किसी को वास्तव में औसत या औसत संभावना वाले लोगों की आवश्यकता होती है, चूँकि - निश्चित रूप से - परिणाम सामान्यतः बड़ी संख्या में अवयव वाली प्रणाली के लिए समान होते हैं, जैसा कि सांख्यिकीय यांत्रिकी में होता है। माध्य संभाव्यता के साथ वितरण फलन प्राप्त करने की विधि सी.जी. डार्विन और आर.एच. फाउलर द्वारा विकसित की गई है और इसलिए इसे डार्विन-फाउलर विधि के रूप में जाना जाता है। यह विधि सांख्यिकीय वितरण फलन प्राप्त करने के लिए सबसे विश्वसनीय सामान्य प्रक्रिया है। चूंकि विधि चयनकर्ता चर (गिनती प्रक्रिया की अनुमति देने के लिए प्रत्येक अवयव के लिए प्रस्तुत किया गया कारक) को नियोजित करती है, इसलिए विधि को चयनकर्ता चर की डार्विन-फाउलर विधि के रूप में भी जाना जाता है। ध्यान दें कि वितरण फ़ंक्शन प्रायिकता के समान नहीं है - सीएफ। मैक्सवेल-बोल्ट्ज़मैन वितरण, बोस-आइंस्टीन वितरण, फर्मी-डिराक वितरण। यह भी ध्यान दें कि वितरण फ़ंक्शन $$f_i$$ जो उन अवस्थाओं के अंश का माप है जो वास्तव में अवयव द्वारा व्याप्त हैं, $$f_i = n_i/g_i $$ या $$ n_i= f_ig_i$$ द्वारा दिया गया है, जहाँ $$g_i$$ ऊर्जा स्तर $$i$$ की गिरावट है उर्जा से $$\varepsilon_i$$ और $$n_i$$ इस स्तर पर उपस्थित अवयव की संख्या है (उदाहरण के लिए फर्मी-डिराक आंकड़ों में 0 या 1)। कुल ऊर्जा $$E$$ और अवयव की कुल संख्या $$N$$ को फिर $$E = \sum_i n_i\varepsilon_i$$ और $$N = \sum n_i$$ द्वारा दिए जाते हैं।

डार्विन-फाउलर पद्धति का उपचार ई. श्रोडिंगर, फाउलर और फाउलर और ई. ए गुगेनहेम, के. हुआंग, और एच. जे. डब्ल्यू. मुलर-कर्स्टनके ग्रंथों में किया गया है। आर.बी. डिंगल की पुस्तक में बोस-आइंस्टीन संघनन की व्युत्पत्ति के लिए इस विधि पर भी चर्चा की गई है और इसका उपयोग किया गया है।

मौलिक आँकड़े
$$N=\sum_in_i$$ स्वतंत्र अवयव के लिए $$n_i$$ ऊर्जा $$\varepsilon_i$$ के स्तर पर और $$E=\sum_in_i\varepsilon_i$$ विहित प्रणाली के लिए तापमान $$T$$ के साथ ताप स्नान हम निर्धारित करते हैं
 * $$ Z = \sum_\text{arrangements}e^{-E/kT} = \sum_\text{arrangements}\prod_iz_i^{n_i}, \;\;\; z_i = e^{-\varepsilon_i/kT}.$$

सभी व्यवस्थाओं का औसत, माध्य व्यवसाय संख्या है
 * $$ (n_i)_\text{av} = \frac{\sum_jn_jZ}{Z} = z_j\frac{\partial}{\partial z_j}\ln Z. $$

एक चयनकर्ता चर $$\omega$$ व्यवस्थित करके सम्मिलित करें
 * $$ Z_\omega = \sum \prod_i(\omega z_i)^{n_i}.$$

मौलिक सांख्यिकी में $$N$$ अवयव (a) अलग-अलग हैं और इन्हें पैकेट $$n_i$$ के साथ व्यवस्थित किया जा सकता है स्तर $$\varepsilon_i$$ पर अवयव जिनकी संख्या है
 * $$ \frac{N!}{\prod_in_i!}, $$

जिससे इस स्थिति में
 * $$ Z_\omega = N!\sum_{n_i}\prod_i\frac{(\omega z_i)^{n_i}}{n_i!}.$$

(बी) स्तर $$\varepsilon_i$$ की अधोगति $$g_i$$ के लिए अनुमति देते हुए यह अभिव्यक्ति बन जाती है
 * $$ Z_\omega = N!\prod_{i=1}^{\infty}\left(\sum_{n_i=0,1,2,\ldots}\frac{(\omega z_i)^{n_i}}{n_i!}\right)^{g_i} = N!e^{\omega\sum_ig_iz_i}. $$

चयनकर्ता चर $$\omega$$ किसी को $$\omega^N$$ गुणांक निकालने की अनुमति देता है जो की $$Z$$ है। इस प्रकार
 * $$ Z = \left(\sum_ig_iz_i\right)^N, $$

और इसलिए
 * $$ (n_j)_\text{av} = z_j\frac{\partial}{\partial z_j}\ln Z = N\frac{g_je^{-\varepsilon_j/kT}}{\sum_ig_ie^{-\varepsilon_i/kT}}.$$

यह परिणाम जो अधिकतमीकरण द्वारा प्राप्त सबसे संभावित मूल्य से सहमत है, इसमें एक भी सन्निकटन सम्मिलित नहीं है और इसलिए यह स्पष्ट है, और इस प्रकार इस डार्विन-फाउलर विधि की शक्ति को प्रदर्शित करता है।

क्वांटम आँकड़े
हमारे पास उपरोक्तानुसार है
 * $$ Z_{\omega}=\sum\prod (\omega z_i)^{n_i}, \;\; z_i=e^{-\varepsilon_i/kT}, $$

जहाँ $$n_i$$ ऊर्जा स्तर $$\varepsilon_i$$ में अवयव की संख्या है। चूंकि क्वांटम सांख्यिकी में अवयव अप्रभेद्य हैं, इसलिए अवयव को पैकेटों में विभाजित करने की विधियों की संख्या की कोई प्रारंभिक गणना नहीं की गई है $$n_1, n_2, n_3, ...$$ आवश्यक है। इसलिए योग $$\sum$$ केवल संभावित मानों $$n_i$$ के योग को संदर्भित करता है।

फर्मी-डिराक आँकड़ों की स्थिति में हमारे पास है
 * $$ n_i=0$$ या $$ n_i=1 $$

प्रति अवस्था. ऊर्जा स्तर $$g_i$$ ऊर्जा स्तर के लिए $$\varepsilon_i$$ स्थितियाँ हैं। इसलिए हमारे पास है
 * $$ Z_\omega=(1+\omega z_1)^{g_1}(1+\omega z_2)^{g_2}\cdots=\prod(1+\omega z_i)^{g_i}.$$

बोस-आइंस्टीन सांख्यिकी की स्थिति में हमारे पास है
 * $$ n_i=0,1,2,3, \ldots \infty.$$

पहले जैसी ही प्रक्रिया से हम वर्तमान स्थिति में प्राप्त करते हैं
 * $$Z_{\omega}=(1+\omega z_1+(\omega z_1)^2 + (\omega z_1)^3 + \cdots)^{g_1}(1+\omega z_2 + (\omega z_2)^2 + \cdots)^{g_2} \cdots.$$

किंतु
 * $$1 + \omega z_1 + (\omega z_1)^2 + \cdots = \frac{1}{(1 - \omega z_1)}.$$

इसलिए
 * $$ Z_\omega=\prod_i(1-\omega z_i)^{-g_i}.$$

दोनों स्थितियों को सारांशित करना और $$Z$$ की परिभाषा को याद करते हुए, हम पाते हैं कि $$Z$$ में $$\omega^N$$ का गुणांक है
 * $$ Z_\omega=\prod_i(1\pm \omega z_i)^{\pm g_i},$$

जहां ऊपरी संकेत फर्मी-डिराक सांख्यिकी पर प्रयुक्त होते हैं, और निचले संकेत बोस-आइंस्टीन सांख्यिकी पर प्रयुक्त होते हैं।

आगे हमें फ़ंक्शन $$\omega^N$$ के स्थिति में $$Z_\omega.$$ में $$\phi(\omega)$$ के गुणांक का मूल्यांकन करना होगा जिसे इस प्रकार विस्तारित किया जा सकता है
 * $$ \phi(\omega) = a_0 + a_1\omega + a_2\omega^2 + \cdots, $$

$$\omega^N$$ का गुणांक कॉची के अवशेष प्रमेय की सहायता से है,
 * $$ a_N = \frac{1}{2\pi i}\oint \frac{\phi(\omega)d\omega}{\omega^{N+1}}.$$

हम ध्यान दें कि इसी प्रकार गुणांक $$Z$$ उपरोक्त के रूप में प्राप्त किया जा सकता है
 * $$ Z=\frac{1}{2\pi i}\oint\frac{Z_{\omega}}{\omega^{N+1}}d\omega\equiv \frac{1}{2\pi i}\int e^{f(\omega)}d\omega,$$

जहाँ
 * $$ f(\omega)=\pm\sum_ig_i\ln (1\pm \omega z_i)-(N+1)\ln\omega.$$

अंतर करने से एक प्राप्त होता है
 * $$ f'(\omega) = \frac{1}{\omega}\left[\sum_i\frac{g_i}{(\omega z_i)^{-1}\pm 1}-(N+1)\right],$$

और
 * $$ f''(\omega) = \frac{N+1}{\omega^2}\mp \frac{1}{\omega^2}\sum_i\frac{g_i}{[(\omega z_i)^{-1}\pm 1]^2}.$$

अब कोई स्थिर बिंदु $$f(\omega)$$ पर $$\omega_0$$ के पहले और दूसरे डेरिवेटिव का मूल्यांकन करता है जिस पर $$f'(\omega_0)=0.$$ सैडल बिंदु $$Z$$ के आसपास $$\omega_0$$ के मूल्यांकन की इस पद्धति को तीव्रतम अवतरण की विधि के रूप में जाना जाता है। तब कोई एक प्राप्त करता है
 * $$ Z = \frac{e^{f(\omega_0)}}{\sqrt{2\pi f''(\omega_0)}}.$$

हमारे पास $$f'(\omega_0) = 0 $$ है और इसलिए
 * $$(N+1) = \sum_i\frac{g_i}{(\omega_0z_i)^{-1}\pm 1}$$

(तब से +1 नगण्य है $$N$$ बड़ी है)। हम क्षण में देखेंगे कि यह अंतिम संबंध केवल सूत्र है
 * $$ N = \sum_in_i.$$

हमें मूल्यांकन करके माध्य व्यवसाय संख्या $$(n_i)_{av}$$ प्राप्त होती है
 * $$(n_j)_{av} = z_j\frac{d}{dz_j}\ln Z = \frac{g_j}{(\omega_0z_j)^{-1}\pm 1} = \frac{g_j}{e^{(\varepsilon_j-\mu)/kT} \pm 1}, \quad e^{\mu/kT}= \omega_0.$$

यह अभिव्यक्ति आयतन $$N$$ में कुल $$V$$ के अवयव की औसत संख्या देती है जो तापमान $$T$$ पर 1-कण स्तर $$\varepsilon_j$$ पर अवनति $$g_j$$ के साथ व्याप्त है (उदाहरण के लिए एक प्राथमिक संभावना देखें)। संबंध के विश्वसनीय होने के लिए यह जांचना चाहिए कि उच्च क्रम के योगदान प्रारंभ में परिमाण में कम हो रहे हैं जिससे सैडल बिंदु के आसपास का विस्तार वास्तव में एक स्पर्शोन्मुख विस्तार उत्पन्न कर सके।