Plasmon-enhanced Raman scattering can push single-molecule vibrational spectroscopy beyond a regime addressable by classical electrodynamics. We employ a quantum electrodynamics (QED) description of the coherent interaction of plasmons and molecular vibrations that reveal the emergence of nonlinearities in the inelastic response of the system. For realistic situations, we predict the onset of phonon-stimulated Raman scattering and a counterintuitive dependence of the anti-Stokes emission on the frequency of excitation. We further show that this QED framework opens a venue to analyze the correlations of photons emitted from a plasmonic cavity.Surface Enhanced Raman Scattering (SERS) is a spectroscopic technique in which the inelastic scattering from a molecule is increased by placing it in a hotspot of a plasmonic cavity, where the electric fields associated with the incident and the scattered photons are strongly enhanced (see the schematic in Fig. 1(a)). 1 The difference between the energy of those two photons provides a fingerprint of the molecule, i.e., detailed chemical information about its vibrational structure. Since the initial observation of Raman scattering from single molecules, 2,3 the use of a variety of plasmonic structures that act as effective optical nanoantennas over the last decades has allowed a tremendous advance of this molecular spectroscopy 4 . Metallic particles such as nanoshells 5,6 , nanorings 7 , nanorods 8 , nanowires 9 , or nano stars 10 , as well as plasmonic nanogap structures formed in particle dimers [11][12][13] , nanoparticle-on-a-mirror morphologies [14][15][16][17] , or nanoclusters 18 , are among the variety of structures that offer huge and controllable enhancements of the field intensity in their hotspots, boosting the inherently weak Raman scattering intensity, 1,19 and ultimately enabling the chemical identification and imaging of particular vibrational modes of a molecule with subnanometer resolution. 20 These results suggest that some experiments might have reached the regime where the quantum-mechanical nature of both the molecular vibrations and the plasmonic cavity emerges, 21 and call for an adequate theoretical description that goes beyond the classical treatment of the electric fields produced in plasmonic cavities. 1,22,23 In this work we address the underlying quantummechanical nature of Raman scattering processes by quantizing as bosonic excitations both the vibrations of the molecule and the electromagnetic field of a plasmonic cavity. The description of the vibrations through bosonic operators can be justified by considering the harmonic approximation to the energy landscape of the molecule along a generalized atomic coordinate ( Fig. 1(b)), such as the length of a molecular bond, e.g., C=O. 21,22 These vibrations interact with the cavity photons through a nonlinear Hamiltonian, reminiscent of that found in optomechanical systems. 24 In this description, the large enhancement of the Raman scattering from a molecule in the plasmonic cavity occurs thanks to ...