Overview of Electron Heat-Flow Effects

• Electron heat-flow effect is the phenomenon where electrons regulate thermal energy through mechanisms such as spin-dependent transport, quantum confinement, and kinetic moment interactions.
• It demonstrates that heat transport can occur without net charge flow, relying on context-specific variables like Lorenz ratios, field-aligned fluxes, and transition-state currents.
• This effect has broad applications from nanoscale spin caloritronic devices and quantum-dot heat diodes to plasma, planetary, and astrophysical systems.

Electron heat-flow effect denotes a class of phenomena in which electrons transport, redistribute, rectify, or regulate thermal energy, rather than merely carrying charge. In the cited literature, the effect appears in spin-caloritronic ferromagnets, quantum dots and ballistic quantum circuits, disordered and hydrodynamic conductors, electron–phonon-coupled nanowires, electrons on liquid helium, laboratory plasmas, planetary ionospheres, and the solar wind [(Flipse et al., 2011); (Ruokola et al., 2011); (Sivre et al., 2018); (Yuan et al., 12 Sep 2025); (Bale et al., 2013)]. Across these settings, electron heat flow is governed not by a single constitutive law but by context-specific objects such as spin-dependent Peltier coefficients, Lorenz ratios, field-aligned heat fluxes, transition-state heat currents, and kinetic moment equations.

1. Foundational descriptions and transport variables

The operative heat-flow quantity depends on the physical setting. In spin caloritronics, the relevant object is the spin-dependent Peltier coefficient,

$\Pi_s=\Pi_\uparrow-\Pi_\downarrow .$

In mesoscopic conductors, the comparison between electrical and thermal transport is often phrased through the Wiedemann–Franz relation,

$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$

In space and plasma physics, the central observable is the field-aligned electron heat flux,

$q_\parallel=\vec q\cdot \hat B,$

often normalized by a free-streaming or saturation scale. In donor–acceptor electron transfer across a thermal gradient, the relevant quantity is a bath-resolved heat current that can remain finite even when the net electron flux vanishes [(Flipse et al., 2011); (Majidi et al., 2021); (Bale et al., 2013); (Craven et al., 2016)].

Context	Heat-flow quantity	Representative relation
Spin caloritronics (Flipse et al., 2011)	Spin Peltier coefficient	$\Pi_s=\Pi_\uparrow-\Pi_\downarrow$
Quantum-dot transport (Majidi et al., 2021)	Lorenz ratio	$\kappa_e=L_0TG$, $L=\kappa_e/(GT)$
Solar-wind transport (Bale et al., 2013)	Field-aligned heat flux	$q_\parallel=\vec q\cdot\hat B$
Bithermal electron transfer (Craven et al., 2016)	Steady-state heat current	$$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$

A central conceptual result is that electron heat flow is not synonymous with ordinary electrical conduction. In the donor–acceptor theory of electron transfer across a thermal gradient, an open electron-transfer channel contributes to enhanced heat transport between sites even when they are in electronic equilibrium. At steady state with zero net charge current but nonzero forward and backward hopping fluxes, the net heat current is

$$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$

so heat flows from hot to cold even though the net electron flux vanishes (Craven et al., 2016). This establishes, in a particularly explicit form, that charge neutrality of the steady state does not imply thermal inactivity.

2. Spin-dependent and magnetically controlled heat flow

The most direct magnetically controllable electron heat-flow effect in the cited literature is the spin Peltier effect. In a ferromagnet, spin-up and spin-down electrons carry different amounts of heat because $\Pi_\uparrow\neq \Pi_\downarrow$; a spin current therefore generates a net heat current, whereas in a nonmagnetic metal the two channels are equivalent and their heat currents cancel for a pure spin current. The total heat current is written as

$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$0

and, neglecting heat entering or leaving the stack and Joule heating, the spin-dependent temperature change at a ferromagnet/nonmagnet interface is

$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$1

The physical driver of the extra heat flow is the spin accumulation $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$2 at the interface (Flipse et al., 2011).

Experimentally, the effect was observed in a permalloy/copper/permalloy spin-valve pillar with a platinum bottom contact and a gold top contact. Temperature was measured by a constantan–platinum thermocouple placed on top and electrically isolated from the bottom contact by an $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$3 nm aluminum-oxide layer. An ac lock-in technique separated the first harmonic response $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$4, which scales with $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$5 and captures Peltier-related temperature changes, from the second harmonic response $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$6, which scales with $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$7 and reflects Joule heating. In the parallel configuration, the temperature profile is mainly the conventional Peltier “zigzag” pattern across interfaces; in the antiparallel configuration, spin accumulation at the interfaces produces an additional temperature difference $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$8 at each interface, and the difference between the two magnetic states is $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$9. From thermocouple measurements and 3-D finite-element modeling, the extracted permalloy spin Peltier coefficient is approximately $q_\parallel=\vec q\cdot \hat B,$0 to $q_\parallel=\vec q\cdot \hat B,$1 mV (Flipse et al., 2011).

The significance of this result is the demonstration of magnetic control of heat flow. The reported effect is modest in ferromagnetic metals such as permalloy, but the authors note that it could be enhanced in materials with larger thermopower, such as ferromagnetic oxides, making the spin Peltier effect a potential route to local, programmable cooling in nanoscale electronics (Flipse et al., 2011).

3. Quantum-confined rectification, filtering, and gate control

In quantum-dot and island devices, electron heat flow becomes strongly energy selective. A single-electron heat diode was proposed using two quantum dots or two metallic islands, each tunnel-coupled to its own reservoir but only capacitively coupled to each other. Charge transport through the structure is forbidden, yet energy transfer remains possible through a four-step sequential tunneling cycle that moves one quantum of heat,

$q_\parallel=\vec q\cdot \hat B,$2

In forward bias, $q_\parallel=\vec q\cdot \hat B,$3 and $q_\parallel=\vec q\cdot \hat B,$4, sequential tunneling transfers heat from right to left; in reverse bias, $q_\parallel=\vec q\cdot \hat B,$5 and $q_\parallel=\vec q\cdot \hat B,$6, the cold right reservoir does not have enough thermal energy to activate the necessary transitions, and the cycle is exponentially suppressed by asymmetric Coulomb blockade. The paper states that a proper diode requires $q_\parallel=\vec q\cdot \hat B,$7 to be at least of order $q_\parallel=\vec q\cdot \hat B,$8, that strong rectification typically requires $q_\parallel=\vec q\cdot \hat B,$9, and that for a rectification ratio above $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$0, often $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$1 if the forward current is to remain reasonably large. In favorable cases the rectification ratio can reach $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$2 (Ruokola et al., 2011).

A related but experimentally realized device is the single-quantum-dot heat valve. In a thermally biased single-dot junction, electron temperature maps in the source electrode display Coulomb-diamond-like patterns as functions of gate voltage $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$3 and bias voltage $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$4. Away from charge degeneracy and at zero bias, essentially no heat current passes through the dot; at charge degeneracy and zero bias, the dot transmits energy between hot source and cold drain even when net charge current is zero; at finite bias, Joule heating can dominate the cooling effect near degeneracy. The maximum heat transfer occurs right at the charge degeneracy point, and the observed temperature map shows an ellipsoidal cooling region around that point, with a crossover from cooling to heating as $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$5 or $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$6 is changed (Dutta et al., 2020).

Quantum confinement can also suppress electronic heat flow below the Wiedemann–Franz prediction. In an InAs nanowire transistor with a quantum dot formed spontaneously near pinch-off, the measured Lorenz ratio near an isolated resonance was

$\Pi_s=\Pi_\uparrow-\Pi_\downarrow$7

so the heat conductance was about $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$8 below the Wiedemann–Franz expectation. The effect was attributed to energy-selective transport through a narrow transmission window. At larger gate voltage, where the tunnel couplings were larger by about a factor of $\Pi_s=\Pi_\uparrow-\Pi_\downarrow$9, the Wiedemann–Franz law was nearly restored, with reported Lorenz ratios $\kappa_e=L_0TG$0, $\kappa_e=L_0TG$1, $\kappa_e=L_0TG$2, and $\kappa_e=L_0TG$3, each with uncertainty about $\kappa_e=L_0TG$4 (Majidi et al., 2021).

Taken together, these results show that quantum confinement allows heat current to be rectified, switched, or suppressed independently of the corresponding electrical conductance. This suggests that electron heat-flow effects in nanostructures are governed as much by spectral selectivity and interaction-induced activation thresholds as by the mere existence of a conducting path.

4. Many-body, ballistic, and fluctuation-mediated heat transport

Electron heat flow in mesoscopic circuits is not determined solely by average conductance. In a floating metallic island connected by ballistic channels, a heat Coulomb blockade was observed: electrical conductance remained at the ballistic quantum limit,

$\kappa_e=L_0TG$5

while thermal conductance was reduced to

$\kappa_e=L_0TG$6

Equivalently, the low-temperature electronic heat current became $\kappa_e=L_0TG$7 instead of $\kappa_e=L_0TG$8. The mechanism is the suppression of the single symmetric charge mode by local charge conservation in the floating node, while the $\kappa_e=L_0TG$9 neutral modes remain fully transmitted (Sivre et al., 2018).

A complementary fluctuation-based result was obtained in a quantum circuit consisting of a heated metallic node connected through one tunable quantum point contact and a resistance $L=\kappa_e/(GT)$0. There, thermal shot noise or “delta-$L=\kappa_e/(GT)$1 noise” was measured and directly validated. The generic-channel noise obeyed

$L=\kappa_e/(GT)$2

and the same partition factor $L=\kappa_e/(GT)$3 produced an additional electronic heat-flow contribution when combined with the island’s charging physics. At low temperatures,

$L=\kappa_e/(GT)$4

showing that fluctuations themselves carry and redistribute heat in the circuit (Sivre et al., 2020).

In quantum Hall edge-channel junctions, electron heating changes equilibration thresholds. Two co-propagating cyclotron-split edge channels were brought into interaction over a tunable length $L=\kappa_e/(GT)$5, and the threshold voltage for the onset of radiative inter-edge transitions was found always below the nominal cyclotron gap,

$L=\kappa_e/(GT)$6

with $L=\kappa_e/(GT)$7 meV, and decreasing as $L=\kappa_e/(GT)$8 increased. The dominant mechanism was modeled as electron heating caused by hot-carrier injection through elastic scattering. The local temperature obeyed

$L=\kappa_e/(GT)$9

so longer interaction paths produced stronger heating, which broadened the Fermi distribution and lowered the apparent emission threshold (Paradiso et al., 2011).

Ballistic heat transport through extended one-dimensional networks also departs from single-channel intuition. In a GaAs/AlGaAs 1D waveguide network, the drain temperature increase $q_\parallel=\vec q\cdot\hat B$0 was proportional to $q_\parallel=\vec q\cdot\hat B$1, as expected from Joule’s law, and no temperature increase was observed when the transmitting waveguide was closed, implying negligible reservoir-to-reservoir heat transfer through electron–phonon interaction below $q_\parallel=\vec q\cdot\hat B$2 K. Yet $q_\parallel=\vec q\cdot\hat B$3 was not proportional to the number of populated subbands $q_\parallel=\vec q\cdot\hat B$4, because heat split among multiple branches of the network rather than remaining confined to a single path (Riha et al., 2014).

These studies collectively reject the common simplification that electron heat flow is a direct thermal analogue of dc charge transport. In ballistic and interacting circuits, neutral modes, partition noise, floating-node charge conservation, and path splitting alter thermal transport without necessarily altering electrical conductance in the same way.

5. Electron–phonon coupling, hydrodynamics, and ultrafast nonequilibrium matter

In low-dimensional solids, electron heat flow is often governed by its coupling to phonons and by the structure of disorder. In disordered graphene, Keldysh theory with impurity averaging showed that the low-temperature heat flux $q_\parallel=\vec q\cdot\hat B$5 depends qualitatively on whether the electron–phonon coupling is deformation-potential or vector-potential-like. For weak screening, disorder enhances the low-temperature deformation-potential heat flux and changes the associated power law from $q_\parallel=\vec q\cdot\hat B$6 to $q_\parallel=\vec q\cdot\hat B$7. By contrast, vector-potential coupling is suppressed by disorder, changing $q_\parallel=\vec q\cdot\hat B$8 to $q_\parallel=\vec q\cdot\hat B$9. For strong screening, both couplings yield a $$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$0 low-temperature law. The distinction arises because the charge sector has a diffusion pole, whereas the pseudospin sector does not (Chen et al., 2012).

In InAs/InP heterostructure nanowires containing a quantum dot, the dominant heat path was not direct electronic diffusion through the highly resistive dot but a phonon-mediated sequence: $$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$1 The coupled finite-element model used

$$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$2

with

$$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$3

The measured drain-side electron temperature rise was unexpectedly large even though the quantum-dot resistance was $$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$4, and the authors concluded that electron and phonon temperatures were highly coupled even at temperatures as low as $$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$5 K. They also inferred an apparent Lorenz number

$$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$6

which implies a breakdown of the Wiedemann–Franz law in the InAs portion of the nanowire (Matthews et al., 2012).

Hydrodynamic electron flow introduces another nonlocal structure. For an inhomogeneous electron liquid with a weakly non-uniform momentum relaxation time in a spherical constriction, the absence of viscosity produces a Landauer-dipole-like temperature distribution, asymmetrically deformed along the current by inelastic electron–phonon scattering. The corresponding asymmetry survives in $$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$7, which the paper identifies as a universal feature of inhomogeneous hydrodynamic electron flow. When viscosity is included, it suppresses the thermal Landauer dipole and leads to a hot spot exactly at the center of the constriction (Zhang et al., 2021).

Strongly nonequilibrium electron gases can mimic thermal flow while retaining nonthermal signatures. Electrons floating on liquid $$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$8He and driven by cyclotron-resonance microwaves showed a density-profile-dependent sign reversal of the central density change $$\left(d\mathcal Q_1/dt\right)_{ss}=-\left(d\mathcal Q_2/dt\right)_{ss}$$9: in plateau profiles the central density decreased, while in caldera profiles it increased, consistent with the Poisson–Boltzmann heating picture. Best-fit temperatures were about $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$0 K at $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$1 dBm and about $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$2 K at $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$3 dBm, vastly above the helium bath temperature of $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$4 mK. However, deviations from the thermal model, especially for positive $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$5, suggested that other physical mechanisms can also provide a measurable contribution (Chepelianskii et al., 2018).

Ultrafast nanocluster experiments further separated intrinsic from extrinsic heat flow. In Au$$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$6 nanoclusters on thin-film substrates, intrinsic heat flow referred to energy transfer from hot cluster electrons to the Au lattice, whereas extrinsic heat flow referred to electronic or vibrational coupling between cluster and substrate. A four-temperature model described the coupled dynamics of cluster electrons, cluster lattice, substrate electrons, and substrate lattice. The key structural result was that reversible diffraction-peak broadening appeared only when lattice heating of the nanoclusters was dominated by intrinsic heat flow; it was absent when heat was injected as hot substrate phonons. The disordering rose with a time constant $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$7 ps and scaled as $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$8 with $$\left(\frac{d\mathcal Q_1}{dt}\right)_{ss} = -\left(\frac{d\mathcal Q_2}{dt}\right)_{ss} = \mathcal J_{ss}\, \frac{2E_{\mathrm{R}1}E_{\mathrm{R}2}(T_2-T_1)} {E_{\mathrm{R}1}T_1+E_{\mathrm{R}2}T_2},$$9, which the authors interpreted as hot electrons modifying the potential energy surface and activating surface diffusion even below the equilibrium threshold for surface pre-melting (Vasileiadis et al., 2018).

6. Plasma, planetary, and astrophysical manifestations

In partially ionized planetary ionospheres, electron heat flow modifies the ambipolar electric field. A diffusion theory based on the eight-moment approximation showed that the standard ambipolar field derived from the electron pressure gradient alone omits a heat-flow term arising when electron temperature varies along magnetic field lines. The heat-flow-inclusive field was denoted $\Pi_\uparrow\neq \Pi_\downarrow$0, the classical field $\Pi_\uparrow\neq \Pi_\downarrow$1, and the comparison quantified by $\Pi_\uparrow\neq \Pi_\downarrow$2. A critical parallel electron temperature gradient was identified,

$\Pi_\uparrow\neq \Pi_\downarrow$3

Comparison with Endurance sounding-rocket measurements below $\Pi_\uparrow\neq \Pi_\downarrow$4 km gave potential-drop slopes of about $\Pi_\uparrow\neq \Pi_\downarrow$5 on ascent and $\Pi_\uparrow\neq \Pi_\downarrow$6 on descent. Classical ambipolar theory without heat flow systematically underestimated the measured potential drop, especially during descent, whereas the heat-flow-inclusive solution aligned much better with the data. The same framework was then used to address Venus’s electric potential drop anomaly: for $\Pi_\uparrow\neq \Pi_\downarrow$7 K and $\Pi_\uparrow\neq \Pi_\downarrow$8, the reported ambipolar field was about $\Pi_\uparrow\neq \Pi_\downarrow$9 without heat flow and $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$00 with heat flow for $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$01 (Yuan et al., 12 Sep 2025).

In the solar wind at $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$02 AU, the field-aligned electron heat flux exhibits a collisional-to-collisionless transition. Spitzer–Härm theory predicts

$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$03

with $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$04. Measurements showed that

$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$05

equivalently $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$06, representing about $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$07 of the data. For $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$08, the flux saturated at

$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$09

independent of further increases in $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$10; about $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$11 of the data lay in this regime. The authors found that the collisionless regulation was not obviously consistent with a whistler heat-flux instability and was more compatible with magnetosonic regulation, wave–particle scattering, an interplanetary electric potential, or a general kinetic flux-limitation mechanism (Bale et al., 2013).

Near the Sun, Parker Solar Probe measurements between $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$12 and $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$13 AU showed that the absolute heat flux was anticorrelated with solar-wind speed, while the normalized net heat flux $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$14 was anticorrelated with plasma beta on all orbits. In high-$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$15 regions near the heliospheric current sheet, the net heat flux often decreased even though the omnidirectional suprathermal electron flux remained comparable or increased. The paper emphasized that many such dropouts were therefore not disconnections from the Sun; a true disconnection would require both low net heat flux and a drop in omnidirectional suprathermal flux. The observations were reported as inconsistent with regulation primarily by collisional mechanisms near the Sun and consistent instead with theoretical instability thresholds associated with oblique whistler and magnetosonic modes (Halekas et al., 2020).

In fusion-relevant plasmas, electron heat flow can be strongly non-diffusive. In KSTAR tokamak core plasmas, fast electron heat transport events were identified as non-diffusive avalanche-like processes. The radial propagation speed was about $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$16 in L-mode and about $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$17 in weak-ITB plasmas, with an inferred escape time of about $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$18 ms, roughly $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$19 times faster than the energy confinement time. The fluctuation spectrum obeyed $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$20, the Hurst exponents were $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$21 and $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$22, and the observed electron-temperature profile corrugation had width roughly $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$23, interpreted as a mesoscale $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$24 shear-flow structure. The long-range avalanche-like events occurred when the profile corrugation was destroyed (Choi et al., 2018).

A different plasma manifestation appears in magnetized hohlraums, where kinetic modeling shows that electron heat flow, inverse bremsstrahlung, and Nernst advection are tightly coupled. The generalized heat flux contains not only the classical conduction term $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$25 but also an anomalous term $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$26, identified as anomalous heat flow up a density gradient. The Nernst velocity was written as

$\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$27

linking magnetic-field advection directly to electron heat flow. For an initial $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$28 T field, the field on axis grew to about $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$29 T within $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$30 ns, the wall field reached nearly $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$31 T, and magnetic-field amplification in the wall was about a factor of $\kappa_e=L_0TG,\qquad L\equiv \frac{\kappa_e}{GT}.$32 (Joglekar et al., 2015).

The comparative literature therefore shows that electron heat-flow effects are not a single phenomenon but a recurrent transport motif. Heat can be carried by spin imbalance, energy-selective tunneling, neutral modes, partition noise, phonon-assisted pathways, hydrodynamic convection, avalanche-like turbulence, or collisionless kinetic populations. This suggests a general principle: whenever electron distributions are structured by spin, confinement, interactions, disorder, or weak collisionality, thermal transport can decouple from the simplest charge-transport intuition and acquire its own control variables, thresholds, and instabilities.