Negative superfluid density and spatial instabilities in driven superconductors

We consider excitation of Higgs modes via the modulation of the BCS coupling within the Migdal-Eliashberg-Keldysh theory of time-dependent superconductivity. Despite the presence of phonons, which break integrability, we observe Higgs amplitude oscillations reminiscent of the integrable case. The dynamics of quasiparticles follows from the effective Bogolyubov-de Gennes equations, which represent a Floquet problem for the Bogoliubov quasiparticles. We find that when the Floquet-Bogoliubov bands overlap, the homogeneous solution formally leads to a negative superfluid density, which is no longer proportional to the amplitude of the order parameter. This result indicates an instability, which we explore using spatially-resolved BdG equations. Spontaneous appearance of spatial inhomogeneities in the order parameter is observed and they first occur when the superfluid density becomes unphysical. We conclude that the homogeneous solution to time-dependent superconductivity is generally unstable and breaks up into a complicated spatial landscape via an avalanche of topological excitations.