---
title: Kubo-Bastin Formalism in Quantum Transport
url: https://www.emergentmind.com/topics/kubo-bastin-formalism
type: topic
---

# Kubo-Bastin Formalism in Quantum Transport

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 [2005.04678, 1602.04864, 1202.1960]:
\[
\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 [1602.04864, 1202.1960]. 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:
\[
\sigma_{kl} = \sigma^I_{kl} + \sigma^{II}_{kl}
\]
with
\[
\begin{aligned}
\sigma^I_{kl} &= \frac{\hbar}{4\pi} \int d\epsilon\, \partial_\epsilon f(\epsilon)\, \mathrm{Tr}\left\{ [j_k\,G^r\,j_l - j_l\,G^a\,j_k](G^r - G^a) \right\} \\
\sigma^{II}_{kl} &= \frac{\hbar}{4\pi} \int d\epsilon\, f(\epsilon)\, \mathrm{Tr}\left\{ j_k\,G^r\,j_l\,\partial_\epsilon G^r - j_k\,\partial_\epsilon G^r\,j_l\,G^r + j_l\,\partial_\epsilon G^a\,j_k\,G^a - j_l\,G^a\,j_k\,\partial_\epsilon G^a \right\}
\end{aligned}
\]
However, algebraic analysis establishes that both terms contain a common “overlap” term, $\sigma^{\mathrm{ol}}$, such that [2005.04678]:
\[
\sigma^I = \sigma^{\mathrm{surf}} + \sigma^{\mathrm{ol}}, \quad \sigma^{II} = \sigma^{\mathrm{sea}} - \sigma^{\mathrm{ol}}
\]
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 $j_k \leftrightarrow j_l$ exchange. This “symmetrized” form unambiguously projects:
\[
\sigma_{kl} = \sigma^{\mathrm{surf}}_{kl} + \sigma^{\mathrm{sea}}_{kl}
\]
with
\[
\begin{aligned}
\sigma^{\mathrm{surf}}_{kl} &= \frac{\hbar}{4\pi} \int d\epsilon\, \partial_\epsilon f(\epsilon)\, \mathrm{Tr}\left\{ j_k (G^r - G^a) j_l (G^r - G^a) \right\} \\
\sigma^{\mathrm{sea}}_{kl} &= -\frac{\hbar}{4\pi} \int d\epsilon\, f(\epsilon)\, \mathrm{Tr}\left\{ [j_k (\partial_\epsilon G^r + \partial_\epsilon G^a) j_l - j_l (\partial_\epsilon G^r + \partial_\epsilon G^a) j_k] (G^r - G^a) \right\}
\end{aligned}
\]
By construction, these have no overlap [2005.04678, 2408.16611].

## 3. Physical Meaning: Intrinsic and Extrinsic Contributions

The decomposition delineates intrinsic (Berry-curvature) and extrinsic (scattering-driven) mechanisms:
- $\sigma^{\mathrm{surf}}$ (“Fermi-surface”): Proportional to $\partial_\epsilon f$, vanishes in the clean ($\eta \rightarrow 0$) limit for current operators of the form $j_k = -\frac{ie}{\hbar}[G^{-1},x_k]$. This term thus purely represents extrinsic, scattering-induced conductivity, and in the appropriate limit reduces to the Kubo-Greenwood formula for diagonal conductivity.
- $\sigma^{\mathrm{sea}}$ (“Fermi-sea”): Weighted by $f(\epsilon)$, 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, $\sigma^{\mathrm{sea}}$ reduces to
\[
\sigma^{\mathrm{sea}}_{kl} \rightarrow -\frac{e^2}{2\pi\hbar} \int d\epsilon\, f(\epsilon)\, \mathrm{Tr}\{ [G^r - G^a][x_k, x_l] \}
\]
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 $[-1,1]$ and Green's functions and spectral operators are expanded up to order $M$, with a dampling kernel (e.g., Jackson) to suppress Gibbs oscillations. Moments
\[
\mu_{mn}^{\alpha\beta} = \frac{g_m\,g_n}{(1+\delta_{m0})(1+\delta_{n0})}\, \mathrm{Tr}\left[ v_\alpha\,T_m(\tilde H)\,v_\beta\,T_n(\tilde H) \right]
\]
are evaluated stochastically via random-phase vector sampling, ensuring statistical errors scale as $1/\sqrt{SR}$ with the number of vectors $R$ and disorder realizations $S$ [1602.04864]. Disorder averaging and energy integrations are efficiently managed with adaptive quadrature or complex-contour methods [2408.16611].

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 [2408.16611]. This generates an intrinsic level broadening, making the single-particle spectrum continuous without artificial $\eta$.

## 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 $\sigma_{xy} = \nu e^2/h$, with plateaus and vanishing longitudinal conductivity between bands [1202.1960].
- 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 [1602.04864].
- 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 [2005.04678, 2408.16611].

## 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 [2408.16611]. 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 | $\sigma^{II} = 0$           | $\sigma^{\mathrm{sea}} = \sigma^{\mathrm{ol}}$    |
| Lattice Rashba ferromagnet| Strong extrinsic oscillations| $\sigma^{\mathrm{sea}}$ yields Berry peaks        |
| Transition-metal bilayer  | Overlap in $\sigma^I,\sigma^{II}$| $\sigma^{\mathrm{sea}}$ recovers smooth Berry signal|
| Kagome antiferromagnet    | In-plane vs out-plane split  | $\sigma^{\mathrm{surf}}$ (extrinsic), $\sigma^{\mathrm{sea}}$ (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 [2005.04678].

## 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 [2005.04678, 2408.16611].

Source: https://www.emergentmind.com/topics/kubo-bastin-formalism