Quantitative comparison of TDDFT-calculated HHG yields in ring-shaped organic molecules

Stephanie N. Armond,Kyle A. Hamer,Ravi Bhardwaj,Francois Mauger,Kenneth Lopata,Kenneth J. Schafer,Mette B. Gaarde
2024-09-18
Abstract:We compare the high-harmonic-generation (HHG) yield driven by a mid-infrared laser in three organic ring-shaped molecules, calculated using time-dependent density-functional theory (TDDFT). We average the yield over the relative orientation of the molecules and the linearly-polarized, 1825 nm driving laser pulse in order to compare to experimental spectra obtained by Alharbi et al., Phys. Rev. A 92, 041801 (2015). We find that the raw TDDFT-calculated HHG yield in cyclohexane (CHA) is strongly overestimated compared to those of benzene and cyclohexene, and that this can be attributed to unphysically large contributions from CHA orbitals lying well below the highest-occupied molecular orbital. We show that implementing a simple orbital-resolved scaling factor, which corrects the yield of the tunneling ionization contribution to the first step in the HHG process, leads to much better comparisons with experimental results. Our results are encouraging for the use of TDDFT in systematic computations of HHG in large molecules.
Atomic Physics
What problem does this paper attempt to address?
### Problems Addressed by the Paper The paper primarily addresses the following issues: 1. **Quantitative Comparison of High Harmonic Generation**: Using Time-Dependent Density Functional Theory (TDDFT), the high harmonic generation (HHG) yields of three cyclic organic molecules under mid-infrared laser fields were calculated and quantitatively compared. 2. **Comparison with Experimental Results**: The calculated HHG yields were compared with the experimental results published by Alharbi et al. in 2015 to validate the effectiveness of the theoretical calculations. 3. **Correction of Ionization Potentials**: It was found that the original TDDFT calculations for cyclohexane (CHA) significantly overestimated the HHG yields due to the underestimation of the ionization potentials of low-energy orbitals, especially HOMO-2. A simple orbital-resolved scaling factor based on the ADK ratio was proposed to correct this issue, thereby better matching the experimental results. ### Main Findings - **Differences in HHG Yields of Different Molecules**: There are significant differences in the HHG yields of benzene (BNZ), cyclohexene (CHE), and cyclohexane (CHA), with the original calculations for CHA being particularly high. - **Orbital Contribution Analysis**: By analyzing the contributions of different orbitals to the HHG yields, it was found that the ionization potential of the HOMO-2 orbital in CHA was underestimated, leading to an overestimation of its HHG yield. - **Correction Method**: A simple ADK ratio correction factor was introduced to adjust the ionization steps in the HHG process, making the calculated results more consistent with the experimental data, especially in suppressing the yield of CHA. ### Significance This study not only validates the effectiveness of TDDFT in systematically calculating HHG for large molecules but also proposes a simple and effective correction method that helps improve computational accuracy and theoretical predictive capability.