An entropic barrier to equilibration in the discrete nonlinear Schrödinger equation: Rigorous bounds for the microcanonical Brownian motion

技術要約：離散非線形シュレディンガー方程式におけるエントロピー障壁による平衡化

問題提起
本論文は、離散非線形シュレディンガー方程式（DNLSE）の緩和時間スケールを調査し、特に高エネルギー密度においてシステムが「凝縮相」へと入る現象に焦点を当てている。エネルギー密度 $\varepsilon$ が臨界閾値 $\varepsilon_c = 2$ を超えると、システムが劇的な平衡化の遅延を示すことが知られている。このとき、エネルギーは長期間にわたって単一の振動子に集中する（凝縮体）。数値シミュレーションでは、平衡化時間スケールがシステムサイズ $N$ に対して急速に発散することが示唆されているが、平衡化時間 $\tau_{eq}$ に関する厳密な解析的境界はこれまで欠如していた。著者らは、ハミルトニアン動力学をマイクロカノニカル面上のブラウン運動に置き換えることで、$\tau_{eq}$ の厳密な上界および下界を提供することを目的としている。

手法
著者らは、固定された全質量 $Q=N$ とエネルギー $E=N\varepsilon$ によって定義されるマイクロカノニカル面 $M_{N,\varepsilon}$ 上のシステムをモデル化している。動力学は、反射境界条件を持つ多様体上のブラウン運動によって近似される。平衡化時間は、スペクトルギャップ $\lambda_1$（ラプラシアン $-\Delta_M$ の最初の非ゼロのノイマン固有値）を介して $\tau_{eq} = 2/\lambda_1$ という関係によって支配される。

分析は、エネルギー密度 $\varepsilon$ に基づいて2つのフェーズに分けられる：

1. 流体相 ($\varepsilon < 2$): システムは非局在状態にある。
2. 凝縮相 ($\varepsilon > 2$): システムは局在を示し、一つのサイトがマクロな量のエネルギーを保持する。

$\lambda_1$ を抑えるために、著者らは2つの異なる数学的手法を用いている：

• $\lambda_1$ の上界（$\tau_{eq}$ の下界）: レイリー商を用いて、$\lambda_1$ を上から推定するための試行関数を構成する。流体相では、単一変数の試行関数で十分である。凝縮相では、凝縮体がサイト間を遷移することを捉えるために、2変数試行関数を利用し、実質的にシステムを構成空間内の「首（ネック）」によって接続された二重井戸問題としてモデル化する。
• $\lambda_1$ の下界（$\tau_{eq}$ の上界）: $\lambda_1$ を下から抑えるために、多様体上でボホナー公式を積分した幾何学的テクニックを適用する。これには、マイクロカノニカル面の境界面の凹性を補償する特定のベクトル場の構成が必要となる。$\varepsilon \ge 2$ の場合、この議論は「深いキャップ」（一つの座標が非常に大きい領域）に局所化され、その後、ポアンカレ不等式とオーバーラップ推定を用いたグルーイング（接合）手順を通じてグローバルな多様体へと拡張される。

主要な貢献と結果

1. マイクロカノニカル面の幾何学
本論文は、$M_{N,\varepsilon}$ の幾何学を厳密に特徴付けている。$\varepsilon > 2$ のとき、この曲面は、特定のサイトが凝縮体を保持する $N$ 個の異なる「チャンバー（室）」であり、それらがペアごとに狭い「チューブ」によって接続されている。これらのチューブの最も狭い部分である「ネック」は、エントロピー的なボトルネックを表している。ネックの面積は、チャンバーの面積よりも指数関数的に小さく、スケールは $e^{-b(\varepsilon)\sqrt{N}}$ となる。

2. $\lambda_1$ の上界（平衡化時間の低界）

• 流体相 ($\varepsilon < 2$): 著者らは、$\tau_{eq} = O(1)$ であることを証明している。具体的には、単一サイトの周辺密度に関連するシュトゥルム＝リウヴィル問題から導かれる通り、$\lambda_1$ は $\varepsilon$ に依存する定数によって下から抑えられる。
• 凝縮相 ($\varepsilon > 2$): 平衡化時間は $N$ に対して伸長指数関数的（stretched exponential）であることが示される。$\lambda_1$ の上界は、2つのチャンバー間のネックを横切って符号が変化する試行関数を用いて導出される。結果は以下の通りである：
  $$\lambda_1 \le \exp\left( -(\sqrt{2}-1)\sqrt{\varepsilon-2}\sqrt{N} + O(\ln N) \right)$$
  これは、$\tau_{eq} \sim \exp(\Theta(\sqrt{N}))$ であることを意味する。

3. $\lambda_1$ の下界（平衡化時間の上界）
著者らは、ボホナー–レイリー基準を用いて、$\lambda_1$ の厳密な下界（したがって $\tau_{eq}$ の上界）を確立している。

• $\varepsilon < 2$ の場合、グローバルなベクトル場が $\lambda_1$ に対する一様な下界を与える。
• $\varepsilon > 2$ の場合、境界は凝縮体が存在する深いキャップに議論を局所化し、局所的なスペクトルギャップが大きいことを証明した後、これらの局所的な推定値をグルーイングすることで得られる。この結果は、伸長指数関数的なスケーリングを裏付けている：
  $$\lambda_1 \ge c_\varepsilon N^{-C_\varepsilon} \exp\left( -\sqrt{N(\varepsilon-2)} - O(N^{1/4}) \right)$$
  したがって、$\varepsilon > 2$ のとき：
  $$(\sqrt{2}-1)\sqrt{\varepsilon-2} \le \liminf_{N\to\infty} \frac{\ln \tau_{eq}}{\sqrt{N}} \le \limsup_{N\to\infty} \frac{\ln \tau_{eq}}{\sqrt{N}} \le \sqrt{\varepsilon-2}$$

意義と主張
本論文は、サイズとともに発散する、DNLSEマイクロカノニカル面上の動力学の緩和時間に関する初の厳密な境界を提供したと主張している。

• エントロピー的起源: 著者らは、凝縮相における遅い動力学をエントロピー障壁に帰している。システムは、構成空間内で凝縮体をサイト間で移動させるために、低体積の「ネック」を通過しなければならない。このプロセスは、エネルギー障壁ではなく、マイクロカノニカル面の幾何学によって支配されている。
• 他のモデルとの区別: 本結果は、凝縮体の移動が多項式時間スケールで行われる他の確率粒子系（ゼロ範囲過程など）とDNLSEを区別するものである。ここでは、質量とエネルギーの両方の保存がボトルネックを生み出し、伸長指数関数的な時間スケールをもたらしている。
• ハミルトニアン動力学との比較: 著者らは、結果がブラウン運動に対して導出されたものであるが、バス（熱浴）と結合したハミルトニアンDNLSEの数値的研究とも、「奇妙な一致」があることを指摘している。数値研究もまた、同様の $\sqrt{N}$ 依存性を持つブリーザーの指数関数的な寿命を示している。しかし、著者らはメカニズムが異なること（断熱不変量 vs エントロピー障壁）を強調しており、これらの厳密な境界を元のハミルトニアン動力学へ転用できるかどうかは依然として未解決の問題であるとしている。
• 方法論的影響: 本論文は、高次元多様体上の幾何学的解析（ボホナー–レイリー議論）と変分法を組み合わせることで、エルゴード性の破れのタイムスケールを制御できることを示しており、この手法は他の確率的問題にも適用可能であると示唆している。

著者らは、自身の境界の鋭さについては謙虚な姿勢を保っており、$\lambda_1$ の下界（$\tau_{eq}$ の上界）は鋭くない可能性が高いと考えている一方で、$\lambda_1$ の上界（$\tau_{eq}$ の下界）はおそらく鋭いと考えている。また、前係数はブラウン運動に選択された特定の計量に依存するが、指数部分は平衡測度の特性であることを述べている。