1. Pendahuluan: Keajaiban Zat Padat yang Menghubungkan Panas dan Listrik Tanpa Bagian Bergerak

Diluncurkan pada akhir musim panas 1977, wahana antariksa kembar NASA Voyager 1 dan Voyager 2 telah melesat menembus ruang antarbintang—lebih dari 24 miliar kilometer dari Matahari—selama lebih dari 47 tahun sambil terus mengirimkan telemetri ke Bumi. Pada jarak sejauh itu, radiasi matahari turun hingga kurang dari seperseribu intensitas di Bumi sehingga panel surya tidak berguna, sementara mesin kalor mekanis dengan piston atau turbin pasti sudah lama rusak akibat keausan.

Jantung yang menghidupi Voyager selama setengah abad adalah Generator Termoelektrik Radioisotop (RTG): sumber daya zat padat yang mengubah perbedaan suhu antara panas peluruhan alfa Plutonium-238 ($^{238}\text{Pu}$, $T_H \approx 1273\,\text{K}$) dan dinginnya antariksa ($T_C \approx 573\,\text{K}$) secara langsung menjadi listrik arus searah (DC) menggunakan pasangan semikonduktor silikon-germanium ($\text{SiGe}$). Di Bumi, hukum fisika yang sama ketika dijalankan terbalik menggerakkan elemen Peltier (Pendingin Termoelektrik, TEC): tanpa setetes pun refrigeran atau kompresor, arus listrik membekukan satu sisi keramik, menstabilkan laser telekomunikasi serat optik DWDM hingga $\pm 0{,}01\,\text{K}$, menggerakkan mesin PCR, dan mendinginkan sensor CCD astronomi.

Mengapa gradien suhu menghasilkan tegangan listrik, dan arus listrik memompa panas? Mengapa generator termoelektrik belum menggantikan turbin uap di pembangkit listrik? Jawabannya terletak pada dilema paling mendasar dalam fisika benda terkondensasi: triad saling bertentangan (trade-off fundamental) antara koefisien Seebeck $S$, konduktivitas listrik $\sigma$, dan konduktivitas termal $\kappa$.

2. Menyatukan Tiga Efek Termoelektrik Klasik dan Resiprositas Onsager

Termoelektrik berdiri di atas tiga fenomena abad ke-19: efek Seebeck (1821), efek Peltier (1834), dan efek Thomson (1851).

2.1 Efek Seebeck: Pelebaran Termal Fermi-Dirac dan Kesetimbangan Difusi-Hanyutan

Ketika konduktor atau semikonduktor rangkaian terbuka ($\mathbf{J} = 0$) dikenai gradien suhu $\nabla T$, medan elektrostatik internal $\mathbf{E} = -\nabla V$ (tepatnya $-\nabla \tilde{\mu}/e$) muncul secara spontan sebanding dengan $\nabla T$:

Konstanta $S$ (dalam $\mu\text{V/K}$ atau $\text{V/K}$) adalah koefisien Seebeck atau daya termoelektrik. Secara mikroskopis, efek ini berasal dari fungsi distribusi Fermi-Dirac:

1. Pelebaran Termal Permukaan Fermi di Ujung Panas: Pada $T > 0$, tepi Fermi melebar pada jendela energi $\sim 2k_B T$. Ujung panas ($T_H$) memiliki populasi pembawa muatan berenergi tinggi di atas $E_F$ jauh lebih banyak daripada ujung dingin ($T_C$).
2. Kesetimbangan Dinamis antara Difusi Termal dan Hanyutan Elektrostatik: Pembawa muatan panas berdifusi menuju ujung dingin. Dalam semikonduktor tipe-n, elektron ($q = -e$) menumpuk di ujung dingin ($S_n < 0$); dalam semikonduktor tipe-p, lubang/hole ($q = +e$) menumpuk di ujung dingin ($S_p > 0$). Medan listrik internal yang terbentuk mendorong arus hanyutan (drift) berlawanan hingga arus difusi termal dan arus hanyutan elektrostatik saling meniadakan secara tepat ($\mathbf{J}_{\text{diff}} + \mathbf{J}_{\text{drift}} = 0$). Dengan menghubungkan kaki tipe-p dan tipe-n secara termal paralel dan secara elektrik seri (pasangan-$\pi$), tegangan rangkaian terbuka saling menjumlah: $V_{\text{oc}} = (S_p - S_n)(T_H - T_C)$.

2.2 Rumus Mott dan Relasi Pisarenko

Melalui ekspansi Sommerfeld dari persamaan transpor Boltzmann, koefisien Seebeck semikonduktor terdegenerasi memenuhi rumus Mott:

di mana $\sigma(E) = \frac{1}{3} e^2 g(E) v^2(E) \tau(E)$ adalah konduktivitas spektral dan $g(E)$ adalah rapat keadaan elektronik (DOS). Koefisien Seebeck mengukur asimetri energi konduksi di sekitar $E_F$. Untuk pita parabolik 3D isotropik dengan hamburan fonon akustik, diperoleh relasi Pisarenko:

Mencapai koefisien Seebeck yang besar membutuhkan (1) massa efektif rapat keadaan $m^*$ yang berat dan (2) konsentrasi pembawa muatan $n$ yang rendah ($S \propto n^{-2/3}$).

2.3 Efek Peltier, Efek Thomson, dan Relasi Resiprokal Kelvin-Onsager

Ketika arus searah $I$ melewati sambungan isotermal antara dua material A dan B, kalor diserap atau dilepaskan secara reversibel dengan laju:

Secara mikroskopis, ketika elektron melintas dari logam ($E \approx E_F$) ke pita konduksi semikonduktor tipe-n ($E_{\text{trans}} > E_C$) di sambungan dingin, ia harus mendaki penghalang energi $\Delta E = E_{\text{trans}} - E_F > 0$ dengan menyerap fonon dari kisi kristal sehingga mendinginkan sambungan; di sambungan panas, elektron turun kembali ke $E_F$ dan membuang kelebihan energinya sebagai fonon. Demikian pula, efek Thomson (1851) menjelaskan pertukaran kalor reversibel $\dot{q}_{\text{total}} = \frac{\mathbf{J}^2}{\sigma} - \mathcal{K}_T (\mathbf{J} \cdot \nabla T)$ sepanjang konduktor homogen dengan gradien $\nabla T$.

Dalam termodinamika non-kesetimbangan linier Lars Onsager (1931, Nobel Kimia 1968), rapat arus listrik $\mathbf{J}_e$ dan fluks kalor $\mathbf{J}_Q$ memenuhi:

Invariansi pembalikan waktu mikroskopis menuntut resiprositas Onsager $L_{12} = L_{21}$, yang langsung menghasilkan relasi Kelvin:

Ini mengungkap bahwa koefisien Seebeck $S = \frac{\Pi}{T}$ adalah tepat entropi yang diangkut per satuan muatan listrik ($S = S_{\text{pembawa}} / q$)!

  flowchart TD
    subgraph TEG["1. Mode Pembangkitan Seebeck (Generator Termoelektrik: TEG)"]
        TH1["Reservoir Panas T_H (Masukan kalor Q_H dari gas buang / radioisotop)"] --> N1["Kaki Semikonduktor n: Elektron e⁻ berdifusi termal ke sisi dingin"]
        TH1 --> P1["Kaki Semikonduktor p: Hole h⁺ berdifusi termal ke sisi dingin"]
        N1 --> TC1["Pendingin T_C (Pembuangan kalor Q_C)"]
        P1 --> TC1
        TC1 --> LOAD["Beban Eksternal R_L: GGL V = (S_p - S_n)(T_H - T_C) menghasilkan daya P"]
    end

    subgraph TEC["2. Mode Pendinginan Peltier (Pendingin Termoelektrik: TEC)"]
        DC["Catu Daya DC Eksternal mengalirkan arus I"] --> N2["Kaki n: Elektron e⁻ mendaki penghalang antarmuka (E_C - E_F)"]
        DC --> P2["Kaki p: Hole h⁺ masuk ke pita valensi E_V"]
        N2 --> COLD["Sambungan Dingin T_C: Fonon kisi diserap Q_c = (S_p - S_n)T_C I"]
        P2 --> COLD
        COLD --> HOT["Sambungan Panas T_H: Kelebihan energi dibuang ke heatsink Q_H"]
    end

3. Penurunan Termodinamika Angka Mérito $ZT$, Efisiensi Maksimum, dan Kapasitas Pendinginan

3.1 Efisiensi Konversi Maksimum $\eta_{\max}$ pada Generator Termoelektrik (TEG)

Untuk pasangan termoelektrik antara $T_H$ dan $T_C$ ($\Delta T = T_H - T_C$) dengan hambatan internal $R$, konduktansi termal $K$, dan koefisien Seebeck $S = S_p - S_n$, laju kalor yang diserap di sisi panas $Q_H$ dan daya listrik $P$ pada beban $R_L$ adalah:

Dengan meminimalkan hasil kali geometri $RK$, seluruh sifat material menyatu ke dalam angka mérito tak berdimensi (Figure of Merit) $ZT$ yang dirumuskan oleh Abram F. Ioffe (1949):

di mana $PF = S^2 \sigma$ adalah faktor daya dan $\kappa = \kappa_e + \kappa_L$ adalah jumlah konduktivitas termal elektronik $\kappa_e$ dan kisi (fonon) $\kappa_L$. Memaksimalkan $\eta(m) = P / Q_H$ terhadap rasio hambatan $m = R_L / R$ memberikan $m_{\text{opt}} = \sqrt{1 + Z\bar{T}}$ ($\bar{T} = \frac{T_H + T_C}{2}$) serta efisiensi konversi maksimum:

3.2 Daya Pendinginan $Q_c$, Perbedaan Suhu Maksimum $\Delta T_{\max}$, dan COP pada Mode Peltier

Pada pendingin Peltier (TEC), neraca kalor di sisi dingin ($T_C$) adalah:

Karena pemanasan Joule tumbuh secara kuadratik ($\propto I^2$), terdapat arus optimal $I_{\max} = \frac{S T_C}{R}$ ($\frac{\partial Q_c}{\partial I} = 0$). Pada beban kalor nol ($Q_c = 0$), penurunan suhu maksimum $\Delta T_{\max}$ dan koefisien kinerja maksimum $\text{COP}_{\max}$ adalah:

4. Dilema Fundamental Zat Padat: Triad Saling Bertentangan $S$, $\sigma$, dan $\kappa$

Mengapa begitu sulit melampaui $ZT \approx 1$? Karena ketiga parameter saling berkonflik melalui konsentrasi pembawa muatan $n$:

1. Konflik $S$ vs $\sigma$ (Kompromi Pisarenko): Menaikkan $n$ meningkatkan $\sigma = n e \mu$, tetapi menjatuhkan $S \propto m^* n^{-2/3}$. Memperbesar massa efektif $m^*$ menurunkan mobilitas $\mu = \frac{e\tau}{m^*}$.
2. Konflik $\sigma$ vs $\kappa_e$ (Hukum Wiedemann-Franz): Elektron yang sama mengangkut muatan sekaligus kalor menurut:

di mana $L$ adalah bilangan Lorenz. Oleh karena itu, titik optimal berada pada semikonduktor terdegenerasi terdoping berat ($n \approx 10^{19}\text{–}10^{20}\,\text{cm}^{-3}$) dengan celah pita sesuai kriteria Chasmar-Stratton ($E_g \approx 6\text{–}10\,k_B T$) guna menekan konduksi bipolar.

5. Menembus Batas $ZT \approx 1$: Rekayasa Pita dan Rekayasa Fonon Modern

  flowchart TD
    subgraph TRIAD["Kemacetan Klasik (Triad Saling Bertentangan: Batas ZT ≈ 1)"]
        N_UP["Naikkan Konsentrasi Pembawa n"] --> SIG_UP["Konduktivitas Listrik σ = neμ Naik (Baik)"]
        N_UP --> S_DOWN["Koefisien Seebeck S ∝ n⁻²/³ Turun via Pisarenko (Buruk)"]
        SIG_UP --> KE_UP["Konduktivitas Termal Elektronik κ_e = LσT Naik via Wiedemann-Franz (Buruk)"]
    end

    subgraph MODERN["Strategi Dekopling Modern (Menembus ZT > 2,0–2,6)"]
        BAND["1. Rekayasa Struktur Pita (Memaksimalkan S²σ)"] --> QD["Kurungan Kuantum (2D/1D/0D) & Doping Resonansi (Tl:PbTe): Mempercuram d ln g(E)/dE"]
        BAND --> CONV["Konvergensi Pita: Menyelaraskan lembah L dan Σ untuk melipatgandakan N_v tanpa kehilangan mobilitas μ"]
        PHONON["2. Rekayasa Fonon (Menekan κ_L ke Batas Kaca Amorf)"] --> PGEC["PGEC & Ion 'Rattler' (Skutterudite Terisi): Hamburan Fonon Resonansi"]
        PHONON --> ALLSCALE["Arsitektur Hierarkis Semua-Skala: Cacat Titik + Nanopresipitat Endotaksial + Batas Butir"]
        PHONON --> ANHARM["Anharmonisitas Intrinsik Raksasa (Kristal Tunggal SnSe): Grüneisen γ ≈ 4,1 Hamburan Umklapp"]
    end

    TRIAD ==> MODERN

1. Konsep PGEC (“Phonon-Glass Electron-Crystal”, Glen Slack, 1995): Dalam skutterudite terisi ($\text{CoSb}_3$) dan klatrat, ion tamu di dalam sangkar atom bergetar secara anharmonik (“rattling”) dan menghamburkan fonon akustik tanpa menghambat elektron.
2. Rekayasa Pita: Kurungan kuantum berdimensi rendah (Hicks & Dresselhaus, 1993), doping tingkat resonansi ($\text{Tl}:\text{PbTe}$, Heremans et al., 2008), dan konvergensi pita ($m^* = N_v^{2/3} m^*_b$, menyelaraskan lembah $L$ dan $\Sigma$ hingga $N_v = 16$) meningkatkan $S$ tanpa mengorbankan mobilitas $\mu$.
3. Nanostrukturisasi Hierarkis Semua-Skala dan $\text{SnSe}$: Menggabungkan cacat titik atomik, nanopresipitat endotaksial koheren ($\text{SrTe}$), dan batas butir mesoskopis (Kanatzidis et al., 2012) menghamburkan seluruh spektrum fonon ($ZT \approx 2{,}2$). Pada 2014, kristal tunggal $\text{SnSe}$ mencapai rekor $ZT = 2{,}6 \pm 0{,}3$ pada $923\,\text{K}$ berkat anharmonisitas kisi intrinsik raksasa ($\gamma \approx 4{,}1$, $\kappa_L \approx 0{,}23\,\text{W/(m}\cdot\text{K)}$).

6. Aplikasi Nyata dan Masa Depan: Dari Wahana Antariksa hingga Spinkaloritronik

1. Eksplorasi Antariksa Dalam (RTG & MMRTG): Dari termokopel $\text{SiGe}$ pada Voyager 1 & 2 dan New Horizons hingga MMRTG ($\text{PbTe/TAGS}$) pada rover Mars Curiosity dan Perseverance serta oktokopter nuklir Dragonfly di Titan.
2. Manajemen Termal Presisi (TEC): Stabilisasi suhu $\pm 0{,}01\,\text{K}$ pada laser telekomunikasi DWDM, mesin siklus termal PCR, dan pendinginan bebas getaran untuk detektor CCD dan inframerah.
3. Pemulihan Panas Buang, Termoelektrik Wearable, dan Spinkaloritronik: Sensor IIoT mandiri, semikonduktor anorganik ulet suhu ruang ($\text{Ag}_2\text{S}$) dan polimer konduktif (PEDOT:PSS) pemanen panas tubuh, serta efek Spin Seebeck (SSE) dan efek Nernst Anomali (ANE, $\mathbf{E}_{\text{ANE}} \propto \mathbf{M} \times \nabla T$) pada semimetal Weyl magnetik topologis.