Equation-of-motion orbital-optimized coupled-cluster doubles method with the density-fitting approximation: An efficient implementation

Orbital-optimized coupled-cluster methods are very helpful for theoretical predictions of the molecular properties of challenging chemical systems, such as excited states. In this research, an efficient implementation of the equation-of-motion orbital-optimized coupled-cluster doubles method with the density-fitting (DF) approach, denoted by DF-EOM-OCCD, is presented. The computational cost of the DF-EOM-OCCD method for excitation energies is compared with that of the conventional EOM-OCCD method. Our results demonstrate that DF-EOM-OCCD excitation energies are dramatically accelerated compared to EOM-OCCD. There are almost 17-fold reductions for the ${\text{C}}_5{\text{H}}_{12}$ molecule in an aug-cc-pVTZ basis set with the RHF reference. This dramatic performance improvement comes from the reduced cost of integral transformation with the DF approach and the efficient evaluation of the particle-particle ladder (PPL) term, which is the most expensive term to evaluate. Further, our results show that the DF-EOM-OCCD approach is very helpful for the computation of excitation energies in open-shell molecular systems. Overall, we conclude that our new DF-EOM-OCCD implementation is very promising for the study of excited states in large-sized challenging chemical systems.