How do Fermi-Dirac statistics contribute to fermions in the Pauli Exclusion Principle? 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
 A: 
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:

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.

A: 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
A: 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.
