---
title: Electron Transpiration Cooling (ETC)
url: https://www.emergentmind.com/topics/electron-transpiration-cooling-etc
type: topic
---

# Electron Transpiration Cooling (ETC)

Electron transpiration cooling (ETC) denotes electron-mediated heat removal by preferentially extracting energetic electrons from a hot region. In the cited literature, the term spans two distinct physical realizations. In cryogenic semiconductor devices, electron transpiration is the energy-selective tunnelling of hot electrons into a superconductor, reducing the electron-gas temperature in a semiconductor island [1405.4119]. In hypersonic aerothermodynamics, ETC refers to thermionic emission from a hot leading edge, with cooling determined by whether emitted electrons escape, are recollected, or redistribute heat over the surface [2502.07387; 2508.05900]. Across both settings, ETC is governed by an energy filter, a return path for parasitic heat, and a global transport constraint that determines whether nominal emission translates into net cooling.

## 1. Definitions and physical regimes

In the low-temperature solid-state formulation, ETC is implemented with superconducting tunnel contacts that selectively remove electrons with energies above the Fermi level by more than the superconducting gap $\Delta$. Prest et al. describe this as electron cooling in silicon using platinum silicide as a superconductor contact, where hot electrons tunnel out of the semiconductor and the electron gas cools as a result [1405.4119].

In the hypersonic formulation, ETC is identified with thermionic cooling, also described as electron transpiration cooling, for leading edges that can reach temperatures exceeding $2000\,^\circ\mathrm{C}$. There, the emitted electron population is set by thermionic emission and, in some models, augmented by photoemission. Cooling is not determined by emission alone: reflected electrons, incident plasma electrons, ion effects, sheath structure, and downstream collection all enter the energy balance [2502.07387; 2508.05900].

A useful unifying interpretation is that ETC always depends on energy-selective charge removal from the emitting surface or electron gas, followed by a transport problem. In the superconductor–semiconductor case, the filter is the superconducting density of states. In hypersonic ETC, the filter is the work function plus sheath or virtual-cathode barrier. This suggests that ETC is less a single device class than a family of nonequilibrium electron-transport cooling mechanisms.

## 2. Cryogenic ETC in superconductor–semiconductor junctions

Prest et al. formulate cooling power per tunnel junction as
$$
\dot Q(V,T_e,T_s)=\frac{1}{e^2R_N}\int_{-\infty}^{\infty}(E-eV)\,n_S(E)\,[f_N(E-eV,T_e)-f_S(E,T_s)]\,dE,
$$
where $R_N$ is the normal-state junction resistance, $V$ the bias per junction, and $f_N$ and $f_S$ the Fermi functions in the normal and superconducting electrodes at temperatures $T_e$ and $T_s$ respectively [1405.4119]. Physically, electrons with energies above the Fermi level by more than $\Delta$ tunnel out, removing heat.

The superconducting density of states is written in Dynes form,
$$
n_S(E)=\mathrm{Re}\left[\frac{E+i\Gamma}{\sqrt{(E+i\Gamma)^2-\Delta^2(T)}}\right],
$$
with $\Gamma$ the Dynes sub-gap leakage parameter and $\Delta(T)$ following the BCS temperature dependence, $\sim\Delta_0\tanh[1.74\sqrt{T_c/T-1}]$, with $\Delta_0\simeq1.764\,k_B T_c$ at $T\ll T_c$ [1405.4119]. The role of $\Gamma$ is central: sub-gap leakage enables unwanted low-energy tunnelling and directly degrades cooling.

The dominant parasitic load is electron–phonon coupling,
$$
P_{e-\mathrm{ph}}=\Sigma V(T_e^5-T_{\mathrm{lat}}^5),
$$
where $\Sigma$ is the material-specific coupling constant, $V$ the electron-gas volume, and $T_{\mathrm{lat}}$ the lattice temperature [1405.4119]. Because the dependence is $T^5$, the lattice heat load falls rapidly at sub-kelvin temperatures, enabling net cooling. Heat balance for an S–Sm–S device is then written as
$$
\dot Q(V,T_e,T_s)+P_{e-\mathrm{ph}}+P_J=0,
$$
with $P_J=I^2R_s$ the residual Joule heating in series resistances [1405.4119].

The device realization uses a 10 nm PtSi thin film on silicon-on-insulator with a 140 nm buried oxide. The central island is implanted with As to $\sim8\times10^{19}\,\mathrm{cm^{-3}}$ and has sheet resistance $\simeq100\,\Omega/\square$. Each PtSi–Si contact has area $2.5\,\mu\mathrm{m}\times5\,\mu\mathrm{m}$, and two junctions in series form an S–Sm–S cooler. The measured normal-state resistance is $R_N\simeq300\,\Omega$ per junction, with series island resistance $R_s\simeq320\,\Omega$. Aluminium pads overlap PtSi for four-point I–V measurements, eliminating contact-resistance artefacts [1405.4119].

## 3. Thin-film PtSi realization and measured cooling

The PtSi implementation is motivated by gap engineering. Bulk PtSi has $T_c\simeq1.0\,\mathrm{K}$, whereas a 10 nm film suppresses $T_c$ to $\simeq0.79\,\mathrm{K}$, corresponding to $\Delta_0\simeq70\,\mu\mathrm{eV}$ [1405.4119]. A smaller $\Delta$ reduces sub-gap leakage currents $\sim\Gamma/\Delta$ and shifts the optimum cooling bias to lower voltages, improving performance at $T_{\mathrm{bath}}\lesssim100\,\mathrm{mK}$. Thin films also crystallize better, yielding smoother PtSi–Si interfaces and higher-quality tunnel barriers [1405.4119].

At $T_{\mathrm{bath}}=100\,\mathrm{mK}$, the measured I–V characteristics show strong suppression of sub-gap current near $V\simeq\pm0.1\,\mathrm{mV}$, while $dI/dV$ exhibits the classic U-shape of a superconducting tunnel junction with rounded peaks at $\pm2\Delta/e$ [1405.4119]. Fitting the isotherm model with fixed $T_e=100\,\mathrm{mK}$ gives $\Delta_0\simeq70\,\mu\mathrm{eV}$ and $\Gamma/\Delta\simeq8\times10^{-3}$. When the heat-balance model is incorporated, the bias-dependent electron temperature drops from $100\,\mathrm{mK}$ to $\sim50\,\mathrm{mK}$ at $V\simeq0.12\,\mathrm{mV}$ [1405.4119].

Peak cooling power at $100\,\mathrm{mK}$ is $P_c\simeq0.72\,\mathrm{pW}$, or $\sim29\,\mathrm{fW/\mu m^2}$, and the coefficient of performance is described as typically a few percent at optimum bias [1405.4119]. The lower $\Delta$ of PtSi allows operation down to $50\,\mathrm{mK}$, below the regime easily accessed with Al, for which $\Delta_{\mathrm{Al}}\simeq180\,\mu\mathrm{eV}$ [1405.4119].

The limiting factors are explicit. Dynes sub-gap leakage of $\Gamma/\Delta\simeq8\times10^{-3}$ is higher than in state-of-the-art NIS coolers, where $10^{-5}$–$10^{-4}$ is typical. Series resistance narrows the usable bias window and introduces Joule heating, while residual thermal coupling to the lattice and environmental heat leaks set a floor to the minimum achievable $T_e$ [1405.4119]. The same Schottky-barrier tunnel junctions can in principle operate as thermometers as well as refrigerators, and PtSi/Si coolers could be monolithically integrated with silicon qubits, single-electron pumps, or ultra-low-noise detectors [1405.4119].

## 4. Thermionic and photoemission ETC for hypersonic surfaces

For hypersonic leading edges, the fundamental thermionic-emission law is the Richardson–Dushman relation
$$
J_{\rm th}(T)=A\,T^2\exp\!\left[-\frac{\Phi}{k_B T}\right],
$$
with local cooling power density estimated as
$$
q_{\rm cool}^{\rm th}(T)=J_{\rm th}(T)\,(\Phi+2k_B T)
$$
when each escaping electron carries away an average energy on the order of $\Phi+2k_B T$ [2502.07387]. Boyer and Fisher write the same ideal thermionic cooling contribution per unit area as
$$
q_{\rm emit}=J_t\,(q\phi_{WF}+2k_B T_w),
$$
with $J_t$ the actual emission current, emphasizing that the useful cooling rate depends on the emission that survives sheath constraints rather than the Richardson–Dushman value alone [2508.05900].

Ghosh and Fisher extend the formulation to photoemission using a modified Fowler–DuBridge description. The spectral electron emission density including photon energy $\hbar\omega$ is written as $J_{\rm PEED}^{\rm work}(u,\theta_p,T)$ and, with a virtual-cathode potential included, as $J_{\rm PEED}^{\rm work+VC}(u,\theta_p,T)$. The reflected current density per spectral bin is then
$$
J_{\rm PEED}^{\rm refl}(u,\theta_p,T)=\tilde n_{\rm norm}(u,\theta_p,T)\Bigl[J_{\rm PEED}^{\rm work}(u,\theta_p,T)-J_{\rm PEED}^{\rm work+VC}(u,\theta_p,T)\Bigr],
$$
where $\tilde n_{\rm norm}$ is the normalized 3D random-energy probability density [2502.07387]. In this framework, net cooling is the difference between emitted and reflected electron energy fluxes.

The sheath or virtual-cathode structure is prescribed in Cartesian, cylindrical, and spherical coordinate systems. The reported sheath depths and extents are $-0.175\,\mathrm{V}$ and $\sim21\,\mu\mathrm{m}$ for Cartesian, $-0.146\,\mathrm{V}$ and $\sim25\,\mu\mathrm{m}$ for cylindrical, and $-0.168\,\mathrm{V}$ and $\sim30\,\mu\mathrm{m}$ for spherical geometry [2502.07387]. Electron trajectories are launched from a random-energy distribution,
$$
\tilde n_{3D}(u,\theta_p,T)=T\ln\!\Bigl[1+\exp\bigl[(E_F-\tfrac12m_eu^2)/(k_B T)\bigr]\Bigr]\sin\theta_p,
$$
and propagated under
$$
\dot z=v_z,\qquad \ddot z=-\frac{q_e}{m_e}\frac{d\phi(z)}{dz}
$$
using a Runge–Kutta 5(4) ODE integrator [2502.07387].

The surface geometry in that study is a planar disk of radius $50\,\mu\mathrm{m}$, representing a local stagnation-point region. Two temperature profiles are examined: a step function with $T=2000\,\mathrm{K}$ for $0\le r\le20\,\mu\mathrm{m}$ and $T=1500\,\mathrm{K}$ for $r>20\,\mu\mathrm{m}$, and a Gaussian profile $T(r)=2000-500\exp[-(r/15\,\mu\mathrm{m})^2]\,\mathrm{K}$ [2502.07387]. The characteristic gradient length is stated to be comparable to the electron lateral-travel distance, making ETC partly a heat-spreading mechanism rather than only a local heat-removal mechanism.

Photoemission contributes materially in the visible-light model with $\omega=3.73\times10^{15}\,\mathrm{rad/s}$ and $\hbar\omega\approx2.5\,\mathrm{eV}$, where total emitted current can increase by $10$–$20\%$ over pure thermionic values at $2000\,\mathrm{K}$, and preliminary tests show $J_{\rm emitted}$ up to $1.2\times$ larger when photoemission is included [2502.07387]. A plausible implication is that ETC optimization for hypersonic surfaces may require co-design of temperature field, work function, and illumination spectrum.

## 5. Kinetic transport limits, backflow, and misleading cooling metrics

Zhang et al. model ETC-relevant transport with a one-dimensional-in-space, three-dimensional-in-velocity electrostatic PIC–MCC plasma diode spanning a full cathode–anode gap [2605.07083]. The domain length is $L=0.10\,\mathrm{m}$ with $N_z=4096$ cells, $\Delta z\approx24.4\,\mu\mathrm{m}\approx0.04\,\lambda_D$, time step chosen so that $\omega_{pe}\Delta t\approx0.008$, and $N_{\rm ppc}=1000$ particles per cell [2605.07083]. Poisson’s equation is solved with $\phi(0)=0$ and $\phi(L)=40\,\mathrm{V}$; thermionic emission is injected at the cathode as a prescribed flux $\Gamma_{\rm emit}$ sampled from a half-Maxwellian at $T_c=0.2\,\mathrm{eV}$ [2605.07083].

The diagnostic structure distinguishes emitted flux, reflected flux, net transport, and anode collection:
$$
\Phi_{\rm emit}\equiv\Gamma_{\rm emit},\qquad
\beta=\frac{\Gamma_{\rm ref}}{\Gamma_{\rm emit}},\qquad
\eta_{\rm net}=1-\beta,\qquad
\eta_{\rm anode}=\frac{\Gamma_{\rm anode}}{\Gamma_{\rm emit}}.
$$
Cathode-side cooling metrics are defined as
$$
q_{\rm cooling,1}=q_{\rm emit,theory}-q_{\rm back,em,cath}-q_{\rm back,p,cath}-q_{\rm ion,cath},
$$
and
$$
q_{\rm cooling,2}=q_{\rm cooling,1}-q_{\rm ion,cath,recomb},
$$
with the ion-recombination term written as $q_{\rm ion,cath,recomb}=\Gamma_{\rm ion,cath}\,e\,E_{\rm ion}$ and $E_{\rm ion}=24.59\,\mathrm{eV}$ [2605.07083].

The principal numerical result is a sharp transition from weak-backflow transport to backflow-limited transport. Below $7.0\times10^{19}\,\mathrm{m^{-2}s^{-1}}$, the backflow ratio remains below $13\%$ and $\eta_{\rm net}\approx1$. At $7.25\times10^{19}\,\mathrm{m^{-2}s^{-1}}$, $\Gamma_{\rm ref}=3.92\times10^{19}\,\mathrm{m^{-2}s^{-1}}$, $\beta\approx54.0\%$, and $\eta_{\rm net}\approx46.0\%$. Above $7.5\times10^{19}\,\mathrm{m^{-2}s^{-1}}$, further emission produces over-compensation: $\beta$ rises to $\sim85\%$ while $\Phi_{\rm net}$ and $\Phi_{\rm anode}$ decrease [2605.07083].

The transition is associated with full-gap potential restructuring. Near $7.25\times10^{19}\,\mathrm{m^{-2}s^{-1}}$, the interior potential collapses and the rise is concentrated near the anode, weakening the near-cathode barrier and promoting backflow [2605.07083]. A purely local virtual-cathode estimate,
$$
\eta_{\rm VC}=\exp\!\Bigl[\frac{\phi_{\rm min}-\phi_w}{T_{\rm emit}}\Bigr],
$$
would give $\eta_{\rm VC}\approx0.99$ for the shallow local dip $\phi_{\rm min}\approx-2.2\times10^{-3}\,\mathrm{V}$, whereas the full-gap effective barrier
$$
\phi_{\rm vc,eff}=T_{\rm emit}\ln\!\Bigl(\frac{\Gamma_{\rm net}}{\Gamma_{\rm emit}}\Bigr)\approx-0.156\,\mathrm{V}
$$
shows that the loss is dominated by global sheath and transport collapse rather than a single local dip [2605.07083].

This result directly addresses a common misconception: stronger emission does not necessarily improve ETC. Zhang et al. show that cathode-side cooling metrics may continue to grow after the transition because they count energy removed at emission, even when emitted electrons fail to escape the full gap [2605.07083]. In that regime, boundary-emission diagnostics cease to be reliable proxies for useful ETC transport.

## 6. Passive ETC energetics, system-level constraints, and design window

Boyer and Fisher analyze passive ETC for a discretized $1\,\mathrm{cm}$-radius tungsten leading edge with wall temperature near $2000\,\mathrm{K}$ at stagnation, embedding a one-dimensional collisionless plasma sheath model into a leading-edge framework [2508.05900]. The sheath model uses
$$
\frac{d^2\phi(x)}{dx^2}=-\frac{\rho(x)}{\varepsilon_0},
$$
with wall and sheath-edge boundary conditions, flowfield electron flux
$$
J_{fe}=n_{e0}\sqrt{\frac{k_B T_{e0}}{2\pi m_e}}\exp\!\Bigl(\frac{q\phi_w}{k_B T_{e0}}\Bigr),
$$
and ion flux
$$
J_i=n_{i0}u_B,\qquad u_B=\sqrt{\frac{k_B T_{e0}}{m_i}}.
$$
The effective work function is reduced by the Schottky effect:
$$
\phi_{WF,\rm eff}=\phi_{WF}-\frac{1}{q}\sqrt{\frac{q^3|E_w|}{4\pi\varepsilon_0}}.
$$
The maximum floating-emitter current with cold ions is given as
$$
J_{c,\rm floating}=\bigl[C_{\rm crit}(1-C_{\rm crit})\bigr]J_i,\qquad
C_{\rm crit}=1-8.3\sqrt{\frac{m_e}{m_i}},
$$
and the actual emission current must be found iteratively with the sheath solution [2508.05900].

The complete surface energy balance includes convective–kinetic heating, radiative heating, thermionic cooling, flowfield electron and ion collection heating, and Joule heating:
$$
q_{\rm net}=q_{\rm conv}+q_{\rm rad}-q_{\rm emit}+q_{\rm collect},
$$
or, more explicitly,
$$
\begin{aligned}
q_{\rm net} &=
J_{fe}\bigl[q\phi_{WF}+2k_B T_{e0}\bigr]
-J_i\Bigl[q\phi_{WF}-\bigl(q\phi_w+\tfrac12 m_i u_B^2+I_0\bigr)\Bigr] \\
&\quad
-J_t\bigl[q\phi_{WF}+2k_B T_w\bigr]
+Q_j .
\end{aligned}
$$
This formulation makes the central systems result explicit: ETC cooling must be compared against the heating caused by recollected flowfield electrons and ions, not only against emission-side energy removal [2508.05900].

The parametric study reveals a narrow operating window. Defining $Y_0\approx2.8\times10^6$ as the minimum seeding needed to just prevent a virtual cathode, corresponding to $n_{i0}\approx9.1\times10^{19}\,\mathrm{m^{-3}}$, the Schottky-enhanced Richardson model predicts $\sim4.1\,\mathrm{MW/m^2}$ cooling at $Y=Y_0$, while the full ETC circuit yields only $\sim3.0\,\mathrm{MW/m^2}$ [2508.05900]. At $Y=10Y_0$, Schottky theory increases cooling by $\sim25\%$, yet the full circuit loses $44\%$ of cooling, to $\sim1.7\,\mathrm{MW/m^2}$. At $Y\gtrsim50Y_0$, the circuit reverses and heats the leading edge by about $+1.9\,\mathrm{MW/m^2}$ [2508.05900].

That reversal is attributed to the floating potential at the stagnation point becoming less positive than downstream, so charge conservation drives electrons away from the stagnation region and deposits heat instead of extracting it [2508.05900]. Raising the work function from $2.5\,\mathrm{eV}$ to $3.5\,\mathrm{eV}$ at $Y=Y_0$ reduces emission until the ETC circuit produces net heating instead of cooling, and above $\sim4.5\,\mathrm{eV}$ heating grows rapidly [2508.05900].

The same study identifies blunt radii and dielectric coatings as protective but ETC-suppressing. Blunter radii of $5$–$20\,\mathrm{cm}$ behave more like floating surfaces, so ETC is essentially off. An insulating coating forces $\phi_w=\phi_f$ everywhere, reducing adverse overcompensation heating by up to $\sim70\%$ but also blocking net cooling at stagnation [2508.05900]. Reducing wall resistivity by two orders of magnitude recovers $\sim40\%$ of lost cooling, after which diminishing returns set in; collector lengths beyond $\sim1\,\mathrm{m}$ likewise show diminishing returns, while even $0.1\,\mathrm{m}$ retains $\sim45\%$ of nominal cooling [2508.05900].

The design guidance emerging from these results is internally constrained rather than monotonic. Passive ETC should operate in a narrow window near exact or slight under-compensation, stated as $Y\in[Y_0,10Y_0]$ for the examples studied; use low-$\phi_{WF}$ materials or surface treatments in the range $2.0$–$2.5\,\mathrm{eV}$; maintain high electrical conductivity without driving the system into severe overcompensation; and use sharp leading-edge radii near $1\,\mathrm{cm}$ if ETC is to remain active [2508.05900]. This suggests that ETC for hypersonic vehicles is not limited by raw thermionic capability but by a coupled electrothermal circuit condition in which collection heating, sheath stability, and geometry can nullify or reverse nominal cooling.

Source: https://www.emergentmind.com/topics/electron-transpiration-cooling-etc