Kubo-Bastin Formalism in Quantum Transport

• Kubo-Bastin formalism is a rigorous analytic and computational approach for determining linear response transport coefficients in non-interacting quantum systems via single-particle Green’s functions.
• It decomposes transport responses into intrinsic (Berry-curvature) and extrinsic (scattering-driven) contributions, clarifying physical mechanisms behind charge and spin conductivities.
• Widely applied in quantum Hall systems, graphene, and spin-orbit coupled materials, it leverages numerical methods like Chebyshev polynomial expansion for large-scale simulations.

The Kubo-Bastin formalism is a foundational analytic and computational approach for evaluating linear response transport coefficients in non-interacting quantum systems. It rigorously encodes the response of an observable—such as (spin) current density—to an external perturbation, under thermal equilibrium, by expressing the response coefficients in terms of single-particle Green's functions and operator traces. The framework generalizes the more elementary Kubo formula and is especially instrumental for computing both intrinsic (Berry-curvature) and extrinsic (scattering-driven) transport contributions. It provides the basis for quantitative calculations of conductivity tensors, spin Hall responses, and related nonequilibrium observables in a wide class of electronic materials and nanostructures.

1. General Structure of the Kubo-Bastin Formula

The Kubo-Bastin expression for a linear response coefficient $\sigma_{kl}$ relates the induced observable $B$ due to an external perturbation $A$ in a non-interacting system with Hamiltonian $\mathcal H_0$. The retarded/advanced Green's functions are defined as $$G^{r(a)}(\epsilon) = [\epsilon - \mathcal H_0 \pm i\eta]^{-1}$$ in the $\eta\rightarrow 0^+$ limit. The central formula reads (Bonbien et al., 2020, Garcia et al., 2016, Dutta et al., 2012): $$\sigma_{kl} = -\frac{\hbar}{2\pi} \int d\epsilon\, f(\epsilon)\, \mathrm{Tr}\left\{ \left[ j_k\,\partial_\epsilon G^r\, j_l - j_l\,\partial_\epsilon G^a\, j_k \right] (G^r - G^a) \right\}$$ where $j_k,\,j_l$ are current operators, $f(\epsilon)$ is the Fermi-Dirac distribution, and $\mathrm{Tr}$ denotes the trace over single-particle states.

Alternative but equivalent forms, including those using spectral operators and velocity matrix elements, can be derived for model-specific implementations (Garcia et al., 2016, Dutta et al., 2012). The formalism applies for computing both charge and spin conductivities, and, via operator substitutions, any observable linear in the applied field.

2. Decompositions: Smrčka-Středa, Overlap, and Permutation Symmetry

The standard route to physical interpretation and practical separation of contributions involves decomposing the Kubo-Bastin formula. The Smrčka-Středa method splits the integral by parts into so-called “Fermi-surface” and “Fermi-sea” terms: $B$0 with

$B$1

However, algebraic analysis establishes that both terms contain a common “overlap” term, $B$2, such that (Bonbien et al., 2020): $B$3 leading to ambiguities and misattribution between intrinsic and extrinsic transport parts, especially near flat bands in multiband systems.

To resolve this, permutation-based decomposition exploits the symmetry or antisymmetry of terms under $B$4 exchange. This “symmetrized” form unambiguously projects: $B$5 with

$B$6

By construction, these have no overlap (Bonbien et al., 2020, Joao et al., 2024).

3. Physical Meaning: Intrinsic and Extrinsic Contributions

The decomposition delineates intrinsic (Berry-curvature) and extrinsic (scattering-driven) mechanisms:

• $B$7 (“Fermi-surface”): Proportional to $B$8, vanishes in the clean ($B$9) limit for current operators of the form $A$0. This term thus purely represents extrinsic, scattering-induced conductivity, and in the appropriate limit reduces to the Kubo-Greenwood formula for diagonal conductivity.
• $A$1 (“Fermi-sea”): Weighted by $A$2, survives in the clean limit, and encodes geometric/Berry-curvature effects. It is the rigorous quantification of intrinsic transport, fully determined by the band geometry and independent of the disorder or scattering rate.

Explicitly, in the clean limit, $A$3 reduces to

$A$4

recovering the modern theory of the anomalous/spin Hall effect and the geometric torque.

4. Numerical Methods and Computational Schemes

Chebyshev polynomial expansion provides a tractable approach to implementing the Kubo-Bastin formula for large-scale systems. Hamiltonians are rescaled to $A$5 and Green's functions and spectral operators are expanded up to order $A$6, with a dampling kernel (e.g., Jackson) to suppress Gibbs oscillations. Moments

$A$7

are evaluated stochastically via random-phase vector sampling, ensuring statistical errors scale as $A$8 with the number of vectors $A$9 and disorder realizations $\mathcal H_0$0 (Garcia et al., 2016). Disorder averaging and energy integrations are efficiently managed with adaptive quadrature or complex-contour methods (Joao et al., 2024).

For systems with open boundaries or contacts (i.e., two-terminal geometries), self-energy terms from semi-infinite leads are included in the Green's functions (Joao et al., 2024). This generates an intrinsic level broadening, making the single-particle spectrum continuous without artificial $\mathcal H_0$1.

5. Practical Implementations and Model Systems

The Kubo-Bastin formalism is widely used for quantifying transport phenomena in models of quantum Hall systems, spin Hall insulators, and heterostructures:

• In integer quantum Hall systems, the formalism recovers quantization of the Hall conductance as $\mathcal H_0$2, with plateaus and vanishing longitudinal conductivity between bands (Dutta et al., 2012).
• In graphene with random spin-orbit impurities, the Chebyshev-Kubo-Bastin approach captures spin Hall and longitudinal conductivities as functions of impurity concentrations and SOC strength, and clarifies how intrinsic SOC gaps and Rashba SOC generate quantized and sign-changing Hall plateaux (Garcia et al., 2016).
• In proximity-induced graphene with spin-orbit or exchange fields, the overlap-free permutation decomposition cleanly separates Berry-driven and scattering-driven currents and torques, resolving ambiguities in prior approaches (Bonbien et al., 2020, Joao et al., 2024).

6. Relationship to Other Linear Response Frameworks

Comparison with the Keldysh formalism, especially in two-terminal device settings, reveals a numerically exact equivalence for any observable when formulated with the same system partitioning and consideration of both Fermi-surface and Fermi-sea terms (Joao et al., 2024). The Kubo-Bastin sea term in the Keldysh approach must capture the nonlocal voltage drop across the device to reproduce gauge invariance and maintain physical consistency in nonequilibrium responses. This resolution reconciles previously noted discrepancies in the evaluation of occupation-dependent observables and highlights the necessity of correct spectrally resolved integration schemes.

7. Illustrative Examples and Applications

Model Physical System	Clean Limit Signature	Permutation Decomposition Outcome
Two-band Dirac/Rashba gas	$\mathcal H_0$3	$\mathcal H_0$4
Lattice Rashba ferromagnet	Strong extrinsic oscillations	$\mathcal H_0$5 yields Berry peaks
Transition-metal bilayer	Overlap in $\mathcal H_0$6	$\mathcal H_0$7 recovers smooth Berry signal
Kagome antiferromagnet	In-plane vs out-plane split	$\mathcal H_0$8 (extrinsic), $\mathcal H_0$9 (intrinsic)

In all these contexts, the permutation (symmetrized) decomposition isolates the intrinsic (Berry curvature) contributions and eliminates false attribution of extrinsic features (e.g. from trivial flat bands) to geometric responses (Bonbien et al., 2020).

8. Significance and Resolution of Ambiguities

The permutation-based Kubo-Bastin decomposition establishes a rigorous, overlap-free separation between intrinsic (band-geometry/Berry-curvature) and extrinsic (scattering or disorder-driven) contributions to any linear transport coefficient in solids. Numerical equivalence with the Keldysh approach, when both are implemented with physically grounded device partitioning and open boundary conditions, demonstrates the physical completeness and universality of the formalism. This framework underpins modern theoretical and computational descriptions of quantum transport phenomena, anomalous and spin Hall effects, and nonequilibrium spin-orbit torques across diverse material systems (Bonbien et al., 2020, Joao et al., 2024).