Improving convergence of experimental and computed vertical ionization energies using the ionization potential equation-of-motion coupled-cluster with singles and doubles method

Comparison of statistically evaluated experimental vertical ionization energies (IEs) for 53 medium-sized molecules (6-34 atoms) with ionization potential equation-of-motion coupled-cluster with singles and doubles (IP-EOMCCSD) computations shows that discrepancies between computed and experimental results can be accounted for with a combination of experimental and theoretical contributions. Discrepancies can be minimized by extrapolating computations to the complete basis set limit and correcting for vibrational zero-point energy (ZPE) while comparing with experimental IEs calculated as the intensity-weighted mean band position to account for band asymmetries. This procedure reduced the average discrepancy for ethylene, (E)-2-butene, 2,5-dihydrofuran, and pyrrole from 0.25 to 0.05 eV. Agreement between reported vertical IEs and computations without either making adjustments as described in this paper or using complete simulation of the ionization spectrum should be considered fortuitous. The comparisons made in this work show that estimates of vertical and adiabatic IE made using IP-EOMCCSD extrapolated to the complete basis set limit and corrected for vibrational ZPE can be used with reasonable confidence when experimental values are not available.