- スターリングの公式の異なる2つの表式の関係．
統計力学でよく出てくるStirlingの公式について考察します．
スターリングの公式の2つの表式
スターリングの公式（Stirling's approximation/Stirling's formula）は，ガンマ関数（階乗）に関する公式です．物理の教科書では\begin{align}
&\log N! \\
&=N(\log N-1) +O(\log N) \tag{1}\label{eq:st_phys}\\
&\qquad\qquad\qquad\qquad\qquad(N\rightarrow\infty)
\end{align}と書かれることが多く，数学の教科書では
\begin{align}
\lim_{s\to \infty}\frac{\Gamma(s+1)}{\sqrt{2\pi} s^{s+1/2}e^{-s}}
=1
\tag{2}\label{eq:st_math}
\end{align}と書かれることが多いです．
ここで，式 (\ref{eq:st_phys})の$O(\cdot)$はLandauの記号と呼ばれるものです ．詳しくは次の記事を参照して下さい：
式 (\ref{eq:st_phys})の導出
式 (\ref{eq:st_phys})を示しましょう．以下に再掲します：&\log N! \\
&=N(\log N-1) +O(\log N) \\
&\qquad\qquad\qquad\qquad\qquad(N\rightarrow\infty)
\end{align}
【証明】
よく$\log x$のグラフを図示して説明されるように，以下の不等式
\begin{align}
& \int_1^N\log x\,\mathrm{d}x \\
& \qquad \leq\log N!=\sum_{n=1}^N\log n \\
& \qquad\qquad \leq\int_1^N\log x\,\mathrm{d}x+\log N
\end{align}が成立します．
これは，（よく図を書いてされる説明からわかるように）ほとんどRiemann積分の定義から示されます：
区間$I$上でRiemann可積分な関数$f(x)$は
\begin{align}
& \sum_{k} \inf\{f(x)|x\in I_k\} \\
& \qquad \leq \int_I f(x)\,\mathrm{d}x
\leq \sum_{k} \sup\{f(x)|x\in I_k\}
\end{align}を満たします（$\{I_k\}_k$は$I$の分割）．これを$f(x)=\log x$に適用すれば，最初の不等式より
\begin{align}
\int_1^N \log x\,\mathrm{d}x
\geq \sum_{n=1}^{N-1} \min_{[n,n+1]}\log x
=\sum_{n=1}^{N-1}\log n
\end{align}が，２つ目の不等式より
\begin{align}
\int_1^N \log x\,\mathrm{d}x
\leq \sum_{n=1}^{N-1} \max_{[n,n+1]}\log x
=\sum_{n=1}^{N}\log n
\end{align}がわかります．
部分積分により
\begin{align}
\int_1^N\log x\,\mathrm{d}x
=N(\log N-1)+1
\end{align}なので，上の不等式は
\begin{align}
&\frac{1}{\log N} \\
&\qquad \leq\frac{\log N!-N(\log N-1)}{\log N} \\
&\qquad\qquad \leq\frac{1}{\log N}+1
\end{align}と変形できます．
この式から，任意の$M>0$をfixするとき，$N>M$において
\begin{align}
\left|\frac{\log N!-N(\log N-1)}{\log N}\right|
<\frac{1}{\log M}+1
\end{align}が成立することがわかりました．
よって，$O(\cdot)$の定義（ランダウの記号の使い方 - Notes_JP参照）から，式 (\ref{eq:st_phys})が成立します．//
2式の関係について
式 (\ref{eq:st_phys})と式 (\ref{eq:st_math})の関係について調べてみましょう．式 (\ref{eq:st_math})は
\begin{align}
\Gamma(s+1)
\sim \sqrt{2\pi} s^{s+1/2}e^{-s}
\qquad(s\rightarrow\infty)
\end{align}と表すことができますが，この式で$\log$をとった式である
\begin{align}
& \log \Gamma(s+1) \\
&\qquad \sim\log \left[\sqrt{2\pi} s^{s+1/2}e^{-s}\right]
\qquad(s\rightarrow\infty)
\end{align}が成立することがわかります．
【証明】
$f(s)=\Gamma(s+1)$, $g(s)=\sqrt{2\pi} s^{s+1/2}e^{-s}$とします．
$\displaystyle\lim_{s\to\infty} \left(f(s)/g(s)\right)=1$と，$\displaystyle\lim_{s\to\infty}|\log g(s)|=\infty$により
\begin{align}
&\left|\frac{\log f(s)}{\log g(s)}-1\right| \\
&=\left|\frac{\log \left|f(s)/g(s)\right|}{\log g(s)}\right|
\rightarrow 0 \qquad(s\rightarrow\infty).
\end{align}よって，$\log f(s)\sim \log g(s)$が言えました．//
示した式から，$s=N$（自然数）の場合に式 (\ref{eq:st_phys})に帰着することがわかります．