Escape over a saddle by coloured noise: theory and numerics

तकनीकी सारांश: रंगीन शोर (Coloured Noise) द्वारा सैडल के ऊपर से पलायन: सिद्धांत और संख्यात्मक विश्लेषण

समस्या विवरण
स्टोकेस्टिक (stochastic) गतिशील प्रणालियाँ यादृच्छिक विक्षोभों द्वारा प्रेरित संक्रमणों को मॉडल करने के लिए आवश्यक हैं, विशेष रूप से एक सैडल बिंदु (saddle point) के स्थिर मैनिफोल्ड (stable manifold) के पार आकर्षक साम्यावस्था (attracting equilibria) से पलायन के मामले में। जबकि छोटे-शोर (small-noise) सीमा के अंतर्गत फ्राइडलिन-वेंट्ज़ेल लार्ज डेविएशन प्रिंसिपल (Freidlin-Wentzell LDP) के माध्यम से यह सैद्धांतिक रूप से अच्छी तरह से समझा जा चुका है, मौजूदा कम्प्यूटेशनल विधियों को महत्वपूर्ण सीमाओं का सामना करना पड़ता है। वर्तमान दृष्टिकोण आमतौर पर निम्नलिखित समस्याओं से जूझते हैं:

1. डीजेनरेट शोर (Degenerate noise): ऐसी स्थितियाँ जहाँ शोर पूरे फेज़ स्पेस (phase space) की दिशाओं में कार्य नहीं करता है (जैसे, यांत्रिक प्रणालियों में जहाँ शोर केवल मोमेंटम समीकरणों के माध्यम से प्रवेश करता है)।
2. रंगीन शोर (Coloured noise): परिमित सहसंबंध समय (finite correlation times) वाला यथार्थवादी बल, जिसे अक्सर आदर्श श्वेत शोर (white noise) के बजाय श्वेत शोर पर रैखिक फिल्टर के माध्यम से मॉडल किया जाता है।
3. अनबाउंडेड समय-क्षितिज (Unbounded time-horizons): यह भौतिक वास्तविकता कि पलायन की घटनाएँ अनंत समय अंतराल पर होती हैं, जबकि संख्यात्मक विधियाँ अक्सर परिमित अंतरालों पर निर्भर करती हैं, जिससे अक्षमता और समय-अनुवाद विसंगतियाँ (time-translation degeneracies) उत्पन्न होती हैं।

अनंत समय-क्षितिज पर बैरियर ऊँचाई (barrier heights) की गणना करने वाली मौजूदा विधियाँ शोर सहप्रसरण मैट्रिक्स (noise covariance matrix) की व्युत्क्रमणीयता (invertibility) पर निर्भर करती हैं, जो उन्हें डीजेनरेट शोर मामलों के लिए अनुपयुक्त बनाती है। इसके अतिरिक्त, दुर्लभ घटनाओं का सीधा नमूनाकरण (direct sampling) कम्प्यूटेशनल रूप से अत्यधिक खर्चीला है, और अनंत-समय प्रक्षेप पथों (trajectories) का मानक विविक्तीकरण (discretization) उन क्षेत्रों में संसाधनों को बर्बाद करता है जहाँ गतिकी लगभग रैखिक होती है।

कार्यप्रणाली (Methodology)
लेखक इन चुनौतियों का समाधान करने के लिए तीन विशिष्ट पद्धतिगत नवाचारों को संयोजित करने वाला एक एकीकृत ढांचा प्रस्तावित करते हैं:

1. डीजेनरेट और रंगीन शोर के लिए हैमिल्टनियन इष्टतम नियंत्रण (Hamiltonian Optimal Control - HOC):

   • डीजेनरेट शोर (जहाँ सहप्रसरण मैट्रिक्स $C = \sigma\sigma^T$ विलुप्त/सिंगुलर है) को संभालने के लिए, यह शोध पत्र हैमिल्टनियन इष्टतम नियंत्रण विधि को अपनाता है। एक्शन फंक्शनल (action functional) जिसमें $C^{-1}$ शामिल है (जो डीजेनरेट शोर के लिए अपरिभाषित है) को न्यूनतम करने के बजाय, इस समस्या को एक विस्तारित फेज़ स्पेस में एक नियतकालिक परिवर्तनशील कार्य (deterministic variational task) के रूप में पुनर्गठित किया गया है।
   • प्रत्येक अवस्था घटक के लिए संयुग्मी चरों (conjugate variables/adjoints) को पेश करके, यह विधि हैमिल्टन के गति के समीकरणों को व्युत्पन्न करती है। यह दृष्टिकोण शोर सहप्रसरण मैट्रिक्स के व्युत्क्रमण (inversion) से पूरी तरह बचता है, जो बिना किसी स्पष्ट प्री-कंडीशनिंग या सख्त हाइपोएलिप्टिसिटी (hypoellipticity) धारणा के स्वाभाविक रूप से डीजेनरेट शोर को समायोजित करता है।
   • रंगीन शोर को फिल्टर चरों (filter variables) को विस्तारित फेज़ स्पेस में जोड़कर, जो ऑर्नस्टीन-उलेनबैक (Ornstein-Uhlenbeck) प्रक्रियाओं द्वारा नियंत्रित होते हैं, मॉडल किया जाता है। यह रंगीन शोर की समस्या को उच्च-आयामी अवस्था वेक्टर वाले एक मानक SDE ढांचे में बदल देता है, जहाँ शोर केवल फिल्टर चरों पर कार्य करता है, जिससे मूल प्रणाली के लिए स्वतः ही एक डीजेनरेट शोर संरचना प्राप्त होती है।

2. विभाजन की विधि (Method of Division - MOD) अनंत समय-क्षितिज के लिए:

   • अनंत समय-अंतरालों के अनुकरण की कम्प्यूटेशनल अव्यवहार्यता और समय-अनुवाद विसंगति के कारण होने वाली संख्यात्मक अस्थिरता को संबोधित करने के लिए, लेखक विभाजन की विधि प्रस्तुत करते हैं।
   • संक्रमण प्रक्षेप पथ (transition trajectory) को तीन अलग-अलग चरणों में विभाजित किया गया है:
     • प्रारंभिक पलायन (Initial Escape): प्रणाली $(-\infty, 0]$ के दौरान स्थिर साम्यावस्था (sink) से बाहर निकलती है। इस चरण को सिंक के आसपास रैखिक गतिकी का उपयोग करके विश्लेषणात्मक रूप से माना जाता है। एक्शन की गणना कंट्रोलेबिलिटी ग्रैमियन (controllability Gramian) के माध्यम से की जाती है, जो एक लियापुनोव समीकरण (Lyapunov equation) को हल करता है।
     • संक्रमण चरण (Transitional Phase): मुख्य गैर-रैखिक संक्रमण एक परिमित अंतराल $[0, T]$ पर होता है। यह एकमात्र भाग है जिसे HOC विधि का उपयोग करके संख्यात्मक रूप से हल किया जाता है।
     • अंतिम दृष्टिकोण (Final Approach): प्रणाली $[T, \infty)$ के दौरान सैडल (saddle) की ओर बढ़ती है। इस चरण को भी सैडल के आसपास गतिकी को रैखिकीकृत करके विश्लेषणात्मक रूप से माना जाता है। अस्थिर दिशा (unstable direction) को स्पष्ट रूप से संभाला जाता है, जबकि स्थिर दिशाएँ स्वाभाविक रूप से सैडल की ओर प्रवाहित होती हैं।
   • समय-अनुवाद विसंगति को तोड़ना (Breaking Time-Translation Degeneracy): यह सुनिश्चित करने के लिए कि इष्टतम प्रक्षेप पथ समय अक्ष के साथ स्वतंत्र रूप से फिसल न सकें, एक टर्मिनल बाधा (terminal constraint) लागू की जाती है। प्रक्षेप पथ को समय $T$ पर सैडल के चारों ओर एक पूर्व-निर्धारित अनुप्रस्थ दीर्घवृत्त (transverse ellipsoid) को पार करना आवश्यक है। यह बाधा समरूपता (symmetry) को तोड़ती है, जिससे एक सुव्यवस्थित संख्यात्मक समस्या सुनिश्चित होती है।

3. संख्यात्मक कार्यान्वयन (Numerical Implementation):

   • समग्र अनुकूलन समस्या (optimization problem) प्रक्षेप पथ और इष्टतम एंडपॉइंट्स ($\phi_0$ सिंक के पास और $\phi_T$ दीर्घवृत्त पर) के सापेक्ष एक्शन के योग को न्यूनतम करती है।
   • बाधाओं को हल करने के लिए ऑगमेंटेड लैग्रेंजियन विधि (Augmented Lagrangian method) का उपयोग किया जाता है। यह दृष्टिकोण ऑब्जेक्टिव फंक्शन को कॉनवेक्सिफाई (convexify) करता है और लैग्रेंज मल्टीप्लायर्स तथा पेनल्टी पैरामीटर्स को टर्मिनल दीर्घवृत्त बाधा को संतुष्ट करने के लिए पुनरावृत्ति (iteratively) से अपडेट करता है, जिससे अनंत पेनल्टी पैरामीटर्स की आवश्यकता के बिना इल-कंडीशनिंग से बचा जा सके।
   • एल्गोरिदम कुशल अभिसरण (convergence) के लिए नॉन-लीनियर कंजुगेट ग्रेडिएंट विधियों का उपयोग करता है।

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

• डीजेनरेट शोर का विस्तार: यह शोध पत्र डीजेनरेट शोर वाले सिस्टम के लिए संक्रमण दरों और सर्वाधिक संभावित पथों (most probable paths) की गणना करने के लिए एक सुदृढ़ कम्प्यूटेशनल ढांचा प्रदान करता है, जो एक ऐसी स्थिति है जहाँ पहले सहप्रसरण मैट्रिक्स व्युत्क्रमण वाली विधियाँ अक्षम थीं।
• रंगीन शोर को संभालना: स्टेट स्पेस में रैखिक फिल्टरों को एकीकृत करके और HOC का उपयोग करके, यह विधि श्वेत शोर अनुमानों पर निर्भर हुए बिना भौतिक रूप से यथार्थवादी रंगीन शोर को प्रभावी ढंग से संभालती है।
• अनंत समय-क्षितिज समाधान: विभाजन की विधि (Method of Division) अनंत समय के लिए दुर्लभ घटनाओं के लिए एक नवीन कम्प्यूटेशनल रणनीति प्रदान करती है। यह कम्प्यूटेशनल रूप से कठिन गैर-रैखिक संक्रमण क्षेत्र को अलग करती है और साम्यावस्था के आसपास के स्पर्शोन्मुखी व्यवहार (asymptotic behavior) को विश्लेषणात्मक रूप से हल करती है, जिससे कम्प्यूटेशनल दक्षता और अभिसरण में उल्लेखनीय वृद्धि होती है।
• समय-अनुवाद विसंगति को हटाना: एक अनुप्रस्थ टर्मिनल बाधा (transverse terminal constraint) का परिचय अनंत-समय प्रक्षेप पथों के लिए संख्यात्मक समाधानों को स्थिर करने के लिए एक कठोर तंत्र प्रदान करता है।

परिणाम और अनुप्रयोग
प्रस्तावित ढांचे की प्रभावशीलता को दो विशिष्ट केस स्टडीज के माध्यम से प्रदर्शित किया गया है:

1. इनवर्टेड डबल-वेल पोटेंशियल (Inverted Double-Well Potential): एक सरलीकृत मॉडल जिसका उपयोग डीजेनरेट और रंगीन शोर को संभालने की पद्धति की सैद्धांतिक क्षमताओं को प्रदर्शित करने के लिए किया गया है।
2. शिप कैपसाइज मॉडल (Ship Capsize Model): एक सरलीकृत रोल-हीव मॉडल जो समुद्री इंजीनियरिंग में एक महत्वपूर्ण अनुप्रयोग का प्रतिनिधित्व करता है। पद्धति सफलतापूर्वक स्टोकैस्टिक वेव फोर्सिंग के तहत एक जहाज के लिए सबसे संभावित पलायन पथों और संक्रमण दरों के घातांकीय स्केलिंग नियमों की गणना करती है, जो कैपसाइज दर में डोमिनेंट एक्सपोनेंशियल योगदान को पकड़ती है।

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