考える波 第 27 回 ／ 第 XI 部・仮説に決着をつける
前回 残した最後の謎に、答えが出ます
対数一様は、どこから来るのか ──
自然が用意していたのは、指数関数だった
1/f を作るのに、べき乗分布を用意する必要はありません。指数の中に「ふつうの量」が入っていれば、それだけでいい。
★ しかも \(\beta-1\) は二つの桁数の比なので、指数は勝手に 1 に釘づけされます。
そして予想が一つ、外れました ── 離散スケール不変性は 3 dB の痕跡を残します。
前回の結論は「分数階が現れるのは、粗視化した自由度がスケール不変に分布しているときである」でした。
そして 1/f 雑音だけが、形を変えて残りました ── ★ なぜ自然には、対数一様な時定数分布がこれほど多いのか。
本回で答えます。
01まず確認 ── 対数一様から、本当に 1/f が出るか
このときの積分は手で解けます
数値で確かめます（\(\tau\in[10^{-6},10^{6}]\) 秒 ＝ 12 桁 を対数等分）。
| \(\omega\) [rad/s] | \(S(\omega)\) | 傾き [/dec] | 予言 |
|---|---|---|---|
| \(10^{-4}\) | 3.1216e+4 | ── | 帯の外 |
| \(10^{0}\) | 3.141589 | −1.0000 | ★ \(\pi A/\omega=3.141593\)（誤差 1.3e−6） |
| \(10^{3}\) | 3.1396e−3 | −1.0002 | −1 |
| \(10^{4}\) | 3.1216e−4 | −1.0025 | 帯の外 |
01節の中心
対数一様 ⇒ 1/f は確かです（帯の中で傾きちょうど −1、振幅も \(\pi A/\omega\) と 6 桁 一致）。
★ そして帯の外では平らに寝ます ── 1/f には必ず上下の折れ点があります。無限には続きません。
02★ 本題 ── 対数一様はどこから来るのか
\(g(\tau)\propto1/\tau\) は、\(\tau\) を直接 見ると特別な形に見えます。ところが ── 変数を一つ替えるだけで、いちばんありふれた形になります。
\(\ln\tau\) は障壁 \(E\) の一次関数です
\(D\) が平らなら、そのまま \(1/\tau\)
02節の中心
★★ \(E\) が一様 ⟺ \(\ln\tau\) が一様 ⟺ \(g(\tau)\propto1/\tau\)
べき乗分布を用意する必要はありません。「指数の中身が平らである」だけでいい ── これはごく普通の状況です。
どれくらい桁が稼げるか（\(\tau_0=10^{-13}\) s、室温 300 K）。1 桁 あたりの障壁は \(kT\ln10=0.05953\) eV。
| 障壁 \(E\) の幅 [eV] | \(\tau\) の広がり [桁] | \(\tau_{\max}\) |
|---|---|---|
| 0.10 | 1.68 | 5.2e−7 s |
| 0.36 | ★ 6.0 | ── |
| 0.70 | 11.76 | 6.3e+3 s ＝ 1.75 時間 |
| 1.00 | 16.80 | 6.9e+8 s |
★ 1/f が 6 桁 続くのに必要なのは、0.36 eV ぶんの障壁のばらつきだけ。安すぎます。
そして \(E_{\max}=1.0\) eV なら \(\tau_{\max}=1.75\) 時間 ── 1/f は \(2.5\times10^{-5}\) Hz 以下で折れるはず。実験で \(10^{-6}\) Hz まで折れないなら、障壁は 1.3 eV 以上まで伸びています。
03なぜ「1」なのか ── 指数は二つの桁数の比で決まる
\(D(E)\) が平らでなかったらどうなるか。\(D(E)\propto e^{cE}\) として、帯の中で \(D\) が 10 倍・1000 倍 変わる場合を計算します。
| \(D(E)\) の変化 | 測った \(\beta\) | 予言 \(1+kT\,d\ln D/dE\) |
|---|---|---|
| 平ら | 0.99998 | 1.00000 |
| 10 倍 | 1.08499 | 1.08504 |
| 1000 倍 | 1.25485 | 1.25511 |
03節の中心
この例では \(\tau\) が 11.8 桁 広がるので ── \(D\) が 10 倍 変わっても \(\beta-1=1/11.8=0.085\)、1000 倍 でも 0.255。
★★ 分母が大きいので、\(\beta\) は勝手に 1 の近くに釘づけされます。
実測の 1/f が \(\beta=0.8\)〜\(1.4\) に散らばるのは、まさにこの比です。── ★ 1/f は微調整の結果ではありません。避けるほうが難しい。
04実例① \(\alpha\) 崩壊 ── 指数写像のいちばん極端な例
第 19 回で「崩壊定数は虚数の周波数」と見ました。その崩壊定数こそ、指数写像の教科書的な例です ── \(\lambda\propto e^{-G}\)、\(G\) はガモフ因子 \(\propto Z/\sqrt Q\)。
| 核種 | \(Q\) [MeV] | 半減期 | \(\log_{10}T_{1/2}\) | 直線の予言 |
|---|---|---|---|---|
| \(^{232}\)Th | 4.0816 | 4.43e+17 s | 17.647 | 17.784 |
| \(^{238}\)U | 4.2700 | 1.41e+17 s | 17.149 | 17.778 |
| \(^{226}\)Ra | 4.8710 | 5.05e+10 s | 10.703 | 10.867 |
| \(^{222}\)Rn | 5.5900 | 3.30e+5 s | 5.519 | 5.688 |
| \(^{214}\)Po | 7.8340 | 164.3 µs | −3.784 | −3.700 |
| \(^{212}\)Po | 8.9540 | 0.299 µs | −6.524 | −6.552 |
残差の標準偏差 0.575 桁（24 桁 の範囲に対して）
04節の中心
★★ \(Q\) は 2.19 倍 しか違わないのに、半減期は 24.2 桁 違う。
★ 第 19 回との接続：あそこでは「崩壊定数は虚数の周波数だから観測窓の壁が効かない」と書きました。本回ではもう一段 深く ── 崩壊定数が宇宙年齢から µ 秒まで連続に埋まっているのは、\(Q\) という「ふつうの量」が指数の中に入っているからです。
05実例② 離散スケール不変性 ── ★ ここで予想が外れました
対数一様には、もう一つの作り方があります ── 等比数列 \(\tau_n=\tau_0\lambda^n\)。第 13 回のエフィモフ状態（\(\lambda=e^{\pi/s_0}=22.69\)）がまさにこれ。1/f は出るのか、出るならどう違うのか。
| 比 \(\lambda\) | 1 桁 あたりの段数 | リップル [dB p-p] | 予言 \(2/\cosh(\pi^2/\ln\lambda)\) |
|---|---|---|---|
| 2 | 3.32 | 0.00004 | 0.00005 |
| 3 | 2.10 | 0.00871 | 0.00872 |
| 10 | 1.00 | 0.95608 | 0.95662 |
| 22.69（エフィモフ） | 0.74 | ★ 2.951 | 2.967 |
| 100 | 0.50 | 8.301 | 8.700 |
どの \(\lambda\) でも 1/f 自体は出ます。違いは対数周期のリップルで、\(\lambda\) が大きいほど大きい ── そして予言の式と \(\lambda\le22.7\) では 3 桁 で一致します。
05節の中心 ── ★ 予想が外れた
エフィモフ比のリップルは 2.95 dB。第 13 回の「対数周期性は見つけにくい」とは逆の結論です。
★ 第 13 回で消えていたのは 虚数階フィルタのインパルス応答（\(e^{-\beta\pi}\) で潰れる）。
★ 本回で残るのは 離散した時定数の梯子がスペクトルに刻むリップル。
── 「離散スケール不変性」と「虚数階」は、別の見え方をします。
これは検証可能な予言になります ── 実測の 1/f が対数周期リップルを見せないなら、梯子の比は小さいはず。
| リップルの検出限界 | 許される最大の \(\lambda\) |
|---|---|
| 1.00 dB 以下 | 10.24 |
| 0.10 dB 以下 | ★ 4.52 |
| 0.03 dB 以下 | 3.57 |
★★ 1/f が 0.1 dB の平坦さで測れているなら、時定数の梯子の比は 4.5 倍 以下。── つまり 1/f の背後にあるのは連続（密）な分布であって、エフィモフ型の粗い梯子ではありません。02 節の活性化過程と整合します。
06対照実験 ── 条件はもっと緩かった
\(\ln\tau\) がガウス分布（幅 \(\sigma\) 桁）の場合。測るのは「傾き」ではありません ── \(\omega\tau=1\) では単一のローレンツでも傾きはちょうど −1 になってしまう（肩の点）。測るべきは幅：傾きが \(-1\pm0.1\) に収まる \(\omega\) の範囲が何桁あるか。
| \(\sigma\) [桁] | 1/f に見える幅 [桁] | 判定 |
|---|---|---|
| ≈0（単一 \(\tau\)） | 0.10 | ★ 肩の一点だけ。1/f ではない |
| 0.5 | 0.26 | 中間 |
| 1.0 | 0.62 | 中間 |
| 2.0 | 2.02 | ★ 1/f と区別がつかない |
| 4.0 | 7.58 | ★ 1/f と区別がつかない |
06節の中心
\(\sigma\) を広げると、1/f に見える幅はほぼ \(\sigma\) に比例して伸びます（\(\sigma=4\) 桁 で 7.6 桁）。
── ★ 「対数一様」である必要すらありません。対数軸で広ければ、それで 1/f に見える。自然が用意すべきものは、さらに緩かった。
07仕組みを並べて採点する
| 仕組み | 形 | 必要な条件 | 対数一様が出るか |
|---|---|---|---|
| 活性化（アレニウス） | \(\tau=\tau_0e^{E/kT}\) | 障壁 \(E\) が平らに散らばる | ◎ ほぼ自動 |
| トンネル効果 | \(\tau\propto e^{2\kappa d}\) | 距離 \(d\) が平らに散らばる | ◎ ほぼ自動 |
| \(\alpha\) 崩壊（第 19 回） | \(\lambda\propto e^{-G}\) | \(Q\) が並ぶだけ | ◎ 実測（24 桁） |
| 掛け算的ランダム過程 | 対数正規 | \(\sigma\) が 2 桁 超で十分 | ○ 条件つき |
| 離散スケール不変性 | \(\tau_n=\tau_0\lambda^n\) | 第 13 回の梯子 | ◎ ただし 3 dB のリップル |
| 臨界現象 | \(\tau\propto\xi^z\) | 相関長が発散するとき | ○ ★ 微調整が要る |
| 一様なばね鎖（第 26 回） | \(\tau\) が一定 | ── | × 整数階のまま |
07節の中心
上から三つは、どれも「指数の中に、ふつうの量が入っている」だけ。べき乗分布を用意する必要も、臨界点に微調整する必要もありません。
★ 対して臨界現象は温度をぴったり \(T_c\) に合わせないとスケール不変になりません ── 自己組織化臨界（SOC）が注目されたのは、この微調整を系が勝手にやってくれる仕組みだったからです。
でも 1/f の説明としては、活性化のほうがずっと安上がりです。
練習問題
- 1/f 雑音の \(\beta\) が温度を上げると 1 に近づく系と、遠ざかる系がある。本回の式でどう読むか。
答えを見る
\(\beta-1=kT\,d\ln D/dE\) の両側が温度で動くからです。分子の \(kT\) は温度に比例して増えるので、\(D(E)\) が傾いていれば \(\beta\) は 1 から遠ざかります。
一方、温度を上げると観測周波数に効く障壁の窓 \(E^\ast=kT\ln(1/\omega\tau_0)\) が高いほうへ移動するので、そこで \(D\) がたまたま平らなら \(\beta\) は 1 に近づきます。
★ つまり \(\beta(T)\) を測ることは、\(D(E)\) を走査していることになります ── これがダッタ＝ホーンの解析の考え方です。本稿はその形を数値で確認しただけで、実測データの解析はしていません。 - 1/f が \(10^{-6}\) Hz まで折れずに続く、という報告から何が言えるか。
答えを見る
障壁の上端が 1.3 eV 以上まで伸びていると言えます。\(\tau_{\max}=\tau_0e^{E_{\max}/kT}\) で \(\tau_0=10^{-13}\) s、300 K なら \(E=1.0\) eV で \(\tau=1.75\) 時間（\(2.5\times10^{-5}\) Hz）、\(E=1.3\) eV で \(6.9\times10^{8}\) s。
★ そして「無限に続く 1/f」は存在しません ── 総パワー \(\int(1/f)df\) は対数発散するので、必ずどこかで折れます。折れ点が見えないのは、測定時間が足りないだけです。 - 第 13 回で「対数周期性は見つけにくい」と結論したのに、本回で 3 dB のリップルが出たのはなぜか。矛盾か。
答えを見る
矛盾ではなく、別のものを測っています。第 13 回で潰れていたのは虚数階フィルタのインパルス応答で、振幅が \(e^{-\beta\pi}\) で落ちる（\(\beta=3\) で \(10^{-6}\)）。
本回で残るのは離散した時定数の梯子がスペクトルに刻む周期的なうねりで、大きさは \(2/\cosh(\pi^2/\ln\lambda)\)。
★ 前者は「連続な虚数階が作る対数振動」、後者は「離散した実数階が作る対数振動」 ── 同じ \(\ln\omega\) 軸の上の話でも、出どころが違います。そして後者のほうがずっと見つけやすい。
まとめ 自然が用意していたのは、べき乗分布ではなく指数写像だった
第 26 回が残した問い ──「なぜ自然には、対数一様な時定数分布がこれほど多いのか」に答えました。
まず対数一様 ⇒ 1/f を確認（帯の中で傾きちょうど −1、振幅も \(\pi A/\omega\) と 6 桁 一致）。ただし帯の外は平ら ── 1/f には必ず上下の折れ点があります。
★ そして答えは指数写像でした。\(\tau=\tau_0e^{E/kT}\) なら \(\ln\tau\) は \(E\) の一次関数なので、\(g(\tau)=kT\,D(E)/\tau\) ── \(D\) が平らなだけで \(1/\tau\) が出ます。べき乗分布を用意する必要はありません。1 桁 あたりの障壁は室温で 0.0595 eV なので、0.36 eV のばらつきで 6 桁 の 1/f。
★★ 指数が 1 に釘づけされる理由も出ました ── \(\beta-1\) ＝（\(D\) が変化する桁数）÷（\(\tau\) が広がる桁数）。\(D\) が 1000 倍 変わっても \(\beta-1=0.25\)。実測の 1/f が 0.8〜1.4 に散らばるのは、まさにこの比です。
\(\alpha\) 崩壊はこの指数写像の極端な例で、\(Q\) は 2.19 倍 しか違わないのに半減期は 24.2 桁（ガイガー＝ナタル、\(R^2=0.9955\)）── 第 19 回の「崩壊定数が連続に埋まっている」理由がこれです。
★ そして予想が一つ外れました。離散スケール不変性（\(\tau_n=\tau_0\lambda^n\)）は 1/f を作りますが、エフィモフ比 22.69 では 2.95 dB の対数周期リップルを残します（予言 \(2/\cosh(\pi^2/\ln\lambda)\) と 3 桁 一致）。第 13 回で潰れていたのは虚数階のインパルス応答で、これは別物でした。逆にこれは検証可能な予言になります ── 1/f が 0.1 dB 平坦なら、梯子の比は 4.5 倍 以下、つまり背後は連続な分布です。
★ 1/f は微調整の結果ではなく、避けるほうが難しい現象だった。
速さが「ふつうの量の指数」で決まる系なら、勝手に対数一様になる。
★★ そして第 26 回の扉が閉じます ── 第 26 回「分数階 ⇐ スケール不変な自由度」、第 27 回「スケール不変な自由度 ⇐ 指数写像」。分数階の起源は、指数関数だった。
★ 次回は、その指数がどこから来るのかを追います ── トンネル効果も、熱活性化も、ボルツマン因子も、どれも \(e^{-(\text{何か})}\) の形。なぜ自然はこんなに指数関数だらけなのか。
★ 波として見ると答えは一つです ── 減衰する波はすべて指数だから。第 6 回のエバネッセント波が、ここで戻ってきます。
この文書は「考える波」シリーズ第 27 回です。計算は kenshou/taisuu.js で行っています。本回の数値について：スペクトルの重ね合わせ・障壁分布からの \(\beta\)・梯子のリップル・対数正規の幅は、すべて本稿で数値計算したものです。\(\alpha\) 崩壊の 9 核種の \(Q\) 値と半減期は文献値（PDG／NNDC の代表値）で、最小二乗の傾き 1.5065 と \(R^2=0.9955\) が本稿の計算です。ガイガー＝ナタルの法則も、\(g(\tau)=kT\,D(E)/\tau\) という活性化からの導出（マクウォーターの模型）も、ダッタ＝ホーンの \(\beta(T)\) 解析も、いずれもよく知られた結果で本稿の発見ではありません ── 本稿がしたのは、これらを第 26 回の残した問いへの答えとして並べ、\(\beta-1\) が桁数の比になることとリップルの大きさを数値で確かめたことです。リップルの予言式 \(2/\cosh(\pi^2/\ln\lambda)\) はメリン再和の \(k=1\) 項で、\(\lambda=100\) で 5 % ずれるのは高次項（\(k\ge2\)）のためです。「1/f の起源」は依然として系ごとに議論がある問題で、本稿は「活性化過程で説明できる系が多い」と述べているだけです。 ── 印刷する場合はブラウザの「印刷」から「PDF に保存」を。