Skip to content

Number Operators

A number operator measures the occupation of a mode. For a bosonic mode ii,

Ni=ai†ai.N_i = a_i^\dagger a_i.

For a fermionic mode ii,

Ni=ci†ci.N_i = c_i^\dagger c_i.

In both cases, number states are eigenstates:

Ni∣n1,n2,…⟩=ni∣n1,n2,…⟩.N_i \lvert n_1,n_2,\ldots\rangle = n_i \lvert n_1,n_2,\ldots\rangle.

The same symbol NiN_i is often used for bosons and fermions because the interpretation is the same: it counts how many particles or excitations occupy mode ii. The allowed eigenvalues depend on the statistics.

Creation and annihilation operators raise and lower occupation. The number operator combines them so that the net effect is to return to the same occupation state with a numerical eigenvalue.

For bosons,

ai∣…,ni,…⟩B=ni ∣…,ni−1,…⟩B,a_i \lvert\ldots,n_i,\ldots\rangle_B = \sqrt{n_i}\, \lvert\ldots,n_i-1,\ldots\rangle_B,

and

ai†∣…,ni−1,…⟩B=ni ∣…,ni,…⟩B.a_i^\dagger \lvert\ldots,n_i-1,\ldots\rangle_B = \sqrt{n_i}\, \lvert\ldots,n_i,\ldots\rangle_B.

Therefore

ai†ai∣…,ni,…⟩B=ni∣…,ni,…⟩B.a_i^\dagger a_i \lvert\ldots,n_i,\ldots\rangle_B = n_i \lvert\ldots,n_i,\ldots\rangle_B.

The fermionic formula has the same eigenvalue statement, with ni=0n_i=0 or 11.

For a bosonic mode,

ni=0,1,2,….n_i=0,1,2,\ldots .

The spectrum of Ni=ai†aiN_i=a_i^\dagger a_i is the set of nonnegative integers. On one mode,

N∣n⟩B=n∣n⟩B.N\lvert n\rangle_B = n\lvert n\rangle_B.

For many modes,

Ni∣n1,…,ni,…⟩B=ni∣n1,…,ni,…⟩B.N_i \lvert n_1,\ldots,n_i,\ldots\rangle_B = n_i \lvert n_1,\ldots,n_i,\ldots\rangle_B.

Since bosons can repeatedly occupy the same mode, NiN_i is not a projector. For example,

Ni2∣…,ni,…⟩B=ni2∣…,ni,…⟩B,N_i^2 \lvert\ldots,n_i,\ldots\rangle_B = n_i^2 \lvert\ldots,n_i,\ldots\rangle_B,

which is generally not the same as NiN_i unless ni=0n_i=0 or 11.

For a fermionic mode,

ni∈{0,1}.n_i\in\{0,1\}.

Thus

Ni∣…,0i,…⟩F=0,N_i\lvert\ldots,0_i,\ldots\rangle_F = 0,

and

Ni∣…,1i,…⟩F=∣…,1i,…⟩F.N_i\lvert\ldots,1_i,\ldots\rangle_F = \lvert\ldots,1_i,\ldots\rangle_F.

The fermionic number operator is a projector:

Ni2=Ni.N_i^2=N_i.

Using the anticommutation relation {ci,ci†}=I\{c_i,c_i^\dagger\}=I,

Ni2=ci†cici†ci=ci†(I−ci†ci)ci=ci†ci−ci†ci†cici=Ni.\begin{aligned} N_i^2 &= c_i^\dagger c_i c_i^\dagger c_i\\ &= c_i^\dagger(I-c_i^\dagger c_i)c_i\\ &= c_i^\dagger c_i - c_i^\dagger c_i^\dagger c_i c_i\\ &= N_i. \end{aligned}

The last term vanishes because (ci†)2=0(c_i^\dagger)^2=0 and ci2=0c_i^2=0.

The total number operator is the sum over modes:

Ntot=∑iNi.N_{\mathrm{tot}} = \sum_i N_i.

On a number state,

Ntot∣n1,n2,…⟩=(∑ini)∣n1,n2,…⟩.N_{\mathrm{tot}} \lvert n_1,n_2,\ldots\rangle = \left(\sum_i n_i\right) \lvert n_1,n_2,\ldots\rangle.

For fixed-particle-number sectors, this eigenvalue is constant. For a general Fock-space state

∣Ψ⟩=∑nCn∣n⟩,\lvert\Psi\rangle = \sum_{\mathbf n} C_{\mathbf n}\lvert\mathbf n\rangle,

the state has definite total number only if all basis states with nonzero CnC_{\mathbf n} have the same value of ∑ini\sum_i n_i.

Commutators with Creation and Annihilation

Section titled “Commutators with Creation and Annihilation”

For bosonic modes,

[Ni,aj†]=δijaj†,[Ni,aj]=−δijaj.[N_i,a_j^\dagger] = \delta_{ij}a_j^\dagger, \qquad [N_i,a_j] = -\delta_{ij}a_j.

For fermionic modes, the same commutators with the number operator hold:

[Ni,cj†]=δijcj†,[Ni,cj]=−δijcj.[N_i,c_j^\dagger] = \delta_{ij}c_j^\dagger, \qquad [N_i,c_j] = -\delta_{ij}c_j.

The algebra used to prove the formulas differs for bosons and fermions, but the meaning is identical:

  • creating in mode jj raises NjN_j by one;
  • annihilating in mode jj lowers NjN_j by one;
  • occupations of other modes are unchanged by that single-mode operation.

For the total number operator,

[Ntot,aj†]=aj†,[Ntot,aj]=−aj,[N_{\mathrm{tot}},a_j^\dagger] = a_j^\dagger, \qquad [N_{\mathrm{tot}},a_j] = -a_j,

and similarly for fermionic cj†,cjc_j^\dagger,c_j.

Let did_i denote either bosonic or fermionic annihilation operators when only number-counting identities are being discussed. A bilinear operator

di†djd_i^\dagger d_j

moves one quantum from mode jj to mode ii. It changes individual mode occupations but preserves total number:

[Ntot,di†dj]=0.[N_{\mathrm{tot}},d_i^\dagger d_j]=0.

For a particular mode kk,

[Nk,di†dj]=(δki−δkj)di†dj.[N_k,d_i^\dagger d_j] = (\delta_{ki}-\delta_{kj})d_i^\dagger d_j.

This says exactly what the bilinear does: it raises the occupation of mode ii and lowers the occupation of mode jj.

The one-body Hamiltonian

H0=∑iϵiNiH_0 = \sum_i \epsilon_i N_i

is diagonal in the occupation basis. More general number-conserving one-body operators have the form

H1=∑ijhijdi†dj.H_1 = \sum_{ij} h_{ij} d_i^\dagger d_j.

They can mix modes while preserving total particle number.

If a Hamiltonian has no explicit time dependence, the Heisenberg equation gives

ddtNtot(t)=iℏ[H,Ntot](t).\frac{d}{dt}N_{\mathrm{tot}}(t) = \frac{i}{\hbar} [H,N_{\mathrm{tot}}](t).

Thus total particle number is conserved when

[H,Ntot]=0.[H,N_{\mathrm{tot}}]=0.

Number-conserving many-particle Hamiltonians do not connect sectors with different total occupation. If a state begins in the N=3N=3 sector, it stays in that sector.

Terms such as

di†dj†ordidjd_i^\dagger d_j^\dagger \quad \text{or} \quad d_i d_j

change total number by two and do not commute with NtotN_{\mathrm{tot}}. They appear in effective descriptions such as pairing Hamiltonians or driven bosonic systems, but then particle number is not conserved in the simple mode-counting sense.

The operator NiN_i counts occupation of mode ii. If the one-particle mode basis is changed, the meaning of NiN_i changes too. A state with definite occupation in one basis can be a superposition of number states in another basis.

The total number operator is basis-independent under unitary changes of the one-particle mode basis, but individual mode occupations are not. This is why one must state the modes before assigning physical meaning to a number operator.

In a position-space basis, the same counting idea becomes the local density ψ†(x)ψ(x)\psi^\dagger(\mathbf x)\psi(\mathbf x) introduced in Field Operators; the basis expansion connecting NiN_i and field notation is treated in Mode Expansions.

  • Treating NiN_i as a particle label rather than a mode-occupation operator.
  • Forgetting that bosonic NiN_i can have eigenvalues 2,3,…2,3,\ldots.
  • Forgetting that fermionic NiN_i is a projector with eigenvalues 00 and 11.
  • Assuming mode occupation is conserved just because total number is conserved.
  • Writing a Hamiltonian with pair-creation terms and still claiming particle number is conserved.
  • Forgetting that individual number operators depend on the chosen mode basis.
  • P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press, 1958.
  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, 1971.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014.
  • A. Altland and B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press, 2010.
  1. Compute Ntot∣2,0,1⟩BN_{\mathrm{tot}}\lvert2,0,1\rangle_B.
Solution

The total occupation is 2+0+1=32+0+1=3, so

Ntot∣2,0,1⟩B=3∣2,0,1⟩B.N_{\mathrm{tot}}\lvert2,0,1\rangle_B = 3\lvert2,0,1\rangle_B.
  1. Use the bosonic action on number states to show that N∣n⟩B=n∣n⟩BN\lvert n\rangle_B=n\lvert n\rangle_B.
Solution

For one bosonic mode,

a∣n⟩B=n ∣n−1⟩B.a\lvert n\rangle_B = \sqrt n\,\lvert n-1\rangle_B.

Then

a†a∣n⟩B=a†n ∣n−1⟩B=nn ∣n⟩B=n∣n⟩B.a^\dagger a\lvert n\rangle_B = a^\dagger\sqrt n\,\lvert n-1\rangle_B = \sqrt n\sqrt n\,\lvert n\rangle_B = n\lvert n\rangle_B.
  1. Prove that a fermionic number operator is a projector.
Solution

Let Ni=ci†ciN_i=c_i^\dagger c_i. Then

Ni2=ci†cici†ci=ci†(I−ci†ci)ci=ci†ci−ci†ci†cici=Ni.\begin{aligned} N_i^2 &= c_i^\dagger c_i c_i^\dagger c_i\\ &= c_i^\dagger(I-c_i^\dagger c_i)c_i\\ &= c_i^\dagger c_i - c_i^\dagger c_i^\dagger c_i c_i\\ &= N_i. \end{aligned}

The last term vanishes because (ci†)2=0(c_i^\dagger)^2=0.

  1. Show that di†djd_i^\dagger d_j preserves total number.
Solution

Use

[Ntot,di†]=di†,[Ntot,dj]=−dj.[N_{\mathrm{tot}},d_i^\dagger]=d_i^\dagger, \qquad [N_{\mathrm{tot}},d_j]=-d_j.

Then

[Ntot,di†dj]=[Ntot,di†]dj+di†[Ntot,dj]=di†dj−di†dj=0.\begin{aligned} [N_{\mathrm{tot}},d_i^\dagger d_j] &= [N_{\mathrm{tot}},d_i^\dagger]d_j + d_i^\dagger[N_{\mathrm{tot}},d_j]\\ &= d_i^\dagger d_j - d_i^\dagger d_j\\ &= 0. \end{aligned}
  1. Does the term di†dj†d_i^\dagger d_j^\dagger conserve total number?
Solution

No. Using the same commutator rule,

[Ntot,di†dj†]=[Ntot,di†]dj†+di†[Ntot,dj†]=di†dj†+di†dj†=2di†dj†.\begin{aligned} [N_{\mathrm{tot}},d_i^\dagger d_j^\dagger] &= [N_{\mathrm{tot}},d_i^\dagger]d_j^\dagger + d_i^\dagger[N_{\mathrm{tot}},d_j^\dagger]\\ &= d_i^\dagger d_j^\dagger + d_i^\dagger d_j^\dagger\\ &= 2d_i^\dagger d_j^\dagger. \end{aligned}

The operator creates two quanta, so it changes total number by two.