1. مقدمة: معجزة الحالة الصلبة التي تربط الحرارة بالكهرباء دون أي أجزاء متحركة
منذ إطلاقتهما في أواخر صيف عام 1977، يواصل المسباران الفضائيان التوأمان التابعان لوكالة ناسا «فوياجر 1» (Voyager 1) و«فوياجر 2» (Voyager 2) رحلتهما عبر الفضاء بين النجمي — على بعد أكثر من 24 مليار كيلومتر من الشمس — لأكثر من 47 عاماً مع الاستمرار في إرسال البيانات إلى الأرض. عند هذه المسافة الهائلة، تنخفض شدة الإشعاع الشمسي إلى أقل من واحد من عشرة آلاف من قيمتها عند الأرض، مما يجعل الألواح الشمسية عديمة الجدوى؛ كما أن أي محرك حراري ميكانيكي يعتمد على المكابس أو التوربينات كان سيتعطل منذ عقود بفعل التآكل.
إن القلب النابض الذي يمد مسباري فوياجر بالطاقة طوال نصف قرن هو المولد الكهرحراري بالنظائر المشعة (RTG): مصدر طاقة من الحالة الصلبة يحول فرق درجات الحرارة بين حرارة اضمحلال ألفا للبلوتونيوم-238 ($^{238}\text{Pu}$، $T_H \approx 1273\,\text{K}$) وبرودة الفضاء السحيق ($T_C \approx 573\,\text{K}$) مباشرة إلى تيار كهربائي مستمر باستخدام مزدوجات أشباه موصلات السيليكون-الجرمانيوم ($\text{SiGe}$). وعلى الأرض، يعمل القانون الفيزيائي ذاته في الاتجاه المعاكس داخل عناصر بلتييه (المبردات الكهرحرارية TEC): فبدون قطرة واحدة من غاز التبريد أو ضاغط ميكانيكي، يؤدي مرور تيار كهربائي إلى تجميد أحد الوجهين الخزفيين، وتثبيت الطول الموجي لليزر اتصالات الألياف البصرية DWDM بدقة $\pm 0.01\,\text{K}$، وتشغيل أجهزة تفاعل البوليميراز المتسلسل (PCR)، وتبريد مستشعرات التصوير الفلكي CCD.
لماذا يولد التدرج الحراري فرق جهد كهربائي، ولماذا يضخ التيار الكهربائي الحرارة؟ ولماذا لم تحل المولدات الكهرحرارية محل التوربينات البخارية في محطات الطاقة؟ تكمن الإجابة في واحدة من أعقد معضلات فيزياء المادة المكثفة: الثلاثية المتعارضة (المقايضة الأساسية) بين معامل سيبيك $S$، والموصلية الكهربائية $\sigma$، والموصلية الحرارية $\kappa$.
2. توحيد التأثيرات الكهرحرارية الكلاسيكية الثلاثة وعلاقات أونساغر التبادلية
ترتكز الكهرحرارية على ثلاث ظواهر اكتُشفت في القرن التاسع عشر: تأثير سيبيك (1821)، وتأثير بلتييه (1834)، وتأثير طومسون (1851).
2.1 تأثير سيبيك: التوسع الحراري لتوزيع فيرمي-ديراك وتوازن الانتشار والانجراف
عندما يخضع موصل أو شبه موصل في دائرة مفتوحة ($\mathbf{J} = 0$) لتدرج حراري مكاني $\nabla T$، ينشأ تلقائياً مجال كهروستاتيكي داخلي $\mathbf{E} = -\nabla V$ (أو بدقة $-\nabla \tilde{\mu}/e$) يتناسب طردياً مع $\nabla T$:
$$\nabla V = -S \nabla T \quad \Longleftrightarrow \quad S = -\frac{\Delta V}{\Delta T}$$يُعرف الثابت $S$ (بوحدة $\mu\text{V/K}$ أو $\text{V/K}$) باسم معامل سيبيك (Seebeck coefficient) أو القوة الدافعة الكهرحرارية. ويُفسَّر مجهرياً عبر دالة توزيع فيرمي-ديراك:
$$f(E, T) = \frac{1}{\exp\left(\frac{E - \mu}{k_B T}\right) + 1}$$- التوسع الحراري لسطح فيرمي عند الطرف الساخن: عند درجة حرارة $T > 0$، تتسع حافة فيرمي ضمن نافذة طاقية $\sim 2k_B T$. ويحتوي الطرف الساخن ($T_H$) على كثافة أعلى بكثير من حاملات الشحنة عالية الطاقة فوق مستوى فيرمي $E_F$ مقارنة بالطرف البارد ($T_C$).
- التوازن الديناميكي بين الانتشار الحراري والانجراف الكهروستاتيكي: تنتشر الحاملات الساخنة نحو الطرف البارد. ففي شبه الموصل من النوع n تتراكم الإلكترونات السالبة ($q = -e$) عند الطرف البارد ($S_n < 0$)؛ وفي شبه الموصل من النوع p تتراكم الفجوات الموجبة ($q = +e$) عند الطرف البارد ($S_p > 0$). يولّد هذا الفصل الشحني مجالاً كهربائياً داخلياً يدفع تيار انجراف معاكساً حتى يتعادل تيار الانتشار الحراري مع تيار الانجراف الكهروستاتيكي تماماً ($\mathbf{J}_{\text{diff}} + \mathbf{J}_{\text{drift}} = 0$). وعند توصيل فرعي p وn حرارياً على التوازي وكهربائياً على التوالي (مزدوجة $\pi$)، يتجمع جهد الدائرة المفتوحة: $V_{\text{oc}} = (S_p - S_n)(T_H - T_C)$.
2.2 صيغة موت وعلاقة بيسارينكو
باستخدام مفكوك سومرفيلد لمعادلة بولتزمان للنقل، يخضع معامل سيبيك في أشباه الموصلات المنحلة لصيغة موت (Mott Formula):
$$S = -\frac{\pi^2 k_B^2 T}{3e} \left. \frac{d \ln \sigma(E)}{dE} \right|_{E = E_F}$$حيث $\sigma(E) = \frac{1}{3} e^2 g(E) v^2(E) \tau(E)$ هي الموصلية الطيفية و$g(E)$ هي كثافة الحالات الإلكترونية (DOS). يقيس معامل سيبيك إذن عدم التناظر الطاقي للتوصيل حول مستوى فيرمي $E_F$. ولحزمة طاقة مكافئة ثلاثية الأبعاد مع تشتت الفونونات الصوتية، تؤول صيغة موت إلى علاقة بيسارينكو (Pisarenko Relation):
$$S = \frac{8 \pi^2 k_B^2}{3 e h^2} m^* T \left( \frac{\pi}{3n} \right)^{2/3}$$يتطلب الحصول على معامل سيبيك كبير (1) كتلة فعالة ثقيلة لكثافة الحالات $m^*$ و(2) تركيز حاملات شحنة منخفضاً $n$ ($S \propto n^{-2/3}$).
2.3 تأثير بلتييه، وتأثير طومسون، وعلاقات كلفن-أونساغر التبادلية
عندما يمر تيار مستمر $I$ عبر وصلة متساوية الحرارة بين مادتين مختلفتين A وB، تُمتص الحرارة أو تنبعث بشكل عكوس بمعدل:
$$\dot{Q}_P = (\Pi_A - \Pi_B) I = \Pi_{AB} I$$مجهرياً، عندما ينتقل إلكترون من المعدن ($E \approx E_F$) إلى حزمة التوصيل في شبه موصل من النوع n ($E_{\text{trans}} > E_C$) عند الوصلة الباردة، يتعين عليه تسلق حاجز طاقي $\Delta E = E_{\text{trans}} - E_F > 0$ عبر امتصاص الفونونات من الشبكة البلورية مما يبرّد الوصلة؛ وعند الوصلة الساخنة يهبط عائداً إلى $E_F$ مطلقاً فائض طاقته كفونونات. وبالمثل يصف تأثير طومسون (1851) التبادل الحراري العكوس $\dot{q}_{\text{total}} = \frac{\mathbf{J}^2}{\sigma} - \mathcal{K}_T (\mathbf{J} \cdot \nabla T)$ في موصل متجانس تحت تدرج حراري $\nabla T$.
وفي الديناميكا الحرارية الخطية لغير الاتزان التي صاغها لارس أونساغر (1931، جائزة نوبل 1968)، يرتبط تيار الشحنة $\mathbf{J}_e$ وتدفق الحرارة $\mathbf{J}_Q$ بالمصفوفة الحركية:
$$\begin{pmatrix} \mathbf{J}_e \\ \mathbf{J}_Q \end{pmatrix} = \begin{pmatrix} L_{11} & L_{12} \\ L_{21} & L_{22} \end{pmatrix} \begin{pmatrix} -\frac{1}{e}\nabla \tilde{\mu} \\ \nabla \left(\frac{1}{T}\right) \end{pmatrix}$$يفرض تناظر انعكاس الزمن المجهري علاقة أونساغر التبادلية $L_{12} = L_{21}$، ومنها تنتج مباشرة علاقتا كلفن:
$$\Pi = S \cdot T, \qquad \mathcal{K}_T = T \frac{dS}{dT}$$وهذا يكشف المعنى الفيزيائي العميق لمعامل سيبيك: إن $S = \frac{\Pi}{T}$ هو بالضبط الإنتروبيا المنقولة لكل وحدة شحنة كهربائية ($S = S_{\text{carrier}} / q$)!
flowchart TD
subgraph TEG["1. نمط توليد الطاقة لسيبيك (المولد الكهرحراري: TEG)"]
TH1["المستودع الساخن T_H (حرارة داخلة Q_H من العادم / النظير المشع)"] --> N1["فرع شبه موصل n: انتشار حراري للإلكترونات e⁻ نحو البارد"]
TH1 --> P1["فرع شبه موصل p: انتشار حراري للفجوات h⁺ نحو البارد"]
N1 --> TC1["المشتت البارد T_C (طرد الحرارة Q_C)"]
P1 --> TC1
TC1 --> LOAD["حمل خارجي R_L: القوة الدافعة V = (S_p - S_n)(T_H - T_C) تولد القدرة P"]
end
subgraph TEC["2. نمط تبريد بلتييه (المبرد الكهرحراري: TEC)"]
DC["مصدر تيار مستمر خارجي يمرر التيار I"] --> N2["فرع n: الإلكترونات e⁻ تتسلق حاجز السطح البيني (E_C - E_F)"]
DC --> P2["فرع p: الفجوات h⁺ تدخل حزمة التكافؤ E_V"]
N2 --> COLD["الوصلة الباردة T_C: امتصاص فونونات الشبكة Q_c = (S_p - S_n)T_C I"]
P2 --> COLD
COLD --> HOT["الوصلة الساخنة T_H: تفريغ الطاقة الفائضة للمشتت الحراري Q_H"]
end
3. الاشتقاق الديناميكي الحراري لمعامل الجدارة $ZT$ والكفاءة القصوى وقدرة التبريد
3.1 الكفاءة القصوى $\eta_{\max}$ للمولد الكهرحراري (TEG)
لمزدوجة كهرحرارية تعمل بين $T_H$ و$T_C$ ($\Delta T = T_H - T_C$) بمقاومة داخلية $R$ ومواصلة حرارية $K$ ومعامل سيبيك $S = S_p - S_n$، تُعطى الحرارة المسحوبة من المصدر الساخن $Q_H$ والقدرة الكهربائية $P$ المسلمة للحمل $R_L$ بالعلاقتين:
$$Q_H = S T_H I + K (T_H - T_C) - \frac{1}{2} I^2 R, \qquad P = I^2 R_L = S(T_H - T_C)I - I^2 R$$وبتحسين النسبة الهندسية لتقليل الجداء $RK$، تتجمع خصائص المادة في معامل الجدارة اللابعدي $ZT$ الذي قدمه أبرام يوفي (1949):
$$ZT = \frac{S^2 \sigma T}{\kappa} = \frac{S^2 \sigma T}{\kappa_e + \kappa_L}$$حيث $PF = S^2 \sigma$ هو معامل القدرة و$\kappa = \kappa_e + \kappa_L$ هو مجموع الموصلية الحرارية الإلكترونية $\kappa_e$ والشبكية (الفونونية) $\kappa_L$. وبتعظيم الكفاءة $\eta(m) = P / Q_H$ بالنسبة إلى $m = R_L / R$ نجد $m_{\text{opt}} = \sqrt{1 + Z\bar{T}}$ ($\bar{T} = \frac{T_H + T_C}{2}$) ومعادلة الكفاءة القصوى:
$$\eta_{\max} = \underbrace{\frac{T_H - T_C}{T_H}}_{\text{كفاءة كارنو }\eta_C} \cdot \underbrace{\frac{\sqrt{1 + Z\bar{T}} - 1}{\sqrt{1 + Z\bar{T}} + \frac{T_C}{T_H}}}_{\text{معامل المادة}}$$3.2 قدرة التبريد الصافية $Q_c$ وأقصى فرق حراري $\Delta T_{\max}$ ومعامل الأداء COP
في مبرد بلتييه (TEC)، يُعطى التوازن الحراري عند الوجه البارد ($T_C$) بالمعادلة:
$$Q_c = \underbrace{(S_p - S_n) T_C I}_{\text{(1) تبريد بلتييه}} - \underbrace{\frac{1}{2} I^2 R}_{\text{(2) التسخين العكسي لجول}} - \underbrace{K (T_H - T_C)}_{\text{(3) التوصيل الحراري العكسي لفورييه}}$$وبما أن تسخين جول ينمو تربيعياً ($\propto I^2$)، يوجد تيار أمثل $I_{\max} = \frac{S T_C}{R}$ ($\frac{\partial Q_c}{\partial I} = 0$). وعند انعدام الحمل الحراري الخارجي ($Q_c = 0$)، نحصل على أقصى فرق في درجات الحرارة $\Delta T_{\max}$ وأقصى معامل أداء $\text{COP}_{\max}$:
$$\Delta T_{\max} = \frac{1}{2} Z T_C^2, \qquad \text{COP}_{\max} = \frac{T_C}{T_H - T_C} \cdot \frac{\sqrt{1 + Z\bar{T}} - \frac{T_H}{T_C}}{\sqrt{1 + Z\bar{T}} + 1}$$4. معضلة الحالة الصلبة الأساسية: الثلاثية المتعارضة لـ $S$ و$\sigma$ و$\kappa$
لماذا يصعب تجاوز حاجز $ZT \approx 1$؟ لأن المعاملات الثلاثة تتصارع عبر تركيز حاملات الشحنة $n$:
- تعارض $S$ مع $\sigma$ (مقايضة بيسارينكو): زيادة $n$ ترفع $\sigma = n e \mu$، لكنها تخفض $S \propto m^* n^{-2/3}$ بشدة. أما زيادة الكتلة الفعالة $m^*$ فتخفض الحركية $\mu = \frac{e\tau}{m^*}$.
- تعارض $\sigma$ مع $\kappa_e$ (قانون فيدمان-فرانز): الإلكترونات نفسها تنقل الشحنة والحرارة معاً وفق:
حيث $L$ هو عدد لورنتز. ولذلك تقع المنطقة المثلى في أشباه الموصلات المنحلة عالية الإشابة ($n \approx 10^{19}\text{–}10^{20}\,\text{cm}^{-3}$) مع فجوة طاقة تحقق معيار تشاسمار-ستراتون ($E_g \approx 6\text{–}10\,k_B T$) لكبح التوصيل ثنائي القطب.
| النطاق الحراري | نظام المواد | الحرارة المثلى $T_{\text{opt}}$ | فجوة الطاقة $E_g$ | ذروة $ZT$ النموذجية | البنية البلورية، الفيزياء الأساسية والتطبيقات |
|---|---|---|---|---|---|
| التبريد الفائق إلى الغرفة ($< 450\,\text{K}$) | سبائك $\text{Bi}_2\text{Te}_3$ | $250\text{–}400\,\text{K}$ | $\sim 0.15\,\text{eV}$ | $0.8\text{–}1.4$ | بنية طبقية معينية، تعدد أودية $N_v=6$. المعيار العالمي لمبردات بلتييه (TEC). |
| الحرارة المتوسطة I ($500\text{–}850\,\text{K}$) | $\text{PbTe}$ وTAGS-85 | $600\text{–}800\,\text{K}$ | $\sim 0.30\,\text{eV}$ | $1.2\text{–}2.2$ | بنية مكعبة. تقارب الحزم L–$\Sigma$ ($N_v = 4 \to 16$)، إشابة رنينية بـ $\text{Tl}$ ورواسب نانوية $\text{SrTe}$. عربات المريخ MMRTG. |
| الحرارة المتوسطة II ($550\text{–}900\,\text{K}$) | السكوتيروديت المملوء ($\text{Ba}_x\text{La}_y\text{Yb}_z\text{Co}_4\text{Sb}_{12}$) | $700\text{–}850\,\text{K}$ | $\sim 0.3\,\text{eV}$ | $1.2\text{–}1.8$ | بنية أقفاص ($\text{CoSb}_3$). الأيونات الضيفة «تقعقع» لا توافقياً (آلية PGEC) فتشتت الفونونات. عوادم السيارات. |
| المتوسطة–العالية ($650\text{–}1000\,\text{K}$) | بلورة $\text{SnSe}$ الأحادية | $750\text{–}923\,\text{K}$ | $\sim 0.86\,\text{eV}$ | $2.2\text{–}2.6+$ | بنية طبقية مطوية ($Pnma \to Cmcm$). لاتوافقية شبكية ذاتية عملاقة ($\gamma \approx 4.1$) تعطي $\kappa_L \approx 0.23\,\text{W/(m}\cdot\text{K)}$. |
| الحرارة العالية ($700\text{–}1100\,\text{K}$) | سبائك Half-Heusler ($\text{ZrNiSn}$, $\text{NbFeSb}$) | $800\text{–}1050\,\text{K}$ | $\sim 0.5\,\text{eV}$ | $1.0\text{–}1.5$ | متانة ميكانيكية فائقة، غير سامة، معامل قدرة ضخم ($PF > 50\,\mu\text{W cm}^{-1}\text{K}^{-2}$). أفران الصلب والصناعة. |
| الحرارة فوق العالية ($900\text{–}1300\,\text{K}$) | سبائك $\text{SiGe}$ ($\text{Si}_{0.8}\text{Ge}_{0.2}$) | $1000\text{–}1273\,\text{K}$ | $\sim 0.8\,\text{eV}$ | $0.8\text{–}1.3$ | شبكة تساهمية ماسية مستقرة عند $1000\,^\circ\text{C}$. تغذي مولدات RTG لمسباري فوياجر ونيو هورايزونز. |
5. كسر سقف $ZT \approx 1$: هندسة الحزم وهندسة الفونونات الحديثة
flowchart TD
subgraph TRIAD["عنق الزجاجة الكلاسيكي (الثلاثية المتعارضة: سقف ZT ≈ 1)"]
N_UP["زيادة تركيز الحاملات n"] --> SIG_UP["ارتفاع الموصلية الكهربائية σ = neμ (مفيد)"]
N_UP --> S_DOWN["هبوط معامل سيبيك S ∝ n⁻²/³ وفق بيسارينكو (ضار)"]
SIG_UP --> KE_UP["ارتفاع الموصلية الحرارية الإلكترونية κ_e = LσT وفق فيدمان-فرانز (ضار)"]
end
subgraph MODERN["استراتيجيات فك الارتباط الحديثة (الوصول إلى ZT > 2.0–2.6)"]
BAND["1. هندسة بنية الحزم (تعظيم معامل القدرة S²σ)"] --> QD["الحصر الكمي (2D/1D/0D) والإشابة الرنينية (Tl:PbTe): زيادة انحدار d ln g(E)/dE"]
BAND --> CONV["تقارب الحزم: محاذاة أودية L وΣ لمضاعفة N_v دون فقدان الحركية μ"]
PHONON["2. هندسة الفونونات (خفض κ_L إلى الحد الزجاجي)"] --> PGEC["مفهوم PGEC والأيونات المقعقعة (السكوتيروديت المملوء): تشتيت رنيني للفونونات"]
PHONON --> ALLSCALE["البنية الهرمية شاملة المقاييس: عيوب نقطية + رواسب نانوية متوافقة + حدود حبيبية"]
PHONON --> ANHARM["اللاتوافقية الذاتية العملاقة (بلورة SnSe الأحادية): معامل غرونايزن γ ≈ 4.1"]
end
TRIAD ==> MODERN
- مفهوم PGEC («زجاج الفونونات – بلورة الإلكترونات»، غلين سلاك، 1995): في السكوتيروديت المملوء ($\text{CoSb}_3$) والكلathrates، تهتز الأيونات الضيفة داخل الأقفاص الذرية بشكل لا توافقي («rattling») فتشتت الفونونات الصوتية دون إعاقة الإلكترونات.
- هندسة بنية الحزم: يتيح الحصر الكمي منخفض الأبعاد (هيكس ودريسلهوس، 1993)، والإشابة بالمستويات الرنينية ($\text{Tl}:\text{PbTe}$، هيريمانز وآخرون، 2008)، وتقارب الحزم ($m^* = N_v^{2/3} m^*_b$، محاذاة حزمتي $L$ و$\Sigma$ لرفع التعدد إلى $N_v = 16$) رفع $S$ دون التضحية بالحركية $\mu$.
- الهندسة الفونونية الهرمية شاملة المقاييس وبلورة $\text{SnSe}$: عبر دمج العيوب النقطية الذرية والرواسب النانوية المتوافقة ($\text{SrTe}$) والحدود الحبيبية الميزوسكوبية (كاناتزيديس وآخرون، 2012) جرى تشتيت كامل طيف الفونونات ($ZT \approx 2.2$). وفي عام 2014، حققت البلورة الأحادية $\text{SnSe}$ رقماً قياسياً بلغ $ZT = 2.6 \pm 0.3$ عند $923\,\text{K}$ بفضل اللاتوافقية الشبكية الذاتية العملاقة ($\gamma \approx 4.1$، $\kappa_L \approx 0.23\,\text{W/(m}\cdot\text{K)}$).
6. التطبيقات العملية وآفاق المستقبل: من مسابير الفضاء إلى الكهرحرارية الدورانية (Spin Caloritronics)
- استكشاف الفضاء العميق (RTG وMMRTG): من مزدوجات $\text{SiGe}$ في فوياجر 1 و2 ونيو هورايزونز إلى مولدات MMRTG ($\text{PbTe/TAGS}$) في عربتي المريخ كيوريوسيتي وبيرسيفيرانس ومروحية دراغون فلاي النووية المتجهة إلى تيتان.
- الإدارة الحرارية الدقيقة (TEC): تثبيت حرارة ليزر اتصالات DWDM في حدود $\pm 0.01\,\text{K}$، وأجهزة PCR الطبية، والتبريد الخالي من الاهتزاز لمستشعرات الفلك والأشعة تحت الحمراء.
- استرجاع الحرارة المهدرة، والكهرحرارية القابلة للارتداء، والكهرحرارية الدورانية: مستشعرات إنترنت الأشياء الصناعي (IIoT) ذاتية التغذية، وأشباه الموصلات غير العضوية المطيلة ($\text{Ag}_2\text{S}$) والبوليمرات الموصلة (PEDOT:PSS) لحصاد حرارة الجسم، إضافة إلى تأثير سيبيك المغزلي (SSE) وتأثير نيرنست الشاذ (ANE، $\mathbf{E}_{\text{ANE}} \propto \mathbf{M} \times \nabla T$) في أشباه معادن فايل المغناطيسية الطوبولوجية.