First passage time for an underdamped harmonic oscillator and application to the power of an information engine

तकनीकी सारांश: अनडैम्पड हार्मोनिक ऑसिलेटर के लिए फर्स्ट पैसेज टाइम और सूचना इंजन में इसका अनुप्रयोग

समस्या विवरण
फर्स्ट पैसेज टाइम (FPT) वितरण का अनुमान लगाना सांख्यिकीय भौतिकी की एक मौलिक समस्या है, जिसके अनुप्रयोग रासायनिक अभिक्रिया दरों से लेकर सूचना इंजनों (information engines) तक विस्तृत हैं। जबकि लैंज़िवन समीकरणों द्वारा वर्णित ओवरडैम्पड प्रणालियों के लिए FPT वितरण अच्छी तरह से स्थापित हैं, अनडैम्पड शासन (underdamped regime) चुनौतीपूर्ण बना हुआ है। अनडैम्पड प्रणालियों में, जड़त्व (inertia) विस्थापन के समीकरणों में एक द्वितीय अवकलज (second derivative) पेश करता है, जो प्रक्रिया को नॉन-मार्कोवियन (non-Markovian) बना देता है। फलस्वरूप, FPT संभाव्यता की गणना करना ऐतिहासिक रूप से एक "कठिन कार्य" रहा है, क्योंकि बंद-रूप समाधान (closed-form solutions) सामान्यतः केवल विशिष्ट सीमाओं (जैसे कि दीर्घ-समय सीमा में मुक्त कण) के लिए ही उपलब्ध होते हैं। यह शोध पत्र एक ब्राउनियन हार्मोनिक ऑसिलेटर के लिए FPT वितरण, $P(t_{fp})$ की गणना को संबोधित करता है, जो एक परिमित गुणवत्ता कारक (quality factor) $Q$ के साथ एक स्थिति थ्रेशोल्ड $x_B$ को पार करता है।

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

1. PDF का अपघटन (Decomposition): कुल संभाव्यता घनत्व फलन $P(t_{fp})$ को तीन योगदानों में विभाजित किया गया है:

   • $P_I \delta(t_{fp})$: तात्कालिक ट्रिगरिंग (instantaneous triggering) को दर्शाता है जहाँ प्रारंभिक स्थिति $x(0) \geq x_B$ है।
   • $P_{II}(t_{fp})$: पहले दोलन काल के भीतर की गतिशीलता ($0 < t_{fp} \leq 2\pi$) का वर्णन करता है।
   • $P_{III}(t_{fp})$: दीर्घ-समय व्यवहार ($t_{fp} > 2\pi$) को कैप्चर करता है।

2. अल्प-समय शासन ($t_{fp} \leq 2\pi$): कम डैम्पिंग ($Q \gg 1$) के लिए, लेखक एक हैमिल्टोनियन सन्निकटन (Hamiltonian approximation) का उपयोग करते हैं। वे मानते हैं कि एक अवधि के दौरान कुल ऊर्जा $E$ लगभग स्थिर रहती है।

   • तात्कालिक ट्रिगर ($P_I$): प्रारंभिक स्थिति के बोल्ट्ज़मैन वितरण का उपयोग करके गणना की गई है, जिससे $P_I = \frac{1}{2}\text{erfc}(\sqrt{B})$ प्राप्त होता है, जहाँ $B = \frac{1}{2}x_B^2$ ऊर्जा अवरोध (energy barrier) है।
   • प्रथम अवधि ($P_{II}$): गतिशीलता को नियतकालिक साइनसोइडल गति $x(t) = x_E \cos(\theta_0 + t)$ के रूप में माना जाता है। संभाव्यता को प्रारंभिक ऊर्जा $E$ और चरण $\theta_0$ के संयुक्त वितरण पर एकीकरण करके प्राप्त किया जाता है। यह एक प्लेटो (plateau) और उसके बाद एक तीव्र क्षय को दर्शाता है, जो समय $t$ पर सुलभ अधिकतम ऊर्जा $E^\dagger(t)$ द्वारा निर्धारित होता है।

3. दीर्घ-समय शासन ($t_{fp} > 2\pi$): एक दोलन काल से अधिक समय के लिए, प्रणाली का व्यवहार थर्मल नॉइज़ के कारण ऊर्जा प्रसार (energy diffusion) द्वारा नियंत्रित होता है।

   • उच्च $Q$ सीमा: समस्या को ऊर्जा अवरोध के पार क्रैमरर एस्केप (Kramers escape) समस्या के रूप में मैप किया गया है। पलायन दर $\Gamma(t)$ की गणना समय-निर्भर ऊर्जा थ्रेशोल्ड $E^\dagger(t)$ के आधार पर की जाती है, जो अंततः एक स्थिर दर $\Gamma_B$ की ओर अभिसरित होती है।
   • परिमित $Q$: लेखक $x_B$ पर एक अवशोषक सीमा स्थिति (absorbing boundary condition) के अधीन, क्रैमरर डिफरेंशियल ऑपरेटर से जुड़े सबसे धीमे आइजनवैल्यू (eigenvalue) $-\lambda_1(x_B)$ की गणना करके पलायन दर $\Gamma$ निर्धारित करते हैं।
   • परिणामी टेल (tail) $P_{III}$, घातांकीय क्षय $e^{-\Gamma_B t_{fp}}$ का अनुसरण करती है।

मुख्य परिणाम
शोध पत्र एक बंद-रूप विश्लेषणात्मक अभिव्यक्ति (समीकरण 2) प्रस्तुत करता है जो केवल ऊर्जा अवरोध $B$ और गुणवत्ता कारक $Q$ द्वारा पैरामीटराइज्ड है।

• प्रायोगिक सत्यापन: सैद्धांतिक भविष्यवाणियों का परीक्षण एक अनडैम्पड माइक्रो-कैंटिलीवर ($Q \approx 7$) का उपयोग करके किया गया जो कमरे के तापमान पर कार्य कर रहा था। स्थिति को 2 MHz सैंपलिंग दर पर पिकोमीटर रिज़ॉल्यूशन के साथ मापा गया था।
• सहमति: प्रायोगिक डेटा ने तीन अलग-अलग शासन के तहत सैद्धांतिक मॉडल के साथ "उत्कृष्ट सहमति" दिखाई:
  1. तात्कालिक ट्रिगर्स की संभाव्यता ($P_I$)।
  2. $0 < t < \pi$ के लिए वितरण का प्लेटो मान ($P_{II}$)।
  3. दीर्घ-समय टेल के लिए घातांकीय क्षय दर ($P_{III}$)।
• सूचना इंजन अनुप्रयोग: व्युत्पन्न FPT वितरण का उपयोग एक सूचना इंजन की शक्ति (power) का अनुमान लगाने के लिए किया गया है। इस सेटअप में, जब ऑसिलेटर एक थ्रेशोल्ड $x_B$ को पार करता है, तो क्षमता कुण्ड (potential well) खिसक जाती है और कार्य निकाला जाता है। औसत कार्य के योग और ड्वेल टाइम (dwell time) के अनुपात के रूप में औसत शक्ति $\langle P \rangle$ की गणना की गई थी। सैद्धांतिक पावर कर्व प्रायोगिक डेटा से मेल खाता है, जो एक इष्टतम ऑपरेटिंग पॉइंट की पहचान करता है जहाँ थ्रेशोल्ड $x_B \approx 2$ और शिफ्ट $L \approx 1$ है।

महत्व और दावे
लेखक दावा करते हैं कि यह कार्य अनडैम्पड प्रणालियों के लिए FPT समस्या में "एक नई अंतर्दृष्टि" प्रदान करता है, जो लघु-समय हैमिल्टोनियन गतिशीलता और दीर्घ-समय विसरणात्मक पलायन (diffusive escape) के बीच के अंतर को पाटते हुए, गुणवत्ता कारकों की पूरी सीमा को कवर करने वाला समाधान प्रदान करता है।

• सैद्धांतिक प्रगति: शोध पत्र का दावा है कि व्युत्पन्न बंद-रूप समाधान (समीकरण 2), अनडैम्पड प्रणालियों के लिए FPT सांख्यिकी में एक प्रगति का प्रतिनिधित्व करता है। यह उल्लेख करता है कि पूर्व में, अनडैम्पड प्रणाली के लिए ज्ञात एकमात्र बंद-रूप समाधान दीर्घ-समय सीमा में मुक्त कण (रैंडम एक्सेलरेशन प्रोसेस) के लिए था।
• क्रैमरर थ्योरी से भिन्नता: लेखक इस बात पर जोर देते हैं कि उनका विवरण मानक क्रैमरर दर सिद्धांत से भिन्न है, जो आमतौर पर बड़े अवरोधों के विसरणात्मक क्रॉसिंग (diffusive crossing) पर लागू होता है। उनके वितरण में, घातांकीय टेल केवल एक दोलन काल से काफी अधिक समय के लिए उभरती है, जबकि अल्प-समय व्यवहार जड़त्वीय दोलनों (inertial oscillations) द्वारा हावी होता है।
• व्यावहारिक उपयोगिता: यह विधि जड़त्व वाले सूचना इंजनों में शक्ति के सटीक अनुमान को सक्षम बनाती है, जो एक ऐसा पैरामीटर है जिसे पूर्ण वितरण ज्ञान के बिना मापना कठिन था। वे कहते हैं कि यह अनडैम्पड गतिशीलता वाले अन्य क्षेत्रों, जैसे स्टोकेस्टिक थर्मोडायनामिक्स में FPT सांख्यिकी को लागू करने का मार्ग प्रशस्त करता है।

शोध पत्र इस बात के साथ समाप्त होता है कि गुणवत्ता कारक के प्रभाव के संबंध में गणनाओं और संख्यात्मक सिमुलेशन का पूर्ण विवरण एक संलग्न शोध पत्र [22] में प्रदान किया गया है।