Skip to content

Core Objects and Notation

Many-body notation compresses large physical structures into a few symbols. That compression is useful only when every index, normalization, ensemble, operator ordering, and limiting convention is recoverable from the page.

This page establishes the default notation for this volume. A specialized article may choose another convention, but it must declare the translation explicitly.

Every calculation should identify:

  1. the retained degrees of freedom;
  2. the Hilbert space or Fock-space sector;
  3. particle, site, mode, spin, and internal indices;
  4. the state or statistical ensemble;
  5. the Hamiltonian and conserved quantities;
  6. boundary conditions and geometry;
  7. Fourier-transform normalization;
  8. operator ordering and time ordering;
  9. units and constants set to one;
  10. the finite-size or thermodynamic limits being taken.

Symbols alone do not supply this information. For example, NN may mean particle number, number of sites, an internal symmetry rank, or the size parameter in a large-NN approximation.

The default visual distinctions are:

ObjectTypical notationComment
scalar numberNN, TT, μ\mu, ttMeaning and units must be stated
spatial vectorr\mathbf r, k\mathbf k, q\mathbf qBoldface distinguishes vectors
abstract operatorHH, NN, SαS^\alphaHats are used when ambiguity remains
state vector$\Psi\rangle$
density operatorρ\rhoPositive, trace one for a normalized state
local sitei,j,ℓi,j,\ellOften integer lattice labels
one-particle modep,q,kp,q,kMay denote momentum or a generic mode
internal componentσ,α,β\sigma,\alpha,\betaSpin, flavor, band, or species
spatial componenta,ba,bCartesian or tensor component
eigenstate labeln,mn,mNot automatically a particle label
imaginary timeτ\tauUsually 0≤τ<βℏ0\le\tau<\beta\hbar

The same Greek letter should not serve simultaneously as a spin index, a coupling, and a critical exponent on one page.

Common size symbols are:

SymbolDefault meaning
NNnumber of particles
NsN_snumber of lattice sites
MMnumber of retained modes or orbitals
VVphysical volume
LLlinear system size
ddspatial dimension
aalattice spacing
n=N/Vn=N/Vparticle-number density
ν=N/Ns\nu=N/N_sfilling per lattice site

For a hypercubic lattice with NsN_s sites,

Ns=(La)dN_s = \left( \frac{L}{a} \right)^d

when every direction has length LL and the endpoints are counted according to a compatible periodic convention.

For an anisotropic system, list L1,…,LdL_1,\ldots,L_d rather than hiding the shape in one LL.

Particle number and site number are independent. The filling

ν=NNs\nu = \frac{N}{N_s}

can be less than, equal to, or greater than one for bosons. For spinless fermions with one orbital per site,

0≤ν≤1,0\le\nu\le1,

while spin or orbital degeneracy changes the maximum occupancy.

A normalized pure state is denoted ∣Ψ⟩|\Psi\rangle, with

⟨Ψ∣Ψ⟩=1.\langle\Psi|\Psi\rangle=1.

A normalized mixed state is a density operator ρ\rho satisfying

ρ≥0,Tr⁡ρ=1.\rho\ge0, \qquad \operatorname{Tr}\rho=1.

Expectation values are

⟨A⟩Ψ=⟨Ψ∣A∣Ψ⟩\langle A\rangle_\Psi = \langle\Psi|A|\Psi\rangle

or

⟨A⟩ρ=Tr⁡(ρA).\langle A\rangle_\rho = \operatorname{Tr}(\rho A).

When the state is clear, the subscript is omitted:

⟨A⟩.\langle A\rangle.

The omission must not conceal an ensemble change.

For a subsystem AA of a joint state ρAB\rho_{AB},

ρA=Tr⁡BρAB.\rho_A = \operatorname{Tr}_B\rho_{AB}.

The label AA may denote a spatial region, a set of sites, a species, or a mode collection. State the partition, especially for identical particles where particle and mode partitions are not interchangeable.

The canonical finite-system definitions live in Reduced States and Partial Trace.

The principal equilibrium ensembles are:

EnsembleFixed macroscopic dataDensity operator
microcanonicalenergy window, particle number, volumenormalized projector onto the allowed shell
canonicaltemperature, particle number, volumee−βH/Ze^{-\beta H}/Z
grand canonicaltemperature, chemical potential, volumee−β(H−μN^)/Ξe^{-\beta(H-\mu\hat N)}/\Xi

The inverse temperature is

β=1kBT.\beta = \frac{1}{k_{\mathrm B}T}.

The canonical partition function is

Z(β)=Tr⁡(e−βH),Z(\beta) = \operatorname{Tr} \left( e^{-\beta H} \right),

and the canonical state is

ρβ=e−βHZ(β).\rho_\beta = \frac{e^{-\beta H}} {Z(\beta)}.

The grand partition function is

Ξ(β,μ)=Tr⁡[e−β(H−μN^)],\Xi(\beta,\mu) = \operatorname{Tr} \left[ e^{-\beta(H-\mu\hat N)} \right],

with state

ρβ,μ=e−β(H−μN^)Ξ(β,μ).\rho_{\beta,\mu} = \frac{ e^{-\beta(H-\mu\hat N)} }{ \Xi(\beta,\mu) }.

Use N^\hat N for the number operator when it could be confused with its expectation or eigenvalue. Once the distinction is established, many-body texts often write NN for both; this volume keeps the hat where ambiguity matters.

Common symbols are:

SymbolQuantityDefining relation
UU or EEinternal energy⟨H⟩\langle H\rangle
SSthermodynamic entropy−kBTr⁡(ρln⁡ρ)-k_{\mathrm B}\operatorname{Tr}(\rho\ln\rho)
FFHelmholtz free energy−kBTln⁡Z-k_{\mathrm B}T\ln Z
Ω\Omegagrand potential−kBTln⁡Ξ-k_{\mathrm B}T\ln\Xi
μ\muchemical potentialconjugate to particle number
PPpressure−(∂F/∂V)T,N-(\partial F/\partial V)_{T,N}
CVC_Vheat capacity(∂U/∂T)V,N(\partial U/\partial T)_{V,N}

The symbol SS may also denote spin, an action, or a structure factor. Context or a subscript should distinguish these meanings.

For a canonical ensemble,

U=−∂∂βln⁡ZU = -\frac{\partial}{\partial\beta} \ln Z

when HH has no implicit β\beta dependence.

For a grand-canonical ensemble,

⟨N^⟩=∂∂(βμ)ln⁡Ξ\langle\hat N\rangle = \frac{\partial}{\partial(\beta\mu)} \ln\Xi

when β\beta is held fixed in that derivative.

Equivalent formulas written as derivatives with respect to μ\mu must carry the factor of β\beta correctly.

Bosonic annihilation and creation operators are denoted

bi,bi†,b_i, \qquad b_i^\dagger,

and satisfy

[bi,bj†]=δij,[bi,bj]=0,[bi†,bj†]=0.\begin{aligned} [b_i,b_j^\dagger] &= \delta_{ij}, \\ [b_i,b_j] &=0, \\ [b_i^\dagger,b_j^\dagger] &=0. \end{aligned}

Fermionic operators are denoted

ci,ci†,c_i, \qquad c_i^\dagger,

and satisfy

{ci,cj†}=δij,{ci,cj}=0,{ci†,cj†}=0.\begin{aligned} \{c_i,c_j^\dagger\} &= \delta_{ij}, \\ \{c_i,c_j\} &=0, \\ \{c_i^\dagger,c_j^\dagger\} &=0. \end{aligned}

The generic symbol aia_i may denote either statistics only after the algebra has been stated.

The mode occupation operators are

niB=bi†bin_i^{\mathrm B} = b_i^\dagger b_i

and

niF=ci†ci.n_i^{\mathrm F} = c_i^\dagger c_i.

The total number operator is

N^=∑ini.\hat N = \sum_i n_i.

Occupation numbers are eigenvalues of nin_i, not probabilities. Bosonic occupations are nonnegative integers; a single spinless fermionic mode has occupation zero or one.

The operator construction belongs to Creation, Annihilation, and Second Quantization.

For spin-1/21/2 degrees of freedom,

Siα=ℏ2σiα,α∈{x,y,z}.S_i^\alpha = \frac{\hbar}{2} \sigma_i^\alpha, \qquad \alpha\in\{x,y,z\}.

The Pauli matrices σα\sigma^\alpha are dimensionless. The spin operators SαS^\alpha carry units of angular momentum unless ℏ=1\hbar=1 has been declared.

The raising and lowering operators may be defined as

Si±=Six±iSiy,S_i^\pm = S_i^x \pm iS_i^y,

or in dimensionless form using σi±\sigma_i^\pm. A page must not switch between S±S^\pm and σ±\sigma^\pm without the corresponding factor of ℏ\hbar.

See Spin and Pauli-Matrix Conventions.

Continuum bosonic or fermionic fields are denoted

ψσ(r),ψσ†(r),\psi_\sigma(\mathbf r), \qquad \psi_\sigma^\dagger(\mathbf r),

where σ\sigma labels spin, species, flavor, or another internal component.

At equal time, bosonic fields satisfy

[ψσ(r),ψσ′†(r′)]=δσσ′δ(d)(r−r′),\left[ \psi_\sigma(\mathbf r), \psi_{\sigma'}^\dagger(\mathbf r') \right] = \delta_{\sigma\sigma'} \delta^{(d)} (\mathbf r-\mathbf r'),

while fermionic fields satisfy

{ψσ(r),ψσ′†(r′)}=δσσ′δ(d)(r−r′).\left\{ \psi_\sigma(\mathbf r), \psi_{\sigma'}^\dagger(\mathbf r') \right\} = \delta_{\sigma\sigma'} \delta^{(d)} (\mathbf r-\mathbf r').

The local number density is

n(r)=∑σψσ†(r)ψσ(r).n(\mathbf r) = \sum_\sigma \psi_\sigma^\dagger(\mathbf r) \psi_\sigma(\mathbf r).

The field operators are operator-valued distributions. Expressions at coincident points may require regularization or a lattice cutoff.

The canonical definitions and mode-expansion construction live in Field Operators. Their use in continuum and lattice many-body models is developed in Field Operators in Many-Body Models.

In a discrete one-particle basis, a number-conserving one-body operator is

H1=∑ijhijai†aj.H_1 = \sum_{ij} h_{ij} a_i^\dagger a_j.

A two-body interaction is written

H2=12∑ijklVijklai†aj†alak.H_2 = \frac{1}{2} \sum_{ijkl} V_{ijkl} a_i^\dagger a_j^\dagger a_l a_k.

The factor 1/21/2 and index ordering depend on how the matrix elements are defined. Antisymmetrized fermionic matrix elements can move numerical factors between VijklV_{ijkl} and the operator sum.

In the continuum, a common form is

H=∫ddr ψ†(r)h0ψ(r)+12∫ddr ddr′ ψ†(r)ψ†(r′)×V(r−r′)ψ(r′)ψ(r).\begin{aligned} H = {}& \int d^dr\, \psi^\dagger(\mathbf r) h_0 \psi(\mathbf r) \\ &+ \frac{1}{2} \int d^dr\,d^dr'\, \psi^\dagger(\mathbf r) \psi^\dagger(\mathbf r') \\ &\qquad\times V(\mathbf r-\mathbf r') \psi(\mathbf r') \psi(\mathbf r). \end{aligned}

Spin or species sums are implicit only when the page says so.

On a lattice, i,ji,j label sites and ⟨i,j⟩\langle i,j\rangle usually denotes a set of neighboring unordered pairs. A hopping Hamiltonian may be written

Ht=−t∑⟨i,j⟩,σ(ciσ†cjσ+cjσ†ciσ).H_t = -t \sum_{\langle i,j\rangle,\sigma} \left( c_{i\sigma}^\dagger c_{j\sigma} + c_{j\sigma}^\dagger c_{i\sigma} \right).

If the sum instead runs over oriented pairs, the Hermitian-conjugate term and counting convention change. The notation ⟨i,j⟩\langle i,j\rangle is not self-defining.

For the Hubbard interaction,

HU=U∑ini↑ni↓.H_U = U \sum_i n_{i\uparrow}n_{i\downarrow}.

The sign of tt can sometimes be changed by a gauge transformation on a bipartite lattice, but not in every geometry or flux sector.

For NsN_s sites at positions rj\mathbf r_j with periodic boundary conditions, this volume uses the symmetric discrete normalization

ckσ=1Ns∑je−ik⋅rjcjσ,c_{\mathbf k\sigma} = \frac{1}{\sqrt{N_s}} \sum_j e^{-i\mathbf k\cdot\mathbf r_j} c_{j\sigma},

with inverse

cjσ=1Ns∑keik⋅rjckσ.c_{j\sigma} = \frac{1}{\sqrt{N_s}} \sum_{\mathbf k} e^{i\mathbf k\cdot\mathbf r_j} c_{\mathbf k\sigma}.

The allowed k\mathbf k values depend on lattice vectors and boundary conditions. With these conventions,

{ckσ,ck′σ′†}=δkk′δσσ′.\{c_{\mathbf k\sigma}, c_{\mathbf k'\sigma'}^\dagger\} = \delta_{\mathbf k\mathbf k'} \delta_{\sigma\sigma'}.

Some references put the entire 1/Ns1/N_s factor in one direction. Formulas for momentum sums and Kronecker deltas must be translated consistently.

For a continuum field in dd dimensions, a symmetric wave-number convention is

ψσ(r)=∫ddk(2π)d/2eik⋅rakσ,\psi_\sigma(\mathbf r) = \int \frac{d^dk} {(2\pi)^{d/2}} e^{i\mathbf k\cdot\mathbf r} a_{\mathbf k\sigma},

with inverse

akσ=∫ddr(2π)d/2e−ik⋅rψσ(r).a_{\mathbf k\sigma} = \int \frac{d^dr} {(2\pi)^{d/2}} e^{-i\mathbf k\cdot\mathbf r} \psi_\sigma(\mathbf r).

This uses wave number k\mathbf k. Momentum is

p=ℏk.\mathbf p = \hbar\mathbf k.

The site-wide Fourier-Transform Conventions page uses momentum-normalized wavefunctions. Translate the measure and normalization when moving between p\mathbf p and k\mathbf k.

Finite-volume box normalization replaces integrals by sums. State the replacement, for example

1V∑k⟶∫ddk(2π)d\frac{1}{V} \sum_{\mathbf k} \longrightarrow \int \frac{d^dk}{(2\pi)^d}

in a periodic continuum limit.

For Heisenberg-picture operators,

A(t)=eiHt/ℏAe−iHt/ℏ.A(t) = e^{iHt/\hbar} A e^{-iHt/\hbar}.

An unordered two-point function is

CAB(t)=⟨A(t)B(0)⟩.C_{AB}(t) = \langle A(t)B(0)\rangle.

The connected function is

CABc(t)=⟨A(t)B(0)⟩−⟨A(t)⟩⟨B(0)⟩.C_{AB}^{\mathrm c}(t) = \langle A(t)B(0)\rangle - \langle A(t)\rangle \langle B(0)\rangle.

For a stationary state commuting with HH, the correlator depends only on the time difference.

The superscript c\mathrm c means connected here. It should not be confused with a creation operator or complex conjugation.

A retarded response function is commonly defined by

χABR(t)=−iℏθ(t)⟨[A(t),B(0)]⟩.\chi_{AB}^{\mathrm R}(t) = -\frac{i}{\hbar} \theta(t) \langle[A(t),B(0)]\rangle.

Some fields use the symbol GABRG^{\mathrm R}_{AB} for the same structure and reserve χ\chi for a specific physical response. The perturbation convention determines possible signs.

A time-ordered Green function for fermionic fields is often

G(1,2)=−i⟨Tψ(1)ψ†(2)⟩G(1,2) = -i \langle \mathcal T \psi(1) \psi^\dagger(2) \rangle

when ℏ=1\hbar=1. Other conventions include 1/ℏ1/\hbar or change the sign for bosons. Every Green-function page must state the definition before using a diagram or Fourier transform.

The volume will use:

  • GG for Green functions;
  • χ\chi for response or susceptibility;
  • CC for generic correlations;
  • S(q,ω)S(\mathbf q,\omega) for a dynamic structure factor.

For local density fluctuations

δnj=nj−⟨nj⟩,\delta n_j = n_j-\langle n_j\rangle,

a finite-lattice static structure factor may be defined as

S(q)=1Ns∑jℓe−iq⋅(rj−rℓ)⟨δnjδnℓ⟩.S(\mathbf q) = \frac{1}{N_s} \sum_{j\ell} e^{-i\mathbf q\cdot (\mathbf r_j-\mathbf r_\ell)} \langle \delta n_j \delta n_\ell \rangle.

The prefactor may instead involve particle number NN. State the normalization because it changes the scaling of Bragg peaks and sum rules.

A dynamic structure factor adds a time transform and an operator-ordering convention:

S(q,ω)∝∫dt eiωt⟨δnq(t)δn−q(0)⟩.S(\mathbf q,\omega) \propto \int dt\, e^{i\omega t} \langle \delta n_{\mathbf q}(t) \delta n_{-\mathbf q}(0) \rangle.

The proportionality sign is intentional here: the exact 2π2\pi, volume, and normalization factors belong on the dedicated definition page.

Imaginary time is denoted τ\tau. A thermal Heisenberg operator is

A(τ)=eτHAe−τHA(\tau) = e^{\tau H} A e^{-\tau H}

when ℏ=1\hbar=1 and τ\tau has inverse-energy units.

If τ\tau instead has units of time, the exponent is

eτH/ℏ.e^{\tau H/\hbar}.

Do not mix these conventions.

Imaginary Time gives the operator-semigroup, heat-kernel, projection, and thermal-trace dictionary with ℏ\hbar explicit.

For inverse-temperature units with ℏ=1\hbar=1, bosonic and fermionic Matsubara frequencies are

ωnB=2πnβ,ωnF=(2n+1)πβ.\begin{aligned} \omega_n^{\mathrm B} &= \frac{2\pi n}{\beta}, \\ \omega_n^{\mathrm F} &= \frac{(2n+1)\pi}{\beta}. \end{aligned}

With explicit ℏ\hbar and time-valued τ\tau, the corresponding angular frequencies carry 1/ℏ1/\hbar relative to the energy convention.

Thermodynamic entropy is

Sth=−kBTr⁡(ρln⁡ρ).S_{\mathrm{th}} = -k_{\mathrm B} \operatorname{Tr} \left( \rho\ln\rho \right).

The dimensionless von Neumann entropy is

SvN(ρ)=−Tr⁡(ρln⁡ρ).S_{\mathrm{vN}}(\rho) = -\operatorname{Tr} \left( \rho\ln\rho \right).

For a reduced state ρA\rho_A, entanglement entropy of a global pure state is

SA=−Tr⁡(ρAln⁡ρA).S_A = -\operatorname{Tr} \left( \rho_A\ln\rho_A \right).

Thus

Sth=kBSvNS_{\mathrm{th}} = k_{\mathrm B} S_{\mathrm{vN}}

for the same density operator, but thermodynamic and entanglement interpretations depend on which state and partition are used.

The logarithm base is natural unless stated. Base-two entropy should be labeled explicitly and is measured in bits.

The thermodynamic limit is written

N,V→∞,NV=nfixed.N,V\to\infty, \qquad \frac{N}{V}=n \quad\text{fixed}.

For a lattice,

Ns→∞N_s\to\infty

with filling, geometry, and coupling scaling specified.

The symbols have distinct meanings:

SymbolMeaning
f∼gf\sim gratio approaches a stated constant, often one
f∝gf\propto gproportional with an unspecified constant
f≃gf\simeq gapproximate equality in a stated regime
f=O(g)f=O(g)bounded in magnitude by a constant times gg asymptotically
f=o(g)f=o(g)f/g→0f/g\to0
Δ\Deltagap, difference, or fluctuation, as locally defined

Do not use ∼\sim merely to mean “roughly similar.”

An extensive quantity XX scales with system size:

X(λN,λV)≃λX(N,V)X(\lambda N,\lambda V) \simeq \lambda X(N,V)

in an appropriate thermodynamic regime.

The associated density

x=XVx = \frac{X}{V}

is intensive. Boundary contributions can scale as Ld−1L^{d-1} and finite-size corrections can violate exact proportionality.

For long-range interactions, the energy may be superextensive unless the coupling is normalized with system size. Extensive and Intensive Quantities is the canonical scaling audit for limiting densities, additivity, boundary terms, fluctuations, and Kac normalization.

This volume keeps ℏ\hbar and kBk_{\mathrm B} explicit on orientation and derivation pages unless setting them to one materially clarifies advanced notation.

If natural many-body units are used, state them near the first formula:

ℏ=1,kB=1,a=1\hbar=1, \qquad k_{\mathrm B}=1, \qquad a=1

are three independent choices.

Setting a=1a=1 makes momentum dimensionless in lattice units but does not make physical momentum dimensionless after units are restored. Setting kB=1k_{\mathrm B}=1 makes temperature carry energy units.

The combinations in an equilibrium exponent must be dimensionless:

βH,βμN^.\beta H, \qquad \beta\mu\hat N.

Use:

  • OBC for open boundary conditions;
  • PBC for periodic boundary conditions;
  • APBC for antiperiodic boundary conditions;
  • twisted boundary conditions with an explicit phase θ\theta.

For a one-dimensional periodic lattice,

cj+Ns=cj.c_{j+N_s} = c_j.

For a twist,

cj+Ns=eiθcj.c_{j+N_s} = e^{i\theta}c_j.

Boundary conditions determine the momentum grid and can change finite-size degeneracies, topology probes, and fermionic sign conventions.

SymbolPossible meaningsPreferred repair
NNparticles, sites, flavors, large-NN rankuse NsN_s, NfN_f, or a sentence
LLlength, number of sites, Liouvillianuse NsN_s for sites and L\mathcal L for Liouvillian
SSentropy, spin, action, structure factoruse SthS_{\mathrm{th}}, S\mathbf S, S\mathcal S, S(q)S(\mathbf q)
Ω\Omegagrand potential, angular frequency, volumeuse ω\omega for frequency and VV for volume
Δ\Deltagap, detuning, order parameter, finite differenceadd a descriptive subscript
Γ\Gammarate, linewidth, vertexdefine units and role
χ\chisusceptibility, response, spinoruse arguments and indices
GGGreen function, conductance, Gibbs free energyuse GRG^{\mathrm R}, G\mathcal G, or descriptive text
TTtemperature, time ordering, transposeuse T\mathcal T for time ordering
σ\sigmaPauli matrix, spin index, conductivitydistinguish σ\boldsymbol\sigma, index position, or σ(ω)\sigma(\omega)

Overloading is acceptable across different pages when local definitions are unambiguous. It is dangerous within one derivation.

Consider

H=−t∑⟨i,j⟩,σ(ciσ†cjσ+h.c.)+U∑ini↑ni↓.\begin{aligned} H = {}& -t \sum_{\langle i,j\rangle,\sigma} \left( c_{i\sigma}^\dagger c_{j\sigma} + \mathrm{h.c.} \right) \\ &+ U \sum_i n_{i\uparrow}n_{i\downarrow}. \end{aligned}

The notation contract is:

  • i,ji,j label lattice sites;
  • σ∈{↑,↓}\sigma\in\{\uparrow,\downarrow\} labels spin;
  • ⟨i,j⟩\langle i,j\rangle is an unordered nearest-neighbor bond set;
  • tt has energy units and is taken real here;
  • UU is the onsite interaction energy;
  • ciσc_{i\sigma} obeys fermionic anticommutation relations;
  • niσ=ciσ†ciσn_{i\sigma}=c_{i\sigma}^\dagger c_{i\sigma};
  • boundary conditions and lattice geometry remain to be stated;
  • particle number
N^=∑i,σniσ\hat N = \sum_{i,\sigma} n_{i\sigma}

is conserved.

For a translationally invariant one-dimensional chain with lattice spacing aa and PBC,

ckσ=1Ns∑je−ikajcjσ.c_{k\sigma} = \frac{1}{\sqrt{N_s}} \sum_j e^{-ikaj} c_{j\sigma}.

The kinetic term becomes diagonal:

Ht=∑k,σϵkckσ†ckσ,H_t = \sum_{k,\sigma} \epsilon_k c_{k\sigma}^\dagger c_{k\sigma},

with

ϵk=−2tcos⁡(ka).\epsilon_k = -2t\cos(ka).

The interaction is local in real space but couples several momenta in momentum space. Representation choice exposes one structure while obscuring another.

For the Hilbert space, parameter conventions, exact limits, symmetries, and physical diagnostics behind this notation, see Hubbard Model.

Before using a many-body formula, ask:

  1. Are indices particles, sites, modes, bands, spin, or species?
  2. Is the state pure, canonical, grand canonical, or nonequilibrium?
  3. Is particle number fixed or fluctuating?
  4. Are operators in Schrödinger, Heisenberg, interaction, or imaginary-time form?
  5. Are Fourier variables momentum or wave number?
  6. Are sums normalized by NN, NsN_s, or VV?
  7. Is a correlation connected, ordered, symmetrized, retarded, or time ordered?
  8. Are ℏ\hbar, kBk_{\mathrm B}, and aa explicit?
  9. Which boundary conditions determine the mode grid?
  10. Which limit is taken first?

This page owns the default notation contract for the Many-Body and Quantum Statistical Mechanics volume.

  • Using NN for both sites and particles in the same equation.
  • Leaving spin, band, or species sums implicit without warning.
  • Confusing an occupation eigenvalue with a probability.
  • Omitting the factor relating SαS^\alpha and σα\sigma^\alpha.
  • Moving between wave number and momentum without ℏ\hbar.
  • Combining Fourier formulas with incompatible normalizations.
  • Forgetting whether a bond sum is oriented or unordered.
  • Double counting a two-body interaction.
  • Using a Green function before defining its sign and ordering convention.
  • Comparing structure factors normalized by NN and by NsN_s as if identical.
  • Mixing energy-valued and time-valued imaginary time.
  • Taking a thermodynamic derivative while an implicit parameter also varies.
  • Setting ℏ\hbar, kBk_{\mathrm B}, or aa to one silently.
  • Writing a thermodynamic limit without density, geometry, or coupling scaling.

Show that β(H−μN^)\beta(H-\mu\hat N) is dimensionless and state the units of β\beta and μ\mu when kBk_{\mathrm B} is explicit.

Solution

The operator HH has energy units. The number operator N^\hat N is dimensionless, so μ\mu has energy units and μN^\mu\hat N also has energy units.

Since

β=1kBT,\beta = \frac{1}{k_{\mathrm B}T},

β\beta has inverse-energy units. Therefore

β(H−μN^)\beta(H-\mu\hat N)

is dimensionless, as required inside an exponential.

Using

ckσ=1Ns∑je−ikrjcjσ,c_{k\sigma} = \frac{1}{\sqrt{N_s}} \sum_j e^{-ik r_j} c_{j\sigma},

show that the momentum-space operators obey the canonical anticommutator for the allowed discrete momenta.

Solution

Compute

{ckσ,ck′σ′†}=1Ns∑jℓe−ikrjeik′rℓ×{cjσ,cℓσ′†}=δσσ′Ns∑jei(k′−k)rj.\begin{aligned} \{c_{k\sigma},c_{k'\sigma'}^\dagger\} &= \frac{1}{N_s} \sum_{j\ell} e^{-ikr_j} e^{ik'r_\ell} \\ &\qquad{}\times \{c_{j\sigma},c_{\ell\sigma'}^\dagger\} \\ &= \frac{\delta_{\sigma\sigma'}}{N_s} \sum_j e^{i(k'-k)r_j}. \end{aligned}

Discrete Fourier orthogonality gives

1Ns∑jei(k′−k)rj=δkk′.\frac{1}{N_s} \sum_j e^{i(k'-k)r_j} = \delta_{kk'}.

Hence

{ckσ,ck′σ′†}=δkk′δσσ′.\{c_{k\sigma},c_{k'\sigma'}^\dagger\} = \delta_{kk'} \delta_{\sigma\sigma'}.

For fermions, show that

[N^,cj]=−cj.[\hat N,c_j] = -c_j.
Solution

With

N^=∑ici†ci,\hat N = \sum_i c_i^\dagger c_i,

only the i=ji=j term contributes. Using the fermionic algebra,

[cj†cj,cj]=cj†cjcj−cjcj†cj=0−(1−cj†cj)cj=−cj.\begin{aligned} [c_j^\dagger c_j,c_j] &= c_j^\dagger c_jc_j - c_jc_j^\dagger c_j \\ &= 0 - (1-c_j^\dagger c_j)c_j \\ &= -c_j. \end{aligned}

Thus annihilation lowers total particle number by one. Similarly,

[N^,cj†]=cj†.[\hat N,c_j^\dagger] = c_j^\dagger.

Given

Ω=−1βln⁡Ξ,\Omega = -\frac{1}{\beta} \ln\Xi,

show that

⟨N^⟩=−(∂Ω∂μ)β,V.\langle\hat N\rangle = -\left( \frac{\partial\Omega}{\partial\mu} \right)_{\beta,V}.
Solution

At fixed β\beta,

∂ln⁡Ξ∂μ=β⟨N^⟩.\frac{\partial\ln\Xi}{\partial\mu} = \beta \langle\hat N\rangle.

Therefore

∂Ω∂μ=−1β∂ln⁡Ξ∂μ=−⟨N^⟩.\frac{\partial\Omega}{\partial\mu} = -\frac{1}{\beta} \frac{\partial\ln\Xi}{\partial\mu} = -\langle\hat N\rangle.

Rearranging gives the result.

Suppose ⟨A⟩=2\langle A\rangle=2, ⟨B⟩=3\langle B\rangle=3, and ⟨AB⟩=7\langle AB\rangle=7. Find the connected correlation and explain what the disconnected contribution is.

Solution

The connected correlation is

⟨AB⟩c=7−(2)(3)=1.\langle AB\rangle_{\mathrm c} = 7-(2)(3) =1.

The disconnected contribution is

⟨A⟩⟨B⟩=6.\langle A\rangle \langle B\rangle =6.

It is the part already explained by the separate one-point expectations. The remaining value measures failure of the two-point expectation to factorize.

  • A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, Dover (2003).
  • G. D. Mahan, Many-Particle Physics, 3rd ed., Springer (2000).
  • P. Coleman, Introduction to Many-Body Physics, Cambridge University Press (2015).
  • A. Altland and B. Simons, Condensed Matter Field Theory, 2nd ed., Cambridge University Press (2010).
  • G. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid, Cambridge University Press (2005).
  • R. K. Pathria and P. D. Beale, Statistical Mechanics, 3rd ed., Elsevier (2011).
  • M. Le Bellac, F. Mortessagne, and G. G. Batrouni, Equilibrium and Non-Equilibrium Statistical Thermodynamics, Cambridge University Press (2004).
  • J. W. Negele and H. Orland, Quantum Many-Particle Systems, Westview Press (1998).