WGQED1D: a microscopic N-atom waveguide-QED toolkit for single-photon scattering and weak-measurement experiments

WGQED1D に関する技術要約

問題提起
原子アンサンブルルの中を通る弱光の伝搬は、量子光学における基本的な問題であり、電磁誘導透過（EIT）、スローライト伝搬、および負の事後選択的滞在時間などの現象に関連している。マクロな線形応答理論（例：ビア・ランベルト吸収およびクラマース・クローニッヒ分散）は確立されているが、$N$ 個の個々の原子による微視的なダイナミクスと、雲全体の統合された光学応答との間の溝を埋めるには、慎重な入力・出力処理が必要である。著者らは、標準的な公式を単純に実装すると、しばしばビア・ランベルト限界に対してオーダー1の程度の定量的な不一致が生じることを指摘している。これらの不一致は、導波路結合の双方向構造に関する微妙な慣習の選択、および光学的厚さ（OD）と原子あたりの結合効率との間の較正に起因している。

手法
本論文は、1次元双方向導波路における $N$ 個の2準位原子による単一光子散乱の微視的シミュレーションのために設計された、検証済みの軽量なPythonツールキット WGQED1D を紹介する。その手法は以下のコアコンポーネントに基づいている：

1. 物理モデル： システムは単一励起部分空間内でモデル化されており、これは単一光子入力と低飽和によって正当化される。ダイナミクスは、$N$ 個の複素振幅 ($b_j$) に作用する標準的な非エルミート有効ハミルトニアン ($H_{\text{eff}}$) によって支配される。このハミルトニアンには以下が含まれる：
   • 遅延フォトンの交換（$e^{ik_0|z_j-z_k|}$）によって媒介される集団的な双極子－双極子相互作用。
   • 自由空間への対角的な損失。
   • 入力ガウス波束を表す駆動項。
2. 入力－出力関係： ツールキットは、雲の境界に対する原子の位相蓄積を考慮した特定の入力－出力関係を用いて、出力電場 ($\phi_{\text{out}}$) を再構成する。
3. 決定的な較正： 本手法の主要な貢献は、原子あたりの導波路結合効率（$\beta = \Gamma_{1D}/\Gamma$）と目標とする光学的厚さ（$OD$）の間の特定の較正の特定と検証である。著者らは、$\beta = OD/(2N)$ である場合に、自由空間のビア・ランベルト限界と定量的一致が得られることを確立した。この2倍の因子は、原子が前方向と後方向の両方に等しい確率で放出するという双方向性の性質から生じるものであり、これは一方向（カイラル）な理想化とは対照的である。
4. 数値実装： ツールキットは、明示的なルンゲ＝クッタ 4(5) 法を用いた決定論的な発展を用いる。これは $O(N^2)$ のメモリおよびステップあたりの演算量で動作し、密な線形代数を利用する。また、Wiseman–Steinberg の弱値公式を用いて事後選択的な原子励起時間 ($\tau_T$) を抽出するためのスペクトル解析ツールも備えている。

主な貢献

• WGQED1D ツールキット： 微視的モデル、入力－出力再構成、およびスペクトル弱値解析を実装した、再現可能なオープンソースのPythonモジュール（約400行）。NumPy、SciPy、Matplotlib といった標準的な科学ライブラリのみに依存している。
• 微視的較正： $\beta = OD/(2N)$ という関係の明示的な導出と数値的検証を行い、微視的シミュレーションとマクロなビア・ランベルト予測との間の系統的な定量的不一致を解決した。
• 検証インフラストラクチャ： 本論文は、3つの異なる検証を文書化している：
  1. 光学的厚さ $OD \in [0.25, 3]$ における解析的なパルス平均ビア・ランベルト透過との一致（誤差 $<0.1\%$）。
  2. 事後選択的な励起時間に関する Wiseman–Steinberg の弱値公式の検証。
  3. 負の原子励起時間に関する Angulo et al. (2026) の実験データとの直接比較。

結果

• ビア・ランベルトの一致： $N=500$ 個の原子を用いたシミュレーションは、テストされた光学的厚さの範囲全体にわたって、解析的なビア・ランベルト透過率を 0.1% 未満の精度で再現しており、希薄限界における $\beta = OD/(2N)$ の較正の妥当性を確認している。
• 弱値の抽出： ツールキットは、解析的な弱値公式と一致する事後選択的な励起時間 ($\tau_T$) を正常に抽出する。不一致は 10% 以内であり、これらは有限のスペクトル分解能および高光学的厚さにおける集団的な有限 $N$ 補正に起因する。
• 実験との比較： 冷たいルビジウム雲における負の事後選択的滞在時間に関する Angulo et al. のデータと比較した際、ツールキットは 6 つのデータポイントのうち 3 つを $1\sigma$ 誤差範囲内で再現している。残りのポイント（主に $OD=4$ の場合）については、シミュレーションは定性的な傾向を追跡しており、場合によっては線形応答の解析公式よりも実験的な測定値に近い結果を示している。
• パフォーマンス： ツールキットは標準的なハードウェア上で効率的に動作し、現実的なアンサンブル（$N \sim 200–1000$）に対して数秒から数分で収束した結果を生成する。

意義と主張
本論文は、WGQED1D を、光学的に厚い原子雲を通る単一光子透過における弱測定、異常分散、群遅延、および事後選択的観測量を研究するための「再現可能な数値プラットフォーム」として位置づけている。著者らは、このツールキットが生のパフォーマンスよりも数値的な信頼性と検証可能性を優先しており、支配方程式への直接的な行単位の一致を提供することを強調している。

本研究の意義は、双方向の導波路結合の微妙な差異を正しく扱う検証済みのツールを提供することにあり、これにより、微視的なシミュレーションをマクロな限界や最近の実験的観測と定量的に一致させることが可能になる。著者らは、より大きなシステムサイズ（$N \gtrsim 10^4$）を必要とするユーザーに対しては、アルゴリズム構造を疎行列代数または GPU バックエンドに拡張可能であると明記しているが、これらは現在実装されていない。ツールキットは、オープンな科学的探究を促進するために MIT ライセンスの下で公開されている。