Can the Brownian diffusion coefficient be reconstructed from Lyapunov exponents?

तकनीकी सारांश: लियापुनोव एक्सपोनेंट्स (Lyapunov Exponents) से ब्राउनियन डिफ्यूजन का पुनर्निर्माण

समस्या विवरण
यह शोध पत्र गैर-संतुलन सांख्यिकीय यांत्रिकी (non-equilibrium statistical mechanics) में दो भिन्न परिमाणात्मकों (quantifiers) के बीच संबंध की जांच करता है: प्रसार गुणांक ($D$), जो स्टोकेस्टिक प्रणालियों में मैक्रोस्कोपिक परिवहन को अभिलक्षणित करता है, और अधिकतम लियापुनोव एक्सपोनेंट ($\lambda_1$), जो नियतात्मक (deterministic) गतिशील प्रणालियों के अराजक गुणों को अभिलक्षणित करता है। सामान्यतः, इन दोनों मात्राओं के बीच कोई सीधा संबंध नहीं होता है; अराजक नियतात्मक प्रणालियाँ अक्सर प्रसार प्रदर्शित करती हैं, जबकि गैर-अराजक प्रणालियाँ आमतौर पर ऐसा नहीं करती हैं, और स्टोकेस्टिक प्रणालियों में उनके नियतात्मक समकक्षों के विपरीत अनंत कोलमोगोरोव-सिनाई एंट्रॉपी होती है। लेखक एक विशिष्ट गैर-संतुलन प्रणाली—एक स्थानिक रूप से आवर्ती, सममित विभव (potential) में ac-ड्रिवन कण—पर ध्यान केंद्रित करते हैं, जहाँ हालिया अध्ययनों से पता चला है कि प्रसार गुणांक ड्राइविंग आयाम (driving amplitude) के एक अर्ध-आवधिक (quasiperiodic) फलन के रूप में व्यवहार करता है। केंद्रीय प्रश्न यह है कि क्या इस जटिल, तापमान-निर्भर प्रसार व्यवहार को संगत शून्य-तापमान नियतात्मक प्रणाली के गुणों से पुनर्गठित किया जा सकता है।

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

1. स्टोकेस्टिक मॉडल: प्रणाली को एक आयामहीन लैंग्विन समीकरण (Langevin equation - Eq. 1) द्वारा वर्णित किया गया है जो घर्षण गुणांक $\gamma$ वाले एक ब्राउनियन कण का प्रतिनिधित्व करता है, जो विभव $U(x) = -\cos x$ में $a \sin(\omega t + \phi_0)$ द्वारा संचालित है, और तीव्रता $Q \propto T$ वाले थर्मल नॉइज़ से प्रभावित है। प्रसार गुणांक $D$ की गणना माध्य वर्ग विस्थापन (mean square displacement) की दीर्घकालिक सीमा के माध्यम से की जाती है।
2. नियतात्मक समकक्ष: शून्य-तापमान सीमा ($Q=0$) एक नियतात्मक समीकरण (Eq. 2) प्रदान करती है। लेखक इस प्रणाली के अधिकतम लियापुनोव एक्सपोनेंट $\lambda_1$ का विश्लेषण करते हैं, जो विचाराधीन स्ट्रॉन्ग-फ्रिक्शन शासन में गैर-अराजक (ऋणात्मक $\lambda_1$) है।
3. विश्लेषणात्मक सन्निकटन (Analytical Approximation): इस अंतर को पाटने के लिए, लेखक एक अनुमानित सिद्धांत विकसित करने हेतु "वाइब्रेशनल मैकेनिक्स" का उपयोग करते हैं। वे गति को धीमी और तेज़ घटकों में विभाजित करते हैं, जिससे एक प्रभावी डैम्प्ड हार्मोनिक ऑसिलेटर समीकरण (Eq. 19) प्राप्त होता है जिसका पुनर्सामान्यीकृत विभव कठोरता (renormalized potential stiffness) $k = J_0(\psi)$ है। इस रैखिकीकृत प्रणाली की रिलैक्सेशन दरों को लियापुनोव एक्सपोनेंट्स के बराबर माना जाता है।
4. पैरामीटर स्पेस: विश्लेषण स्ट्रॉन्ग-फ्रिक्शन सेक्टर ($\gamma = 3$) पर केंद्रित है जहाँ नियत प्रणाली केवल "लॉक्ड" (आवर्ती) प्रक्षेपवक्र (trajectories) प्रदर्शित करती है, जिसमें आवृत्ति $\omega = 1.59$ और तापमान $Q = 0.1$ को स्थिर रखते हुए ड्राइविंग आयाम $a$ को बदला जाता है। मजबूती (robustness) का परीक्षण $\gamma$ और $Q$ के परिवर्तनों के विरुद्ध किया जाता है।

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

• मजबूत सहसंबंध: प्राथमिक निष्कर्ष शोर युक्त प्रणाली में प्रसार गुणांक $D(a)$ और नियतात्मक प्रणाली में अधिकतम लियापुनोव एक्सपोनेंट $\lambda_1(a)$ के बीच एक आश्चर्यजनक गुणात्मक और मात्रात्मक सहसंबंध है। बावजूद इसके कि नियत प्रणाली गैर-अराजक और गैर-विसरणात्मक (non-diffusive) है, इसका $\lambda_1$ प्रसार गुणांक $D$ के दोलनी व्यवहार की नकल करता है।
• दोलन का तंत्र:
  • प्रसार ($D$): $D$ का अर्ध-आवधिक व्यवहार नियतात्मक अट्रैक्टर संरचना में होने वाले द्विभाजन (bifurcations) द्वारा संचालित होता है। जैसे-जैसे $a$ बढ़ता है, प्रणाली "सामान्य अवस्थाओं" (विभव के न्यूनतम बिंदुओं के आसपास दोलन) और "सममिति-भंग अवस्थाओं" (विभव के अधिकतम बिंदुओं के आसपास दोलन) के बीच संक्रमण करती है। प्रसार गुणांक तब चरम पर होता है जब कण का प्रक्षेपवक्र विभव के दोनों न्यूनतम और अधिकतम बिंदुओं को समाहित कर लेता है, जिससे कुओं (wells) के बीच थर्मल जंप सुलभ हो जाते हैं।
  • लियापुनोव एक्सपोनेंट ($\lambda_1$): $\lambda_1$ में परिवर्तन आवर्ती अट्रैक्टर्स की ओर प्रक्षेपवक्रों के रिलैक्सेशन समय से जुड़े हैं। द्विभाजन बिंदुओं (क्रिटिकल एम्पलीट्यूड $a_c, a_1, a_2, \dots$) के पास, प्रणाली "क्रिटिकल स्लोइंग डाउन" प्रदर्शित करती है, जहाँ $\lambda_1 \to 0$ (ऋणात्मक दिशा से शून्य की ओर बढ़ना) होता है, जो रिलैक्सेशन समय में विचलन के अनुरूप है।
• पुनर्गठन सूत्र (Reconstruction Formula): लेखक प्रसार गुणांक को अधिकतम लियापुनोव एक्सपोनेंट से पुनर्गठित करने के लिए एक अनुमानित सूत्र (Eq. 24) प्रस्तावित करते हैं:
  $$D = \frac{Q}{\gamma I_0^2(-\lambda_1(\lambda_1 + \gamma)/Q)}$$
  जहाँ $I_0$ प्रथम प्रकार का मॉडिफाइड बेसेल फंक्शन (modified Bessel function of the first kind) है। यह सूत्र नियतात्मक प्रणाली से प्राप्त $\lambda_1$ के सटीक संख्यात्मक मान का उपयोग करके स्टोकेस्टिक प्रणाली में $D$ की भविष्यवाणी करता है।
• सहमति और विसंगतियां: प्रस्तावित सूत्र अधिकांश पैरामीटर रेंज में लैंग्विन समीकरण के संख्यात्मक सिमुलेशन के साथ "अत्यंत संतोषजनक" सहमति दिखाता है। हालांकि, प्रसार गुणांक के स्थानीय मैक्सिमा के निकट (जहाँ द्विभाजन होते हैं) विसंगतियां देखी गई हैं, जहाँ रिलैक्सेशन समय का रैखिक सन्निकटन विफल हो जाता है।
• मजबूती (Robustness): $D$ और $\lambda_1$ के बीच सहसंबंध घर्षण गुणांक $\gamma$ और तापमान $Q$ में छोटे परिवर्तनों के विरुद्ध मजबूत है। हालांकि $\gamma$ के साथ $\lambda_1(a)$ का विशिष्ट आकार बदलता है, लेकिन व्युत्क्रम संबंध (जहाँ $\lambda_1$ का शून्य के करीब पहुँचना बढ़ते $D$ के अनुरूप है) बना रहता है। यह सहसंबंध तब भी बना रहता है जब तापमान को पांच गुना बढ़ा दिया जाता है, हालांकि $D > D_0$ (आइंस्टीन डिफ्यूजन) का क्षेत्र चौड़ा हो जाता है।

महत्व और दावे
यह शोध पत्र दावा करता है कि विशिष्ट अनुकूलित पैरामीटर व्यवस्थाओं में, एक स्टोकेस्टिक प्रणाली के जटिल परिवहन गुणों (प्रसार) को उसकी नियतात्मक, शून्य-तापमान समकक्ष (लियापुनोव एक्सपोनेंट्स) के गतिशील गुणों से सटीक रूप से पुनर्गठित किया जा सकता है। यह महत्वपूर्ण है क्योंकि यह एक ऐसे शासन में सूक्ष्म नियतात्मक गतिकी और मैक्रोस्कोपिक परिवहन के बीच एक लिंक प्रकट करता है जहाँ नियत प्रणाली गैर-अराजक और गैर-विसरणात्मक है।

लेखक इस बात पर जोर देते हैं कि यह सभी प्रणालियों के लिए एक सामान्य नियम नहीं है, बल्कि स्ट्रॉन्ग-फ्रिक्शन सीमा में प्रणाली की नियमित अट्रैक्टर संरचना से उत्पन्न होने वाली एक विशिष्ट घटना है। वे अनुमानित सूत्र को पूर्ण स्टोकेस्टिक सिमुलेशन के बिना प्रसार गुणांक का अनुमान लगाने के लिए एक व्यावहारिक उपकरण के रूप में प्रस्तावित करते हैं, बशर्ते कि नियत लियापुनोव एक्सपोनेंट ज्ञात हो। यह कार्य इस बात को रेखांकित करता है कि शोर की उपस्थिति में "धुंधले" (blurred) नियत ऑर्बिट्स अभी भी प्रसार तंत्र को निर्देशित करते हैं, और इन ऑर्बिट्स की स्थिरता (जिसे $\lambda_1$ द्वारा मापा जाता है) थर्मल एस्केप और उसके बाद के प्रसार की सुगमता के लिए एक प्रॉक्सी (proxy) के रूप में कार्य करती है।