[論文レビュー] Theory of strong-field ionization of aligned CO2

Mahmoud Abu-samha, Lars Bojer Madsen|ArXiv.org|Jun 4, 2009

Advanced Chemical Physics Studies被引用数 4

ひとこと要約

本稿では、一電子近似を用いたab initio時間発展シュレーディンガー方程式（TDSE）法を用いて、整列したCO2分子の強場電離をモデル化する。グリッドベースの数値解法により、完全な分子ポテンシャルおよび励起電子状態を組み込むことで、実験的電離収率と前例のない一致を達成した。この結果、電離ダイナミクスは電子密度ではなく中間的共鳴状態によって支配されていることが示され、分子トンネル理論や強場近似といった半解析的理論が失敗する理由が解明された。

ABSTRACT

A theoretical framework for studying strong-field ionization of aligned molecules is presented, and alignment-dependent ionization yields are computed for CO2. Our calculations are in unprecedented agreement with recent experiments. We find that the ionization process is affected by intermediate resonance states, and the alignment-dependent ionization yields do not follow the electron density of the initial states. The theory explains the breakdown of semi-analytical theories, like the molecular tunneling theory and strong-field approximation, were excited electronic structure is neglected.

研究の動機と目的

• H2およびH2+を越える多原子分子の整列状態における強場電離の理論的枠組みを構築すること。
• CO2における整列依存電離の実験的測定値と、半解析的理論（例：分子トンネル理論）との乖離を解消すること。
• 特にトンネル電離と多光子電離の遷移領域において、励起電子状態が強場電離ダイナミクスに与える影響を調査すること。
• 標準的な量子化学コードと互換性を持つグリッドベースのab initio手法を提供し、複雑な分子系への広範な応用を可能とすること。

提案手法

• 速度ゲージにおける有限差分グリッド法を用いて、時間発展シュレーディンガー方程式（TDSE）を数値的に解く。
• GAMESSで計算されたDFT-LDAポテンシャルの角運動量統合を用いて、分子ポテンシャルを部分波展開する。
• 時間発展波動関数をグリッドベースのスペクトル法で射影し、電離収率および運動量分布を抽出する。
• 10周期、800 nmのレーザーパルスを用い、ピーク強度が5.6×10^13および1.1×10^14 W/cm²の条件で実験条件を模擬する。
• 理論的電離収率を実験的整列分布関数と畳み込み処理して、測定データと一致させる。
• pavicic:243001の実験データと比較することで、異なる強度における角度依存電離収率をもとに結果の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 なぜCO2における実験的整列依存電離収率が、初期状態の電子密度に基づく予測から逸脱するのか？
• RQ2 励起電子状態は、CO2のような多原子分子における強場電離ダイナミクスにどのように寄与するのか？
• RQ3 分子トンネル理論や強場近似といった半解析的モデルが、整列CO2の電離を記述する際にどの程度失敗するのか？
• RQ4 励起状態ダイナミクスを含むab initio TDSEシミュレーションは、トンネル電離と多光子電離の遷移領域における実験的電離収率を再現できるか？
• RQ5 中間的共鳴状態は、CO2における電離収率の角度依存性を決定づける役割を果たすのか？

主な発見

• 1.1×10^14 W/cm²における理論的電離収率は、整列分布関数と畳み込まれた後、実験測定値と良好に一致し、モデルの正確性を確認した。
• 1.1×10^14 W/cm²における電離収率は、0°で7.35×10^-4、45°で2.58×10^-3、80°で2.43×10^-4、90°で5.31×10^-5であり、電子密度では予測できない強い角度依存性を示した。
• トンネル理論は実験データから大きく外れた電離収率を予測しており、本領域ではその理論の失敗を示した。
• HOMOの電子がHOMO-1（πu, 2p）軌道に遷移し、全電離収率の約10%を占めることが判明し、励起状態の重要な寄与を示した。
• ピーク強度におけるエネルギースペクトルは、明確に束縛状態を特定しており、πu（HOMO-1）、πg（HOMO）、πu（LUMO）、およびπu（C 3p）が識別された。これは電子構造の正確な記述を裏付けた。
• 本手法により、電離ダイナミクスは初期状態の電子密度ではなく、中間的共鳴状態および励起状態の多重準位によって支配されていることが明らかになった。これにより、半解析的理論の破綻が説明された。