Skip to content

Mode Decompositions

A mode decomposition is a choice of independent one-particle waveforms, orbitals, field patterns, or internal degrees of freedom used to describe a system by mode occupations. It turns the abstract question “what are the subsystems?” into a concrete choice of modes.

Modes are not usually particles. They are basis elements or independently addressable degrees of freedom. A single particle can be in a superposition of modes, many bosons can occupy the same mode, and a fermionic mode can be either empty or occupied. Entanglement between modes is therefore entanglement relative to a chosen mode decomposition.

This page explains the main kinds of modes used in continuous-variable and field-like systems: spatial modes, momentum modes, frequency modes, polarization modes, and mode bases related by unitary transformations. The detailed algebra of creation and annihilation operators lives in the mode occupations and mode expansions pages.

Start with a one-particle Hilbert space h\mathcal h. A discrete mode decomposition is often an orthonormal basis

{∣φi⟩}i∈I,⟨φi∣φj⟩=δij.\{\lvert\varphi_i\rangle\}_{i\in I}, \qquad \langle\varphi_i\vert\varphi_j\rangle = \delta_{ij}.

Each basis vector defines a mode. In a bosonic Fock space, the occupation basis is built from states

∣n1,n2,…⟩,\lvert n_1,n_2,\ldots\rangle,

where nin_i is the number of excitations in mode φi\varphi_i. For fermions, nin_i can only be 00 or 11, and a fixed ordering of modes is needed to define signs.

The word “mode” can also refer to a normalized wavepacket, a cavity field pattern, a lattice orbital, a spin-orbital, a polarization component, a frequency bin, or a waveguide channel. The common feature is that there is an operator that creates or annihilates an excitation in that degree of freedom.

If ∣f⟩\lvert f\rangle is a normalized one-particle state, the corresponding bosonic creation operator may be written

af†=∑i⟨φi∣f⟩ai†a_f^\dagger = \sum_i \langle\varphi_i\vert f\rangle a_i^\dagger

in a discrete mode basis. The state af†∣0⟩a_f^\dagger\lvert0\rangle is one excitation in the wavepacket mode ff.

Spatial modes are one-particle wavefunctions localized in different regions, channels, traps, lattice sites, or beam paths. For example, two localized wavepackets fL(x)f_L(x) and fR(x)f_R(x) may define left and right modes if

∫dx fL∗(x)fR(x)≈0.\int dx\, f_L^*(x)f_R(x) \approx 0.

When the overlap is exactly zero, the corresponding bosonic creation operators obey

[aL,aR†]=0.[a_L,a_R^\dagger]=0.

When the overlap is nonzero,

[af,ag†]=⟨f∣g⟩,[a_{f},a_g^\dagger] = \langle f\vert g\rangle,

so the two wavefunctions are not independent orthonormal modes. One should orthonormalize the modes or keep the overlap explicitly.

Spatial modes are common in double wells, interferometers, waveguides, and lattice models. They are not the same as exact position eigenstates. A physical spatial mode is a normalizable wavepacket or orbital, while ∣x⟩\lvert x\rangle is a generalized eigenket.

Momentum modes diagonalize translation-invariant one-particle Hamiltonians. In a finite box, a convenient orthonormal basis is

φk(x)=1Ωeik⋅x,\varphi_{\mathbf k}(\mathbf x) = \frac{1}{\sqrt\Omega} e^{i\mathbf k\cdot\mathbf x},

where Ω\Omega is the volume and the allowed k\mathbf k values depend on boundary conditions.

The annihilation field can then be expanded as

ψ(x)=1Ω∑keik⋅xdk.\psi(\mathbf x) = \frac{1}{\sqrt\Omega} \sum_{\mathbf k} e^{i\mathbf k\cdot\mathbf x} d_{\mathbf k}.

Momentum modes are natural for free particles, weakly interacting gases, phonons, photons in homogeneous media, and scattering calculations. Interactions or boundaries can mix them, so a definite momentum occupation is not always conserved.

In infinite volume, sums become integrals and Kronecker deltas become Dirac deltas. The precise powers of 2π2\pi depend on the Fourier convention, so the normalization must be stated.

Frequency modes appear when the system is an oscillator field or when a signal is decomposed spectrally. A cavity field has discrete resonant modes with frequencies ωn\omega_n; a traveling pulse may be described by a continuum of frequency modes.

For a single oscillator mode,

H=ℏω(a†a+12).H = \hbar\omega \left( a^\dagger a+\frac12 \right).

For several independent modes,

H=∑jℏωj(aj†aj+12),H = \sum_j \hbar\omega_j \left( a_j^\dagger a_j+\frac12 \right),

before interactions or couplings are included. Frequency modes are especially useful in quantum optics, spectroscopy, and input-output descriptions of fields.

A frequency label is not enough by itself. A complete electromagnetic mode also includes spatial structure and polarization. Likewise, a wavepacket with a finite duration cannot have a perfectly sharp frequency; time-frequency mode decompositions always involve a resolution tradeoff.

Polarization is an internal mode label for fields such as light. If a fixed spatial-temporal mode supports two orthogonal polarizations, one may use creation operators

aH†,aV†a_H^\dagger, \qquad a_V^\dagger

for horizontal and vertical polarization, or

aR†,aL†a_R^\dagger, \qquad a_L^\dagger

for right- and left-circular polarization. The polarization basis is a mode basis inside a two-dimensional internal space.

The full mode label may combine several pieces:

mode=(spatial profile,frequency,polarization).\text{mode} = (\text{spatial profile},\text{frequency},\text{polarization}).

This matters because a phrase such as “the photon is horizontally polarized” is incomplete unless the spatial and spectral mode structure is also controlled or irrelevant. Entanglement can occur between polarization modes, path modes, frequency modes, or combinations of them.

Two complete orthonormal mode bases of the same one-particle Hilbert space are related by a unitary transformation. If

∣χα⟩=∑iUiα∣φi⟩,\lvert\chi_\alpha\rangle = \sum_i U_{i\alpha}\lvert\varphi_i\rangle,

then the corresponding creation operators transform as

bα†=∑iUiαai†.b_\alpha^\dagger = \sum_i U_{i\alpha}a_i^\dagger.

Unitarity of UU preserves the canonical commutation relations for bosons:

[bα,bβ†]=δαβ.[b_\alpha,b_\beta^\dagger] = \delta_{\alpha\beta}.

The same unitary change preserves fermionic anticommutation relations, but fermionic many-mode states still require a consistent mode ordering.

A two-mode example makes the point visible. Define

b+†=aA†+aB†2,b−†=aA†−aB†2.b_+^\dagger = \frac{a_A^\dagger+a_B^\dagger}{\sqrt2}, \qquad b_-^\dagger = \frac{a_A^\dagger-a_B^\dagger}{\sqrt2}.

The one-excitation state in the ++ mode is

b+†∣0⟩=12(∣1A,0B⟩+∣0A,1B⟩).b_+^\dagger\lvert0\rangle = \frac{1}{\sqrt2} \left( \lvert1_A,0_B\rangle + \lvert0_A,1_B\rangle \right).

With respect to the A,BA,B mode split, this is a single excitation delocalized over two modes. With respect to the +,−+,- mode split, it is simply

∣1+,0−⟩.\lvert1_+,0_-\rangle.

The state has not physically changed. The decomposition used to ask the entanglement question has changed.

Mode entanglement is precise only after the mode decomposition and accessible operations have been specified. Several caveats are standard:

  • A mode basis change can change whether a state looks product or entangled.
  • Superselection rules may restrict which mode superpositions are operationally accessible.
  • Fermionic mode entanglement requires parity and sign conventions.
  • Spatial-region decompositions are not always equivalent to particle or mode decompositions.
  • Continuum mode labels require distributions, wavepackets, or finite-volume regularization.

None of these caveats makes mode entanglement meaningless. They say that the physical question must name the modes, the algebra of observables, and the allowed local operations.

Mode decompositions are useful because experiments and Hamiltonians often single out natural modes:

  • cavities select standing-wave spatial and frequency modes;
  • beam splitters select input and output path modes;
  • optical fibers and waveguides select transverse and polarization modes;
  • lattices select site or Bloch modes;
  • traps select oscillator or orbital modes;
  • scattering setups select incoming and outgoing channels.

The best mode decomposition is usually the one in which preparations, measurements, or dynamics take their simplest form. It is a modeling choice constrained by the physical apparatus, not a matter of taste.

  • Treating a mode label as an intrinsic particle identity.
  • Calling nonorthogonal wavepackets independent modes without accounting for their overlap.
  • Forgetting that polarization alone is not a complete optical mode.
  • Assuming momentum modes are natural when boundaries, traps, or interactions strongly mix them.
  • Confusing a basis change in one-particle Hilbert space with a physical operation on a fixed subsystem split.
  • Ignoring fermionic mode ordering when translating between occupation strings.
  • Treating continuum labels as normalizable modes rather than delta-normalized idealizations.
  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, 1971.
  • D. F. Walls and G. J. Milburn, Quantum Optics, 2nd ed., Springer, 2008.
  • M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge University Press, 1997.
  • C. C. Gerry and P. L. Knight, Introductory Quantum Optics, Cambridge University Press, 2005.
  • S. L. Braunstein and P. van Loock, “Quantum information with continuous variables”, Reviews of Modern Physics 77, 513-577, 2005.
  • C. Weedbrook, S. Pirandola, R. Garcia-Patron, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, “Gaussian quantum information”, Reviews of Modern Physics 84, 621-669, 2012.
  1. Commutator under a mode change. Let bα†=∑iUiαai†b_\alpha^\dagger=\sum_i U_{i\alpha}a_i^\dagger, where UU is unitary and [ai,aj†]=δij[a_i,a_j^\dagger]=\delta_{ij}. Show that [bα,bβ†]=δαβ[b_\alpha,b_\beta^\dagger]=\delta_{\alpha\beta}.
Solution

The annihilation operator is

bα=∑iUiα∗ai.b_\alpha = \sum_i U_{i\alpha}^*a_i.

Then

[bα,bβ†]=∑ijUiα∗Ujβ[ai,aj†]=∑iUiα∗Uiβ=(U†U)αβ=δαβ.\begin{aligned} [b_\alpha,b_\beta^\dagger] &= \sum_{ij} U_{i\alpha}^*U_{j\beta} [a_i,a_j^\dagger] \\ &= \sum_i U_{i\alpha}^*U_{i\beta} = (U^\dagger U)_{\alpha\beta} = \delta_{\alpha\beta}. \end{aligned}
  1. One excitation in two bases. Using b+†=(aA†+aB†)/2b_+^\dagger=(a_A^\dagger+a_B^\dagger)/\sqrt2, show that b+†∣0⟩b_+^\dagger\lvert0\rangle equals a superposition of one excitation in modes AA and BB.
Solution

Apply the definition:

b+†∣0⟩=12(aA†∣0⟩+aB†∣0⟩).b_+^\dagger\lvert0\rangle = \frac{1}{\sqrt2} \left( a_A^\dagger\lvert0\rangle + a_B^\dagger\lvert0\rangle \right).

Since

aA†∣0⟩=∣1A,0B⟩,aB†∣0⟩=∣0A,1B⟩,a_A^\dagger\lvert0\rangle=\lvert1_A,0_B\rangle, \qquad a_B^\dagger\lvert0\rangle=\lvert0_A,1_B\rangle,

one obtains

b+†∣0⟩=12(∣1A,0B⟩+∣0A,1B⟩).b_+^\dagger\lvert0\rangle = \frac{1}{\sqrt2} \left( \lvert1_A,0_B\rangle + \lvert0_A,1_B\rangle \right).
  1. Nonorthogonal wavepackets. Let af†a_f^\dagger and ag†a_g^\dagger create bosons in normalized wavepackets ff and gg. Show that [af,ag†]=⟨f∣g⟩[a_f,a_g^\dagger]=\langle f\vert g\rangle.
Solution

Expand in an orthonormal basis:

af†=∑i⟨φi∣f⟩ai†,ag†=∑j⟨φj∣g⟩aj†.a_f^\dagger = \sum_i \langle\varphi_i\vert f\rangle a_i^\dagger, \qquad a_g^\dagger = \sum_j \langle\varphi_j\vert g\rangle a_j^\dagger.

Then

af=∑i⟨f∣φi⟩ai.a_f = \sum_i \langle f\vert\varphi_i\rangle a_i.

Using [ai,aj†]=δij[a_i,a_j^\dagger]=\delta_{ij},

[af,ag†]=∑ij⟨f∣φi⟩⟨φj∣g⟩[ai,aj†]=∑i⟨f∣φi⟩⟨φi∣g⟩=⟨f∣g⟩.\begin{aligned} [a_f,a_g^\dagger] &= \sum_{ij} \langle f\vert\varphi_i\rangle \langle\varphi_j\vert g\rangle [a_i,a_j^\dagger] \\ &= \sum_i \langle f\vert\varphi_i\rangle \langle\varphi_i\vert g\rangle = \langle f\vert g\rangle. \end{aligned}
  1. Complete optical mode labels. Why is “horizontal polarization” usually not a complete mode label for a photon?
Solution

Polarization is only one part of the mode label. A complete optical mode also needs spatial structure and spectral or temporal structure, at least to the accuracy relevant for the experiment. Two photons with horizontal polarization but different wavepackets, frequencies, paths, or transverse profiles occupy different modes.

  1. Mode entanglement caveat. Explain why the state b+†∣0⟩b_+^\dagger\lvert0\rangle can look entangled in the A,BA,B mode basis but product in the +,−+,- mode basis.
Solution

The state is

b+†∣0⟩=∣1+,0−⟩b_+^\dagger\lvert0\rangle = \lvert1_+,0_-\rangle

in the +,−+,- decomposition, so it is a product of the two new mode factors. In the A,BA,B decomposition, the same vector is

12(∣1A,0B⟩+∣0A,1B⟩).\frac{1}{\sqrt2} \left( \lvert1_A,0_B\rangle + \lvert0_A,1_B\rangle \right).

The vector is the same; the tensor-product structure used to ask the entanglement question is different.