数学・物理通信
7巻3号

編集

2017 年 6 月

新関章三・矢野2017 年 6 月 1 日

忠1目次 (Contents)
1. 変数を含むコタンジェントの冪乗の有限和
中西

襄

2

世戸憲治

12

伊藤仁之

20

矢野

忠

22

新関章三

30

Noboru NAKANISHI

2

2. 量子力学における周期ポテンシャル問題 (5)3. 超弦理論と Dirac のハサミ

4. n 次元の球の体積 1

5. 編集後記

1. Finite Sums of Powers of the Cotangents Involving a Variable

2. A Periodic Potential Problem in Quantum Mechanics (5)Kenji SETO

12

Hitoshi ITO

20

Tadashi YANO

22

Shozo NIIZEKI

30

3. Superstring Theory and Dirac’s Scissors

4. The Volume of n-dimensinal Hypersphere 1

5. Editorial Comments

変数を含むコタンジェントの冪乗の有限和
Finite Sums of Powers of the Cotangents Involving a Variable中西襄1

Noboru NAKANISHI

1

はじめに「コサインの正負冪乗の有限和とチェビシェフ多項式」（「数学・物理通信」6-7）で，変数 x を含むコサインの正負冪乗の有限和を求める公式を与えた．この公式から，変数 x を含むタンジェント（もしくはコタンジェント）の 2 乗以上の冪乗の有限和を求める公式を導くことが可能である．しかし，その「追記」（「数学・物理
通信」7-2）で注意したように，各項の r 依存性は，タンジェントの周期 π に対応した (r/N )π ではなく，コサインの周期 2π に対応した (r/N )(2π) であった．
そこで今回の論文では，周期 π に対応した級数の和公式を求めることを考える．タンジェントについて考える代わりに，コタンジェントの和公式を主にして定式化した．そのほうが少しばかりきれいに扱えるからである 2 ．タンジェントに関する公式は，(1.5) に与えるように，コタンジェントに関するそれからすぐに導ける．
出発点は「コサインの正負冪乗の有限和とチェビシェフ多項式」において例外とした 1 乗の場合の和公式である．そしてコタンジェントのよく知られた微分公式 (d/dx) cot x = −(1 + cot2 x) に注目して，一般冪乗の級数和に対する漸化式を導く．この両者を用いて和公式を計算するのである．
以下で考察する級数は，

Wn (N, x) ≡

N
−1 [
∑
r=0

(
rπ )]ncot x +
N

(1.1)

である（r の前に因子 2 がないことに注意）．コタンジェントの周期性 cot(θ + π) = cot θ を使うと，定義から明らかに

Wn (N, x) = Wn (N, x + π)

(1.2)

である．また cot(−θ) = − cot θ と (1.2) を使い，和の取り順を逆にすると，Wn (N, −x) = (−1)n Wn (N, x)(1.3)

を得る．すなわち Wn (N, x) は，n が偶数ならば x の偶関数，奇数ならば奇関数である．
タンジェントの冪乗和との関係は，tan θ = (cot θ)−1 により，(
rπ )]ntan x +
N

(1.4)

(
π)
W−n (N, x) = (−1)n Wn N, x +
2

(1.5)

W−n (N, x) =

N
−1 [
∑
r=0

であるが，cot(θ + (π/2)) = − tan θ により，となるので，タンジェントに関する和公式は改めて計算しなくてもよい．
1 京都大学名誉教授, nbr-nak@trio.plala.or.jp
2 (2.2) と (6.3) とを比較せよ．

2

2

n = 1 の場合
(1.1) の n = 0 の場合は明らかに
W0 (N, x) = N(2.1)

W1 (N, x) = N cot(N x)

(2.2)

である．

(1.1) の n = 1 の場合は，で与えられる 3 ．証明は「コサインの正負冪乗の有限和とチェビシェフ多項式」（「数学・物理通信」6-7）で与えたタンジェントの場合の証明を少し修正すればよい．
証明： 円分方程式 X N − 1 = 0 の N 個の解は X = e2πi(r/N ) (r = 0, 1, · · · , N − 1) であるから，その左辺はXN − 1 =

N
−1
∏

(
)
X − e(2r/N )πi

(2.3)

r=0

のように因数分解ができる．X = e−2xi を代入すると，

e−2N xi − 1 =N−1∏( −2xi)e− e(2r/N )πir=0= (e−xi )N (e(1/N )πi )N (N −1)/2N−1∏( −(x+(r/N )π)i)e
− e(x+(r/N )π)i

(2.4)

r=0

= (−2i)N e−(N x−(N −1)π/2)iとなるから，(rπ )sin x +
N
r=0

N
−1
∏

(
rπ )
1
sin x +
= N −1 sin(N x)
N
2
r=0

N
−1
∏

(2.5)

を得る．ただしサインのオイラー公式を用いた．(2.5) の対数をとりそれを微分すれば，(
2rπ )cot x +
= N cot(N x)
N
r=0

N
−1
∑

(2.6)

が得られる．□

3

コタンジェントの冪乗の級数和
cot θ の微分が − csc2 θ = −(1 + cot2 θ) であることから，
Wn′ (N, x) = −n[Wn−1 (N, x) + Wn+1 (N, x)](3.1)

である．ただしプライムは x に関する微分を表す．n ≧ 1 のとき，(3.1) から漸化式

1
Wn+1 (N, x) = − Wn′ (N, x) − Wn−1 (N, x)n3 この式は森口繁一他著「数学公式 II」p.18 に載っている．なお，(2.5) は p.26 にある．

3

(3.2)

が成立することがわかる．(2.1) と (2.2) より W0 (N, x) = N および W1 (N, x) = w であるから，(3.2) を使え
ば逐次 Wn (N, x) を計算することができる．ただしw ≡ N cot(N x)(3.3)

w′ = −(w2 + N 2 )(3.4)

とおいた．このときである．数学的帰納法により，Wn (N, x) は w と N に関する多項式で，その各項の w に関する次数と N に関する次数との和は n を超えないことがわかる．
以下において，混乱を生ずるおそれがない場合は Wn (N, x) を Wn と略記する．n = 5 までの計算結果は，

W0 = N,W1 = w,
W2 = w2 + N (N − 1),
W3 = w3 + (N 2 − 1)w,
1
4
W4 = w4 + (N 2 − 1)w2 + N (N − 1)(N 2 + N − 3),
3
3
5 2
1
5
3
W5 = w + (N − 1)w + (N 2 − 1)(2N 2 − 3)w
3
3

(3.5)

となる．さらに n = 8 まで進めると，

W6 = w6 + 2(N 2 − 1)w4 +1(N 2 − 1)(17N 2 − 23)w2
15

1
N (N − 1)[2N (N + 1)(N 2 − 4) + 15],15
7
7
W7 = w7 + (N 2 − 1)w5 + (N 2 − 1)(11N 2 − 14)w3
345
1+ (N 2 − 1)(17N 4 − 53N 2 + 45)w,
45
84
W8 = w8 + (N 2 − 1)w6 + (N 2 − 1)(9N 2 − 11)w4
3
15
8+(N 2 − 1)(31N 4 − 88N 2 + 66)w2
315
1+
N (N − 1)[N (N + 1)(17N 4 − 95N 2 + 213) − 315].
315
+

(3.6)

これらの式から Wn の一般形を推測するのは不可能に近い．だがこれをW2 − wW1 = (N − 1)W0 ,W3 − wW2 = (N − 1)W1 ,1
1
W4 − wW3 = (N 2 − 1)W2 + (N − 1)(N − 2)W0 ,
3
3
1 2
1
W5 − wW4 = (N − 1)W3 + (N − 1)(N − 2)W1 ,
3
3

4

(3.7)

および

1 2
1
(N − 1)W4 + (N 2 − 1)(N 2 − 4)W2
345
1
+ (N − 1)(N − 2)(N 2 + 3N − 13)W0 ,45
1
1
W7 − wW6 = (N 2 − 1)W5 + (N 2 − 1)(N 2 − 4)W3
345
1+ (N − 1)(N − 2)(N 2 + 3N − 13)W1 ,45
1
1
W8 − wW7 = (N 2 − 1)W6 + (N 2 − 1)(N 2 − 4)W4
3
45
1
2
2
(N − 1)(N − 4)(2N 2 − 11)W2+
945
1(N − 1)(N − 2)(2N 4 + 6N 3 − 28N 2 − 96N + 251)W0+945
W6 − wW5 =(3.8)のように書き直してみると，幾分規則性が見えてくる．たとえば，W2m − wW2m−1 と W2m+1 − wW2m の展開係数はすべて一致するようである．Wn−2 の係数は，n ≧ 4 で (N 2 − 1)/3 である．さらに Wn−4 の係数は，

n ≧ 6 で (N 2 − 1)(N 2 − 4)/45 であるようだ．このような性質は，Wn−2k に対し，n ≧ 2k + 2 で成立するのかもしれない．しかし，やはり Wn − wWn−1 の一般形の推定は難しい．ただし，N = 1 のときに右辺が 0 にな
ること，および N = 2 のときに右辺が Wn−2 に等しくなることは，次節で示すように一般的に証明できる．

4

N = 2k の場合
N = 1 のときは，定義から明らかに
Wn = cotn x = wn

(4.1)

Wn − wWn−1 = Wn−2

(4.2)

である．

N = 2 のときはが成立する．
証明：

簡単のため，sin x ≡ s，cos x ≡ c と書くと，Wn = cotn x + cotn (x + (π/2)) = cotn x + (−1)n tann x
=

(4.3)

c2n + (−1)n s2n
sn cn

なるゆえ，

c2 − s2 c2n−2 + (−1)n−1 s2n−2·
sc
sn−1 cn−12n
2 2n−2c −s c+ (−1)n (s2n − c2 s2n−2 )=sn cn2(n−2)c+ (−1)n−2 s2(n−2)= Wn −
sn−2 cn−2= Wn − Wn−2

wWn−1 = W1 Wn−1 =となる．これを移項すればよい．□

5

(4.4)

N = 4 のとき，(3.5)，(3.6) から
W2 = w2 + 12,
W3 = w3 + 15w,
W4 = w4 + 20w2 + 68,
W5 = w5 + 25w3 + 145w,(4.5)

W6 = w6 + 30w4 + 249w2 + 396,
W7 = w7 + 35w5 + 378w3 + 1183w,
W8 = w8 + 40w6 + 532w4 + 2512w2 + 2308を得る．また，(3.7)，(3.8) からは

W2 − wW1 = 3W0 ,W3 − wW2 = 3W1 ,W4 − wW3 = 5W2 + 2W0 ,W5 − wW4 = 5W3 + 2W1 ,(4.6)

W6 − wW5 = 5W4 + 4W2 + 2W0 ,W7 − wW6 = 5W5 + 4W3 + 2W1 ,W8 − wW7 = 5W6 + 4W4 + 4W2 + 2W0
を得る．(4.6) における係数は驚くべくシンプルである．

N = 8 のとき，(3.5)，(3.6) から，
W2 = w2 + 56,
W3 = w3 + 63w,
W4 = w4 + 84w2 + 1288,
W5 = w5 + 105w3 + 2625w,(4.7)

W6 = w6 + 126w4 + 4473w2 + 32312,
W7 = w7 + 147w5 + 6762w3 + 92799w,
W8 = w8 + 168w6 + 9492w4 + 194256w2 + 816136
を得る．(3.7)，(3.8) からは

W2 − wW1 = 7W0 ,W3 − wW2 = 7W1 ,W4 − wW3 = 21W2 + 42W0 ,W5 − wW4 = 21W3 + 42W1 ,W6 − wW5 = 21W4 + 84W2 + 70W0 ,W7 − wW6 = 21W5 + 84W3 + 70W1 ,W8 − wW7 = 21W6 + 84W4 + 468W2 + 398W0
6

(4.8)

を得る．
上の結果から，N が 2 の冪乗であると，いずれの式も係数はすべて正の整数になるようである．もし係数がすべて整数だという性質が，N = 2k のときつねに成立するとすると，一般式に現れる係数の N 依存性はかな
り制約のきついものとなる．

(3.7) の第 3 式以下と (3.8) の各式の右辺第 1 項に現れる (N 2 − 1)/3 は，N = 2k のとき必ず整数になる．じっさい，任意の隣接する 3 つの整数はそのうち 1 個が必ず 3 の倍数だから，2k − 1，2k ，2k + 1 のうちの１つは3 の倍数である．2k は 3 の倍数ではないから，N 2 − 1 = (2k − 1)(2k + 1) は 3 の倍数である．
同様に (3.8) の第 2 項に現れる (N 2 − 1)(N 2 − 4) は，N = 2k （k ≧ 2）のとき 45 の倍数である．なぜなら，
上述により 22k − 1 と 4(22(k−1) − 1) はともに 3 の倍数である．したがって，(N 2 − 1)(N 2 − 4) は 9 の倍数である．隣接する 5 つの整数 2k − 2，2k − 1，2k ，2k + 1，2k + 2 のうちの 1 つは 5 の倍数であるが，2k は 5 の倍数ではないから，(N 2 − 1)(N 2 − 4) はまた 5 の倍数でもある．したがってそれは 45 で整除される．

(3.8) 第 1 式と第 2 式の第 3 項に現れる (N − 1)(N − 2)(N 2 + 3N − 13) の場合は，= (N 2 − 1)(N 2 − 4) −15(N − 1)(N − 2) だから上の場合に帰着する．ほかの場合も同様に扱えると期待されるが，一般的な考察は難しい．

5正冪の級数和と負冪の級数和との関係
1 節で指摘したように，タンジェントの冪乗和はコタンジェントの冪乗和に帰着する．すなわち，冪指数 nが負の場合は，(1.5) によって n が正の場合に還元される．しかし，級数の和の性質を利用すると，x の値を変えない還元公式が得られるのである．面白いことに，この還元公式は，N が偶数か奇数かによって形が異なる．

N が偶数ならば，W−n (N, x) = (−1)n Wn (N, x).

(5.1)

証明：

W−n (N, x) =

N
−1 [
∑
r=0

N/2−1
N
−1 )[
(
(
∑
rπ )]n ( ∑rπ )]n=tan x ++
tan x +
N
N
r=0

(5.2)

r=N/2

において，tan θ = tan(θ + π) を使って右辺第 1 項の r についての和を N から 3N/2 − 1 までに変えれば，
3N/2−1 [W−n (N, x) =

∑

r=N/2

(

N
−1 [
(
∑
rπ )]nrπ π )]ntan x +=+
tan x +
N
N
2
r=0

(5.3)

となる．tan(θ + π/2) = − cot θ なるゆえ，

W−n (N, x) = (−1)n

N
−1 [
∑
r=0

(
rπ )]ncot x += (−1)n Wn (N, x)
N

(5.4)

を得る．□

N が奇数ならば，W−n (N, x) = (−1)n [Wn (2N, x) − Wn (N, x)].

7

(5.5)

証明：

W−n (N, x) =

−1)/2
( (N∑
r=0

=

N
−1
∑

+

(N +1)/2−1)/2( (3N∑

=

r=0

=

N
−1 [
∑
r=0

N
−1
∑

+

r=N
N
−1 [
∑

)[

(N +1)/2(rπ )]ntan x +
N)[(rπ )]ntan x +
N((r + (N + 1)/2)π )]ntan x +
N

(5.6)(
(2r + 1)π π )]n+tan x +2N2

= (−1)n

N
−1 [
∑

(
(2r + 1)π )]ncot x +
.
2N

r=0

他方，

Wn (N, x) =

N
−1 [
∑
r=0

(
2rπ )]n
cot x +
2N

(5.7)

であるから，(−1)n W−n (N, x) + Wn (N, x) = Wn (2N, x)(5.8)となり，(5.5) を得る．□

6

タンジェントの冪乗の級数和負冪の場合，すなわちタンジェントの正冪の級数和

W−n (N, x) =

N
−1 [
∑
r=0

(
rπ )]ntan x +
N

(6.1)

を考察する（n > 0）．

−x + (π/2) ≡ x̃ とおくと，(1.3) と (1.5) から，
W−n (N, x) = (−1)n W−n (N, −x) = Wn (N, x̃)(6.2)

である．つまりタンジェントに関する式は，コタンジェントに関する式の x を x̃ に置き換えれば得られる．と
くに，

(
Nπ )
W−1 (N, x) = w̃ ≡ N cot(N x̃) = −N cot N x −2である．すなわち w̃ は，N が偶数ならば −N cot(N x)，N が奇数ならば N tan(N x) である．

(6.3)

(3.5) の x を x̃ に置き換えると，
W−1 = w̃,W−2 = w̃2 + N (N − 1),W−3 = w̃3 + (N 2 − 1)w̃,
14W−4 = w̃4 + (N 2 − 1)w̃2 + N (N − 1)(N 2 + N − 3),
3
3
1
5 2
3
5
W−5 = w̃ + (N − 1)w̃ + (N 2 − 1)(2N 2 − 3)w̃
3
3
8

(6.4)

がわかる．n ≧ 6 も同様．
さて，N が偶数の場合は，w̃ = −w であるから，上の結果は (5.1) と (1.3) から直ちに従うものである．しかし N が奇数の場合は，(5.5) と上の結果とのコンシステンシーは自明ではない．(5.5) の n = 1 の場合は，コタンジェントの倍角公式 cot(2θ) = 21 (cot θ − tan θ) を用いて，

W−1 (N, x) = −[W1 (2N, x) − W1 (N, x)]= −2N cot(2N x) + N cot(N x)= −N [cot(N x) − tan(N x)] + N cot(N x)(6.5)= N tan(N x)= w̃
のように確かめられる．同様に (5.5) の n = 2, 3 の場合は，W−2 (N, x) = W2 (2N, x) − W2 (N, x)= [2N cot(2N x)]2 + 2N (2N − 1) − [N cot(N x)]2 − N (N − 1)
= N 2 [cot(N x) − tan(N x)]2 − N 2 cot2 (N x) + 3N 2 − N(6.6)= [N tan(N x)]2 + N 2 − N= w̃2 + N (N − 1);W−3 (N, x) = −[W3 (2N, x) − W3 (N, x)]= −[2N cot(2N x)]3 − [(2N )2 − 1][2N cot(2N x)]+ [N cot(N x)]3 + (N 2 − 1)N cot(N x)= −N 3 [cot(N x) − tan(N x)]3 − (4N 2 − 1)N [cot(N x) − tan(N x)](6.7)+ N 3 cot3 (N x) + (N 2 − 1)N cot(N x)= [N tan(N x)]3 + (N 2 − 1)[N tan(N x)]= w̃3 + (N 2 − 1)w̃のように確認できる．
この性質はきつい条件であり，Wn (N, x) の形をほとんど決定する．たとえば，W2 (N, x) を w2 + f (N ) と仮定するとき，(6.6) の計算から，f (N ) は関数方程式

f (2N ) = 2f (N ) + 2N 2

(6.8)

を満たさなければならないはずである．(6.8) の解は

f (N ) = N 2 + cN

(6.9)

で与えられる．c は任意定数だが，(4.1) を使えば f (1) = 0 でなければならないので，c = −1 と決まる．
一般に，関係式 Wn (2N, x) = Wn (N, x)+(−1)n W−n (N, x) （N を奇数に限定しなくてもよい）と「W−n (N, x)の w̃ = N tan(N x) での関数形は Wn (N, x) の w = N cot(N x) での関数形と同形である」という命題に，w のn 次多項式 Wn (N, x) の最高次の項の係数が 1 であることと付加条件 Wn (1, x) = wn をつければ，Wn (N, x) の形は一意的に決まることが証明できる（数学的帰納法と関数方程式の線形性を用いる．
）．

9

7

x の特殊値
特殊値 x = 0 のときを考える．cot 0 は無限大になってしまうので，n < 0，N = 2m + 1 の場合を考える．

このとき w̃ = N tan(N x) = 0 となるから，n = −2k ならば定数項のみが生き残る．したがって (3.5) と (3.6)
から，W−2 (2m + 1, 0) = 2m(2m + 1),
2
W−4 (2m + 1, 0) = m(2m + 1)(4m2 + 6m − 1),
3
2m(2m + 1)(32m4 + 80m3 + 40m2 − 20m + 3),
W−6 (2m + 1, 0) =15
2
W−8 (2m + 1, 0) =m(2m + 1)(1088m6 + 3808m5 + 3920m4315+ 280m3 − 868m2 + 322m − 45)(7.1)となる．他方，x = 0 のとき，(1.4) において r = 0 の項はゼロ，r = m + 1 から r = 2m までの項の和は r = 1
から r = m までの項の和に等しいことに注意すれば，

W−2k (2m + 1, 0) = 2S̃k (m)であることがわかる．ただしS̃k (m) ≡

m
∑

tan2k

r=1

(

rπ )
2m + 1

(7.2)

(7.3)

は「三角関数の冪乗の有限和」（「数学．物理通信」5-4）で用いた記号である．S̃k (m) の具体的な式は「三角関数の冪乗の有限和 2」（「数学．物理通信」5-6）の (2.5)（k = 1, 2, 3 の場合）と付記の第 2 式（k = 4 の場合）に与えた．これらを 2 倍したものはもちろん (7.1) の右辺と完全に一致する．
同様に，特殊値 x = π/(2N ) を考えると，cot(π/2) = 0 なので w = 0 となる．つまり，Wn の定数項のみが生き残る．
まず N が偶数ならば，N = 2m とすると，(
π )W2k 2m,= W2k (2m, x)= 2Uk (m)4mw=0であることがわかる．ただしUk (m) ≡

m−1
∑

cot2k

r=0

( (2r + 1)π )
4m

(7.4)

(7.5)

とする．その具体的な式は，上記「三角関数の冪乗の有限和 2」の (3.6)（k = 1, 2, 3 の場合）と付記の第 3 式（k = 4 の場合）に与えた．
奇数 N = 2m + 1 の場合，W2k の定数項は W−2k のそれと同じだから，k = 1, 2, 3, 4 では (7.1) の右辺と一10致する．じっさい，そのことは直接の変形(
W2k 2m + 1,2m
) ∑( (2r + 1)π )
π
=
cot2k
2(2m + 1)
2(2m + 1)
r=0

=

2m
∑

tan2k

( (2r + 1)π2(2m + 1)

r=0

=

2m
∑

tan2k

r=0
m
∑

=2

r=1

−

π)
2

( (r − m)π )tan2k2m + 1(rπ )2m + 1= 2S̃k (m) = W−2k (2m + 1, 0)
でも確認できる．

11

(7.6)

量子力学における周期ポテンシャル問題 (5)世戸 憲治 ∗A Periodic Potential Problem in Quantum Mechanics (5)Kenji Seto∗

1 はじめに
前回の「量子力学における周期ポテンシャル問題 (4)」(「数学・物理通信」7 巻 1 号) に続いて，ここでは，2次式で表される周期ポテンシャルが負の場合を解析してみる．すなわち，1 次元の Schrödinger 方程式

[

−

]
ℏ2 d2
+ V (x) Ψ = EΨ
2
2m dx

(1.1)

において，ポテンシャル V (x) が周期 2ℓ の周期関数 V (x + 2ℓ) = V (x) で，−ℓ ≤ x < ℓ の範囲で

V (x) = V0

( x )2ℓ,

V0 < 0(1.2)と定義された場合，前回はこの V0 が正の場合を扱ったが，今回は，これが負の場合を解析してみる．途中まで
は，前回のものと同じ処方箋に従うが，この方程式の実関数で表される解を求めることは，かなり煩雑な計算を余儀無くされる．
以下，数式簡素化のため，座標 x, および，ポテンシャルの大きさ V0 , エネルギー E を無次元化し，

x
→ x,
ℓ

2mℓ2
V0 → V0 ,
ℏ2

2mℓ2
E → E
ℏ2

(1.3)

と改めて置き直すことにする．この置き換えで方程式は，

]
+
E
−
V
(x)
Ψ =0
2

[ d2
dx

(1.4)

となり，ポテンシャル V (x) は周期 2 の周期関数 V (x + 2) = V (x) で，−1 ≤ x < 1 の範囲でとなり，この範囲での方程式は

V (x) = V0 x2

(1.5)

]
2
Ψ =0
+
E
−
V
x
0
dx2

(1.6)

[ d2

となる．
∗ 北海学園大学名誉教授

E-mail: seto@pony.ocn.ne.jp12

2 方程式の解法
初めは，この方程式の解法として，前回のものをそのまま踏襲した形で書いておく．すなわち，定数 µ をµ = (4V0 )1/4

(2.1)

と導入し，これを用いて，変数 x から z に，また，エネルギー E から κ に，z = µx,

E
1
−
2
µ
2

(2.2)

1 z2 ]
−
Ψ =0
2
4

(2.3)

κ=

と変換すると，−1 ≤ x < 1 での方程式 (1.6) は，

[ d2
dz

+κ+
2

となる．これは，Weber の微分方程式と呼ばれるもので，その解は，放物柱関数 (Weber 関数) Dκ (z) であり，
これは合流型超幾何関数 1 F1 を用いて，
κ/2Dκ (z) = 2√ −z2 /4 [πe√
( κ 1 z2 )( 1 − κ 3 z 2 )]
1
2z
;
−
, ;
1 F1 − ,1 F1
Γ ((1 − κ)/2)2 2 2
Γ (−κ/2)
2
2 2

(2.4)

と定義される *1 ．ここまでは，前回のものと同じであるが，今回の V0 < 0 の場合は，(2.1) 式で定義された µは，√
ν ≡ 4 4|V0 |µ = eπi/4 ν,(2.5)と複素数になる *2 ．ここで，ν はこの第 2 式で定義する．これから，(2.2) 式で定義された z ，および，κ は，

1
κ = −iλ − ,2z = eπi/4 νx,

λ≡

E
ν2

(2.6)

となる．ただし，λ はこの第 3 式で定義する．以下，関数定義の利便性を考慮して，z に替わる実変数として，

ζ = νx(2.7)

を定義しておく．
ここで，Weber 関数 (2.4) 式について，元の方程式 (2.3) が z の符号反転に対し不変，したがって，ζ の符号反転にも不変であることから，この大括弧をはずしたときの第 1 項目，第 2 項目の定数係数を除いたものを，

1 1 i )
+ , ; ζ22
4 2 2
( iλ 3 3 i )−iζ 2 /4
S2 (ζ) =ζ e
+ , ; ζ2
1 F1
2
4 2 2
S1 (ζ) =e−iζ /4 1 F1
2

( iλ

(2.8)

と定義したとき，これらも，また，元の方程式 (2.3) の解となることに注意する．ここで，S1 (ζ), S2 (ζ) は，それぞれ，偶関数，奇関数である．さらに，元の方程式 (2.3) は実係数の方程式なので，これらの式の複素共役を*1 「数学公式 3」(岩波全書) p.75-78*2 より一般には，µ = e±πi/4 ν,or e±3πi/4 ν も候補となり得るが，ここでは，1 つの候補を求めるだけで十分である．

13

とったものも解になるので，これらの実部，あるいは，虚部をとったものも解である．ここでは，実部をとることにし，[( iλ 1 1 i )]2T1 (ζ) =ℜ e−iζ /4 1 F1
+ , ; ζ22
4 2 2[
( iλ 3 3 i )]−iζ 2 /4
+ , ; ζ2
T2 (ζ) =ℜ ζ e
1 F1
2
4 2 2

(2.9)

を方程式 (1.6) の −1 ≤ x < 1 における 2 個の独立解とする *3 ．なお，S1 , S2 と同じく，T1 (ζ), T2 (ζ) も，それぞれ，ζ の偶関数，奇関数となる．これらの解を，合流型超幾何関数の部分で，ζ の低次の項に展開すると，T1 (ζ) = cos
T2 (ζ) = cos( ζ 2 )(
4
( ζ 2 )(
41−

( ζ2 ) ζ2
λ 2)
ζ + sin
+ ···
2
4 4

( ζ2 ) ζ3
λ )ζ − ζ 3 + sin
+ ···
6
4 4

(2.10)となる．ここで，これら関数の Wronskian

W (ζ) = T1 (ζ)T2′ (ζ) − T1′ (ζ)T2 (ζ)(2.11)

を求める．ここで，プライムは ζ での微分を表す．T1 , T2 が方程式 (1.6) を満たすことから，dW/dζ = 0 となるので，この Wronskian は定数である．そこで，(2.10) 式を利用して，ζ = 0 でこの値を見積もると，W (ζ) ≡ 1

(2.12)となる．
この関数 T1 , T2 を用いて，方程式 (1.6) の一般解を，任意定数 A, B を用いて，

Ψ (x) = AT1 (νx) + BT2 (νx)(2.13)

とおく．これは，あくまでも，範囲 −1 ≤ x < 1 での解であって，ポテンシャルが周期的であっても，この解が周期的に続くわけではない．Floquet の定理によると，波動関数は，周期が 1 つ増えるごとに元の波動関数に，伝播数と呼ばれる定数 K (0 ≤ K ≤ π) を用いた位相 eiK あるいは e−iK を付けたものが次の周期の波動関数になることが知られている．したがって，範囲 1 ≤ x < 3 での波動関数は

[
(
)
(
)]
Ψ (x) = eiK AT1 ν(x − 2) + BT2 ν(x − 2) ,

or

[
(
)
(
)]
Ψ (x) = e−iK AT1 ν(x − 2) + BT2 ν(x − 2) (2.14)となる．この K の前に付く符号は縮退した 2 個の波動関数が存在することを意味するが，ここでは当分のあいだこの第 1 式の方で議論を進める．また，K = 0, π のときは，これら 2 つの波動関数は同じものになってしまうが，そのときは K で微分したものがもう 1 つの解を与えるはずである．ただし，この種の計算は面倒になるので，ここではその場合を例外として扱わない．したがって，以後，K の範囲は 0 < K < π とする．
つぎに成すべきことは，この (2.13) (2.14) 式を，x = 1 で，関数自身とその微係数が連続となるように繋ぐことである．結果は，[]AT1 (ν) + BT2 (ν) = eiK AT1 (ν) − BT2 (ν) ,[]AT1′ (ν) + BT2′ (ν) = eiK − AT1′ (ν) + BT2′ (ν)(2.15)*3 仮に，T1 , T2 のいずれか一方でも虚部を取ってしまうと，以下で述べる Wronskian がゼロとなるので，2 個の独立解としては不適格である．また，「付録 1」に解そのものとしても不適格な理由を示す．

14

となる．ここで，関数 T1 , T2 の偶奇性を用いた．あるいは，これらをまとめて，((1 − eiK )T1 (ν))( )(1 + eiK )T2 (ν)
A(1 + eiK )T1′ (ν)(1 − eiK )T2′ (ν)

B

( )
0
=

(2.16)

0となる．ここで，A, B が両方ともゼロとならないためには，この係数行列式の値がゼロでなければならず，

cos(K) = T1 (ν)T2′ (ν) + T1′ (ν)T2 (ν)(2.17)を得る．ここで，Wronskian の式 (2.11) (2.12) を用いた．この式から，伝播数 K が求まるためには，右辺の絶対値が 1 以下でなければならず，これから，関数 T1 , T2 に含まれる λ，あるいは，(2.6) 式からエネルギー

E の取り得る値の範囲が決まり，いわゆるエネルギーのバンド構造が決まることになる．この意味で，この式を固有値方程式と呼ぶことにする．
ここで，A, B の値を (2.15) の第 1 式が成り立つように，

B = −(1 − eiK )T1 (ν)

A = (1 + eiK )T2 (ν),(2.18)と選ぶことにする．これを用いて，一般の x に対する波動関数は，n を整数として，2n − 1 ≤ x < 2n + 1 の範
囲で，

[
(
)
(
)]
Ψ (x) = eiKn (1 + eiK )T2 (ν)T1 ν(x − 2n) − (1 − eiK )T1 (ν)T2 ν(x − 2n)(2.19)と表される．もちろん，これは，まだ，規格化されたものではない．

3 波動関数の規格化
以下では，エネルギー依存性を明示するため，波動関数 Ψ (x)，および，関数 T1 (ζ), T2 (ζ) を，それぞれ，

Ψ (x, E), T1 (ζ, λ), T2 (ζ, λ) と記すことにする．ここで，波動関数の規格化のため，一般に異なる 2 個のエネル
ギー E, E ′ を考え，これに属する (2.6) 式の λ を λ, λ′ ，また，これに対応して， (2.17) 式から決まる伝播数

K を K, K ′ とする．波動関数の規格化の式は，(2.19) 式を考慮すると，∫ ∞∞ ∫ 2n+1∑′

Ψ (x, E)Ψ (x, E )dx =
−∞

n=−∞

′

Ψ (x, E)Ψ (x, E )dx =

∞
( ∑

2n−1

−i(K−K ′ )n

)∫ 1

e

n=−∞

Ψ (x, E)Ψ (x, E ′ )dx

−1

(3.1)
となる．ここで，上付傍線は複素共役をとることを意味する．この和の部分は，超関数公式

lim

M →∞

M
∑

′

e−i(K−K )n = 2πδ(K − K ′ ),0 < K, K ′ < π

(3.2)

n=−M

が使え，デルタ関数を用いた直交性がでる．つぎに，積分の部分を実行するには，E, E ′ に対応する (1.6) 式を，

[ d2
dx

]
2
+
E
−
V
x
Ψ (x, E) = 0,
0
2

]
′
2
Ψ (x, E ′ ) = 0
+
E
−
V
x
0
dx2

[ d2

(3.3)

と書いておく．ただし，E に対応する方は複素共役を取ったものにしておく．この第 1 式に Ψ (x, E ′ ) を，また，第 2 式には Ψ (x, E) を掛けてから，辺々を引き算すると，]d [dΨ (x, E ′ ) dΨ (x, E)
Ψ (x, E)−Ψ (x, E ′ ) − (E − E ′ )Ψ (x, E)Ψ (x, E ′ ) = 0
dx
dx
dx
15

(3.4)

となり，これを，−1 から 1 まで積分すると，∫ 1

Ψ (x, E)Ψ (x, E ′ )dx =

−1

]1
1 [
dΨ (x, E ′ ) dΨ (x, E)
′
Ψ
(x,
E)
−
Ψ
(x,
E
)
E − E′
dx
dx
−1

(3.5)

となる．ここで，この右辺に n = 0 とした (2.19) 式を代入すると，16 項がでてくるが，関数 T1 , T2 の偶奇性を用いると半分は消え，さらに，Wronskian の (2.11) (2.12) 式，固有値方程式 (2.17) 式を用いると，うまくま
とめることができ，∫ 1

Ψ (x, E)Ψ (x, E ′ )dx =

−1

[
′
′ ]
2ν [
T1 (ν, λ)T2 (ν, λ) (1 + e−i(K−K ) ) cos(K ′ ) − (e−iK + eiK )
′
E−E
]
[
′
′ ]
− T1 (ν, λ′ )T2 (ν, λ′ ) (1 + e−i(K−K ) ) cos(K) − (e−iK + eiK )(3.6)となる．これを (3.1) 式に代入し，分母の E, E ′ を (2.6) の第 3 式を用いて，λ, λ′ で表すと，

∫ ∞

[
[
′
′ ]
4π
T1 (ν, λ)T2 (ν, λ) (1 + e−i(K−K ) ) cos(K ′ ) − (e−iK + eiK )′ν(λ − λ )][′′ ]− T1 (ν, λ′ )T2 (ν, λ′ ) (1 + e−i(K−K ) ) cos(K) − (e−iK + eiK ) δ(K − K ′ ) (3.7)Ψ (x, E)Ψ (x, E ′ )dx =−∞という結果になる．ここで，1 つ注意が必要である．λ を決めると固有値方程式 (2.17) から伝播数 K が一意に決まるが，逆に K を決めたときには，λ は一意に決まるわけではなく，各バンドごとの λ が決まるということである．したがって，λ と λ′ が異なるバンドに属するときは，例え，K = K ′ であっても λ ̸= λ′ なので，この式の値はゼロとなり，異なるバンド間の直交性がでる．また，λ と λ′ が同じバンドに属するときは，この式は0/0 の不定形となるので，極限値をとることにする．すなわち，l’Hôpital の定理にしたがって，分母，分子をそれぞれ λ で偏微分してから，λ′ → λ とおく．結果は，∫ ∞−∞Ψ (x, E)Ψ (x, E ′ )dx =8π
dKT1 (ν, λ)T2 (ν, λ) sin(K)δ(K − K ′ )
ν
dλ

(3.8)

となる．この式のデルタ関数を除いた部分は正定値になるべき量であり，0 < K < π で sin(K) は正なので，積T1 (ν, λ)T2 (ν, λ) と dK/dλ は必ず同符号でなければならない．このことの直接証明はここではしないが，次節の数値計算例では確かに満たされている．さらに，この式のデルタ関数の部分をエネルギー E で表すようにす
ると，

∫ ∞
−∞

Ψ (x, E)Ψ (x, E ′ )dx = N 2 (E)δ(E − E ′ ),N 2 (E) = 8πν T1 (ν, λ)T2 (ν, λ) sin(K)(3.9)となる．ここに，N (E) は規格化定数で，Ψ (x, E)/N (E) が規格化された波動関数となる．なお，K の符号を変えた (2.14) の第 2 式を採用した場合は，sin(K), dK/dλ の符号が両方共に反転するので，ここでの規格化の式はそのまま変更なしに使えることを注意する．

4 数値計算例(2.9) 式で定義された関数 T1 (ζ, λ), T2 (ζ, λ) の λ = 0 の場合については，Bessel 関数で表すことができることを「付録 1」に示すが，一般の λ についてはどのような関数なのかは，この式からは想像できない．そこで，
数値的にグラフ化したものを以下の図 1-1, および，図 1-2 に示す．この図は，水平右向きに ζ 軸，右斜め上方

16

向に λ 軸，上方向に T1 (図 1-1)，または，T2 (図 1-2) をとって，0 ≤ ζ ≤ 5, 0 ≤ λ ≤ 10 の範囲で立体的に図示したものである．前回で用いた関数，S1 (z, κ), S2 (z, κ) に比べると，今回のものはかなりおとなしい関数で，T1 , T2 共に ζ 方向にゆっくり振動しているのがわかる．ただし，今回のものは，E < 0 したがって，λ < 0 の領域もあり得るが，この領域では T1 , T2 共に振幅はかなり大きくなるが，図の見易さからこの領域を図示して
いない．

図 1-1

関数 T1 (ζ, λ) のグラフ

図 1-2

関数 T2 (ζ, λ) のグラフ

つぎに，固有値方程式 (2.17) を用いて，ポテンシャルの大きさ V0 を一定値に固定したとき，エネルギー E の値に対し，伝播数 K がどのように決まるかを数値的に求めてみた．V0 = −50 とした場合を，−50 < E < 250の範囲で，以下の図 2 に示す．エネルギーが小さいときは K の値が一定の幅を持ちながら飛び飛びに決まり，
いわゆるエネルギーのバンド構造が現れる．エネルギーが大きくなるにつれ，これらのバンドは癒着してしまい
完全なバンドではなくなるが，これらも含めて，バンドということにする．この図で，曲線に付いている色は，(3.8) 式のところで述べたように，積 T1 (ν, λ)T2 (ν, λ) が正のときは赤，負のときは青となるように付けたものである．この T1 (ν, λ)T2 (ν, λ) の正負は，これら曲線の傾き dK/dE = dK/ν 2 dλ の正負と一致していることが
わかる．

図2

K と E の関係 (V0 = −50)

図 3 は，V0 の値を固定することなく，|V0 | - E 平面を細かなメッシュに分割し，各メッシュ間をスキャンさせながら固有値方程式 (2.17) の右辺の値を求めていき，その絶対値が 1 以上のときは何も印さず， 1 のときは赤17

点で，また，1 より小さくなるにつれ色を連続的に変化させ， −1 のときに青点になるように印したものである．
この図から，全体のバンド構造がどのようになっているかが良く理解できる．

図3|V0 | - E 平面におけるバンド構造

5 おわりに手を替え品を替え，ではなく，手も替えずに，品だけちょっぴり替えた同じレシピで，周期ポテンシャルの問
題を 5 編も書いてしまった．しかも，今回のものは，前回のものと解析的にはほとんど同じである．しかし，前回のものは，実数の範囲での計算だったのに対し，今回のものは，複素数の数値計算が入ってくるので，これにかなり手間取ってしまった．この数値計算は，Visual Basic を使ってやっているが，これには，複素演算のルーチンが備わっていないため，自分でプログラムを組むしか方法がなかったためである．
話は変わるが，最初に，この周期ポテンシャル問題に興味を持ったのは，バンド間の直交性とバンド内におけるエネルギーの直交性がどのように導かれるか，ということが，漠然としていてよく理解できなかったことである．本によっては，i, j をバンド番号とし，その直交性 δi,j と，バンド内におけるエネルギーの直交性

δ(E − E ′ ) の積，δi,j δ(E − E ′ ) が出るように書かれたものがある．しかし，1 本の式から，この Kronecker とDirac の 2 個のデルタ関数が出ること自体，おかしなことである．実際はここで示したように，エネルギーの直交性だけ導出できれば，それで十分であり，これは当然のことながら，異なるバンド間の直交性も意味している．しかも，バンドといっても，完全に離れたものばかりではなく，癒着してしまった不完全バンドもあり得る
ので，バンド間の直交性ということ自体，無意味ではないかと考えられる．これまでに書いた 5 編の論文でこのことを明らかにしてきたつもりである．

6 付録 1: エネルギーがゼロの場合の関数 T1 , T2
「数学公式 3」(岩波全書) P.57 に載っている合流型超幾何関数と Bessel 関数を結び付ける公式

Jν (z) =()(z/2)ν e−iz
1
1 F1 ν + , 2ν + 1; 2izΓ (ν + 1)
2

18

(6.1)

において，ν = −1/4, または，ν = 1/4 とおき，z = ζ 2 /4 とすると，

e−iζ /4 1 F1
2

( ζ2 )(1 1 i )( 3 )( ζ 2 )1/4
J−1/4
, ; ζ2 = Γ
4 2 2
4
8
4

( ζ2 )(3 3 i )( 5 )( ζ 2 )−1/4
J1/4e−iζ /4 1 F1 , ; ζ 2 = Γ
4 2 2
4
8
4

(6.2)

2

という式が得られる．これは，(2.9) 式で定義された関数 T1 (ζ, λ), T2 (ζ, λ) で，λ = 0, したがって，エネルギーを E = 0 としたときの関数が，Bessel 関数を用いて，T1 (ζ, 0) = Γ( 3 )( ζ 2 )1/4
4

8

J−1/4

( ζ2 )
4( 5 ) ( ζ 2 )−1/4( ζ2 )
T2 (ζ, 0) = Γ
ζ
J1/4
4
8
4

(6.3)

と表されることになる．しかも，(6.2) 式の左辺には，虚数の i が含まれるが，右辺は虚数が消えて実数のみとなっている．したがって，(2.9) 式のところで，仮に，虚部を採用したとすると，E = 0 でこれら関数は恒等的にゼロとなり，Wronskian がゼロとなるばかりでなく，解そのものとしても不適格なものになってしまう．

[ 謝辞 ]
今回もこの原稿を書くにあたって，京都大学名誉教授の中西襄先生にご精読いただき，たくさんのコメントをいただきました．ここに，謹んで感謝いたします．

19

超弦理論と Dirac のハサミ
Superstring Theory and Dirac’s Scissors
伊藤仁之
Hitoshi Ito1

はじめに⌈ 数学・物理通信 ⌋ のあまり熱心でない，読者です．近畿大学を退職して，木版画の制作に没入するようになり，時間のゆとりがなくなったからだが，数式を追いかけるのが煩わしくなってきたと言うこともある．この度 (７巻１号)，数式がひとつしかない中西先生の ⌈ 超弦理論はなぜつぶれたのか ⌋ を見つけ，これなら読めるかもしれないと，読み始めたが，こんどは，数学ばかりか，物理的概念もひとつの単語か，数行の説明で登場し，未知のものもたくさんあるし，以前勉強したはずだがなかなか思い出せないというのも多く，正確に理解するのは不可能だった．しかし，超弦理論展開の過程で数々の論理的不整合があるにもかかわらず，希望的観測でそれらの無視を続け，迷路に迷い込んでいったことを先生が指摘していることを読み取ることができたのは幸いだった．というのは，私が超弦理論に関心を持ち始めたのは今年になってからだが，当初からこの理論に疑いを抱いていたからだ．
超弦の超が怪しいと言うことである．中西先生の解説にあるように，超はボゾンとフェルミオンが対になって現れることだが，これを弦理論に当てはめると，弦はスピノルとしても振る舞うと言うことだ．ボゾンは古典力学とのつながりが強い Schrödinger 方程式を相対論化した場の方程式を量子化すると出現するから，点に変わって１次元的に広がった弦を素子とした場を考えるというアイデアもあり得るかもしれないが，フェルミ
オンはそういうわけにはいかない．フェルミオンは Dirac 方程式に従う場を量子化すると出現する．この方程式が内包する排他律が反交換関係による量子化を要求するわけだ．
弦にはスピノル性が潜在している，と言ったのは他ならぬ Dirac だった．Dirac はハサミに取り付けたひも(弦) の他端をどこかに固定し，ハサミを２回転して生じるひものねじれが並進運動で解けることに気づき，弦のスピノル性を示唆した (Dirac 自身が書いたものは残ってないようで，どこまで本気だったか不明)．弦にス
ピノル性が潜在しているなら，弦にフェルミオン的ラグランジアンを乗せて量子化すれば，スピノル性が顕在化してフェルミオンが現れ，超弦理論万々歳となったかもしれないが，弦にはスピノル性などないことを次節で示す．(なお，中西調書では”ボゾン弦理論も完全に棄却される”となっている．)

2

Dirac のハサミ
ひもの一端をどこかに固定し，他端にハサミを取り付ける．ひもを軸として右回りに２回転する．これを斜め上方から見ており，視線の方向を空間固定の回転軸とする．この２回転ねじれは，ひもの中程を持ち上げてハサミがひもの下を通るようにして，並進運動で元の位置に戻すと，ほどける．これはスピノルの (巨視的！)模型だ．というわけだが，これは怪しい．実験幾何学の手法で謎を解いてみよう．Dirac がいう並進運動は回転運動に分解することができる．２回転状態から，ひもの上で左１回転する．この回転を元に戻すとき，ひもを持ち上げ，ひも軸を反転させて右１回転する．この右回転は，ひも軸を基準にすると左回転になる．つまり，
ひも軸を基準にすると，右２回転は左２回転で元に戻るという当たり前のことを言ってるに過ぎない [1]．
この Dirac のハサミは，さすがに，専門誌の論文が引用することはないようだが，専門家向けの本が紹介し，スピノル性のさらなる展開が可能であるかのように扱っている [2, 3]．超弦理論が着想され，大発展する２０
世紀後半は，⌈ 弦にはスピノル性が潜在している ⌋ という思想が蔓延していた時代ではなかったかと推測する

20

が，科学啓蒙誌をあさったくらいでは，確かめようがない．そこで，それを裏付ける雰囲気のようなものに触れることにする．

3超弦理論の現状管見
超弦理論の現状を見極めるには，専門家と批判者を集めた国際会議を開くのが最良の策であろうが，私のようなアマチュアはとうてい参加できない．横目で眺めていようと思っていたら，⌈ 超ひも理論でひらく時空のトビラ ⌋ という講演会が目に飛び込んできた．大阪大学で開催された物理学会の一般聴衆に向けた啓蒙講演会だった．聴講してみた．講演者は漫画チックであることを標榜する面白い先生だったが，超弦の超の説明をスルーしていた．そこで，⌈ Dirac のハサミを枕に，超の意味を説明すれば，よりいっそう漫画チックになって面
白いではないか ⌋ と質問してみたが，取り合ってもらえなかった．講演後，2,3 言葉を交わしたが，”The only

game in town”の人かどうか見極めかねた．
このすこし前に，超弦理論を大まかにでも把握しておこうと，大学院向けの参考書を購入してみた [4]．まえがきの冒頭に，ここで議論する超弦理論の描像を一口に言ってしまえば，⌈ かって超弦理論と呼ばれていた ⌋広がった物体の理論というものに近い と書いてあった．

参考文献[1] 伊藤仁之， Dirac のハサミと Feynman のワインダンス 物理教育 49 巻 (4),2001 訂正記事投稿中[2] C.W. Misner, K.S. Thorne and J.A. Wheeler, GRAVITATION, Chapter 41, Freeman and company, NY, 1973.
[3] R.Penrose and W.Rindler, Spinors and Space-Time Vol.1, pp 42-44, Cambridge U.P., 1984
[4] 太田信義， 超弦理論・ブレイン・M 理論, 丸善出版,201221

n 次元の球の体積 1
矢野忠1

The Volume of n-dimensinal Hypersphere 1
Tadashi YANO2

1

はじめに
すでに 3 回のシリーズとして球の体積を求める一連のエッセイを愛数協（愛媛県数学教育協議会）の機関誌『研究と実践』に掲載した [1][2][3]．その後，いくつかの文献を見る機会があった [4][5][6][7][8]．そのことから前のエッセイを改訂するのが適当と思われた．これが [1] を改訂して「数学・物理通信」に投稿する第一の理由である．第二の理由として，これがより重要な理由だが，以前の記事の読者は小学校，中学校や高校の先生に限られていたが，「数学・物理通信」に掲載すれば，もっと広い読者層が期待できる 3 ．

[1] の改訂にあたって，文章を大幅に書き換えることも考えたが，それを断念した．その代わりにその後に読んだ文献から生じた私の疑問等を示すためにその該当箇所の近くに（問）をつけた．その（問）に対する解答
例は 6 節に述べた．はじめて読むときはこの（問）の箇所は無視して読んでほしい．
（問）の部分はかなり細かな事項になるからである．もっとも問に対する答で示された知識はこのシリーズの後のエッセイで詳しく述べられることも多い．
さて，まえおきはそのくらいにして，何回かに分けて n 次元の球の体積についてその求め方を述べてみたい．
数学には縁のない普通の人にとって n 次元の球について考える必要は一生に一度もないだろう．ましてや，その体積を求めようなんて気になることはないかもしれない．
3
まあ，でも半径 a の円の面積が πa2 であることは小学校で学んだであろうし，半径 a の球の体積が 4π
3 a であることは高等学校で学ぶだろう．それらをどうやって求めたかということは知らないにしても．
また，円は 2 次元の球だし，3 次元の球が普通に球といわれているものであることを知れば，4 次元の球とか 5 次元の球というのはどういう風に数学的に表されるのだろうかと少しは関心が起こってくるかもしれない．
また，1 次元の球だって考えられるのだ．

2球の数学的な表現
半径 a の n 次元の球を数式で表すことはそんなに難しいことではない．いま，n 次元の直交座標系を考えよう．変数として n 個の変数を考える必要があるから，それらを (x1 , x2 , · · · , xn ) としよう．
なんだか，はじめから難しそうだが，これは普通の平面上の直交座標なら，(x1 , x2 ) となる．添字を使わない表し方だと (x, y) と表される．すなわち，x1 = x, x2 = y である．変数の数が多いときには使う文字の数が足らなくなるので，数学では変数に添字をつけて表すのが普通である．
さて，本筋にもどって，n 次元の球は

x21 + x22 + · · · + x2n = a2 ,

a>0

(2.1)

1 元愛媛大学工学部
2 yanotad@earth.ocn.ne.jp
3 n 次元の球の体積の求め方を調べてみたいと思った動機はカオスの話題と関係している．2 次の写像 xn+1 = axn (1 − xn ) で，コントロール・パラメータ a はその値が a = 3.5699456 · · · をこえるとカオス的となる [9]．そのことを山口昌哉先生の愛媛大学での集中講義
で聞いて数年してから，このカオス的な写像を乱数として使うことができないかと考えた（同じ着想の研究として [10][11] を参照）．この 2 次の写像で発生させた数が乱数として使えるのかどうかに関心があったので，卒論ゼミに来た学生にそれをテーマに与えた．それはうまくいかなかったのだが，それを試す一つの手段としてカオス的に発生した数でモンテカルロ法で多重積分を評価してみたらどうかということを考えた．答えのわかっている多重積分の一つとして n 次元の球の体積に関心をもった [12]22と表される．3 次元空間での球と区別するために n 次元の「超球 (hypersphere)」というのが正しいのかもしれないが，今後簡単のために「超」という字を省いて単に球という．a は球の半径である．
低次元の場合を考えてみよう．いま，n = 1, 2, 3 のとき

x21 = a2

(2.2)

x21 + x22 = a2

(2.3)

x21 + x22 + x23 = a2

(2.4)

である．もちろん (2.3) は平面上の円の方程式であり，(2.4) は 3 次元の球の方程式である．

3低次元の球の体積
低次元の球を表す式を (2.2)-(2.4) で与えたので，これらの球の体積を求めてみよう．まず，1 次元の球の体

積 ℓ は −a ≤ x1 ≤ a をみたす線分の長さで

∫ a
ℓ=
−a

で与えられる．つぎに 2 次元の球の 体積

dx1 = 2a

S とは実は円の面積のことであるから，
S = πa2で与えられる．これは 2 重積分を用いて

(3.1)

(3.2)

∫∫

S=

D : x21 + x22 ≤ a2

dx1 d2 ,(3.3)

Dで求められる．直交座標 (x1 , x2 ) は 2 次元の極座標 (r, θ) に

x1 = r cos θ
x2 = r sin θで変換すれば，微小面積要素 dS = dx1 dx2 は 2 次元の極座標では

dx1 dx2 = rdrdθ
となるから，

∫∫
S=

rdrdθ
D
∫ a

∫ 2π
rdr
dθ
0
0
(
)
1 2
= 2π
a
2
=

= πa2

(3.4)

確かに円の面積 S = πa2 が得られる 4 ．

3 次元の球の体積 V は
V =

4π 3
a
3

4 2 次元以上の直交座標系から極座標系への変換ではヤコビアンの計算が必要になる．これは [13] で求められる．

23

(3.5)

で与えられるのは周知のことだが，3 重積分を実際に計算して，(3.5) で与えられた 3 次元の球の体積も導出してみよう．
そのためにまず 3 次元の極座標 (r, θ, ϕ) を次式で導入しよう．

x1 = r cos θ
x2 = r sin θ cos ϕx3 = r sin θ sin ϕ
3 次元の球の体積 V は 3 重積分を用いて∫∫∫V =dx1 dx2 dx3 ,D : x21 + x22 + x23 ≤ a2

(3.6)

D

で定義される．(3.6) の 3 重積分の微小体積要素 dV = dx1 dx2 dx3 は 3 次元の極座標では

dx1 dx2 dx3 = r2 sin θdrdθdϕであるから，∫∫∫
r2 sin θdrdθdϕ

V =
∫ a

D

∫ π
∫ 2π
r2 dr
sin θdθ
dϕ
0
0
0
∫ a
= 4π
r2 dr

=

0

4π 3
a
=
3

(3.7)

ここで，上から 2 行目の dΩ ≡ sin θdθdϕ は微小立体角であるので全立体角 Ω = 4π であることを用いてもよい．

4n 次元の球の体積

3 次元の球の体積を求めた，上の方法を一般化すれば，n 次元の球の体積を求めることは原理的に可能であるが，その計算はかなり厄介である．それで，ここでは便法を考えよう [14]5 ．
そのためにつぎの事実に着目しよう．1 次元の球の体積 ℓ は球の半径 a に比例し，2 次元の球の体積 S は a2に比例し，また 3 次元の球の体積 V は a3 に比例している．したがって，n 次元の球の体積 Vn は an に比例することがわかる．この推論は n 次元の球の体積は長さ L の n 乗であるから次元的にもあっている．
そこで，n 次元の球の体積 Vn を

Vn = Cn an(4.1)と表すことにしよう．ここで，Cn は定数である．

——————————————————————————（問 4.1） [5][6] では n 次元の球の体積 Vn とその球の表面積 Sn を考えて，それらの関係を導いている．この関係を導け．
——————————————————————————
これから定数 Cn を求めよう．そうすれば，n 次元の球の体積 Vn が求められる．そのためにまずつぎの積分を

5 [5][6] はこの [14] にヒントを得ているのではないかと思うが，[5] は文献を引用していないので，自分で考えたのかもしれない．[6]も文献の引用はないが，書かれた文からは明らかに自分の考案した計算法ではないことが察せられる．

24

考えよう．

∫ ∞
In =

···

−∞
(∫ ∞

=

∫ ∞

exp[−(x21 + x22 + · · · + x2n )]dx1 dx2 · · · dxn
)n
−x2
e
dx
−∞

−∞

=π

n/2

(4.2)

ここで

∫ ∞

e−x dx =
2

√
π

−∞
6

であることを用いた ．

——————————————————————————（問 4.2） なぜ唐突に (4.2) の積分を考えたのか．
——————————————————————————
一方，x21 + x22 + · · · + x2n = r2 とおき，上の積分 In を半径 r から r + dr の球殻に分けて積分することを考えれば 7 ，
dx1 dx2 · · · dxn = dVn= Cn drn= Cn nrn−1 dr= Cn [(r + dr)n − rn ]であるから，求める積分 In は

∫ ∞rn−1 e−r dr2

In = Cn n

(4.3)

(4.4)

0

となる．

——————————————————————————（問 4.3） n 次元の極座標を定義して，Vn = Cn rn の Cn を n 次元の極座標の角度積分として表せ．また Cn
は何を意味するか．

——————————————————————————
いま，r2 = s とおけば，2rdr = ds であるから∫n ∞ n/2−1 −s
In = Cn
s
e ds
2 0
(
)
n
n
= Cn Γ
2
2
6 この積分の値は∫ ∞2

e−x dx

I=
−∞

とおけば，

I2 =

∫ ∞ ∫ ∞−∞

22e−(x +y ) dxdy−∞と表される．ここで，2 次元の極座標を用いれば，

I2 =

∫ 2π

∫ ∞
dθ

0

r2 = s とおけば，ds = 2rdr であるから
I2 = π2re−r dr0∫ ∞e−s ds = π0√
したがって，I = π が得られる．
7 (4.4) で球を球殻に分けて積分していることは (4.3) の最後の等式が成り立つことからわかる．

25

(4.5)

ここで，Γ(x) はガンマ関数である．したがって，(4.2) と (4.5) から

Cn =

In
2In
2π n/2==(n/2)Γ(n/2)
nΓ(n/2)nΓ(n/2)(4.6)と求められる．
以上から，n 次元の球の体積 Vn は

Vn = Cn an(4.1)と求められた．そこで，係数 Cn は (4.6) で与えられる．

Cn は (問 4.3) で与えられた角度積分をしても得られる．その角度積分については [17] で述べる．
ところで，この結果は次元数が低いときに前に得られた結果と一致することを確かめておこう．そのためにはガンマ関数 [15] の性質を知っている必要がある．ガンマ関数は整数 n の階乗 n! を任意の実数に拡張したものであり，その定義は∫ ∞Γ(x + 1) =tx e−t dt

(4.7)

0

である．またΓ(x + 1) = xΓ(x)Γ(1) = 1
Γ(2) = 1√
Γ(1/2) = π = π 1/2
n = 0 または自然数

Γ(n + 1) = n!,が成り立つ．まず n = 1 のときには

C1 =

2π 1/2=2
Γ(1/2)

つぎに n = 2 のときには

C2 =

π
=π
Γ(1)

さらに，n = 3 のときには

C3 =4
2π 3/2= π
3Γ(3/2)3が得られ，これらは 3 節で求められた値と一致している．

5

おわりに
n 次元の球の体積を正統な方法ではなく少し頭を使った便法で求めた．計算は面倒になるかもしれないが，もっと正統な方法で求めることも必要であろう．これは [17] で示すことにする．一般の n 次元の球の体積を正統な方法で求めるための準備として具体的な 4, 5 次元の球の体積を [13] で求めることにしよう．

6

問の解答例
——————————————————————————[5][6] では n 次元の球の体積 Vn とその球の表面積 Sn を考えて，それらの関係を導いている．こ（問 4.1）の関係を導け．

26

——————————————————————————（問 4.1）解答例[5][6] によれば，Vn は球の半径 a だけの関数であり，
∫∫
∫
∫ a
Vn =
···
dx1 dx2 · · · dxn =
Sn dr
D

(6.1)

0

ここで領域 D は x21 + x22 + · · · + x2n ≤ a2 をみたす n 次元の球の中である．これは球の表面積に微小な球殻の厚さ dr をかけて積分すれば，球の体積 Vn が得られることを示している．
逆に (6.1) から，球の体積 Vn を球の半径 a で微分すれば，球の表面積

Sn =

dVn
da

(6.2)

が得られる．(4.3) の Cn [(r + dr)n − rn ] は半径 r + dr の球から半径 r の球を取り去ったものであるので，厚さ dr の球殻の部分の体積を表している．もちろん dr が微小なときには dr の 2 次以上の高次の項は無視していいので，(r + dr)n − rn = nrn−1 dr となる．ここで Cn nrn−1 は球の表面積を表す．

——————————————————————————（問 4.2） なぜ唐突に (4.2) の積分を考えたのか．
——————————————————————————（問 4.2）解答例 どうして (4.2) の積分を考えたのかは本当のところはわからない．しかし，
∫ ∞
√
2
I1 =
e−x dx = π

(6.3)

−∞

という積分はよく知られている．

[16] によれば
In = I1n
(∫ ∞
)n
−x2
=
e
dx
−∞
∫ ∞
∫ ∞
=
···
exp[−(x21 + x22 + · · · + x2n )]dx1 dx2 · · · dxn
−∞

(6.4)

−∞

であることはすぐにわかるので，この In の積分の値がすぐにわかる．さらにそれに n 次元の極座標を導入

すれば，In が極座標系を用いて別の方法でも積分できそうである．それで Cn が求められるのではないかという予想が立つ．

——————————————————————————
n 次元の極座標を定義して，Vn = Cn an の Cn を n 次元の極座標の角度積分として表せ．また Cn
は何を意味するか．
——————————————————————————（問 4.3）（問 4.3）の解答例

まず n 次元の極座標を定義しよう．n 次元の直交座標 (x1 , x2 , . . . , xn−1 ) と極座標(r, θ1 , θ2 , . . . , θn−1 ) の関係は次の式で与えられる．
x1 = r cos θ1 ,(6.5)

x2 = r sin θ1 cos θ2 ,(6.6)

x3 = r sin θ1 sin θ2 cos θ3 ,

(6.7)

...

(6.8)

xn−1 = r sin θ1 sin θ2 . . . sin θn−2 cos θn−1 ,(6.9)

xn = r sin θ1 sin θ2 . . . sin θn−2 sin θn−1 .

(6.10)

27

0 ≤ x21 + x22 + · · · + x2n ≤ a2 で表される球の内部にある (x1 , x2 , . . . , , xn ) に対する (r, θ1 , θ2 , · · · , θn−1 ) の変域は
0 ≤ r ≤ a,0 ≤ θ1 , θ2 , · · · , θn−2 ≤ π,

0 ≤ θn−1 ≤ 2πである．
これから体積要素 (volume element) は

dVn = dx1 dx2 · · · dxn−1 dxn= Jdrdθ1 dθ2 dθ3 · · · dθn−1(6.11)であり，ヤコビアン J がわかると n 次元の極座標系での体積要素が求められる．
このヤコビアンの導出は [17] に譲るが，，J = rn−1 sinn−2 θ1 sinn−3 θ2 · · · sin2 θn−3 sin θn−2

(6.12)である．したがって

1
Cn =
n

∫ π

∫ π
n−2

sin
0

∫ π
n−3

θ1 dθ1sinθ2 dθ2 · · ·0∫ 2π
sin θn−2 dθn−2

0

dθn−1

(6.13)0と表される．

n 次元の球の体積 Vn = Cn an で半径 a を 1 とおけば，Vn = Cn となり，Cn は半径 1 の n 次元の球（半径 1の球を単位球という）の体積である [7] [16]．
(2004.6.22) (2009.4.7 改訂) (2017.5.16 三訂)

参考文献[1] 矢野 忠，n 次元の球の体積 1，研究と実践（愛数協），88 号（2005.9) 7-12
[2] 矢野 忠，n 次元の球の体積 2，研究と実践（愛数協），90 号 (2006.5) 1-4[3] 矢野 忠，n 次元の球の体積 3，研究と実践（愛数協），91 号 (2006.7) 1-4[4] スミルノフ，『高等数学教程 II 巻（第 1 分冊）』（共立出版，1958）290-291[5] 芦田正巳，『統計力学を学ぶ人のために』（オーム社, 2006）341-343[6] EMAN の物理学・統計力学・n 次元球の体積,
http://homepage2.nifty.com/eman/statisitic/sphere_vol.html[7] 新関章三，岩崎正春，多次元単位球の体積計算法，数学・物理通信，1 巻 9 号 (2011.9) 13-17[8] 山本哲朗，『行列解析ノート』，（サイエンス社，2013）142-147[9] H. Goldstein, C. Poole, J. Safko（矢野，江沢，淵崎 訳），『古典力学』（吉岡書店，2009）702-711[10] S. C. Phatak and S. S. Rao, Logistic map: A possible random-number generator, Phys. Rev. E, Vol.
51, No. 4, 3670-3678 (1995)[11] 長瀬隆久，ロジスティック写像を用いた正規乱数発生法，情報処理学会論文誌，Vol. 42,No. 6, 1727-1729(2001) とそこに引用された論文を参照せよ．
28

[12] 宮武修，脇本和昌，『乱数とモンテカルロ』（森北出版，1978）55-58，ここには一般の n 次元の球の体積については記載がないが，５次元の多重積分の例がある．

[13] 矢野 忠，n 次元の球の体積 2，数学・物理通信，7 巻 4 号（2017.6) 12-19[14] 久保亮五，『統計力学』（共立出版，1952）46[15] 高木貞治，『解析概論』改訂第三版（岩波書店，1961）248-258[16] http://wasan.hatenablog.com/entry/20110503/1304372922[17] 矢野 忠，n 次元の球の体積 3，数学・物理通信，7 巻 5 号 (2017.6) 掲載予定

29

編集後記
伊予の国，松山は春の終りから夏にかけて，気温の高低の差が激しい日々が続き，毎日の体の調子を具合よく維持するのに少々苦労しております．皆さん，お元気ですか．
さて，今回の投稿論文もそれぞれの分野に応じて興味深い内容を提供して下さっており，編集者としてもうれしいかぎりである．そして，このような状態が今後とも続けばと神に祈りたい気持ちでいっぱいである．このような精神と，いつも頭にある言葉 “継続は力なり”をモット―にして編集にこれからも携わっていきたい．
それにしても編集者の一人が数学関係であるのに数学の論文が少ないのが，気になる．私自身の身体の調子もあり，歯がゆい限りではある．しかし，このサーキュラーの発刊の趣旨に沿った健全な発展のために是も非もなく努力したい．
（新関章三）

30