Shironetsu Blog

@shironetsuのブログ

一般フィボナッチ数列の周期について(2)

(1)の続き.shironetsu.hatenadiary.com
 (1)の記事で実数列としてのフィボナッチ数列について一通り説明した.
ここからは整数列としてのフィボナッチ数列の性質を見ていく.


用語の定義

 そのまえに, 整除性の話を始めるために, ここから成分は特に述べない限り(指数や位数, 添え字を除いて)剰余環Zm=Z/mZで考えることにする. つまり前回定義した一般フィボナッチ数列はZm上の数列に, 行列ΦGはZmに成分を持つ2×2行列の集合M(2;Zm)の元として扱うことになる. 等号はZmでの等号の意味で使うものと実数としての等号が混在することになるが, 文脈からどちらの意味か分かるようになっている(はず).

 これを踏まえ, まず「周期」や「系列」の意味をはっきりさせてしまおう.

定義2.1
 一般フィボナッチ数列{G[n]}の周期cは全てのnに対して
f:id:shironetsu:20151121172351p:plain
を満たす最小の自然数である.

 このように定義すれば当然その存在性が問題になるが, 鳩の巣論法と,Zm上一般フィボナッチ数列の2つ組が前後に一対一であることから明らかだろう. また, 行列でこれと同値な定義を与えるとその意味は一層明確になる.
f:id:shironetsu:20151121172457p:plain

 ΦがZm上可逆であることから結局周期とはこの最後の式を満たす最小の自然数であることが分かる. また, 同様の式を満たす自然数が全てcの倍数であり, その逆も成り立つこともほぼ明らかだろう.


 同様に出現点についても定義しておく.

定義2.2
 出現点eとは, F[e]=0 を満たす最小の自然数である.

 出現点の存在性は周期が存在することから明らかである. ただしこちらにはF[e+1]についての制約が無い.

 フィボナッチ数列では剰余について0を挟んだ両側は等しくなるから,行列を使って次のように言い換えることもできる.

出現点eとは, あるa≠0が存在して
f:id:shironetsu:20151121172526p:plain
となる最小の自然数である.ただしI単位行列.

 注意すべきことは, 周期cは初項に依存するものであるのに対して出現点は(0,1)から始まるフィボナッチ数列にのみ意味を持つことだ. 思い出しておくと, この記事で述べることはその出現点, フィボナッチ数列で最初に0が出てくる番号さえ分かれば他の系列についても大体の性質が分かるということだった.

定義2.3
(a,b)∈Zm×Zmに対して, (a,b)系列*1とは
f:id:shironetsu:20151121172710p:plain
によって定まるZm×Zm上の同値関係"~Φ"による同値類である.

 一応順序をなくした同値類として「系列」の意味を定めたが, 1つ代表元を取ってΦをかければすぐに復元できる. また, この言葉を使うと周期は系列の元の個数ということになる.

 ところで0を含む系列と含まない系列があることを見たのだった. これらをそれぞれ「節有り」と「節無し」と呼び, 1周期の中に0を含む数を「節数」と呼ぼう.

定義2.4
 (0,0)系列でない(a,b)系列Cの節数sとは
f:id:shironetsu:20151121172735p:plain
によって定まる非負整数である. このとき節有り系列は s>0 の系列であり, 節無し系列とは s=0 の系列である.

 この記事で使う言葉についての定義はこれで終わった.

p元体上の一般フィボナッチ数列

 一応一般の自然数に対して以上の用語を定めたが, ここからは専ら法は素数に限ることにする. そうするのはpに対して剰余環Zpがp元体(有限素体ともいう;記号Fpを使う)になって格段に調べやすくなるからに他ならない.

 目標としてここで証明することを先に掲げてしまおう.

以下でpは5以外の奇素数とする.

(I) p≡±1 mod 10のとき F[p-1]≡0, F[p]≡1 mod p
p≡±3 mod 10のとき F[p]≡-1, F[p+1]≡0 mod p

(II)出現点eと節有り系列の節数sについて,
e≡0 mod 4のとき s=2
e≡±1 mod4 となるのは p≡1 mod 4のときに限り, s=4
e≡2 mod 4のときs=1
なお, 節有り系列の節数は全て一致している.

(III)p≡±3 mod10,またはe≡0,1 mod 4のとき
   (0,0)系列を除く全ての系列の周期は(0,1)系列の周期に一致する.
   p≡±1 mod10かつe≡2 mod 4のとき(0,0)系列を除いて
   (0,1)系列の周期cより短い周期の系列が存在し, その周期はc/2

 これらを使えば前回の記事の表が, 出現点を調べるだけで書けてしまう.

 行列を使った表記については既に述べてきたが, ここからの話はp元体Fp上の2次正方行列が「主役」を演じることになる. *2

 「系列」はFFp上の同値関係から定められていたが, その中の元(a,b)と次のFp上の2次正方行列Aとの全単射により, 行列間の同値関係が誘導される.
f:id:shironetsu:20151121173016p:plain
 さらにこの行列をよく見ると,
f:id:shironetsu:20151121173038p:plain
と分解できることが分かる. なお行列Aは以降もこの意味で用いることにする.

 この関係からさらに興味深い構造が見えてくるが, とりあえずは素朴な数論から
示すべきことについて述べていこう.

 まず周期と出現点,節の関係について簡単だが重要な関係を示す.

命題2.1
素数pに対して, (0,1)系列の周期cと出現点e,節数sについて
c=seで, s=1,2,4のいずれか.

(0,1)系列がフィボナッチ数列に対応することに注意.

証明
(0,1)系列の周期は
f:id:shironetsu:20151121173116p:plain
を満たす最小の自然数であり, また出現点eはあるa∈Fpが存在して
f:id:shironetsu:20151121173127p:plain
を満たす最小の自然数である.
cをeで割りc=qe+r (0≦r<e)とすると
f:id:shironetsu:20151121173137p:plain
r≠0だとeの定義に反するためr=0でcはeの倍数. さらにqはa^q=1を満たす最小の自然数になる. F^e=aI, F^2e=a^2I,…,F^qe=Iよりqは節数sになるため, c=se. eとaの関係について, 行列式を比較すると,(-1)^e=a^2
(1)eが奇数のとき
 a^2=-1から a≠1, a^2=-1≠1, a^3≠1, a^4=1よりs=4.
(ii)eが偶数のとき
 a^2=1より a=1またはa=-1. よってs=1,2

(証明終)

 「F[e]=0ならF[4e+1]=1」はもっと素朴に示すこともできる.
フィボナッチ数の「加法定理」から,
f:id:shironetsu:20151121173204p:plain
行列で計算すればすっきりするところをまわりくどくやっているだけではあるが,
フィボナッチ数列の関係式に立ち返ってみるのも悪くはない.

 上で定義した行列A[a,b]の行列式は, 系列を分類する上では本質的な意味を持つ. 同じ系列に含まれるという意味で同値な2つのFFpの元は, 対応する行列式が等しいかちょうど(-1)倍になっている. (a,0)が含まれる系列では, 対応する行列式は±a^2であるから,(b,0)もまた含まれるならb^2=±a^2でなければならない. そのためこれを満たすbはaを含めて高々4個に制限されることになる.

 これで目標に掲げた(II)は部分的に示されたことになる. 先に示してしまったのは, 簡単なため.

 では目標に掲げた(I)を示そう.
カタランの公式を思い出すと, 奇素数pに対してp番目のフィボナッチ数は,
f:id:shironetsu:20151121173242p:plain
ただし今度の等号はp元体の元としての意味である. フェルマーの小定理から2^{p-1}=1, さらに2項係数について
f:id:shironetsu:20151121173254p:plain
の関係を思い出すと,
f:id:shironetsu:20151121173305p:plain

 再びフェルマーの小定理からこの値は±1であることは分かるが, そのどちらだろうか.ここで登場するのがルジャンドル記号である. つまり平方剰余が関わってくることになる.

平方剰余

 数論の基本的な内容ではあるが, 極めて重要なため平方剰余について一通り述べておこう. ここでは合同式の記号を使うことにする. またpは奇素数とする.

定義2.5
 pの倍数でないaに対して x^2≡a mod pを満たす整数xが存在するとき, aを(pに対する)平方剰余, 存在しないとき平方非剰余という.

 たとえば7に対して
1^2≡1, 2^2≡4, 3^2≡2, 4^2≡2, 5^2 ≡4, 6^2=1
から, 1,2,4は平方剰余, 3,5,6は平方非剰余である.

定義2.6
 ルジャンドル記号(a/p)は
f:id:shironetsu:20151121173358p:plain
によって定められる.
法pにおいて平方剰余と平方非剰余はともに(p-1)/2個ずつ存在する.

命題2.2
 ルジャンドル記号に対して次の関係が成立する.ここでa,bはpの倍数でない整数, pとqは相異なる奇素数とする.
f:id:shironetsu:20151121173455p:plain
一番上の法則は
(平方剰余)×(平方剰余)=(平方剰余)
(平方非剰余)×(平方非剰余)=(平方剰余)
(平方剰余)×(平方非剰余)=(平方非剰余)
となることを言っている.

 平方剰余の相互法則はいわゆる初等数論の中でもっとも重要な定理の一つで, それを認識したガウスがそのいくつもの異なる証明を与えたことで有名である.

 平方剰余についてはとりあえずこれで十分. では中断していた証明に戻ろう
オイラーの規準からp元体Fpの元として
f:id:shironetsu:20151121173519p:plain
となるから, 5が平方剰余であるか否かを見ればいい. 5以外の奇素数pについて相互法則から,
f:id:shironetsu:20151121173530p:plain
法10に置き換えられるのはpが奇数のため.

結局, 5以外の奇素数pに対して
f:id:shironetsu:20151121173542p:plain
が示された.

 F[p-1],F[p+1]に対しても同じようにしてpに対する剰余を求められるが, やや煩雑である. そこで行列を使った方法を試みよう. *3

 Φに対して, ケーリー・ハミルトンの法則から,
f:id:shironetsu:20151121173606p:plain
1行目から2行目で「平方完成」, 2行目から3行目で両辺(p-1)/2乗,
3行目から4行目で両辺に(2Φ-I)をかけた.
4行目の左辺を2項展開すると,
f:id:shironetsu:20151121173617p:plain
pについて場合分けすると,

(1)(5/p)=+1 ; 5が平方剰余であるとき
f:id:shironetsu:20151121173637p:plain
(2)(5/p)=-1; 5が平方非剰余であるとき
f:id:shironetsu:20151121173649p:plain
ルジャンドル記号(5/p)を使ってまとめると,
f:id:shironetsu:20151121173707p:plain
とも書ける.

 この結果は重要だ. 5が平方剰余のとき出現点も周期もp-1の約数になる.5が平方非剰余のとき出現点はp+1の約数, 周期は2p+2の(p+1以外の)約数になる.表を見ればこのことが確かめられるだろう.

 実は出現点そのものについて(この記事では)これ以上強い主張ができない. 前の記事の表を再び見ると, たとえばp=83では出現点はp+1=84, 周期は2p+2=168でたっぷりかかるパターンになっている. その一方p=89は11(p-1=88の約数)が出現点で, 周期も44と短い.表には載っていないが, 23番目のフィボナッチ数28657は素数であるため, 出現点23(28657+1を割り切って商は1246),周期92ともっと短い*4.

 この出現点と周期の性質については後で再び別証明が与えられることを見るだろう. では続いて目標に掲げた(II)(III)についてまとめて証明に入る.

対角化, ジョルダン標準形, 体の拡大

 (a,b)系列の周期cは
f:id:shironetsu:20151121173748p:plain
を満たす最小の自然数なのだった.

 行列で表したのだからΦを対角化しよう.固有値をλとして固有方程式を求めると,
f:id:shironetsu:20151121173808p:plain
思い出さねばならないのはこれがp元体Fp上での方程式だということだ. 5が平方剰余のとき2つの異なる解が存在するが, pが平方非剰余のときFp上に解は存在せず, p=5のときは重解になる.

それぞれ場合分けして考える.

(1)5が平方剰余 ; p≡±1 mod 10のとき
 x^2=5はp元体上2つの解を持つからその一方を√5と書いてしまう. このときもう一方は-√5. そして分数を分母の逆元をかける意味で使うと結局実数の時と同じように2つの固有値
f:id:shironetsu:20151121173845p:plain
になる.2行目は解と係数の関係である.
たとえばp=19のときx^2=5の解は9,-9=10の2つで, 固有値は5,15の2つ.

 対角行列Λを次のように定めて, 各固有値に対応する固有ベクトルを並べた行列Pを用いてΦを対角化できる.
f:id:shironetsu:20151121173908p:plain
これを使ってΦのべき乗をとると周期cについて次のように書ける.
f:id:shironetsu:20151121173922p:plain
ここでA[a,b]=aI+bΦという関係があったことを思い出すと,
f:id:shironetsu:20151121173935p:plain
この成分を見てさらに場合分けする.

(イ)a+bα≠0かつa+bβ≠0のとき
(a+bα)(a+bβ)=a^2-ab-b^2≠0と表せて*5, cはα^c=1かつβ^c=1を満たす最小の自然数になる. ゆえにcはαとβの(Fpの乗法群での)位数の最小公約数ということになる.

(ロ)a+bα=0のとき
b=-a/α=βaで, この系列は a,βa,β^2a,β^3a… という「等比数列」になる. このときa+bβ=(2+β)a≠0(a,bがともに0とはなっていないことを仮定)だから, 周期はβの位数.

(ハ)a+bβ=0のとき
(ロ)と同様に周期はαの位数.

 αとβの位数の関係を考える. 位数の大きくないほうをα,その位数をkとする.*6

(i)kが奇数のとき
f:id:shironetsu:20151121174052p:plain
仮定よりβの位数はk以上であるため位数は2kになる. これらの公約数は2k. 結局,
a^2-ab-b^2≠0またはb=βaのとき周期2k,
b=αaのとき周期kと他のちょうど半分の長さになる.

(ii)kが偶数のとき
同様にして
f:id:shironetsu:20151121174113p:plain
仮定からβの位数はk.
a^2-ab-b^2≠0,b=αa,b=βaいずれの場合にも周期はkになる.

 続けて出現点についても考える.出現点eはあるd≠0が存在して次の条件を満たす最小の自然数である.
f:id:shironetsu:20151121174125p:plain

(i)eが奇数のとき
(α^e)^2=-1より-1が平方剰余になるから, 第一補充法則よりp≡1 mod 4
eはこれを満たす最小の自然数だから, αの位数はk=4e. すなわちαの位数が偶数の場合に相当し, 節数は4.

(ii)eが偶数のとき
α^2e=1となるが,α^e=-1, 1, またe≡0,2 mod 4(それぞれ自然数e'によりe=4e', e=4e'+2とおく)のときとでさらに分ける.

  (ii-1) α^e=-1 かつ e≡0 mod 4のとき
   k=2eでαの位数が偶数の場合に相当し, 節数は2.

  (ii-2) α^e=-1 かつ e≡2 mod 4のとき
  f:id:shironetsu:20151121174211p:plain
  よりeの最小性に反す.

  (ii-3) α^e=1 かつ e≡0 mod 4のとき
  f:id:shironetsu:20151121174254p:plain
  よりeの最小性に反す.

  (ii-4)α^e=1 かつ e≡2 mod 4のとき
  f:id:shironetsu:20151121174320p:plain
  よってαの位数はk=2e'+1=e/2で, 奇数の場合に相当し節数は1

 まとめると
eが奇数のときk=4e, すべての系列で周期は等しく4e, 節数4
e≡0 mod4のときk=2e, すべての系列で周期は等しく2e, 節数2
e≡2 mod4のときk=e/2, (a, αa)系列で周期e/2, その他の系列は周期eで節数1

 ちなみにαとβに対してフェルマーの小定理を使うと,
f:id:shironetsu:20151121174353p:plain
が言える.

(2)p=5のとき
 1つしかないわけだからわざわざ行列で見なくてもいいが統一的に扱う*7.
 固有方程式の解は
f:id:shironetsu:20151121174420p:plain
で重根を持ち対角化できないが, ジョルダン標準形には変形できる.
行列PJが存在して, (1)と同様にして(a,b)系列の周期cを次のように表せる.
f:id:shironetsu:20151121174437p:plain

 (イ)a+bλ≠0; b≠3aのとき
  行列の対角成分からλ^c=1. これを代入すると(1,2)成分から
  f:id:shironetsu:20151121174455p:plain
  (「微分のような形」になって指数が肩から降りてきているため合同記号を使った.)
  λの位数は3^2=4, 3^3=2, 3^4=1より4だから, cは4と5の最小公倍数で20.

   出現点eについて考えると, あるdが存在して
  f:id:shironetsu:20151121174512p:plain
  を満たす最小の自然数で, c=4e. 節数4.

 (ロ)a+bλ=0; b=3aのとき
  行列の(1,2)成分からbλ^c=b. cはλの位数で4

……と, 特殊なp=5の場合もジョルダン標準形によって扱うことができた.ちなみにPを実際に求めてΦのべき乗をとるとF[n]=n・3^(n-1)と表せることがわかる.

(3)5が平方非剰余 ; p≡±3 mod 10のとき
 固有方程式はp元体の中に解を持たない. 言い換えると固有多項式Fp係数多項式環Fp[X]において既約である. そこで2次方程式の解が実数の中に存在しなかったとき虚数単位を導入したように, この固有多項式の根を「新しい数」として付け加えてしまうことを考える. 代数の言葉で言うと X^2-X-1 の生成するイデアル(X^2-X-1)による剰余環Fp[X]/(X^2-X-1)$を考えるのだ. これは体になって, Fp上の2次元ベクトル空間, 拡大次数2の拡大体である. Xの代表する類 X+(X^2-X-1) をφとすると, φ^2-φ-1=0でこれが根の1つになる. <1,φ>はベクトル空間としての基底になり, 新しい体はFp+Fpφで表せる.
これはガロア体GF(p, 2)とも呼ばれるものである. この中で固有方程式はもう1つの解 1-φ を持つためこちらは φ* としよう. また, ここでは単純にFp^2とこの体を表すことにしよう.

 まわりくどい導入をしたが, 結局「5の平方根」をp元体上に付け加えるのと同じことである.
f:id:shironetsu:20151121174723p:plain
平方剰余と平方非剰余の積は平方非剰余だったから, 全ての平方非剰余がFp^2の中では異なる根を2つずつ持てることになる.

 これを使うことで, 行列ΦFp^2上の行列として対角化することができる. 手順としては(1)の5が平方非剰余だったときとほとんど同じで, ただ単にαとβをφとφ*で置き換えれば済むので省略し, 周期cについて以下のような表し方を得る.
f:id:shironetsu:20151121174816p:plain
ただし今度はa+bφもa+bφ*も0とはならない. というのも, 仮に0になればaとbはともに0とはならないFpの元であるという仮定に反するためである. 1とφがFp上のベクトル空間の基底として一次独立だったことを思い出そう.

 つまり周期cは全ての系列でφとφ*の, 今度はFp^2の乗法群での位数の最小公倍数ということになる.

 出現点と周期の関係については短い周期の系列が現れない点のみ異なり, ほかは(1)と同様.

 ちなみに「p≡±3 mod 10のとき一般フィボナッチ数列の周期はフィボナッチ数列の周期に一致する」という命題はアメリカの数学者ドナルド・ウォールの名前から「ウォールの定理」と呼ばれるらしい.

 なお以下のようにして命題(II)が再度示される.
Fp^2の乗法群は元の個数が p^2-1 の巡回群であるから,
f:id:shironetsu:20151121174903p:plain
一方,
f:id:shironetsu:20151121174911p:plain
よって,
f:id:shironetsu:20151121174920p:plain
が言える.
これはφとφ*の位数が2(p+1)で割り切れ, (p+1)ではないことも意味する.


結論

 以上の考察から出現点による系列の分類ができた. まとめると, 出現点e, 周期c, (0,1)系列の節数eについて以下の関係が成り立つ. 周期が系列に依存する値であることに注意.

素数pに対して

  • p≠5のとき
    • e ≡ ±1 mod 4 のとき全系列でc=4e, s=4
    • e ≡ 0 mod4 のとき全系列でc=2e, s=2
    • e ≡ 2 mod4のとき
      • e≡±1 mod 10 のとき λ^2-λ-1=0の解がFp上に2つ存在し,位数が小さいほうαに対し(a, αa)系列でc=e/2, その他の系列はc=eでs=1
      • e≡±3 mod 10 のとき 全系列でc=e, s=1
  • p=5のとき

   e=5,s=4で,c=20の系列とc=4の系列が存在.

系列の個数については単にFFpから(0,0)を除いたp^2-1の元を周期で分けていけばいい.

目標に達したのでここで終える.
ここまで来ればどれも当たり前のように思えるが, フィボナッチ数列という単純な再帰数列を通じて整数の色々な性質や代数的な構造が現れてきたことには驚かされる. フィボナッチ数列やはりおもしろい.

 補遺でリュカ数列などについて述べる.

*1:「軌道」と言うほうがグローバルだったかもしれないが, 自分の中で「系列」で定着しすぎてしまったのでこう呼ぶ. しかしあまりローカルな言葉を生やしたくない気持ちもある.

*2:Zmは可換環だからZm上の行列から定義してもいいのだが

*3:この証明は中村滋『フィボナッチ数の小宇宙』(日本評論社, 2002年)p.215に載っている証明を部分的に変えたものである. その他この記事の多くの部分もこの本によるところが大きい.

*4:ちなみにフィボナッチ数列中に無限に素数が存在するか否かは未解決問題である(はず).

*5:これは単にdet(A)≠0を意味する.「行列式が0で無いならフィボナッチ数の周期と一致する」ということ自体はすぐに言えるがそうすると固有値が見えない.

*6:上でx^2=5の解の一方を√5として定めたが, 正負で分ける実数の場合と違ってどちらの解であるかを制限しない.

*7:一般の2階線形回帰数列を考えるときは必要になる.これについては次の記事で述べる.

一般フィボナッチ数列の周期について(1)

 祖先型の目で眺めたこの惑星は, 十分近づいても赤道と極以外に陸地が無いように見える.しかし周回軌道上でよく目を凝らせばその勘違いが, 青色と海を結びつける先入観に起因するものだと気付く. あるところでは群れをなし, またあるところでは毛羽立ち, 赤道では渦巻く白い雲の下にいくらか隠されながらも滑らかで穏やかな表面を見せるのが海洋だ.その一方で, 皺がよったようにわずかずつ異なる「青色」が縦横に走る領域があることが分かるだろう. 海に見えていた部分の3割は青い植物で覆われた大陸なのだ.

 伝統的な分類で主系列のF型星に属すこの系の太陽は我々の故郷のソルよりいくらか重く,その放つ光はやや明るい. そしてスペクトルのピークは短波長側にずれている. つまり少しだけ青い. 植物が持つクロロフィル――と相似な分子――が, もっともエネルギーに富む波長を選んでいないように見えるのはここでもどうやら同じようだ.

 といっても光合成の担い手が全て青いわけではない. やはり地球同様に植物の上陸が一度だけ起こったこの惑星では,現生の陸生青色植物の光合成色素は数億年間経った今でも高い一様性を保っている. しかし水中に目を転ずれば進化の試行錯誤の跡と多様性を見ることができるだろう. いくつもの系統樹の枝に渡って存在する単細胞の藻類は, この惑星を酸素で改造した祖先と微妙に異なる多様な光合成色素を持ち, 淡水海水に限らずあらゆる水圏で食物連鎖を支えていた. カルシウムの骨格を獲得したある藻類のグループはほとんど黒に近い紺色で, 浅瀬の動物に住処を与え生態系にとって重要な役割を果たしていた.

 我々が出会った緋色藻は, 高緯度地域の湖で赤く小さな体を漂わせる,そんな多くの興味深い藻類の中の1つだった.

[ライブラリ>数学>数論>フィボナッチ数列]

 0,1,1,2,3,5,8,13,……言わずと知れたフィボナッチ数列.少し計算すると15番目の610で初めて最後の桁が0になる.さらに根気よく計算を続けていけば最後の桁の数がある周期を持つことに気付くだろう.と言っても21番目の10946に至ってもまだ一巡りしない.そこで最後の桁の数だけで計算を続けてみよう.

0,1,1,2,3,5,8,3,1,4,5,9,4,3,7,0,7,7,4,1,5,6,1,7,8,5,3,8,1,9,
0,9,9,8,7,5,2,7,9,6,5,1,6,7,3,0,3,3,6,9,5,4,9,3,2,5,7,2,9,1,
0,1,1,…

と60番目で0, 61番目で1が出てきて以降は繰り返されることが分かる.すなわち周期は60である.*1.

 さて, 「最後の桁の数」は我々の手が10本(十進数で)の指を備えている事情から
見えているに過ぎない. 他の自然数*2
に対する剰余の周期も気になるところである.

mod 2 : 0,1,1,0,1,1,…
mod 3 : 0,1,1,2,0,2,2,1,0,1,1,….
mod 4 : 0,1,1,2,3,1,0,1,1,…
mod 5 : 0,1,1,2,3,0,3,3,1,4,0,4,4,3,2,0,2,2,4,1,0,1,1,…

こうなると更に気になるのは0,1以外から始めた場合だ.とりあえず法を素数に限るのが常套手段*3.

 上の列を見ると2,3では全ての組が尽くされているが, mod5で出てきていない組があることに気付く.探すとそれは0が出てこない周期4の列,

1,3,4,2,1,3,…

であることが分かる.これら2つでmod5は尽くされるのだ. ただし0,0,0,…は自明だからこの例に限らず除くとする.

 mod7ではどうか.

● 0, 1, 1, 2, 3, 5, 1, 6, 0, 6, 6, 5, 4, 2, 6, 1, 0, 1, …
● 0, 2, 2, 4, 6, 3, 2, 5, 0, 5, 5, 3, 1, 4, 5, 2, 0, 2, …
● 0, 3, 3, 6, 2, 1, 3, 4, 0, 4, 4, 1, 5, 6, 4, 3, 0, 3, …

これで7*7-1個の2つ組は尽くされた.

 mod11では?

● 0, 1, 1, 2, 3, 5, 8, 2, 10, 1, 0, 1, …
● 0, 2, 2, 4, 6, 10, 5, 4, 9, 2, 0, 2, …
● 0, 3, 3, 6, 9, 4, 2, 6, 8, 3, 0, 3, …
● 0, 4, 4, 8, 1, 9, 10, 8, 7, 4, 0, 4,…
● 0, 5, 5 ,10, 4, 3, 7, 10, 6, 5, 0, 5,…
● 0, 6, 6, 1, 7, 8, 4, 1, 5, 6, 0, 6, …
● 0, 7, 7, 3, 10, 2, 1, 3, 4, 7, 0, 7, …
● 0, 8, 8, 5, 2, 7, 9, 5, 3, 8, 0, 8, …
● 0, 9, 9, 7, 5, 1, 6, 7, 2, 9, 0, 9, …
● 0, 10, 10, 9, 8, 6, 3, 9, 1, 10, 0, 10, …

0の入る場合はこれで全てだ. しかし周期10*10=100組だけしか入っていない.
11*11のマス目に印を付けてみると, 出てこない組が分かる.

● 1, 4, 5, 9, 3, 1, 4, …
● 1, 8, 9, 6, 4, 10, 3, 2, 5, 7, 1, 8, …
● 2, 8, 10, 7, 6, 2, 8, …

周期5+10+5=20組. これで全てだ.

 mod13では

● 0, 1, 1, 2, 3, 5, 8, 0, 8, 8, 3, 11, 1, 12, 0, 12, 12, 11, 10, 8, 5, 0, 5, 5, 10, 2, 12, 1, 0, 1, …
● 0, 2, 2, 4, 6, 10, 3, 0, 3, 3, 6, 9, 2, 11, 0, 11, 11, 9, 7, 3, 10, 0, 10, 10, 7, 4, 11, 2, 0, 2, …
● 0, 4, 4, 8, 12, 7, 6, 0, 6, 6, 12, 5, 4, 9, 0, 9, 9, 5, 1, 6, 7, 0, 7, 7, 1, 8, 9, 4, 0, 4, …

0を含む列は周期28が3つ, ここに入る2つ組は84. 0を含まないものは

● 2, 1, 3, 4, 7, 11, 5, 3, 8, 11, 6, 4, 10, 1, 11, 12, 10, 9, 6, 2, 8, 10, 5, 2, 7, 9, 3, 12, 2, 1, …
● 4, 2, 6, 8, 1, 9, 10, 6, 3, 9, 12, 8, 7, 2, 9, 11, 7, 5, 12, 4, 3, 7, 10, 4, 1, 5, 6, 11, 4, 2, …
● 8, 4, 12, 3, 2, 5, 7, 12, 6, 5, 11, 3, 1, 4, 5, 9, 1, 10, 11, 8, 6, 1, 7, 8, 2, 10, 12, 9, 8, 4, …

周期28が3つ. 0を含むものと同じ形だ.

 観察のため, mod17, 19も続けてみる.

mod17

0を含むもの
●0, 1, 1, 2, 3, 5, 8, 13, 4, 0, 4, 4, 8, 12, 3, 15, 1, 16, 0, 16, 16, 15, 14, 12, 9, 4, 13, 0, 13, 13, 9, 5, 14, 2, 16, 1, …
●0, 2, 2, 4, 6, 10, 16, 9, 8, 0, 8, 8, 16, 7, 6, 13, 2, 15, 0, 15, 15, 13, 11, 7, 1, 8, 9, 0, 9, 9, 1, 10, 11, 4, 15, 2, …
●0, 3, 3, 6, 9, 15, 7, 5, 12, 0, 12, 12, 7, 2, 9, 11, 3, 14, 0, 14, 14, 11, 8, 2, 10, 12, 5, 0, 5, 5, 10, 15, 8, 6, 14, 3, …
●0, 6, 6, 12, 1, 13, 14, 10, 7, 0, 7, 7, 14, 4, 1, 5, 6, 11, 0, 11, 11, 5, 16, 4, 3, 7, 10, 0, 10, 10, 3, 13, 16, 12, 11, 6, …
周期36*4列.

0を含まないもの
●1, 3, 4, 7, 11, 1, 12, 13, 8, 4, 12, 16, 11, 10, 4, 14, 1, 15, 16, 14, 13, 10, 6, 16, 5, 4, 9, 13, 5, 1, 6, 7, 13, 3, 16, 2, …
●1, 4, 5, 9, 14, 6, 3, 9, 12, 4, 16, 3, 2, 5, 7, 12, 2, 14, 16, 13, 12, 8, 3, 11, 14, 8, 5, 13, 1, 14, 15, 12, 10, 5, 15, 3, …
●1, 7, 8, 15, 6, 4, 10, 14, 7, 4, 11, 15, 9, 7, 16, 6, 5, 11, 16, 10, 9, 2, 11, 13, 7, 3, 10, 13, 6, 2, 8, 10, 1, 11, 12, 6, …
●1, 9, 10, 2, 12, 14, 9, 6, 15, 4, 2, 6, 8, 14, 5, 2, 7, 9, 16, 8, 7, 15, 5, 3, 8, 11, 2, 13, 15, 11, 9, 3, 12, 15, 10, 8, …
周期36*4列.

mod19

0を含むもの
●0, 1, 1, 2, 3, 5, 8, 13, 2, 15, 17, 13, 11, 5, 16, 2, 18, 1, …
●0, 2, 2, 4, 6, 10, 16, 7, 4, 11, 15, 7, 3, 10, 13, 4, 17, 2, …
●0, 3, 3, 6, 9, 15, 5, 1, 6, 7, 13, 1, 14, 15, 10, 6, 16, 3, …
●0, 4, 4, 8, 12, 1, 13, 14, 8, 3, 11, 14, 6, 1, 7, 8, 15, 4, …
●0, 5, 5, 10, 15, 6, 2, 8, 10, 18, 9, 8, 17, 6, 4, 10, 14, 5, …
●0, 6, 6, 12, 18, 11, 10, 2, 12, 14, 7, 2, 9, 11, 1, 12, 13, 6, …
●0, 7, 7, 14, 2, 16, 18, 15, 14, 10, 5, 15, 1, 16, 17, 14, 12, 7, …
●0, 8, 8, 16, 5, 2, 7, 9, 16, 6, 3, 9, 12, 2, 14, 16, 11, 8, …
●0, 9, 9, 18, 8, 7, 15, 3, 18, 2, 1, 3, 4, 7, 11, 18, 10, 9, …
●0, 10, 10, 1, 11, 12, 4, 16, 1, 17, 18, 16, 15, 12, 8, 1, 9, 10, …
●0, 11, 11, 3, 14, 17, 12, 10, 3, 13, 16, 10, 7, 17, 5, 3, 8, 11, …
●0, 12, 12, 5, 17, 3, 1, 4, 5, 9, 14, 4, 18, 3, 2, 5, 7, 12, …
●0, 13, 13, 7, 1, 8, 9, 17, 7, 5, 12, 17, 10, 8, 18, 7, 6, 13, …
●0, 14, 14, 9, 4, 13, 17, 11, 9, 1, 10, 11, 2, 13, 15, 9, 5, 14, …
●0, 15, 15, 11, 7, 18, 6, 5, 11, 16, 8, 5, 13, 18, 12, 11, 4, 15, …
●0, 16, 16, 13, 10, 4, 14, 18, 13, 12, 6, 18, 5, 4, 9, 13, 3, 16, …
●0, 17, 17, 15, 13, 9, 3, 12, 15, 8, 4, 12, 16, 9, 6, 15, 2, 17, …
●0, 18, 18, 17, 16, 14, 11, 6, 17, 4, 2, 6, 8, 14, 3, 17, 1, 18, …
周期18*18列

0を含まないもの
●1, 5, 6, 11, 17, 9, 7, 16, 4, …
●1, 15, 16, 12, 9, 2, 11, 13, 5, 18, 4, 3, 7, 10, 17, 8, 6, 14, …
●2, 10, 12, 3, 15, 18, 14, 13, 8, …
周期9が2つ, 周期18が1つ.

こうして課題が見えてきた.
- 周期を決めるのは何か?
- 0が初めて現れるのはいつか?
- 0を含まない系列があるのはいかなる場合か?
- 0を含むとすれば循環節中にいくつ含むか?
- 他より短い周期の系列が生じるのはいつか?

 (0,1)から始まる系列において0が初めて出てくる番号は出現点(エントリーポイント)と呼ばれ, 周期と同様に重要である.また, 0を含む系列が1つ存在すれば, 各成分に対してある数をかけることで他の系列が全て得られることはほとんど明らかだろう.このことを踏まえ, 愚直にp*pのマス目にチェックを入れながら各系列について計算する方法で上の
項目についてもっと調べてみよう.

法p 0を含む系列 周期 出現点 0の数 0を含まない系列
2 1 3 3 1 0
3 1 8 4 2 0
5 1 20 5 4 1
7 3 16 8 2 0
11 10 10 10 1 2+1
13 3 28 7 4 3
17 4 36 9 4 4
19 18 18 18 1 2+1
23 11 48 24 2 0
29 28 14 14 1 32+2
31 30 30 30 1 2+1
37 9 76 19 4 9
41 20 40 20 2 22
43 21 88 44 2 0
47 23 32 16 2 46
53 13 108 27 4 13
59 58 58 58 1 2+1
61 15 60 15 4 47
67 33 136 68 2 0
71 70 70 70 1 2+1
73 18 148 37 4 18
79 78 78 78 1 2+1
83 41 168 84 2 0
89 22 44 11 4 158
97 24 196 49 4 24

 「0を含まない系列」の項の"+1","+2"は短い周期のものがその数だけ現れていることを示す. 表には書いていないが,その長さは他の系列のちょうど半分になっている.つまり短い系列はそれほどたくさん現れずどうやらその長さは他の半分らしい.

 ……ただしp=5の場合は他とは異なる. 上で見たように, 0を4つ含む周期20の系列と1,3,4,2,…の周期4の系列が現れる. これは端的に5という数がフィボナッチ数において特殊な役割を果たしていることを表している.

 今から示していくのは, 実は出現点さえ分かればこの表の他の項目は埋められるということだ. p*pのマス目を埋めていく必要はないのだ.

 そのためにフィボナッチ数について定義や重要ないくつかの命題をまとめよう.

定義1.1 フィボナッチ数列{F[n]}とは,
f:id:shironetsu:20151114184328p:plain
によって定まる再帰数列である.

 この定義では添え字は非負整数に限っているが負数番への拡張は容易だ.単にn<0に対してはF[n]=F[n+2]-F[n+1]とすればいい. そして単に偶数番に負号が付くだけであることもすぐに分かる.一応命題として述べておこう.

命題1.1 n>0に対して
f:id:shironetsu:20151114184343p:plain

 しかし初項が(0,1)だけではあまりに窮屈だ.同じ形の再帰数列を, 初項を一般の整数の組にも拡げて「一般フィボナッチ数列」とまとめて呼ぼう.

定義1.2 初項が(a,b)∈Z^2の一般フィボナッチ数列{G[n]}を
f:id:shironetsu:20151114184400p:plain
によって定める.

 後に分かることだが, a=2, b=1の場合はフィボナッチ数列と「対」になると言える重要な数列である.これをリュカ数列*4と呼ぶ.

定義1.3

リュカ数列[L[n]}とは
f:id:shironetsu:20151114184415p:plain
によって定まる数列である.

最初の数項は2,1,3,4,7,11,18,29,…となる.

 さて, 一般フィボナッチ数列も何項か計算してみればすぐにその規則に気付くだろう.

a, b, a+b, a+2b, 2a+3b, 3a+5b, …

係数にフィボナッチ数が現れている. これを命題として述べよう.

命題1.2 初項(a,b)の一般フィボナッチ数列{G[n]}に対して
f:id:shironetsu:20151114184430p:plain
証明は略すが数学的帰納法からすぐに示せる.

 より一般的な形として

命題1.3 
f:id:shironetsu:20151114184441p:plain
さらにここから系として

命題1.4
f:id:shironetsu:20151114184456p:plain
も簡単に示される.

 この関係式から感じられるように, フィボナッチ数列は行列を使うと見通しが良くなる.線型回帰数列であるから当然のことではある.簡単のためこの2つ組を「進める」行列に記号Fを当てておき以後断りなく使う.

命題1.5
行列Φ,横ベクトルg[n]を以下のように定めると次の関係式が成り立つ.
f:id:shironetsu:20151121140728p:plain

さらに行列G[n]を定めると, 
f:id:shironetsu:20151121140738p:plain

特に,
f:id:shironetsu:20151121140805p:plain

ベクトルに対して後ろから行列を書ける形式にしたのは横ベクトルで使いたいことが理由で, 特に深い意味は無い.

行列式を比較することから次の命題が得られる.

命題1.6
f:id:shironetsu:20151114184717p:plain

3番目の等式はシムソンの等式と呼ばれている*5.

 この行列Fの対角化と累乗により次の有名な等式、ビネーの公式*6*7が得られる.

命題1.7(ビネーの公式)

λ^2-λ-1の2つの解をα,βとすると,
f:id:shironetsu:20151114184747p:plain

リュカ数列に対しても同様の等式が得られる.

命題1.8
f:id:shironetsu:20151114184756p:plain

 今でこそ当たり前の公式だが, やはり整数しか出てこない数列の一般項を無理数のべき乗の線形和で表せてしまうというのは驚きである. これが整数になることはもちろん帰納法からすぐに示せるが, 有理数であることなら「直接」示せるだろう. αとβのべき乗を足すと√5が「消える」ことを言えばよいのだ. そのために二項定理を使う.

f:id:shironetsu:20151114184807p:plain

 これはカタランの公式と呼ばれるものである.

命題1.9(カタランの公式)
f:id:shironetsu:20151114184819p:plain

 リュカ数についても同様の式が成り立つが, 特に直接用いることはないため省略.

 ここまでで以後用いるフィボナッチ数についての一般的な定義と命題については一通り述べ終えた. この準備を踏まえ, 次の記事から整除性と周期について述べる.

*1:この観察はラグランジュによるものとしてコクセター『幾何学入門』(ちくま学芸文庫)で紹介されている……もっと前から知られていそうな気はするが

*2:数論の話になるため自然数に0を含めない

*3:法10を導入にしておきながら…

*4:一般の二階線形回帰数列にも同じくフランスの数学者エドゥアール・リュカ(1842-1891)の名前が冠せられていたりするがこの名前で呼ぶ

*5:それ以前にカッシーニ、カタラン達も発見していたようである.

*6:J.P.M.ビネー(仏 1786-1856). この名で呼ばれることが多いこの公式は, しかしビネーによって「再発見」されたものに過ぎないことが今では知られているようだ. 約1世紀前にオイラー, ダニエル・ベルヌーイ, ド・モアブル達が既に得ていたらしい. このような発見者と名称の不一致は歴史上しばしば見られるが(誰が最初に言い出してどうやって広まるのだろう?), ここでは一応ビネーの公式の名で呼ぶことにする.

*7:多分発見当時は行列の対角化自体をしたわけではないと思う.

2次元・3次元ラグランジュ点の位置と安定性

 私を無理矢理平面の肉体に押し込み、案内役も無しにこの宇宙を探索させた〈圧搾者〉も、精神の消失からは私を保護してくれたようだった。復活した直後の安堵と微睡みのなか視野に飛び込んできたのは星空だった。やはり2πラジアンの。しかし今度は赤い太陽は無い。ただしどうやら拠り所の無い空虚、恒星間空間の寂しさを味わうことになるわけではなさそうだ。この領域を重力によって支配している星がどれなのかはすぐに見当がつく。視野中の近い場所に並ぶ、一際明るく輝く2つのそっくりな青い星、連星だ。そして私のいる位置が特殊な場所であることもまもなく分かった。双生児に挟まれた領域の視角度はおよそπ/3ラジアン――ここは正三角形の頂点であり、そして多分連星系のラグランジュ点。

 重力の逆2乗則が成り立つ我々の3次元宇宙では、共通重心の周りを重力によって円運動する2天体の周りに、無視できるほど小さい質量の物体が2天体との相対位置を変えずに留まることのできる点が5つ存在する。2天体からの重力と遠心力が釣り合う点である。これらのことをラグランジュと呼ぶ。一般的に、質量の大きい主星Aと小さい伴星Bに対して、各点は図のようにL1~L5と呼ばれる。ABを結ぶ直線上にL3、L1、L2が並び、ABを1つの辺にもつ2つの正三角形の頂点のうち回転方向前方にL4が、後方にL5が位置している。
f:id:shironetsu:20150831002009p:plain:w400

 太陽-地球系のL1点は常に太陽を観測することができる。現在も運用中のSOHOや、最近地球の前を横切る月を撮影して話題となったDSCOVRはここにある。

 太陽-地球系のL2点は地球の陰によって太陽光が遮られるため宇宙望遠鏡の設置に誂え向きである。実際、宇宙背景放射を観測したWMAPやプランクはここに置かれた。ハッブル宇宙望遠鏡の後継機であるジェームズウェッブ宇宙望遠鏡*1もここに置かれる予定になっているらしい。

 L3点は主星を挟んで伴星の反対側に位置するため、太陽-地球系、地球-月系のどちらでも特に利用されていない。

 太陽-木星系のL4、L5点にはトロヤ群と呼ばれる小惑星の集合がある。なお、より一般には太陽と各惑星の系のL4、L5点に位置する小惑星の集合をトロヤ群と呼び、土星や地球にも存在するらしい。*2

 ところで、L1~L3には天然の物体は存在しないとされている。というのも、ここでは物体は軌道修正無しには安定しないためである*3。一方、L4、L5点にはトロヤ群のような小惑星や塵が存在するとされている。ただしその安定性には主星の質量が伴星の(25+3√69)/2≒24.96より大きいという条件が課せられる。たとえば太陽-地球、太陽-木星、地球-月系ではそれぞれ33万、1000、81倍と条件が満たされる*4。ところが、最近探査機ニューホライズンズの撮影した初めての詳細な地表の画像で話題になった冥王星-カロン系では、冥王星の質量がカロンの9倍程度であるため*5、この系の正三角形ラグランジュ点で軌道は安定しない。

 さて、ここまでは3次元の話だった。では「逆1乗則」の引力が働く2次元ではどうか?実は相変わらず直線上に3つ、正三角形の頂点に2つ、ほとんど同じ位置関係で「ラグランジュ点」が存在する。そしてL1~L3が不安定なのも同じである。しかしL4、L5での振る舞いがやや異なる。なんと2天体の質量比によらず安定するのだ。

 連星系ならここに太陽風や星間塵などに由来する物質が集積してニーヴンの『インテグラル・ツリー』に出てくる「スモークリング」のような地面のないガスの濃い領域ができるかも……。何せ2次元では重力は対数ポテンシャルの「井戸」を持つから物質に脱出速度は存在しない。ずっと空の同じ位置にある2つの太陽からの光を得るべく領地争いをする植物、コリオリ力の働きを生来知るジェット噴射推進動物、L4とL5で独立に発達する生態系……色々想像すると楽しい。
 ただ問題があるとすれば自己重力で惑星になってしまうのではないかということ。それはそれで面白そうだが「ふわふわ」した感じが失われてしまう。ついでに言うと安定性も分からないし、そもそも天体の形成過程などは考え出したらかなり危うい部分が多くなってしまう……。でもSFの舞台としてはおもしろい場所だと思う。

 ルンゲ-クッタ法で計算してみると*6実際この点の周りで準周期的な軌道を描く様子を見ることができる。

f:id:shironetsu:20150831003209p:plain
左から順に主星:伴星質量比1:1、9:1、81:1の、逆1乗則重力下でのL4点のまわりでの物体の軌道。2天体間距離が1で重心は原点、yが負の側に重いほうの星がある。

比較用に主星:伴星質量比81:1(地球:月)の、逆2乗則重力下でのL4点のまわりの物体の軌道を示しておく。
f:id:shironetsu:20150831005202p:plain
ちなみにこれらの図に見える楕円は、主星:伴星質量比が大きいときには短軸を延長すると主星にぶつかる。

 いや、むしろ(25+3√69)/2という数の出てくる逆2乗則重力のほうが奇妙と言えるかもしれない。単純な問題設定の中からどうしたらこんな数が出てくるのか?せっかくなので「逆k乗」の重力としてラグランジュ点の導出から安定性まで一緒に計算してしまおう。とはいっても専らk=1,2を想定し、それぞれ2次元、3次元の場合に対応している。なお平面の問題として処理するため3次元の場合のz座標について考慮されないが、z方向はxy方向とは独立に単振動すること、つまり安定であることが知られている。

運動方程式

 xy座標上で、質量(1-q)Mの天体AをxA = (-qR,0)に、質量qMの天体BをxB = ((1-q)R,0)におく。
ただし 0 < q < 1 とする。このとき重心は原点にくる。2天体が重心の周りで円軌道をとるとすると,距離のk乗に逆比例する重力の下で、運動方程式からその角速度ωは
f:id:shironetsu:20150831003644p:plain
である。この系において、位置x=(x,y)にある質量m(Mに比べて無視できるほど小さい)の物体の運動方程式
f:id:shironetsu:20150831003706p:plain
である。ただしJはπ/2回転の行列で
f:id:shironetsu:20150831003844p:plain
である。

 mω^2で割り、時間tによる微分をωtによる微分に変え(ただし引き続きドットを使う)、X=(X,Y)=x/Rによって無次元化すると
f:id:shironetsu:20150831003855p:plain
すなわち
f:id:shironetsu:20150831003903p:plain
簡単のため
f:id:shironetsu:20150831003913p:plain

とおく。

この微分方程式の平衡点を求める。そこがラグランジュ点である。

ラグランジュ点の導出

時間変化が無いことから、
f:id:shironetsu:20150831003923p:plain
が解くべき式となる。第2式から2つの場合に分けられることが分かる。

(i) Y = 0のとき(直線解)
第1式に代入して
f:id:shironetsu:20150831003936p:plain

  • k=1のとき

f:id:shironetsu:20150831003948p:plain
f(X) = X-(1-q)/(X+q)-q/(X-1+q)とおくと、
f:id:shironetsu:20150831004007p:plain
から、X < -q, -q < X < 1-q, 1-q < X にそれぞれ1つずつ解を持つことが分かる。順にX3,X1,X2とし、それぞれに対応する点をL3,L1,L2と呼ぶことにする。なお数字の振り方は一般的なラグランジュ点に対応させた。
ちなみにqが0に近いとき、その値は
f:id:shironetsu:20150831004016p:plain
によって近似できる。

  • k = 2のとき

k = 1のときと同じ区間に同数の解がある。
f:id:shironetsu:20150831004053p:plain
(i) Y ≠ 0のとき(正三角形解)
第2式をYで割って、
f:id:shironetsu:20150831004110p:plain
(X+q)(式20)-(式7)および(式7)-(X+q-1)(式20)より、rA = rB = 1が分かる。これより、(X,Y) = (1/2-q,±√3/2)。Yが正の点をL4、Yが負の点をL5としておく。


安定性の解析

 線形近似によって平衡点の周りでの安定性を調べる。平衡点(X0,Y0)からわずかにずれた点を(X0+ξ,Y0+η)とする。

 ここで予めf(z,w) = z/(z^2+w^2)^((k+1)/2)の1次の項までのテイラー展開を求めておく。
f:id:shironetsu:20150831004138p:plain
より、
f:id:shironetsu:20150831004150p:plain
これを用いると、
f:id:shironetsu:20150831004200p:plain

(i) 直線解
f:id:shironetsu:20150831004212p:plain
γ^2 = (1-q)/rA^(k+1)+q/rB^(k+1) (γ>0)と置くと、この2元2階線型連立微分方程式固有値λは
f:id:shironetsu:20150831004227p:plain
ここで、式(10)より、
f:id:shironetsu:20150831004238p:plain
を導いておく。

  • k = 1のとき

f:id:shironetsu:20150831004251p:plain
L2では X0 < -q < 0, rA < rB より、γ^2-1 > 0
L3では X0 > 1-q > 0,rA > rB より、γ^2-1 > 0
L1での安定性を調べるために、0 < p < 1として X1 = p(1-q)+(1-p)(-q) = p-q(AとBをp:(1-p)に内分する点)とおき、式(10)を用いてqを消去すると、
f:id:shironetsu:20150831004304p:plain
ゆえに γ^2-1 > 0。なお q = 1/2 のとき X1 = 0 となるが、このとき γ^2-1 = 3 > 0 で同様。

以上より、直線解では必ず正の実数の固有値を持つため不安定。

  • k = 2のとき

k = 1のときよりやや煩雑だが必ず正の固有値を持つことが示される。

(ii) 正三角形解
(X0,Y0) = (1/2-q,±√(3)/2)
f:id:shironetsu:20150831004319p:plain
固有値λは、
f:id:shironetsu:20150831004329p:plain
から、
f:id:shironetsu:20150831004341p:plain

  • k = 1のとき

f:id:shironetsu:20150831004352p:plain
0 < q < 1で0 < q(1-q) < 1/4であり、-2 < λ^2 < -3/2, -1/2 < λ^2 < 0となり、λは必ず4つの異なる純虚数となる。従って軌道は安定。

  • k = 2のとき

f:id:shironetsu:20150831004403p:plain
0 < q < (9-√69)/18, (9+√69)/18 < q < 1 のとき、λは必ず4つの異なる純虚数となって安定。

(9-√69)/18 < q < (9+√69)/18 のとき、1-27q(1-q) < 0 より、
f:id:shironetsu:20150831004414p:plain
 このときλは純虚数でない複素数を含む。式(43)から、4つの固有値の和は0(∵解と係数の関係)であるため、その実部が全て負となることはない。ゆえに不安定渦状点になる。

 結局、2つの天体の質量比 r = q/(1-q)について、
・ 0 < r < (25-3√69)/2, (25+3√69)/2 < r のとき安定
・ (25-3√69)/2 < r < (25+3√69)/2のとき不安定となることが分かる。

 復活によって得た新しい肉体は赤い太陽を周回していたときのそれとは異なっていることに気付いた。自分の姿をはっきりと見ることはできないが、おそらく対角線の長さが2対1程度の菱形をいくらか丸くした形だ。こうして星と自分の体を眺める眼は両方の鈍角の頂点の周囲にあるように思われた。この肉体はここにいるのが自然であるような姿なのだろうか?
つまりこのラグランジュ点を利用して生きる住民の? 他の点よりも集積しやすいとはいえ、物質はやはり希薄に見える。ここだけで暮らすのはあまり魅力的には思えない。しかしもし中継基地があるならこの連星系から抜け出して旅立つには重要な場所だろう。

 星と星を渡る――定義はともかく文明を持たないという条件付きで――生物は進化の過程を考えるとほとんどありそうに思えない。時間、エネルギー、物質のどれもが不足する。だがそれはあくまでも故郷の宇宙での話だ。現に私は赤い太陽の軌道上の生物圏を見てきた。物理定数同士の関係が変わり、例えば化学が相対的により激しいエネルギーの交換を伴うものであれば、生物のあり方も全く異なるものにもなりうるだろう。

 私は2つの青い「目」に見つめ続けられながら、正三角形の頂点で力を抜いて待つことにした。

*1:予算の増大と完成時期の延期で計画は色々と難航しているらしい。広報ツイッターアカウント→ @NASAWebbTelescp https://twitter.com/NASAWebbTelescp

*2:L4よりL5のほうが小惑星が多いという謎についてクラークの『神の鉄槌』で触れている箇所があった気がする。しかし観測誤差によるものという話もあって詳しいところはよく分からない。

*3:視聴当初は気付かなかったのだが、今放送している戦姫絶唱シンフォギアGX第1話「奇跡の殺戮者」冒頭で出てきたスペースシャトルラグランジュ点で活動していたとか(ちらっと風鳴司令が口にしている)。ナスターシャ教授はどこのラグランジュ点に漂っていたのだろう。

*4: 太陽系で最大の惑星は木星であるため、他の惑星も太陽とのラグランジュ点を持つ。  また、地球 : 月は太陽系の惑星-衛星では最大の質量比を持つ。ゆえに他の7個の惑星ではいずれも条件が満たされるが、単純な制限三体問題として考えられるのか(たとえば木星には質量の近いガリレオ衛星が周回していて、他の衛星からの影響を無視できない)、安定性なのかどうか、よく分からない。

*5:冥王星が惑星から「降格」されたはまさにこの「衛星」が「惑星」に対して大きすぎるという事実が一つの理由になっている。

*6:Perlでプログラムを書いて出てきた数値をExcelにプロットさせた。gnuplotのような専用ソフトにやらせたいところだったのだが…。

燠火(解説)

 グレッグ・イーガン『白熱光』、ハル・クレメント『重力の使命』*1ロバート・L・フォワード『竜の卵』、スティーブン・バクスター『天の筏』、ラリイ・ニーヴン『インテグラル・ツリー』、そしてアーサー・C・クラーク野尻抱介の多くの短編……等々多くのSFで扱われる「重力ネタ」*2。身近でありながら少し違った環境に出ると驚くべき振る舞いをする、もっともSF的な自然法則の一つであり、特にハードSFの題材として親しまれている。

 「燠火」(おきび)はそういったSFにあこがれて書いた。「作業仮説」よりはSFとしての体をなすように目指したが、結局(特に最後の方)長々とした説明にならざるを得なくなって残念。

 舞台は2次元宇宙。相変わらず『ディアスポラ』の5次元宇宙物理やA・K・デュードニー『プラニバース』に影響を受けている。作中の世界の構築に興味深い科学的考察を使ったわけではないが、語り手がぶつぶつと考える物理は本物である。高度が下がれば加速し、上がれば減速するというのは我々の宇宙と同じ。というか衛星の挙動には、少なくともこの作中では面白い現象は起きない。ただしその推進方法が問題となる。

 実は当初、語り手が最初に予想したように「磁場中の浮力」を使って高度を変えて推進する方法について考えていた。舞台もそれをもとにして構築した。しかし閉じた定常電流は内部に均一な磁場を作るのみで、外部には磁場を作らなかったのだ!つまり環状電流同士は磁気的な相互作用を持たない。ざっくり言うと2次元空間に磁力は存在しない*3。これには困った。磁気による衛星の推進といえば藤井太洋『オービタル・クラウド』が記憶に新しいが、あの手が使えないのである。というわけで白状するとその推進機構については投げた。捕食者に食われてしまった語り手が再び目を覚ますときが来れば語ることになるだろう(?)。

 語り手が考えた問題はもう一つある。重力そのものの由来である。こちらには一応の答えを出した。自然に考えれば重力は距離に逆比例しそうだが、一般相対論から静的重力は現れない。リーマン曲率テンソルとエネルギー運動量テンソルの独立成分の数がともに6個となってしまうことが原因、と説明されるがこれについてはよく知らない。しかし3次元でニュートン万有引力の法則を再現するときと同じように線形近似で計算するとそれが分かる。作中でも述べたとおり、質量ではなく運動エネルギーが重力の源として卓越する。一般相対論由来の万有引力の法則は存在しない。ではどうすればいいか、というと最も素朴には電磁気の法則と同じものを成り立たせればいい。クーロンの法則と万有引力の法則が完全に対応するようにするのである。ただし同符号同士で引き合うように。その電荷の対応物はまさしく「重力質量」である。慣性質量と重力質量が同じものであることを言う等価原理を無視してしまう。マクスウェル方程式型の重力法則は相対論的共変性を持つため数学的には問題は無い*4……が、負のエネルギーが現れてしまうのが問題らしい。この点については正直よく分からない。

 2次元の電磁気と重力の法則についての面白そうな話題を入れた短いストーリーを書くことができて満足した。もっと考察ができれば、たとえば運動エネルギーによる重力の支配する世界とか、回転する帯電円盤内の磁場に満たされた空間を描くことができそう。2次元世界の話というのはポピュラーでありながら、ハードSF的考察の加えられたものは案外少ない。面白い話題だと思うので、そういうSFはもっと読みたいし、色々と考えたりもしてみたい。

*1:原題は"Mission of Gravity"。『重力の使命』は早川書房版で、東京創元社での邦訳タイトルは『重力への挑戦』。半世紀以上前に書かれた素朴な味わいのハードSFの名作として有名だが、数年前に創元SF文庫から一時期復刊されたのみで今は中古でしか手に入らない。どっちでもいいから復刊して!

*2:書いていて気付いたが真っ赤な霧って『天の筏』の世界の視覚的イメージそのまんまじゃないか……。パロディとして許してほしい。

*3:これは実は『プラニバース』にもさらっと書かれているのだが真意を理解していなかった。

*4:一般相対論は重力を説明するものではなくなるが消えてしまうわけではない

燠火(後半)

 同類――私と同じくプレスされた人間?この世界についての有益な情報を得るためか、単に共感を求めるためか、危険の可能性を省みることなど忘れた私は一心不乱にコンタクトの手段をこの体の中から見つけるべく、再度まだ使っていない感覚を探った。まずはこちらが名乗って挨拶し、名前を教えもらおう。それまで彼のことは1号と呼ぶ。

 ここが真空であることには今やほとんど疑いを持っていない。音波はだめだ。では電磁波は?体を縁取る線に並ぶ点の中にクラゲの如く発光するものがあり、それが意思疎通に用いられるものなら受容器も備わっていなくてはならない。そしてその「言語」を解する能力が――たとえ人間の言語並かそれ以上に複雑なものであろうと――〈圧搾者〉は私の脳に必ずそれをインストールしているはずだ。自分の目指す姿に玩具の光るコマを想像した私はその滑稽さを感じながらも、あるともないとも知れない感覚を操ろうと必死になった。しかし一向に応答は見られない。1号は無視を決め込んでいるのか、眠っているのか、あるいは死んでいるのか。ここでふと彼がずっと動かないでいてくれるなら重力の「逆比例則」を検証することができることに気付いた。
CW解は4つの定数A、θ0、c0、c1によって
f:id:shironetsu:20150817185224p:plain
とも表せるがこちらのほうが意味は明確だ。「太陽」からの平均軌道半径は、目標物体のそれをRとするとR+2c0で表せる。近ければ角度方向正の向きへ離れてゆき、遠ければ負の向きへ離れていく。しかし平均半径が等しければ円軌道をとる目標物体の周りを、公転と反対の向きへ縦横比1 : 2αの楕円状に周期2π/αωで回転することになる。公転周期はT=2π/ωであるから、その1/α倍だ。今もしk=2なら1/√2倍――すなわち1回の公転の間に目標物体の周りを約1.4回巡ることになる。目印は必ずしも円軌道でなくていい。平均距離さえ等しければお互いの位置関係は周期T/αで元に戻る。

 これを利用すれば1号を利用してαを測れる。まず彼に追いつき、「太陽」からの平均距離の等しい軌道に乗る。そして公転周期を測るため空に目印となる星を見つけ、1回公転する間の彼との上下関係の変化回数を測る。1回半弱であればその計測はk=2を支持する結果となる。

 問題は平均距離の等しい軌道に乗る方法、公転周期を同期させることだ。公転軌道の扁平率が1に近いほどそのずれの測定に及ぼす影響は大きくなる。動的な修正によって公転周期を彼と同期させることは可能だろうか?

 そんなふうに逆比例則の検証方法について考えてぼんやりとしていると、1号が突然動き出す気配を感じた――波だ。彼が眠りから覚め空間に波を立てたのだ。それは人間の器官だけを持っていたときには知らなかった感覚だった。存在に気付き襲い掛かってくることを警戒した私はすぐに動き出せるように身構えた。心臓があれば速くなった鼓動が聞こえるところだろう。

 緊張の瞬間はすぐに彼によって破られた。彼はゆっくりと斜め下方に向かって動き出した。私をおそれて逃げている?しかしここで見失ってしまうわけにはいかない。私は彼を追跡することにした。

 1号は私が経験したのより明らかに速い速度で降下していった。新参者にとってはなかなかのスリルを伴う体験だ。体感で――ここでは主観的時間以外の基準に意味はない――数分間の降下の後、星が瞬き始めていることに気付いた。太陽の最上部層を潜ったのだ。これが故郷の太陽; ソルで私が通常の有機物の体を持っていれば、とっくに焼け焦げるどころかプラズマと化して跡形もなくなっているところだ。だがこの体にはその心配は無さそうだった。私の懸念は別のところにあった。大気のある層に潜れば必然的に抵抗による減速を受ける。自ら体勢を立て直すことができなくなるほど減速してしまえばあとは落下するのみだ。この体でそんな事態を回避することができる高度を気付かないうちに下回ってしまうことはありえる。1号にこのままついていくことに不安を感じ始めたとき、ようやく彼が減速するのを感じた。

 そこには赤い霞がかかっていた。人間の色彩感覚の残滓を備えていた私には熱さを感じさせる風景だったが、実際には熱は感じていなかった。そもそも温感が無い可能性も無視できないが。その代わりに私が感じたのは原始的な衝動だった。1号がこの層に降りてきた理由もおそらく同じだ。このクラゲの肉体は栄養源としてこの赤い霞を摂取するようにできているのだ。この4回対称の体に口があるとすれば腕と腕の間か、そこにあるとなら消化器も袋状のものが4つずつか、分析的な目で理性を保てるふりをしながら私は特に味覚を刺激しない赤い霞が空腹感を満たしていくのを感じた。これではまるで、その世界の食べ物を口にした瞬間に元の世界との繋がりを断たれ、そこに囚われることになる古い物語の類型をおそれる子供のようではないか。自分がこのクラゲになりきってしまうことへの抵抗を忘れないようにしながら、私は霞の濃い方向へと絶え間なく移動を続けた。

 この層にやってきた目的を思い出したのは空腹感が満たされ、理性が主導権を取り戻した後だった。いつの間にか1号は視界から消え、周囲には一段と濃くなった赤い霞が拡がっていた。食後の倦怠感に包まれた私は、そこで追跡が失敗に終わったとみなすことに決めた。ここに食料があるなら同類にはまた出会えるだろう。

 訳も分からず口にした赤い霞で満ちたこの空間を改めて見回すと、今更ながらに気味の悪さを感じることになった。赤は何よりもまず血を思い出させる。アンタレスがサソリの心臓であるように、この太陽も心臓ならさしずめ私は臓器を蝕む寄生虫だ。

 ここで私は疑問を抱いた。アンタレスの表面にもこんな霞はかかっているのだろうか。多少太陽に近づいたとは言え、この層は明らかにまだまだ希薄な領域だ。「表面温度」を測るには多分まだ冷たすぎる。

 減速の恐怖から高度を維持しようと努めて漂っていると、しばらくの後にヒントは現れた。赤い綿だ。私の5分の1程度の大きさの扁平な楕円、ちょうど霞を固めて雲にしたような印象を抱かせる。それと同時に私が想像したのは植物プランクトンの群体だった。小さな粒子状の単細胞の藻が集まって形成する群体――赤色光が光合成色素の色なら。となると〈圧搾者〉が私に赤色を見せた理由についての仮説は一旦捨てねばならなくなる。この赤色が星の表面温度では無く植物の発する色であるなら、緑色を見せていればもっと早くに気付けたはずだ。

 プランクトンを食べる私のこの体はもしかするとクジラなのかもしれない。それとももう少し単純な軟体動物か。いずれにせよこの赤い藻類より複雑な体制を持っているのは確かだ。であれば植物の方ももっと高等なものがより下層、大気のいくらか濃い領域に繁茂しているに違いない。そして太陽の本当の表面も見ることができるだろう。私は再び降下を開始した。

 未開の地に足を踏み入れた博物学者の気分で「森」を目指す探検は、しかし「波」によって中断を余儀なくされることとなった。何かが迫ってくる。さっきと同じ同類ではない。本能はサイレンを鳴らし、私の感覚を鋭敏にさせた。本能的な恐怖――それは天敵、捕食者の接近を意味する。予想しておくべきことだった。ここがよほど原始的な生態系の発展途上にない限り、草食動物がいるなら当然肉食動物もいる。そして私のこの体は草食動物、被食者のもの。捕食者に狙われていることは警戒してしかるべきだ。

 公転方向前方下側! 「波」を捉える感覚が捕食者の接近方向を察知させた。姿の見えない捕食者は、その高度を上げれば角速度を減じて私とランデブーすることが出来る。そうなれば私の負け。しかし同時に私も高度を上げれば公転後方に後退して逃げることができる。

 それ以上軌道力学を検討する暇など無い。すぐに宇宙空間への上昇に転じた。逃げ切れることを祈ろう。
結局私はこの一心不乱の鬼ごっこに勝利することができた。おそらく最初にいたところと同じくらいの高度に戻ってきたのだろう。星は再び瞬きのない点となり、視界の半分以上が暗黒の天球に占められた。

 ジャングル探検をするにはまだ私は未熟すぎた。もう少しここで星を眺めていよう。星を睨みながら、公転前方と後方で星の色に変化は見られないか観察した。もし公転速度が光速に比べられる程度に大きければ前方は短波長側に、後方は長波長側にドップラー偏移を起こすだろう。しばらく待ったが変化は見られなかった。無論暖色灯のもとでも色の認識が変わらないように、わずかな波長の変化を脳が処理して補正している可能性もあるが、さしあたってここにある物体が亜光速に達しているとみなす必要はないだろう。

 この観察は相対論を考えるための第一歩だ。私の関心は太陽の照らすもとに縛り付けている重力だった。
重力を説明するのはアインシュタイン方程式だ。
f:id:shironetsu:20150817185236p:plain
定数κは故郷の3+1次元宇宙でなら8πG/c^4と表されるだろう。しかし万有引力定数はいわば派生的な定数だ。ちょうど真空の透磁率とクーロン定数の関係のように。n+1次元時空ではκを定数にとるほうがいい。

 そう、弱場近似の仮定の下でアインシュタイン方程式に線形近似を施すと一般相対論の方程式は電磁気の方程式との類似で語りうるものになるのだ。ミンコフスキー時空からのずれが微小であるとするとは
f:id:shironetsu:20150817185243p:plain
ととることである。エネルギー運動量テンソルとして完全流体のそれを取ると、「ローレンス条件」のもとでダランベール演算子□を用いて計量の方程式として
f:id:shironetsu:20150817185252p:plain
が得られる。ρは各点の質量密度、vは質量の流れの速度で、完全流体の式に従う。そこでの質量m、速度Vの質点の運動方程式は測地線の方程式の近似として
f:id:shironetsu:20150817185305p:plain
が得られる。実際には軌道はmに依存しない。ここで場に時間変化は無いという仮定を加えて次のように置き換えると電磁気との類似は明らかになる;
f:id:shironetsu:20150817185312p:plain
いわば「重力の電磁ポテンシャル」と「重力の電磁場」、そして「重力のローレンツ力」だ。ここから「重力のスカラーポテンシャル」についてポアソン方程式を解くとニュートン万有引力の法則が導かれる――が、それは3+1次元以上の宇宙でしか通用しない。n=2では明らかに事情が異なる。式を見れば分かるように、(n-2)の因子が入っていることによって、静止した質量は「重力電場」を作り出せないのだ。「重力電場」の源としては質量の流れの速度の2乗に比例する項が支配的になる。
f:id:shironetsu:20150817185318p:plain
 エネルギーと質量の等価性を考えれば自然なことだ。2+1次元宇宙では運動エネルギーが重力場を作る。

 違いはそれだけではない。静止質量が卓越しなくなることで、「重力磁場」も同程度の働きを担うことになる。普通の磁場は平行電流同士で引き合い、反平行なら反発する。アンペールの法則だ。しかし重力磁場は平行な質量の流れ同士で反発し、反平行で引き合う。ちょうど万有引力クーロン力で符号が逆転するように。

 私はここまで考えて赤い太陽を眺めた。ここに来たときと違い、今ではあの赤さが光合成色素によるものだと考えている。2+1次元宇宙版ヘルツシュプリング・ラッセル図にあてはめるのはしばらくやめておこう。もしかすると天球に分布する星々のあの色も星本来の色ではないのかもしれない。

 そんな奇妙な重力のもとで、あの太陽はまとまっていられるのだろうか?核融合を進行させる圧力は?しかし事実としてこの太陽は存在するのだ。そこを回る私や1号、捕食者も。

 重力を電磁場のアナロジーで語る前に、2次元電磁場についての理解を深めておかねば。視覚がある以上――偏執的に脳による再構築を疑わない限り――ほとんど間違いなく電磁波は存在する。

 では電磁場はどのような力を及ぼすか?基本はローレンツ力だが、場にかかる応力を電気力線のように捉えるためには先にも用いたマクスウェル応力テンソルを考えるのが便利だ。
行列形式で書くと、
f:id:shironetsu:20150817185325p:plain
となる。ただし2行目でBはB12成分を表す。

 2次元において磁場の独立な成分はB唯一つしか存在しない。これを擬スカラーに対応――3次元において2階反対称テンソルの磁場を擬ベクトルに対応させるのと同じく空間の向き付けに依存――させると、それは「濃さ」としてイメージできるだろう。中性が灰色なら、正は白色に、負は黒色とすればよい。同符号なら強め合って色は濃くなり、異符号なら弱めあって灰色に近づく。平面内に磁場のグラデーションができることになる。便利な見方だ。

 注意を引いたのは磁場の応力テンソルだった。この形は等方的な応力、圧力のそれだ。符号を考えるとBの「濃い」ところほど―白黒どちらでも――圧力が強くなることを意味する。

 磁場のおよぼす力は3次元では磁力線を持つ弾性体のアナロジーで捉えられてきた。ここ2次元では弾性体が完全流体に置き換わるということだ。圧力B^2/2μ0の流体――磁場を作るものは磁場の中を「浮き沈み」する。

 そのイメージが閃きをもたらした。磁場の白黒は電流の左回り・右回りに対応するだろう。ならば磁場に関して環状電流はモノポールのように振舞うはずだ。そしてその相互作用は圧力の考え方を使うと、まさしく異符号――右回りと左回り――同士で引き合い、同符号同士で反発することになる!というのも、異符号同士だとその間にできる磁場は弱められ、圧力も小さくなり、「浮力」が働くからだ。同符号同士ではその逆のことが起こる。

 気分転換に太陽の方へと近づいてみたり、遠ざかったりしてみた。私はこれを飛行として捉えていた。しかしもし磁場の中を浮き沈みしているのなら泳ぎにたとえるほうが近いかもしれない。太陽が磁場をもっているなら、体の中の電流の向きを変えることで浮き沈みすることが可能になる。

 仮説を立てたなら計算しよう。円電流同士の相互作用を求めるため、今度は電磁ポテンシャルで電磁場の方程式から始める。
f:id:shironetsu:20150817185331p:plain
時間変化が無いと仮定して磁場についてのみ考える。原点から遠方で十分早く減衰する境界条件の下でポアソン方程式を解くと
f:id:shironetsu:20150817185338p:plain
から多重極展開が得られる。「点磁荷」が存在しないことや多重極子の対称性から、高次の項を無視すると
f:id:shironetsu:20150817185345p:plain
Mが磁気双極子モーメント、Tが磁気四極子モーメントを意味する。
双極子モーメントの作る磁場を考えよう。B12成分のみを考えればよいから、
f:id:shironetsu:20150817185351p:plain
――不意をつかれた。ゼロだ。磁気双極子モーメントは2次元では磁場を作らないということだ。

 3次元においては磁石の源は(少なくとも古典論的には)微小な環電流の作る磁気双極子モーメントとして説明される。いわば磁場の源の最小単位だ。私がこの世界に来る前に高次元の磁場の源として計算したのも磁気双極子で、確かに磁場は形成されていた。その磁気双極子が、2次元でも存在するのに磁場に寄与しないのだ。

 では四極子を計算しよう。四極子は各成分が先の関係式によって結ばれているため、独立成分は2つしか存在せず、
f:id:shironetsu:20150817185358p:plain
と表せる。これならどうか?さっきよりいくらか複雑な計算の結果はまたもゼロだった。この時点で私は悟った。おそらく2次元空間で閉じた定常電流は遠方に磁場を作らない。

 私は悩んだ。ヒントは素朴にアンペールの法則を思い出すことで得られた。磁場と電流の関係は
f:id:shironetsu:20150817185405p:plain
となるのだった。これはナブラを使った表記では
f:id:shironetsu:20150817185411p:plain
を意味する。この式の意味するところは、「擬スカラー」Bの勾配が電流をπ/2回転した方向を向くということだ!となるとそのイメージは簡単だ。円形に正の向きに電流が流れる回路を考えよう。電流をπ/2回転させると全て中心を向く。「デルタ関数的」な勾配を考えると、そこにできる磁場は回路に囲まれ3次元方向に張り出した円柱として想像できるだろう。内側、外側でそれぞれ一定値をとることになるのだ。そして自然に考えると外側でその値はゼロになる。私はここでソレノイドコイルのつくる磁場を思い出した。無限に長いソレノイドコイルは内外それぞれに均一な磁場を作る。ここ2次元世界では均一な磁場はもっとありふれたものになる。先の磁場の濃さの考え方を使うと、回路内部で均一な白色または黒色で、外部は中性の灰色になる。これをグラデーションとは呼べない。ある意味で回路がモノポールになるという見方は正しかった。しかしそれは内部にしか磁場をもたないものだった。

 ここまでくればマクスウェル方程式全体をより直感的な形に書き換えてしまおう。π/2回転の行列をKで表すと、
f:id:shironetsu:20150817185422p:plain
が電磁場の方程式となり、ローレンツ力は
f:id:shironetsu:20150817185428p:plain
になる。

 私は途方に暮れた。磁気双極子同士の相互作用として太陽への浮沈を説明できると考えていたが、それは間違いだった。太陽と私以外を真空と考える以上、磁場は力を与えない。

 では何を仮定すればいいのだろう。実はここを真空をみなすのは間違っている?あるいは輻射圧はこの体を浮き沈みさせる程度に大きい?それとも全く別の力?

 もっと考えるための材料が要る。私はここに来てもっとも単純な可能性としてのいくつかの物理法則を仮定して考えることを試しているに過ぎない。

 何から手をつけようかと思案しているときだった。捕食者はこの高度にもいたのだ――後方上側から「波」がくるのを感じた。下側に逃げるほか無い。2度目の逃走の中、私は重力について考えた。この重力が計量によるものであることを支持する現象は未だ見ていない。しかし星の様子はあまりにも太陽に似ている。場が2つのマクスウェル方程式を持つことを禁じる法則はあるのだろうか?もし無ければ、重力質量とそれについてのマクスウェル方程式を考えてみるほうが、運動エネルギーが重力源となると考えるよりも、私を太陽の方向へ引っ張り太陽自身をまとめる力の可能性としては単純かもしれない。

 捕食者から逃げる中でそんなことを考える余裕を持つためには、彼らについて私は無知すぎた。私の目が捉えたのは、先の割れた楔形をした赤い斑――というか破線――の物体だった。表面の艶は甲虫のそれを思わせる。彼らは狼のように群れを成して狩りをするための連携をとれる賢さがあったのだ。まんまと罠に追い込まれた被食者は、その裂けた楔形の顎に捕らえられ、身を潜めていた他の仲間が3,4と集まってくるのをもがきながら眺めた。

 やるべきことはまだある、と言うよりまだ何もやっていない。〈圧搾者〉の目的が何であれ、わざわざ呼び寄せた以上むざむざとここで殺してしまいはしないだろう。私は意識が途切れ途切れになってゆく中で次に期待した。



================


(この短いSF(?)中で語り手が一人考えることは、2+1次元宇宙の電磁場についての考察に基づいている。2次元、「フラットランド」の電磁場である。作中いきなり出てくる2次元マクスウェル方程式の妥当性については、一般次元(n+1次元時空)電磁場についての考察そのものを扱った「作業仮説」
「作業仮説」というSFもどきを書いた - Shironetsu Blog
を読んでほしい。)

(内容に疑わしい点も多い。指摘、感想等あればコメントでください。)

燠火(前半)

 目覚めた私が捉えたのは2πラジアンの視野だった――4πステラジアンではなく。そのうち半分よりいくらか少ない領域はいくらか眩しさを感じる赤い光の線が占め、もう半分にはわずかに赤みがかったり青みがかったりしている弱く白い光点が所々に混ざった暗い線が見える。注意すれば、赤い光の線と黒い線のぼやけた2つの境界の一方から小さい光点が現れ、もう一方には同じ速さで消えていくのが分かる。私はそのゆっくりとした流れにオルゴールのホイールを想起しながら、夢の続きの場面のようなこの光景を理解することを始めた。


 夢の中で展開されるのと同じ、目に見えたものを都合よく理解しようとする論理が働き、私はそれを自分が2次元の世界にいることの自然な結果だと理由付けして納得した。あるいはまだ夢の中なのかもしれない。しかし、意識がはっきりとしてくるにつれ、視界に入るものの異常さを認識するだけの常識と自分の置かれた状況を検討するための材料を記憶の中から集める能力を取り戻していった。夢なら視界が明瞭すぎる。目に見えるものについて判断する前に理解が訪れる夢の感覚は無かった。当然陥るべき狂乱状態を経過することなく自らの身に起こったことの分析を始めた冷静さ、あるいは鈍さにはむしろ混乱させられた。だがすぐにその理由に思い至った。この世界、2+1次元宇宙へと私は来るべくして来たのだ。

 眠りに着く前の記憶は鮮明だった。私は一般次元の磁気双極子同士の相互作用の計算に悪戦苦闘していたのだ。

 3+1次元において、磁気双極子は軸性ベクトル量と対応させれば電気双極子が作る電場と同じ形の磁場を作る。これに加え、電場の応力テンソルと磁場の応力テンソルも同じ形になることから、磁気双極子同士に働く引力・斥力は電気双極子同士のそれと同じ形に表せることが分かる。すなわち磁気双極子は、その相互作用について、絶対値が等しく異符号の"磁荷"が両端に付いたダンベルと等価になるのだ――ただし"磁荷"間の引力がクーロンの法則に従うとして。

 より一般のn+1次元ではそううまくはいかない。磁場は2階の反対称テンソルとしてn(n-1)/2個の独立成分を持ち、磁気双極子も同じく2階の反対称テンソルと同じように振舞う。3+1次元でうまくいくのは、空間の体積形式を介した双対によって2階反対称テンソルとベクトルとの間に一対一の対応関係――ただし空間の向き付けに依存する――が存在するからだ。これは平面とその垂線との一対一関係があることと同じ事情によると言える。

 私が試みたのは、もっとも単純な状況として、2次元平面の定常的な円電流の作る磁気双極子モーメント同士の相互作用を計算することだった。それには2つの間に挟まれたn-1次元超平面の各点の応力を足し合わせればよい。ちょうど3次元と同じように。ところが手順は正しいに違いないのにガンマ関数やベータ関数が次々と現れる積分計算を経て得た答えはnに3を代入したとき3次元の場合と整合しなかった。計算過程を点検しても間違いは見つけられなかった……というか乱雑に書き散らしたためにその作業は意味を成さなかった。

 心が折れた私は布団を被った。人間が一日に取り組むことのできる積分の回数はその体力による。目が覚めて気力が回復していたら自分に一からやり直しさせよう。

だがそもそも「向かい合わせ」の双極子間の引力、斥力を計算したところでその相互作用の定性的な理解には遠く及ばないのだ。一般次元の「磁石」間に働く力を理解する術を探しながら私は眠りについた。

ところでこの過程では意図的に2次元の場合が除外されていた。「一般次元」と称しながら、である。ポテンシャルが対数になることや成分が少なくなることからくる例外的な扱いを後回しにしていたのだ。

そして目覚めたのがこの2次元空間だった。誰が私をぺちゃんこにしたのだろう?多分私は2次元への軽視によって2次元世界を尊ぶ何者かの怒りを買ったのだ。その〈圧搾者〉はこのフラットランドへ閉じ込めることで私を改心させることを目指しているのだろう。

私はいくらか催眠効果を持っていそうな光点のゆったりとした流れを眺めながらそこまでの推論を進めた。自分が嫌っているはずの、根拠の薄い悟りを土台にした論理展開に自嘲的な気分になったが、突然2次元に押し潰されて無事でいることに比べれば遥かに許しやすい荒唐無稽さだった。

 そろそろこの世界と向き合わねばならない。〈圧搾者〉はどういう技術によってか、私の精神に破壊的なダメージを与えることなく2次元のこの肉体に転写していた。削ぎ落とされたであろう無数の器官と同時に省略された神経の不在による異常は今のところ感じられなかった。残された神経はいくらか簡略化されているとは言えそれを平面に転写することは可能だろうか?ひょっとするとそんなことはせず精神は3次元に取り残されているのかもしれない。だとすればこの体も、2次元空間内に埋め込んで全てを完結させるという手の込んだパズルなどせずに、「外部」の空間に飛び出した部品をつぎはぎして集めた塊にすぎないのかも――あるいはもっと単純に全てが表面だけの見せかけでありさえするかも――しれない。しかし私はそうは信じられなかった。これは「現実」だ――あらゆる細部にわたって自己完結した法則に基づいて動いているという意味で。この世界に対して、破綻をきたすことなく整合性を保証するためにあてがわれるもっとも単純な方法としての数学的仕組みとして物理法則が働いていると考えるのは決して間の抜けた妄想ではないだろう。その確信は〈圧搾者〉が力ずくの方法として脳を薬漬けにしたことによるものでもありえた。だがもしそうなら、2次元を尊重させるために嘘の舞台と心理への直接の干渉のどっちつかずの方法を取ることにどれだけ意味があるか?いずれにせよ確信とそれを支える感情は自己肯定的だった。何一つとして信用に足りる証拠を自分に示すことはできなかった。しかしこれ以上の疑いで、この世界を解釈するためのもっとも簡単な道具を自分から奪うつもりにもなれなかった。

 2次元への単なる埋め込み以上におそろしささえ感じる巧妙な技術は、この2次元の肉体への神経の調整と適合にあった。視覚は2πラジアンの視野を、もともとそうであったかのように自然に処理していた。ただし3次元的視野も忘れることなく。

 少し力を入れると4本の「腕」が動くのが視野に入った。それは気味が悪く驚きをもたらす体験であったには違いないが、半ば諦めのようなものを感じながらすぐに受け入れる他なかった。付いていたのは4本ともほとんど同じ形状、大きさの、わずかに透き通った肌色を帯びる鰭のようなものだった。いや、鰭と言えるかどうかは分からない。少なくとも、1次元的情報から奥行きを捉える感覚――3次元人の立体視の対応物――はその一本線を、鋭角二等辺三角形の辺を膨らませた形に認識した。
この鰭の生える部位に眼は付いているはずが無いから、人間が2つの目からの情報を統合するのと同様の方法で隙間は埋められているはずだ。五感の他の感覚由来の情報は、眼の周囲にこれらの感覚器が備わっているとしても今は感じられなかった。鰭はおそらく胴と同じかそれ以上の長さに渡って延びていたが、力を入れても互いに触れさせることはならず、触覚についても分からない。自分の体のうち見ることができるのはこの4つの鰭だけだった。

 ここまで徹底的な神経の再配線が施されていることが明らかであれば、自らの連続性についての確信は揺らいで当然のはずだった。脳をかき混ぜられたと信じているのに自分の過去への信用を喪失しないのは矛盾している。だがまたしても混乱は訪れず、新しい肉体への理解は自分を再発見する作業であるかのように淡々と進んでいった。

 私はここまでで知ったことから自分の姿として4本の鰭の生えた円盤を想像した。鰭を不器用にじたばたと左右に反らす憐れな円盤だ。ちらっとミズクラゲを連想したが、すぐに思いなおした。
この2次元世界においてこの形態はもっとも厚みのある姿だ。「平ら」に思われるのは、存在しない"z軸"方向から3次元的な視点で想像しているからに過ぎない。その仮想の"z軸"、いわば真横の方向から見れば内臓まで全て見えるだろう。ただし管状構造は存在しないはずだ。そういう意味ではやはりクラゲに近いのかもしれない。それ以前に、消化器の存在や代謝の仕組みの有無など、この肉体がこの世界で生存可能なつくりになっていることへも疑いを向けなければならない。何を食べればいいというのだ?

 眼は少なくともその鰭同士の間に4つは配置されているはずだがもっと多いかもしれない。その眼から延びる神経のつながりについては何も分からない。交叉は最小限にとどめられているには違いないが全く無いというわけにはいかないだろう。直線的な電気的信号伝達とは異なる、例えば複数種の化学物質の濃度勾配を積極的に利用するような巧妙な仕組みが存在しているのかもしれない。

 その眼によって再び周囲に注意を向けると、光点は変わりなく流れ続けていた。〈圧搾者〉はこの憐れなクラゲに、親切心からか、3次元の類似物への連想につながりやすい色覚を備えさせていた。視界の半分に広がるのは「天球」、暗闇を背景とした光点はその一つ一つが星なのだ。その色が恒星の表面温度に対応するものなのか、可視光のスペクトルと色の関係までは分からない。しかし少なくともそのどれもが見慣れた星の色だった。暗闇が黒くある必然性さえ無いのだから、この類似は恣意的なものであったが、周囲の理解の簡単化には大きく貢献していた。

 ではもう半分に広がる赤い光は何か?間違いない――太陽だ。この赤さが「本物」かどうかは分からない。しかし〈圧搾者〉がもしこの色付けに意味を与えているならきっとヘルツシュプルング・ラッセル図の右側、天球に輝く白い星より穏やかに「燃えて」いることを示しているのだろう。この炭火のような光に絶え間なく片面を焼かれて無事でいられる肉体は一体何でできているのか、考えを進めようとしたところで諦めた。全てに対応関係が付くわけではないだろう。

 つまり、今の私は太陽の軌道上を衛星、いや、「惑星」として自由落下しているのだ。太陽の見え方にほとんど変わりは無い。角運動量の保存則に従って円軌道を描き続ける無力な質点だ。

私はここで電磁波を自明視していることに気付いた。2+1次元時空での電磁場を公平に扱うことこそプレスされた原因なのだった。手始めに一般次元のマクスウェル方程式を思い描いた――
f:id:shironetsu:20150628133839p:plain
順にガウスの法則、磁荷の非存在、ファラデーの電磁誘導の法則、アンペール-マクスウェルの法則と、3+1次元で呼ばれるものに相当する。独立な式の数は1+nC3+nC2+n個だ。

2+1次元を例外的なものにするのはまず磁場の成分の数だといえる。n+1次元においてその独立成分の数がn(n-1)/2個になるため、2次元ではたった1成分のみ、すなわち磁場は"擬"スカラーになる。1,2成分をBで表せば,マクスウェル方程式
f:id:shironetsu:20150628133922p:plain
とたった4つの式からなることになる。磁荷の非存在に相当する式はなくなる。

再び一般次元を考えると、簡単な式変形から、波動方程式
f:id:shironetsu:20150628133945p:plain
を得る.

――私は赤い光を浴びながら、戯れにシュテファン-ボルツマンの法則を導出することを試みた。――

 エネルギー密度、ポインティングベクトル、マクスウェル応力テンソルはそれぞれ
f:id:shironetsu:20150628134006p:plain
であり、電磁場が持つ等方的な圧力Pは、応力テンソルのトレースから
f:id:shironetsu:20150628134016p:plain
となるべきものである。――ここでn=3においては直ちにP=u/3が得られること、n=2では電場の、n=4では磁場の応力テンソルのトレースが0となることに気付いた。これが3次元の特殊性の1つと
言えるものであるのは確かだが3次元を有利にしているものでありえるだろうか?――電磁ポテンシャルを介して自由電磁場の波動方程式を解けば、エネルギー密度のうち電場と磁場の担うものが相等しくなることが3次元と同様に分かる。ここから、電磁波の状態方程式として
f:id:shironetsu:20150628134058p:plain
を得る。後は熱力学の問題だ。uが温度Tのみの関数であるという仮定の下で
f:id:shironetsu:20150628134136p:plain
と、エネルギーが温度の(n+1)乗に比例することが分かる。――2次元なら温度の3乗だ。しかしそれを知ったところでこの太陽について分かることは何もない。放射強度、比例定数、半径、必要な数値はどれもこの体では到底得られそうにない。もっとも、得る必要もないからこれ以上考えるのはやめた。

 この結果を導いて満足した段になって私は数式の変形がとんとんと進んだことにようやく気付いた。紙と鉛筆と手を奪ったことへのささやかな補償として強化が施されたのだろうか?それよりは彼らが計算によってこの世界を理解することを重視したために憐れな脳を改造したという見方のほうが、もっともらしかった。困惑することに、私を一番動揺させたのは注意を向けるまで気付かずにいたこの内部の変化だった。眠りにつく前の自分との連続性への確信が揺るぐほど私は哲学に浸るつもりはなかったが、〈圧搾者〉が頭の中で複雑に入り組んでいるに違いない機能とその周辺に傷を付けずにこの変化をもたらしたと能天気に考えるのは難しいことだった。とは言え、配線の方法への不信を上回る魅力的な贈り物であることは認めざるを得ず、疑いは単に自分が正気を保って屈服しきったわけではないことを確認するための防衛的な心理であったのかもしれない。

 どこまでできるか試すべく、計算の対象を探そうと再び周囲を眺めた私の眼が捉えたのはこの世界に来て初めて認めた影だった。「下」側――上下を意識することに近い感覚はあったが、きっとこの状況では胃があれば吐き気を引き起こしたに違いない人間の平衡感覚とはまた違ったものだった――、太陽の上の黒い点だ。多分気付いていなかっただけでしばらく前から見える位置にあったのだろう。私は初めそれを目印に太陽の自転を測れることを期待したが、また新たに気付いた遠近の感覚がその距離がもっと近いものであることを告げた――あれは私と同じ「惑星」だ。

 正体不明の物体に近づくことへの恐れを脇に追いやった私は、その方向へ近付いていく方法を考えはじめ、この鰭は役に立つものだろうかと思案しつつ眺めた。既に気付いていたことだがこの高度に大気はない。あったとしてもこの鰭に抵抗を与えられないほど希薄だ。代替案として、この鰭がもし太陽帆なら光圧から運動量を得られることを思いついた私はどの角度なら望んだ方向に進めるか鰭を動かして頭を悩ませる一方で目標を睨んだ。加速に気付いたのはそのすぐ後だった。向かっているのは正しく目標の方向だ。鰭が太陽帆としての役に立った?違う、目標方向の反対側に引かれるのを感じながら私は思いなおした。それなら今まで鰭を動かしたときにも加速は感じたはずだ。別の何かが体内で働いていた。

 自転車を見たこともないのにサドルにまたがってすぐに運転できたかのような感覚を味わいながら、半ば無意識的な操舵で体が動くのに任せて私は軌道上の物体のランデブーをモデル化したヒルの方程式を思い出した。――
f:id:shironetsu:20150628134349p:plain
原点は円運動する目標物体、xは動径方向にとった軸、yは角度方向の軸を意味する。ωは目標物体の角速度だ。x方向に進めば遠心力がより引き離す方向へ働き、速度を出せば進行方向右手にコリオリ力が加わる。今見えている惑星を目指してまっすぐ進めば、余分な角運動量がその手前方向に軌道を逸らすことになるのだ。しかし軌道力学を知っていようと必ず犯すであろう失敗を回避する本能が
この体には組み込まれていたようだった。目標への接近は着々と進む。自ら身につけた技量ではないとはいえ、全自動のプログラムに任せきるよりはいくらか操縦の感覚があるこの体験を楽しんだ。

 kが4未満ならα=√(4-k)として、力fがかからない状況で解として
f:id:shironetsu:20150628135024p:plain
を持つ。k=3にz方向の単振動を加えればCW解として知られる解になる。

 問題はkだった。kは中心力;重力の距離との依存性を表す量だ。kが4を超えればわずかな摂動が指数関数的な落下か離脱を引き起こす。こうして安定的にのんきに回り続けてきた以上、kが4未満であることは確かだ。フラックスのアナロジーから、その値はn次元においてnになる――逆(n-1)乗則――のが確からしいが、根拠に欠ける。ポアソン方程式を考えることが必要だった。

ここに至って2次元の重力の特殊性を思い出したが(私に測地線を辿らせてきたこの重力は何が源だ?)、いつの間にか視直径が数度に膨れていた、ストーブの前に置かれているかのように赤く照らされた目標の姿を片側の眼で認めたことによる驚きがそれをかき消した。光沢のある線とその間に挟まれて並んだ暗い破線、所々に見えるやや明るい点――ゆっくりと回転するその輪郭を、2次元遠近感覚は中心の膨らんだ4つの刃を持つ手裏剣として翻訳した。私は同類に出会ったのだ。

(「燠火」(前編)終)

ゼンデギ (ハヤカワ文庫SF)

ゼンデギ (ハヤカワ文庫SF)

金属があることから

 先日の記事で
「これを記した筆者は重力理論の検証を兼ねた宇宙開発を目指している。しかしこの星には鉱物資源が無いため苦労することになるかもしれない。」
と書いた。
「作業仮説」スポイラー補足解説 - Shironetsu Blog

 当初はこのネタをもっと膨らませようと思っていた。
 人間の文明の成立にとって金属が必要不可欠であったことは言うまでもない。しかし、もし金属が使えない天体に知的生命が生まれたら?なかなかに絶望的であることは確か。高度な機械は作れない。医療を発達させられない。普通の方法ではコンピュータを作れない。そんな環境で宇宙進出は望めるだろうか?知性を強化させたり精神の電子化へ進む術はあるか?

 ところで、この宇宙の恒星は世代によって分けられる。大雑把には種族Iが我々の太陽を含む新しい星で、過去に存在した恒星による元素合成で重元素を含み、その惑星は地球のように(過去の超新星爆発などによって)色々な元素を持てる。種族IIは古い星で、軽い元素で構成される。*1

 種族IIの星は惑星を持ちにくいらしいが、この宇宙でも現実に重元素に乏しい天体は存在するためそういう問題は起こりうる(あるいは起こりえた)と思う。もっと身近に、カール・セーガン的ガス惑星生物*2が知性を持てたとして、彼らが自力で肉体に束縛された運命から逃れる術を手に入れられるだろうか。

 好ましい性質をもつ有機物の素材を合成する生物を品種改良で生み出すという生物工学的な手段や、自らの進化という気が遠くなる"優生学的"方法によっていくつかの問題を解決したり、長大な種族的計画を立てて金属が合成されるまで耐え忍んだり、SFとしてそういう可能性を考えるのは面白そう。

 話を戻すと、この話でそういうことに踏み込むのをあきらめたのは、重力的な振る舞いがこの宇宙と大きく異なる天体での元素合成や「第二世代以降」の天体の形成の問題を考えるための知識が全く無かったため(他の点についてはその知識が足りているのかというとそんなことはないが……)。『ディアスポラ』でも触れられていることだが多分ブラックホールだらけになる。そういうわけで金属の問題についてはあまり立ち入らないことにした。

────────────

 こういうありえた可能性について考えるにつけ、人間原理的な「後付け」であるとはいえ、この宇宙とこの惑星が好ましい性質を持っていてよかった、とほっとする。

 知性が発現するだけでは足りない。金属が無ければその先に進めない。知能はあるのに金属が足りない、という悲劇的な状況に立たされることなく済んだことにありがたみを感じねば……。

 とここまで来ると人間にも何かしら「それ以上はどうしても先に進めなくなる」というような軛が別に存在しているのではないかと未知の可能性に不安にもなる。そういった「どうしようもない困難」に気付いてしまう日は来るだろうか。

*1:「第一世代星」である種族IIIというのもあるらしい。

*2:『コスモス』などで木星の気球型生物やそれを狩る生物の可能性が考案されている。アーサー・クラークの短編「メデューサとの出会い」はこのアイデアを受けたもの。『2010年宇宙の旅』では同様の生物がモノリスに見切りをつけられてしまい、エウロパ生物の繁栄のために、第二の太陽の中にくべられることになる。