Temperature and mean axial momentum vs. laser intensity of electrons released from O$_2$ by an 800 nm ultrashort pulsed laser

技術要約：800 nm 超短パルスレーザーによって酸素から放出された電子の温度および平均軸方向運動量とレーザー強度の関係

問題提起
本論文は、空気中におけるレーザーフィラメントの電気力学的流体シミュレーションに不可欠な初期条件、具体的には熱平衡状態の電子温度（$T$）および平均ポスト光学軸方向運動量（$\langle p_{fz} \rangle$）を決定する必要性に対処している。先行研究では、フィラメントの導電率測定に基づき、電子密度と内部レーザー強度は決定されているが、これらパラメータと、集束された 800 nm 超短パルスレーザー（USPL）が酸素（$O_2$）と相互作用する際のピークレーザー強度（$I_0$）を結びつける堅牢なモデルが欠けていた。課題は、理論的な強電場近似（SFA）モデルと実験的な電子運動エネルギースペクトルの間を整合させることにあり、実験データは、再散乱や再結合といった複雑な電子・イオン相互作用に起因する特徴（低エネルギーのサージや高エネルギーのプラトーなど）を示しており、標準的な非断熱 SFA 理論では予測できないものとなっている。

手法
著者らは、Li らおよび Luo らによって定式化された、再定式化された非断熱強電場近似（SFA0）に基づいた半経験的モデルを提示している。このモデルは、時間依存性はあるが定常的なレーザーパルスエンベロープを仮定し、特定の瞬時電離時刻 $t_0$ に対する最も確率の高いトンネル経路のみを考慮して、電離過程を準古典的に扱う。

SFA0 理論と実験データの間の乖離を埋めるために、著者らは理論的な電子運動エネルギースペクトル $S_0(U)$ に対して、2 つの現象論的な修正を導入している。

1. SFA1 (シーリング・モデル): 実験データには存在しない低運動エネルギー領域（$U < U_c$）における電子集団の理論的なサージに対処するため、臨界エネルギー $U_c = 1.60$ eV でスペクトルに「天井（シーリング）」を設けている。これは、親イオンへの再結合による低エネルギー電子の抑制を考慮したものである。
2. SFA2 (エネルギー・リスケーリング): $4\pi$ ステラジアン平均データで観察されるより高い熱平衡温度に適合させるため、SFA1 スペクトルに対してエネルギー乗数 $\zeta = 2.067$ を適用している。これは、電子が 1 光周期内に親イオンからの弾性再散乱によって受けるエネルギーブーストを考慮したものである。

平均軸方向運動量 $\langle p_{fz} \rangle$ は、再散乱の古典的な運動学を用いて推定される。モデルは、電子が親イオンに戻る時刻（$t_1$）と、その後の角度 $\phi$ での弾性散乱後の運動量成分を計算する。全運動量は、異なる方向で記録された経験的な温度データ（$T_x, T_+, T_*$）を用いて角度分布に重み付けを行い、パルスプロファイルと立体角について積分することで導出される。

主な貢献

• 精緻化されたスペクトルモデル: 著者らは、標準的な SFA0 を拡張し、SFA1 および SFA2 を導入することで、さまざまな ケルディッシュ・パラメータ（$\gamma_0$）において 800 nm 線形偏光に曝露された $O_2$ の実験スペクトルに適合させることに成功した。
• 熱平衡温度 ($T$) と強度: 著者らは、ピーク強度 $I_0$（または $\gamma_0$）の関数としての熱平衡電子温度 $T$ の半経験的な関係を示す。モデルは、$x$ 方向放出から導かれた温度（$T_x$）、$2\pi$ 平面平均（$T_+$）、および $4\pi$ ステラジアン平均（$T_*$）を区別している。
• 平均軸方向運動量 ($\langle p_{fz} \rangle$) の推定: SFA2 熱平衡温度と古典的な再散乱運動学を組み合わせることで、$\langle p_{fz} \rangle$ を推定する手法が開発された。これには、レーザーパルス強度の時空間平均を考慮するための補正係数 $\alpha(\gamma_0)$ が含まれており、ピーク強度 $I_0$ を持つパルスの特性と、特定の強度における固有の特性を区別している。
• 偏光解析: モデルは線形偏光（$\epsilon=0$）と円偏光（$\epsilon=1$）を比較している。SFA0 は円偏光の方が高い温度を予測するが、線形偏光に関する実験データは、再散折効果によって線形の場合の方が同等または高い実効温度をもたらすことを示唆しており、マイクロ波放射における円偏光の期待される利点を打ち消している。

結果

• スペクトル適合: SFA2 モデル（$S_2(U)$）は、$0.82 \leq \gamma_0 \leq 1.30$ の範囲において、$4\pi$ ステラジアンで平均された実験スペクトルに対して良好な適合を示す。このモデルは、再散乱に起因するとされる実験で観察された高エネルギー・テールの「プラトー」をうまく再現している。
• 温度の傾向: SFA2 から得られる熱平衡温度 $T_2$ は、$\gamma_0 = 0.47$ まで有用な近似であることが示されたが、多光子電離領域が支配的となり SFA の仮定が崩れる $\gamma_0 > 1.40$ では精度が低下する。
• 運動量の計算: 計算された平均軸方向運動量 $\langle p_{fz} \rangle$ は、再散乱の影響を大きく受けることが示されている。$\gamma_0 = 1$ のとき、再散乱による寄与（$\langle p_{fz} \rangle_b$）は約 $0.80 \times 10^{-27}$ kg-m/s であり、非再散乱による寄与（$\langle p_{fz} \rangle_a$）は $1.66 \times 10^{-27}$ kg-m/s である。本論文は、$\gamma_0 < 1$ の場合、中性 $O_2$ の顕著な電離枯渇が $T_x$ データから導出される運動量推定に誤差をもたらす可能性があると指摘している。
• 補正係数: 解析により、ピーク強度における固有の運動量は、パルス平均の運動量よりも高いことが明らかになった。$\gamma_0 = 1$ のときの補正係数は $\alpha(\gamma_0) \approx 0.76$ であり、$\gamma_0 > 1$ ではさらに減少する。

意義および主張
本研究の主要な意義は、空気中のレーザーフィラメントの電気力学的流体シミュレーションのための初期条件（$T$ および $\langle p_{fz} \rangle$）を提供することであると著者らは主張している。これらの条件は、公開された実験データと、電子の再散乱および再結合を考慮した修正 SFA 理論を統合することで導出されている。

著者らは、自らの貢献を、強電場理論における根本的なブレイクスルーではなく、データ解析およびシミュレーションの初期化のためのツールとして控えめに位置づけている。また、以下の制限事項も認めている。

• 本モデルは $0.82 \leq \gamma_0 \leq 1.30$ で最も正確である。
• $\epsilon=1$（円偏光）の扱いは、スペクトルデータが不足しているため予備的なものである。
• クーロン集束、$O_2$ 内の 2 つの原子間の量子干渉、および高 $\gamma$ における原子共鳴などの効果は含まれていない。
• $\gamma_0 < 1$ の場合、現在の運動量推定では十分に捉えられていない中性分子の枯渇効果により、重大な誤差が生じる可能性がある。

本論文は、導電率測定を電子密度および温度に結びつけることを目的とした、USPL によって形成される空気フィラメントの物理的特性を決定するための、より広範な研究（[1] および [2] として参照）の補完的な構成要素として機能している。