ऊष्मा स्थानांतरण भौतिकी

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

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

परिचय
ऊष्मा कणों की तापमान-निर्भर गति से जुड़ी तापीय ऊर्जा है। ऊष्मा अंतरण विश्लेषण में प्रयुक्त अतिसूक्ष्म आयतन के लिए मैक्रोस्कोपिक ऊर्जा समीकरण है $$\nabla \cdot \mathbf{q} = -\rho c_p \frac{\partial T}{\partial t} + \sum_{i,j} \dot s_{i-j},$$ जहाँ $q$ ऊष्मा प्रवाह वेक्टर है, $−ρc_{p}(∂T/∂t)$ आंतरिक ऊर्जा ($ρ$ घनत्व है, $c_{p}$ स्थिर दबाव पर विशिष्ट ताप क्षमता है, $T$ तापमान है और $t$ समय है) का अस्थायी परिवर्तन है, और $$\dot s$$ थर्मल ऊर्जा ($i$ और $j$ प्रमुख ऊर्जा वाहकों के लिए हैं) से ऊर्जा रूपांतरण है। इसलिए ये शब्द ऊर्जा परिवहन, संचयन और परिवर्तन का प्रतिनिधित्व करते हैं। ऊष्मा प्रवाह वेक्टर $q$ तीन मैक्रोस्कोपिक मौलिक मोड से बना है, जो थर्मल चालन ($q_{k} = −k∇T$, $k$: तापीय चालकता), संवहन ($q_{u} = ρc_{p}uT$, $u$: वेग), और विकिरण ($ \mathbf q_r = 2\pi \int_{0}^{\infty} \int_{0}^{\pi} \mathbf s I_{ph,\omega} \sin(\theta) d\theta \, d\omega$, $ω$: कोणीय आवृत्ति, $θ$ : ध्रुवीय कोण, $I_{ph,ω}$: वर्णक्रमीय, दिशात्मक विकिरण तीव्रता, $s$: यूनिट वेक्टर) है। अर्थात्, $q = q_{k} + q_{u} + q_{r}$.

एक बार ऊर्जा रूपांतरण और थर्मोफिजिकल गुणों की स्थिति और गतिकी ज्ञात हो जाने पर गर्मी हस्तांतरण के भाग्य का वर्णन उपरोक्त समीकरण द्वारा किया जाता है। इन परमाणु-स्तर के तंत्रों और गतिकी को ऊष्मा स्थानांतरण भौतिकी में संबोधित किया जाता है। सूक्ष्म तापीय ऊर्जा को प्रमुख ऊर्जा वाहक फोनन (p), इलेक्ट्रॉन (e), द्रव कण (f), और फोटॉन (ph) द्वारा संग्रहीत, परिवहन और परिवर्तित किया जाता है।

लंबाई और समय का पैमाना
पदार्थ के थर्मोफिजिकल गुण और प्रमुख वाहकों के मध्य परस्पर क्रिया और ऊर्जा विनिमय की गतिशीलता परमाणु-स्तर के विन्यास और अंतःक्रिया पर आधारित होती है। तापीय चालकता जैसे परिवहन गुणों की गणना पारंपरिक और क्वांटम यांत्रिकी का उपयोग करके इन परमाणु-स्तर के गुणों से की जाती है। प्रमुख वाहकों की क्वांटम अवस्थाएँ (उदाहरण के लिए संवेग, ऊर्जा) श्रोडिंगर समीकरण (जिसे प्रथम सिद्धांत या एबी इनिटियो कहा जाता है) से प्राप्त की जाती हैं और इंटरैक्शन दर (कैनेटिक्स के लिए) की गणना क्वांटम अवस्थाओं और क्वांटम पर्टर्बेशन सिद्धांत (क्वांटम यांत्रिकी) का उपयोग करके की जाती है। फर्मी स्वर्ण नियम के रूप में तैयार किया गया)। एब इनिटियो (शुरुआत से लैटिन) सॉल्वर (सॉफ्टवेयर) की विविधता मौजूद है (उदाहरण के लिए, एबिनिट, कैस्टेप, गाऊसी (सॉफ्टवेयर), क्यू केम,  एस्प्रेसो जितना , सिएस्टा (कंप्यूटर प्रोग्राम), वीएएसपी, WIEN2k)। आंतरिक कोश (कोर) में इलेक्ट्रॉन गर्मी हस्तांतरण में शामिल नहीं होते हैं, और आंतरिक-कोश इलेक्ट्रॉनों के बारे में उचित अनुमान से गणना बहुत कम हो जाती है। क्वांटम उपचार, जिसमें संतुलन और नॉनक्विलिब्रियम एब इनिटियो आणविक गतिशीलता (एमडी) शामिल हैं, जिसमें बड़ी लंबाई और समय शामिल है, गणना संसाधनों द्वारा सीमित हैं, इसलिए सरलीकृत मान्यताओं के साथ विभिन्न वैकल्पिक उपचारों और बल गतिकी का उपयोग किया गया है। पारंपरिक (न्यूटोनियन) एमडी में, परमाणु या अणु (कण) की गति अनुभवजन्य या प्रभावी अंतःक्रिया क्षमता पर आधारित होती है, जो बदले में एबी इनिटियो गणना के वक्र-फिट या थर्मोफिजिकल गुणों के वक्र-फिट पर आधारित हो सकती है। अनुरूपित कणों के समुच्चय से, स्थैतिक या गतिशीलता थर्मल गुण या बिखरने की दर प्राप्त होती है। अभी भी बड़े लंबाई के पैमाने पर (मेसोस्केल, जिसमें कई माध्य मुक्त पथ शामिल हैं)बोल्ट्ज़मैन परिवहन समीकरण समीकरण (बीटीई) लागू किया जाता है जो पारंपरिक हैमिल्टनियन-सांख्यिकीय यांत्रिकी पर आधारित है। बीटीई स्थिति और गति वैक्टर (एक्स, पी) के संदर्भ में कण राज्यों पर विचार करता है और इसे राज्य व्यवसाय संभावना के रूप में दर्शाया जाता है। व्यवसाय में संतुलन वितरण (ज्ञात बोसॉन, फ़र्मियन और मैक्सवेल-बोल्ट्ज़मैन कण) हैं और ऊर्जा (गर्मी) का परिवहन किसी भी संतुलन (प्रेरक बल या क्षमता के कारण) के कारण होता है। परिवहन के केंद्र में बिखरने की भूमिका है जो वितरण को संतुलन की ओर मोड़ती है। प्रकीर्णन संबंध समय या माध्य मुक्त पथ द्वारा प्रस्तुत किया जाता है। विश्राम का समय (या इसका व्युत्क्रम जो अंतःक्रिया दर है) अन्य गणनाओं (अब इनिटियो या एमडी) या अनुभवजन्य रूप से पाया जाता है। बीटीई को मोंटे कार्लो विधि आदि से संख्यात्मक रूप से हल किया जा सकता है। लंबाई और समय के पैमाने के आधार पर, उपचार का उचित स्तर (ab initio, MD, या BTE) चुना जाता है। ऊष्मा स्थानांतरण भौतिकी विश्लेषण में थर्मल ऊर्जा संचयन, परिवहन और परिवर्तन से संबंधित राज्यों और गतिज के साथ कई पैमाने शामिल हो सकते हैं (उदाहरण के लिए, एबी इनिटियो या पारंपरिक एमडी से इंटरैक्शन दर का उपयोग करके बीटीई)।

तो, ऊष्मा स्थानांतरण भौतिकी पारंपरिक और क्वांटम यांत्रिक दृष्टिकोण से चार प्रमुख ऊर्जा वहन और उनकी गतिकी को कवर करती है। यह निम्न-आयामीता और आकार प्रभावों सहित मल्टीस्केल (एबी इनिटियो, एमडी, बीटीई और मैक्रोस्केल) विश्लेषण को सक्षम बनाता है।

फ़ोनोन
फोनन (मात्राबद्ध जाली कंपन तरंग) केंद्रीय थर्मल ऊर्जा वाहक है जो गर्मी क्षमता (समझदार गर्मी संचयन) और संघनित चरण में प्रवाहकीय गर्मी हस्तांतरण में योगदान देता है, और थर्मल ऊर्जा रूपांतरण में बहुत महत्वपूर्ण भूमिका निभाता है। इसके परिवहन गुणों को फोनन चालकता टेंसर K द्वारा दर्शाया जाता हैp (डब्ल्यू/एम-के, फूरियर कानून क्यू सेk,p = -केp⋅∇ टी) थोक सामग्री के लिए, और फोनन सीमा प्रतिरोध एआरp,b[के/(डब्ल्यू/एम2)] ठोस इंटरफेस के लिए, जहां ए इंटरफ़ेस क्षेत्र है। फोनन विशिष्ट ऊष्मा क्षमता cv,p(J/kg-K) में क्वांटम प्रभाव शामिल है। फोनन से जुड़ी तापीय ऊर्जा रूपांतरण दर इसमें शामिल है $$\dot{s}_{i\mbox{-}j}$$. ऊष्मा स्थानांतरण भौतिकी वर्णन और भविष्यवाणी करती है, सीv,p, 'क'p, आरp,b(या संचालन जीp,b) और $$\dot{s}_{i\mbox{-}j}$$, परमाणु-स्तर के गुणों पर आधारित।

संतुलन क्षमता के लिए ⟨φ⟩o एन परमाणुओं वाले सिस्टम में, कुल क्षमता ⟨φ⟩ संतुलन पर टेलर श्रृंखला विस्तार द्वारा पाई जाती है और इसे दूसरे डेरिवेटिव (हार्मोनिक सन्निकटन) द्वारा अनुमानित किया जा सकता है $$\begin{align} \langle\varphi\rangle &= \langle\varphi\rangle_\mathrm{o} + \left.\sum_i\sum_\alpha\frac{\partial\langle\varphi\rangle}{\partial d_{i\alpha}}\right|_\mathrm{o}d_{i\alpha} + \left.\frac{1}{2}\sum_{i,j}\sum_{\alpha,\beta}\frac{\partial^2\langle\varphi\rangle}{\partial d_{i\alpha}\partial d_{j\beta}}\right|_\mathrm{o}d_{i\alpha}d_{j\beta} + \left.\frac{1}{6}\sum_{i,j,k}\sum_{\alpha,\beta,\gamma}\frac{\partial^3\langle\varphi\rangle}{\partial d_{i\alpha}\partial d_{j\beta}\partial d_{k\gamma}}\right|_\mathrm{o} d_{i\alpha}d_{j\beta}d_{k\gamma}+ \cdots \\ &\approx \langle\varphi\rangle_\mathrm{o} + \frac{1}{2}\sum_{i,j}\sum_{\alpha,\beta}\Gamma_{\alpha\beta}d_{i\alpha}d_{j\beta}, \end{align} $$ जहां घi परमाणु i का विस्थापन वेक्टर है, और Γ विभव के दूसरे क्रम के व्युत्पन्न के रूप में स्प्रिंग (या बल) स्थिरांक है। परमाणुओं के विस्थापन के संदर्भ में जाली कंपन के लिए गति का समीकरण ['डी'(जेएल,टी): समय टी पर एल-वें इकाई सेल में जे-वें परमाणु का विस्थापन वेक्टर] है $$m_j\frac{d^2\mathbf{d}(jl,t)}{dt^2} = -\sum_{j'l'} \boldsymbol{\Gamma} \binom{j \ j^\prime}{l \ l'}\cdot \mathbf{d} (j' l', T), $$ जहां m परमाणु द्रव्यमान है और 'Γ' बल स्थिरांक टेंसर है। परमाणु विस्थापन सामान्य मोड का योग है ['s'α: मोड α, ω का यूनिट वेक्टरp: तरंग की कोणीय आवृत्ति, और 'κ'p: तरंग वेक्टर]। इस समतल-तरंग विस्थापन का उपयोग करते हुए, गति का समीकरण आइगेनवैल्यू समीकरण बन जाता है $$\mathbf{M} \omega_p^2 (\boldsymbol{\kappa}_p,\alpha) \mathbf{s}_\alpha(\boldsymbol{\kappa}_p) = \mathbf{D} (\boldsymbol{\kappa}_p) \mathbf{s}_\alpha(\boldsymbol{\kappa}_p), $$ जहां M विकर्ण द्रव्यमान मैट्रिक्स है और D हार्मोनिक डायनेमिक मैट्रिक्स है। इस eigenvalue समीकरण को हल करने से कोणीय आवृत्ति ω के मध्य संबंध मिलता हैpऔर तरंग वेक्टर 'κ'p, और इस संबंध को फ़ोनन फ़ोनन#विक्षेपण संबंध कहा जाता है। इस प्रकार, फोनन फैलाव संबंध मैट्रिक्स एम और डी द्वारा निर्धारित किया जाता है, जो परमाणु संरचना और घटक परमाणुओं के मध्य बातचीत की ताकत पर निर्भर करता है (इंटरेक्शन जितना मजबूत होगा और परमाणु जितने हल्के होंगे, फोनन आवृत्ति उतनी ही अधिक होगी और बड़ी होगी) ढलान ''dωp/डीएमp). हार्मोनिक सन्निकटन के साथ फोनन प्रणाली का हैमिल्टनियन है $$\mathrm{H}_p = \sum_x \frac{1}{2m} \mathbf{p}^2(\mathbf{x}) + \frac{1}{2}\sum_{\mathbf{x},\mathbf{x}'}\mathbf{d}_i(\mathbf{x})D_{ij}(\mathbf{x}-\mathbf{x}')\mathbf{d}_j(\mathbf{x}'),$$ जहां घijपरमाणुओं i और j, और 'd' के मध्य गतिशील मैट्रिक्स तत्व हैi (डीj) i (j) परमाणु का विस्थापन है, और 'p' संवेग है। इससे और फैलाव संबंध के समाधान से, क्वांटम उपचार के लिए फोनन विनाश ऑपरेटर को परिभाषित किया गया है $$b_{\kappa,\alpha} = \frac{1}{N^{1/2}}\sum_{\kappa_p,\alpha} e^{-i(\boldsymbol{\kappa}_p\cdot\mathbf{x})}\mathbf{s}_\alpha(\boldsymbol{\kappa}_p)\cdot \left[\left(\frac{m\omega_{p,\alpha}}{2\hbar}\right)^{1/2}\mathbf{d}(\mathbf{x}) + i\left(\frac{1}{2\hbar m\omega_{p,\alpha}}\right)^{1/2}\mathbf{p}(\mathbf{x})\right],$$ जहां N, α द्वारा विभाजित सामान्य मोड की संख्या है और ħ कम प्लैंक स्थिरांक है। सृजन संचालिका संहार संचालिका का सहायक है, $$ b_{\kappa,\alpha}^\dagger = \frac{1}{N^{1/2}}\sum_{\kappa_p,\alpha} e^{i(\boldsymbol{\kappa}_p\cdot\mathbf{x})}\mathbf{s}_\alpha(\boldsymbol{\kappa}_p)\cdot \left[\left(\frac{m\omega_{p,\alpha}}{2\hbar} \right)^{1/2}\mathbf{d}(\mathbf{x})-i\left(\frac{1}{2\hbar m\omega_{p,\alpha}}\right)^{1/2} \mathbf{p}(\mathbf{x})\right].$$ बी के संदर्भ में हैमिल्टनियनκ,α†और बीκ,αएच हैp= एसκ,αभाईp,α[बीκ,α†बीκ,α+ 1/2] और बीκ,α†बीκ,αफोनन नंबर ऑपरेटर है. क्वांटम-हार्मोनिक ऑसिलेटर की ऊर्जा E हैp= एसκ,α [एफp(के,ए) + 1/2]सीओp,α('क'p), और इस प्रकार फोनन ऊर्जा की मात्रा ħωp.

फ़ोनन फैलाव संबंध ब्रिलोइन जोन (पारस्परिक स्थान में आदिम कोशिका के अन्दर का क्षेत्र) और राज्यों डी के फ़ोनन घनत्व के अन्दर सभी संभावित फ़ोनन मोड देता हैp(संभावित फोनन मोड की संख्या घनत्व)। फ़ोनन समूह वेग यूp,gफैलाव वक्र का ढलान है, dωp/डी'के'p. चूँकि फोनन बोसोन कण है, इसका अधिभोग बोस-आइंस्टीन वितरण {fpओ = [exp(ħωp/कBटी)-1]−1, केB: बोल्ट्ज़मान स्थिरांक}. राज्यों के फोनन घनत्व और इस अधिभोग वितरण का उपयोग करते हुए, फोनन ऊर्जा ई हैp(टी) = '∫'डीp(ओहp)एफp(ओहp,T)ħωpdωp, और फोनन घनत्व n हैp(टी) = '∫'डीp(ओहp)एफp(ओहp,T)dωp. फोनन ताप क्षमता सीv,p(ठोस सी मेंv,p= सीp,p, सीv,p: स्थिर-आयतन ताप क्षमता, सीp,p: स्थिर-दबाव ताप क्षमता) डेबी मॉडल (रैखिक फैलाव मॉडल) के लिए फोनन ऊर्जा का तापमान व्युत्पन्न है, है $$c_{v,p} = \left.\frac{dE_p}{dT}\right|_v = \frac{9k_\mathrm{B}}{m} \left(\frac{T}{T_D} \right)^3 n \int_0^{T_D/T} \frac{x^4 e^x}{\left(e^x - 1 \right)^2} dx \qquad (x = \frac{\hbar\omega}{k_\mathrm{B}T}),$$ जहां टीD डेबी मॉडल है, एम परमाणु द्रव्यमान है, और एन परमाणु संख्या घनत्व है (क्रिस्टल 3एन के लिए फोनन मोड की संख्या घनत्व)। इससे डेबाई टी3 नियम मिलता है|डेबाई टी3कम तापमान पर नियम और उच्च तापमान पर डुलोंग-पेटिट नियम।

गैसों के गतिज सिद्धांत से, प्रमुख वाहक की तापीय चालकता i (p, e, f और ph) है $$ k_i = \frac{1}{3} n_i c_{v,i}u_i\lambda_i,$$ कहां एनiवाहक घनत्व है और ताप क्षमता प्रति वाहक है, यूiवाहक गति और λ हैiमाध्य मुक्त पथ है (प्रकीर्णन घटना से पहले वाहक द्वारा तय की गई दूरी)। इस प्रकार, वाहक घनत्व, ताप क्षमता और गति जितनी अधिक होगी और प्रकीर्णन जितना कम होगा, चालकता उतनी ही अधिक होगी। फोनन के लिए λpफ़ोनों की अंतःक्रिया (प्रकीर्णन) गतिकी का प्रतिनिधित्व करता है और प्रकीर्णन विश्राम समय से संबंधित है τpया दर (= 1/τp) λ के माध्यम सेp= यूpτp. फोनन अन्य फोनन के साथ, और इलेक्ट्रॉनों, सीमाओं, अशुद्धियों आदि के साथ बातचीत करते हैं, और λpमैथिसेन के नियम के माध्यम से इन अंतःक्रिया तंत्रों को जोड़ता है। कम तापमान पर, सीमाओं द्वारा प्रकीर्णन प्रमुख होता है और तापमान में वृद्धि के साथ अशुद्धियों, इलेक्ट्रॉन और अन्य फ़ोनों के साथ अंतःक्रिया की दर महत्वपूर्ण हो जाती है, और अंत में T > 0.2T के लिए फ़ोनन-फ़ोनन प्रकीर्णन प्रमुख हो जाता है।D. में इंटरैक्शन दरों की समीक्षा की जाती है और इसमें क्वांटम गड़बड़ी सिद्धांत और एमडी शामिल हैं।

फैलाव और λ के संबंध में अनुमान के साथ कई चालकता मॉडल उपलब्ध हैंp. एकल-मोड विश्राम समय सन्निकटन (∂f) का उपयोग करनाp“/∂t|s = −fp“/tp) और गैस गतिज सिद्धांत, कैलावे फोनन (जाली) चालकता मॉडल के रूप में $$ k_{p,\mathbf{s}} = \frac{1}{8\pi^3}\sum_{\alpha}\int c_{v,p}\tau_p(\mathbf{u}_{p,g}\cdot\mathbf{s})^2d\kappa \ \ \ \ \ \text{ for component along } \mathbf{s},$$ $$ k_p = \frac{1}{6\pi^3}\sum_{\alpha}\int c_{v,p}\tau_p {u}_{p,g}^2\kappa^2d\kappa \ \ \ \ \ \ \ \ \text{for isotropic conductivity}.$$ डेबी मॉडल के साथ (एकल समूह वेग यूp,g, और ऊपर गणना की गई विशिष्ट ताप क्षमता), यह बन जाती है $$ k_p = \left(48\pi^2\right)^{1/3} \frac{k_\mathrm{B}^3 T^3}{a h_\mathrm{P}^2 T_\mathrm{D}} \int_0^{T/T_\mathrm{D}}\tau_p \frac{x^4 e^x}{\left(e^x-1\right) ^2}dx,$$ जहाँ a जालक स्थिरांक a = n है−1/3घन जाली के लिए, और n परमाणु क्रमांक घनत्व है। सुस्त फोनन चालकता मॉडल मुख्य रूप से ध्वनिक फोनन बिखरने (तीन-फोनन इंटरैक्शन) पर विचार करते हुए दिया गया है $$ k_p = k_{p,S} = \frac{3.1\times10^{12}\langle M\rangle V_a^{1/3}T_{D,\infty}^3}{T\langle\gamma_G^2\rangle N_o^{2/3}}\qquad \text{ high temperatures } ( T > 0.2 T_D,\text{ phonon-phonon scattering only)},$$ जहाँ $⟨M⟩$ आदिम कोशिका में परमाणुओं का औसत परमाणु भार है, Va=1/एन प्रति परमाणु औसत आयतन है, टीD,∞उच्च तापमान डिबाई तापमान है, टी तापमान है, एनo आदिम कोशिका में परमाणुओं की संख्या है, और ⟨γ2G⟩ उच्च तापमान पर ग्रुनेसेन स्थिरांक या पैरामीटर का मोड-औसत वर्ग है। इस मॉडल का शुद्ध गैर-धातु क्रिस्टल के साथ व्यापक रूप से परीक्षण किया गया है, और समग्र समझौता अच्छा है, यहां तक ​​कि जटिल क्रिस्टल के लिए भी।

बल गतिकी और परमाणु संरचना विचार के आधार पर, उच्च क्रिस्टलीय और मजबूत इंटरैक्शन वाली सामग्री, जो हल्के परमाणुओं (जैसे हीरे और ग्राफीन) से बनी होती है, में बड़ी फोनन चालकता होने की उम्मीद है। जाली का प्रतिनिधित्व करने वाली सबसे छोटी इकाई कोशिका में से अधिक परमाणु वाले ठोस में दो प्रकार के फोनन होते हैं, अर्थात् ध्वनिक और ऑप्टिकल। (ध्वनिक फोनन अपने संतुलन की स्थिति के बारे में परमाणुओं के चरण-चरण आंदोलन हैं, जबकि ऑप्टिकल फोनन जाली में आसन्न परमाणुओं के चरण-बाहर आंदोलन हैं।) ऑप्टिकल फोनन में उच्च ऊर्जा (आवृत्ति) होती है, लेकिन चालन गर्मी हस्तांतरण में छोटा योगदान होता है, उनके छोटे समूह वेग और अधिभोग के कारण।

हेटरो-संरचना सीमाओं के पार फोनन परिवहन (आर के साथ दर्शाया गया है)।p,b, इंटरफेशियल थर्मल प्रतिरोध) सीमा बिखरने के अनुमान के अनुसार ध्वनिक और फैलाना बेमेल मॉडल के रूप में तैयार किया गया है। बड़ा फोनन ट्रांसमिशन (छोटा आरp,b) उन सीमाओं पर होता है जहां सामग्री जोड़े में समान फ़ोनन गुण होते हैं (यूp, डीp, आदि), और अनुबंध में बड़े आरp,bतब होता है जब कुछ सामग्री दूसरे की तुलना में नरम (कम कट-ऑफ फोनन आवृत्ति) होती है।

इलेक्ट्रॉन
इलेक्ट्रॉन के लिए क्वांटम इलेक्ट्रॉन ऊर्जा अवस्थाएं इलेक्ट्रॉन क्वांटम हैमिल्टनियन का उपयोग करके पाई जाती हैं, जो आम तौर पर गतिज (-ħ) से बना होता है2∇2/2me) और संभावित ऊर्जा शब्द (φ)।e). परमाणु कक्षक, फ़ंक्शन (गणित) जो परमाणु में इलेक्ट्रॉन या इलेक्ट्रॉनों की जोड़ी के तरंग-जैसे व्यवहार का वर्णन करता है, इस इलेक्ट्रॉन हैमिल्टनियन के साथ श्रोडिंगर समीकरण से पाया जा सकता है। हाइड्रोजन जैसे परमाणु (नाभिक और इलेक्ट्रॉन) इलेक्ट्रोस्टैटिक क्षमता (कूलम्ब कानून) के साथ श्रोडिंगर समीकरण के बंद-रूप समाधान की अनुमति देते हैं। से अधिक इलेक्ट्रॉन वाले परमाणुओं या परमाणु आयनों के श्रोडिंगर समीकरण को इलेक्ट्रॉनों के मध्य कूलम्ब इंटरैक्शन के कारण विश्लेषणात्मक रूप से हल नहीं किया गया है। इस प्रकार, संख्यात्मक तकनीकों का उपयोग किया जाता है, और इलेक्ट्रॉन विन्यास को सरल हाइड्रोजन-जैसे परमाणु ऑर्बिटल्स (पृथक इलेक्ट्रॉन ऑर्बिटल्स) के उत्पाद के रूप में अनुमानित किया जाता है। एकाधिक परमाणुओं (नाभिक और उनके इलेक्ट्रॉन) वाले अणुओं में आणविक कक्षीय (एमओ, अणु में ऋणावेशित सूक्ष्म अणु का विन्यास तरंग-जैसे व्यवहार के लिए गणितीय कार्य) होता है, और परमाणु कक्षाओं के रैखिक संयोजन (एलसीएओ) जैसी सरलीकृत समाधान तकनीकों से प्राप्त होते हैं।. आणविक कक्षक का उपयोग रासायनिक और भौतिक गुणों की भविष्यवाणी करने के लिए किया जाता है, और उच्चतम व्याप्त आणविक कक्षक (HOMO/LUMO) और सबसे कम रिक्त आणविक कक्षक (HOMO/LUMO) के मध्य का अंतर अणुओं की उत्तेजित अवस्था का माप है।

धात्विक ठोसों की क्रिस्टल संरचना में, मुक्त इलेक्ट्रॉन मॉडल (शून्य क्षमता, φe= 0) संयोजकता इलेक्ट्रॉनों के व्यवहार के लिए प्रयोग किया जाता है। हालाँकि, क्रिस्टल संरचना | आवधिक जाली (क्रिस्टल) में, आवधिक क्रिस्टल क्षमता होती है, इसलिए इलेक्ट्रॉन हैमिल्टनियन बन जाता है $$ \mathrm{H}_e = - \frac{\hbar^2}{2m_e}\nabla^2 + \varphi_c(\mathbf{x}),$$ कहां एमeइलेक्ट्रॉन द्रव्यमान है, और आवधिक क्षमता φ के रूप में व्यक्त की जाती हैc(एक्स) = एसg φgexp[i('g'∙'x')] ('g': व्युत्क्रम जाली वेक्टर)। इस हैमिल्टनियन के साथ समय-स्वतंत्र श्रोडिंगर समीकरण (आइजेनवैल्यू समीकरण) के रूप में दिया गया है $$ \mathrm{H}_e \psi_{e,\mathbf{x}}(\mathbf{x}) = E_e(\boldsymbol{\kappa}_e) \psi_{e,\mathbf{x}}(\mathbf{x}),$$ जहां eigenfunction ψe,κइलेक्ट्रॉन तरंग फ़ंक्शन है, और eigenvalue Ee('क'e), इलेक्ट्रॉन ऊर्जा है (κe: इलेक्ट्रॉन वेववेक्टर)। वेववेक्टर, κ के मध्य संबंधe और ऊर्जा ईeइलेक्ट्रॉनिक बैंड संरचना प्रदान करता है। व्यवहार में, अनेक-निकाय समस्या के रूप में जाली | अनेक-निकाय प्रणालियों में क्षमता में इलेक्ट्रॉनों और नाभिक के मध्य परस्पर क्रिया शामिल होती है, लेकिन यह गणना बहुत जटिल हो सकती है। इस प्रकार, कई अनुमानित तकनीकों का सुझाव दिया गया है और उनमें से है घनत्व कार्यात्मक सिद्धांत (डीएफटी), पूर्ण इंटरैक्शन के बजाय स्थानिक रूप से निर्भर इलेक्ट्रॉनिक घनत्व के कार्यात्मक का उपयोग करता है। डीएफटी का व्यापक रूप से एबी इनिटियो सॉफ्टवेयर (ABINIT, CASTEP, क्वांटम एस्प्रेसो, SIESTA (कंप्यूटर प्रोग्राम), VASP, WIEN2k, आदि) में उपयोग किया जाता है। इलेक्ट्रॉन विशिष्ट ऊष्मा ऊर्जा अवस्थाओं और अधिभोग वितरण (फ़र्मी-डिराक आँकड़े) पर आधारित है। सामान्य तौर पर, इलेक्ट्रॉन की ताप क्षमता बहुत उच्च तापमान को छोड़कर छोटी होती है जब वे फोनन (जाली) के साथ थर्मल संतुलन में होते हैं। इलेक्ट्रॉन ठोस में, विशेष रूप से धातुओं में, ताप संचालन (आवेश वहन के अलावा) में योगदान करते हैं। ठोस में तापीय चालकता टेंसर विद्युत और फोनन तापीय चालकता टेंसरों का योग है 'K' = 'K'e + केp.

इलेक्ट्रॉन दो थर्मोडायनामिक बलों से प्रभावित होते हैं [आवेश से, ∇(EF/यह हैc) जहां ईF फर्मी स्तर और ई हैcप्राथमिक आवेश और तापमान प्रवणता है, ∇(1/T)] क्योंकि उनमें आवेश और तापीय ऊर्जा दोनों होती है, और इस प्रकार विद्युत धारा 'जे' होती है।e और ताप प्रवाह q को थर्मोइलेक्ट्रिक टेंसर (ए) के साथ वर्णित किया गया हैee, एet, एte, और एtt) ऑनसागर पारस्परिक संबंधों से जैसा $$ \mathbf{j}_e = \mathbf{A}_{ee}\cdot\nabla\frac{E_\mathrm{F}}{e_c} + \mathbf{A}_{et}\cdot\nabla\frac{1}{T} ,\ \ \text{and}$$ $$ \mathbf{q}= \mathbf{A}_{te}\cdot\nabla\frac{E_\mathrm{F}}{e_c} + \mathbf{A}_{tt}\cdot\nabla\frac{1}{T}.$$ इन समीकरणों को j में परिवर्तित करनाe विद्युत क्षेत्र के संदर्भ में समीकरण ईe और ∇T और 'q' समीकरण 'j' के साथe और ∇T, (आइसोट्रोपिक परिवहन के लिए अदिश गुणांक का उपयोग करते हुए, αee, एet, एte, और αttके बजाय एक'ee, एet, एte, और एtt) $$ \mathbf{j}_e = \alpha_{ee}\mathbf{e}_e - \frac{\alpha_{et}}{T^2}\nabla T \qquad (\mathbf{e}_e = \alpha_{ee}^{-1}\mathbf{j}_e+\frac{\alpha_{ee}^{-1}\alpha_{et}}{T^2}\nabla T), $$ $$ \mathbf{q}= \alpha_{te}\alpha_{ee}^{-1}\mathbf{j}_e-\frac{\alpha_{tt}-\alpha_{te}\alpha_{ee}^{-1}\alpha_{et}}{T^2}\nabla T.$$ विद्युत चालकता/प्रतिरोधकता σe(ओह−1m−1)/ पीe (Ω-m), विद्युत तापीय चालकता ke(डब्ल्यू/एम-के) और सीबेक/पेल्टियर गुणांक αS (वी/के)/एP (वी) को इस प्रकार परिभाषित किया गया है, $$ \sigma_e = \frac{1}{\rho_e}=\alpha_{ee}, \ \ k_e = \frac{\alpha_{tt}-\alpha_{te}\alpha_{ee}^{-1}\alpha_{et}}{T^2},\mathrm{and} \ \alpha_\mathrm{S} = \frac{\alpha_{et}\alpha_{ee}^{-1}}{T^2} \ \ (\alpha_\mathrm{S} = \alpha_\mathrm{P}T). $$ विभिन्न वाहक (इलेक्ट्रॉन, मैग्नन, फोनन और पोलरॉन) और उनकी परस्पर क्रियाएं सीबेक गुणांक को काफी हद तक प्रभावित करती हैं। सीबेक गुणांक को दो योगदानों, α के साथ विघटित किया जा सकता हैS = एS,pres + एS,trans, कहां αS,pres वाहक-प्रेरित एन्ट्रापी परिवर्तन में योगदान का योग है, अर्थात, αS,pres = एS,mix + एS,spin + एS,vib (एS,mix: मिश्रण की एन्ट्रॉपी, αS,spin: स्पिन एन्ट्रापी, और αS,vib: कंपन एन्ट्रापी)। अन्य योगदान αS,trans किसी वाहक को हिलाने में हस्तांतरित शुद्ध ऊर्जा को qT (q: वाहक आवेश) से विभाजित किया जाता है। सीबेक गुणांक में इलेक्ट्रॉन का योगदान अधिकतर α में होता हैS,pres. αS,mix आमतौर पर हल्के डोप किए गए अर्धचालकों में प्रमुख होता है। किसी प्रणाली में इलेक्ट्रॉन जोड़ने पर मिश्रण की एन्ट्रापी में परिवर्तन तथाकथित हेइक्स सूत्र है $$ \alpha_\mathrm{S,mix} = \frac{1}{q} \frac{\partial S_\mathrm{mix}}{\partial N} = \frac{k_\mathrm{B}}{q}\ln\left(\frac{1 - f_e^\mathrm{o}}{f_e^\mathrm{o}}\right),$$ जहाँ एफeओ = एन/एनaसाइटों (वाहक एकाग्रता) के लिए इलेक्ट्रॉनों का अनुपात है। रासायनिक क्षमता (μ) का उपयोग करते हुए, तापीय ऊर्जा (kBटी) और फर्मी फ़ंक्शन, उपरोक्त समीकरण को वैकल्पिक रूप, α में व्यक्त किया जा सकता हैS,mix = (केB/क्यू)[(ईe- μ)/(kBटी)]। सीबेक प्रभाव को स्पिन तक विस्तारित करते हुए, लौहचुंबकीय मिश्र धातु अच्छा उदाहरण हो सकता है। सीबेक गुणांक में योगदान, जो सिस्टम की स्पिन एन्ट्रापी को बदलने वाले इलेक्ट्रॉनों की उपस्थिति के परिणामस्वरूप होता है, α द्वारा दिया जाता हैS,spin = एसspin/क्यू = (केB/q)ln[(2s + 1)/(2s0 +1)], जहां एस0 और एस क्रमशः वाहक की अनुपस्थिति और उपस्थिति में चुंबकीय स्थल के शुद्ध स्पिन हैं। इलेक्ट्रॉनों के साथ कई कंपन प्रभाव भी सीबेक गुणांक में योगदान करते हैं। कंपन आवृत्तियों का नरम होना कंपन एन्ट्रापी में परिवर्तन उत्पन्न करता है, इसका उदाहरण है। कंपन एन्ट्रापी मुक्त ऊर्जा का नकारात्मक व्युत्पन्न है, अर्थात, $$ S_\mathrm{vib} = -\frac{\partial F_\mathrm{mix}}{\partial T} = 3Nk_\mathrm{B}T\int_0^\omega \left\{\frac{\hbar\omega}{2k_\mathrm{B}T}\coth\left(\frac{\hbar\omega}{2k_\mathrm{B}T}\right) - \ln \left[2\sinh\left(\frac{\hbar\omega}{2k_\mathrm{B}T}\right)\right] \right\}D_p(\omega)d\omega,$$ जहां घp(ω) संरचना के लिए फ़ोनन घनत्व की स्थिति है। उच्च तापमान सीमा और अतिशयोक्तिपूर्ण कार्यों की श्रृंखला विस्तार के लिए, उपरोक्त को α के रूप में सरल बनाया गया हैS,vib = (ΔSvib/क्यू) = (केB/क्यू)एसi(-देखनाi/ओi).

उपरोक्त ऑनसेगर फॉर्मूलेशन में प्राप्त सीबेक गुणांक मिश्रण घटक α हैS,mix, जो अधिकांश अर्धचालकों पर हावी है। हाई-बैंड गैप सामग्री जैसे बी में कंपन घटक13C2 बहुत महत्वपूर्ण है.

सूक्ष्म परिवहन को ध्यान में रखते हुए (परिवहन किसी संतुलन का परिणाम नहीं है), $$ \mathbf{j}_e = -\frac{e_c}{\hbar^3}\sum_p\mathbf{u}_e f_e^\prime = -\frac{e_c}{\hbar^3k_\mathrm{B}T}\sum_p\mathbf{u}_e\tau_e \left(-\frac{\partial f_e^\mathrm{o}}{\partial E_e}\right)(\mathbf{u}_e\cdot\mathbf{F}_{te}),$$ $$ \mathbf{q}=\frac{1}{\hbar^3}\sum_p(E_e-E_\mathrm{F})\mathbf{u}_ef_e^\prime = \frac{1}{\hbar^3k_\mathrm{B}T}\sum_p \mathbf{u}_e \tau_e \left(-\frac{\partial f_e^\mathrm{o}}{\partial E_e}\right)(E_e-E_\mathrm{F})(\mathbf{u}_e\cdot\mathbf{F}_{te}),$$ जहां तुमe इलेक्ट्रॉन वेग वेक्टर है, एफe(एफeo) इलेक्ट्रॉन नोक्विलिब्रियम (संतुलन) वितरण है, τeइलेक्ट्रॉन प्रकीर्णन समय है, ईeइलेक्ट्रॉन ऊर्जा है, और 'एफ'te ∇(E) से विद्युत और तापीय बल हैF/यह हैc) और ∇(1/T). जे के लिए थर्मोइलेक्ट्रिक गुणांक को सूक्ष्म परिवहन समीकरणों से संबंधित करनाeऔर क्यू, थर्मल, इलेक्ट्रिक और थर्मोइलेक्ट्रिक गुणों की गणना की जाती है। इस प्रकार, केeविद्युत चालकता σe और तापमान T के साथ बढ़ती है, जैसा कि विडेमैन-फ्रांज कानून प्रस्तुत करता है [ke/(पीeTe) = (1/3)(πkB/यह हैc)2= $0 W-Ω/K^{2}$]. इलेक्ट्रॉन परिवहन (σ के रूप में दर्शाया गया हैe) वाहक घनत्व n का फलन हैe,cऔर इलेक्ट्रॉन गतिशीलता μe(पीe= औरcne,cμe). एमeइलेक्ट्रॉन प्रकीर्णन दर द्वारा निर्धारित होता है $$\dot{\gamma}_e$$ (या विश्राम का समय, $$\tau_e = 1/\dot{\gamma}_e $$) अन्य इलेक्ट्रॉनों, फोनन, अशुद्धियों और सीमाओं के साथ बातचीत सहित विभिन्न इंटरैक्शन तंत्रों में।

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

द्रव कण
द्रव कण किसी भी रासायनिक बंधन को तोड़े बिना द्रव चरण (गैस, तरल या प्लाज्मा) में सबसे छोटी इकाई (परमाणु या अणु) है। द्रव कण की ऊर्जा को संभावित, इलेक्ट्रॉनिक, ट्रांसलेशनल, कंपनात्मक और घूर्णी ऊर्जा में विभाजित किया गया है। द्रव कण में ऊष्मा (थर्मल) ऊर्जा का संचयन तापमान पर निर्भर कण गति (अनुवादात्मक, कंपनात्मक और घूर्णी ऊर्जा) के माध्यम से होता है। इलेक्ट्रॉनिक ऊर्जा को केवल तभी शामिल किया जाता है जब तापमान तरल कणों को आयनित करने या भिन्न करने या अन्य इलेक्ट्रॉनिक संक्रमणों को शामिल करने के लिए पर्याप्त उच्च हो। द्रव कणों की ये क्वांटम ऊर्जा अवस्थाएँ उनके संबंधित क्वांटम हैमिल्टनियन का उपयोग करके पाई जाती हैं। ये हैं एचf,t = −(एच2/2m)∇2, एचf,v= −(एच2/2m)∇2 + Γx2/2 और एचf,r = −(एच2/2If)∇2ट्रांसलेशनल, वाइब्रेशनल और रोटेशनल मोड के लिए। (Γ: हुक का नियम, If: अणु के लिए जड़ता का क्षण)। हैमिल्टनियन से, परिमाणित द्रव कण ऊर्जा अवस्था ईfऔर विभाजन फलन (सांख्यिकीय यांत्रिकी) Zf[मैक्सवेल-बोल्ट्ज़मैन आँकड़ों के साथ|मैक्सवेल-बोल्ट्ज़मैन (एमबी) अधिभोग वितरण] के रूप में पाए जाते हैं यहाँ, जीfअध:पतन है, n, l, और j संक्रमणकालीन, कंपनात्मक और घूर्णी क्वांटम संख्याएँ हैं, Tf,vकंपन के लिए विशिष्ट तापमान है (= ħωf,v/कB, : कंपन आवृत्ति), और टीf,rघूर्णी तापमान है [= ħ2/(2आईfkB)]. औसत विशिष्ट आंतरिक ऊर्जा Z के माध्यम से विभाजन फ़ंक्शन से संबंधित हैf, $$ e_f = (k_\mathrm{B}T^2/m)(\partial \mathrm{ln}Z_f/\partial T)|_{N,V}.$$ ऊर्जा अवस्थाओं और विभाजन फ़ंक्शन के साथ, द्रव कण विशिष्ट ताप क्षमता cv,fविभिन्न गतिज ऊर्जाओं के योगदान का योग है (गैर-आदर्श गैस के लिए संभावित ऊर्जा भी जोड़ी जाती है)। क्योंकि अणुओं में स्वतंत्रता की कुल डिग्री परमाणु विन्यास द्वारा निर्धारित होती है, cv,fकॉन्फ़िगरेशन के आधार पर भिन्न-भिन्न सूत्र हैं,
 * अनुवादात्मक $$ E_{f,t,n} = \frac{\pi^2\hbar^2}{2m} \left(\frac{n_x^2}{L^2}+\frac{n_y^2}{L^2}+\frac{n_z^2}{L^2}\right) \ \ \ \text{and} \ \ \ Z_{f,t} \sum_{i = 0}^\infty g_{f,t,i} \exp \left(-\frac{E_{f,t,i}}{k_\mathrm{B}T}\right) = V \left(\frac{m k_\mathrm{B}T}{2\pi\hbar^2}\right)^{3/2},$$
 * कंपनात्मक $$ E_{f,v,l} = \hbar\omega_{f,v}\left(1 + \frac{1}{2}\right) \ \ \text{and} \ \ Z_{f,v}\sum_{j = 0}^\infty \exp\left[-\left(l+\frac{1}{2}\right)\frac{\hbar\omega_{f,v}}{k_\mathrm{B}T}\right] = \frac{\exp(-T_{f,v}/2T)}{1-\exp(-T_{f,v}/T)},$$
 * घूर्णी $$ E_{f,r,j} = \frac{\hbar^2}{2I_f} \ \ \text{and} \ \ Z_{f,r}\sum_{j = 0}^\infty (2j+1)\exp \left[-\frac{-\hbar^2j(j+1)}{2I_f k_\mathrm{B}T}\right] \approx \frac{T}{T_{f,r}} \left(1 + \frac{T_{f,r}}{3T} + \frac{T_{f,r}^2}{15T^2}+ \cdots\right),$$
 * कुल $$ E_{f} = \sum_i E_{f,i} = E_{f,t} + E_{f,v} + E_{f,r} + \dots \ \ \text{and} \ \ Z_{f}=\prod_{i}Z_{f,i} = Z_{f,t}Z_{f,v}Z_{f,r}\dots .$$

जहां आरgगैस स्थिरांक है (= NAkB, एनA: एवोगैड्रो स्थिरांक) और एम आणविक द्रव्यमान (किलो/किलोमीटर) है। (बहुपरमाणुक आदर्श गैस के लिए, एनo अणु में परमाणुओं की संख्या है।) गैस में, स्थिर दबाव विशिष्ट ताप क्षमता सीp,fइसका मान बड़ा है और अंतर तापमान T, वॉल्यूमेट्रिक थर्मल विस्तार गुणांक β और इज़ोटेर्मल संपीड़ितता κ [c पर निर्भर करता है।p,f- सीv,f= टीβ2/(आरfके), आरf: द्रव घनत्व]। सघन तरल पदार्थों के लिए कणों के मध्य परस्पर क्रिया (वैन डेर वाल्स इंटरेक्शन) को शामिल किया जाना चाहिए, और सीv,fऔर सीp,fतदनुसार परिवर्तन होगा. कणों की शुद्ध गति (गुरुत्वाकर्षण या बाहरी दबाव के तहत) संवहन ऊष्मा प्रवाह 'q' को जन्म देती हैu = पीfcp,fमेंfटी. चालन ताप प्रवाह 'क्यू'kआदर्श गैस के लिए गैस गतिज सिद्धांत या बोल्ट्ज़मैन परिवहन समीकरणों से प्राप्त किया जाता है, और तापीय चालकता होती है $$ k_{f} = \tfrac{1}{3}n_f c_{p,f}\langle u_f^2\rangle\tau_{f\mbox{-}f},$$ जहां तुमf2⟩1/2आरएमएस (मूल माध्य वर्ग) थर्मल वेग (3k) हैBएमबी वितरण फ़ंक्शन से टी/एम, एम: परमाणु द्रव्यमान) और τf-fविश्राम का समय है (या अंतर्टकराव समय अवधि) [(21/2π डी2nf⟨मेंf⟩)−1गैस गतिज सिद्धांत से, ⟨uf⟩: औसत तापीय गति (8kBटी/πm)1/2, d: द्रव कण (परमाणु या अणु) का टकराव व्यास, nf: द्रव संख्या घनत्व]।
 * मोनोआटोमिक आदर्श गैस $$ c_{v,f} = \left.\frac{\partial e_f}{\partial T}\right|_V = \frac{3R_g}{2M},$$
 * द्विपरमाणुक आदर्श गैस $$ c_{v,f} = \frac{R_g}{M} \left\{ \frac{3}{2} + \left(\frac{T_{f,v}}{T}\right)^2 \frac{\exp(T_{f,v,i}/T)}{[\exp(T_{f,v,i}/T)-1]^2} + 1 + \frac{2}{15} \left(\frac{T_{f,v}}{T}\right)^2 \right\},$$
 * अरैखिक, बहुपरमाणुक आदर्श गैस $$ c_{v,f} = \frac{R_g}{M} \left\{3+ \sum_{j=1}^{3N_o-6} \left(\frac{T_{f,v}}{T}\right)^2 \frac{\exp(T_{f,v,i}/T)}{[\exp(T_{f,v,i}/T)-1]^2} \right\} .$$

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

द्रव कण अन्य प्रमुख कणों के साथ परस्पर क्रिया कर सकते हैं। कंपन या घूर्णी मोड, जिनमें अपेक्षाकृत उच्च ऊर्जा होती है, फोटॉन के साथ बातचीत के माध्यम से उत्तेजित या क्षय होते हैं। गैस लेजर द्रव कणों और फोटॉन के मध्य इंटरेक्शन बल गतिकी को नियोजित करते हैं, और सीओ में लेजर कूलिंग पर भी विचार किया गया है2 गैस लेजर. इसके अलावा, तरल पदार्थ के कण ठोस सतहों (फिसिसोरेशन और केमिसोरेशन) पर सोख सकते हैं, और सोखने वाले (द्रव कण) में कुंठित कंपन मोड ई बनाकर क्षय हो जाते हैं−-ज+जोड़े या फ़ोनन। इन अंतःक्रिया दरों की गणना द्रव कण और फर्मी गोल्डन नियम पर एब इनिटियो गणना के माध्यम से भी की जाती है।

फोटॉन
फोटॉन विद्युतचुंबकीय विकिरण का क्वांटा है|विद्युतचुंबकीय (ईएम) विकिरण और थर्मल विकिरण के लिए ऊर्जा वाहक है। ईएम तरंग पारंपरिक मैक्सवेल के समीकरणों द्वारा नियंत्रित होती है, और विद्युत चुम्बकीय क्षेत्र की मात्रा का उपयोग ब्लैक-बॉडी विकिरण (विशेष रूप से पराबैंगनी आपदा को समझाने के लिए) जैसी घटनाओं के लिए किया जाता है। कोणीय आवृत्ति ω की क्वांटा ईएम तरंग (फोटॉन) ऊर्जाphई हैph = hωph, और बोस-आइंस्टीन वितरण फ़ंक्शन (एफ) का अनुसरण करता हैph). परिमाणित विकिरण क्षेत्र (द्वितीय परिमाणीकरण) के लिए फोटॉन हैमिल्टनियन है $$ \mathrm{H}_{ph} = \frac{1}{2} \int \left(\varepsilon_\mathrm{o}\mathbf{e}_e^2 + \mu_\mathrm{o}^{-1}\mathbf{b}_e^2\right)dV = \sum_\alpha \hbar \omega_{ph,\alpha} \left(c_\alpha^\dagger c_\alpha + \frac{1}{2}\right),$$ कहां ईe और बीe ईएम विकिरण के विद्युत और चुंबकीय क्षेत्र हैं, εo और μo मुक्त-स्थान पारगम्यता और पारगम्यता हैं, वी इंटरैक्शन वॉल्यूम है, ωph,αα मोड और c के लिए फोटॉन कोणीय आवृत्ति हैα†और सीαफोटॉन निर्माण और विनाश संचालक हैं। वेक्टर क्षमता 'ए'e ईएम क्षेत्रों की (उदाe = −∂ae/∂t और 'बी'e = ∇×ae) है $$ \mathbf{a}_{e} (\mathbf{x},t) = \sum_\alpha \left(\frac{\hbar}{2\varepsilon_\mathrm{o}\omega_{ph,\alpha}V}\right)^{1/2} \mathbf{s}_{ph,\alpha} \left(c_\alpha e^{i \boldsymbol{\kappa}_\alpha \cdot \mathbf{x}} + c_\alpha^\dagger e^{-i\boldsymbol{\kappa}_\alpha\cdot\mathbf{x}}\right), $$ जहाँ एसph,α इकाई ध्रुवीकरण वेक्टर है, κα तरंग सदिश है.

विभिन्न प्रकार के फोटॉन उत्सर्जन के मध्य ब्लैकबॉडी विकिरण इंटरफोटॉन इंटरैक्शन के बिना थर्मल ऊर्जा वितरण के साथ फोटॉन गैस मॉडल को नियोजित करता है। रैखिक फैलाव संबंध (अर्थात्, फैलाव रहित) से, चरण और समूह गति बराबर हैं (यू)।ph= डी ωph/dk = ωph/के, यूph: फोटॉन गति) और डिबाई (फैलाव रहित फोटॉन के लिए प्रयुक्त) राज्यों का घनत्व डी हैph,b,ωdω = ωph2dωph/पी2uph3. डी के साथph,b,ωऔर संतुलन वितरण एफph, फोटॉन ऊर्जा वर्णक्रमीय वितरण डी.आईb,ωया डी.आईb,λ(एलph: तरंग दैर्ध्य) और कुल उत्सर्जक शक्ति ईbके रूप में व्युत्पन्न हैं $$ dI_{b,\omega} = \frac{D_{ph,b,\omega}f_{ph}u_{ph}d\omega_{ph}}{4\pi} =\frac{\hbar\omega_{ph}^3}{4\pi^3u_{ph}^2} \frac{1}{e^{\hbar\omega_{ph}/k_\mathrm{B}T}-1} d\omega_{ph} \ \text{or} \ d I_{b,\lambda} = \frac{4\pi\hbar u_{ph}^2 d \lambda_{ph}}{\lambda_{ph}^5(e^{2\pi\hbar u_{ph} / \lambda_{ph}k_\mathrm{B}T}-1)} $$ (प्लैंक का नियम), $$ E_b = \int_0^\infty d E_{b,\lambda} = \sigma_\mathrm{SB}T^4\ \text{, where} \ \sigma_\mathrm{SB} = \frac{\pi^2 k_\mathrm{B}^4}{60 \hbar^3 u_{ph}^2} $$ (स्टीफ़न-बोल्ट्ज़मैन कानून)।

ब्लैकबॉडी विकिरण की तुलना में, लेजर उत्सर्जन में उच्च दिशात्मकता (छोटा ठोस कोण ΔΩ) और वर्णक्रमीय शुद्धता (संकीर्ण बैंड Δω) होती है। इलेक्ट्रॉनिक ऊर्जा अवस्थाओं के मध्य गुंजयमान संक्रमण (उत्तेजित उत्सर्जन) के आधार पर लेज़रों की रेंज दूर-अवरक्त से लेकर एक्स-रे/γ-किरणों तक होती है। निकट-क्षेत्र विकिरण ताप स्थानांतरण|ऊष्मीय रूप से उत्तेजित द्विध्रुवों और अन्य विद्युत/चुंबकीय संक्रमणों से निकट-क्षेत्र विकिरण उत्सर्जन स्थलों से कम दूरी (तरंग दैर्ध्य के क्रम) के अन्दर बहुत प्रभावी होता है। फोटॉन कण गति के लिए बीटीई पीph = hωphएस/यूphदिशा के साथ-साथ अवशोषण/उत्सर्जन का अनुभव हो रहा है $$ \textstyle \dot{s}_{f,ph-e}\ $$ (=मेंphσph,ω[एफph(ओहph,टी) - एफph('एस')], पीph,ω: वर्णक्रमीय अवशोषण गुणांक), और पीढ़ी/निष्कासन $$ \textstyle \dot{s}_{f,ph,i}$$, है $$ \frac{\partial f_{ph}}{\partial t} + u_{ph}\mathbf{s}\cdot\nabla f_{ph} = \left.\frac{\partial f_{ph}}{\partial t}\right|_s + u_{ph}\sigma_{ph,\omega}[f_{ph}(\omega_{ph},T)-f_{ph}(\mathbf{s})]+ \dot{s}_{f,ph,i}. $$ विकिरण की तीव्रता के संदर्भ में (Iph,ω= यूphfphभाईphDph,ω/4पी, डीph,ω: राज्यों का फोटॉन घनत्व), इसे विकिरण हस्तांतरण (ईआरटी) का समीकरण कहा जाता है $$ \frac{\partial I_{ph,\omega}(\omega_{ph}, \mathbf{s})}{u_{ph} \partial t} + \mathbf{s}\cdot\nabla I_{ph,\omega} (\omega_{ph},\mathbf{s}) = \left.\frac{\partial I_{ph,\omega}(\omega_{ph}, \mathbf{s})}{u_{ph}\partial t}\right|_s + \sigma_{ph,\omega}[I_{ph,\omega}(\omega_{ph},T)-I_{ph}(\omega_{ph},\mathbf{s})]+ \dot{s}_{ph,i}. $$शुद्ध विकिरणीय ऊष्मा प्रवाह वेक्टर है $ \mathbf{q}_r = \mathbf{q}_{ph} = \int_0^\infty\int_{4\pi}\mathbf{s} I_{ph,\omega}d \Omega d\omega.$ आइंस्टीन गुणांक से, वर्णक्रमीय अवशोषण गुणांक σph,ωईआरटी में है, $$ \sigma_{ph,\omega} = \frac{\hbar\omega\dot{\gamma}_{ph,a}n_e}{u_{ph}},$$ जहाँ $$\dot{\gamma}_{ph,a}$$ अंतःक्रिया संभाव्यता (अवशोषण) दर या परमाणु वर्णक्रमीय रेखा बी है12(जे−1m3s−1), जो विकिरण क्षेत्र की प्रति इकाई वर्णक्रमीय ऊर्जा घनत्व (1: जमीनी अवस्था, 2: उत्तेजित अवस्था), और n प्रति इकाई समय की संभावना देता हैeइलेक्ट्रॉन घनत्व (जमीनी अवस्था में) है। इसे संक्रमण द्विध्रुव आघूर्ण 'μ' का उपयोग करके प्राप्त किया जा सकता हैeएफजीआर और आइंस्टीन गुणांक के मध्य संबंध के साथ। औसत σph,ωω से अधिक औसत फोटॉन अवशोषण गुणांक σ देता हैph.

L लंबाई के वैकल्पिक रूप से मोटे माध्यम के मामले में, अर्थात्, σphएल >> 1, और गैस गतिज सिद्धांत का उपयोग करते हुए, फोटॉन चालकता kph16σ हैSBT3/3σph(पीSB: स्टीफ़न-बोल्ट्ज़मान स्थिरांक, σph: औसत फोटॉन अवशोषण), और फोटॉन ताप क्षमता एनphcv,ph16σ हैSBT3/uph.

फोटॉन में ऊर्जा की सबसे बड़ी श्रृंखला होती है और यह विभिन्न प्रकार के ऊर्जा रूपांतरणों में केंद्रीय होता है। फोटॉन विद्युत और चुंबकीय संस्थाओं के साथ परस्पर क्रिया करते हैं। उदाहरण के लिए, विद्युत द्विध्रुव जो बदले में ऑप्टिकल फोनन या द्रव कण कंपन, या इलेक्ट्रॉनिक संक्रमण के संक्रमण द्विध्रुव क्षणों से उत्तेजित होते हैं। ऊष्मा स्थानांतरण भौतिकी में, फोनन के इंटरेक्शन बल गतिकी का इलाज परटर्बेशन सिद्धांत (फर्मी गोल्डन रूल) और इंटरेक्शन हैमिल्टनियन का उपयोग करके किया जाता है। फोटॉन-इलेक्ट्रॉन इंटरैक्शन है $$ \mathrm{H}_{ph-e} = -\frac{e_c}{m_e} \left(a + a^\dagger\right)\mathbf{a}_e\cdot\mathbf{p}_e = -\left(\frac{\hbar\omega_{ph,\alpha}}{2\varepsilon_o V}\right)^{1/2} (\mathbf{s}_{ph,\alpha}\cdot e_c \mathbf{x}_e)\left(a + a^\dagger\right)\left(ce^{i\mathrm{\kappa}\cdot\mathrm{x}}+c^\dagger e^{-i\mathrm{\kappa}\cdot\mathrm{x}}\right), $$ जहां पीe द्विध्रुव आघूर्ण सदिश है और a†और ए इलेक्ट्रॉन की आंतरिक गति का निर्माण और विनाश है। फोटॉन टर्नरी इंटरैक्शन में भी भाग लेते हैं, उदाहरण के लिए, फोनन-सहायता वाले फोटॉन अवशोषण/उत्सर्जन (इलेक्ट्रॉन ऊर्जा स्तर का संक्रमण)। द्रव कणों में कंपन मोड फोटॉन उत्सर्जित या अवशोषित करके क्षय या उत्तेजित हो सकता है। उदाहरण ठोस और आणविक गैस लेजर शीतलन हैं। ईएम सिद्धांत के साथ पहले सिद्धांतों के आधार पर एबी इनिटियो गणनाओं का उपयोग करते हुए, विभिन्न विकिरण गुण जैसे कि ढांकता हुआ फ़ंक्शन (विद्युत पारगम्यता, ε)e,ω), वर्णक्रमीय अवशोषण गुणांक (σph,ω), और जटिल अपवर्तन सूचकांक (एमω), पदार्थ में फोटॉन और विद्युत/चुंबकीय संस्थाओं के मध्य विभिन्न इंटरैक्शन के लिए गणना की जाती है। उदाहरण के लिए, काल्पनिक भाग (εe,c,ω) जटिल ढांकता हुआ फ़ंक्शन (εe,ω= ईe,r,ω+ मैं ईe,c,ω) बैंडगैप में इलेक्ट्रॉनिक संक्रमण के लिए है $$ \varepsilon_{e,c,\omega} = \frac{4\pi^2}{\omega^2V}\sum_{i\isin \mathrm{VB},j\isin \mathrm{CB}}\sum_{\kappa} w_\kappa |p_{ij}|^2 \delta(E_{\kappa,j}-E_{\kappa,i}-\hbar\omega), $$ जहां V इकाई-सेल आयतन है, VB और CB वैलेंस और चालन बैंड को दर्शाते हैं, wκκ-बिंदु और पी से जुड़ा वजन हैijसंक्रमण गति मैट्रिक्स तत्व है। वास्तविक भाग ε हैe,r,ωε से प्राप्त होता हैe,c,ωक्रेमर्स-क्रोनिग संबंध का उपयोग करना $$ \varepsilon_{e,r,\omega} = 1 + \frac{4}{\pi}\mathbb{P}\int_{0}^\infty \mathrm{d}\omega'\frac{\omega'\varepsilon_{e,c,\omega'}}{\omega'^2-\omega^2}.$$ यहाँ, $$\mathbb{P}$$ कॉची प्रमुख मूल्य को दर्शाता है।

अन्य उदाहरण में, सुदूर आईआर क्षेत्रों के लिए जहां ऑप्टिकल फोनन शामिल हैं, ढांकता हुआ फ़ंक्शन (εe,ω) के रूप में गणना की जाती है $$ \frac{\varepsilon_{e,\omega}}{\varepsilon_{e,\infty}} = 1 + \sum_j\frac{\omega_{\mathrm{LO},j}^2 - \omega_{\mathrm{TO},j}^2}{\omega_{\mathrm{TO},j}^2 - \omega^2 - i\gamma\omega} ,$$ जहां LO और TO अनुदैर्ध्य और अनुप्रस्थ ऑप्टिकल फोनन मोड को दर्शाते हैं, j सभी IR-सक्रिय मोड हैं, और γ ऑसिलेटर मॉडल में तापमान-निर्भर भिगोना शब्द है। εe,∞उच्च आवृत्ति ढांकता हुआ पारगम्यता है, जिसकी गणना डीएफटी गणना की जा सकती है जब आयनों को बाहरी क्षमता के रूप में माना जाता है।

इन ढांकता हुआ फ़ंक्शन से (εe,ω) गणना (उदाहरण के लिए, एबिनिट, वीएएसपी, आदि), जटिल अपवर्तक सूचकांक एमω(=एनω+ मैं श्रीमानω, एनω: अपवर्तन सूचकांक और κω: विलुप्ति सूचकांक) पाया जाता है, अर्थात्, एमω2=ईe,ω= ईe,r,ω+ मैं ईe,c,ω). निर्वात या वायु से सामान्य आपतित आदर्श सतह का सतह परावर्तन R इस प्रकार दिया गया है आर = [(एनω- 1)2+श्रीω2]/[(एनω+ 1)2+श्रीω2]. फिर वर्णक्रमीय अवशोषण गुणांक σ से पाया जाता हैph,ω= 2o कω/मेंph. विभिन्न विद्युत संस्थाओं के लिए वर्णक्रमीय अवशोषण गुणांक नीचे दी गई तालिका में सूचीबद्ध हैं।

यह भी देखें

 * ऊर्जा अंतरण
 * दूरी बदलना
 * ऊर्जा परिवर्तन|ऊर्जा परिवर्तन (ऊर्जा रूपांतरण)
 * थर्मल भौतिकी
 * ताप विज्ञान
 * थर्मल इंजीनियरिंग