Cluster perturbation theory. XI. Excitation-energy series using a variational excitation-energy function

Traditionally, excitation energies in coupled-cluster (CC) theory have been calculated by solving the CC Jacobian eigenvalue equation. However, based on our recent work [Jørgensen et al., Sci. Adv. 10, eadn3454 (2024)], we propose a reformulation of the calculation of excitation energies where excitation energies are determined as a conventional molecular property. To this end, we introduce an excitation-energy function that depends on the CC Jacobian and the right and left eigenvectors for the Jacobian eigenvalue problem. This excitation-energy function is variational with respect to the right and left eigenvectors but not with respect to the cluster amplitudes. Instead, the cluster amplitudes satisfy the cluster-amplitude equations, and we set up an excitation-energy Lagrangian by adding to the excitation-energy function the cluster-amplitude equations with an undetermined multiplier for each cluster-amplitude constraint. The excitation-energy Lagrangian is variational in all its parameters. Based on the variational property of the Lagrangian, we have determined two quadratically convergent excitation-energy series: the total-order cluster-perturbation (tCP) and variational cluster-perturbation (vCP) excitation-energy series. Calculations of the excitation energies of three small molecules have shown that the vCP series is to be preferred over the tCP series. The test calculations have been carried out for CPS(D) expansions [targeting the CC singles-and-doubles (CCSD) wave function from the CC singles wave function] and the CPSD(T) expansion [targeting the CC singles-doubles-triples (CCSDT) wave function from the CCSD wave function]. For the S(D) and SD(T) orbital excitation space calculations, we obtain in the second vCP iteration excitation energies with a mean deviation from CCSD excitation energies of about 0.04 eV for the S(D) orbital spaces, and for the SD(T) orbital space calculation, we obtain a mean deviation from the CCSDT excitation energies of 0.001 eV.