Fermion doubling

The fermion doubling problem is a problem that is encountered when naively trying to put fermionic fields on a lattice. It consists in the appearance of spurious states, such that one ends up having 2^d fermionic particles (with d the number of discretized dimensions) for each original fermion. In order to solve this problem, several strategies are in use, such as and staggered fermions.

Mathematical overview[]

The action of a free Dirac fermion in d dimensions,^{[note 1]} of mass m, and in the continuum (i.e. without discretization) is commonly given as

${\displaystyle S=\int d^{d}x\;{\bar {\psi }}(i\partial \!\!\!/-m)\psi \;.}$

Here, the Feynman slash notation was used to write

${\displaystyle \gamma _{\mu }a^{\mu }=a\!\!\!/}$

where ${\displaystyle \gamma _{\mu }}$ are the gamma matrices. When this action is discretized on a cubic lattice, the fermion field ${\displaystyle \psi (x)}$ is replaced with a discretized version ${\displaystyle \psi _{x}}$, where ${\displaystyle x}$ now denotes the lattice site. The derivative is replaced by the finite difference. The resulting action is now:^[1]

${\displaystyle S=a^{d}\sum _{x,\mu }{\frac {i}{2a}}({\bar {\psi }}_{x}\gamma _{\mu }\psi _{x+{\hat {\mu }}}-{\bar {\psi }}_{x+{\hat {\mu }}}\gamma _{\mu }\psi _{x})-a^{d}\sum _{x}m{\bar {\psi }}_{x}\psi _{x}\;,}$

where ${\displaystyle a}$ is the lattice spacing and ${\displaystyle {\hat {\mu }}}$ is the vector of length ${\displaystyle a}$ in the ${\displaystyle \mu }$ direction. When one computes the inverse fermion propagator in momentum space, one readily finds:^[1]

${\displaystyle S^{-1}(p)=m+{\frac {i}{a}}\sum _{\mu }\gamma _{\mu }\sin(p_{\mu }a)\;.}$

Due to the finite lattice spacing the momenta ${\displaystyle p_{\mu }}$ have to be inside the (first) Brillouin zone, which is typically taken to be the interval ${\displaystyle [-\pi /a,\pi /a[}$.

When simply taking the limit ${\displaystyle a\to 0}$ in the above inverse propagator, one recovers the correct continuum result. However, when instead expanding this expression around a value of ${\displaystyle p_{\mu }}$ where one or more of the components are at the corners of the Brillouin zone (i.e. equal to ${\displaystyle \pi /a}$), one finds the same continuum form again, although the sign in front of the gamma matrix can change.^{[2][note 2]} This means that, when one of the components of the momentum is near ${\displaystyle \pi /a}$, the discretized fermion field will again behave like a continuum fermion. This can happen with all d components of the momentum, leading to —with the original fermion with momentum near the origin included— 2^d different "tastes" (in analogy to flavor).^{[note 3]}

The Nielsen–Ninomiya theorem[]

Nielsen and Ninomiya proved a theorem^[3] stating that a local, real, free fermion lattice action, having chiral and translational invariance, necessarily has fermion doubling. The only way to get rid of the doublers is by violating one of the presuppositions of the theorem —for example:

• explicitly violate chiral symmetry, giving an infinitely high mass to the doublers which then decouple.
• So-called "" have a nonlocal action.
• Staggered fermions
• Interacting fermions ^{[4][5][6][7][8]}

See also[]

• Staggered fermions: a way to reduce the number of doublers
• Acoustic and optical phonons: a similar phenomenon in solid state crystals

Notes and references[]

Notes

References