A counterexample to Fermi isospectral rigidity for two dimensional discrete periodic Schrödinger operators

この論文は、数値的検証を用いて、2 次元離散周期シュレーディンガー作用素においてフェルミ等スペクトル性が剛性を持たないことを示す非自明な実数値周期ポテンシャルの存在を証明し、2 次元におけるフェルミ多様体の既約性に関する 1990 年代の予想を反証したものである。

要約

この論文は、数学の難しい分野（量子力学や格子理論）における「常識」を覆す、とても面白い発見について書かれています。専門用語を避け、日常の例えを使って分かりやすく解説しましょう。

🎵 音楽と「隠れたリズム」の話

まず、この研究の舞台は**「格子（グリッド）」と呼ばれる、無限に広がる点の並びです。そこには「ポテンシャル（可能性の場）」という、目に見えない「地形」や「重み」が乗っています。これを「シュレーディンガー演算子」と呼びますが、イメージとしては「無限に続くピアノの鍵盤」や「広大な森の地形」**だと思ってください。

この「地形」がどうなっているかによって、電子（音）がどう動くかが決まります。

1. 「フェルミ・アイソスペクトラル」とは？

この論文の核心は**「フェルミ・アイソスペクトラル」という概念です。
これを「音の響き」**に例えてみましょう。

• ゼロ・ポテンシャル（0）: 何もない平坦な地面、あるいは何もない静かな部屋。ここでの「音の響き（スペクトル）」は単純で、規則正しいものです。
• あるポテンシャル（V）: 地面に小さな段差や、壁が置かれている状態。

通常、地面に段差（ポテンシャル）があれば、音の響き（スペクトル）は変わります。だから、「音が同じなら、地面は平坦でなければならない（段差はない）」というのが、これまでの**「常識（rigidity/剛性）」**でした。

しかし、この論文は**「実は、段差（ポテンシャル）があっても、特定の条件下では『音の響き』がゼロ（平坦）と全く同じになることがある！」**と証明しました。

2. 2 次元の「魔法の地形」

これまでの研究では、3 次元以上の世界では「音が同じなら、地形は必ず平坦（ゼロ）である」ということが証明されていました。つまり、3 次元以上では「音の響き」から地形を完全に特定できる（剛性がある）のです。

でも、2 次元（平面）の世界ではどうだろう？
これが、この論文で解決された疑問です。

著者たちは、**「3 行 5 列の小さなパターンが無限に繰り返される、2 次元の魔法の地形」**を見つけ出しました。

• この地形は、平坦ではありません（段差があります）。
• でも、不思議なことに、この地形の上を電子が動く時の「特定のエネルギーレベル（音の高さ）」における「音の響き」は、何もない平坦な地面と完全に一致するのです。

まるで、**「複雑な模様が描かれた壁紙を貼っても、部屋全体の『広がり感（音の響き）』が、何もない白い壁と全く同じになる」**ような現象です。

3. 昔の「予想」は間違っていた

1990 年代、有名な数学者たちは**「2 次元の世界でも、実在する（物理的に意味のある）地形なら、音が同じになるのは『平坦な地面』だけだ」**と予想していました（Gieseker-Knorrer-Trubowitz 予想）。

しかし、この論文は**「いいえ、違います！平坦じゃない地形でも、音が同じになる『魔法の地形』が存在します！」**と、その予想を否定しました。

4. どうやって証明したの？（コンピュータの力）

この「魔法の地形」を見つけるのは、人間の頭脳だけで計算するのは不可能なほど複雑な方程式です。

著者たちは、**「数値証明（Numerical Certification）」**という高度な技術を使いました。

• ステップ 1: コンピュータに「音が同じになるような地形」を探させました。すると、コンピュータは「たぶん、この数値の組み合わせが正解に近いよ」という答え（近似解）を見つけました。
• ステップ 2: ここが重要ですが、ただの近似では「もしかしたら誤差で間違っているかも？」という疑念が残ります。そこで、**「クラウツキの方法（Krawczyk's method）」**という、数学的に「間違いなく、この範囲に正解が存在する」と証明する厳密なアルゴリズムを使いました。
• 結果: 2 つの異なるコンピュータプログラム（Macaulay2 と Julia）で計算し直し、**「間違いなく、この地形が存在する」**と数学的に証明しました。

🌟 まとめ：この発見が意味すること

1. 2 次元の世界は特別: 3 次元以上では「音から地形がわかる」のが当たり前ですが、2 次元では「音が同じでも、地形は違う」ということがあり得ることがわかりました。
2. 予想の崩壊: 昔の天才数学者たちの「2 次元ではそんなことはあり得ない」という予想が、実は間違っていたことが証明されました。
3. 新しい可能性: この発見は、物質科学や量子力学において、**「同じ性質を持つのに、内部構造が全く違う新しい材料」**が存在する可能性を示唆しています。

つまり、この論文は**「2 次元の世界には、目には見えない『魔法の地形』が隠れていて、それによって物理的な性質が欺かれることがある」**という、驚くべき事実を、コンピュータの力を使って数学的に証明した物語なのです。

技術概要

論文概要

この論文は、2 次元離散周期シュレーディンガー演算子における「フェルミ同スペクトル（Fermi isospectrality）」の剛性（rigidity）に関する長年の未解決問題に対し、数値的証明（numerical certification）を用いて反例を構成し、否定答を与えることを目的としています。また、フェルミ多様体の既約性に関する 1990 年代の予想（Gieseker-Knörrer-Trubowitz 予想）も同時に反証しています。

1. 研究背景と問題設定

• 対象: 格子 $\mathbb{Z}^d$ 上の離散シュレーディンガー演算子 $H = \Delta + V$。ここで $\Delta$ は離散ラプラシアン、$V$ は $\Gamma$-周期ポテンシャルです。
• フェルミ同スペクトル性: 2 つの周期ポテンシャル $X, Y$ が、あるエネルギーレベル $\lambda_0$ においてフェルミ多様体（Fermi variety）$F_{\lambda_0}(X)$ と $F_{\lambda_0}(Y)$ が一致する場合、これらは「フェルミ同スペクトル」と呼ばれます。
• 既知の結果（$d \ge 3$）: 3 次元以上では、ゼロポテンシャルとフェルミ同スペクトルである実ポテンシャルはゼロポテンシャルのみであることが証明されています（Theorem 1.1）。
• 未解決の問題（$d = 2$）: 2 次元の場合、ゼロポテンシャルとフェルミ同スペクトルである実周期ポテンシャル $V$ は、必ず $V=0$ となるか？（Question 1）。
  • 直感的には 3 次元と同様に剛性が成り立つと考えられていましたが、証明されていませんでした。
• 関連する予想: Gieseker, Knörrer, Trubowitz (1990 年代) は、「2 次元における非定数実周期ポテンシャルに対して、フェルミ多様体はすべてのエネルギーレベルで既約である」と予想していました（Conjecture 1）。

2. 主要な結果

この論文は以下の 2 つの主要な定理を証明しています。

1. 定理 1.2（反例の存在）:
   2 次元において、ゼロポテンシャルとフェルミ同スペクトルである非自明な実周期ポテンシャルが存在します。具体的には、周期 $3\mathbb{Z} \times 5\mathbb{Z}$のポテンシャル$V$が見つかり、エネルギー$\lambda=0$において$F_0(V) = F_0(0)$ が成り立ちます。

   • これにより、2 次元におけるフェルミ同スペクトル剛性は成り立たないことが示されました（Question 1 への否定答）。

2. 定理 1.4（既約性の反証）:
   2 次元において、あるエネルギーレベルでフェルミ多様体が**可約（reducible）**となる非定数実周期ポテンシャルが存在します。

   • これにより、Gieseker-Knörrer-Trubowitz 予想が一般には偽であることが示されました。

3. 手法と数値的証明の枠組み

この研究の核心は、代数的な問題設定と、厳密な数値的証明（Numerical Certification）の組み合わせにあります。

3.1. 代数的定式化

• フロケ行列（Floquet Matrix）: 周期 $q_1 \times q_2$ のポテンシャル $V$ に対し、フロケ境界条件を課した行列 $D_V(z)$ を構成します。ここで $z = (z_1, z_2)$ はフロケ乗数です。
• 多項式方程式系への帰着:
  • フェルミ同スペクトル条件 $F_0(V) = F_0(0)$ は、行列式 $\det(D_V(z)) = \det(D_0(z))$ と同値です（Lemma 2.1）。
  • 周期 $(3, 5)$ の場合、$V$ の値を変数 $v_{i,j}$ とみなし、$\det(D_V(z)) - \det(D_0(z)) = 0$ を展開して係数を比較することで、14 個の多項式 $f_1(v) = \dots = f_{14}(v) = 0$ からなる連立方程式系 $F=0$ を構築しました。
  • この系は 15 変数（ポテンシャルの値）と 14 方程式ですが、平均値が 0 になるなどの条件を加えて解の存在を調べます。

3.2. 近似解の探索

• 多項式系の次数が高く、ベズーの定理に基づく解の総数は約 45 億に達するため、全解を数値的に求めるのは非現実的です。
• ホモトピー法（Homotopy Continuation）: HomotopyContinuation.jl パッケージを用い、ランダムな超平面切片と交差する解の軌跡を 1 つずつ追跡する「lazy-evaluation」アプローチを採用しました。
• 約 400 分の計算後、142,295 番目の軌跡で実数解（近似解 $V^*$）を発見しました。

3.3. 厳密な数値的証明（Krawczyk 法）

• 単なる数値計算では誤差の存在が証明として不十分です。そこで、Krawczyk 法（区間算術を用いた不動点定理の応用）を用いて、近似解 $V^*$ の近傍に厳密な実数解が存在することを証明しました。
• 手順:
  1. 近似解 $V^*$ を中心とする区間ボックス $B$ を定義します。
  2. ニュートン作用素 $N_F$ に対して、区間算術を用いて $N_F(B) \subset \text{interior}(B)$ となることを厳密に検証します。
  3. バナッハの不動点定理により、ボックス内に唯一の不動点（すなわち方程式の解）が存在し、かつそれが実数解であることを保証します。
• 実装: この証明を 2 つの独立した環境で実行し、結果の一致を確認しました。
  • Macaulay2: NumericalCertification パッケージを使用。
  • Julia: HomotopyContinuation.jl 内の Krawczyk 法実装を使用。
• 両方の環境で、非特異な実数解の存在が約 1 分〜5 分で証明されました。

4. 結論と意義

• 剛性の破れ: 2 次元離散周期シュレーディンガー演算子において、ゼロポテンシャルとフェルミ同スペクトルである非自明な実ポテンシャルが存在することが初めて証明されました。これは、高次元（$d \ge 3$）で成り立つ剛性が、2 次元では成り立たないことを示す決定的な反例です。
• フェルミ多様体の構造: 非定数実ポテンシャルであっても、フェルミ多様体が可約になる場合があることを示しました。これにより、1990 年代の Gieseker-Knörrer-Trubowitz 予想が一般には偽であることが判明しました。
• 計算機数学の応用: 代数的幾何学と数値的証明（Computer-Assisted Proof）を組み合わせることで、従来の解析的手法では到達できなかった複雑な非線形問題の解決に成功した点に大きな意義があります。特に、Macaulay2 と Julia の両方で独立した検証を行ったことは、結果の信頼性を高めています。

この研究は、逆問題（inverse problems）や可積分系、固体物理学におけるバンド構造の理解に対して、新たな視点と数学的厳密性を提供しています。