Eigenvalue-by-Eigenvalue Comparison of a Sierra--Rodríguez-Laguna-Type Spectrum with the Riemann Zeros

技術的要約：シエラ–ロドリゲス＝ラグナ型スペクトルとリーマンの零点との固有値ごとの比較

問題の定義
本論文は、リーマンゼータ関数の非自明な零点（$1/2 + i\gamma_n$）が自己共役作用素の固有値に対応するというヒルベルト–ポリア予想の妥当性を調査している。具体的には、著者らはベリー–キーティングのハミルトニアン（$H_{BK} = xp$）の修正モデルであるシエラ–ロドリゲス＝ラグナ（SR）モデルを分析している。$H_{BK}$はリーマンの零点の平均密度を再現するが、連続スペクトルを持つ。$x = \ell_x$に硬い壁を持つSRモデル（$H_{SR} = x(p + \ell_p^2/p)$）は、離散スペクトルをもたらす。中心となる問いは、SRハミルトニアンの個々の固有値（$E_n$）がリーマンの零点の始点（$\gamma_n$）と一対一で一致させることができるか、また、もし可能であれば、$n \to \infty$における差 $E_n - \gamma_n$ の漸近挙動はどうなるか、ということである。

手法
著者らは、数値計算と厳密な漸近解析の組み合わせを用いている：

1. 数値的スペクトルの生成： SRハミルトニアンの自己共役拡張は、固有値に関する超越方程式 $\text{Re}[e^{-i\theta/2}K_{1/2+iE/2}(2\pi)] = 0$ を導く（ここで $K_\nu$ は変形ベッセル関数であり、$\theta$ は自由パラメータである）。著者らは、$\theta = 1.417\pi$ においてこの方程式を数値的に解き、最初の606個の固有値（$E_1, \dots, E_{606}$）を得た。
2. 漸近解析（最急降下法）： $E_n$ の閉形式の予測を導出するために、著者らは、大きな $E$ に対する変形ベッセル関数 $K_{1/2+iE/2}(2\pi)$ の積分表示に対して最急降下法（鞍点法）を適用した。彼らは支配的な鞍点を特定し、主要項の位相挙動を導出した。
3. 厳密な漸近式による検証： 隠れた位相項の欠落がないことを確実にするため、著者らは、自身の鞍点法による導出を、大きな虚数次数を持つ $K_\nu(x)$ に関する既知の厳密な漸近公式（Temme, Balogh, Booker, Strömbergsson, Then, および Tseng を引用）と照合・検証した。
4. 量子化とランベルトW関数による解法： 量子数 $n$ を用いて $E_n$ を明示的に表すために、量子化条件をランベルトW関数を用いて解いた。
5. リーマン–ヴァン・マンゴルト公式との比較： 導出した $E_n$ をリーマン・ヴァン・マンゴルトの計数関数 $N(T)$ と比較した。彼らは、$N(T)$ の跳躍不連続点における揺らぎ項 $S(T)$ の扱い、特に用いられている慣習を批判的に分析した。
6. 経験的検証： 理論的な予測を、A. Odlyzko によって表形式にまとめられた最初の105個のリーマンの零点に対してテストした。

主な貢献と結果

• 漸近的一致： 著者らは、最急降下法によって導かれた変形ベッセル関数の主要項の位相が、厳密な漸近結果と正確に一致することを実証した。これにより、ベッセル関数自体がSRスペクトルとリーマンの零点の間の数値的な不一致の原因である可能性を排除した。
• 位相シフトの解決： 標準的な予測（$N(\gamma_n) = n$ と仮定するもの）と数値データの間に不一致が特定された。著者らは、これを計数関数の跳躍における $S(T)$ の取り扱いに起因すると突き止めた。彼らは、跳躍における $S(T)$ の「中間点慣習」を用いる標準的なアプローチが、有効な量子数を $1/2$ シフトさせることを主張している。
  • 標準的なアプローチ： $\gamma_n \sim g(n - 7/8)$
  • 修正されたアプローチ： $\gamma_n \sim g(n - 11/8)$
• 閉形式の予測： 修正された計数関数を用いて、著者らは $E_n - \gamma_n$ の閉形式の予測式を導出した：
  $$E_n - \gamma_n \sim \frac{\theta + 15\pi/4}{1 + W(n/e)}$$
  ここで $W$ はランベルトW関数であり、これは差が $O(1/\ln n)$ で減衰することを予測している。
• 数値的一致： 最初の105個の零点を比較した際、修正された予測（図2の黒線）は、測定されたデータに対して 0.03% 以内の精度で一致した。対照的に、修正されていない「標準的な」予測は、約24%の偏差を示した。データの残差の散らつきは、$S(T)$ に関するセルベルグの中央極限定理と一致している。
• アーティファクトの特定： 著者らは、$E \approx 947.5$ 付近にある偽の根（元は $E_{607}$ とラベル付けされていたもの）を特定し、除去した。これは、2つの根が合体した数値的なアーティファクトであり、異常な準位間隔比によって確認された。

意義
本論文は、SRハミルトニアンの固有値がリーマンの零点を追跡しているかという問いに対し、満足のいく解決策を提供したと主張している。その意義は以下の点にある：

1. 精度： SRモデルが、単に統計的にではなく、予測可能な誤差項（$n \to \infty$ で消失する）を伴って、固有値ごとにリーマンの零点を再現することを確立した。
2. 理論的明晰さ： 以前観察されていた不一致は、ハミルトニアンの物理学やベッセル関数の性質によるものではなく、リーマンの計数関数 $N(T)$ の定義における微妙な慣習によるものであることを明らかにした。
3. モデルの妥当性： 本研究は、SRハミルトニアンがヒルベルト–ポリア予想の背後にある作用素の候補としての地位を強化するものであり、その主要な漸近挙動は堅牢であり、残りの不一致はすべて理解可能で修正可能であることを示している。

著者らは控えめなトーンを維持しており、一致は強力であり解決策も満足できるものであるものの、その議論は $S(T)$ の特定の慣習に依存しており、これらは文献における標準的なものであるが、さらなる精査の余地があることを述べている。彼らはリーマン仮説を証明したと主張しているのではなく、むしろ特定の量子モデルと既知の零点との間の精密な対応関係を確立したのである。