Dimensional reduction for optical beams with thermal nonlocal nonlinearity

技術要約：熱的非局所非線形性を伴う光ビームの次元削減

問題提起
熱屈折効果に由来する光学的な非線形性は、熱拡散と吸収によって支配される長距離の材料応答をもたらす。勾配屈折率（GRIN）媒体において、この非局所性は、局所的な3次非線形媒質では不安定化して崩壊してしまうはずの渦ビームやリング・ソリトンを安定化させるなど、モード間の相互作用を根本的に変化させる。しかし、熱拡散方程式に結合した完全な空間非線形シュレディンガー方程式（NLSE）から、これらのダイナミクスを正確にモデリングすることは計算負荷が高い。局所的なカー非線形性に対しては既存の簡約化された次元モデル（非多項式シュレディンガー方程式、NPSEなど）が存在するが、これらは熱的応答の非局所的な性質を考慮できていない。著者らは、モード次数および非局所性の程度への依存性を捉えつつ、熱的非局所媒体におけるビームダイナミクスの扱いやすい簡約化された記述の必要性に取り組んでいる。

手法
著者らは、NPSEの枠組みを熱的に媒介される非局所非線形性を組み込む形へと拡張している。そのアプローチには以下のステップが含まれる：

1. 支配方程式： 系は、光電場エンベロープ $\Phi$ に対するパラキシャル空間NLSEと、温度誘起の屈折率変化に対する定常状態の熱方程式を結合させることでモデル化される。弱い吸収と特定の境界条件の仮定の下で、熱方程式は2次元の遮蔽ポアソン方程式へと簡約化される。結果として得られる光学伝搬方程式には、修正ベッセル関数カーネル（$K_0$）との畳み込みによって定義される非局所演算子が含まれる。
2. 変分アンザッツ： 次元削減を行うために、著者らは損失のない系のラグランジアン定式化に基づく変分法を用いる。横方向の電場プロファイルは、任意の方位角次数 $m$（ここで $S=|m|$）およびゼロの径方向指数 $p=0$ を持つ単一のラグエール・ガウス（LG）モードであると仮定するが、この定式化は一般的な $p$ も許容する。極めて重要な点として、アンザッツには径方向の位相チャープ項 $b(z)r^2/2$ が含まれており、著者らはこれが非局所媒体における回折を適切に記述するために必要であると特定している。
3. 有効モデルの導出： ラグランジアン密度を横平面にわたって積分することにより、著者らは有効な1次元作用を導出する。このプロセスにより、モード次数と変分的な横幅 $\sigma(z)$ に依存する非局所的な寄与を含む、明示的な解析表現が得られる。非局所項は、強度プロファイルのハンケル変換を用いて評価され、基本ガウスモード（$p=0, m=0$）に対する閉じた形式の表現、および高次モードに対する一般的な積分形式を与える。
4. 運動方程式： 有効作用から導かれるオイラー＝ラグランジュ方程式は、以下の結合系をもたらす：
   • $\sigma$ に依存する非線形項を含む、軸方向振幅 $A(z)$ の発展方程式。
   • 局所的な回折、GRINポテンシャル、および非局所的な熱圧のバランスに基づいて、横幅 $\sigma(z)$ を決定する代数方程式。
   • 径方向位相チャープ $b(z)$ に関する方程式。

主な結果

• 解析表現： 本論文は、熱的相互作用の強さを定量化する非局所積分 $I_{pS}(\sigma)$ の明示的な表現を導出している。基本ガウスモードについては、これは指数積分関数 $Ei$ を用いて表現される。
• 径方向チャープの役割： 本研究は、標準的なNPSEアンザッツ（径方向位相項を持たないもの）は非局所媒体には不十分であることを示している。ビームの幅の呼吸振動（ブリージング振動）を正しく捉えるためには、チャープパラメータ $b(z)$ の包含が不可欠である。
• 定常状態： 著者らは、ビーム幅を支配する有効ポテンシャル $V_{eff}(\sigma)$ を分析している。退色（デフォーカス）領域において、非局所性長 $L_{eff}$ を増加させると、ポテンシャルの最小値がより大きな幅の方へシフトし、結果としてより幅の広い定常ビームが得られることを示している。集束（フォーカス）領域では、非局所性は小さな幅においてより深いポテンシャルの井戸を形成する。
• 数値検証： 簡約モデルは、パラキシャル回折方程式の完全な2次元スプリットステップフーリエシミュレーションによって検証されている。
  • 短距離非局所性の領域（$L_{eff} = 0.5$）では、変分アンザッツは様々な渦次数に対して厳密な定常径方向強度プロファイルと極めて良好に一致している。
  • 強非局所性の領域（$L_{eff} = 5.0$）では、厳密なビームプロファイルは仮定されたLG関数形から逸脱する。その結果、変分モデルは定常強度プロファイルの細部において不一致を示す。
  • これらの静的なプロファイルの不一致にもかかわらず、簡約モデルは、吸収によるパワーの減衰や横幅の減衰を伴う呼吸振動を含む、ビームの動的な発展を正確に捉えている（$L_{eff} = 40$ の場合でも）。

意義と主張
本論文は、GRINファイバーにおける熱的に媒介される相互作用のための、一般化された枠組みを提供することを主張している。その主要な貢献は、モード次数および横幅への依存性を明示的に組み込んだ、簡約化された1次元モデルの導出である。著者らは、モデルが高度に非局所的な領域（ビーム形状がLGアンザッツから歪む領域）において厳密な横方向プロファイルを完全に再現することはできないものの、安定性や伝搬発展といった全体的なビームダイナミクスを理解するための強力なツールであり続けることを強調している。本研究は、遮蔽ポアソン型の非局所性を伴うシステム（これは熱媒体だけでなく、液晶、電気透析効果、プラズマ相互作用にも関連する）への変分次元削減手法の適用範囲を広げるものである。著者らは、現在のモデルは連続波および保存系（吸収は現象論的に追加されている）に限定されており、将来的な一般化においては、マルチモード場、飽和応答、およびパルス系における時間依存の熱緩和に対処できる可能性があると述べている。