Equation generator for equation-of-motion coupled cluster assisted by computer algebra system

We present an equation generator algorithm that utilizes second-quantized operators in normal order with respect to a correlated or non-correlated reference and the corresponding Wick theorem. The algorithm proposed here, written with Mathematica, enables the generation of non-redundant strings of second-quantized operators that, after classification, are directly assigned to many-body term quantities used to construct the many-body Hamiltonian. We demonstrate the capabilities of the algorithm by computing the coupled-cluster amplitude equations and various blocks of the equation-of-motion many-body Hamiltonian. A comprehensive description of this four-step algorithm is provided alongside concrete examples.