Stirling's approximation

Stirling's approximation is named after the Scottish mathematician James Stirling (1692-1770).

$\left.\ln N!\right.=\ln 1+\ln 2+\ln 3+...+\ln N=\sum _{k=1}^{N}\ln k.$

using Euler-MacLaurin formula one has

$\sum _{k=1}^{N}\ln k=\int _{1}^{N}\ln x\,dx+\sum _{k=1}^{p}{\frac {B_{2k}}{2k(2k-1)}}\left({\frac {1}{n^{2k-1}}}-1\right)+R,$

where B₁ = −1/2, B₂ = 1/6, B₃ = 0, B₄ = −1/30, B₅ = 0, B₆ = 1/42, B₇ = 0, B₈ = −1/30, ... are the Bernoulli numbers, and R is an error term which is normally small for suitable values of p.

Then, for large N,

$\ln N!\sim \int _{1}^{N}\ln x\,dx\sim N\ln N-N.$

after some further manipulation one arrives at (apparently Stirling's contribution was the prefactor of ${\sqrt {2\pi }}$)

$N!={\sqrt {2\pi N}}\;N^{N}e^{-N}e^{\lambda _{N}}$

where

${\frac {1}{12N+1}}<\lambda _{N}<{\frac {1}{12N}}.$

For example:

N	N! (exact)	N! (Stirling)	Error (%)
1	1	0.92213700	8.44
2	2	1.91900435	4.22
3	6	5.83620959	2.81
4	24	23.5061751	2.10
5	120	118.019168	1.67
6	720	710.078185	1.40
7	5040	4980.39583	1.20
8	40320	39902.3955	1.05
9	362880	359536.873	0.93
10	3628800	3598695.62	0.84

When one is dealing with numbers of the order of the Avogadro constant ($10^{23}$) this formula is essentially exact. In computer simulations the number of atoms or molecules (N) is invariably greater than 100; for N=100 the percentage error is approximately 0.083%.

Gosper’s formula

Gosper’s formula :

$n!\approx {\sqrt {2\pi \left(n+{\frac {1}{6}}\right)}}\;\left({\frac {n}{e}}\right)^{n}$

Which results in:

N	N! (exact)	N! (Gosper)
1	1	0.99602180
2	2	1.99736305
3	6	5.99613535
4	24	23.9908895
5	120	119.970030
6	720	719.872829
7	5040	5039.33747
8	40320	40315.9028
9	362880	362850.646
10	3628800	3628560.82

Applications in statistical mechanics

• Ideal gas Helmholtz energy function