Violation of local equilibrium thermodynamics in one-dimensional Hamiltonian-Potts model

तकनीकी सारांश: एक आयामी हैमिल्टोनियन-पोट्स मॉडल में स्थानीय संतुलन ऊष्मागतिकी का उल्लंघन

समस्या विवरण
हैमिल्टोनियन प्रणालियों की गतिशीलता को आमतौर पर स्थानीय संतुलन ऊष्मागतिकी (local equilibrium thermodynamics) मानते हुए हाइड्रोडायनामिक समीकरणों द्वारा वर्णित किया जाता है। हालाँकि, यह सन्निकटन तीव्र ऊष्मा धाराओं (heat currents) और तीव्र तापमान प्रवणता (temperature gradients) की उपस्थिति में, विशेष रूप से प्रथम-क्रम के प्रावस्था संक्रमण (first-order phase transitions) के दौरान विफल हो जाता है। एक प्रमुख प्रश्न यह है कि क्या एक गैर-संतुलन स्थिर अवस्था (nonequilibrium steady state - NESS) में सह-अस्तित्व वाले क्रमित (ordered) और अव्यवdt (disordered) चरणों के बीच का तापमान, साम्यावस्था संक्रमण तापमान ($T_c$) के बराबर होता है। जबकि वैश्विक ऊष्मागतिकी (global thermodynamics) ऊष्मा धाराओं द्वारा मेटास्टेबल अवस्थाओं के स्थिरीकरण के कारण $T_c$ से विचलन की भविष्यवाणी करती है, दो-आयामी हैमिल्टोनियन-पोट्स मॉडलों का पिछले संख्यात्मक सत्यापन महत्वपूर्ण परिमित-आकार (finite-size) और परिमित-खुरदरापन (finite-roughness) प्रभावों से ग्रस्त रहा है। ये प्रभाव, जो मोटे स्थानिक विविक्तिकरण (coarse spatial discretization) और सीमित प्रणाली आकार से उत्पन्न होते हैं, अक्सर स्थानीय संतुलन के उल्लंघन को छिपा देते हैं, जिससे थर्मोडायनामिक सीमा तक पहुँचने के लिए गणनात्मक रूप से अत्यधिक कठिन सिमुलेशन की आवश्यकता होती है।

कार्यप्रणाली
इन सीमाओं को दूर करने के लिए, लेखक एक न्यूनतम, नियंत्रित संख्यात्मक मॉडल प्रस्तावित करते हैं: एक आयामी हैमिल्टोनियन-पोट्स मॉडल जिसमें भिन्नात्मक-क्रम (fractional-order) स्थानिक डेरिवेटिव शामिल हैं।

• मॉडल निर्माण: लघु-सीमा अंतःक्रियाओं (short-range interactions) वाला मानक एक-आयामी मॉडल प्रावस्था संक्रमण का समर्थन नहीं करता है। $n=5$ अवस्थाओं के लिए प्रथम-क्रम का प्रावस्था संक्रमण प्रेरित करने हेतु, लेखक एक गैर-स्थानीय स्थानिक डेरिवेटिव ऑपरेटर $D$ पेश करते हैं। इस ऑपरेटर को फूरियर प्रतिनिधित्व के माध्यम से परिभाषित किया गया है जहाँ निम्न तरंग संख्या (low wavenumbers) पर स्पेक्ट्रल घनत्व $\rho(k) = \text{const}$ के रूप में स्केल करता है, जो प्रभावी रूप से एक मानक दो-आयामी मॉडल के निम्न-तरंग संख्या घनत्व की नकल करता है। यह आयामी न्यूनीकरण (dimensional reduction) एक ऐसे एक-आयामी परिवेश में प्रावस्था सह-अस्तित्व के अध्ययन की अनुमति देता है जो अनुप्रस्थ ज्यामितीय बाधाओं से मुक्त है।
• हैमिल्टोनियन: प्रणाली एक हैमिल्टनियन $H = \int dx \, h$ द्वारा शासित है, जहाँ गतिज (kinetic) और ग्रेडिएंट पद $D$ के भिन्नात्मक डेरिवेटिव का उपयोग करते हैं। पोटेंशियल पद को एक नियमित सिम्प्लेक्स (simplex) के शीर्षों के रूप में व्यवस्थित $n$ समान मिनिमा (minima) बनाने के लिए निर्मित किया गया है।
• सिमुलेशन सेटअप: सीमाओं पर तापमान $T_1$ और $T_2$ (जो $T_1 < T_c < T_2$ को संतुष्ट करते हैं) वाले हीट बाथ लगाकर गैर-संतुलन स्थिर अवस्थाएँ (NESS) प्राप्त की जाती हैं। सीमा क्षेत्रों को लैंगव डायनामिक्स (Langevin dynamics) का उपयोग करके मॉडल किया गया है, जबकि केंद्रीय क्षेत्र ऊर्जा-संरक्षण हैमिल्टोनियन समीकरणों के माध्यम से विकसित होता है।
• संख्यात्मक कार्यान्वयन: हैमिल्टोनियन डायनामिक्स के लिए एक द्वितीय-क्रम सिम्प्लेक्टिक स्कीम (second-order symplectic scheme) और लैंगव डायनामिक्स के लिए एक द्वितीय-क्रम स्टोकेस्टिक रनगे-कुट्टा (stochastic Runge-Kutta) स्कीम का उपयोग करके गति के समीकरणों को एकीकृत किया जाता है। प्रणाली का आकार ($L_x = 256$) और जाली अंतराल (lattice spacing $\Delta x = 1/4$) विविक्तिकरण प्रभावों को न्यूनतम करने और गणनात्मक व्यवहार्यता बनाए रखने के लिए चुना गया है।

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

1. 1D में प्रावस्था सह-अस्तित्व की प्राप्ति: यह अध्ययन सफलतापूर्वक प्रदर्शित करता है कि भिन्नात्मक डेरिवेटिव का परिचय एक आयामी प्रणाली को प्रथम-क्रम का प्रावस्था संक्रमण और स्थिर प्रावस्था सह-अस्तित्व प्रदर्शित करने में सक्षम बनाता है, जो कि आमतौर पर उच्च आयामों या दीर्घ-परासी अंतःक्रियाओं (long-range interactions) तक सीमित एक घटना है।
2. स्थानीय संतुलन का उल्लंघन: संख्यात्मक सिमुलेशन प्रकट करते हैं कि स्थिर इंटरफेस (stationary interface) पर तापमान ($\theta$) साम्यावस्था संक्रमण तापमान ($T_c \approx 0.634$) से व्यवस्थित रूप से विचलित होता है। विशेष रूप से, इंटरफेस तापमान $T_c$ से अधिक पाया गया है, जो एक 'सुपरहीटेड' (superheated) क्रमित अवस्था के स्थिरीकरण का संकेत देता है। यह विचलन सांख्यिकीय अनिश्चितता से कहीं अधिक बना रहता है, जो स्थानीय संतुलन के टूटने का स्पष्ट प्रमाण प्रदान करता है।
3. वैश्विक ऊष्मागतिकी के साथ मात्रात्मक समझौता: देखे गए इंटरफेस तापमान वैश्विक ऊष्मागतिकी के भविष्यवाणियों के साथ मात्रात्मक रूप से सुसंगत हैं। वैश्विक तापमान $\tilde{T}$, जिसे स्थानीय तापमान के घनत्व-भारित स्थानिक औसत के रूप में परिभाषित किया गया है, और इंटरफेस तापमान विचलन के लिए सैद्धांतिक सूत्र (एक वेरिएशनल सिद्धांत से व्युत्पन्न) सिमुलेशन डेटा का सटीक वर्णन करते हैं।
4. तापीय चालकता विषमता (Thermal Conductivity Asymmetry): परिणाम दर्शाते हैं कि क्रमित चरण में तापीय चालकता ($\kappa_o$), अव्यक्त (disordered) चरण की तुलना में ($\kappa_d$) कम है, जो इस सैद्धांतिक भविष्यवाणी के अनुरूप है कि इंटरफेस तापमान उस चरण की ओर झुकता है जिसकी तापीय चालकता कम होती है।

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