The accurate prediction of ionization potentials (IPs) is central to understanding molecular reactivity, redox behavior, and spectroscopic properties. While vertical IPs can be accessed directly from electronic excitations at fixed nuclear geometries, the computation of adiabatic IPs requires nuclear gradients of the ionized states, posing a major theoretical and computational challenge, especially within correlated frameworks. Among the most promising approaches for IP calculations is the many-body Green's function GW method, which provides a balanced compromise between accuracy and computational efficiency. Furthermore, it is applicable to both finite and extended systems. Recent work has established formal connections between GW and coupled-cluster doubles (CCD) theory, leading to the first derivation of analytic GW nuclear gradients via a unitary CCD framework. In this work, we present an alternative, fully analytic formulation of GW nuclear gradients based on a modified version of the traditional equation-of-motion CCD formalism, enabling the inclusion of missing correlation effects in the traditional CCD methods.
Kitsaras et al. (Tue,) studied this question.