Chebyshev self-imaging from inverse-sampled angular spectra

技術要約：逆サンプリング角スペクトルによるチェビシェフ自己像形成

問題提起
タルボット効果に代表される従来の自己像形成は、傍作用領域における空間周波数の二次位相蓄積に依存している。この領域では、横方向に周期的な場は、固定された距離（タルボット長 $z_T$）およびその整数倍において自己再構成を行う。この再生距離は、調和位相（$\propto m^2$）が可換であるため、場の帯域幅には依存しない。すなわち、すべてのモードが単一の周期の整数倍において同時に再位相化する。これまでの研究では、ガウス和を評価したり整数を干渉法によって因数分解したりするために、この二次位相カーネルを利用してきたが、その根底にある伝搬算術は固定された二次的なものである。著者らは、自己像形成の算術自体を、単に固定されたカーネルに対して数論的構造を読み取るのではなく、数論的構造を直接的に露呈または計算するように設計できるのではないかと問いかけている。

手法
著者らは、波動の角スペクトルをエンジニアリングすることで、自己像形成の根本的な算術を変更する手法を提案している。空間周波数の標準的な線形間隔ではなく、モード指数 $m$ の逆数冪で角スペクトルをサンプリングする。

場は $M$ 個の横平面波の重ね合わせとしてモデル化される：
$$E(x,z) = \sum_{m=1}^{M} a_m e^{ik_x(m)x} e^{-ik_x^2(m)z/2k_o}$$
ここで、横方向の空間周波数は $k_x(m) = k_L m^{s/2}$ と定義される。

• 従来の場合 ($s > 0$): $s=2$（タルボット）または $s=1$（モンゴメリー効果）の場合、再位相化の周期は単位区間内にネストし、$M$ に依存しない普遍的な再生距離を持つ。
• 逆サンプリングの場合 ($s = -r$): 著者らは $s = -r$（$r$ は正の整数）とし、その結果 $k_x(m) = k_L m^{-r/2}$ と設定した。この構成では、モード $m$ の位相は $\exp(-2\pi i \zeta / m^r)$ として進化する（ここで $\zeta = z/z_T$）。

この逆サンプリング下では、モード $m$ の再位相化周期は $m^r$ に比例する。厳密な自己像形成（$E(x,z) = E(x,0)$）は、伝搬距離 $\zeta$ がすべての個別周期 $\{m^r\}_{m=1}^M$ の公倍数である場合にのみ発生する。したがって、最初の厳密な再生距離は、モード指数の $r$ 乗の最小公倍数（LCM）によって決定される：
$$z_{rev}(M) = z_T \cdot \text{lcm}(1, 2, \dots, M)^r$$

主要な貢献と結果
本論文は、このエンジニアリングが、再生距離を固定された定数から、整数の乗法的構造によって支配されるプログラム可能な量へと変貌させることを示している。

1. チェビシェフ再生と素数冪階段関数:
   再生座標は $z_{rev}/z_T = e^{r\psi(M)}$ と表される。ここで $\psi(M)$ は第二チェビシェフ関数（$\psi(M) = \ln(\text{lcm}(1, \dots, M))$）である。著者らは、モードが追加されても再生距離は一定に保たれるが、モード指数 $M$ が新しい素数冪 $p^a$ に達したときに、距離が $p^r$ の因子分ジャンプすることを観察している。これにより、「素数冪階段」が伝搬距離の中に形成され、自由空間の回折を整数の乗法的構造の物理的な読み出し装置へと変えている。

2. 約数再生:
   厳密な再生の間、場は「約数再生」を示す。中間的な整数平面 $\zeta$ において、インデックスが $\zeta$ を割り切るモードのサブセットが再位相化し、部分的な自己像を形成する。軸上の強度 $F(\zeta)$ は、整数 $M$ までの約数格子に対応する階層的なピークを表示する。最初の厳密な再生は、これら部分的なピークとは異なり、すべてのモードが同時に再位相化する最初の点として区別される。

3. 実験的実証:
   連続波レーザー（$\lambda = 635$ nm）と位相のみを制御する空間光変調器（SLM）を用い、逆サンプリング・スペクトル（$r=1$）をプログラムし、軸方向強度を測定した。

• 連続するスペクトル $S = \{2, \dots, M\}$ について、$M=3, 4, 5, 6$ のとき $\zeta = 6, 12, 60, 60$ で最初の再生を観察した。$M=5$ と $M=6$ の間のプラトーは、合成数（6）を追加しても LCM が変化しない一方で、素数（5）を追加すると5倍の増加をもたらすことを裏付けている。
• 疎なスペクトル（例：$S=\{4,5,6,7\}$）は、フルセット $\{1, \dots, 7\}$ と同じ再生距離を達成できることが示され、これは、再生が単なる帯域幅のサイズではなく、占有された帯域幅の素数冪含有量に基づいてプログラム可能であることを示している。
• $r=2$（逆二乗サンプリング）の測定も提示されており、$p^2$ のジャンプを示し、二乗された約数格子の構造を確認している。

意義
本論文は、波動の再生距離が固有の固定された特性ではなく、プログラム可能な算術量であることを確立している。モード位相率を成形することで、著者らは、タルボット効果の加法的可換性から、整数の乗法的構造へと支配原理を転換させている。

その意義は主に以下の3点にある：

• 基礎物理学: 自己像形成の算術を反転させ、再生を普遍的な定数ではなく、モード周期のLCMへと遅延させる。
• 数論的読み出し: 伝搬距離自体が、モード指数に対する素数および素数冪のテストを実行する。発射から再生までの強度の「カーペット」は、整数の約数格子を物理的に描き出す。
• 量子論的アナロジー: $r=2$ のケース（$k_x \propto m^{-1}$）は、エネルギー・スペクトル（$E_n \propto -1/n^2$）に対応する。著者らは、これらの光学的な再生が、量子波動束の研究で通常用いられる近似的なテイラー展開の代わりに、厳密な算術限界を用いるリドバーグ波動束の古典的エミュレータとして機能することを示唆している。

著者らは、再生距離が帯域幅に対して指数関数的に増大するため（$z_{rev} \sim e^{rM}$）、自由空間でアクセス可能な $M$ には限界があるものの、この構成はプラットフォームに依存しないと述べている。分散型ファイバー光学を用いれば、$\sim 10^4$ までの整数のLCMを伝搬長にエンコードできる可能性があり、これによってアクセス可能な範囲を拡張できることを示唆している。本研究は、有限のスペクトル遮断を、単なる打ち切り誤差としてではなく、プログラム可能な算術レジスタとして再定義している。