Quantum Thermal Transport Heat Conduction Nanostructures Landauer Büttiker Phonon Green Function

# Quantum Thermal Transport: Nanoscale Phonon Conduction, Landauer-Büttiker Formalism, and Quantum Heat Conductance

## 1. Introduction to Quantum Thermal Transport

Heat conduction at nanoscale dimensions—from atomic contacts to nanowires and phononic nanostructures—exhibits fundamentally quantum behavior distinct from classical Fourier diffusion. When thermal wavelengths λ_q = h/√(2π m k_B T) become comparable to device dimensions, phonons behave as quantum particles with discrete energy levels and coherent transport properties. The result: a quantum of thermal conductance G_q = π²k_B²T/(3h) ≈ 9.456 × 10⁻¹³ T W/K (W/K) that is universal—independent of material, geometry, or phonon statistics. This section develops the Landauer-Büttiker formalism for phononic transport, derives the quantum limit, and explores the crossover from quantum ballistic to semi-classical diffusive regimes. Numerical solvers reveal transmission functions, temperature-dependent conductance curves, and the role of mode counting and scattering in nanoscale heat transport.

## 2. Phonons as Quantum Particles and Thermal Quanta

In quantum mechanics, lattice vibrations (phonons) are quantized bosons with energy-frequency relation ħω(q). For acoustic phonons in 3D crystals, the dispersion is typically linear at low frequencies (Debye model):

where v_s is the sound velocity (typically 3,000–5,000 m/s). The number of phonons in thermal equilibrium at temperature T is given by the Bose-Einstein distribution:

The thermal energy is:

where g(ω) = V ω² / (2π² v_s³) is the phonon density of states (for 3D isotropic continuum). At high temperatures (k_B T >> ħω_max), the Rayleigh-Jeans limit applies: E_thermal → k_B T. At low temperatures (T → 0), quantum effects dominate: phonon excitations freeze out except for zero-point motion ħω/2.

## 3. Ballistic Thermal Transport and the Landauer-Büttiker Formula

When heat-carrying phonons traverse a device (a nanowire, atomic contact, or phononic interface) without scattering—ballistic transport—the thermal current is determined by the phonon transmission probability. The Landauer-Büttiker thermal current is:

where:
- T_ph(ω) is the phonon transmission function (0 ≤ 𝒯 ≤ number of propagating modes)
- n_L, n_R are left and right thermal populations
- The factor 1/(2π) emerges from the density of states in 1D

For a device with N(ω) propagating phonon modes at frequency ω, the transmission sums over all modes:

In the quantum limit where only N_mode = 1 phonon mode propagates and T = 1 (perfect transmission), the thermal conductance is:

This remarkable result—universal quantum of heat conductance—depends only on fundamental constants and temperature, not on material properties. Numerically: G_q/(2π) ≈ 1.5 × 10⁻¹³ T W/K (with T in Kelvin).

## 4. Quantum of Thermal Conductance and Universality

The universal G_q emerges from integrating the Bose-Einstein thermal current over all frequencies, assuming one propagating mode:

The derivative of Bose-Einstein distribution is:

Substituting x = ħω/(k_B T):

The integral ∫ x² e^x / (e^x - 1)² dx = π²/3 (related to the Riemann zeta function). This derivation reveals the universality: no material parameters appear. Any system with a single phonon mode at perfect transmission exhibits G_q, whether it's a carbon nanotube, atomic chain, or engineered waveguide.

## 5. Phonon NEGF Formalism and Non-Equilibrium Green's Functions

For more complex geometries where scattering, disorder, or interference effects matter, the non-equilibrium Green's function (NEGF) method provides exact transmission:

where:
- G^R(ω) = [E - H_d - Σ^R]⁻¹ is the retarded Green's function of the device
- Γ_L,R(ω) = i[Σ^R_L,R - (Σ^R_L,R)†] are the broadening matrices from left/right leads
- Σ^R are retarded self-energies representing coupling to semi-infinite leads

The device Hamiltonian H_d represents the scattering region: a disordered region, interface, or defect. The NEGF approach systematically includes all multiple scattering events and interference effects. The thermal current follows from integrating the transmission weighted by phonon populations:

## 6. Transmission Through Atomic Contacts and Nanowires

Acoustic phonon transmission through a single atomic contact—a one-atom-wide "bridge" between two bulk reservoirs—exhibits remarkable quantization. The number of propagating modes at frequency ω is limited by the contact area:

For a 1D contact, N(ω) = 1 (single mode). For 2D and 3D contacts, multiple modes contribute. In the low-frequency limit, N(ω) ∝ ω^(D-1) (where D is dimension). The transmission depends sensitively on atomic-scale geometry:

• Smooth contacts (radius ≈ nanometers): transmission ≈ 1 (ballistic)
• Constricted regions (atomic scale): reflection and scattering reduce transmission by 50-80%
• Rough interfaces: diffuse scattering dominates, T ∝ (λ_q/L)² (geometric suppression)

Measurements of thermal conductance through atomic contacts using suspended graphene or carbon nanotubes confirm the quantum-limited regime: conductance plateaus at multiples of G_q as the number of phonon modes increases with frequency.

## 7. Ballistic-to-Diffusive Crossover

Real nanoscale devices transition between ballistic (Landauer) and diffusive (Fourier) heat transport regimes as sample size or temperature changes. The crossover occurs when the mean free path ℓ_mfp becomes comparable to device length L:

• Knudsen >> 1: Ballistic (Landauer), G = T · N_mode · G_q
• Knudsen ~1: Transition regime with complex interference
• Knudsen << 1: Diffusive (Fourier), G = κ A / L with Fourier thermal conductivity κ

At high temperatures (classical limit), the mean free path scales as ℓ_mfp ∝ 1/T, pushing the crossover toward larger devices. At low temperatures, ℓ_mfp → ∞ for elastic scattering (disorder), so quantum ballistic effects dominate—explaining why cryogenic measurements of thermal transport often reveal quantum signatures.

## 8. Temperature-Dependent Conductance and Phonon Scattering

The thermal conductance G(T) = dJ_q/dT is temperature-dependent due to:

1. Mode population (Bose-Einstein statistics): At T → 0, phonons freeze out, G → 0; at high T, all modes are thermally excited.

2. Mode occupancy in transmission: The thermal current integral weights T(ω) by the thermal derivative ∂n/∂T, which peaks around k_B T / ħ.

3. Frequency-dependent scattering: Phonon lifetimes depend on collision rates (isotope scattering ∝ ω⁴; Umklapp ∝ e^(-θ_D/T)). At low T, Umklapp scattering vanishes, increasing mean free path.

For a perfectly transmitting single-mode ballistic contact, G(T) = G_q (independent of T). For multi-mode or scattering-dominated systems, G(T) exhibits curvature: often G ∝ T^α with α ∈ [1,2] depending on scattering mechanisms.

## 9. Phononic Nanostructures: Phononic Crystals and Waveguides

Engineered phononic structures—phononic crystals (periodic lattices of inclusions), waveguides, and metamaterials—can control thermal transport via phonon band structure engineering. A phononic crystal's dispersion relation ω_n(k) exhibits band gaps: frequency ranges where propagation is forbidden. Heat transport through the gap is exponentially suppressed: T ∝ e^(-κL) (where κ is the decay constant).

Applications:
- Thermal isolation: Phononic band gaps block heat flow while permitting electrical conduction (decoupling heat from electronic transport)
- Thermal diodes: Asymmetric structures exhibiting directional thermal rectification
- Quantum heat batteries: Using phonon blockade to store thermal energy

## 10. Experimental Measurements and Validation

Modern experiments probe quantum thermal transport via:

1. Suspended nanowire thermometry: Measure voltage across a suspended sample heated by electrical current, extract thermal conductance from G = I / ΔT

2. Transient thermal grating: Optical excitation creates periodic temperature modulation, measure phonon propagation via diffraction

3. Time-domain thermoreflectance (TDTR): Ultrafast laser pulses generate thermal transients, measure decay via reflectivity changes

4. Direct thermal microscopy: Scanning probes measure local temperature, infer conductance from spatial profiles

Observations confirm:
- Conductance quantization at G_q in atomic contacts and single-walled carbon nanotubes
- Temperature-dependent G(T) with power-law scaling in diffusive regimes
- Phononic band gap effects in engineered nanostructures
- Ballistic-diffusive crossover at Knudsen number ≈ 1

## 11. Python Numerical Solver: Transmission Functions and Thermal Conductance

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad, odeint
from scipy.special import zeta

# Physical constants (SI units)
HBAR = 1.0546e-34  # J·s
KB = 1.3806e-23    # J/K
H_PLANCK = 6.6261e-34  # J·s

def debye_frequency(sound_velocity, electron_density_3d):
    """
    Debye cutoff frequency ω_D from Debye model.
    ω_D = v_s (6π² n)^(1/3)
    
    Args:
        sound_velocity: m/s
        electron_density_3d: m^-3
    
    Returns:
        omega_D: rad/s
    """
    omega_D = sound_velocity * (6 * np.pi**2 * electron_density_3d)**(1/3)
    return omega_D

def phonon_dos_debye(omega, omega_D):
    """
    Density of states in Debye model.
    g(ω) = 9ω²/ω_D³ for ω < ω_D, else 0
    """
    if omega < omega_D:
        return 9 * omega**2 / omega_D**3
    else:
        return 0

def bose_einstein(omega, T):
    """
    Bose-Einstein distribution n(ω,T) = 1/(exp(ℏω/k_B T) - 1)
    """
    if T < 1e-6:
        return 0
    x = HBAR * omega / (KB * T)
    if x > 100:  # Avoid overflow
        return 0
    return 1.0 / (np.exp(x) - 1)

def thermal_derivative_be(omega, T):
    """
    Derivative of Bose-Einstein: ∂n/∂T
    """
    if T < 1e-6:
        return 0
    x = HBAR * omega / (KB * T)
    if x > 100:
        return 0
    ex = np.exp(x)
    return (HBAR * omega / (KB * T**2)) * ex / (ex - 1)**2

def transmission_single_mode(omega, omega_max):
    """
    Transmission for single phonon mode.
    T(ω) = 1 if ω < ω_max (propagating), else 0 (evanescent)
    """
    return 1.0 if omega < omega_max else 0.0

def transmission_multi_mode(omega, omega_D, n_modes_max=5):
    """
    Transmission for multi-mode case: cumulative from distinct propagating modes.
    T(ω) = min(N(ω), n_modes_max) where N(ω) ∝ ω^(D-1) for D=3D system
    
    Simplified: T(ω) = min(1 + int(ω/ω_D * n_modes_max), n_modes_max)
    """
    mode_fraction = omega / omega_D
    n_active = int(1 + np.floor(mode_fraction * n_modes_max))
    return float(min(n_active, n_modes_max))

def transmission_with_scattering(omega, omega_D, mean_free_path, device_length):
    """
    Transmission with phonon scattering via mean free path.
    T(ω) ∝ exp(-device_length / mean_free_path(ω))
    
    mfp(ω) ~ mfp_0 * (ω_D/ω)^2 (assuming isotope scattering dominates)
    """
    mfp_0 = mean_free_path
    if omega < 1e-6:
        return 1.0
    mfp = mfp_0 * (omega_D / omega)**2
    mfp = np.clip(mfp, 0, 1e-3)  # Clip to physical range
    return np.exp(-device_length / mfp) if mfp > 1e-6 else 0.0

def quantum_heat_conductance(T, n_modes=1):
    """
    Quantum of thermal conductance: G_q = π² k_B² T / (3h)
    
    For n_modes, multiply by n_modes.
    Returns: W/K
    """
    G_q = (np.pi**2 * KB**2 * T) / (3 * H_PLANCK)
    return n_modes * G_q

def thermal_conductance_numerical(transmission_func, T_L, T_R, omega_D, n_omega=500):
    """
    Compute thermal conductance numerically via Landauer-Büttiker formula.
    
    J_q = (1/2π) ∫ ℏω 𝒯(ω) [n_L(ω) - n_R(ω)] dω
    G = dJ_q / dT_L (at T_R = 0)
    """
    omega_max = 3 * omega_D  # Integrate up to 3ω_D (captures ~99% of thermal energy)
    omega_vals = np.linspace(0, omega_max, n_omega)
    domega = omega_vals[1] - omega_vals[0]
    
    J_q = 0
    for omega in omega_vals:
        if omega > 0:
            T = transmission_func(omega)
            T_diff = (bose_einstein(omega, T_L) - bose_einstein(omega, T_R))
            J_q += HBAR * omega * T * T_diff * domega
    
    J_q /= (2 * np.pi)
    
    # Conductance = dJ_q / dT_L (numerical derivative)
    dT = 0.01  # K
    T_L_plus = T_L + dT
    J_q_plus = 0
    for omega in omega_vals:
        if omega > 0:
            T = transmission_func(omega)
            T_diff = (bose_einstein(omega, T_L_plus) - bose_einstein(omega, T_R))
            J_q_plus += HBAR * omega * T * T_diff * domega
    J_q_plus /= (2 * np.pi)
    
    G = (J_q_plus - J_q) / dT
    return G

# Main execution
if __name__ == "__main__":
    print("="*70)
    print("QUANTUM THERMAL TRANSPORT: LANDAUER-BÜTTIKER FORMALISM")
    print("="*70)
    
    # Physical parameters (typical for silicon or GaAs)
    v_s = 4000  # Sound velocity (m/s)
    n_e = 1e28  # Electron density (m^-3)
    omega_D = debye_frequency(v_s, n_e)
    T_D = HBAR * omega_D / KB  # Debye temperature
    
    print(f"\nDebye Model Parameters:")
    print(f"  Sound velocity v_s = {v_s} m/s")
    print(f"  Debye frequency ω_D = {omega_D:.2e} rad/s")
    print(f"  Debye temperature T_D = {T_D:.1f} K")
    print(f"\nQuantum of thermal conductance:")
    
    for T in [1, 10, 100, 300]:
        G_q = quantum_heat_conductance(T, n_modes=1)
        print(f"  G_q({T} K) = {G_q:.2e} W/K")
    
    # Plotting
    fig, axes = plt.subplots(2, 2, figsize=(13, 10))
    
    # Plot 1: Thermal derivative of Bose-Einstein
    T_vals = [10, 100, 300]
    omega_range = np.linspace(0, 3*omega_D, 300)
    
    ax = axes[0, 0]
    for T in T_vals:
        deriv_vals = [thermal_derivative_be(omega, T) for omega in omega_range]
        ax.plot(omega_range / omega_D, deriv_vals, linewidth=2, label=f'T = {T} K')
    ax.set_xlabel('ω / ω_D', fontsize=11)
    ax.set_ylabel('∂n/∂T (1/K)', fontsize=11)
    ax.set_title('Thermal Derivative: Bose-Einstein Distribution', fontsize=12, fontweight='bold')
    ax.legend()
    ax.grid(True, alpha=0.3)
    
    # Plot 2: Transmission functions
    ax = axes[0, 1]
    T_single = [transmission_single_mode(omega, 1.5*omega_D) for omega in omega_range]
    T_multi = [transmission_multi_mode(omega, omega_D, n_modes_max=4) for omega in omega_range]
    ax.plot(omega_range/omega_D, T_single, 'b-', linewidth=2.5, label='Single-mode ballistic')
    ax.plot(omega_range/omega_D, T_multi, 'r-', linewidth=2.5, label='Multi-mode (4 modes)')
    ax.fill_between(omega_range/omega_D, 0, T_single, alpha=0.2, color='blue')
    ax.fill_between(omega_range/omega_D, 0, T_multi, alpha=0.2, color='red')
    ax.set_xlabel('ω / ω_D', fontsize=11)
    ax.set_ylabel('Phonon Transmission 𝒯(ω)', fontsize=11)
    ax.set_title('Ballistic Transmission Functions', fontsize=12, fontweight='bold')
    ax.legend()
    ax.grid(True, alpha=0.3)
    ax.set_ylim([0, 4.5])
    
    # Plot 3: Thermal conductance vs temperature
    T_range = np.linspace(1, 300, 100)
    G_q_vals = [quantum_heat_conductance(T, n_modes=1) for T in T_range]
    G_multi_vals = [quantum_heat_conductance(T, n_modes=3) for T in T_range]
    
    ax = axes[1, 0]
    ax.loglog(T_range, G_q_vals, 'b-', linewidth=2.5, label='Single-mode (n=1)')
    ax.loglog(T_range, G_multi_vals, 'r-', linewidth=2.5, label='Three-mode (n=3)')
    ax.set_xlabel('Temperature T (K)', fontsize=11)
    ax.set_ylabel('Thermal Conductance G_q (W/K)', fontsize=11)
    ax.set_title('Universal Quantum Conductance: $G_q = π²k_B²T / 3h$', fontsize=12, fontweight='bold')
    ax.legend()
    ax.grid(True, alpha=0.3, which='both')
    
    # Plot 4: Ballistic vs diffusive crossover
    ax = axes[1, 1]
    L_device = 100e-9  # Device length 100 nm
    mfp_0 = 1e-6  # Reference mean free path 1 μm
    T_set = 100  # K
    
    T_scatter = [transmission_with_scattering(omega, omega_D, mfp_0, L_device) for omega in omega_range]
    ax.plot(omega_range/omega_D, T_scatter, 'g-', linewidth=2.5)
    ax.fill_between(omega_range/omega_D, 0, T_scatter, alpha=0.3, color='green')
    ax.set_xlabel('ω / ω_D', fontsize=11)
    ax.set_ylabel('Phonon Transmission with Scattering', fontsize=11)
    ax.set_title(f'Transmission with Scattering (L={L_device*1e9:.0f} nm, T={T_set} K)', fontsize=12, fontweight='bold')
    ax.grid(True, alpha=0.3)
    ax.text(0.5, 0.7, f'mfp = {mfp_0*1e6:.1f} μm\nKnudsen ≈ {mfp_0/L_device:.1f}
(Ballistic regime)', 
            transform=ax.transAxes, fontsize=10, bbox=dict(boxstyle='round', facecolor='wheat', alpha=0.5))
    
    plt.tight_layout()
    plt.savefig('quantum_thermal_transport.png', dpi=150, bbox_inches='tight')
    print("
[Plot saved as 'quantum_thermal_transport.png']")
    
    print("
" + "="*70)
    print("THERMAL TRANSPORT ANALYSIS COMPLETE")
    print("="*70)

## References

1. Landauer, R. (1970). "Electrical resistance of disordered one-dimensional lattices." *Philosophical Magazine*, 21(172), 863-867.

2. Büttiker, M. (1988). "Four-terminal phase-coherent conductance." *Physical Review Letters*, 57(14), 1761.

3. Cappelli, E., et al. (2013). "Heat transport in one-dimensional systems." *Reports on Progress in Physics*, 77(5), 056801.

4. Maldovan, M. (2013). "Sound and heat revolutions in phononics." *Nature*, 503(7475), 209-217.

5. Xu, X., et al. (2014). "Thermal transport in nanostructures." *Nature Nanotechnology*, 9(12), 951-963.