Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Statistical Mechanics

Fermi-Dirac Distribution

The grand canonical potential for non-interacting fermions:

Ω[β,μ]=i1βlog(N=01eβ(NϵiNμ))=\Omega[\beta, \mu] = -\sum_i {1\over\beta} \log\left(\sum_{N=0}^1 e^{-\beta\left(N\epsilon_i - N\mu\right)}\right) =
=i1βlog(1+eβ(ϵiμ))== -\sum_i {1\over\beta} \log\left(1 + e^{-\beta\left(\epsilon_i - \mu\right)}\right) =
=iΩi,= \sum_i \Omega_i\,,

Where Ωi\Omega_i is the grand canonical potential for a single particle state:

Ωi=1βlog(1+eβ(ϵiμ)).\Omega_i = - {1\over\beta} \log\left(1 + e^{-\beta\left(\epsilon_i - \mu\right)}\right)\,.

The total (average) number of particles in the state ii is:

Ni=(Ωiμ)T,V=N_i = -\left(\partial\Omega_i\over\partial\mu\right)_{T,V} =

=(1βlog(1+eβ(ϵiμ)))μ== -{\partial \left(- {1\over\beta} \log\left(1 + e^{-\beta\left(\epsilon_i - \mu\right)}\right) \right)\over\partial\mu} =

=1β11+eβ(ϵiμ)eβ(ϵiμ)β== {1\over\beta} {1\over 1 + e^{-\beta(\epsilon_i - \mu)}} e^{-\beta(\epsilon_i - \mu)} \beta =

=1eβ(ϵiμ)+1.= {1\over e^{\beta(\epsilon_i - \mu)} + 1} \,.

The total (average) number of particles in the full system is:

N=(Ωμ)T,V=(iΩiμ)T,V=i((Ωiμ)T,V)=N = -\left(\partial\Omega\over\partial\mu\right)_{T,V} = -\left(\partial\sum_i\Omega_i\over\partial\mu\right)_{T,V} = \sum_i \left(-\left(\partial\Omega_i\over\partial\mu\right)_{T,V}\right) =

=iNi=i1eβ(ϵiμ)+1.= \sum_i N_i = \sum_i {1\over e^{\beta(\epsilon_i - \mu)} + 1} \,.

Bose-Einstein Distribution

The grand canonical potential for non-interacting bosons:

Ω[β,μ]=i1βlog(1eβ(ϵiμ))=iΩi,\Omega[\beta, \mu] = \sum_i {1\over\beta} \log\left(1 - e^{-\beta\left(\epsilon_i - \mu\right)}\right) = \sum_i \Omega_i\,,

Where Ωi\Omega_i is the grand canonical potential for a single particle state:

Ωi=1βlog(1eβ(ϵiμ)).\Omega_i = {1\over\beta} \log\left(1 - e^{-\beta\left(\epsilon_i - \mu\right)}\right)\,.

The total (average) number of particles in the state ii is:

Ni=(Ωiμ)T,V=N_i = -\left(\partial\Omega_i\over\partial\mu\right)_{T,V} =

=(1βlog(1eβ(ϵiμ)))μ== -{\partial \left({1\over\beta} \log\left(1 - e^{-\beta\left(\epsilon_i - \mu\right)}\right) \right)\over\partial\mu} =

=1β11eβ(ϵiμ)(eβ(ϵiμ))β== -{1\over\beta} {1\over 1 - e^{-\beta(\epsilon_i - \mu)}} (-e^{-\beta(\epsilon_i - \mu)}) \beta =

=1eβ(ϵiμ)1.= {1\over e^{\beta(\epsilon_i - \mu)} - 1} \,.

The total (average) number of particles in the full system is:

N=(Ωμ)T,V=(iΩiμ)T,V=i((Ωiμ)T,V)=N = -\left(\partial\Omega\over\partial\mu\right)_{T,V} = -\left(\partial\sum_i\Omega_i\over\partial\mu\right)_{T,V} = \sum_i \left(-\left(\partial\Omega_i\over\partial\mu\right)_{T,V}\right) =

=iNi=i1eβ(ϵiμ)1.= \sum_i N_i = \sum_i {1\over e^{\beta(\epsilon_i - \mu)} - 1} \,.