Skip to content

One-Body Operators

A one-body operator is an observable or Hamiltonian term obtained by applying the same single-particle operator to each particle. If AA is an operator on the one-particle Hilbert space h\mathcal h, then its fixed-NN action is

A^(N)=∑α=1NA(α),\widehat A^{(N)} = \sum_{\alpha=1}^N A^{(\alpha)},

where A(α)A^{(\alpha)} acts as AA on the α\alphath tensor slot and as the identity on the others. For identical particles this sum is symmetric under relabeling of the formal slots, so it maps the bosonic and fermionic subspaces into themselves.

The companion many-body application guide develops reduced-density-matrix expectations, transition rules, dynamics, truncation, and effective-operator cautions.

In occupation-number language, the same operator has the compact Fock-space form

A^=∑ijAij di†dj,Aij=⟨φi∣A∣φj⟩.\widehat A = \sum_{ij} A_{ij}\,d_i^\dagger d_j, \qquad A_{ij} = \langle\varphi_i\vert A\vert\varphi_j\rangle.

Here {∣φi⟩}\{\lvert\varphi_i\rangle\} is an orthonormal one-particle mode basis, and di,di†d_i,d_i^\dagger denote either bosonic mode operators or fermionic mode operators with the appropriate algebra. The bilinear di†djd_i^\dagger d_j removes one particle from mode jj and creates one particle in mode ii, so the net particle number is unchanged.

For distinguishable particles, the tensor-slot notation is literal. For three particles, for example,

A^(3)=A⊗I⊗I+I⊗A⊗I+I⊗I⊗A.\widehat A^{(3)} = A\otimes I\otimes I + I\otimes A\otimes I + I\otimes I\otimes A.

For identical particles, the slots are bookkeeping devices rather than physical labels. The same sum is nevertheless meaningful because it treats every slot identically. This is the standard way to lift a single-particle observable to the NN-particle sector.

If A=hA=h is a one-particle Hamiltonian, then h^(N)\widehat h^{(N)} is the part of a many-particle Hamiltonian in which each particle moves independently in the same external one-particle operator. Interactions between pairs of particles are not one-body terms; they require two annihilation operators and two creation operators.

Choose an orthonormal basis {∣φi⟩}\{\lvert\varphi_i\rangle\} of h\mathcal h. The one-particle operator can be written as

A=∑ijAij∣φi⟩⟨φj∣,A = \sum_{ij} A_{ij} \lvert\varphi_i\rangle \langle\varphi_j\rvert,

with matrix elements

Aij=⟨φi∣A∣φj⟩.A_{ij} = \langle\varphi_i\vert A\vert\varphi_j\rangle.

If AA is Hermitian, then

Aij=Aji∗.A_{ij} = A_{ji}^*.

Diagonal matrix elements weight occupations. Off-diagonal matrix elements transfer occupation between modes. This distinction is one of the main reasons second quantization is useful: it separates “how many particles are in each mode” from “how the operator mixes modes.”

The Fock-space lift of AA is

A^=∑ijAij di†dj.\widehat A = \sum_{ij} A_{ij}\,d_i^\dagger d_j.

This formula is often denoted dΓ(A)d\Gamma(A) in mathematical treatments. It acts on every particle in a variable-particle-number state while preserving each fixed-NN sector.

On the one-particle sector,

di†dj∣φk⟩=δjk∣φi⟩,d_i^\dagger d_j \lvert\varphi_k\rangle = \delta_{jk} \lvert\varphi_i\rangle,

so

A^∣φk⟩=∑iAik∣φi⟩,\widehat A \lvert\varphi_k\rangle = \sum_i A_{ik} \lvert\varphi_i\rangle,

which is exactly the one-particle action of AA.

On the vacuum,

A^∣0⟩=0,\widehat A\lvert0\rangle=0,

because a one-body operator needs a particle to act on. On an NN-particle sector, the same formula reproduces the first-quantized sum ∑α=1NA(α)\sum_{\alpha=1}^N A^{(\alpha)}.

For bosons and i≠ji\ne j,

ai†aj∣…,ni,…,nj,…⟩B=(ni+1)nj ∣…,ni+1,…,nj−1,…⟩B.a_i^\dagger a_j \lvert\ldots,n_i,\ldots,n_j,\ldots\rangle_B = \sqrt{(n_i+1)n_j}\, \lvert\ldots,n_i+1,\ldots,n_j-1,\ldots\rangle_B.

For i=ji=j the same bilinear is the number operator:

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

For fermions, ci†cjc_i^\dagger c_j moves one fermion from mode jj to mode ii when mode jj is occupied and mode ii is empty. The result is zero if jj is empty or if ii is already occupied. A sign may appear from the chosen ordering of fermionic modes:

ci†cj∣n1,n2,…⟩F=±∣n1′,n2′,…⟩F,c_i^\dagger c_j \lvert n_1,n_2,\ldots\rangle_F = \pm \lvert n_1',n_2',\ldots\rangle_F,

where nj′=0n_j'=0, ni′=1n_i'=1, and all other occupations are unchanged. The sign is not a new physical rule; it is the occupation-number representation of fermionic antisymmetry.

Take

A=∣φr⟩⟨φr∣.A = \lvert\varphi_r\rangle \langle\varphi_r\rvert.

Then Aij=δirδjrA_{ij}=\delta_{ir}\delta_{jr}, so

A^=dr†dr=Nr.\widehat A = d_r^\dagger d_r = N_r.

A mode-occupation operator is therefore a one-body operator associated with a one-particle projector.

If the one-particle Hamiltonian is diagonal,

h∣φi⟩=ϵi∣φi⟩,h\lvert\varphi_i\rangle = \epsilon_i\lvert\varphi_i\rangle,

then

H^1=∑iϵi Ni.\widehat H_1 = \sum_i \epsilon_i\,N_i.

This is the standard free Hamiltonian for independent particles or independent excitations in modes with energies ϵi\epsilon_i.

For two modes with

h=(ϵ1tt∗ϵ2),h = \begin{pmatrix} \epsilon_1 & t \\ t^* & \epsilon_2 \end{pmatrix},

the second-quantized Hamiltonian is

H^1=ϵ1d1†d1+ϵ2d2†d2+t d1†d2+t∗d2†d1.\widehat H_1 = \epsilon_1 d_1^\dagger d_1 +\epsilon_2 d_2^\dagger d_2 +t\,d_1^\dagger d_2 +t^* d_2^\dagger d_1.

The diagonal terms count occupations. The off-diagonal terms transfer occupation between modes.

If the one-particle basis includes spin states ∣s⟩\lvert s\rangle, then a spin component SaS^a is lifted as

S^a=∑ss′⟨s∣Sa∣s′⟩ds†ds′.\widehat S^a = \sum_{ss'} \langle s\vert S^a\vert s'\rangle d_s^\dagger d_{s'}.

If the particle also has orbital modes, the mode label must include both orbital and spin quantum numbers. For example, dμs†d_{\mu s}^\dagger creates a particle in orbital mode μ\mu and spin state ss.

For a continuum position basis, introduce nonrelativistic field operators through the mode expansion:

ψ(x)=∑iφi(x)di,ψ†(x)=∑iφi∗(x)di†.\psi(\mathbf x) = \sum_i \varphi_i(\mathbf x)d_i, \qquad \psi^\dagger(\mathbf x) = \sum_i \varphi_i^*(\mathbf x)d_i^\dagger.

The operator meaning of ψ(x)\psi(\mathbf x) and ψ†(x)\psi^\dagger(\mathbf x) is treated in Field Operators.

If AxA_x is the position-space differential or multiplicative operator representing AA, then

Aij=∫d3x φi∗(x)Axφj(x).A_{ij} = \int d^3x\, \varphi_i^*(\mathbf x) A_x\varphi_j(\mathbf x).

Substituting the mode expansion into the bilinear gives

A^=∫d3x ψ†(x)Axψ(x),\widehat A = \int d^3x\, \psi^\dagger(\mathbf x) A_x \psi(\mathbf x),

with the usual domain and boundary-condition assumptions needed for differential operators.

For a single species of nonrelativistic particles in an external potential,

H^1=∫d3x ψ†(x)(−ℏ22m∇2+U(x))ψ(x).\widehat H_1 = \int d^3x\, \psi^\dagger(\mathbf x) \left( -\frac{\hbar^2}{2m}\nabla^2 +U(\mathbf x) \right) \psi(\mathbf x).

The potential part can also be written using the density operator n(x)=ψ†(x)ψ(x)n(\mathbf x)=\psi^\dagger(\mathbf x)\psi(\mathbf x):

U^=∫d3x U(x)n(x).\widehat U = \int d^3x\, U(\mathbf x)n(\mathbf x).

The formula

A^=∑ijAijdi†dj\widehat A = \sum_{ij} A_{ij}d_i^\dagger d_j

has the same visible shape for bosons and fermions. What changes is the algebra of the dd operators and the allowed occupation numbers.

For bosons,

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

For fermions,

{ci,cj†}=δij,{ci,cj}=0.\{c_i,c_j^\dagger\}=\delta_{ij}, \qquad \{c_i,c_j\}=0.

In both cases, the ordinary commutator with the total number operator vanishes:

[Ntot,A^]=0.[N_{\mathrm{tot}},\widehat A]=0.

This is because every term di†djd_i^\dagger d_j creates and annihilates exactly one particle. A one-body operator can move particles among modes, rotate spin, or mix internal states, but it cannot change total particle number.

The coefficients AijA_{ij} depend on the chosen one-particle basis. The Fock-space operator A^\widehat A does not.

If

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

then the transformed matrix is

A~αβ=∑ijUiα∗AijUjβ,\widetilde A_{\alpha\beta} = \sum_{ij} U_{i\alpha}^* A_{ij}U_{j\beta},

and the transformed annihilation operators satisfy

d~α=∑iUiα∗di.\widetilde d_\alpha = \sum_i U_{i\alpha}^*d_i.

The bilinear form is invariant:

∑αβA~αβd~α†d~β=∑ijAijdi†dj.\sum_{\alpha\beta} \widetilde A_{\alpha\beta} \widetilde d_\alpha^\dagger \widetilde d_\beta = \sum_{ij} A_{ij}d_i^\dagger d_j.

Thus a mode basis is needed to write the formula, but the physical operator is basis-independent.

  • Writing only one tensor-slot term in first quantization instead of summing over all particles.
  • Replacing a non-diagonal one-body operator by ∑iAiiNi\sum_i A_{ii}N_i and losing the mode-mixing terms.
  • Thinking di†djd_i^\dagger d_j changes particle number because it contains a creation operator.
  • Forgetting that fermionic bilinears carry signs through the chosen mode ordering.
  • Treating spin as an extra decoration rather than part of the complete one-particle mode label.
  • Confusing one-body terms with interactions; pair interactions require two-body operators.
  • Assuming the field-operator formula is relativistic. Here ψ(x)\psi(\mathbf x) is a nonrelativistic annihilation field.
  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, McGraw-Hill, 1971.
  • J. W. Negele and H. Orland, Quantum Many-Particle Systems, Addison-Wesley, 1988.
  • A. Altland and B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press, 2010.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014.
  1. Show that A^=∑ijAijdi†dj\widehat A=\sum_{ij}A_{ij}d_i^\dagger d_j reproduces the one-particle action of AA on ∣φk⟩\lvert\varphi_k\rangle.
Solution

Use dj∣φk⟩=δjk∣0⟩d_j\lvert\varphi_k\rangle=\delta_{jk}\lvert0\rangle and di†∣0⟩=∣φi⟩d_i^\dagger\lvert0\rangle=\lvert\varphi_i\rangle. Then

A^∣φk⟩=∑ijAijdi†dj∣φk⟩=∑iAik∣φi⟩.\widehat A\lvert\varphi_k\rangle = \sum_{ij} A_{ij} d_i^\dagger d_j \lvert\varphi_k\rangle = \sum_i A_{ik} \lvert\varphi_i\rangle.

This is the expansion of A∣φk⟩A\lvert\varphi_k\rangle in the same basis.

  1. Let hh be a two-mode one-particle Hamiltonian with h11=ϵ1h_{11}=\epsilon_1, h22=ϵ2h_{22}=\epsilon_2, h12=th_{12}=t, and h21=t∗h_{21}=t^*. Write H^1\widehat H_1.
Solution

Insert the matrix elements into the one-body formula:

H^1=ϵ1d1†d1+ϵ2d2†d2+t d1†d2+t∗d2†d1.\widehat H_1 = \epsilon_1 d_1^\dagger d_1 +\epsilon_2 d_2^\dagger d_2 +t\,d_1^\dagger d_2 +t^*d_2^\dagger d_1.

The first two terms count occupations; the last two move occupation between the two modes.

  1. Compute a1†a2∣1,3⟩Ba_1^\dagger a_2\lvert1,3\rangle_B.
Solution

First annihilate one boson in mode 22, then create one in mode 11:

a1†a2∣1,3⟩B=32 ∣2,2⟩B=6 ∣2,2⟩B.a_1^\dagger a_2 \lvert1,3\rangle_B = \sqrt{3}\sqrt{2}\, \lvert2,2\rangle_B = \sqrt6\, \lvert2,2\rangle_B.
  1. Show that every one-body bilinear di†djd_i^\dagger d_j commutes with NtotN_{\mathrm{tot}}.
Solution

For both bosons and fermions,

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

Therefore

[Ntot,di†dj]=[Ntot,di†]dj+di†[Ntot,dj]=di†dj−di†dj=0.[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.
  1. In position space, write the second-quantized form of a one-particle potential U(x)U(\mathbf x).
Solution

The one-particle operator is multiplication by U(x)U(\mathbf x), so

U^=∫d3x ψ†(x)U(x)ψ(x).\widehat U = \int d^3x\, \psi^\dagger(\mathbf x) U(\mathbf x) \psi(\mathbf x).

Equivalently, using n(x)=ψ†(x)ψ(x)n(\mathbf x)=\psi^\dagger(\mathbf x)\psi(\mathbf x),

U^=∫d3x U(x)n(x).\widehat U = \int d^3x\, U(\mathbf x)n(\mathbf x).