わかる c·t=一定 第 1 回 / スペック表を作るところから始める

メモリと演算を割り算すると、状態方程式が出てきます

宇宙は、1ビットに
0.035回しか計算していない 「宇宙は有限リソースの計算機だ」── この見方を、比喩ではなく
スペック表として最後まで書き切ります。最初の一枚で、もう妙なことが起きます。

必要な道具:割り算、プランク単位、対数\(\Omega/N=\dfrac{\ln 2}{3\pi^2(1+w)}\)

前シリーズ「わかる共形変換」は、\(c\cdot t=\text{一定}\) を 語り方としては完璧、モデルとしては失格 と判定して閉じました。この新シリーズは、そこで扱いきれなかった一点だけを掘ります ── そもそも \(c\cdot t=\text{一定}\) は、物理を情報理論で見た世界だったという一点です。まずは計算機としてのスペックを取ります。メモリ、クロック、演算数。三つ揃えて割り算すると、宇宙論のいちばん基本的な量が、まったく別の顔で出てきます。

01メモリ ── 地平面に何ビット書けるか

計算機の第一のスペックはメモリです。宇宙の場合、これはホログラフィック原理が答えを持っています ── 情報は体積ではなく面積に書かれる。地平面の面積 \(A\) を、プランク面積で割って \(4\ln 2\) で割ると、ビット数になります。

計算 ── メモリ

地平面のビット数

$$N=\frac{A}{4\ell_P^2\ln 2}=\frac{\pi}{\ln 2}\left(\frac{R_H}{\ell_P}\right)^2$$

数を入れる(\(R_H=c\,t_0=1.305\times10^{26}\) m、\(\ell_P=1.616\times10^{-35}\) m)

$$\frac{R_H}{\ell_P}=8.075\times10^{60}\qquad\Longrightarrow\qquad N=2.956\times10^{122}\ \text{bit}$$

これが宇宙の全メモリです。地球上の全ストレージが \(10^{22}\) ビット程度なので、その \(10^{100}\) 倍。

02クロック ── 何回ティックしたか

第二のスペックはクロック。プランク時間を1ティックと数えます。

計算 ── クロック $$\frac{t_0}{t_P}=\frac{4.354\times10^{17}\,\mathrm{s}}{5.391\times10^{-44}\,\mathrm{s}}=8.075\times10^{60}\ \text{tick}$$

対数で数えると

$$\ln\frac{t_0}{t_P}=140.2$$

メモリの \(R_H/\ell_P\) とクロックの \(t_0/t_P\) が、同じ \(8.075\times10^{60}\) になっているのに気づいたと思います。これが前シリーズ番外編③の 「1ティックに1セル」 ── \(dR_H/dt=c\)、つまり \(c\cdot t=\text{一定}\) そのものです。今回はここには触れず、先へ進みます。

対数のほうが、本当のクロックかもしれない 生の 8×10⁶⁰ という数より、\(\ln(t_0/t_P)=140\) のほうが計算機としては意味があります。宇宙の履歴は「10⁶⁰回の同じ処理」ではなく「140回の倍々」だからです。実際、1対数ティックごとにメモリは \(e^2=7.39\) 倍に増えます(\(N\propto t^2\) なので)。140 回ぶん掛けると \(e^{280}\approx10^{122}\) ── ちゃんとメモリの値が再現します。宇宙は、140 ステップしか走っていない機械です。

03演算数 ── 何回、状態を変えたか

第三のスペックが演算数です。「1回の演算」に物理的な下限を与えるのが マルゴラス=レヴィティン限界 で、エネルギー \(E\) の系が直交する状態へ移るのに最低 \(\pi\hbar/2E\) の時間がかかる、と言います。つまり演算のレートは \(2E/\pi\hbar\) が上限。

地平面の中の全エネルギーは、平坦 FLRW ならいつでも \(E=Mc^2=c^4R_H/2G\) です(番外編③06節の恒等式 ── ディラックの大数仮説の正体だったもの)。膨張則を \(a\propto t^{p}\) とすると \(R_H=ct/p\) なので、レートは時間に比例して増えます。

計算 ── 演算数(三行)

① レート

$$\frac{d\Omega}{dt}=\frac{2E}{\pi\hbar}=\frac{c^4R_H}{\pi\hbar G}=\frac{c^5\,t}{\pi\hbar G\,p}$$

② 積分する

$$\Omega=\frac{c^5t^2}{2\pi\hbar G\,p}$$

③ \(\ell_P^2=\hbar G/c^3\) を使うと、単位が全部消える

$$\Omega=\frac{1}{2\pi p}\left(\frac{ct}{\ell_P}\right)^2$$

\(p=1\)(\(c\cdot t=\text{一定}\))で数を入れると \(\Omega=1.038\times10^{121}\) 回。ロイド(2002)の有名な見積もり「宇宙は \(10^{120}\) 回の演算をした」と、同じ桁です。

◇ ◇ ◇

04核心 ── メモリを演算で割る

スペックが三つ揃いました。計算機を評価するとき、いちばん見たい数はこれです ── 1ビットあたり、何回演算したか。メモリばかり大きくて演算しない機械なのか、その逆なのか。

やってみると、驚くことが起きます。

計算 ── 割るだけ $$\frac{\Omega}{N}=\frac{1}{2\pi p}\left(\frac{ct}{\ell_P}\right)^{2}\ \Bigg/\ \frac{\pi}{\ln 2}\left(\frac{ct}{p\,\ell_P}\right)^{2} =\frac{p\,\ln 2}{2\pi^{2}}$$

\(t\) が、きれいに消えました。 メモリも演算数もどちらも \(t^2\) で増えるので、比は宇宙の年齢に依りません。宇宙が生まれて1秒でも、いまでも、\(10^{100}\) 年後でも、1ビットあたりの演算回数はまったく同じ

そして残った \(p\) は、膨張則そのものです。状態方程式 \(w\) との関係 \(p=2/3(1+w)\) を入れると ──

この回の結論

$$\boxed{\ \frac{\Omega}{N}=\frac{\ln 2}{3\pi^{2}\,(1+w)}\ }$$

1ビットあたりの演算回数は、状態方程式だけで決まる純粋な数。宇宙の年齢にも、大きさにも、中身の詳細にも依らない。

宇宙論でいちばん基本的な量(状態方程式)と、計算機科学でいちばん基本的な比(演算数 ÷ メモリ)が、同じ一つの数だった、ということです。

判定手続きにかける 前シリーズ最終回の作法どおり、この量を二列に分けます。\(\Omega\) も \(N\) もただの個数で、次元がありません。だから共形変換(=物差しの取り替え)では動かない ── 右の列、物理です。しかも動かないどころか、\(w\) を直接読んでいる。「宇宙の膨張則を測る」ことと「宇宙の演算効率を測る」ことが、同じ作業だったということになります。

05四人を並べる

宇宙\(w\)\(p\)\(\Omega/N\)1回あたり
stiff(キネーション)\(+1\)1/30.0117185.4 ビットに1回
放射\(+1/3\)1/20.0175657.0 ビットに1回
物質\(0\)2/30.0234142.7 ビットに1回
\(c\cdot t=\)一定\(-1/3\)10.0351228.5 ビットに1回
ド・ジッター\(-1\)発散時間とともに増える

どれも 1 に遠く届きません。宇宙は 138 億年かけて、一度も「1ビットあたり1回」を計算したことがない。 しかもそれは「まだ若いから」ではありません。比は時間に依らないので、いつの時代もそうだったのです。

068つ目の特徴づけ

ここで \(c\cdot t=\text{一定}\) の行を、もう少し正確に言い直します。\(p=ct/R_H\) ですから、結論の式はこう書けます。

同じ式を、比で書く $$\frac{\Omega}{N}=\frac{\ln 2}{2\pi^{2}}\cdot\frac{ct}{R_H}$$

よって

$$R_h=ct\qquad\Longleftrightarrow\qquad \frac{\Omega}{N}=\frac{\ln 2}{2\pi^{2}}=0.0351152\cdots$$

前シリーズ番外編③は、\(a\propto t\) の特徴づけを七つ並べて「これだけの性質が一点に集まるのは、それ自体が驚くべきこと」と書きました。八つ目を足せます。

ポテンシャルが純粋な二次形式\(U=-\phi^2/2t_0^2\)(共形変換 第5回)
共形時間が \((-\infty,\infty)\) 全体ミンコフスキーをまるごと覆う(第6回)
加速がゼロ\(\ddot a=0\)。地平線問題の分水嶺(第5回)
共動ハッブル半径が一定アドレス空間が動かない(番外編②)
粒子的地平線と事象の地平線が両方無限全メモリに読み書き両方のアクセス
\(dR_H/dt=c\)、セル数=ティック数1ティックに1セル(番外編③)
強いエネルギー条件の飽和\(\rho+3p=0\)。特異点定理の境界(番外編③)
1ビットあたりの演算回数が \(\ln 2/2\pi^2\)演算とメモリの比が、光が走った距離と地平半径の比に一致する(今回)

⑥と⑧は、同じ事実の別の顔です。⑥が「長さと時間」で言っているものを、⑧は「演算とビット」で言っている ── 次元付きの言い方と、無次元の言い方。前シリーズの判定手続きで言えば、⑧のほうが右の列にいます。

07では、「1」に届く宇宙はあるのか

結論の式を逆に解いてみましょう。\(\Omega/N=1\) を要求すると ──

逆に解く $$1+w=\frac{\ln 2}{3\pi^{2}}=0.02341\qquad\Longrightarrow\qquad \boxed{\,w=-0.97659\,}$$

ほとんど \(-1\)。つまり 宇宙定数のすぐ手前です。そして観測されている暗黒エネルギーの状態方程式は \(w=-1.03\pm0.03\) ── すでにその向こう側にいます

いま宇宙は暗黒エネルギー優勢期に入ったところです。ド・ジッター的な膨張では \(N\) が一定に飽和する一方で演算レートは落ちないので、比は時間とともに増え続けます。

ド・ジッターでの比 $$\frac{\Omega}{N}=\frac{\ln 2}{\pi^{2}}\,Ht\qquad\Longrightarrow\qquad \frac{\Omega}{N}=1\ \text{になるのは}\ Ht=\frac{\pi^{2}}{\ln 2}=14.24$$

\(H_\Lambda=H_0\sqrt{\Omega_\Lambda}=1.81\times10^{-18}\,\mathrm{s^{-1}}\) を入れると

$$t\approx2.5\times10^{11}\ \text{年}\qquad(\text{約 2500 億年後})$$

今回いちばん言いたいこと

宇宙は 138 億年、一度も「1ビットあたり1回」を計算していない
そしていま、初めてその境界を越える側の膨張則に入った。

図:横軸は宇宙の対数年齢(\(\ln(t/t_P)\)、今日が 140.2)。冪則の宇宙では 1ビットあたりの演算回数は完全に水平で、\(w\) のツマミを動かすと線の高さだけが変わる。銅の破線が「1 演算/ビット」── 冪則の宇宙がここに届くのは \(w<-0.977\)、つまり宇宙定数のすぐ手前に入ってからだけ

w = −0.33
1ビットあたりの演算回数(冪則) 1 演算/ビット Λ優勢の未来

ツマミを左端(\(w\to-1\))に寄せていくと、水平線がゆっくり持ち上がります。それでも銅の破線に触れるのは \(w=-0.977\) を越えてから ── 宇宙定数のすぐ手前でしか、計算機は「割に合わない」状態を脱しない。右へ振るほど(物質・放射・stiff)線は下がる一方です。灰色の破線が、今日から先の Λ 優勢の未来。ここだけが水平ではありません。

◇ ◇ ◇

08種明かし ── なぜ両方 \(t^2\) だったのか

比が時間に依らなかったのは、\(N\) も \(\Omega\) も \(t^2\) だったからです。これは偶然ではありません。

なぜ \(t^2\) か
メモリ \(N\)面積だから。\(N\propto R_H^2\)、\(R_H\propto t\)。ホログラフィック原理そのもの
演算 \(\Omega\)レートが \(E\propto R_H\propto t\) に比例し、それを積分するから。\(E=c^2R_H/2G\times c^2\) という恒等式が効いている

後者の恒等式は、前シリーズ番外編③06節で「ディラックの大数仮説はただの恒等式だった」と片付けたあの式です。ハッブル半径は、いつでも中身のシュヴァルツシルト半径にぴったり等しい。 だから宇宙のエネルギーは半径に比例し、演算レートも半径に比例する。

まとめると、こういう連鎖です。

連鎖 $$\underbrace{N\propto R_H^{2}}_{\text{ホログラフィー}}\qquad \underbrace{\frac{d\Omega}{dt}\propto E\propto R_H}_{\text{大数の恒等式}}\qquad\Longrightarrow\qquad \frac{\Omega}{N}\propto\frac{\int R_H\,dt}{R_H^{2}}=\text{(膨張則だけの関数)}$$

「演算数がメモリに比例する」のは、ホログラフィー と 大数の恒等式 の積でした。どちらも \(c\cdot t=\text{一定}\) に固有の性質ではなく、平坦 FLRW ならいつでも成り立つ ── だから比は、膨張則の情報しか残さない。

09おまけ ── メモリ増設のほうが速い

同じ数を、レートで読み直しておきます。

今日の値(\(p=1\)) $$\frac{dN}{dt}=\frac{2N}{t_0}=1.36\times10^{105}\ \mathrm{bit/s},\qquad \frac{d\Omega}{dt}=\frac{2\Omega}{t_0}=4.77\times10^{103}\ \mathrm{ops/s}$$ $$\frac{dN/dt}{d\Omega/dt}=\frac{2\pi^{2}}{p\ln 2}=28.5$$

宇宙は、1回演算するあいだに 28.5 ビットのメモリを増設している。 計算機としては異常な設計です ── メモリを積み続けるだけで、そこに書かれたものを処理する速度がまるで追いつかない。前シリーズ番外編②が「メモリ \(10^{122}\)、クロック 140 回」を極端にメモリ過多・クロック過少と呼んだのは、この比のことでした。今回はその「過多」に、\(2\pi^2/\ln2\) という値が付いたことになります。

しかも、ほとんど空 番外編②の計算では、今日の宇宙が実際に使っているエントロピーは容量の \(1.5\times10^{-18}\) しかありません。つまり ── 28.5 ビットに 1 回しか演算しない機械が、そのメモリの \(10^{-18}\) しか使っていない。書き込みも計算も、限界からはるかに遠い。宇宙は資源が有限な計算機ではありますが、その有限な資源すら、まったく使い切っていないのです。この「空っぽさ」が何を意味するかは、第3回で扱います。
正直な線

「演算」の定義は、マルゴラス=レヴィティン限界による 直交状態への遷移回数の上限 です。意味のある計算をしている、という主張ではありません。宇宙が実際にその上限で走っているという証拠もない ── これはスペック表であって、ベンチマークではない

係数は規約に敏感です。レートを \(2E/\pi\hbar\) と取り、\(E\) を地平面内の全エネルギー \(Mc^2=c^4R_H/2G\) と取ると \(\Omega=1.04\times10^{121}\) になります。ロイド(2002)の \(\sim10^{120}\) との差は、この取り方の違いです。ビット数を \(\ln2\) で割るか(bit)割らないか(nat)でも \(\Omega/N\) は 1.44 倍動きます。したがって「0.035」という数字そのものより、「時間に依らず、\(w\) だけで決まる」という構造のほうが本体です。

\(w\) は宇宙全体の実効的な状態方程式であって、単一成分の話ではありません。実際の宇宙は放射→物質→Λ と乗り換えていくので、\(\Omega/N\) も区間ごとに違う値を取り、正確には \(\int R_H dt/R_H^2\) を数値で追う必要があります。表の値は「その \(w\) がずっと続いたら」という理想化です。

07節の「約 2500 億年後」は、今日以降を純粋なド・ジッターと近似した見積もりです。また \(w=-1.03\pm0.03\) は \(w\) を定数と仮定したフィットの値で、\(w<-1\) 側(ファントム)は将来の測定で動きうる領域です。「すでに向こう側」は中心値についての言い方であって、確定した事実ではありません。

練習問題(今回の式だけで解けます)

  1. 放射優勢期(\(w=1/3\))の宇宙では、1ビットあたり何回演算するか。
    答えを見る
    \(\Omega/N=\ln2/(3\pi^2\cdot4/3)=\ln2/(4\pi^2)=0.01756\)。57 ビットに 1 回。\(c\cdot t=\text{一定}\) のちょうど半分で、放射のほうがさらに「計算しない」宇宙です。
  2. なぜ \(\Omega/N\) は宇宙の年齢に依らないのか。一行で。
    答えを見る
    \(N\propto R_H^2\propto t^2\)(ホログラフィー)で、\(\Omega\propto\int E\,dt\propto\int R_H\,dt\propto t^2\)(\(E=c^4R_H/2G\) の恒等式)だから。両方が同じ \(t^2\) なので、比には膨張則の情報しか残りません。
  3. \(\Omega/N\) は帳簿か、物理か。
    答えを見る
    物理。\(\Omega\) も \(N\) もただの個数で無次元なので、共形変換で動きません。しかも \(w\) を一意に決めるので、動くなら本当に別の宇宙です ── 前シリーズ最終回の判定手続きで、右の列に落ちます。
  4. \(\Omega/N=1\) を満たす \(w\) を求め、それが何を意味するか述べよ。
    答えを見る
    \(1+w=\ln2/3\pi^2=0.0234\) より \(w=-0.9766\)。宇宙定数のすぐ手前です。冪則で膨張する宇宙が「1ビットあたり1回」に達するには、ほとんど加速膨張しきっている必要がある ── 逆に言えば、物質や放射が支配する宇宙は原理的に計算機として割に合いません。
  5. (やや難)⑥「1ティックに1セル」と⑧「\(\Omega/N=\ln2/2\pi^2\)」が同じ事実であることを示せ。
    答えを見る
    \(\Omega/N=(\ln2/2\pi^2)(ct/R_H)\) と書ける。⑥は \(dR_H/dt=c\)、これを積分すると \(R_H=ct\) すなわち \(ct/R_H=1\)。代入すれば \(\Omega/N=\ln2/2\pi^2\)。逆も同じ。「地平面が光と同じ速さで進む」と「演算とメモリの比が \(\ln2/2\pi^2\)」は、\(ct/R_H=1\) という一つの式の二つの読み方です。前者は次元付きの言い方、後者は無次元の言い方。

まとめ 割り算したら、状態方程式が出た

宇宙を計算機としてスペック表にすると、メモリ \(N=\pi(R_H/\ell_P)^2/\ln2=2.96\times10^{122}\) ビット、クロック \(8.08\times10^{60}\) ティック(対数では 140 回)、演算 \(\Omega=(ct/\ell_P)^2/2\pi p=1.04\times10^{121}\) 回。ここまでは既知の見積もりの再確認です。面白いのは、メモリと演算を割ったときでした。

両方が \(t^2\) で増えるので時間が消え、残るのは膨張則だけ ── \(\Omega/N=\ln2/3\pi^2(1+w)\)。1ビットあたりの演算回数は、状態方程式そのものです。放射で 1/57、物質で 1/42.7、\(c\cdot t=\text{一定}\) で 1/28.5宇宙は 138 億年かけて、一度も「1ビットあたり1回」に達していないし、比は時間に依らないので、いつの時代もそうでした。

そして \(\Omega/N=(\ln2/2\pi^2)(ct/R_H)\) と書けるので、\(R_h=ct\) はこの比が \(\ln2/2\pi^2\) にちょうど等しいこと ── 番外編③の七つの特徴づけに、八つ目が足せます。⑥「1ティックに1セル」の、無次元版です。「1」に届くには \(w=-0.977\) が必要で、観測されている暗黒エネルギーはその向こう側。純ド・ジッター近似なら、宇宙が初めて 1 演算/ビットを超えるのは約 2500 億年後です。

種明かしは二つの恒等式でした ── メモリが \(t^2\) なのはホログラフィー、演算が \(t^2\) なのはハッブル半径=シュヴァルツシルト半径(ディラックの大数)。どちらも平坦 FLRW なら常に成り立つので、割り算の結果には膨張則しか残らなかった。おまけに、メモリの増設は演算より \(2\pi^2/\ln2=28.5\) 倍速く、しかもそのメモリは \(10^{-18}\) しか使われていません。

この文書は「わかる c·t=一定」シリーズ第1回、物理好きの高校生・大学生向け読み物です。地平面エントロピー \(S=k_BA/4\ell_P^2\)、マルゴラス=レヴィティン限界(直交状態への遷移レートの上限 \(2E/\pi\hbar\))、および宇宙の総演算数が \(\sim10^{120}\) 程度になることは Lloyd (2002, PRL 88, 237901) 以来の標準的な見積もりです。平坦 FLRW における \(M=c^2R_H/2G\) は恒等式で、Dirac の大数仮説がこれと同型であることは本シリーズ前作の番外編③で示した通りです。本稿の \(\Omega=(ct/\ell_P)^2/2\pi p\)、\(\Omega/N=p\ln2/2\pi^2=\ln2/3\pi^2(1+w)\)、\(\Omega/N=1\) を与える \(w=-0.97659\)、およびド・ジッターでの \(\Omega/N=(\ln2/\pi^2)Ht\)(\(Ht=\pi^2/\ln2=14.24\) で 1 に達する)は、いずれも本稿での導出・数値計算です。数値は \(t_0=4.3536\times10^{17}\) s、\(\ell_P=1.6163\times10^{-35}\) m、\(t_P=5.3912\times10^{-44}\) s を用いました。「演算」はエネルギーが許す遷移回数の上限であり、意味のある計算を指すものではありません。 レートの規約(\(2E/\pi\hbar\) か \(4E/h\) か)、\(E\) の取り方、bit と nat の別によって \(\Omega/N\) の数値は 1.5 倍程度動きます ── 時間に依存せず \(w\) のみで決まるという構造が本稿の主張であり、数値そのものではありません。\(w\) は宇宙全体の実効的な状態方程式で、表の値は「その \(w\) が全時代続いたら」という理想化です。約 2500 億年という見積もりは純ド・ジッター近似(\(H_\Lambda=H_0\sqrt{\Omega_\Lambda}\)、\(H_0=67.66\) km/s/Mpc、\(\Omega_\Lambda=0.685\))によります。暗黒エネルギーの \(w=-1.03\pm0.03\) は定数 \(w\) を仮定したフィットの値です。エントロピー使用率 \(1.5\times10^{-18}\) は Egan & Lineweaver (2010, ApJ 710, 1825) の \(S_{\rm obs}=3.1\times10^{104}k_B\) によります。線形膨張(\(c\cdot t=\)一定、\(R_h=ct\))は検証途上の少数派モデルで、初期宇宙まで外挿した場合は元素合成と矛盾します(Lewis, Barnes & Kaushik 2016, MNRAS 460, 291)。学術的な標準はインフレーションを含む \(\Lambda\)CDM モデルです。 ── 印刷する場合はブラウザの「印刷」から「PDF に保存」を(印刷版ではスライダーと解答は静止・非表示になります)。

印刷 / PDF 化:⌘+P(Windows は Ctrl+P)。画面ではスライダーで状態方程式を変え、1ビットあたりの演算回数が上下する様子が見えます。「答えを見る」で解答が開きます。