1. Introducción: El milagro del estado sólido que une calor y electricidad sin piezas móviles

Lanzadas a finales del verano de 1977, las sondas espaciales gemelas Voyager 1 y Voyager 2 de la NASA llevan más de 47 años surcando el vacío interestelar —a más de 24.000 millones de kilómetros del Sol— enviando telemetría a la Tierra. A esa distancia, la irradiancia solar cae a menos de una diezmilésima parte de su valor en la Tierra, inutilizando los paneles solares, y cualquier motor térmico con pistones o turbinas habría fallado hace décadas por desgaste mecánico.

El corazón que mantiene vivas a las Voyager es el Generador Termoeléctrico de Radioisótopos (RTG): una fuente de estado sólido que convierte directamente la diferencia de temperatura entre el calor de desintegración alfa del plutonio-238 ($^{238}\text{Pu}$, $T_H \approx 1273\,\text{K}$) y el frío del espacio ($T_C \approx 573\,\text{K}$) en corriente continua mediante termopares semiconductores de silicio-germanio ($\text{SiGe}$). En la Tierra, la misma ley física invertida impulsa las celdas Peltier (refrigeradores termoeléctricos, TEC): sin refrigerantes ni compresores, una corriente eléctrica congela una cara cerámica, estabiliza láseres de fibra óptica DWDM con precisión de $\pm 0{,}01\,\text{K}$, acciona termocicladores PCR y enfría sensores CCD astronómicos.

¿Por qué un gradiente térmico genera voltaje y una corriente bombea calor? ¿Y por qué los módulos termoeléctricos aún no han sustituido a las turbinas de vapor? La respuesta reside en uno de los dilemas más profundos de la física de la materia condensada: la tríada contraindicada (compromiso fundamental) entre el coeficiente Seebeck $S$, la conductividad eléctrica $\sigma$ y la conductividad térmica $\kappa$.

2. Unificación de los tres efectos termoeléctricos clásicos y reciprocidad de Onsager

La termoelectricidad descansa sobre tres fenómenos descubiertos en el siglo XIX: el efecto Seebeck (1821), el efecto Peltier (1834) y el efecto Thomson (1851).

2.1 El efecto Seebeck: Ensanchamiento térmico de Fermi-Dirac y equilibrio difusión-deriva

Cuando un conductor o semiconductor en circuito abierto ($\mathbf{J} = 0$) se somete a un gradiente de temperatura $\nabla T$, surge espontáneamente un campo electrostático interno $\mathbf{E} = -\nabla V$ (más exactamente $-\nabla \tilde{\mu}/e$) proporcional a $\nabla T$:

La constante $S$ (en $\mu\text{V/K}$ o $\text{V/K}$) es el coeficiente Seebeck o poder termoeléctrico. Su origen microscópico se explica mediante la función de distribución de Fermi-Dirac:

1. Ensanchamiento térmico de la superficie de Fermi en el extremo caliente: A $T > 0$, el escalón de Fermi se difumina en una ventana $\sim 2k_B T$. El extremo caliente ($T_H$) posee una población mucho mayor de portadores de alta energía por encima de $E_F$ que el extremo frío ($T_C$).
2. Equilibrio dinámico entre difusión térmica y deriva electrostática: Los portadores calientes se difunden hacia el extremo frío. En un semiconductor tipo n, los electrones ($q = -e$) se acumulan en el lado frío ($S_n < 0$); en un semiconductor tipo p, los huecos ($q = +e$) se acumulan en el lado frío ($S_p > 0$). El campo eléctrico interno resultante impulsa una corriente de deriva opuesta hasta que la corriente de difusión térmica y la corriente de deriva electrostática se cancelan exactamente ($\mathbf{J}_{\text{diff}} + \mathbf{J}_{\text{drift}} = 0$). Conectando ramas p y n térmicamente en paralelo y eléctricamente en serie (unicouple en $\pi$), el voltaje en circuito abierto se suma: $V_{\text{oc}} = (S_p - S_n)(T_H - T_C)$.

2.2 Fórmula de Mott y relación de Pisarenko

Mediante la expansión de Sommerfeld de la ecuación de transporte de Boltzmann, el coeficiente Seebeck de un semiconductor degenerado sigue la fórmula de Mott:

donde $\sigma(E) = \frac{1}{3} e^2 g(E) v^2(E) \tau(E)$ es la conductividad espectral y $g(E)$ la densidad de estados electrónicos (DOS). El coeficiente Seebeck mide la asimetría energética de la conducción alrededor de $E_F$. Para una banda parabólica 3D isótropa con dispersión por fonones acústicos, se obtiene la relación de Pisarenko:

Lograr un alto coeficiente Seebeck requiere (1) una gran masa efectiva de densidad de estados $m^*$ y (2) una baja concentración de portadores $n$ ($S \propto n^{-2/3}$).

2.3 Efecto Peltier, efecto Thomson y relaciones recíprocas de Kelvin-Onsager

Cuando una corriente continua $I$ atraviesa la unión isoterma entre dos materiales A y B, se absorbe o libera calor reversiblemente a razón de:

Microscópicamente, cuando un electrón pasa del metal ($E \approx E_F$) a la banda de conducción de un semiconductor tipo n ($E_{\text{trans}} > E_C$) en la unión fría, debe escalar una barrera de energía $\Delta E = E_{\text{trans}} - E_F > 0$ absorbiendo fonones de la red cristalina y enfriando la unión; en la unión caliente, cae de regreso a $E_F$ liberando su exceso de energía en forma de fonones. Asimismo, el efecto Thomson (1851) describe el intercambio térmico reversible $\dot{q}_{\text{total}} = \frac{\mathbf{J}^2}{\sigma} - \mathcal{K}_T (\mathbf{J} \cdot \nabla T)$ en un conductor homogéneo con un gradiente $\nabla T$.

En la termodinámica lineal fuera del equilibrio de Lars Onsager (1931, Premio Nobel 1968), los flujos eléctrico $\mathbf{J}_e$ y térmico $\mathbf{J}_Q$ cumplen:

La invariancia bajo inversión temporal impone la reciprocidad de Onsager $L_{12} = L_{21}$, de donde se deducen las relaciones de Kelvin:

Esto revela que el coeficiente Seebeck $S = \frac{\Pi}{T}$ es exactamente la entropía transportada por unidad de carga eléctrica ($S = S_{\text{portador}} / q$).

  flowchart TD
    subgraph TEG["1. Modo Generación Seebeck (Generador Termoeléctrico: TEG)"]
        TH1["Fuente caliente T_H (Entrada de calor Q_H de escape / radioisótopo)"] --> N1["Rama semiconductora n: Electrones e⁻ se difunden térmicamente al lado frío"]
        TH1 --> P1["Rama semiconductora p: Huecos h⁺ se difunden térmicamente al lado frío"]
        N1 --> TC1["Sumidero frío T_C (Rechazo de calor Q_C)"]
        P1 --> TC1
        TC1 --> LOAD["Carga externa R_L: La f.e.m. V = (S_p - S_n)(T_H - T_C) genera potencia P"]
    end

    subgraph TEC["2. Modo Refrigeración Peltier (Enfriador Termoeléctrico: TEC)"]
        DC["Fuente CC externa impulsa corriente I"] --> N2["Rama n: Electrones e⁻ escalan la barrera interfacial (E_C - E_F)"]
        DC --> P2["Rama p: Huecos h⁺ entran en la banda de valencia E_V"]
        N2 --> COLD["Unión fría T_C: Absorción de fonones de la red Q_c = (S_p - S_n)T_C I"]
        P2 --> COLD
        COLD --> HOT["Unión caliente T_H: Exceso de energía disipado al sumidero Q_H"]
    end

3. Derivación termodinámica de la figura de mérito $ZT$, la eficiencia máxima y la capacidad de enfriamiento

3.1 Eficiencia máxima $\eta_{\max}$ de un generador termoeléctrico (TEG)

Para un par termoeléctrico entre $T_H$ y $T_C$ ($\Delta T = T_H - T_C$) con resistencia interna $R$, conductancia térmica $K$ y coeficiente Seebeck $S = S_p - S_n$, el flujo de calor absorbido en el lado caliente $Q_H$ y la potencia eléctrica entregada a la carga $R_L$ son:

Al optimizar la geometría para minimizar $RK$, las propiedades del material se agrupan en la figura de mérito adimensional $ZT$ introducida por Abram F. Ioffe (1949):

donde $PF = S^2 \sigma$ es el factor de potencia y $\kappa = \kappa_e + \kappa_L$ la suma de la conductividad térmica electrónica $\kappa_e$ y de red (fonones) $\kappa_L$. Maximizando $\eta(m) = P / Q_H$ respecto a $m = R_L / R$ se obtiene el acople óptimo $m_{\text{opt}} = \sqrt{1 + Z\bar{T}}$ ($\bar{T} = \frac{T_H + T_C}{2}$) y la eficiencia máxima de conversión:

3.2 Potencia de enfriamiento $Q_c$, diferencia máxima de temperatura $\Delta T_{\max}$ y COP en celdas Peltier

En un enfriador Peltier (TEC), el balance térmico en la cara fría ($T_C$) es:

Como el calentamiento de Joule crece cuadráticamente ($\propto I^2$), existe una corriente óptima $I_{\max} = \frac{S T_C}{R}$ ($\frac{\partial Q_c}{\partial I} = 0$). Con carga térmica nula ($Q_c = 0$), la caída máxima de temperatura $\Delta T_{\max}$ y el coeficiente de desempeño máximo $\text{COP}_{\max}$ son:

4. El dilema fundamental del estado sólido: La tríada contraindicada de $S$, $\sigma$ y $\kappa$

¿Por qué es tan difícil superar $ZT \approx 1$? Porque los tres parámetros compiten entre sí a través de la concentración de portadores $n$:

1. Conflicto $S$ vs $\sigma$ (compromiso de Pisarenko): Aumentar $n$ eleva $\sigma = n e \mu$, pero desploma $S \propto m^* n^{-2/3}$. Aumentar la masa efectiva $m^*$ reduce la movilidad $\mu = \frac{e\tau}{m^*}$.
2. Conflicto $\sigma$ vs $\kappa_e$ (ley de Wiedemann-Franz): Los mismos electrones que transportan carga transportan calor según:

donde $L$ es el número de Lorenz. Por ello, el punto óptimo se encuentra en los semiconductores degenerados fuertemente dopados ($n \approx 10^{19}\text{–}10^{20}\,\text{cm}^{-3}$) con brecha energética acorde al criterio de Chasmar-Stratton ($E_g \approx 6\text{–}10\,k_B T$) para suprimir la conducción bipolar.

5. Rompiendo el techo de $ZT \approx 1$: Ingeniería de bandas e ingeniería de fonones

  flowchart TD
    subgraph TRIAD["Cuello de botella clásico (Tríada contraindicada: Límite ZT ≈ 1)"]
        N_UP["Aumentar concentración de portadores n"] --> SIG_UP["Conductividad eléctrica σ = neμ sube (Bueno)"]
        N_UP --> S_DOWN["Coeficiente Seebeck S ∝ n⁻²/³ cae vía Pisarenko (Malo)"]
        SIG_UP --> KE_UP["Conductividad térmica electrónica κ_e = LσT sube vía Wiedemann-Franz (Malo)"]
    end

    subgraph MODERN["Estrategias modernas de desacoplamiento (Ruptura hacia ZT > 2,0–2,6)"]
        BAND["1. Ingeniería de estructura de bandas (Maximizar S²σ)"] --> QD["Confinamiento cuántico (2D/1D/0D) y dopaje resonante (Tl:PbTe): Empinar d ln g(E)/dE"]
        BAND --> CONV["Convergencia de bandas: Alinear valles L y Σ para elevar N_v sin perder movilidad μ"]
        PHONON["2. Ingeniería de fonones (Reducir κ_L al límite vítreo)"] --> PGEC["PGEC e iones 'rattlers' (Skutteruditas rellenas): Dispersión resonante de fonones"]
        PHONON --> ALLSCALE["Arquitectura jerárquica multiescala: Defectos puntuales + nanoprecipitados endotaxiales + granos"]
        PHONON --> ANHARM["Anarmonicidad intrínseca gigante (Monocristal SnSe): Grüneisen γ ≈ 4,1 dispersión Umklapp"]
    end

    TRIAD ==> MODERN

1. El concepto PGEC («Vidrio de fonones – Cristal de electrones», Glen Slack, 1995): En skutteruditas rellenas ($\text{CoSb}_3$) y clatratos, los iones huésped dentro de las jaulas atómicas oscilan de forma anarmónica («rattling») y dispersan los fonones acústicos sin frenar los electrones.
2. Ingeniería de bandas: El confinamiento cuántico (Hicks y Dresselhaus, 1993), el dopaje de nivel resonante ($\text{Tl}:\text{PbTe}$, Heremans et al., 2008) y la convergencia de bandas ($m^* = N_v^{2/3} m^*_b$, alineando los valles $L$ y $\Sigma$ hasta $N_v = 16$) elevan $S$ sin sacrificar la movilidad $\mu$.
3. Nanoestructuración jerárquica multiescala y $\text{SnSe}$: Combinando defectos puntuales atómicos, nanoprecipitados endotaxiales coherentes ($\text{SrTe}$) y bordes de grano mesoscópicos, Kanatzidis et al. (2012) lograron dispersar todo el espectro de fonones ($ZT \approx 2{,}2$). En 2014, el monocristal $\text{SnSe}$ alcanzó $ZT = 2{,}6 \pm 0{,}3$ a $923\,\text{K}$ gracias a su anarmonicidad reticular intrínseca gigante ($\gamma \approx 4{,}1$, $\kappa_L \approx 0{,}23\,\text{W/(m}\cdot\text{K)}$).

6. Aplicaciones reales y fronteras futuras: De las sondas espaciales a la espincaloritrónica

1. Exploración del espacio profundo (RTG y MMRTG): Desde los termopares de $\text{SiGe}$ en Voyager 1 y 2 y New Horizons hasta los MMRTG ($\text{PbTe/TAGS}$) de los rovers marcianos Curiosity y Perseverance y el octocóptero nuclear Dragonfly en Titán.
2. Gestión térmica de precisión (TEC): Estabilización en $\pm 0{,}01\,\text{K}$ de láseres de telecomunicaciones DWDM, termocicladores PCR y enfriamiento sin vibraciones de detectores CCD e infrarrojos.
3. Recuperación de calor residual, termoelectricidad vestible y espincaloritrónica: Sensores IIoT autónomos, semiconductores inorgánicos dúctiles ($\text{Ag}_2\text{S}$) y polímeros conductores (PEDOT:PSS) para recolectar calor corporal, así como el efecto Seebeck de espín (SSE) y el efecto Nernst anómalo (ANE, $\mathbf{E}_{\text{ANE}} \propto \mathbf{M} \times \nabla T$) en semimetales de Weyl magnéticos topológicos.