Skip to content

Matsubara Frequency Table

Canonical treatment: Bosonic and Fermionic Matsubara Frequencies owns the derivation, transform conventions, routing arguments, numerical discussion, exercises, and references. This page is the quick lookup table.

Let β=1/(kBT)\beta=1/(k_{\mathrm B}T) and let the imaginary-time circumference be Lτ=βℏL_\tau=\beta\hbar.

ChannelBoundary conditionEnergy gridAngular-frequency gridZero modeIndex reflection
BosonicF(τ+Lτ)=F(τ)F(\tau+L_\tau)=F(\tau)Ωm=2πm/β\Omega_m=2\pi m/\betaωm=2πm/(βℏ)\omega_m=2\pi m/(\beta\hbar)m=0m=0m↦−mm\mapsto-m
FermionicF(τ+Lτ)=−F(τ)F(\tau+L_\tau)=-F(\tau)εn=(2n+1)π/β\varepsilon_n=(2n+1)\pi/\betaωn=(2n+1)π/(βℏ)\omega_n=(2n+1)\pi/(\beta\hbar)nonen↦−n−1n\mapsto-n-1

The energy variables have units of energy; the angular frequencies have units of inverse time and obey Ωm=ℏωm\Omega_m=\hbar\omega_m. In natural units ℏ=kB=1\hbar=k_{\mathrm B}=1, the two notations look identical, so the convention must still be declared.

One consistent energy convention is

F(τ)=1β∑ne−iεnτ/ℏF(iεn),F(iεn)=1ℏ∫0βℏdτ eiεnτ/ℏF(τ).F(\tau) = \frac1\beta\sum_n e^{-i\varepsilon_n\tau/\hbar}F(i\varepsilon_n), \qquad F(i\varepsilon_n) = \frac1\hbar \int_0^{\beta\hbar}d\tau\, e^{i\varepsilon_n\tau/\hbar}F(\tau).

If angular frequencies are used instead, the inverse-transform prefactor is 1/(βℏ)1/(\beta\hbar). Do not combine the exponent from one convention with the sum measure from the other.

  • boson ±\pm boson is bosonic;
  • fermion −- fermion is bosonic;
  • fermion ++ boson is fermionic;
  • a composite operator uses the grid fixed by its thermal boundary condition, not simply by the statistics of its microscopic constituents.

Chemical potential normally appears in the propagator denominator, for example iεn+μ−ϵki\varepsilon_n+\mu-\epsilon_{\mathbf k}, rather than by redefining the standard grid. Twisted boundary conditions are an intentional generalization and must state their phase convention.

  • Preserve paired indices: mm with −m-m, and nn with −n−1-n-1.
  • Treat the bosonic zero mode separately when infrared behavior is important.
  • The smallest nonzero bosonic energy is 2π/β2\pi/\beta; the smallest fermionic magnitude is π/β\pi/\beta.
  • A Matsubara cutoff is a numerical regulator, not a physical bandwidth.
  • Analytic continuation relates imaginary-axis data to a specified retarded, advanced, or other real-frequency object; it is not the substitution iεn→ωi\varepsilon_n\to\omega without boundary data.