1. Введение: твердотельное чудо, связывающее тепло и электричество без единой движущейся детали

Запущенные в конце лета 1977 года космические зонды-близнецы НАСА «Вояджер-1» (Voyager 1) и «Вояджер-2» (Voyager 2) уже более 47 лет летят сквозь межзвёздное пространство — на расстоянии свыше 24 миллиардов километров от Солнца — и продолжают передавать радиотелеметрию на Землю. На таком расстоянии солнечное излучение падает менее чем до одной десятитысячной от земного, делая солнечные батареи бесполезными, а любой механический тепловой двигатель с поршнями или турбинами давно вышел бы из строя из-за износа.

Сердцем, питающим «Вояджеры» полвека подряд, служит радиоизотопный термоэлектрический генератор (РИТЭГ, 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}$, управляет ПЦР-термоциклерами и охлаждает астрономические ПЗС-матрицы.

Почему градиент температуры рождает напряжение, а ток перекачивает тепло? И почему термоэлектрогенераторы ещё не вытеснили паровые турбины? Ответ кроется в фундаментальной дилемме физики конденсированного состояния — взаимоисключающей триаде между коэффициентом Зеебека $S$, электропроводностью $\sigma$ и теплопроводностью $\kappa$.

2. Объединение трёх классических термоэлектрических эффектов и взаимность Онсагера

Термоэлектричество опирается на три явления XIX века: эффект Зеебека (1821), эффект Пельтье (1834) и эффект Томсона (1851).

2.1 Эффект Зеебека: тепловое размытие Ферми — Дирака и баланс диффузии и дрейфа

Когда разомкнутый проводник или полупроводник ($\mathbf{J} = 0$) помещён в градиент температуры $\nabla T$, внутри него самопроизвольно возникает электростатическое поле $\mathbf{E} = -\nabla V$ (точнее, $-\nabla \tilde{\mu}/e$), пропорциональное $\nabla T$:

Коэффициент $S$ (в $\text{мкВ/К}$ или $\text{В/К}$) называется коэффициентом Зеебека или термо-ЭДС. Микроскопически он следует из распределения Ферми — Дирака:

1. Тепловое размытие поверхности Ферми на горячем конце: При $T > 0$ ступенька Ферми размывается в окне $\sim 2k_B T$. На горячем конце ($T_H$) концентрация высокоэнергетических носителей выше уровня Ферми $E_F$ значительно больше, чем на холодном ($T_C$).
2. Динамическое равновесие термодиффузии и электростатического дрейфа: Горячие носители диффундируют к холодному концу. В полупроводнике 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 Формула Мотта и соотношение Писаренко

Разложение Зоммерфельда кинетического уравнения Больцмана даёт для вырожденных полупроводников формулу Мотта:

где $\sigma(E) = \frac{1}{3} e^2 g(E) v^2(E) \tau(E)$ — спектральная проводимость, а $g(E)$ — плотность электронных состояний (DOS). Коэффициент Зеебека измеряет энергетическую асимметрию проводимости вблизи $E_F$. Для 3D параболической зоны при рассеянии на акустических фононах получаем соотношение Писаренко (А. Ф. Иоффе, Н. Л. Писаренко, 1940):

Высокий коэффициент Зеебека требует (1) большой эффективной массы плотности состояний $m^*$ и (2) низкой концентрации носителей $n$ ($S \propto n^{-2/3}$).

2.3 Эффект Пельтье, эффект Томсона и соотношения взаимности Кельвина — Онсагера

При пропускании постоянного тока $I$ через изотермический контакт двух разнородных материалов A и B выделяется или поглощается обратимая тепловая мощность:

Микроскопически при переходе электрона из металла ($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$ потоки связаны кинетической матрицей:

Обратимость времени на микроуровне требует соотношения взаимности Онсагера $L_{12} = L_{21}$, откуда напрямую следуют соотношения Кельвина:

Таким образом, коэффициент Зеебека $S = \frac{\Pi}{T}$ есть в точности энтропия, переносимая единичным электрическим зарядом ($S = S_{\text{носителя}} / 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$ равны:

Минимизация произведения $RK$ сводит все свойства материала к введённому академиком А. Ф. Иоффе (1949) безразмерному параметру — термоэлектрической добротности $ZT$:

где $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}$) и формулу Иоффе для максимального КПД:

3.2 Холодопроизводительность $Q_c$, максимальный перепад температур $\Delta T_{\max}$ и холодильный коэффициент COP

В охладителе Пельтье (TEC) баланс тепла на холодном спае ($T_C$) имеет вид:

Поскольку джоулево тепло растёт как $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}$ равны:

4. Фундаментальная дилемма твёрдого тела: взаимоисключающая триада $S$, $\sigma$ и $\kappa$

Почему так трудно преодолеть рубеж $ZT \approx 1$? Потому что все три параметра жёстко связаны через концентрацию носителей $n$:

1. Конфликт $S$ и $\sigma$ (компромисс Писаренко): Рост $n$ повышает $\sigma = n e \mu$, но обрушивает $S \propto m^* n^{-2/3}$. Увеличение эффективной массы $m^*$ снижает подвижность $\mu = \frac{e\tau}{m^*}$.
2. Конфликт $\sigma$ и $\kappa_e$ (закон Видемана — Франца): Одни и те же электроны переносят и заряд, и тепло:

где $L$ — число Лоренца. Поэтому оптимум достигается в сильнолегированных вырожденных полупроводниках ($n \approx 10^{19}\text{–}10^{20}\,\text{см}^{-3}$) с шириной запрещённой зоны по критерию Часмара — Страттона ($E_g \approx 6\text{–}10\,k_B T$) для подавления биполярной теплопроводности.

5. Прорыв потолка $ZT \approx 1$: зонная и фононная инженерия

  flowchart TD
    subgraph TRIAD["Классический твердотельный тупик (Триада S-σ-κ: предел 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

1. Концепция PGEC («Фононное стекло — электронный кристалл», Глен Слэк, 1995): В заполненных скуттерудитах ($\text{CoSb}_3$) и клатратах гостевые ионы внутри атомных клеток совершают ангармонические колебания («rattling»), рассеивая акустические фононы без рассеяния электронов.
2. Зонная инженерия: Низкоразмерное квантовое ограничение (Хикс и Дрессельхаус, 1993), резонансные уровни ($\text{Tl}:\text{PbTe}$, Хереманс и др., 2008) и конвергенция зон ($m^* = N_v^{2/3} m^*_b$, выравнивание долин $L$ и $\Sigma$ до $N_v = 16$) повышают $S$ без потери подвижности $\mu$.
3. Всемасштабная иерархическая фононная инженерия и $\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{Вт/(м}\cdot\text{К)}$).

6. Практические применения и будущее: от дальнего космоса до спин-калоритроники

1. Космические РИТЭГи (RTG и MMRTG): От $\text{SiGe}$-термопар «Вояджеров» до MMRTG ($\text{PbTe/TAGS}$) марсоходов Curiosity и Perseverance и ядерного октокоптера Dragonfly для полётов на Титане.
2. Прецизионное терморегулирование (TEC): Стабилизация температуры оптоволоконных DWDM-лазеров с точностью $\pm 0{,}01\,\text{K}$, ПЦР-амплификаторы ДНК и безвибрационное охлаждение ИК- и ПЗС-детекторов.
3. Рекуперация сбросного тепла, носимая термоэлектрика и спин-калоритроника: Автономные датчики IIoT, пластичные неорганические полупроводники ($\text{Ag}_2\text{S}$) и проводящие полимеры (PEDOT:PSS) для генерации энергии от тепла тела, а также спиновый эффект Зеебека (SSE) и аномальный эффект Нернста (ANE, $\mathbf{E}_{\text{ANE}} \propto \mathbf{M} \times \nabla T$) в топологических магнитных полуметаллах Вейля.