Stationary $1/f^α$ noise in discrete models of the Kardar-Parisi-Zhang class

離散モデルにおけるカルダール・パリジ・ザン（KPZ）普遍性クラスの成長界面の高さ変動を解析した結果、有限系において定常状態に達し、時間相関関数が動的スケーリングに従い、パワースペクトルが非自明な周波数領域で指数 5/3 の $1/f^\alpha$ 型スケーリングを示すことが確認され、これにより変動が広義定常であることが立証された。

要約

この論文は、**「成長する表面の『ざらつき』が、実はある法則に従って静かに揺らいでいる」**という発見について書かれています。

専門用語を避け、身近な例え話を使って解説しますね。

1. 物語の舞台：砂山と雪だるまの成長

まず、想像してみてください。砂を一つずつ積み上げて「砂山」を作っている場面を。
あるいは、雪だるまに雪玉を転がして大きくしていく場面です。

• KPZ 普遍性クラス（カーダール・パリー・ザング）： これは、砂山や雪だるま、あるいは液体の表面など、「不規則に成長するもの」に共通するルールの名前です。物理学者たちは、これらが実は同じような「成長の癖」を持っていることに気づきました。
• ざらつき（粗さ）： 表面は平らではなく、ボコボコしています。このボコボコの「高さ」が時間とともにどう変わるかを観察するのがこの研究です。

2. 過去の謎：「1/f ノイズ」という不思議な音

以前、科学者たちはこの「ざらつき」の高さの変化を聴覚に例えて分析しました。
すると、**「1/f ノイズ（ワン・エフ・ノイズ）」**という不思議な音が聞こえてきたのです。

• 1/f ノイズとは？
  • 音楽や自然の音（川のせせらぎ、風の音、心拍など）によくあるリズムです。
  • 「低い音（ゆっくりした変化）」と「高い音（速い変化）」が、ある決まったバランスで混ざり合っています。
  • これまで、この「ざらつき」の成長は**「永遠に落ち着かない（定常的ではない）」**と考えられていました。つまり、雪だるまがいつまで経っても大きくなり続け、永遠に「成長途中」の状態だと思われていたのです。

3. この論文の発見：「実は、小さな箱の中では落ち着いている！」

著者たちは、**「もし、その砂山を『小さな箱』の中で育てたらどうなるか？」**と考えました。

• 大きな箱（無限の大きさ）の場合：
  砂山は無限に大きくなり、いつまで経っても「成長中」の状態が続きます。これは「定常的（落ち着いている）」とは言えません。
• 小さな箱（有限の大きさ）の場合：
  箱の壁にぶつかるまで成長すると、もうこれ以上大きくはなりません。そこで、**「高さは一定の範囲で揺れ動くだけ」という「定常状態（落ち着き）」**に達します。

今回の重要な発見：
「小さなシステム（箱）の中では、実は**『定常状態』**に達している！」
つまり、このざらつきは「永遠に成長し続ける不安定なもの」ではなく、「ある一定の範囲で揺れている安定した状態」だったのです。

4. 2 つの証拠：「記憶」と「音」

著者たちは、この「落ち着き」を証明するために 2 つのテストを行いました。

1. 「記憶」のテスト（自己相関関数）：

   • 「今、砂山の高さが高いなら、10 分後にはどうなっている？」という**「過去の記憶」**を調べました。
   • 結果：「ある一定の時間（箱の大きさによる）が過ぎれば、過去の記憶は消え去る」ことがわかりました。これは、システムが「定常的」である証拠です。
   • また、この記憶の消え方は単純な指数関数ではなく、**「複雑な曲線」**を描いていました。

2. 「音」のテスト（パワースペクトル）：

   • 先ほどの「1/f ノイズ」の音（パワースペクトル）を詳しく聞きました。
   • 結果：低い音（ゆっくりした変化）の部分で、**「ある一定の音量以下には下がらない」**という「下限（カットオフ）」が見つかりました。
   • これは、**「システムが有限の大きさを持っている」**ことを示すサインです。無限に大きければ、この下限は存在しないはずだからです。

5. 結論：ウィーナー・チンヒンの定理が使える！

物理学には**「ウィーナー・チンヒンの定理」**という有名なルールがあります。

「『過去の記憶（相関）』と『音の成分（スペクトル）』は、実は表裏一体の関係にある」

これまで、このルールは「定常的なもの」にしか適用できないとされていました。
しかし、今回の研究では、**「小さな箱の中で成長するざらつきは、実は定常的だった」と証明したため、「このルール（ウィーナー・チンヒンの定理）が、この現象にも当てはまる！」**と結論づけました。

まとめ：何が変わったの？

• 以前の考え方： 「成長する表面のざらつきは、永遠に不安定で、定常的ではない（だから複雑な理論が必要だ）。」
• 今回の発見： 「小さなシステム（箱）の中では、実は落ち着いて（定常的に）揺れている。だから、昔からある『記憶と音の関係』というシンプルな法則が使える！」

イメージ：
まるで、**「無限に広がる海（無限のシステム）」は常に荒れていて落ち着かないけれど、「小さなプール（有限のシステム）」**に入れた水は、波が壁に反射して一定のリズムで揺れているだけ、という状態です。

この発見は、実験で観測される「ざらつき」のデータを、よりシンプルで確実な物理法則を使って理解できる道を開いたのです。

技術概要

以下は、提示された論文「Stationary 1/f α noise in discrete models of the Kardar-Parisi-Zhang class（カーダール・パリジ・チャンクラス離散モデルにおける定常な 1/f α ノイズ）」の技術的サマリーです。

1. 研究の背景と問題設定

非平衡系における「1/f ノイズ（ピンクノイズ）」は、電子デバイスから生物学的活動まで多様なシステムで観測される普遍的な現象ですが、その一般的な説明は未だ確立されていません。特に、カーダール・パリジ・チャン（KPZ）普遍性クラスに属する成長する粗い界面（Kinetic Rough Interfaces）において、界面の高さ変動が示す 1/f ノイズの性質について、以下の対立する見解が存在していました。

• 従来の見解（非定常性）: 大規模な系（$L \to \infty$）や有限の観測時間において、パワースペクトル密度（PSD）に下限の遮断周波数（ローカットオフ）が存在しないため、信号は「非定常（non-stationary）」であり、ウィーナー・ヒンチンの定理（相関関数とパワースペクトルの関係）が適用できないとされてきました。また、時間相関関数が時間並進対称性を満たさず、老化（aging）挙動を示すと報告されていました。
• 本研究の問い: 有限サイズの系において、十分に長い観測時間（相関時間よりも十分長い時間）を設けることで、系が「定常状態（stationary state）」に到達しうるか？そして、その場合、ウィーナー・ヒンチンの定理は有効か？

2. 手法とモデル

本研究では、KPZ 普遍性クラスに属する 4 つの代表的な離散モデル（セル・オートマトン）をシミュレーションし、固定されたサイトにおける高さ変動 $\delta h_L(t)$ の時間的相関とパワースペクトルを独立して解析しました。

• 対象モデル:
  1. Kim-Kosterlitz (KK) モデル
  2. 弾性衝突堆積（Ballistic Deposition: BD）モデル
  3. エッチング（Etching: EM）モデル（堆積バージョン）
  4. シングルステップ（Single-Step: SS）モデル
• シミュレーション条件:
  • 1 次元格子（サイズ $L$）、周期的境界条件。
  • 時間単位はモノレイヤー単位（$L$ 回の堆積試行を 1 ステップとする）。
  • 過渡状態（transient）を除外し、十分長い観測時間 $T$（$T \gg L^z$、ここで $z$ は動的指数）までシミュレーションを行い、定常状態への到達を確認。
• 解析手法:
  • 二時間自己相関関数: $C(\tau, L) = \langle \delta h_L(t) \delta h_L(t+\tau) \rangle$ を計算。
  • パワースペクトル密度（PSD）: フーリエ変換を用いて $S(f, L)$ を計算。
  • 有限サイズスケーリング（FSS）: 異なる系サイズ $L$ に対するデータの崩れ（data collapse）を確認し、臨界指数を決定。

3. 主要な結果

A. 相関関数の性質と定常性の確認

• 非指数関数的減衰: 相関関数は指数関数的に減衰せず、$C(\tau) \sim \tau^{1/z}$（$\tau \ll L^z$）のべき乗則に従うことが確認されました。
• 相関時間の発散: 相関時間 $\tau_c$ は系サイズとともに発散し、$\tau_c \sim L^z$ となります。
• 定常性の立証: 観測時間 $T$ が相関時間 $L^z$ よりも十分に長い場合、相関関数は時間差 $\tau$ のみの関数となり、時間並進対称性が回復します。これは、有限サイズの系において**広義定常（wide-sense stationary）**な状態が達成可能であることを示しています。
• スケーリング則: 正規化された相関関数は $C(\tau, L)/C(0, L) = F_C(\tau/L^z)$ とスケーリングします。

B. パワースペクトルと 1/f ノイズ

• 下限カットオフの存在: 定常状態では、PSD は相関時間の逆数に相当する周波数 $f_c \sim L^{-z}$ 以下で一定の値（フラット）を示します。これは、有限サイズの系では相関時間が有限であるため、低周波領域でパワーが飽和することを意味します。
• 1/f^α スケーリング: 非自明な周波数領域（$L^{-z} \ll f \ll 1/2$）において、PSD は $S(f) \sim f^{-\alpha}$ のべき乗則に従います。
• スペクトル指数の決定:
  • 数値解析により、動的指数 $z \approx 1.51(2)$（理論値 $3/2$ と一致）、スペクトル指数 $\alpha \approx 5/3$ を得ました。
  • 理論的関係式 $\alpha = 1 + 2\chi/z$（KPZ 普遍性クラスでは $\chi=1/2, z=3/2$）から導かれる $\alpha = 5/3$ と完全に一致します。
• モデル普遍性: KK, BD, EM, SS のすべてのモデルで同じスケーリング挙動が観測され、これが KPZ 普遍性クラスに固有の性質であることが確認されました。

4. 理論的枠組みとウィーナー・ヒンチンの定理

本研究の核心的な貢献は、スケーリング理論的な議論によって以下の関係を確立した点にあります。

• ウィーナー・ヒンチンの定理の適用: 系が定常状態にある場合、PSD は相関関数のフーリエ変換として導出可能です。
• スケーリング関数の対応: 相関関数のスケーリング関数 $F_C(u)$ と PSD のスケーリング関数 $F_S(u)$ の間には、フーリエ変換の性質に基づき、$F_S(u) \sim u^{-(1+1/z)}$（$u \gg 1$）という関係が導かれます。これにより、観測された $\alpha = 1 + 1/z = 5/3$ という値が理論的に裏付けられました。
• 低周波領域の振る舞い: 低周波（$f \to 0$）におけるパワーは $S(f \to 0, L) \sim L^{\alpha z}$ と系サイズに依存して増加し、総パワーは $P(L) \sim L^{(\alpha-1)z}$ とスケーリングします。

5. 結論と意義

• 定常状態の到達可能性: 無限大の系では非定常（老化）挙動を示す KPZ 過程ですが、有限サイズの系において十分な観測時間を設ければ、広義定常状態に到達できることを実証しました。
• 理論と実験の架け橋: 従来の研究では「非定常」と見なされていた KPZ 界面の高さ変動が、実際には定常的な 1/f^α ノイズ（$\alpha=5/3$）を生成しており、ウィーナー・ヒンチンの定理が有効であることを示しました。
• 実験的限界の解明: 実験において低周波カットオフが観測されにくい理由は、相関時間 $L^z$ が系サイズとともに発散するため、有限の観測時間では定常状態に到達する前に観測が終了してしまう（非定常と誤認される）ことに起因すると結論付けました。

本研究は、KPZ 普遍性クラスにおけるノイズの性質を再定義し、非平衡統計力学における定常性とスケーリング理論の整合性を明確にした重要な業績です。