Sub-Kelvin Thermal Conductivity

• Sub-Kelvin thermal conductivity is the metric that quantifies heat transport in materials below 1 K, essential for quantum device operation and low-temperature studies.
• Steady-state techniques in dilution refrigerators and numerical differentiation reveal distinct power-law dependencies, highlighting mechanisms like boundary and TLS scattering.
• Material-specific properties from HR/LR Si, sapphire, Nb films, and amorphous polymers guide cryogenic system design and optimize thermal management in quantum technologies.

Sub-Kelvin thermal conductivity quantifies the ability of materials to conduct heat at temperatures below 1 kelvin, a regime central to quantum device operation, ultra-sensitive calorimetry, and fundamental studies of phonon, magnon, and electron thermal transport. Heat conduction mechanisms at sub-Kelvin temperatures differ fundamentally from higher-temperature regimes, becoming dominated by boundary scattering, disorder-induced mechanisms (e.g., two-level systems in glasses), and, in certain cases, quantum topological excitations. Experimental results span a variety of conductors, insulators, superconductors, and low-dimensional systems, revealing material-specific power-law temperature dependencies and scattering phenomena that are critical for thermal management and device design at millikelvin temperatures.

1. Experimental Methodologies for Sub-Kelvin Thermal Conductivity

Sub-Kelvin κ(T) is typically extracted via steady-state techniques in dilution refrigerators, where controlled heater power is applied at one end of a high-purity sample, and micromachined or thin-film resistive thermometers (e.g., RuO₂, NIS junctions) monitor temperature differences. The "integrated thermal conductivity method" measures the applied heater power $P$ versus induced temperature difference ΔT across known sample geometry, with κ(T) deduced by numerical differentiation:

$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$

where $L$ is the distance between thermometers and $A$ is the cross-sectional area. Systematic uncertainties—from calibration, parasitic heat leaks, wiring shunts, and contact resistances—are minimized by careful thermalization, calibration against primary thermometers (e.g., $^{60}$Co nuclear orientation), and statistical analysis of measurement repeatability, often propagating errors to yield χ²/ndf = 1 in reported fits (Drobizhev et al., 2016, Erhart et al., 2023, Bon-Mardion et al., 7 May 2026, Feshchenko et al., 2016).

Special considerations pertain to on-chip and low-dimensional systems. In nanoscale devices, power levels below 10 nW are required to keep temperature gradients within measurement linearity. In 2D conductors, local electronic temperatures are measured via temperature-dependent universal conductance fluctuations (UCF), with superconducting contacts employed to block out-diffusion from the sample and isolate intrinsic cooling mechanisms (Draelos et al., 2018).

2. Material-Dependent Thermal Conductivity and Power Laws

Sub-Kelvin κ(T) is highly material- and microstructure-dependent, reflecting underlying scattering mechanisms. Representative power-law fits and scaling exponents are summarized as follows:

Material	κ(T) W/m·K	Power Law Exponent n	Dominant Scattering Mechanism	Reference
High-resistivity Si	$5\times10^{-2}$	$n=2.4$	Boundary scattering (boundary-limited Casimir)	(Bon-Mardion et al., 7 May 2026)
Low-resistivity Si	$8\times10^{-4}$	$n=2.2$	Phonon–hole (impurity) scattering	(Bon-Mardion et al., 7 May 2026)
Sapphire	$2\times10^{-3}$	$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$0	Grain boundary/structural disorder	(Bon-Mardion et al., 7 May 2026)
Borosilicate glass	$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$1	$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$2	Two-level system (TLS, amorphous)	(Bon-Mardion et al., 7 May 2026)
PTFE (amorphous)	see text	$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$3	Tunneling two-level systems	(Drobizhev et al., 2016)
Al₂O₃ ceramic	$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$4 (@0.7 K)	$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$5–$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$6 (above 0.3 K), $\kappa(T) = \frac{L}{A} \frac{dP}{dT}$7 (below)	Crystallite boundary + defect scattering	(Drobizhev et al., 2016)
Nb thin films	$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$8–$\kappa(T) = \frac{L}{A} \frac{dP}{dT}$9 (@0.3-0.6K)	$L$0	Residual conduction in the superconducting state	(Feshchenko et al., 2016)
Polystyrene (UPS-923A)	$L$1	$L$2	TLS/scattering in amorphous polymer	(Erhart et al., 2023)
Graphene (electronic)	$L$3	$L$4 (WF)	Electronic diffusion + e-ph coupling	(Draelos et al., 2018)

Power-law exponents generally reflect the dominant phonon scattering regime: $L$5 for Debye-boundary-limited crystals, $L$6 for TLS-limited amorphous solids, and non-integer exponents in disordered or doped systems where multiple mechanisms compete (Drobizhev et al., 2016, Bon-Mardion et al., 7 May 2026, Erhart et al., 2023).

3. Scattering Mechanisms and Theoretical Models

Sub-Kelvin heat conduction is generally governed by acoustic phonons, with the mean free path $L$7 set by extrinsic (boundary, grain, impurity) and intrinsic (TLS, e-ph interaction) scattering:

• Crystalline materials (e.g., sapphire, HR Si): For HR Si, diffusive boundary scattering dominates ($L$8 mm scale, Casimir), yielding slightly sub-cubic κ(T) ($L$9), attributed to partial mode conversion at boundaries. Sapphire shows $A$0 50 μm, T-indep., consistent with grain-boundary scattering from sub-μm crystallites, enforcing $A$1 but with reduced magnitude (Bon-Mardion et al., 7 May 2026).
• Doped/defected semiconductors (e.g., LR Si): Phonon–hole scattering rates $A$2 produce $A$3, and typical $A$4 values are much shorter (~30 μm at 0.5 K) than in pure silicon (Bon-Mardion et al., 7 May 2026).
• Amorphous solids and polymers (PTFE, borosilicate, polystyrene): TLS models predict $A$5, verified in PTFE, borosilicate, and polystyrene-based scintillators, where measured n is within experimental uncertainty of 2 (Erhart et al., 2023, Drobizhev et al., 2016).
• Ceramics (Al₂O₃): Sintered alumina samples exhibit mixed scattering. Above 0.3 K, $A$6 (long-range Debye phonons between crystallites), while below 0.3 K a linear-in-T term emerges, reflecting elastic defect scattering localized within the microstructure (Drobizhev et al., 2016).
• 2D systems (graphene): In-plane electronic thermal conduction is governed by the Wiedemann-Franz law, $A$7, while electron-phonon cooling is described by $A$8, with measured $A$9 and $^{60}$0 at carrier densities $^{60}$1 cm$^{60}$2 (Draelos et al., 2018).
• Superconductors (Nb thin-film): Observed $^{60}$3 between 0.1 and 0.6 K exceeds the BCS-predicted exponential suppression by orders of magnitude. $^{60}$4 is two orders smaller than the normal-state Wiedemann-Franz conductance and two orders lower than superconducting Al, making thin Nb an effective thermal insulator for sub-Kelvin circuitry (Feshchenko et al., 2016).

4. Hall Thermal Transport and Exotic Contributions

Thermal transport at sub-Kelvin temperatures can be highly nontrivial in materials supporting topological or correlated excitations. In the frustrated pyrochlore magnet Yb₂Ti₂O₇, a "gigantic" thermal Hall effect (κ{xy}) is observed below 1 K. The thermal Hall angle (|κ{xy}|/κ{xx}) reaches nearly 2%, with a strong dependence on the field and temperature regime—peaking well above the transition temperature to canted ferromagnetic order (T{CFM} ≈ 0.275 K). Unlike conventional phonon skew scattering, the behavior is attributed to the magnon Hall effect (Berry curvature of magnon bands), with a crossover to negative κ{xy} as field-aligned phases set in and the magnon gap opens. This non-vanishing κ{xy} at T < 0.3 K demonstrates the sensitivity of sub-Kelvin thermal transport to collective quantum excitations and topological effects (Hirschberger et al., 2019).

5. Application to Quantum Technologies and Cryogenic System Design

Sub-Kelvin thermal conductivity is a critical consideration in the architecture and scaling of quantum devices:

• Substrate selection: HR Si provides maximum thermal evacuation for control/readout stages (κ(300 mK) = $^{60}$5 W/m·K), while LR Si, sapphire, and borosilicate offer lower conductance suitable for qubit thermal isolation (κ(300 mK) = $^{60}$6–$^{60}$7 W/m·K) (Bon-Mardion et al., 7 May 2026).
• On-chip routing: Superconducting Nb lines increase in-plane conductance by up to a factor of four on thinned LR-Si-based SOCs at 0.5 K, yet the substrate remains the dominant path for heat—an important design parameter for large-scale integration (Bon-Mardion et al., 7 May 2026).
• Integration strategies: To prevent unwanted thermal cross-talk, function partitioning using 3D integration (e.g., SIP, stacked architectures) is essential, with separate substrates or interposers for qubits and control electronics, and vertical heat extraction routes that minimize interface/boundary resistances (Bon-Mardion et al., 7 May 2026).
• Thermal budgets: Measured conductance values show that per-chip heating must be limited to $^{60}$810–50 nW at 0.3 K to keep ΔT within manageable bounds, constraining the integration of dissipative cryo-CMOS and reinforcing the trend toward functionally segregated architectures (Bon-Mardion et al., 7 May 2026).
• Polymers and detector design: Amorphous polymers such as PTFE and polystyrene represent limiting cases for low thermal conductivity, important for applications where thermal isolation is essential (e.g., support/structural elements for low-background detectors) (Drobizhev et al., 2016, Erhart et al., 2023).
• Superconducting wiring: Thin Nb can be used to bridge islands electrically while providing strong thermal isolation due to its steep (T^{4.5}) suppression of G even above 0.3 K; relevant for constructing calorimeters, bolometric detectors, or nanocircuit thermal switches (Feshchenko et al., 2016).

6. Summary Table: Selected Sub-Kelvin Thermal Conductivity Power Laws

Material/System	κ(T) [W/(m·K)]	Temperature Range (K)	Reference
PTFE (Teflon)	$^{60}$9	0.17–0.43	(Drobizhev et al., 2016)
Al₂O₃ (Sample D)	$5\times10^{-2}$0	0.10–0.72	(Drobizhev et al., 2016)
UPS-923A (Polystyrene)	$5\times10^{-2}$1	0.10–0.70	(Erhart et al., 2023)
HR Si (substrate)	$5\times10^{-2}$2	0.1–1.0	(Bon-Mardion et al., 7 May 2026)
LR Si (substrate)	$5\times10^{-2}$3	0.1–1.0	(Bon-Mardion et al., 7 May 2026)
Borosilicate glass	$5\times10^{-2}$4	0.1–1.0	(Bon-Mardion et al., 7 May 2026)
Sapphire	$5\times10^{-2}$5	0.1–1.0	(Bon-Mardion et al., 7 May 2026)
Nb thin film (5 μm)	$5\times10^{-2}$6	0.1–0.6	(Feshchenko et al., 2016)

These experimentally determined κ(T) dependences are critical inputs for precise thermal simulation and cryogenic engineering in quantum computing, detector design, and low-temperature physics.

7. Outlook and Open Issues

The sub-Kelvin regime continues to expose deviations from semiclassical models. Anomalously high κ in superconductors above BCS predictions, pronounced topology-driven Hall effects, and dramatic material-to-material variation in phonon mean free path all reflect subtle scattering and excitation processes unique to low temperatures. Accurate characterization of G(T) in complex substrates, interfaces, nanostructures, and hybrid systems remains an active area, informing both device design and the fundamental understanding of emergent quantum phenomena in condensed matter at the lowest temperatures accessible in the laboratory.

---

References:

(Feshchenko et al., 2016, Drobizhev et al., 2016, Draelos et al., 2018, Hirschberger et al., 2019, Erhart et al., 2023, Bon-Mardion et al., 7 May 2026)