[论文解读] Theoretical description of the ionization processes with a discrete basis set representation of the electronic continuum

Mateusz S. Szczygieł, Michał Lesiuk+1 more|arXiv (Cornell University)|Sep 17, 2019

Laser-Matter Interactions and Applications参考文献 6被引用 7

一句话总结

该论文提出了一种新颖的离散基组方法，利用高斯基组中的库仑型轨道来近似电子连续态波函数，实现了对光电子截面和光电子角分布的精确计算。该方法在高达60 eV的光子能量范围内，对H、He和H₂的实验数据表现出极佳的一致性，证明了其与标准量子化学软件的兼容性以及在长度表象和速度表象下的规范不变性。

ABSTRACT

In this work we present a new method of approximating the continuum wavefunctions with a discrete basis set. This method aims to be at least compatible with other well known methods of the electronic structure theory to describe processes in the electronic continuum. As an example of the application we present the results of the calculations of the photoionization cross sections for the hydrogen atom and hydrogen molecule and for the helium atom. We also obtained the photoelectron angular distribution for the hydrogen molecule. The agreement of our results with the experimental data is very good in the energy range from the ionization threshold up to photon energies of about 60 eV.

研究动机与目标

• 开发一种实用且计算兼容的电子连续态表示方法，采用离散基组。
• 通过使用高斯型轨道，确保与现有量子化学软件的兼容性。
• 准确计算原子和分子体系的光致电离截面和光电子角分布（MFPAD）。
• 通过与实验数据和既有的理论方法对比，验证该方法在长度表象和速度表象下的有效性。
• 为未来在分子光致电离中扩展至时间依赖和关联电子动力学研究奠定基础。

提出的方法

• 该方法通过从库仑波函数导出的离散基组来近似连续态波函数，适配于高斯型轨道（GTO）框架。
• 连续轨道通过一种变换生成，将精确的库仑波解映射为适合分子哈密顿量的有限基表示。
• 采用长度表象和速度表象计算偶极矩阵元，以确保规范不变性并保持与实验可观测量的一致性。
• 光致电离截面通过频率依赖的极化率 α⁻(ω) 计算，关系式为 σ = (4πω/c) Im[α⁻(ω)]。
• 光电子角分布（MFPAD）通过展开式 ∂σ/∂Ωₖ = Σₗₘ Aₗₘ Yₗₘ(θₖ, φₖ) 计算，其中系数 Aₗₘ 由连续轨道振幅导出。
• 该方法采用完整组态相互作用（CI）波函数表示初始态，并通过插值离散角度点生成平滑的MFPAD曲线，以与实验结果对比。

实验结果

研究问题

• RQ1 能否以一种与标准量子化学软件兼容的方式，利用高斯型轨道构建电子连续态的离散基组表示？
• RQ2 该方法在从电离阈值到60 eV的能量范围内，对氢原子、氦原子和氢分子的光致电离截面预测精度如何？
• RQ3 该方法在长度表象和速度表象下，对实验光电子角分布（MFPAD）的重现程度如何？
• RQ4 该方法在计算可观测量（尤其是偶极矩阵元）时是否保持规范不变性？
• RQ5 该方法能否捕捉实验中观察到的更高阶角动量贡献（如d波特征）在光电子分布中的表现？

主要发现

• 该方法在从电离阈值到60 eV的光子能量范围内，对氢原子、氦原子和氢分子的实验光致电离截面均表现出极佳的一致性。
• 对于H₂，长度表象和速度表象下的计算MFPAD与Hikosaka和Eland（2003）的实验数据高度吻合，尤其在21.2、26.9、40.8、44.0、46.0、60.0和76.0 eV能量点，仅在40.8 eV的平行极化情况下存在微小偏差。
• 在40.8 eV平行极化情况下观测到的d波形状贡献与Semenov和Cherepkov（2003）的先前理论结果一致，验证了该方法捕捉高阶角动量特征的能力。
• 该方法对H₂的2pπᵤ和2sσ₊g离子态（Ito等，2000）的MFPAD实验数据的匹配，证实了其在复杂角分布中的高精度。
• 长度表象和速度表象的计算结果几乎完全一致，证明了方法实现的规范不变性和鲁棒性。
• 该方法成功再现了MFPAD的圆柱对称性和镜像对称性，插值曲线通过对称操作完整反映了全角度依赖性。