Skip to content

Conventions for Quantum Matter

This page defines the default notation contract for lattices, Bloch states, band Hamiltonians, Berry geometry, electromagnetic coupling, response functions, and material observables. It supplements the sitewide Conventions Overview, Fourier Transform Conventions, and Many-Body Symbols and Conventions.

For magnetic claims, Magnetism and Spin Systems supplies the moment, normalization, exchange, ordering-wavevector, field-history, and probe ledger; the convention definitions themselves remain canonical here and in their specialist owners.

Choosing a Model for Quantum Matter owns the preceding choice of retained description. This page begins after that choice and fixes the signs, bases, normalizations, units, and limits needed to compare calculations. Superfluidity and Superconductivity owns the chapter-level choice among pairing, stiffness, electromagnetic response, defects, weak links, and proximity; this page supplies the charge, gauge, Fourier, and response conventions used along those routes. It does not choose the local crystalline learning branch; use the Lattices, Reciprocal Space, and Bloch Electrons gateway for that dependency audit, then return here before comparing formulas.

When the calculation produces band data, pair this notation contract with the Band Theory and Electronic Structure gateway. Its band ledger records the description and parameter provenance, basis and symmetry content, energy reference, filling and state, direct or indirect gap claim, temperature and broadening, target observable, and validity window. This page owns the conventions used to state those fields; the gateway owns the downstream interpretation route.

When the calculation instead produces a dynamical matrix, mode spectrum, dielectric pole, or dressed material excitation, pair this notation contract with the Lattice Vibrations and Collective Modes gateway. This page fixes cell phases, angular-frequency and energy variables, normalizations, response signs, and orders of limits; the gateway owns the status-aware material route.

A specialized article may use another standard convention when that choice clarifies the subject. It must declare the override near the first affected formula and translate back when comparing results. A silent change of sign, normalization, basis embedding, or tensor order is an error.

A reusable quantum-matter calculation must make the following information recoverable:

ItemRequired data
latticeprimitive vectors, basis positions, dimension, orientation
finite geometrynumber of cells, boundary conditions, sample shape
reciprocal spacereciprocal vectors, Brillouin-zone representative, momentum grid
basisorbital, sublattice, layer, spin, band, and gauge conventions
Fourier transformphases, normalization, and whether basis positions enter
Hamiltonianhopping orientation, interaction factors, charge signs, units
statefilling, chemical potential, temperature, preparation
responseperturbation sign, current or density convention, Fourier variable
geometryfield, current, polarization, surface, and tensor-axis directions
limitsthermodynamic, clean, dc, uniform, zero-temperature, and their order

The compact rule is:

A formula is defined only together with its basis, normalization, orientation, charge, units, and limits.

SymbolDefault meaning
r\mathbf rcontinuum position
R\mathbf RBravais-lattice vector
τα\boldsymbol\tau_\alphaposition of basis object α\alpha within the chosen cell
rRα\mathbf r_{\mathbf R\alpha}embedded position R+τα\mathbf R+\boldsymbol\tau_\alpha
ai\mathbf a_iprimitive real-space vector
bi\mathbf b_iprimitive reciprocal-lattice vector
k\mathbf kcrystal wavevector or Bloch label
q\mathbf qtransferred wavevector
G\mathbf Greciprocal-lattice vector
i,ji,jsites or Cartesian components, declared locally
α,β\alpha,\betaorbital, sublattice, layer, or basis indices
n,mn,mband, eigenstate, integer, or Matsubara indices
σ,σ′\sigma,\sigma'physical spin unless explicitly called pseudospin
NcN_cnumber of primitive cells
NorbN_{\mathrm{orb}}number of retained orbitals per primitive cell
Ωc\Omega_cprimitive-cell volume or area
VVphysical sample volume or area in the modeled dimension

Repeated Cartesian indices are not summed unless an article declares an Einstein convention. Lattice, orbital, band, and spin sums are normally explicit.

Boldface denotes geometric vectors. It does not denote a generic column vector in orbital space; write unα(k)u_{n\alpha}(\mathbf k) or declare a bold orbital-space vector locally.

In dd dimensions, a Bravais vector is

R=∑i=1dniai,ni∈Z.\mathbf R = \sum_{i=1}^{d} n_i\mathbf a_i, \qquad n_i\in\mathbb Z.

The primitive vectors are ordered. In three dimensions the default orientation is right handed:

Ωc=a1⋅(a2×a3)>0.\Omega_c = \mathbf a_1\cdot \left( \mathbf a_2\times\mathbf a_3 \right) > 0.

In two dimensions embedded in the xyxy plane,

Ωc=(a1×a2)⋅z^>0.\Omega_c = \left( \mathbf a_1\times\mathbf a_2 \right)\cdot\hat{\mathbf z} > 0.

If a source uses a left-handed ordering, reverse the orientation before comparing Berry-curvature or Hall signs.

A crystal site or localized orbital is embedded at

rRα=R+τα.\mathbf r_{\mathbf R\alpha} = \mathbf R + \boldsymbol\tau_\alpha.

τα\boldsymbol\tau_\alpha is defined modulo a Bravais vector. Moving a basis representative into a neighboring cell changes Fourier phases and matrix representations but not physical positions.

An orbital is not always centered exactly at an atomic coordinate. Wannier centers, bond orbitals, and effective degrees of freedom may have different embeddings. Record the position used by the position operator and by electromagnetic phases.

For a periodic finite lattice with NiN_i cells along ai\mathbf a_i,

Nc=∏i=1dNi.N_c = \prod_{i=1}^{d}N_i.

When the modeled sample is exactly tiled by primitive cells,

V=NcΩc.V = N_c\Omega_c.

Use NcN_c for the number of primitive cells rather than the ambiguous NsN_s. The number of localized modes is NcNorbN_cN_{\mathrm{orb}} before spin or other internal multiplicities.

The reciprocal basis is defined by

ai⋅bj=2πδij.\mathbf a_i\cdot\mathbf b_j = 2\pi\delta_{ij}.

A reciprocal vector is

G=∑i=1dmibi,mi∈Z.\mathbf G = \sum_{i=1}^{d} m_i\mathbf b_i, \qquad m_i\in\mathbb Z.

With the right-handed three-dimensional convention,

b1=2πa2×a3Ωc,b2=2πa3×a1Ωc,b3=2πa1×a2Ωc.\begin{aligned} \mathbf b_1 &= 2\pi \frac{\mathbf a_2\times\mathbf a_3}{\Omega_c}, \\ \mathbf b_2 &= 2\pi \frac{\mathbf a_3\times\mathbf a_1}{\Omega_c}, \\ \mathbf b_3 &= 2\pi \frac{\mathbf a_1\times\mathbf a_2}{\Omega_c}. \end{aligned}

The reciprocal-cell volume is

ΩBZ=(2π)dΩc.\Omega_{\mathrm{BZ}} = \frac{(2\pi)^d}{\Omega_c}.

Here ΩBZ\Omega_{\mathrm{BZ}} means the measure of any primitive reciprocal cell, including the first Brillouin zone.

The default first Brillouin zone is the Wigner–Seitz cell of the reciprocal lattice around k=0\mathbf k=\mathbf0. Any primitive reciprocal cell contains one representative of each equivalence class

k∼k+G.\mathbf k \sim \mathbf k+\mathbf G.

Boundary points have multiple equivalent representatives. A numerical mesh must assign each boundary state once. Integrals are insensitive to that measure-zero choice, while discrete sums are not.

For periodic boundary conditions over NiaiN_i\mathbf a_i, use

km=∑i=1dmiNibi,mi=0,…,Ni−1.\mathbf k_{\mathbf m} = \sum_{i=1}^{d} \frac{m_i}{N_i} \mathbf b_i, \qquad m_i=0,\ldots,N_i-1.

The listed grid is a parallelepiped representative. Points may be folded into the first Brillouin zone for plotting.

With a boundary twist θi\theta_i,

km(θ)=∑i=1dmi+θi/(2π)Nibi.\mathbf k_{\mathbf m}(\boldsymbol\theta) = \sum_{i=1}^{d} \frac{m_i+\theta_i/(2\pi)}{N_i} \mathbf b_i.

State whether θ\boldsymbol\theta is a boundary phase, a physical flux, or an auxiliary parameter used for many-body geometry.

k\mathbf k has units of inverse length. In a continuum free-particle problem,

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

is physical momentum. In a crystal, ℏk\hbar\mathbf k is commonly called crystal momentum, but it is defined modulo ℏG\hbar\mathbf G and need not equal mechanical momentum or mvm\mathbf v.

The band velocity of an isolated differentiable band is

vn(k)=1ℏ∇kεn(k).\mathbf v_n(\mathbf k) = \frac{1}{\hbar} \boldsymbol\nabla_{\mathbf k} \varepsilon_n(\mathbf k).

The effective inverse-mass tensor is

[(mn∗)−1]ij=1ℏ2∂2εn∂ki∂kj.\left[ \left( m_n^\ast \right)^{-1} \right]_{ij} = \frac{1}{\hbar^2} \frac{\partial^2\varepsilon_n} {\partial k_i\partial k_j}.

Neither equation defines a global particle mass. Both describe local band response under the assumptions of semiclassical wave-packet dynamics.

For an annihilation operator cRασc_{\mathbf R\alpha\sigma}, the default lattice transform is

ckασ=1Nc∑Re−ik⋅RcRασ,cRασ=1Nc∑keik⋅Rckασ.\begin{aligned} c_{\mathbf k\alpha\sigma} &= \frac{1}{\sqrt{N_c}} \sum_{\mathbf R} e^{-i\mathbf k\cdot\mathbf R} c_{\mathbf R\alpha\sigma}, \\ c_{\mathbf R\alpha\sigma} &= \frac{1}{\sqrt{N_c}} \sum_{\mathbf k} e^{i\mathbf k\cdot\mathbf R} c_{\mathbf k\alpha\sigma}. \end{aligned}

The basis position τα\boldsymbol\tau_\alpha is not included in the phase. This is called the cell convention.

Discrete orthogonality is

1Nc∑Rei(k−k′)⋅R=δk,k′,1Nc∑keik⋅(R−R′)=δR,R′.\begin{aligned} \frac{1}{N_c} \sum_{\mathbf R} e^{i(\mathbf k-\mathbf k')\cdot\mathbf R} &= \delta_{\mathbf k,\mathbf k'}, \\ \frac{1}{N_c} \sum_{\mathbf k} e^{i\mathbf k\cdot(\mathbf R-\mathbf R')} &= \delta_{\mathbf R,\mathbf R'}. \end{aligned}

The symmetric normalization preserves canonical anticommutators or commutators without extra cell factors.

Some sources define

cˉkασ=1Nc∑Re−ik⋅(R+τα)cRασ.\bar c_{\mathbf k\alpha\sigma} = \frac{1}{\sqrt{N_c}} \sum_{\mathbf R} e^{-i\mathbf k\cdot (\mathbf R+\boldsymbol\tau_\alpha)} c_{\mathbf R\alpha\sigma}.

The two conventions are related by

cˉkασ=e−ik⋅ταckασ.\bar c_{\mathbf k\alpha\sigma} = e^{-i\mathbf k\cdot\boldsymbol\tau_\alpha} c_{\mathbf k\alpha\sigma}.

Define the diagonal embedding matrix

Dαβ(k)=δαβe−ik⋅τα.D_{\alpha\beta}(\mathbf k) = \delta_{\alpha\beta} e^{-i\mathbf k\cdot\boldsymbol\tau_\alpha}.

Then

cˉk=D(k)ck.\bar{\mathbf c}_{\mathbf k} = D(\mathbf k)\mathbf c_{\mathbf k}.

The two bases give identical observables when Hamiltonians, eigenvectors, position operators, and derivatives are transformed together.

In the thermodynamic limit,

1Nc∑k⟶Ωc(2π)d∫BZddk.\frac{1}{N_c} \sum_{\mathbf k} \longrightarrow \frac{\Omega_c}{(2\pi)^d} \int_{\mathrm{BZ}} d^dk.

Equivalently,

∑k⟶V(2π)d∫BZddk.\sum_{\mathbf k} \longrightarrow \frac{V}{(2\pi)^d} \int_{\mathrm{BZ}} d^dk.

The second form assumes V=NcΩcV=N_c\Omega_c. A formula per cell uses the first normalization; a formula per physical volume uses the second.

Define

H0=∑R,Δ∑α,βtαβ(Δ)cRα†cR+Δ,β.H_0 = \sum_{\mathbf R,\boldsymbol\Delta} \sum_{\alpha,\beta} t_{\alpha\beta}(\boldsymbol\Delta) c_{\mathbf R\alpha}^\dagger c_{\mathbf R+\boldsymbol\Delta,\beta}.

tαβ(Δ)t_{\alpha\beta}(\boldsymbol\Delta) is the amplitude multiplying a hop from (R+Δ,β)(\mathbf R+\boldsymbol\Delta,\beta) to (R,α)(\mathbf R,\alpha) in this operator order. Hermiticity requires

tαβ(Δ)=tβα∗(−Δ).t_{\alpha\beta}(\boldsymbol\Delta) = t_{\beta\alpha}^\ast(-\boldsymbol\Delta).

With the default cell transform,

H0=∑kck†h(k)ck,H_0 = \sum_{\mathbf k} \mathbf c_{\mathbf k}^\dagger h(\mathbf k) \mathbf c_{\mathbf k},

where

hαβ(k)=∑Δtαβ(Δ)eik⋅Δ.h_{\alpha\beta}(\mathbf k) = \sum_{\boldsymbol\Delta} t_{\alpha\beta}(\boldsymbol\Delta) e^{i\mathbf k\cdot\boldsymbol\Delta}.

For a one-dimensional chain written as

H=−t∑j(cj†cj+1+cj+1†cj),H = -t \sum_j \left( c_j^\dagger c_{j+1} +c_{j+1}^\dagger c_j \right),

the dispersion is

ε(k)=−2tcos⁡(ka).\varepsilon(k) = -2t\cos(ka).

The minus sign belongs to the Hamiltonian definition. Some sources absorb it into a signed hopping parameter.

In the cell convention,

h(k+G)=h(k).h(\mathbf k+\mathbf G) = h(\mathbf k).

In the orbital-position convention,

hˉ(k)=D(k)h(k)D†(k).\bar h(\mathbf k) = D(\mathbf k) h(\mathbf k) D^\dagger(\mathbf k).

Since

D(k+G)=DGD(k),(DG)αβ=δαβe−iG⋅τα,D(\mathbf k+\mathbf G) = D_{\mathbf G}D(\mathbf k), \qquad \left(D_{\mathbf G}\right)_{\alpha\beta} = \delta_{\alpha\beta} e^{-i\mathbf G\cdot\boldsymbol\tau_\alpha},

the embedded Hamiltonian is generally quasi-periodic:

hˉ(k+G)=DGhˉ(k)DG†.\bar h(\mathbf k+\mathbf G) = D_{\mathbf G} \bar h(\mathbf k) D_{\mathbf G}^\dagger.

Its eigenvalues remain periodic. Confusing matrix periodicity with spectral periodicity is a common source of apparent disagreement.

The active translation operator acts as

⟨r∣TR∣ψ⟩=ψ(r−R).\langle\mathbf r| T_{\mathbf R} |\psi\rangle = \psi(\mathbf r-\mathbf R).

A Bloch eigenstate satisfies

TR∣ψnk⟩=e−ik⋅R∣ψnk⟩.T_{\mathbf R} |\psi_{n\mathbf k}\rangle = e^{-i\mathbf k\cdot\mathbf R} |\psi_{n\mathbf k}\rangle.

Equivalently,

ψnk(r+R)=eik⋅Rψnk(r).\psi_{n\mathbf k}(\mathbf r+\mathbf R) = e^{i\mathbf k\cdot\mathbf R} \psi_{n\mathbf k}(\mathbf r).

Write

ψnk(r)=eik⋅runk(r),\psi_{n\mathbf k}(\mathbf r) = e^{i\mathbf k\cdot\mathbf r} u_{n\mathbf k}(\mathbf r),

where

unk(r+R)=unk(r).u_{n\mathbf k}(\mathbf r+\mathbf R) = u_{n\mathbf k}(\mathbf r).

This structural form suppresses any overall finite-volume normalization. Bloch’s Theorem gives the proof and the distinction between finite-volume eigenvectors and infinite-crystal generalized eigenfunctions.

For a finite periodic crystal, use

⟨ψnk∣ψmk′⟩=δnmδk,k′.\langle \psi_{n\mathbf k} | \psi_{m\mathbf k'} \rangle = \delta_{nm} \delta_{\mathbf k,\mathbf k'}.

For cell-periodic functions, define

⟨u∣v⟩c=1Ωc∫cellddr u∗(r)v(r),\langle u|v\rangle_c = \frac{1}{\Omega_c} \int_{\mathrm{cell}} d^dr\, u^\ast(\mathbf r)v(\mathbf r),

and normalize

⟨unk∣umk⟩c=δnm.\langle u_{n\mathbf k} | u_{m\mathbf k} \rangle_c = \delta_{nm}.

With this cell-average convention and V=NcΩcV=N_c\Omega_c, the corresponding globally normalized wavefunction is

ψnknorm(r)=1NcΩceik⋅runk(r).\psi_{n\mathbf k}^{\mathrm{norm}}(\mathbf r) = \frac{1}{\sqrt{N_c\Omega_c}} e^{i\mathbf k\cdot\mathbf r} u_{n\mathbf k}(\mathbf r).

The factor 1/NcΩc1/\sqrt{N_c\Omega_c} is omitted only when the structural Bloch form or a generalized infinite-volume state is intended.

An article using ∫cell∣u∣2=1\int_{\mathrm{cell}}|u|^2=1 has moved a factor of Ωc\Omega_c into the cell inner product. Translate all matrix elements consistently.

n,mn,m label energy bands after diagonalization. α,β\alpha,\beta label orbitals or sublattices before diagonalization:

∑βhαβ(k)unβ(k)=εn(k)unα(k).\sum_\beta h_{\alpha\beta}(\mathbf k) u_{n\beta}(\mathbf k) = \varepsilon_n(\mathbf k) u_{n\alpha}(\mathbf k).

At degeneracies, individual band eigenvectors are not unique. Use a projector onto the isolated band subspace when a gauge-invariant statement does not require choosing a basis inside it.

For an isolated band,

∣unk⟩⟶eiχn(k)∣unk⟩|u_{n\mathbf k}\rangle \longrightarrow e^{i\chi_n(\mathbf k)} |u_{n\mathbf k}\rangle

does not change the state ray or energy. This is the Bloch eigenvector gauge.

For a degenerate isolated subspace,

∣uak⟩⟶∑b∣ubk⟩Uba(k)|u_{a\mathbf k}\rangle \longrightarrow \sum_b |u_{b\mathbf k}\rangle U_{ba}(\mathbf k)

with unitary U(k)U(\mathbf k) is the corresponding nonabelian gauge freedom.

A k\mathbf k-dependent unitary change of orbital basis also transforms h(k)h(\mathbf k). This is distinct from multiplying a solved eigenvector by a phase. Derivatives with respect to k\mathbf k must include the basis connection when the orbital basis itself depends on k\mathbf k.

In a compatible k\mathbf k-independent cell basis, the velocity matrix is

vi(k)=1ℏ∂h(k)∂ki.v_i(\mathbf k) = \frac{1}{\hbar} \frac{\partial h(\mathbf k)}{\partial k_i}.

After a k\mathbf k-dependent basis transformation, differentiating only the transformed matrix can omit connection terms. State the basis before using derivative formulas for velocity, Berry curvature, or orbital response.

Berry Connection, Curvature, and Orientation

Section titled “Berry Connection, Curvature, and Orientation”

For a normalized cell-periodic eigenvector, use

An,i(k)=i⟨unk∣∂kiunk⟩c.\mathcal A_{n,i}(\mathbf k) = i \langle u_{n\mathbf k} | \partial_{k_i} u_{n\mathbf k} \rangle_c.

Under

∣unk⟩⟶eiχn(k)∣unk⟩,|u_{n\mathbf k}\rangle \longrightarrow e^{i\chi_n(\mathbf k)} |u_{n\mathbf k}\rangle,

the connection transforms as

An,i⟶An,i−∂kiχn.\mathcal A_{n,i} \longrightarrow \mathcal A_{n,i} - \partial_{k_i}\chi_n.

The curvature is

Ωn,ij=∂kiAn,j−∂kjAn,i.\Omega_{n,ij} = \partial_{k_i}\mathcal A_{n,j} - \partial_{k_j}\mathcal A_{n,i}.

In three-vector notation,

Ωn=∇k×An.\boldsymbol\Omega_n = \boldsymbol\nabla_{\mathbf k} \times \boldsymbol{\mathcal A}_n.

This sign agrees with Berry Connection and the Berry Curvature formula card.

For an oriented two-dimensional Brillouin zone with coordinates (kx,ky)(k_x,k_y),

Cn=12π∫BZdkx dky Ωn,xy(k).C_n = \frac{1}{2\pi} \int_{\mathrm{BZ}} dk_x\,dk_y\, \Omega_{n,xy}(\mathbf k).

Reversing the coordinate orientation changes the sign of both the area form and reported Chern number. Always state the real-space and reciprocal-space axis orientation before comparing signs.

A periodic gauge would satisfy

∣un,k+G⟩=∣unk⟩|u_{n,\mathbf k+\mathbf G}\rangle = |u_{n\mathbf k}\rangle

in a periodic orbital basis. Such a gauge may not be globally smooth for a topologically nontrivial isolated band. Numerical calculations should use gauge-covariant overlaps, projectors, or patching rather than enforce an impossible global phase choice.

Charge, Electromagnetic Coupling, and Peierls Phase

Section titled “Charge, Electromagnetic Coupling, and Peierls Phase”

e>0e>0 denotes the elementary charge magnitude. The electron charge is

qe=−e.q_{\mathrm e} = -e.

For a particle of charge qq, minimal coupling is

p⟶p−qA.\mathbf p \longrightarrow \mathbf p-q\mathbf A.

For an electron,

p−qeA=p+eA.\mathbf p-q_{\mathrm e}\mathbf A = \mathbf p+e\mathbf A.

The electromagnetic potentials obey

B=∇×A,E=−∇ϕ−∂A∂t.\begin{aligned} \mathbf B &= \boldsymbol\nabla\times\mathbf A, \\ \mathbf E &= - \boldsymbol\nabla\phi - \frac{\partial\mathbf A}{\partial t}. \end{aligned}

Use

A′=A+∇Λ,ϕ′=ϕ−∂Λ∂t,ψ′=eiqΛ/ℏψ.\begin{aligned} \mathbf A' &= \mathbf A + \boldsymbol\nabla\Lambda, \\ \phi' &= \phi - \frac{\partial\Lambda}{\partial t}, \\ \psi' &= e^{iq\Lambda/\hbar} \psi. \end{aligned}

For an electron, the wavefunction phase is e−ieΛ/ℏe^{-ie\Lambda/\hbar}.

For a hopping term tijci†cjt_{ij}c_i^\dagger c_j, which annihilates at jj and creates at ii, use

tij⟶tijexp⁡[iqℏ∫rjriA⋅dℓ].t_{ij} \longrightarrow t_{ij} \exp \left[ \frac{iq}{\hbar} \int_{\mathbf r_j}^{\mathbf r_i} \mathbf A\cdot d\boldsymbol\ell \right].

For an electron,

tij⟶tijexp⁡[−ieℏ∫rjriA⋅dℓ].t_{ij} \longrightarrow t_{ij} \exp \left[ - \frac{ie}{\hbar} \int_{\mathbf r_j}^{\mathbf r_i} \mathbf A\cdot d\boldsymbol\ell \right].

Reversing the hopping orientation complex conjugates the phase. Products around a closed oriented loop encode the magnetic flux with the same charge and orientation convention.

The Peierls substitution is an approximation tied to localized orbitals and slowly varying fields. Intra-cell dipoles, orbital magnetic moments, and field-induced changes of the basis may require additional terms.

The default is to keep ℏ\hbar and kBk_B explicit. SI units are used for dimensionful electromagnetic formulas unless an article declares Gaussian, atomic, natural, or lattice units.

Use:

  • EE or ε\varepsilon for energy;
  • ω\omega for angular frequency;
  • ν\nu for ordinary frequency only when needed;
  • β=1/(kBT)\beta=1/(k_BT) for inverse energy.

The conversions are

E=ℏω=hν,ω=2πν.E = \hbar\omega = h\nu, \qquad \omega = 2\pi\nu.

A spectrum labeled in electronvolts is energy resolved even if the horizontal variable is traditionally called ω\omega in the literature. State once whether ω\omega is angular frequency or an energy-valued variable used with ℏ=1\hbar=1.

μ\mu denotes chemical potential. EFE_F denotes the zero-temperature chemical potential when that limit is well defined:

EF=lim⁡T→0μ(T).E_F = \lim_{T\to0}\mu(T).

At finite temperature, prefer μ\mu. A band energy crossing μ\mu need not imply a sharp quasiparticle if interactions or disorder broaden the spectrum.

For fermions,

f(E)=1eβ(E−μ)+1.f(E) = \frac{1} {e^{\beta(E-\mu)}+1}.

For bosons with an allowed chemical potential,

nB(E)=1eβ(E−μ)−1.n_B(E) = \frac{1} {e^{\beta(E-\mu)}-1}.

The energy argument and chemical-potential convention must match the Hamiltonian. If time evolution uses the grand Hamiltonian K=H−μN^K=H-\mu\hat N, say so explicitly.

The default filling is electrons per primitive cell:

ν=NeNc.\nu = \frac{N_{\mathrm e}}{N_c}.

This differs from:

n=NeV,n = \frac{N_{\mathrm e}}{V},

the physical number density, and from occupation per orbital,

nˉ=NeNcNorb.\bar n = \frac{N_{\mathrm e}} {N_cN_{\mathrm{orb}}}.

Spin, valley, layer, and orbital degeneracies are not silently included. A phrase such as “half filling” must state the denominator. One electron per site can mean half filling for one spinful orbital, full filling for one spinless orbital, or another fraction in a multiorbital basis.

Use doping with an explicit reference:

δν=ν−νref.\delta\nu = \nu-\nu_{\mathrm{ref}}.

State whether positive δν\delta\nu means electron doping. Hole doping is often denoted p>0p>0 with ν=νref−p\nu=\nu_{\mathrm{ref}}-p; never assume the sign from the letter alone.

For an ordinary time-dependent quantity,

A(ω)=∫−∞∞dt eiωtA(t),A(t)=∫−∞∞dω2πe−iωtA(ω).\begin{aligned} A(\omega) &= \int_{-\infty}^{\infty} dt\, e^{i\omega t} A(t), \\ A(t) &= \int_{-\infty}^{\infty} \frac{d\omega}{2\pi} e^{-i\omega t} A(\omega). \end{aligned}

The Heisenberg operator is

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

With a perturbation

δH(t)=−f(t)B,\delta H(t) = -f(t)B,

the retarded response is

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

Then

δ⟨A⟩(ω)=χABR(ω)f(ω).\delta\langle A\rangle(\omega) = \chi_{AB}^{\mathrm R}(\omega) f(\omega).

If a source uses δH=+fB\delta H=+fB, its susceptibility or response equation carries the compensating sign.

For fermions, define

GαβR(k,t)=−iθ(t)⟨{ckα(t),ckβ†(0)}⟩.G_{\alpha\beta}^{\mathrm R}(\mathbf k,t) = -i\theta(t) \left\langle \left\{ c_{\mathbf k\alpha}(t), c_{\mathbf k\beta}^\dagger(0) \right\} \right\rangle.

When the spectral variable EE has energy units, use

GR(k,E)=1ℏ∫−∞∞dt eiEt/ℏGR(k,t).G^{\mathrm R}(\mathbf k,E) = \frac{1}{\hbar} \int_{-\infty}^{\infty} dt\, e^{iEt/\hbar} G^{\mathrm R}(\mathbf k,t).

The spectral function is

A(k,E)=−1πIm⁡GR(k,E).A(\mathbf k,E) = - \frac{1}{\pi} \operatorname{Im} G^{\mathrm R}(\mathbf k,E).

With canonical normalization,

∫−∞∞dE Aαα(k,E)=1.\int_{-\infty}^{\infty} dE\, A_{\alpha\alpha}(\mathbf k,E) = 1.

A source using energy-valued ω\omega has renamed EE and usually set ℏ=1\hbar=1. Check dimensions before comparing self-energies or linewidths.

For a local operator ORO_{\mathbf R}, define

Oq=1Nc∑Re−iq⋅ROR.O_{\mathbf q} = \frac{1}{\sqrt{N_c}} \sum_{\mathbf R} e^{-i\mathbf q\cdot\mathbf R} O_{\mathbf R}.

The default dynamical structure factor is

SOO(q,ω)=12π∫−∞∞dt eiωt⟨Oq(t)O−q(0)⟩.S_{OO}(\mathbf q,\omega) = \frac{1}{2\pi} \int_{-\infty}^{\infty} dt\, e^{i\omega t} \left\langle O_{\mathbf q}(t) O_{-\mathbf q}(0) \right\rangle.

It satisfies the equal-time sum rule

∫−∞∞dω SOO(q,ω)=⟨OqO−q⟩.\int_{-\infty}^{\infty} d\omega\, S_{OO}(\mathbf q,\omega) = \left\langle O_{\mathbf q} O_{-\mathbf q} \right\rangle.

Experimental cross sections can include 2π2\pi, ℏ\hbar, form-factor, polarization, unit-cell, and detailed-balance factors. Quote the measured normalization separately from the theoretical correlator.

When the operator has basis components, include embedding phases or form factors explicitly:

Oq=1Nc∑R,αe−iq⋅(R+τα)Fα(q)ORα.O_{\mathbf q} = \frac{1}{\sqrt{N_c}} \sum_{\mathbf R,\alpha} e^{-i\mathbf q\cdot (\mathbf R+\boldsymbol\tau_\alpha)} F_\alpha(\mathbf q) O_{\mathbf R\alpha}.

This physical probe phase is not removed merely because the Hamiltonian uses the cell Fourier convention.

Transport, Response, and Optics routes these shared conventions to the appropriate material or terminal framework.

Drude Theory and Boltzmann Transport apply the Fourier, signed-charge, tensor-order, electrochemical-field, and field-orientation conventions below to dc, thermal, optical, and Hall response.

Define conductivity by

Ji(ω)=∑jσij(ω)Ej(ω),J_i(\omega) = \sum_j \sigma_{ij}(\omega) E_j(\omega),

where J\mathbf J is current density. Thus σxy\sigma_{xy} means current along xx in response to electric field along yy.

Resistivity is the matrix inverse:

ρ=σ−1,\rho = \sigma^{-1},

not the componentwise reciprocal. In an isotropic two-dimensional system,

σ=(σxxσxy−σxyσxx).\sigma = \begin{pmatrix} \sigma_{xx} & \sigma_{xy}\\ -\sigma_{xy} & \sigma_{xx} \end{pmatrix}.

Then

ρxy=−σxyσxx2+σxy2.\rho_{xy} = - \frac{\sigma_{xy}} {\sigma_{xx}^2+\sigma_{xy}^2}.

This minus sign is a frequent source of apparent Hall-sign contradictions.

Use a right-handed laboratory frame:

x^×y^=z^.\hat{\mathbf x} \times \hat{\mathbf y} = \hat{\mathbf z}.

State the direction of magnetic field, current, voltage leads, magnetization, and sample normal. A reported sign without these directions is not portable.

With:

  • electron charge qe=−eq_{\mathrm e}=-e;
  • A=i⟨u∣∇ku⟩\mathcal A=i\langle u|\boldsymbol\nabla_{\mathbf k}u\rangle;
  • C=(2π)−1∫dkxdky ΩxyC=(2\pi)^{-1}\int dk_xdk_y\,\Omega_{xy};
  • σxy=Jx/Ey\sigma_{xy}=J_x/E_y;

the filled-band intrinsic contribution is

σxy=−Ce2h.\sigma_{xy} = - C \frac{e^2}{h}.

Equivalently, before a complete occupied band has been integrated,

σxyint=−e2ℏ∑n∫BZd2k(2π)2fn(k) Ωn,xy(k).\sigma_{xy}^{\mathrm{int}} = -\frac{e^2}{\hbar} \sum_n \int_{\mathrm{BZ}} \frac{d^2k}{(2\pi)^2} f_n(\mathbf k)\, \Omega_{n,xy}(\mathbf k).

The minus sign belongs to this complete convention block: it combines the electron current with the stated Berry-connection sign and tensor order. A source using A=−i⟨u∣∇ku⟩\mathcal A=-i\langle u|\boldsymbol\nabla_{\mathbf k}u\rangle reverses both CC and the written Chern-to-Hall bridge, leaving the physical response unchanged.

Equivalently,

σyx=Ce2h.\sigma_{yx} = C \frac{e^2}{h}.

Sources can report the opposite sign by reversing the Berry-connection sign, coordinate orientation, or tensor-index order. Compare the complete convention block, including the fixed physical charge, rather than treating carrier charge as an independent notational toggle.

Bulk electric polarization in a periodic crystal is defined modulo a polarization quantum. Shifting an orbital center by a lattice vector changes a cell dipole while leaving the infinite charge distribution equivalent.

In one dimension, moving charge qq by one lattice vector aa changes polarization by

ΔP=q\Delta P = q

when polarization is expressed as charge per cell boundary in the one-dimensional convention. In dd dimensions, a lattice-vector transport changes polarization by

ΔP=qRΩc.\Delta\mathbf P = \frac{q\mathbf R}{\Omega_c}.

State the origin, ionic contribution, electron-charge sign, occupied-band subspace, and branch. A Berry phase is not an absolute dipole without these data.

The canonical geometric treatment is Zak Phase Preview and the mathematical Berry pages. This conventions page does not choose a material topology claim: enter Topological Quantum Matter to select the appropriate live phase, response, boundary, or computation owner. Berry-Phase Polarization and Charge Pumping is the material owner for polarization branches, ionic-plus-electronic charge, and integrated current over an adiabatic cycle.

Use:

  • ⟨⋯ ⟩\langle\cdots\rangle for a quantum or thermal expectation;
  • (⋯ )‾\overline{(\cdots)} or [⋯ ]dis[\cdots]_{\mathrm{dis}} for disorder averaging;
  • ⟨⋯ ⟩cell\langle\cdots\rangle_{\mathrm{cell}} only after defining a unit-cell average;
  • η>0\eta>0 for a causal infinitesimal;
  • Γ>0\Gamma>0 for a physical linewidth or damping scale.

The operations

⟨A⟩H‾\overline{ \langle A\rangle_H }

and

⟨A⟩H‾\left\langle A \right\rangle_{\overline H}

are generally different. The first solves each disorder realization and then averages; the second replaces the Hamiltonian by an average before solving.

Do not use η\eta as a fitted lifetime. A numerical broadening should be reported with convergence checks and distinguished from physical self-energy.

State boundary conditions separately for each direction:

  • periodic or twisted for bulk calculations;
  • open for edges and finite samples;
  • lead-coupled for transport;
  • slab or semi-infinite for surfaces.

The following limits can fail to commute:

Nc→∞,ω→0,q→0,T→0,Γ→0,t→∞.\begin{gathered} N_c\to\infty, \qquad \omega\to0, \qquad \mathbf q\to0, \\ T\to0, \qquad \Gamma\to0, \qquad t\to\infty. \end{gathered}

For conductivity, specify whether

lim⁡ω→0lim⁡q→0\lim_{\omega\to0} \lim_{\mathbf q\to0}

or

lim⁡q→0lim⁡ω→0\lim_{\mathbf q\to0} \lim_{\omega\to0}

is intended, and whether the thermodynamic and clean limits occur first. The first often represents a spatially uniform time-dependent drive; the second can represent static long-wavelength equilibrium response.

For spontaneous order, a source field hh may require

lim⁡h→0lim⁡Nc→∞⟨M⟩h.\lim_{h\to0} \lim_{N_c\to\infty} \langle M\rangle_h.

Reversing the limits restores the symmetry in each finite system.

Computational Quantum Matter owns the calculation-specific record that applies these conventions to model, solver, convergence, benchmark, and probe validation; this page retains the shared sign, gauge, mesh, unit, and limit contracts.

Fractional reciprocal coordinates mean

k=∑iκibi.\mathbf k = \sum_i \kappa_i\mathbf b_i.

Cartesian coordinates use inverse-length units. Label axes explicitly; a number such as 0.50.5 is meaningless without the basis and units.

A band plot samples a one-dimensional path through a higher-dimensional Brillouin zone. It is not the full band structure or density of states. State:

  • reciprocal coordinates of every labeled point;
  • path order and interpolation;
  • whether horizontal distance is cumulative Euclidean distance;
  • energy zero and chemical potential;
  • spin, orbital, or symmetry projections;
  • whether degeneracies are exact or plotting-tolerance coincidences.

For a uniform mesh with weights wkw_{\mathbf k},

∑kwk=1\sum_{\mathbf k} w_{\mathbf k} = 1

for an average per cell. A Brillouin-zone integral is then approximated by

Ωc(2π)d∫BZddk F(k)≃∑kwkF(k).\frac{\Omega_c}{(2\pi)^d} \int_{\mathrm{BZ}} d^dk\, F(\mathbf k) \simeq \sum_{\mathbf k} w_{\mathbf k} F(\mathbf k).

Symmetry-reduced meshes require weights for the full star. Berry curvature near avoided crossings, Fermi-surface integrals, and singular densities of states need convergence beyond a visually smooth plot.

DifferenceTranslation
cell versus orbital Fourier phasemultiply orbital component by e−ik⋅ταe^{-i\mathbf k\cdot\boldsymbol\tau_\alpha}
e+iωte^{+i\omega t} versus e−iωte^{-i\omega t} forward transformreplace ω→−ω\omega\to-\omega and translate causal prescriptions
angular frequency versus energyE=ℏωE=\hbar\omega and transform measures accordingly
electron charge −e-e versus signed eereplace every Lorentz, minimal-coupling, and current factor consistently
$\mathcal A=+i\langle u\nabla u\rangleversusversus-i\langle u
σxy=Jx/Ey\sigma_{xy}=J_x/E_y versus Jy/ExJ_y/E_xswap tensor indices and account for antisymmetry
per-cell versus per-volume responsemultiply or divide by Ωc\Omega_c with operator normalization
S=ℏσ/2S=\hbar\sigma/2 versus dimensionless S=σ/2S=\sigma/2restore ℏ\hbar in operators and coupling units
h(k)h(\mathbf k) periodic versus quasi-periodicidentify whether basis embedding enters the Fourier phase

A reliable translation checks an invariant:

  • canonical commutator or anticommutator;
  • total state count;
  • Hermiticity;
  • gauge covariance;
  • an equal-time or spectral sum rule;
  • dimensions of the final observable.
  1. Including τα\boldsymbol\tau_\alpha in one Fourier transform but not its inverse.
  2. Calling k\mathbf k momentum while using mv=ℏkm\mathbf v=\hbar\mathbf k in a nonparabolic band.
  3. Comparing Berry signs without comparing reciprocal orientation and tensor order.
  4. Differentiating a k\mathbf k-dependent basis Hamiltonian without connection terms.
  5. Using NsN_s for sites, cells, and orbitals in the same calculation.
  6. Calling a finite-temperature chemical potential the Fermi energy without qualification.
  7. Treating resistivity as 1/σij1/\sigma_{ij} component by component.
  8. Using energy-valued ω\omega in one equation and angular frequency in another.
  9. Applying a Peierls phase with the electron charge sign but the reverse link orientation.
  10. Taking a dc, clean, or thermodynamic limit silently.
  11. Plotting a high-symmetry line and inferring the absence of off-path crossings.
  12. Calling numerical η\eta a physical lifetime without convergence evidence.

Before accepting a calculation, verify:

  • primitive and reciprocal vectors obey ai⋅bj=2πδij\mathbf a_i\cdot\mathbf b_j=2\pi\delta_{ij};
  • real-space and reciprocal orientations are stated;
  • every basis position τα\boldsymbol\tau_\alpha is declared;
  • lattice Fourier phases and normalization form an inverse pair;
  • the Bloch Hamiltonian follows from the declared hopping orientation;
  • state counting gives NcNorbN_cN_{\mathrm{orb}} modes before internal multiplicities;
  • cell and orbital Fourier conventions are not mixed;
  • Bloch eigenvectors use a stated cell inner product;
  • Berry and Hall signs include charge, orientation, and tensor order;
  • energy and angular frequency are dimensionally distinct;
  • response specifies the perturbation sign and current normalization;
  • broadening and disorder averages are defined;
  • every noncommuting limit is written in order;
  • numerical meshes reproduce state counts and sum rules.

Use the default Fourier transform to show that one orbital on each of NcN_c cells produces exactly NcN_c independent crystal-momentum modes.

Solution

The transform matrix is

UkR=1Nce−ik⋅R.U_{\mathbf k\mathbf R} = \frac{1}{\sqrt{N_c}} e^{-i\mathbf k\cdot\mathbf R}.

Discrete orthogonality gives

∑RUkRUk′R∗=δk,k′.\sum_{\mathbf R} U_{\mathbf k\mathbf R} U^\ast_{\mathbf k'\mathbf R} = \delta_{\mathbf k,\mathbf k'}.

The inverse relation gives

∑kUkR∗UkR′=δR,R′.\sum_{\mathbf k} U^\ast_{\mathbf k\mathbf R} U_{\mathbf k\mathbf R'} = \delta_{\mathbf R,\mathbf R'}.

Thus UU is an Nc×NcN_c\times N_c unitary matrix. It maps NcN_c real-space modes to NcN_c momentum modes without changing the Hilbert-space dimension.

Starting from cˉk=D(k)ck\bar{\mathbf c}_{\mathbf k}=D(\mathbf k)\mathbf c_{\mathbf k}, derive the relation between h(k)h(\mathbf k) and hˉ(k)\bar h(\mathbf k).

Solution

The Hamiltonian in the cell basis is

H=∑kck†h(k)ck.H = \sum_{\mathbf k} \mathbf c_{\mathbf k}^\dagger h(\mathbf k) \mathbf c_{\mathbf k}.

Since

ck=D†(k)cˉk,\mathbf c_{\mathbf k} = D^\dagger(\mathbf k) \bar{\mathbf c}_{\mathbf k},

substitution gives

H=∑kcˉk†D(k)h(k)D†(k)cˉk.H = \sum_{\mathbf k} \bar{\mathbf c}_{\mathbf k}^\dagger D(\mathbf k) h(\mathbf k) D^\dagger(\mathbf k) \bar{\mathbf c}_{\mathbf k}.

Therefore

hˉ(k)=D(k)h(k)D†(k).\bar h(\mathbf k) = D(\mathbf k) h(\mathbf k) D^\dagger(\mathbf k).

The spectra agree because the transformation is unitary at each k\mathbf k.

For

H=+t∑j(cj†cj+1+cj+1†cj),H = +t \sum_j \left( c_j^\dagger c_{j+1} +c_{j+1}^\dagger c_j \right),

find the dispersion and compare its band minimum with the convention using −t-t.

Solution

The Fourier transform gives

ε(k)=+2tcos⁡(ka).\varepsilon(k) = +2t\cos(ka).

For t>0t>0, its minimum occurs at k=π/ak=\pi/a, whereas the Hamiltonian with coefficient −t-t gives

ε(k)=−2tcos⁡(ka)\varepsilon(k) = -2t\cos(ka)

and has its minimum at k=0k=0. On a bipartite nearest-neighbor chain the two signs can be related by the gauge transformation cj↦(−1)jcjc_j\mapsto(-1)^jc_j, which shifts momentum by π/a\pi/a. Additional hoppings, boundaries, or flux can make hopping signs physically meaningful.

Derive the transformation of

Ai=i⟨u∣∂iu⟩\mathcal A_i = i\langle u|\partial_i u\rangle

under ∣u′⟩=eiχ∣u⟩|u'\rangle=e^{i\chi}|u\rangle.

Solution

Differentiate:

∂i∣u′⟩=eiχ(i(∂iχ)∣u⟩+∂i∣u⟩).\partial_i|u'\rangle = e^{i\chi} \left( i(\partial_i\chi)|u\rangle + \partial_i|u\rangle \right).

Using ⟨u′∣=⟨u∣e−iχ\langle u'|=\langle u|e^{-i\chi} and ⟨u∣u⟩=1\langle u|u\rangle=1,

Ai′=i⟨u∣[i(∂iχ)∣u⟩+∂i∣u⟩]=−∂iχ+Ai.\begin{aligned} \mathcal A_i' &= i \langle u| \left[ i(\partial_i\chi)|u\rangle + \partial_i|u\rangle \right] \\ &= -\partial_i\chi + \mathcal A_i. \end{aligned}

The curl is unchanged because mixed derivatives commute on a smooth patch:

Ωij′=Ωij.\Omega_{ij}' = \Omega_{ij}.

Invert

σ=(σxxσxy−σxyσxx)\sigma = \begin{pmatrix} \sigma_{xx} & \sigma_{xy}\\ -\sigma_{xy} & \sigma_{xx} \end{pmatrix}

and verify the sign of ρxy\rho_{xy}.

Solution

The determinant is

det⁡σ=σxx2+σxy2.\det\sigma = \sigma_{xx}^2+\sigma_{xy}^2.

Therefore

ρ=σ−1=1σxx2+σxy2(σxx−σxyσxyσxx).\rho = \sigma^{-1} = \frac{1} {\sigma_{xx}^2+\sigma_{xy}^2} \begin{pmatrix} \sigma_{xx} & -\sigma_{xy}\\ \sigma_{xy} & \sigma_{xx} \end{pmatrix}.

Hence

ρxy=−σxyσxx2+σxy2.\rho_{xy} = - \frac{\sigma_{xy}} {\sigma_{xx}^2+\sigma_{xy}^2}.

The sign follows from matrix inversion, not from a new carrier convention.

Show that the Peierls-dressed term

tije(iq/ℏ)∫jiA⋅dℓci†cjt_{ij} e^{(iq/\hbar)\int_j^i\mathbf A\cdot d\boldsymbol\ell} c_i^\dagger c_j

is invariant under the electromagnetic gauge transformation defined above.

Solution

The link phase changes as

exp⁡[iqℏ∫ji(A+∇Λ)⋅dℓ]=eiq(Λi−Λj)/ℏexp⁡[iqℏ∫jiA⋅dℓ].\exp \left[ \frac{iq}{\hbar} \int_j^i (\mathbf A+\boldsymbol\nabla\Lambda) \cdot d\boldsymbol\ell \right] = e^{iq(\Lambda_i-\Lambda_j)/\hbar} \exp \left[ \frac{iq}{\hbar} \int_j^i \mathbf A\cdot d\boldsymbol\ell \right].

The annihilation operator transforms as

cj⟶eiqΛj/ℏcj,c_j \longrightarrow e^{iq\Lambda_j/\hbar}c_j,

so

ci†cj⟶e−iqΛi/ℏeiqΛj/ℏci†cj.c_i^\dagger c_j \longrightarrow e^{-iq\Lambda_i/\hbar} e^{iq\Lambda_j/\hbar} c_i^\dagger c_j.

The operator phase cancels the added link phase. Reversing the link orientation reverses the line integral and complex conjugates the phase.

  • Quantum Matter Overview routes readers to the definition, conceptual map, convention ledger, and next technical chapter.
  • Choosing a Model for Quantum Matter selects the smallest adequate retained description before this page fixes its calculation contract.
  • How to Use This Volume provides goal-based routes and identifies when this convention ledger is needed as a lookup.
  • Quantum Matter Reference and Data applies this page’s convention system to transferable lookup records; this page remains the canonical owner of the lattice, basis, Fourier, Berry, electromagnetic, normalization, and response conventions themselves.
  • Conventions Overview owns sitewide state, operator, unit, and Fourier choices.
  • Many-Body Symbols and Conventions owns generic Fock-space, ensemble, correlator, and limit notation.
  • Phonons applies the cell-phase, angular-frequency, energy, Brillouin-zone, and per-cell normalization contracts to lattice dynamics.
  • Drude Theory applies the signed-charge, Fourier, conductivity-tensor, Hall, and optical conventions.
  • Boltzmann Transport applies the phase-space, degeneracy, collision-integral, electrochemical-field, and thermoelectric conventions.
  • Hall Effect fixes the transverse tensor order, voltage polarity, carrier sign, dimensionality, and resistivity-to-conductivity inversion.
  • Notation Collisions catalogs overloaded symbols.
  • Bloch’s Theorem owns the theorem, proof, finite-size state counting, and Bloch–Floquet qualification.
  • Wannier Functions applies the Bloch-frame, Fourier-phase, orbital-embedding, and reciprocal-sewing conventions fixed here to the construction and localization of band-derived orbitals.
  • Symmetry of Bloch States applies the momentum, reciprocal-sewing, basis, and band-index conventions fixed here to little-group representations and symmetry-constrained degeneracies.
  • Bloch Theorem formula card gives the compact translation result.
  • Berry Curvature formula card gives gauge and projector formulas.
  • Chern Number formula card records normalization and isolation assumptions.
  • Retarded and Advanced Response owns causal analyticity and Fourier structure.
  • Peierls Phase Preview develops lattice gauge covariance.
  • Quantum Matter Map shows where each convention enters model and experiment.
  1. N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston, 1976. Standard lattice, reciprocal-space, Bloch, band, and transport notation.
  2. C. Kittel, Introduction to Solid State Physics, 8th ed., Wiley, 2004. A compact reference for crystal and reciprocal-lattice conventions.
  3. G. D. Mahan, Many-Particle Physics, 3rd ed., Springer, 2000. Green-function, response, spectral, and transport conventions.
  4. A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, Dover, 2003. Operator normalization, propagators, and response functions.
  5. D. Vanderbilt, Berry Phases in Electronic Structure Theory, Cambridge University Press, 2018. Modern Bloch gauge, polarization, Wannier, and Berry conventions.
  6. D. Xiao, M.-C. Chang, and Q. Niu, “Berry Phase Effects on Electronic Properties,” Reviews of Modern Physics 82, 1959–2007 (2010), doi:10.1103/RevModPhys.82.1959. Berry-curvature and semiclassical sign conventions with explicit electron charge.
  7. J. M. Blount, “Formalisms of Band Theory,” Solid State Physics 13, 305–373 (1962), doi:10.1016/S0081-1947(08)60459-2. A foundational treatment of position, derivatives, and gauge in Bloch representations.
  8. R. Resta, “Macroscopic Polarization in Crystalline Dielectrics: The Geometric Phase Approach,” Reviews of Modern Physics 66, 899–915 (1994), doi:10.1103/RevModPhys.66.899. Establishes branch and polarization-quantum conventions.
  9. R. Peierls, Quantum Theory of Solids, Oxford University Press, 1955. Classical source for lattice models and magnetic phases on hopping amplitudes.