Unified Statistical Theory of Heat Conduction in Nonuniform Media

यहाँ यी ज़ेंग और जियानजुन डोंग द्वारा लिखित शोध पत्र "Unified Statistical Theory of Heat Conduction in Nonuniform Media" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

शास्त्रीय ऊष्मा चालन (Classical heat conduction) पारंपरिक रूप से फूरियर के नियम (Fourier's law) द्वारा वर्णित किया जाता है, जो यह मानता है कि ऊष्मा प्रवाह (heat flux) तापमान प्रवणता (temperature gradients) के प्रति तत्काल (instantaneously) और स्थानीय (locally) प्रतिक्रिया देता है। हालांकि मैक्रोस्कोपिक स्तर पर यह सफल है, लेकिन यह ढांचा अत्यंत तीव्र समय पैमाने (सब-टेराहर्ट्ज़) और नैनोमीटर लंबाई के पैमाने (सब-10 nm) पर विफल हो जाता है। इन स्तरों पर, अंतर्निहित फोनोन रिलैक्सेशन समय (phonon relaxation times) और माध्य मुक्त पथ (mean free paths - MFP) प्रयोगात्मक पैमानों के तुल्य हो जाते हैं, जिससे गैर-फूरियर घटनाएं उत्पन्न होती हैं, जैसे:

• क्वासी-बैलिस्टिक परिवहन (Quasi-ballistic transport): जहाँ फोनोन बिना बिखराव (scattering) के यात्रा करते हैं।
• हाइड्रोडायनामिक परिवहन (द्वितीय ध्वनि/Second Sound): तरंग-समान तापमान प्रसार।
• पोइज़ुइल प्रवाह (Poiseuille flow): श्यानता-समान (viscous-like) फोनोन प्रवाह।

मौजूदा सामान्यीकृत मॉडल (जैसे, टेम्पोरल मेमोरी के लिए मैक्सवेल-कैट्टासियो-वर्नोट और स्थानिक गैर-स्थानीयता के लिए गुयेर-क्रुमहांसल) आमतौर पर विशिष्ट स्पर्शोन्मुखी (asymptotic) क्षेत्रों तक सीमित होते हैं। उनमें एक एकीकृत, व्यवस्थित रूप से नियंत्रित ढांचे का अभाव है जो टेम्पोरल मेमोरी (temporal memory), स्थानिक गैर-स्थानीयता (spatial nonlocality), और सामग्री की विषमता (material heterogeneity) (जैसे इंटरफेस, विस्कृत ठोस) को बिना किसी अनुभवजन्य क्लोजर पैरामीटर या क्षेत्र-विशिष्ट मॉडल चयन के, एक साथ संभाल सके।

2. कार्यप्रणाली (Methodology)

लेखक गैर-संतुलन सांख्यिकीय यांत्रिकी (nonequilibrium statistical mechanics) से ज़्वानजिग प्रोजेक्शन-ऑपरेटर फॉर्मलिज्म (Zwanzig projection-operator formalism) पर आधारित एक एकीकृत सैद्धांतिक ढांचा विकसित करते हैं।

• सूक्ष्म अपघटन (Microscopic Decomposition): फोनोन बोल्ट्ज़मैन ट्रांसपोर्ट इक्वेशन (phBTE) को तीन वर्गों के प्रोजेक्टेड वेरिएबल्स में विभाजित किया गया है:
  1. संरक्षित चर (Conserved variables): स्थानीय तापमान क्षेत्र ($T$)।
  2. विषम-पैरिटी (flux) मोड (Odd-parity modes): करंट-कैरिंग मोड ($\eta_\lambda$) जो ऊष्मा प्रवाह का पुनर्निर्माण करते हैं।
  3. सम-पैरिटी (nonlocal) मोड (Even-parity modes): $\vec{q}$-सिमेट्रिक डिस्टॉर्शन मोड ($\theta_\iota$) जो स्थानिक गैर-स्थानीयता को मध्यस्थता प्रदान करते हैं।
• कर्नेल का व्युत्पन्न (Derivation of the Kernel): अप्रासंगिक (तेज) वेरिएबल्स को औपचारिक रूप से हटाकर, लेखक एक कारण संबंधी द्वि-बिंदु स्पैटियोटेम्पोरल कर्नेल (causal two-point spatiotemporal kernel) व्युत्पन्न करते हैं, जिसे $\overleftrightarrow{Z}(\vec{r}, \vec{R}, t)$ के रूप में दर्शाया गया है। यह कर्नेल सभी स्थान ($\vec{R}$) और समय ($t'$) में तापमान प्रवणता के इतिहास को स्थिति $\vec{r}$ और समय $t$ पर ऊष्मा प्रवाह से जोड़ता है।
• सूक्ष्म परिभाषा (Microscopic Definition): कर्नेल को कठोरता से एक स्पैटियोटेम्पोरल ग्रीन-कुबो (stGK) संबंध के माध्यम से परिभाषित किया गया है:
  $$\overleftrightarrow{Z}(\vec{r}, \vec{R}, t) = \frac{1}{k_B T_0^2} \langle \vec{j}(\vec{r}, t) \otimes \vec{j}(\vec{R}, 0) \rangle_{eq}$$
  यह कर्नेल को एक स्थान-विभेदित साम्यावस्था ऊष्मा-प्रवाह समय-सहसंबंध फलन (equilibrium heat-flux time-correlation function) के रूप में परिभाषित करता है, जो क्रिस्टलीय, विस्कृत और विषम प्रणालियों के लिए लागू है और फोनोन-गैस चित्र (phonon-gas picture) पर निर्भर नहीं है।
• **इंटरफेस मॉडलिंग (Interface Modeling):ness Framework को सामग्री इंटरफेस तक विस्तारित करने के लिए इंटरफेस को कर्नेल के भीतर एक स्थानिक रूप से स्थानीयकृत कपलिंग के रूप में माना जाता है, न कि एक सीमा शर्त (boundary condition) के रूप में। यह कपिट्ज़ा प्रतिरोध (Kapitza resistance) को स्थानिक रूप से विभेदित कर्नेल के एक कोर्स-ग्रेन्ड लिमिट के रूप में प्राप्त करने की अनुमति देता है।

3. मुख्य योगदान (Key Contributions)

1. एकीकृत संरचनात्मक वस्तु (Unified Constitutive Object): शोध पत्र स्पैटियोटेम्पोरल कर्नेल $\overleftrightarrow{Z}$ को ऊष्मा परिवहन के प्राथमिक संरचनात्मक डिस्क्रिप्टर के रूप में पहचानता है। यह टेम्पोरल मेमोरी और स्थानिक गैर-स्थानीयता को समान धरातल पर रखता है, जिससे अनुभवजन्य क्लोजर पैरामीटर या क्षेत्र-विशिष्ट मॉडल के बीच चयन करने की आवश्यकता समाप्त हो जाती है।
2. व्यवस्थित स्पर्शोन्मुखी न्यूनीकरण (Systematic Asymptotic Reductions): लेखक प्रदर्शित करते हैं कि शास्त्रीय और गैर-फूरियर मॉडल इस एकल कर्नेल के नियंत्रित सीमाएँ (limits) के रूप में उभरते हैं:
   • लोकल-ट्रांजिएंट (LT): टेम्पोरल मेमोरी बनाए रखता है लेकिन स्थानिक स्थानीयता मानता है (मैक्सवेल-कैट्टासियो-वर्नोट में बदल जाता है)।
   • ज़्वानजिग-डिफ्यूजन (ZD): स्थानिक गैर-स्थानीयता बनाए रखता है लेकिन अर्ध-स्थिर समय विकास मानता है (सामान्यीकृत गैर-स्थानीय प्रसार)।
   • फोनोन-हाइड्रोडायनामिक (PH): तब उभरता है जब कुछ सामूहिक मोड हावी होते हैं (गुयेर-क्रुमहांसल में बदल जाता है)।
   • एटेन्युएटेड-स्ट्रीमिंग (AS): परिमित प्रसार गति और घातांकीय क्षय के साथ क्वासी-बैलिस्टिक परिवहन का वर्णन करता है।
   • हाइब्रिड कलेक्टिव-स्ट्रीमिंग: उन режиनों को पकड़ता है जहाँ सामूहिक और स्ट्रीमिंग मोड सह-अस्तित्व में होते हैं (MFP वितरण वाले पदार्थों में सामान्य)।
3. इंटरफेस के लिए सूक्ष्म आधार (Microscopic Foundation for Interfaces): यह ढांचा इंटरफेशियल थर्मल रेजिस्टेंस (कपिट्ज़ा रेजिस्टेंस) का संरचनात्मक व्युत्पन्न प्रदान करता है। यह दिखाता है कि तापमान जंप एक विच्छेदन (discontinuity) नहीं है, बल्कि कर्नेल के स्थानिक सहसंबंधों के बाउंड्री-लेयर डिस्टॉर्शन का एक कम किया गया डिस्क्रिप्टर है।
4. सामान्य प्रयोज्यता (General Applicability): यह ढांचा केवल क्रिस्टलीय ठोसों तक सीमित नहीं है। विस्कृत हार्मोनिक सॉलिड्स के लिए, यह एलन-फेल्डमैन सीमा (Allen–Feldman limit) को पुनः प्राप्त करता है, जो फोनोन-क्वासिपार्टिक चित्र से परे इसकी प्रयोज्यता को दर्शाता है।

4. परिणाम और संख्यात्मक चित्रण (Results and Numerical Illustration)

लेखकों ने प्रथम-सिद्धांत फोनोन गुणों (RTA के भीतर) का उपयोग करके 300 K पर बल्क सिलिकॉन पर इस सिद्धांत को लागू किया और ट्रांजिएंट थर्मल ग्रेटिंग (TTG) प्रयोगों का अनुकरण किया।

• कर्नेल संरचना (Kernel Structure): सिलिकॉन के लिए RTA कर्नेल MFPs और रिलैक्सेशन समय के एक विस्तृत स्पेक्ट्रम को प्रकट करता है। कोई भी एकल लंबाई या समय पैमाना परिवहन को चरित्रित नहीं करता है; इसके बजाय, यह एक मल्टी-स्केल वितरण द्वारा शासित होता है।
• TTG क्षय गतिकी (TTG Decay Dynamics):
  • स्थानिक गैर-स्थानीयता का प्रभुत्व: अध्ययन में पाया गया कि स्थानिक गैर-स्थालिकता (जो फोनोन MFPs के वितरण से जुड़ी है) फूरियर परिवहन से विचलन का प्राथमिक स्रोत है।
  • टेम्पोरल मेमोरी का दमन: जबकि अकेले टेम्पोरल मेमोरी (LT सीमा में) मजबूत दोलनकारी "सेकंड-साउंड-जैसे" व्यवहार की भविष्यवाणी करता है, पूर्ण कर्नेल में स्थानिक गैर-स्थानीयता को शामिल करने से स्थानिक डीफेजिंग (spatial dephasing) होता है। विभिन्न MFPs वाले मोड अलग-अलग स्थानिक चरणों (phases) को संचित करते हैं, जिससे विनाशकारी हस्तक्षेप (destructive interference) होता है।
  • मात्रात्मक निष्कर्ष: $L=100$ nm के ग्रेटिंग पीरियड पर, LT मॉडल ने प्रारंभिक सिग्नल के $>50\%$ आयाम के साथ दोलन की भविष्यवाणी की। पूर्ण कर्नेल ने इसे $<5\%$ तक कम कर दिया, जो प्रयोगात्मक अवलोकनों के अनुरूप है जो सिलिकॉन में 300 K पर गैर-फूरियर क्षय तो दिखाते हैं लेकिन दोलन नहीं दिखाते।
• सत्यापन (Validation): पूर्ण कर्नेल और उसके ज़्वानजिग-डिफ्यूजन (ZD) न्यूनीकरण ने TTG क्षय दरों और फैलाव संबंधों ($\Gamma(q)$) को सटीक रूप से पुनरुत्पादित किया, जबकि एकल-पैरामीटर हाइड्रोडायनामिक (गुयेर-क्रुमहांसल) मॉडल विशिष्ट पैमानों पर सटीकता खोए बिना सभी लंबाई पैमानों पर व्यवहार को पकड़ने में विफल रहा।

5. महत्व (Significance)

• सैद्धांतिक एकीकरण (Theoretical Unification): यह कार्य सूक्ष्म सांख्यिकीय यांत्रिकी और मैक्रोस्कोपिक निरंतरता परिवहन (continuum transport) के बीच के अंतर को पाटता है। यह परमाणु सिमुलेशन (ग्रीन-कुबो के माध्यम से) और निरंतरता समीकरणों के बीच एक कठोर, पैरामीटर-मुक्त लिंक प्रदान करता है।
• "सेकंड साउंड" विवाद का समाधान: इसके निष्कर्ष स्पष्ट करते हैं कि तरंग-समान तापीय प्रतिक्रियाएं स्वचालित रूप से हाइड्रोडायनामिक सामूहिक व्यवहार का संकेत नहीं देती हैं। सिलिकॉन जैसे पदार्थों में, जो व्यापक MFP वितरण के कारण होने वाली मजबूत स्थानिक डीफेजिंग के कारण, सहसंबद्ध दोलन दब जाते हैं, जिससे परिवहन प्रसार (diffusive) लेकिन गैर-स्थानीय (nonlocal) प्रतीत होता है।
• व्यावहारिक उपयोगिता: यह ढांचा जटिल, गैर-समान प्रणालियों (नैनो-डिवाइस, इंटरफेस, विस्कृत पदार्थ) में ऊष्मा परिवहन को मॉडल करने का एक व्यवस्थित तरीका प्रदान करता है जहाँ पारंपरिक मॉडल विफल हो जाते हैं। यह सुझाव देता है कि "प्रभावी" पैरामीटर (जैसे कि एक एकल गैर-स्थानीय लंबाई) केवल अधिक जटिल अंतर्निखंड कर्नेल के कोर्स-ग्रेन्ड प्रोजेक्शन हैं।
• भविष्य की दिशाएं: यह शोध पत्र विस्कृत प्रणालियों के लिए सीधे साम्यावस्था मॉलिक्यूलर डायनेमिक्स से कर्नेल की गणना करने और गैर-स्थानीय संरचनात्मक संबंधों के लिए सुसंगत सीमा शर्तें विकसित करने का मार्ग प्रशस्त करता है।

निष्कर्षतः, यह शोध पत्र ऊष्मा चालन के लिए स्पैटियोटेम्पोरल हीट-फ्लक्स कोरिलेशन फंक्शन को मौलिक संरचनात्मक नियम के रूप में स्थापित करता है, जो विविध परिवहन क्षेत्रों को एक एकल, कठोर सांख्यिकीय ढांचे के तहत एकीकृत करता है।