量子統計力学に出てくる特殊関数-その1 ガンマ関数の定義と諸性質

本記事は、 統合ゼミのアドベントカレンダー4日目 の一部です。

量子統計力学に出てくる特殊関数-その1 ガンマ関数の定義と諸性質

さて、まずはじめはガンマ関数についてみていきたいと思います！

ガンマ関数の定義と階乗

まずは、ガンマ関数の定義です。

この定義式を部分積分すると、
$$\begin{array}{rl}\mathrm{\Gamma }(n)& ={\int }_0^{\mathrm{\infty }}{e}^{-t}{t}^{n-1}\mathrm{d}t\\ & ={\int }_0^{\mathrm{\infty }}{e}^{-t}{\left(\frac{1}{n}{t}^n\right)}^{\prime }\mathrm{d}t\\ & ={\left[e^{-t}\frac{1}{n}{t}^n\right]}_0^{\mathrm{\infty }}+\frac{1}{n}{\int }_0^{\mathrm{\infty }}{e}^{-t}{t}^n\mathrm{d}t\\ & =\frac{1}{n}\mathrm{\Gamma }(n+1)\end{array}$$
となります。つまり、$\mathrm{\Gamma }(n+1)$と$\mathrm{\Gamma }(n)$の間には次の関係があることがわかります。

$n$が整数の時、これを繰り返し用いると、
$$\begin{array}{rl}\mathrm{\Gamma }(n+1)& =n\mathrm{\Gamma }(n)\\ & =n(n-1)\mathrm{\Gamma }(n-1)\\ & ⋮\\ & =n(n-1)(n-2)\cdots 1\times \mathrm{\Gamma }(1)\end{array}$$
となります。ここで、ガンマ関数の定義式において$n=1$とすると、
$$\begin{array}{rl}\mathrm{\Gamma }(1)& ={\int }_0^{\mathrm{\infty }}{e}^{-t}{t}^{1-1}\mathrm{d}t\\ & ={\int }_0^{\mathrm{\infty }}{e}^{-t}\mathrm{d}t\\ & =-{\left[e^{-t}\right]}_0^{\mathrm{\infty }}\\ & =1\end{array}$$
であるので、
$$\begin{array}{r}\mathrm{\Gamma }(n+1)=n(n-1)(n-1)\cdots 1\times 1.\end{array}$$
よって、$\mathrm{\Gamma }(n+1)$は階乗の定義そのものになります。なんかごつい名前の関数と慣れ親しんだ階乗がこうして結びつくのは感動ものですね。

また、ガンマ関数の定義式を用いると、階乗の積分表示を得ることができます。

半整数値を変数とするガンマ関数

次に半整数値を変数とするガンマ関数についてみていこうと思います！
はじめに、$n$が半整数値のときの【公式1】の初期値$\mathrm{\Gamma }(\frac{1}{2})$を計算してみましょう。
$$\begin{array}{r}\mathrm{\Gamma }\left(\frac{1}{2}\right)={\int }_0^{\mathrm{\infty }}{e}^{-t}{t}^{-\frac{1}{2}}\mathrm{d}t\end{array}$$
積分変数を$t=x^2$と変換すれば、この積分はガウス積分となります。
$$\begin{array}{rl}\mathrm{\Gamma }\left(\frac{1}{2}\right)& ={\int }_0^{\mathrm{\infty }}{e}^{-x^2}{x}^{2(-\frac{1}{2})}\cdot 2x\mathrm{d}x\\ & =2{\int }_0^{\mathrm{\infty }}{e}^{-x^2}\mathrm{d}x\\ & ={\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}{e}^{-x^2}\mathrm{d}x\\ & =\sqrt{\pi }\end{array}$$
これと、【公式1】を繰り返し用いることで、$n$が半整数値のときのガンマ関数を次のように求めることができます。
$$\begin{array}{rl}\mathrm{\Gamma }\left(\frac{3}{2}\right)& =\mathrm{\Gamma }(\frac{1}{2}+1)=\frac{1}{2}\mathrm{\Gamma }\left(\frac{1}{2}\right)=\frac{1}{2}\sqrt{\pi }\\ \mathrm{\Gamma }\left(\frac{5}{2}\right)& =\mathrm{\Gamma }(\frac{3}{2}+1)=\frac{3}{2}\mathrm{\Gamma }\left(\frac{3}{2}\right)=\frac{3}{2}\frac{1}{2}\sqrt{\pi }\\ \mathrm{\Gamma }\left(\frac{7}{2}\right)& =\mathrm{\Gamma }(\frac{5}{2}+1)=\frac{5}{2}\mathrm{\Gamma }\left(\frac{5}{2}\right)=\frac{5}{2}\frac{3}{2}\frac{1}{2}\sqrt{\pi }\\ \mathrm{\Gamma }\left(\frac{9}{2}\right)& =\mathrm{\Gamma }(\frac{7}{2}+1)=\frac{7}{2}\mathrm{\Gamma }\left(\frac{7}{2}\right)=\frac{7}{2}\frac{5}{2}\frac{3}{2}\frac{1}{2}\sqrt{\pi }\\ & ⋮\end{array}$$
一般に、$n=(2l+1)/2(l=0,1,2,3,\cdots )$のときは、
$$\begin{array}{rl}\mathrm{\Gamma }\left(\frac{2l+1}{2}\right)& =\frac{(2l-1)!!}{2^l}\sqrt{\pi }\end{array}$$
となります。ここで、$\{\}!!$という記号は、2重階乗を表しており、整数$m$に対して、次のように定義されています。ただし、$(-1)!!\equiv 1$です。
$$\begin{array}{r}m!!=\{\begin{array}{ll}m(m-2)(m-4)\cdots 3\cdot 1& m\text{:奇数}\\ m(m-2)(m-4)\cdots 4\cdot 2& m\text{:偶数}\end{array}\end{array}$$
通常の階乗だけを用いて、一般の場合を表すことを考えます。2重階乗の定義式において、$m=2l$とすると、$(2l)!!=2^{l}l!$と表すことができるので、これを用いると、
$$\begin{array}{rl}\mathrm{\Gamma }\left(\frac{2l+1}{2}\right)& =\frac{(2l-1)!!}{2^l}\frac{(2l)!!}{2^{l}l!}\sqrt{\pi }\\ & =\frac{(2l)!}{2^{2l}l!}\sqrt{\pi }\end{array}$$
と表すことができました。

$N$次元球の体積と表面積

さて、ガンマ関数がよく出てくる例ってなんですかね？私が思うに$N$次元球の話だと思っています。
そこで、最後は半径$r$の$N$次元球の体積$V_N$と表面積$S_N$を求めてみたいと思います！

体積$V_N$は、どうせ$r^N$に比例することが予想されるので、比例定数を$C_N$とすると、
$$\begin{array}{r}{V}_N=C_N{r}^N\end{array}$$
と表すことができます。また、表面積$S_N$はこれを$r$で微分することにより、
$$\begin{array}{r}{S}_N=\frac{\mathrm{d}{V}_N}{\mathrm{d}r}=C_{N}N{r}^{N-1}\end{array}$$
と表すことができます。

次はこの$C_N$を計算すれば良いのですが。。。おやおや、何やら統合稲荷神社の方から御神託が聞こえて来ましたよ。

かわいいきつねさんが言うので間違いありません、とりあえず$N$重ガウス積分
$$\begin{array}{r}{I}_N={\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}{\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\cdots {\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\exp \{-(r_1^2+r_2^2+\cdots r_N^2)\}\mathrm{d}{r}_1\mathrm{d}{r}_2\cdots \mathrm{d}{r}_N\end{array}$$
を計算していきたいと思います！まずは、直接積分してやりましょう。
被積分関数の指数関数はバラすことができるので、変数分離することができるので、
$$\begin{array}{rl}{I}_N& ={\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\exp (-r_1^2)\mathrm{d}{r}_1{\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\exp (-r_2^2)\mathrm{d}{r}_2\cdots {\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\exp (-r_N^2)\mathrm{d}{r}_N\\ & ={\{{\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\exp (-r^2)\mathrm{d}r\}}^N\\ & ={\left(\mathrm{\Gamma }\left(\frac{1}{2}\right)\right)}^N\\ & ={\pi }^{N/2}\end{array}$$
となります。別の方法として、$N$次元空間の極座標を導入して計算する方法があります。このとき、被積分関数は$r^2=r_1^2+r_2^2\cdots r_N^2$で与えられる$r$のみの関数となります。角度変数に関する積分は、$S_N=\frac{\mathrm{d}{V}_N}{\mathrm{d}r}=C_{N}N{r}^{N-1}$を用いると、
$$\begin{array}{rl}{I}_N& ={\int }_0^{\mathrm{\infty }}\exp (-r^2)S_N\mathrm{d}r\\ & ={\int }_0^{\mathrm{\infty }}\exp (-r^2)C_{N}N{r}^{N-1}\mathrm{d}r\\ & =C_{N}N{\int }_0^{\mathrm{\infty }}\exp (-r^2)r^{N-1}\mathrm{d}r\end{array}$$
となります。ここで、積分変数を$r^2=t$と変換すると、
$$\begin{array}{rl}{I}_N& =C_{N}N{\int }_0^{\mathrm{\infty }}\exp (-r^2)r^{N-1}\mathrm{d}r\\ & =C_{N}N{\int }_0^{\mathrm{\infty }}\exp (-t)t^{(N-2)/2}\frac{\mathrm{d}t}{2}\\ & =C_N\frac{N}{2}{\int }_0^{\mathrm{\infty }}\exp (-t)t^{(N-2)/2}\mathrm{d}t\\ & =C_N\frac{N}{2}\mathrm{\Gamma }\left(\frac{N}{2}\right)\end{array}$$
となります。これらの別々の方法によって得られた結果を等しいとおくと、
$$\begin{array}{r}{C}_N\frac{N}{2}\mathrm{\Gamma }\left(\frac{N}{2}\right)={\pi }^{N-2}\end{array}$$
よって、求めたい比例定数は
$$\begin{array}{r}{C}_N=\frac{{\pi }^{N/2}}{\frac{N}{2}\mathrm{\Gamma }\left(\frac{N}{2}\right)}=\frac{2{\pi }^{N/2}}{N\mathrm{\Gamma }\left(\frac{N}{2}\right)}=\frac{{\pi }^{N/2}}{\mathrm{\Gamma }(\frac{N}{2}+1)}\end{array}$$
となります。（御神託のおかげだなあ（しみじみ）。）
よって、$N$次元球の体積$V_N$と表面積$S_N$は

と表すことができます。

さて、ガンマ関数はまだまだ色々と性質がありますが今回はこの辺で終えておこうと思います。
それでは、次の記事に参りましょう！