■泡細胞の幾何学(その2)

 コクセターは1つの泡に接する泡の数を

 (23+√313)/3=13.56

と計算し,そのアイデアを日記に記しているそうである.

 これは2次方程式

  3x^2−46x+72=0

の解となっていることが見てとれるが,コラム「泡細胞の幾何学」に示した計算法では2次方程式が現れる気配はまったくなかった.どのようにして導出されたものなのだろうか?

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

【1】ヒーウッドの公式

 オイラーの公式:v−e+f=2は単純ですが,要はその使い方というわけで,たとえば,多面体には3角形か4角形面か5角形面が少なくとも1つなければならないことはもっと簡単に証明できます.

 どの領域も少なくとも6つの領域で囲まれていると仮定すると

  6f≦2e

また,このような問題を解くにあたっては,すべての交点で3本の境界線が会している地図だけを考えればよいので,

  3v≦2e

 これらをオイラーの公式に代入すると

  v−e+f≦1/3e−e+2/3e=0≠2

となって矛盾を生じます.したがって,5個以下の隣接領域しかもたない領域が少なくともひとつあることになります.

 平面や球面上に描かれた地図に関するオイラーの公式は

  v−e+f=2

でしたが,トーラス上の地図に関するオイラーの公式は

  v−e+f=0

です.

 トーラスでは6個以下の隣接領域しかもたない領域が少なくともひとつあることを証明するために,どの領域も少なくとも7つの領域で囲まれていると仮定すると

  7f≦2e

また,3v≦2eですから

  v−e+f≦2/7e−e+2/3e=−1/21e≠0

という矛盾を引き出すことができます.

 したがって,トーラスでは6個以下の隣接領域しかもたない領域が少なくともひとつあることになります.このことを利用すると,

  「トーラス上のどんな地図でも7色で塗り分けられる」

ことが証明されます.ヒーウッドは実際に7色を必要とする例もあげています.

 これを証明したヒーウッドはさらにg個の穴があいたトーラス上の地図に関するオイラーの公式

  v−e+f=2−2g

を利用して

(1)2個の穴があいているトーラス上の地図はどれも8色で塗り分けられる

(2)3個の穴があいているトーラス上の地図はどれも9色で塗り分けられる

・・・・・・・・・・・・・・・・・・・・・・・・

(3)10個の穴があいているトーラス上の地図はどれも14色で塗り分けられる

に引き続いて,

(4)g個の穴があいているトーラス上の地図はどれもH(g)色で塗り分けられる

  H(g)=[{7+√(1+48g)}/2]

を証明しました.

 種数gのトーラスの彩色数が,2次不等式

  x^2−7x+12(1−g)≧0

の解に帰着されることについては,コラム「ヒーウッドの公式について」をご参照ください.

 [・]はガウス記号で,

  g:1,2,3, 4, 5, 6, 7, 8, 9,10

  H:7,8,9,10,11,12,12,13,13,14

となるのですが,しかし,ヒーウッドはg≧2に対してそのような地図が実在することを示すことはできませんでした.そのため,この問題は「ヒーウッド予想」と呼ばれることになりました.

 1968年,リンゲルとヤングスは,g個の穴のあいているトーラス上にこれだけの色を必要とする地図が存在することを証明したのですが,ヒーウッド予想(1890年)が最終的に証明されるまでには77年もの歳月が必要だったというわけです.導出法についてはコラム「ヒーウッドの公式について」をご参照願います.

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

【2】1つの泡に接する泡の数(その1)

[Q]3次元では泡細胞は14面体を中心とした分布をなし(面数の平均は<14),すべての泡細胞が14面以上の面をもつことは不可能である.

[A]3次元の空間が多面体により分割されるとき,3個の多面体の面が合して1本の稜線を形成し(内容的には同じことであるが)4本の稜線が1点に集まる.集合体から1個の多面体を分離して考えてみると,分割多面体の幾何学的性格で最も重要なものは,多面体のいずれの頂点にも3本の稜が集まるということである.そしてこのような多面体で空間を充填すれば,1個の頂点は4個の多面体によって共有され,そこには必ず4本の稜が集まる形になる.

 各分割多面体の頂点,辺,面の数をそれぞれvi,ei,fiとする.

  vi−ei+fi=2

分割多面体では1個の頂点に3本の辺が集まり,また1本の辺は2個の頂点を結ぶことから,

  2ei=3vi

また,集合多面体の頂点,辺,面,胞の数をそれぞれV,E,F,Cとすると

  Σ(vi−ei+fi)=2C

 空間分割の面の数などは一義的には決まらず,統計的にしか扱えないので,各多面体の平均頂点数,辺数,面数v1,e1,f1とおくと

  v1C=Σvi,e1C=Σei,f1C=Σfi

  v1C−e1C+f1C=2C

f1≒14が証明したい事柄である.なお,V≠Σvi,E≠Σei,F≠Σfiであることを注意しておく.

 平均的多面体の各面をp角形,各頂点にq面が会するとし,各辺にr個の多面体(p,q)が集まるものとする.境界多面体(p,q)の平均的な頂点数,辺数,面数は(v1,e1,f1)となる.また,頂点に集まる辺の中点を結んでできる多面体はq角形が1つの辺にr面会した多面体(q,r)になっていて,その図形は(v2,e2,f2)で表されるものとする.

 3次元の握手定理は多彩になって

  f1C=2F,v1C=f2V,v2V=2E,e1C=rE=pF=e2V

であるが,仮定により

  q=r=3,v2=4,e2=6,f2=4

であるから

  f1C=2F,v1C=4V,4V=2E,e1C=3E=pF=6V

 これを

  v1C−e1C+f1C=2C

に代入すると

  (2−p/3)F=2C

  f1=2F/C=12/(6−p)

  v1=4V/C=4/(p/6−1)

  e1=6V/C=6/(p/6−1)

 ここで位相幾何学的証明でなく,計量的証明が必要になるのであるが,3次元空間充填であるためには,等式

  cos(π/q)=sin(π/p)sin(π/r)

が成り立たなくてはならない.q=r=3を代入すると

  sin(π/p)=cot(π/3)=√(1/3)

  p=5.1044・・・

  f1=13.398・・・<14

  v1=22.796・・・

となる.

 なお,

  V−E+F−C=0

を用いても,

  f1=2F/C=12/(6−p)

を導き出すことは可能である.厳密にいうと

  V−E+F−C=1

であるが,何千という小さな泡の集合体を想定して,非常に多くの多面体を統計的に扱うので右辺は0としてかまわない.しかし,V−E+F−C=1を用いればもっとf1≒14に近づくであろう.

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

【3】1つの泡に接する泡の数(その2)

 ここでは,別の方法でもっと簡単に1つの泡に接する泡の数を計算してみます.2つの泡がくっついたとき,その境界も球面になり,この3つの球面の接合角度は120°となります.また,互いに120°の角度で交わる石鹸膜の交線は

  arccos(−1/3)=109.471°

で接触します.正四面体の頂点から中心に向かう3枚の膜は互いに120°の角度をなし,中心に集まる4本の線は109.471°(マラルディの角)をなすのです.

 120°と109.471°は石鹸膜が接触するときの基本的な角度ですが,ここで泡が内角109.471°の正多角形からなる正多面体とみなせば,面は

 p=360/109.471=5.104

角形となります.

 また,3v=pf,2e=pfを,オイラーの多面体定理に代入すると

  pf/3−pf/2+f=2

  f=12/(6−p)

  f=13.39,v=22.78,e=34.18

すなわち,泡の平均の姿は22.78個の頂点,34.18本の辺,13.39枚の面からなる面が5.104角形の立体となることがわかります.平均的な泡細胞は正12面体にやや似たものになるというわけです.

 (その1)でもf=12/(6−p)にp=5.104を代入することによってf=13.39が得られましたが,(その1)で行った計算

  sin(π/p)=√(1/3)

はπ/pが正四面体の2面角の半分に等しいことを意味しているので,(その2)の結果とよく符合するのは当然のことです.

 実験的研究から多面体の面数は14面,面の形は五角形がもっとも多いことが知られていますが,理論的にもかなりよく符合する結果が得られ,実験で得られた値を裏付ける1つの根拠を与えてくれるのです.

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

【4】1つの泡に接する泡の数(その3)

  f=(23+√313)/3=13.56

が2次方程式

  ax^2+2bx+c=0

の解だとすると,

  f=(−b+√(b^2−ac))/a

より,

  a=3,b=−23,c=72

  3x^2−46x+72=0

となります.

 また,

  f=(−b+√(b^2−ac))/a

の分子を有理化すると

  f=c/(−b−√(b^2−ac))

より,

  p=6(f−2)/f=6(c+2b+2√(b^2−ac))/c

   =(26+2√313)/12=5.1153

 pが2次方程式

  Ax^2+2Bx+C=0

の解だとすると,

  p=(−B+√(B^2−AC))/A

より,

  A=12,B=−26,C=−48

  12x^2−52x−48=0

となります.

  p=2π/arccos(−1/3)=5.104

からpの近似値を2次方程式の解として求められればよいのですが,そのような方法が思いつきません.

 また,

  f=(23+√313)/3=13.56,p=5.1153

とf=13.39,p=5.104は微妙な食い違いを示していることから,両者は別の理論的模型に由来しているのではないかと推察されます.これ以上の考察は次回以降の宿題としたいと思います.

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