n次元球面の体積(volume of the n-sphere)とは、単位球面 $S^n\subset\mathbb R^{n+1}$ の Riemann 計量による $n$ 次元の体積 $\lvert S^n\rvert$ のことで、$\lvert S^n\rvert=2\pi^{(n+1)/2}/\Gamma(\frac{n+1}2)$ である。Stokes の定理から単位球の体積 $\omega_{n+1}$ との関係 $\lvert S^n\rvert=(n+1)\omega_{n+1}$ が出て、球対称な関数の積分公式を Gauss 積分に使うと値が決まる。単位球の体積は 5 次元で、球面の体積は $S^6$ で最大になり、次元を上げると $0$ に近づく。
前提知識: n次元球面, 変数変換公式(多変数積分), Γ関数, Stokesの定理
単位円周 $S^1$ の長さは $2\pi$、単位球面 $S^2$ の面積は $4\pi$ である。一般の $n$ 次元球面 $S^n\subset\mathbb R^{n+1}$ についても、その $n$ 次元の体積(「表面積」)を Γ 関数で表す公式がある:
$$
\lvert S^n\rvert=\frac{2\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+1}2\bigr)},\qquad \lambda(B^{n+1})=\frac{\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+3}2\bigr)}=\frac{\lvert S^n\rvert}{n+1}.
$$
右の式は、内部の球の体積との関係「体積を半径で微分すると表面積になる」を表す。
この記事では、球面の体積を Riemann 計量から定め、上の公式を 2 つの事実から導く。1 つは Stokes の定理による $\lvert S^n\rvert=(n+1)\lambda(B^{n+1})$ であり、もう 1 つは球対称な関数の積分を 1 変数の積分に直す公式である。後者を Gauss 積分 $\int_{\mathbb R^{n+1}}e^{-\lvert x\rvert^2}dx=\pi^{(n+1)/2}$ に使うと値が決まる。値を並べると、球の体積は 5 次元で、球面の体積は $S^6$ で最大になり、次元が上がると $0$ に近づく。
$n\ge0$ とする。$\lvert x\rvert$ は Euclid ノルム、$\lambda$ は Lebesgue 測度で、
$$
S^n:=\{x\in\mathbb R^{n+1}\mid\lvert x\rvert=1\},\quad B^{n+1}:=\{x\mid\lvert x\rvert<1\},\quad D^{n+1}:=\{x\mid\lvert x\rvert\le1\}
$$
とおく。$\omega_{n+1}:=\lambda(B^{n+1})$ を単位球の体積とする。$r>0$ について $\lambda(rB^{n+1})=r^{n+1}\omega_{n+1}$ である(変数変換公式(多変数積分) の命題「線形写像による体積の変化」を $x\mapsto rx$ に使う)。$D^{n+1}\setminus B^{n+1}=S^n$ は、すべての $\varepsilon\in(0,1)$ について $(1+\varepsilon)B^{n+1}\setminus(1-\varepsilon)B^{n+1}$ に含まれ、その体積 $\bigl((1+\varepsilon)^{n+1}-(1-\varepsilon)^{n+1}\bigr)\omega_{n+1}$ は $\varepsilon\to0$ で $0$ に近づくので、$\lambda(S^n)=0$、$\lambda(D^{n+1})=\omega_{n+1}$ である。
$n\ge1$ のとき、$S^n$ は $\mathbb R^{n+1}$ のコンパクトな $n$ 次元部分多様体であり、$\mathbb R^{n+1}$ の内積から Riemann 計量を受け継ぐ。その Riemann 体積(Riemann多様体 の定義「リーマン体積」)による $S^n$ 全体の体積を $\lvert S^n\rvert$ と書き、$n$ 次元球面の体積 という。$n=0$ では $S^0=\{\pm1\}$ の点の個数をとり、$\lvert S^0\rvert:=2$ とする。
$n=1$ なら $\lvert S^1\rvert$ は単位円周の長さ、$n=2$ なら単位球面の面積である。半径 $r$ の球面 $S^n(r):=rS^n$ は、拡大 $x\mapsto rx$ によって計量を $r^2$ 倍した $S^n$ と同じ Riemann 多様体になるので、$\lvert S^n(r)\rvert=r^n\lvert S^n\rvert$ である(Riemann多様体 の命題「計量の一様な拡大と共形変形」)。
$\mathbb R^{n+1}$ の座標を $x_1,\dots,x_{n+1}$ とし、$\Omega:=dx_1\wedge\cdots\wedge dx_{n+1}$、$E:=\sum_ix_i\,\partial/\partial x_i$ とおく。$\Omega$ に $E$ を代入した $n$ 形式を
$$
\sigma:=\iota_E\Omega=\sum_{i=1}^{n+1}(-1)^{i-1}x_i\,dx_1\wedge\cdots\wedge\widehat{dx_i}\wedge\cdots\wedge dx_{n+1}
$$
とする($\widehat{\ \ }$ はその因子を除くことを表す)。$n=1$ なら $\sigma=x\,dy-y\,dx$ である。
$n\ge0$ について $\lvert S^n\rvert=(n+1)\,\omega_{n+1}$ である。
$n=0$ では $\lvert S^0\rvert=2=\lambda((-1,1))$ である。$n\ge1$ とし、$S^n$ には閉球体 $D^{n+1}$ の境界としての向き(外向きの法線を先頭に置く向き。Stokesの定理 の定義「向き付けられた境界付き多様体と境界の向き」)を入れる。
$d\sigma=(n+1)\Omega$。 $d\bigl(x_i\,dx_1\wedge\cdots\widehat{dx_i}\cdots\wedge dx_{n+1}\bigr)=dx_i\wedge dx_1\wedge\cdots\widehat{dx_i}\cdots\wedge dx_{n+1}=(-1)^{i-1}\Omega$ なので、$d\sigma=\sum_i(-1)^{i-1}(-1)^{i-1}\Omega=(n+1)\Omega$ である。
$S^n$ 上で $\sigma$ は体積形式。 $p\in S^n$ で $E_p=p$ は長さ $1$ の外向き法線である。$T_pS^n$ の正の向きの正規直交基底 $(v_1,\dots,v_n)$ をとると、$(p,v_1,\dots,v_n)$ は $\mathbb R^{n+1}$ の正の向きの正規直交基底なので、$\sigma_p(v_1,\dots,v_n)=\Omega(p,v_1,\dots,v_n)=1$ である。一方、Riemann 体積形式 $\operatorname{vol}=\sqrt{\det(g_{ij})}\,dy^1\wedge\cdots\wedge dy^n$($y$ は正の向きの局所座標)も、$v_k=\sum_ia_{ik}\,\partial/\partial y^i$ と書くと、正規直交性 $a^{\mathsf T}(g_{ij})a=I$ から $(\det a)^2\det(g_{ij})=1$、$\det a>0$ なので $\operatorname{vol}(v_1,\dots,v_n)=\sqrt{\det(g_{ij})}\det a=1$ である。$n$ 次元空間の上の $n$ 次交代形式は 1 つの基底での値で決まるので、$S^n$ への引き戻しについて $\sigma|_{S^n}=\operatorname{vol}$ であり、$\int_{S^n}\sigma=\lvert S^n\rvert$ である。
Stokes の定理。 $D^{n+1}$ はコンパクトな境界付き多様体で、境界は $S^n$ である(Stokesの定理 の節「設定:境界付き多様体・向き・積分」)。Stokesの定理 の定理「Stokes の定理」により
$$
\lvert S^n\rvert=\int_{S^n}\sigma=\int_{D^{n+1}}d\sigma=(n+1)\int_{D^{n+1}}\Omega=(n+1)\,\lambda(D^{n+1})=(n+1)\,\omega_{n+1}
$$
である。$\square$
この命題は「$\lambda(rB^{n+1})=r^{n+1}\omega_{n+1}$ を $r$ で微分すると $r^n\lvert S^n\rvert$」という関係の $r=1$ の場合である。
$n\ge0$ とし、$g\colon(0,\infty)\to[0,\infty]$ を Borel 可測関数とする。このとき
$$
\int_{\mathbb R^{n+1}}g(\lvert x\rvert)\,dx=\lvert S^n\rvert\int_0^\infty g(r)\,r^n\,dr
$$
が成り立つ(両辺は $\infty$ でもよく、左辺の被積分関数は原点で $0$ とする)。$g$ が実数値で、$\int_0^\infty\lvert g(r)\rvert r^n\,dr<\infty$ のときも同じ等式が成り立つ。
$(0,\infty)$ の Borel 集合 $A$ について
$$
\nu_1(A):=\lambda\bigl(\{x\ne0\mid\lvert x\rvert\in A\}\bigr),\qquad\nu_2(A):=(n+1)\,\omega_{n+1}\int_Ar^n\,dr
$$
とおく。$x\mapsto\lvert x\rvert$ は連続なので $\nu_1$ は測度であり、$\nu_2$ も測度である。$0\le a< b$ について $\{a<\lvert x\rvert\le b\}=bD^{n+1}\setminus aD^{n+1}$($a=0$ では原点を除く)なので
$$
\nu_1((a,b])=(b^{n+1}-a^{n+1})\,\omega_{n+1}=\nu_2((a,b])
$$
である。区間 $(a,b]$($0\le a< b$)と空集合の全体は π 系で、$(0,\infty)$ の Borel 集合全体を生成し、増加列 $(0,k]$ は $(0,\infty)$ を覆って $\nu_1((0,k])<\infty$ である。Dynkinのπ-λ定理 の定理「π系の上で一致する測度の一意性」の 2 により $\nu_1=\nu_2$ である。
$g=\chi_A$ なら主張は $\nu_1(A)=\nu_2(A)$ そのもの(prop-nsv-stokes で $(n+1)\omega_{n+1}=\lvert S^n\rvert$)である。両辺は $g$ について線形なので $0$ 以上の単関数でも成り立ち、単関数の増加列 $g_k\uparrow g$ に Lebesgue積分 の定理「単調収束定理」を両辺で使えば一般の $g\ge0$ で成り立つ。実数値の場合は、$\lvert g\rvert$ に使って両辺の可積分性を得たうえで、$g$ の正の部分と負の部分に分けて引けばよい。$\square$
$n\ge0$ について
$$
\lvert S^n\rvert=\frac{2\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+1}2\bigr)},\qquad\omega_{n+1}=\frac{\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+3}2\bigr)}
$$
である。半径 $r$ の球面と球については、それぞれ $r^n$ 倍、$r^{n+1}$ 倍である。
$I:=\int_{\mathbb R^{n+1}}e^{-\lvert x\rvert^2}\,dx$ を 2 通りに計算する。$e^{-\lvert x\rvert^2}=\prod_ie^{-x_i^2}$ なので、Tonelli の定理(Fubiniの定理 の定理「Tonelliの定理」)と 1 変数の値 $\int_{\mathbb R}e^{-t^2}dt=\sqrt\pi$(Gauss積分 の定理「Gauss積分の値」)から $I=\pi^{(n+1)/2}$ である。一方、prop-nsv-radial を $g(r)=e^{-r^2}$ に使い、$s=r^2$ と置換すると($r^n\,dr=\frac12s^{(n-1)/2}\,ds$)、
$$
I=\lvert S^n\rvert\int_0^\infty e^{-r^2}r^n\,dr=\frac{\lvert S^n\rvert}2\int_0^\infty s^{\frac{n+1}2-1}e^{-s}\,ds=\frac{\lvert S^n\rvert}2\,\Gamma\Bigl(\frac{n+1}2\Bigr)
$$
である(Γ関数 の定義「$\Gamma$ 関数」)。両者を比べて $\lvert S^n\rvert$ の式を得る。prop-nsv-stokes と $\Gamma(s+1)=s\Gamma(s)$(Γ関数 の命題「関数等式と階乗」)を $s=\frac{n+1}2$ に使うと
$$
\omega_{n+1}=\frac{\lvert S^n\rvert}{n+1}=\frac{2\pi^{(n+1)/2}}{(n+1)\Gamma\bigl(\frac{n+1}2\bigr)}=\frac{\pi^{(n+1)/2}}{\Gamma\bigl(\frac{n+3}2\bigr)}
$$
である。半径 $r$ の場合は定義の節で見た $r^n$、$r^{n+1}$ 倍である。$\square$
球の体積の式は Gauss積分 の定理「単位球の体積」にもあり、そこでは層ごとの積分 $\int f=\int_0^\infty\lambda(\{f>s\})\,ds$ で示されている。Axl26 5.44(p. 141)は次元を 2 つ下げる漸化式で、同じ値を Γ 関数を使わない形で示している。上の証明はどちらとも違い、球面の体積を先に求めて、球の体積を prop-nsv-stokes から得ている。
$\Gamma(\frac12)=\sqrt\pi$、$\Gamma(1)=1$ と関数等式から、次の値を得る。$\lvert S^n\rvert$ は $\pi^{\lceil n/2\rceil}$ の有理数倍である。
| $n$ | $\lvert S^n\rvert$ | 近似値 | $\omega_{n+1}=\lambda(B^{n+1})$ | 近似値 |
|---|---|---|---|---|
| $0$ | $2$ | $2$ | $2$ | $2$ |
| $1$ | $2\pi$ | $6.2832$ | $\pi$ | $3.1416$ |
| $2$ | $4\pi$ | $12.566$ | $\frac43\pi$ | $4.1888$ |
| $3$ | $2\pi^2$ | $19.739$ | $\frac12\pi^2$ | $4.9348$ |
| $4$ | $\frac83\pi^2$ | $26.319$ | $\frac8{15}\pi^2$ | $5.2638$ |
| $5$ | $\pi^3$ | $31.006$ | $\frac16\pi^3$ | $5.1677$ |
| $6$ | $\frac{16}{15}\pi^3$ | $33.073$ | $\frac{16}{105}\pi^3$ | $4.7248$ |
| $7$ | $\frac13\pi^4$ | $32.470$ | $\frac1{24}\pi^4$ | $4.0587$ |
$k\ge1$ について $\omega_k:=\lambda(B^k)=\pi^{k/2}/\Gamma(\frac k2+1)$ と書き、$\omega_0:=1$ とおく(同じ式で $k=0$ とした値)。
1:主定理の式と $\Gamma(\frac k2+1)=\frac k2\Gamma(\frac k2)$ から $\omega_k/\omega_{k-2}=\pi\,\Gamma(\frac k2)/\Gamma(\frac k2+1)=2\pi/k$ である($k=2$ では $\omega_2=\pi=\frac{2\pi}2\omega_0$)。
2:主定理から $\lvert S^{n+2}\rvert=2\pi^{(n+3)/2}/\Gamma(\frac{n+3}2)=2\pi\,\omega_{n+1}$ である。
3:1 により、$k$ と $k-2$ の比は $2\pi/k$ で、$k\le6$ なら $1$ より大きく、$k\ge7$ なら $1$ より小さい。よって偶数の $k$ では $\omega_0<\omega_2<\omega_4<\omega_6>\omega_8>\cdots$、奇数の $k$ では $\omega_1<\omega_3<\omega_5>\omega_7>\cdots$ である。$\omega_5=\frac8{15}\pi^2>\omega_6=\frac16\pi^3$($\pi<\frac{16}5$ による)なので、最大は $\omega_5$ だけである。2 により $\lvert S^n\rvert=2\pi\,\omega_{n-1}$($n\ge2$)であり、$\lvert S^0\rvert,\lvert S^1\rvert<2\pi\omega_5$ なので、$\lvert S^n\rvert$ の最大は $n-1=5$ のときだけである。
4:$k>4\pi$ なら $\omega_k<\frac12\omega_{k-2}$ なので、$\omega_k\to0$ である。2 により $\lvert S^n\rvert\to0$ である。$\square$
2 は Archimedes の「球面の面積は外接円柱の側面積に等しい」の高次元版にあたる($n=0$ では $\lvert S^2\rvert=2\pi\cdot2=4\pi$)。
$0<\varepsilon<1$ とすると、半径 $1-\varepsilon$ の球と単位球の体積の比は $(1-\varepsilon)^k$ なので、殻 $B^k\setminus(1-\varepsilon)B^k$ が占める割合 $1-(1-\varepsilon)^k$ は $k\to\infty$ で $1$ に近づく。たとえば $\varepsilon=0.01$、$k=1000$ なら割合は $1-0.99^{1000}>0.9999$ である。また単位球を囲む立方体 $[-1,1]^k$ の体積は $2^k$ で、比 $\omega_k/2^k$ は 4 により $0$ に近づく。高次元では、球は立方体の中で無視できるほど小さい。
$\mathbb R^{n+1}\setminus\{0\}$ の上で $\tau:=\sigma/\lvert x\rvert^{n+1}$ とおく。$n=1$ なら $\tau=(x\,dy-y\,dx)/(x^2+y^2)$ は外微分 の例「原点を除いた平面の角度形式」である。
$n\ge1$ とする。$\tau$ は $\mathbb R^{n+1}\setminus\{0\}$ 上の閉形式で、$\int_{S^n}\tau=\lvert S^n\rvert>0$ である。したがって $\tau$ は完全形式でなく、$H^n_{dR}(\mathbb R^{n+1}\setminus\{0\})\ne0$ である。
1 形式 $\alpha$ について $\alpha\wedge\Omega=0$($n+2$ 次)なので、$0=\iota_E(\alpha\wedge\Omega)=\alpha(E)\Omega-\alpha\wedge\sigma$、すなわち $\alpha\wedge\sigma=\alpha(E)\Omega$ である。$\alpha=\sum_jx_j\,dx_j$ とすると $\alpha(E)=\lvert x\rvert^2$ である。$d(\lvert x\rvert^{-(n+1)})=-(n+1)\lvert x\rvert^{-(n+3)}\alpha$ なので、prf-nsv-stokes の $d\sigma=(n+1)\Omega$ とあわせて
$$
d\tau=-(n+1)\lvert x\rvert^{-(n+3)}\lvert x\rvert^2\,\Omega+\lvert x\rvert^{-(n+1)}(n+1)\Omega=0
$$
である。$S^n$ の上では $\lvert x\rvert=1$ なので $\tau=\sigma$ であり、prf-nsv-stokes により $\int_{S^n}\tau=\lvert S^n\rvert>0$ である。$\tau=d\eta$ なら、包含 $\iota\colon S^n\to\mathbb R^{n+1}\setminus\{0\}$ による引き戻しは $\iota^*\tau=d(\iota^*\eta)$ で、Stokesの定理 の系「閉多様体上の積分と完全形式」によりその積分は $0$ になり、矛盾する。$\square$
同じ計算から、$r:=\lvert x\rvert$ について $dr\wedge(\sigma/r)=\Omega$ も成り立つ($dr=\alpha/r$)。$\sigma/r$ を半径 $r$ の球面に制限したものはその体積形式であり、この等式は「体積要素 = 半径方向の長さ × 球面の体積要素」を形式の言葉で述べている。prop-nsv-radial はその測度論での表現である。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| prop-nsv-radial の「被積分関数が $\lvert x\rvert$ だけで決まる」 | $f(x)=\chi_{\{x_1>0\}}(x)\,e^{-\lvert x\rvert^2}$ | $\int f=\lvert S^n\rvert\int_0^\infty f(re_1)r^n\,dr$(左辺 $\frac12\pi^{(n+1)/2}$、右辺 $\pi^{(n+1)/2}$) |
| prop-nsv-radial の「$g\ge0$ または $\int_0^\infty\lvert g\rvert r^n\,dr<\infty$」 | $g(r)=r^{-(n+1)}\sin r$ | 左辺が Lebesgue 積分として定まること(右辺の広義積分は $\lvert S^n\rvert\frac\pi2$ に収束する) |
| ノルムが Euclid ノルム | $\lvert x\rvert_\infty=\max_i\lvert x_i\rvert$ の単位球 $(-1,1)^{n+1}$ | $\omega_{n+1}=\pi^{(n+1)/2}/\Gamma(\frac{n+3}2)$(体積は $2^{n+1}$) |
1 行目:半空間 $\{x_1>0\}$ と $\{x_1<0\}$ は鏡映 $x_1\mapsto-x_1$ で移りあい、境界の超平面は零集合なので、左辺は $I=\pi^{(n+1)/2}$ の半分である。右辺は $f(re_1)=e^{-r^2}$ を使った主定理の証明の計算そのもので $I$ に等しい。
2 行目:$g(r)r^n=\frac{\sin r}r$ で、$\int_0^\infty\frac{\sin r}r\,dr=\frac\pi2$ は広義積分として収束する。しかし $\int_0^\infty\frac{(\sin r)^\pm}{r}\,dr=\infty$($(\sin r)^+$ は $k$ が偶数の区間 $[k\pi,(k+1)\pi]$ で、$(\sin r)^-$ は $k$ が奇数の区間で、それぞれ寄与が $\frac{2}{(k+1)\pi}$ 以上であり、その和は発散する)なので、prop-nsv-radial を $g^\pm$ に使うと $\int g(\lvert x\rvert)^\pm\,dx=\infty$ であり、$g(\lvert x\rvert)$ の Lebesgue 積分は $\infty-\infty$ になって定まらない。
3 行目:$\{\lvert x\rvert_\infty<1\}=(-1,1)^{n+1}$ の体積は $2^{n+1}$ で、$n\ge1$ では $\omega_{n+1}$ と異なる($n=1$ で $4\ne\pi$)。$n=0$ では両方 $2$ である。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する