मल्टीकैनोनिकल एन्सेम्बल

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

प्रेरणा
बड़ी संख्या में स्वतंत्रता की डिग्री वाली प्रणालियों में, जैसे स्पिन (भौतिकी) प्रणालियों में, मोंटे कार्लो एकीकरण की आवश्यकता होती है। इस एकीकरण में, महत्व प्रारूपकरण और विशेष रूप से मेट्रोपोलिस कलन विधि, बहुत ही महत्वपूर्ण प्रौद्योगिकी है। चूँकि, मेट्रोपोलिस कलन विधि प्रारूप के अनुसार बताता है, $$\exp(-\beta E)$$ जहां बीटा तापमान का व्युत्क्रम है। इसका आशय है कि ऊर्जा अवरोध $$\Delta E$$ ऊर्जा स्पेक्ट्रम पर नियंत्रण करना कठिन है। पॉट्स मॉडल जैसे कई स्थानीय ऊर्जा मिनिमा वाले प्रणाली का प्रारूप लेना कठिन हो जाता है, क्योंकि कलन विधि प्रणाली के स्थानीय मिनिमा में फंस जाता है। यह अन्य दृष्टिकोणों, अर्थात् अन्य प्रारूपकरण वितरणों को प्रेरित करता है।

अवलोकन
मल्टीकैनोनिकल एन्सेम्बल मेट्रोपोलिस-हेस्टिंग्स कलन विधि का उपयोग करता है, जो प्रणाली के अवस्था के घनत्व के व्युत्क्रम द्वारा दिए गए प्रारूप वितरण के साथ होता है, जो प्रारूप वितरण के विपरीत है। $$\exp(-\beta E)$$ मेट्रोपोलिस कलन विधि के इस विकल्प के साथ, औसतन, प्रत्येक ऊर्जा पर प्रारूप किए गए अवस्था की संख्या स्थिर होती है, अर्थात यह ऊर्जा पर फ्लैट हिस्टोग्राम के साथ अनुकरण है। यह ऐसे कलन विधि की ओर ले जाता है, जिसके लिए ऊर्जा बाधाओं को दूर करना अब कठिन नहीं है। मेट्रोपोलिस कलन विधि पर अन्य लाभ यह है कि प्रारूपकरण प्रणाली के तापमान से स्वतंत्र है, जिसका अर्थ है कि सिमुलेशन सभी तापमानों के लिए थर्मोडायनामिकल चर के अनुमान की अनुमति प्रदान करता है (इस प्रकार नाम मल्टीकैनोनिकल: कई तापमान)। प्रथम क्रम चरण संक्रमण के अध्ययन में यह बड़ा सुधार है।

मल्टीकैनोनिकल समूह को निष्पादित करने में सबसे बड़ी समस्या यह है कि अवस्था के घनत्व को प्राथमिकता से जानना होगा। मल्टीकैनोनिकल सैंपलिंग में महत्वपूर्ण योगदान वांग और लैंडौ कलन विधि का था, जो अभिसरण तथा अवस्था के घनत्व की गणना करते समय एसिम्प्टोटिक रूप से मल्टीकैनोनिकल समूह में परिवर्तित हो जाता है।

मल्टीकैनोनिकल समूह भौतिक प्रणालियों तक ही सीमित नहीं है। इसे अमूर्त प्रणालियों पर नियोजित किया जा सकता है, जिनमें लागत फलन F होता है। F के संबंध में अवस्था के घनत्व का उपयोग करके, उच्च-आयामी इंटीग्रल की गणना करने या स्थानीय मिनीमा अनुसन्धान के लिए विधि सामान्य हो जाती है।

प्रेरणा
प्रणाली और उसके चरण-स्थान पर विचार करें I $$\Omega$$ विन्यास द्वारा विशेषता $$\boldsymbol{r}$$ में $$\Omega$$ और प्रणाली के चरण-स्थान से एक-आयामी स्थान तक लागत फलन F $$\Gamma$$: $$F(\Omega) = \Gamma = [\Gamma_\min, \Gamma_\max]$$, F का स्पेक्ट्रम है।

औसत मात्रा की गणना $$\langle Q\rangle$$ चरण-स्थान पर अभिन्न के मूल्यांकन की आवश्यकता होती है:


 * $$\langle Q\rangle = \int_\Omega Q(\boldsymbol{r}) P_r(\boldsymbol{r}) \,d\boldsymbol{r}$$

जहाँ $$P_r(\boldsymbol{r})$$ प्रत्येक अवस्था का भार है (उदा. $$P_r(\boldsymbol{r})=1/V$$ समान रूप से वितरित अवस्था के अनुरूप)।

जब Q किसी विशेष अवस्था पर नहीं बल्कि केवल अवस्था के विशेष F के मान पर निर्भर करता है I $$F(\boldsymbol{r}) = F_\boldsymbol{r}$$,के लिए सूत्र $$\langle Q\rangle$$ डायराक डेल्टा फलन जोड़कर F पर एकीकृत किया जा सकता है और इसे इस प्रकार लिखा जा सकता है



\begin{align} \langle Q\rangle & = \int_{\Gamma_\min}^{\Gamma_\max} \int_\Omega Q(F_\boldsymbol{r}) P_r(F_\boldsymbol{r}) \delta(f - F_\boldsymbol{r}) \,d\boldsymbol{r} \,df \\ & = \int_{\Gamma_\min}^{\Gamma_\max} Q(f) \int_\Omega \delta(f - F_\boldsymbol{r}) P_r(F_\boldsymbol{r}) \,d\boldsymbol{r} \, df \\ & = \int_{\Gamma_\min}^{\Gamma_\max} Q(f) P(f) \, df \\ \end{align} $$ जहाँ


 * $$P(f) = \int_\Omega P_r(r)\delta(f - F(\boldsymbol{r})) \,d\boldsymbol{r}$$

F का सीमांत वितरण है I

जब प्रणाली में बड़ी संख्या में स्वतंत्रता की डिग्री होती है, तो इसके लिए विश्लेषणात्मक अभिव्यक्ति होती है I $$\langle Q\rangle$$ प्राप्त करना प्रायः कठिन होता है, और मोंटे कार्लो एकीकरण को सामान्यतः गणना में नियोजित किया जाता है I $$\langle Q\rangle$$ सरलतम फॉर्मूलेशन पर, विधि N समान रूप से वितरित अवस्था का चयन करती है I $$\boldsymbol{r}_i\in \Omega$$, और अनुमानक का उपयोग करता है I


 * $$\overline{Q}_N = \sum_{i = 1}^N Q(\boldsymbol{r}_i) V P_r(\boldsymbol{r}_i)$$

कंप्यूटिंग के लिए $$\langle Q\rangle$$ क्योंकि $$\overline{Q}_N$$ लगभग निश्चित रूप से अभिसरण होता है, $$\langle Q\rangle$$ बड़ी संख्या के नियम द्वारा नियम इस प्रकार है:


 * $$\lim_{N\rightarrow\infty} \overline{Q}_N = \langle Q\rangle.$$

इस अभिसरण की विशिष्ट समस्या यह है कि Q का विचरण बहुत अधिक हो सकता है, जिससे उचित परिणाम प्राप्त करने के लिए उच्च कम्प्यूटेशनल प्रयास करना पड़ता है।

इस अभिसरण को बेहतर बनाने के लिए, मेट्रोपोलिस-हेस्टिंग्स कलन विधि प्रस्तावित किया गया था। आम तौर पर, मोंटे कार्लो पद्धति का विचार अनुमानक के अभिसरण को बेहतर बनाने के लिए महत्व प्रारूप का उपयोग करना है $$\overline{Q}_N$$ एक मनमाना वितरण के अनुसार अवस्था का प्रारूप लेकर $$\pi(\boldsymbol{r})$$, और उपयुक्त अनुमानक का उपयोग करें:


 * $$\overline{Q}_N = \sum_{i = 1}^N Q(\boldsymbol{r}_i) \pi^{-1}(\boldsymbol{r}_i) P_r(\boldsymbol{r}_i)$$.

यह अनुमानक एक मनमाना वितरण से लिए गए नमूनों के माध्य के अनुमानक को सामान्यीकृत करता है। इसलिए, जब $$\pi(\boldsymbol{r})$$ एक समान वितरण है, यह ऊपर एक समान प्रारूप पर उपयोग किए गए वितरण से मेल खाता है।

जब प्रणाली एक भौतिक प्रणाली होता है जो ऊष्मा स्नान के संपर्क में होता है, तो प्रत्येक अवस्था $$\boldsymbol{r}$$ बोल्ट्ज़मान कारक के अनुसार भारित किया जाता है, $$P_r(\boldsymbol{r}_i) \propto \exp(-\beta F_\boldsymbol{r})$$. मोंटे कार्लो में, विहित एन्सेम्बल को चुनकर परिभाषित किया गया है $$\pi(\boldsymbol{r})$$ के आनुपातिक होना $$P_r(\boldsymbol{r}_i)$$. इस स्थिति में, अनुमानक एक साधारण अंकगणितीय औसत से मेल खाता है:
 * $$\overline{Q}_N = \frac{1}{N} \sum_{i = 1}^N Q(\boldsymbol{r}_i)$$

ऐतिहासिक रूप से, यह तेज़ कंप्यूटिंग मशीनों द्वारा अवस्था की गणना के समीकरण के कारण हुआ हीट बाथ के संपर्क में आने वाले प्रणाली पर औसत की गणना करने के लिए मेट्रोपोलिस-हेस्टिंग्स कलन विधि का उपयोग करना था, जहां वजन बोल्ट्जमैन कारक द्वारा दिया जाता है, $$P(\boldsymbol{x}) \propto \exp(-\beta E(\boldsymbol{r}))$$.

जबकि अक्सर ऐसा होता है कि प्रारूप वितरण $$\pi$$ वजन वितरण के लिए चुना गया है $$P_r$$, ऐसा होना आवश्यक नहीं है। एक स्थिति जहां विहित एन्सेम्बल एक कुशल विकल्प नहीं है, वह तब होता है जब इसे एकत्रित होने में मनमाने ढंग से लंबा समय लगता है। एक स्थिति जहां ऐसा होता है जब फलन F में एकाधिक स्थानीय मिनीमा होते हैं। कलन विधि के लिए एक विशिष्ट क्षेत्र को स्थानीय न्यूनतम के साथ छोड़ने की कम्प्यूटेशनल लागत लागत फलन के न्यूनतम मूल्य के साथ तेजी से बढ़ती है। अर्थात्, न्यूनतम जितना गहरा होगा, कलन विधि वहां उतना अधिक समय व्यतीत करेगा, और इसे छोड़ना उतना ही कठिन होगा (स्थानीय न्यूनतम की गहराई के साथ तेजी से बढ़ रहा है)।

लागत फलन के स्थानीय न्यूनतम में फंसने से बचने का एक तरीका प्रारूप प्रौद्योगिकी को स्थानीय न्यूनतम के लिए अदृश्य बनाना है। यह मल्टीकैनोनिकल समूह का आधार है।

मल्टीकैनोनिकल एन्सेम्बल
मल्टीकैनोनिकल समूह को प्रारूप वितरण का चयन करके परिभाषित किया गया है


 * $$\pi(\boldsymbol{r}) \propto \frac{1}{P(F_\boldsymbol{r})}$$

जहाँ $$ P(f)$$ ऊपर परिभाषित F का सीमांत वितरण है। इस विकल्प का परिणाम यह है कि f, m(f) के दिए गए मान के साथ नमूनों की औसत संख्या दी गई है


 * $$ m(f) = \int_\Omega \delta(f - F_\boldsymbol{r}) \pi(\boldsymbol{r}) P_r(\boldsymbol{r})\,d\boldsymbol{r} \propto \int_\Omega \delta(f - F_\boldsymbol{r}) P_r(\boldsymbol{r}) \frac{1}{P(F_\boldsymbol{r})} d\boldsymbol{r} = \frac{1}{P(f)} \int_\Omega \delta(f - F_\boldsymbol{r}) P_r(\boldsymbol{r}) d\boldsymbol{r} = 1$$

अर्थात्, नमूनों की औसत संख्या f पर निर्भर नहीं करती है: सभी लागतें f समान रूप से प्रारूप की जाती हैं, भले ही वे अधिक या कम संभावित हों। यह फ्लैट-हिस्टोग्राम नाम को प्रेरित करता है। हीट बाथ के संपर्क में आने वाली प्रणालियों के लिए, प्रारूपकरण तापमान से स्वतंत्र होता है और एक सिमुलेशन सभी तापमानों का अध्ययन करने की अनुमति देता है।

सुरंग बनाने का समय और महत्वपूर्ण गति धीमी होना
किसी भी अन्य मोंटे कार्लो विधि की तरह, इसमें से लिए गए नमूनों के सहसंबंध होते हैं $$P(\boldsymbol{r})$$. सहसंबंध का एक विशिष्ट माप सुरंग बनाने का समय है। टनलिंग समय को मार्कोव चरणों (मार्कोव श्रृंखला के) की संख्या से परिभाषित किया जाता है, सिमुलेशन को F के न्यूनतम और अधिकतम स्पेक्ट्रम के बीच एक राउंड-ट्रिप करने की आवश्यकता होती है। टनलिंग समय का उपयोग करने के लिए एक प्रेरणा यह है कि जब यह पार हो जाता है स्पेक्ट्रा, यह अवस्था के अधिकतम घनत्व वाले क्षेत्र से होकर गुजरता है, इस प्रकार प्रक्रिया का सहसंबंध समाप्त हो जाता है। दूसरी ओर राउंड-ट्रिप्स का उपयोग यह सुनिश्चित करता है कि प्रणाली सभी स्पेक्ट्रम का दौरा करता है।

क्योंकि हिस्टोग्राम वेरिएबल F पर सपाट है, एक मल्टीकैनोनिक एन्सेम्बल को F मानों की एक-आयामी रेखा पर एक प्रसार प्रक्रिया (यानी एक यादृच्छिक चलना) के रूप में देखा जा सकता है। प्रक्रिया का विस्तृत संतुलन यह निर्देशित करता है कि प्रक्रिया पर कोई स्टोकेस्टिक बहाव नहीं है। इसका तात्पर्य यह है कि सुरंग बनाने का समय, स्थानीय गतिशीलता में, एक प्रसार प्रक्रिया के रूप में स्केल होना चाहिए, और इस प्रकार सुरंग बनाने का समय स्पेक्ट्रम के आकार के साथ चतुष्कोणीय रूप से स्केल होना चाहिए, एन:


 * $$\tau_{tt} \propto N^2$$

चूँकि, कुछ प्रणालियों में (आइज़िंग मॉडल सबसे प्रतिमानात्मक है), स्केलिंग गंभीर रूप से धीमी होने से ग्रस्त है: यह है $$N^{2+z}$$ जहाँ $$z>0$$ विशिष्ट प्रणाली पर निर्भर करता है।

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