So I’m doing some research on Quantum Chromodynamics and I've come to an area in which I’m having some trouble arriving at a conclusive answer. Apparently the Pauli exclusion principle states that no two fermions can have the same quantum numbers, which is supposedly the same as saying that no two fermions can be in the same quantum state, (are those effectively saying the same thing)? And fermions are particles with half integer spins that obey Fermi-Dirac statistics. Please explain what Fermi-Dirac statistics are, I’m having trouble with contradictory or confusing answers in other sites. Help would be appreciated! Very important: The fermions I am talking about are quarks, not electrons

  • $\begingroup$ Fermions have half-integer spin. $\endgroup$
    – ProfRob
    Aug 2, 2021 at 23:59
  • $\begingroup$ This question has nothing to do with the strong force or QCD but ok, maybe people that know more than me can make connections to those specifically. $\endgroup$
    – SuperCiocia
    Aug 3, 2021 at 3:36

3 Answers 3


no two fermions can have the same quantum numbers, which is supposedly the same as saying that no two fermions can be in the same quantum state, (are those effectively saying the same thing)?

Yes they mean the same thing. A quantum state is uniquely determined by a set of quantum numbers.

The fermions I am talking about are quarks, not electrons

Does not matter. Electrons and quarks are both fermions. The Pauli exclusion principle will apply to two identical electrons or to two identical quarks. It will not apply between an electron and a quark as they are distinct particles (distinct mass).

And fermions are particles with half integer spins that obey Fermi-Dirac statistics. Please explain what Fermi-Dirac statistics are,

So there are three levels to this:

1. Bosons vs fermions.

If you have two particles ($1$ and $2$) with wavefunctions $\phi$ and $\psi$, you have to describe the system as a superposition of [case where particle $1$ has wavefunction $\phi$ AND particle $2$ has wavefunction $\psi$] AND vice versa.

Physically, this is because the labelling is done by us the experimenter and thus should not have any effect on the system itself and its evolution. Elementary particles are indistinguishable so we cannot label them once and expect them to "keep" their label.

So the total wavefunction $\Psi(1,2)$ will be: $$ \Psi(1,2) = \phi(1)\psi(2) \pm \phi(2)\psi(1). $$ For identical particles $\psi = \phi$ and hence the $\pm$ determines whether or not the two particles can coexist. The plus sign case means they can co-exist and they are called bosons, while the minus sign gives you a zero and hence they cannot co-exist, and are called fermions.

A neat example is considering two electrons in two states of a square well, and you can see that for the $-$ sign ("anti-symmetric" because $\Psi$ picks up an overall minus sign upon particle exchange $1\leftrightarrow$2) the probability of finding the particles at the same position $x=y$ is always exactly zero:

enter image description here

The argument above is dumbed down a bit but can be rigourously derived from the topology of particle exchange.

The above requirements are usually summarised as commutation and anti-commutation relationships:

\begin{equation} \begin{gathered} \text{bosons} \Leftrightarrow \left [ a,a^\dagger \right ] = 0, \\ \text{fermions} \Leftrightarrow \left \{ a,a^\dagger \right \} = 0. \end{gathered} \end{equation}

2. Spin-statistics theorem.

The spin-statistics theorem, essentially (after lots of maths and proofs and stuff) entails these relations for the creation ($a^\dagger$) and annihilation ($a$) operators of particles (more correctly fields) with angular momentum $j$:

\begin{equation} \begin{gathered} \text{for bosons: } \quad \left [ a, a^\dagger \right ] \propto 1 - (-1)^{2j} , \\ \text{for fermions: } \quad \left \{ a, a^\dagger\right \} \propto 1 - (-1)^{2j}, \end{gathered} \end{equation}

and, in order to match these to the relations above, you get that fermions must satisfy $2j \in \mathbb{Z} \implies j \in \mathbb{Z}/2$. Hence fermions have half-integer spin, because $j$, in the absence of any other (orbital) angular momentum, is just the spin.

3. Particle distribution.

Any single-particle (non interacting) system is described by the following thermal density matrix, that is Boltzmann distributed:

\begin{equation} \mathrm{e}^{-\beta (H-\mu N)} = \sum_i\mathrm{e}^{-\beta (E_i-\mu n_i)}\,|n_i\rangle\langle n_i|, \end{equation}

because it's single-particle, the number operator $\hat N$ gives you $n_i|n_i\rangle$ and the energy operator (Hamiltonian) $\hat H|n_i\rangle$ gives you $E_i |n_i\rangle$.

Now I want to define the average occupation of a state, that is (particles)/(states), $f(E)$: \begin{equation} f = \bar{n} = \frac{1}{\text{total}}\sum_j n_j P(n_j) = \frac{1}{{\sum_j \mathrm{e}^{-n_j(E_j - \mu)/ k_B T}}}\sum_j \langle n_j| \,\mathrm{e}^{-(H - \mu N)/k_B T} n |n_j \rangle = \\ \frac{1}{Z}\sum_j n_j P(n_j), \end{equation} where \begin{equation} P(n) = \frac{1}{Z}\, \mathrm{e}^{-n(E-\mu)/k_B T} \quad \text{with} \quad Z = \sum_n \mathrm{e}^{-n(E-\mu)/k_BT} = \sum_n \left ( \mathrm{e}^{-r} \right )^{n}, \end{equation} $P$ is the probability of finding a particle at energy $E$ (equivalent to number $n$), and $Z$ the partition function.

We said that fermions cannot co-exist, so the sum over the states has to be over $0$ or $1$, because anything $\geq 2$ will give you zero net wavefunction (see first bullet point). This restricts the sum $Z= \sum_{n=0}^1 r^n = 1 + r$ and results in: \begin{equation} \begin{gathered} \bar{n} = \sum_{n=0}^1 n\,P(n) = \frac{1}{Z}\,\mathrm{e}^{-(E-\mu)/k_BT} = \frac{r}{1+r} \\ \Rightarrow \bar{n} = f_{\text{FD}} = \frac{1}{\mathrm{e}^{(E-\mu)/k_B T} + 1}. \end{gathered} \end{equation}

So this is where the Fermi-Dirac (FD) distribution function comes from.

It is derived from the Pauli exclusion principle, and it is hence compatible with it. If you plot it (below, for some non-zero temperature), you see that $f$ is at most one for fermions: you can have maximum $1$ particle per quantum state.

enter image description here


I think the Fermi-Dirac statistic is referring to the Fermi-Dirac distribution function:

$$f(E) = \frac{1}{1+\exp\left(\frac{E-\mu}{kT} \right)},$$

and this function does not exceed one, meaning you never have more than one particle in a certain quantum state, unlike in the case of the Bose-Einstein distribution function.

The quantum state is always defined in terms of quantum numbers, and the spin is always one of the well defined quantum numbers.

Do you want to know how comes that half-integer spin is related to such a distribution function?

In general, half-integer spin -> Pauli exclusion principle -> Fermi-Dirac distribution function.

The relationship half-integer spin -> Pauli exclusion principle is described here: https://en.wikipedia.org/wiki/Spin-statistics_theorem

The link Pauli exclusion principle -> Fermi-Dirac distribution function is described here: https://nanohub.org/resources/5787/download/2009.02.02-ECE606-L9.pdf

  • $\begingroup$ Thanks, I appreciate it! Do you mind particularly explaining what the Fermi-Dirac distribution function dictates or explains? $\endgroup$ Aug 3, 2021 at 0:08
  • $\begingroup$ There are several derivations of the Fermi-Dirac distribution functions from the Pauli exclusion principle. Take a look here nanohub.org/resources/5787/download/2009.02.02-ECE606-L9.pdf $\endgroup$
    – freude
    Aug 3, 2021 at 0:09
  • $\begingroup$ Your questions should be asked other way around: how does the Pauli exclusion principle contribute to the Fermi-Dirac distribution function, since the first one is more fundamental. $\endgroup$
    – freude
    Aug 3, 2021 at 0:11
  • $\begingroup$ Much appreciated, thank you! $\endgroup$ Aug 3, 2021 at 0:13
  • $\begingroup$ half-integer spin -> Pauli exclusion principle -> Fermi-Dirac distribution function $\endgroup$
    – freude
    Aug 3, 2021 at 0:19

This is a really deep and interesting topic! I hope this answer helps to clear up some confusions.

To address your first question, quantum numbers label the different states of the system, so if two particles have the same quantum numbers, they are in the same state.

Fermi-Dirac statistics refer to the phase acquired by the wavefunction when two particles are exchanged. In particular, if you exchange two fermions, the wavefunction should acquire an overall minus sign. Suppose we have a system described by a single quantum number, $n$, and it contains two fermions with position vectors $\mathbf{x}_1$ and $\mathbf{x}_2$ respectively. Let these fermions be in states $n_1$ and $n_2$. Then the wavefunction for the combined system should be: $$\Psi (\mathbf{x}_1,\mathbf{x}_2) = \psi_{n_1}(\mathbf{x}_1)\psi_{n_2}(\mathbf{x}_2) - \psi_{n_1}(\mathbf{x}_2)\psi_{n_2}(\mathbf{x}_1)$$ where the states $\psi_{n_1}(\mathbf{x})$ and $\psi_{n_2}(\mathbf{x})$ are single-particle states. You can see from the two-particle state that $\Psi (\mathbf{x}_1,\mathbf{x}_2) = -\Psi (\mathbf{x}_2,\mathbf{x}_1)$, which encodes the fact that swapping the two fermions results in an overall change of sign. In contrast, two-particle bosonic states obey $\Psi(\mathbf{x}_1,\mathbf{x}_2) = \Psi(\mathbf{x}_2,\mathbf{x}_1)$.

The Pauli exclusion principle states that no two indistinguishable fermions can occupy the same quantum state. We can see that this holds by considering the case that $n_1=n_2$ above. In this situation, $\Psi = 0$ identically. That means there is zero probability for the fermions to be in the same state. Note this is not the case for bosons.

The spin-statistics theorem tells us that the phase acquired by the wavefunction when two particles are exchanged is the same as the phase acquired when one particle undergoes a full $2\pi$ rotation. So fermions pick up a minus sign when you rotate one of them by $2\pi$. This is telling us that they must have half-integer spin, which is not at all trivial! I would suggest reading up more on spin specifically if you're interested in this connection.

Important: The above facts are true for any type of (indistinguishable) fermions. Electrons and quarks obey the same rules in this sense. However, Fermi-Dirac statistics (referring to the minus sign picked up by the wavefunction upon particle exchange) are not the same as the Fermi-Dirac distribution. The FD distribution is obeyed by non-interacting fermions that are in thermal equilibrium with a reservoir with which they can exchange particles. In particular, quarks are very strongly interacting particles and so would not obey the FD distribution. Exchange statistics (and the Pauli-exclusion principle) are a more general, deeper fact.

  • $\begingroup$ Thank you! I appreciate it! $\endgroup$ Aug 3, 2021 at 2:07
  • $\begingroup$ Thank you once again. Two other websites I visited were saying or implying that Fermi-Dirac statistics were simply a WAY particles could be distributed or that the RESULT of Fermi-Dirac statistics was actually Fermi-Dirac distribution. In quarks specifically, you are saying this is incorrect? $\endgroup$ Aug 3, 2021 at 15:29

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.