Divergences in classical and quantum linear response and equation of motion formulations

Calculating molecular properties using quantum devices can be performed through the quantum linear response (qLR) or, equivalently, the quantum equation of motion (qEOM) formulations. Different parameterizations of qLR and qEOM are available, namely naïve, projected, self-consistent, and state-transfer. In the naïve and projected parameterizations, the metric is not the identity, and we show that it depends on redundant orbital rotations. This dependency may lead to divergences in the excitation energies for certain choices of the redundant orbital rotation parameters in an idealized noiseless setting. Furthermore, this leads to a significant variance when calculations include statistical noise from finite quantum sampling.