All-Electron Molecular Tunnel Ionization Based on the Weak-Field Asymptotic Theory in the Integral Representation

Imam S. Wahyutama,Denawakage D. Jayasinghe,Francois Mauger,Kenneth Lopata,Kenneth J. Schafer
2024-08-18
Abstract:Tunnel ionization (TI) underlies many important ultrafast processes, such as high-harmonic generation and strong-field ionization. Among the existing theories for TI, many-electron weak-field asymptotic theory (ME-WFAT) is by design capable of accurately treating many-electron effects in TI. An earlier version of ME-WFAT relied on an accurate representation of the asymptotic tail of the orbitals, which hindered its implementation in Gaussian-basis-set-based quantum chemistry programs. In this work, we reformulate ME-WFAT in the integral representation, which makes the quality of the asymptotic tail much less critical, hence greatly facilitating its implementation in standard quantum chemistry packages. The integral reformulation introduced here is therefore much more robust when applied to molecules with arbitrary geometry. We present several case studies, among which is the CO molecule where some earlier theories disagree with experiments. Here, we find that ME-WFAT produces the largest ionization probability when the field points from C to O, as experiments suggest. An attractive feature of ME-WFAT is that it can be used with various types of multi-electron methods whether of density functional [Phys. Rev. A 106, 052211 (2022)] or multiconfiguration types, which facilitates tunnel ionization calculation in systems exhibiting a strong multireference character, such as transition-metal-containing molecules.
Chemical Physics,Atomic and Molecular Clusters,Computational Physics
What problem does this paper attempt to address?
This paper aims to solve some key problems in the theory of tunnel ionization (TI), especially its application in multi - electron systems. Specifically, the goals of the paper are: 1. **Improve the multi - electron weak - field asymptotic theory (ME - WFAT)**: The existing ME - WFAT depends on the accurate representation of the asymptotic tails of orbital wave functions, which is difficult to implement in quantum chemistry programs using Gaussian basis sets. By reformulating ME - WFAT into an integral form (IR), the paper greatly reduces the requirement for the accuracy of the asymptotic tails, thus making the theory easier to implement in standard quantum chemistry software packages. 2. **Improve computational efficiency and applicability**: By introducing the integral form, the method proposed in the paper is more robust when dealing with molecules of arbitrary geometric structures and can be more widely applied to different molecular systems, especially those with strong multi - reference characteristics, such as molecules containing transition metals. 3. **Verify the validity of the theory**: Through multiple case studies, the paper shows the consistency between the new method and experimental results in calculating the ionization probability of CO molecules. In particular, when the electric field points from the carbon atom to the oxygen atom, the ionization probability predicted by ME - WFAT is the largest, which is in line with experimental observations. 4. **Expand the application range of the theory**: The paper points out that the new IR ME - WFAT method can be combined with various multi - electron methods, whether based on density functional theory (DFT) or multi - configuration methods, which provides more flexibility and accuracy for tunnel ionization calculations of complex systems. In summary, the main contribution of this paper lies in the development of a more efficient and more general multi - electron tunnel ionization theoretical method, enabling it to be applied in a wider range of molecular systems and better handle the effects of multi - electron interactions.