Phase and coherence retrieval from near- and far-field intensities

技術要約：近視野および遠視野の強度からの位相およびコヒーレンス再構成

問題提起
複雑な光学場のコヒーレンス特性を定量化することは、高出力レーザーアレイから量子コヒーレントシステムに至るまで、幅広いアプリケーションにおいて極めて重要である。ゲルシュベルク・サクソン（GS）フレームワークのような反復位相回復アルゴリズムは、複雑な場の再構成における標準的な手法であるが、これらは通常、完全コヒーレントな場を想定している。大規模なレーザーアレイ、X線／電子ビーム、あるいは構造化された量子状態を含む多くの重要な物理系は、部分的な空間コヒーレンスしか示さない。このような場合、場は単純な複素振幅では記述できず、一次空間コヒーレンスを表す4次元の相互強度関数 $J(r, r') = \langle E^*(r)E(r') \rangle$ を必要とする。既存の再構成手法は、多くの場合、制限的な仮定（例：均質なコヒーレンス、既知のモード基底、または低ランク密度行列）に依存しているか、あるいは2つの平面（近視野と遠視野）のみの強度測定から完全な相互強度を回復することに失敗する。

手法
著者らは、2つのフーリエ共役平面における強度測定から、部分コヒーレント場の完全な相互強度を再構成するための新しいパラダイムを提案している。このアプローチは、近視野と遠視野のドメイン間で相互強度を反復的に伝搬および投影するようにGSアルゴリズムを適応させたものである。2つの補完的な実装が開発されている：

1. テンソルGSアルゴリズム： この手法は、4次元の相互強度テンソルに対して交互投影を適用することで、相互強度を直接再構成する。プロセスには以下が含まれる：

   • 4次元フーリエ変換を用いた平面間の相互強度の伝搬。
   • 測定された近視野および遠視野の強度に対して、相互強度の対角要素を投影することによるターゲット強度の強制。
   • 基礎となるコヒーレンス構造を維持しつつ、エルミート対称性と物理的一貫性（$|\tilde{g}_1(k, k')| \leq 1$ の保証）を強制しながら、非対角要素の更新。
   • 収束するまでの反復。

2. モンテカルロ（MC）GSバリアント： 4次元フーリエ変換に伴う計算コストを削減するため、この手法は部分コヒーレント場を、$P$ 個の独立したランダムな複素実現体のアンサンブルとして表現する。

   • アンサンブル内の場は、標準的な2次元フーリエ変換を用いて平面間を伝搬される。
   • 各ステップにおいて、アンサンブル平均の強度がターゲット強度と一致するように、場に共通の関数を乗算する。
   • 相互強度は、$P$ 個の実現体のアンサンブル平均によって推定される。
   • 精度は $1/\sqrt{P}$ に従ってスケールし、制御された近似が可能である。計算複雑性は、テンソルGS法が $M^2 \log(M)$ であるのに対し、MC GS法は $PM \log(M)$ （ここで $M$ は空間モード数）となる。

主な貢献

• 一般的な再構成フレームワーク： 本論文は、定義済みのモード基底、低ランクパラメータ化、または複数の検出平面を必要とせずに、わずか2平面の強度測定から完全な空間相互強度を再構成する一般的な手法を提示している。
• デュアル・アルゴリズム的アプローチ： 高精度なテンソルGS法と計算効率の高いMC GSバリアントの両方の導入により、異なるシステムサイズや精度要件に応じた柔軟性を提供している。
• 不均質性の処理： これらの手法は、空間的に不均質な部分コヒーレント場を扱うように特別に設計されており、従来の逆フーリエ変換（IFT）アプローチが失敗する領域に対応している。

結果
著者らは、数値シミュレーションと実験的実証を通じて両方の手法を検証した：

• 数値シミュレーション：
  • 線形アレイ： ガウス型減衰コヒーレンスを持つ最大 $N=600$ のガウスビームの線形アレイを用いたシミュレーションにおいて、IFT法は短いコヒーレンス長（$\xi \ll Nd$）に対してのみ正確にコヒーレンスを回復するが、テンソルGS法およびMC GS法の両方は、相互強度が強く不均質になる領域（$\xi > Nd$）においても精度を維持することを示した。
  • リングアレイ： リングアレイ（$N=50$）のシミュレーションでは、並進不変性の仮定の崩壊によりIFT法が質的に失敗することを示した。対照的に、提案された手法の両方は、相互強度の正しい構造を正確に回復した。テンソルGS法は高忠実度の再構成を提供し、MC GS法は計算効率の高い近似を提供した。
• 実験的実証：
  • テンソルGS法は、結合とスタッガード・カイラリティが制御された130個の結合レーザーからなる三角形アレイに適用された。
  • このシステムは、非単調な相互コヒーレンス減衰と非自明な位相（スタッガード・アンチボルテックスおよびボルテックス／アンチボルテックスの等価重ね合わせ状態）を生成した。
  • 再構成は理論と極めて良好に一致し、期待される位相巻き（例：$2\pi/3$）およびコヒーレンス振幅のガウス型減衰を再現した。
  • RMS位相誤差は、スタッガード・アンチボルテックス状態で $2\pi/250$、等価重ね合わせ状態で $2\pi/65$ と測定された。

意義と主張
本論文の結果は、このアプローチが大規模な部分コヒーレント光学系の特性評価に実用的であることを示唆していると著者らは主張している。これらの手法は、コヒーレンス構造の事前知識や複数の検出平面を必要としない点で、既存のトモグラフィー的手法を補完するものである。著者らは、2つのフーリエ共渉平面からの相互強度の再構成における理論的な一意性は依然として未解決の問題であるとしつつも、数値的および実験的な結果は、強く不均質なコヒーレンス構造に対して物理的に意味のある解への安定した収束を示すことを強調している。

本研究は、自由空間伝搬を組み込むことで非フーリエ共渉平面からの場を再構成することや、相互強度再構成と密度行列再構成の形式的な類似性を踏まえ、量子トモグラフィーにこのフレームワークを適用することなどの、将来的な拡張を示唆している。しかし、論文は、特定の困難な不均質系における位相とコヒーレンスの再構成能力を実証することに焦点を当て、任意の光学場へと主張を広げすぎないという控えめな立場を維持している。