Collective modes in superconductors including Coulomb repulsion

We numerically study the collective excitations present in isotropic superconductors including a screened Coulomb interaction. By varying the screening strength, we analyze its impact on the system. We use a formulation of the effective phonon-mediated interaction between electrons that depends on the energy transfer between particles, rather than being a constant in a small energy shell around the Fermi edge. This justifies considering also rather large attractive interactions. We compute the system's Green's functions using the iterated equations of motion (iEoM) approach, which ultimately enables a quantitative analysis of collective excitations. For weak couplings, we identify the well-known amplitude (Higgs) mode at the two-particle continuum's lower edge and the phase (Anderson-Bogoliubov) mode at $\omega = 0$ for a neutral system, which shifts to higher energies as the Coulomb interaction is switched on. As the phononic coupling is increased, the Higgs mode separates from the continuum, and additional phase and amplitude modes appear, persisting even in the presence Coulomb interactions.