We investigate the interaction between light and molecular systems modelled as quantum emitters coupled to a multitude of vibrational modes via a Holstein-type interaction. We follow a quantum Langevin equations approach that allows for analytical derivations of absorption and fluorescence profiles of molecules driven by classical fields or coupled to quantized optical modes. We retrieve analytical expressions for the modification of the radiative emission branching ratio in the Purcell regime and for the asymmetric cavity transmission associated with dissipative cross-talk between upper and lower polaritons in the strong coupling regime. We also characterize the Förster resonance energy transfer process between donor-acceptor molecules mediated by the vacuum or by a cavity mode.PACS numbers: 05.60.Gg, 37.30.+i, 81.05.Fb Recent experimental progress [1] has shown that the Purcell enhancement of the zero-phonon line of a single molecule can strongly alter the branching ratio of spontaneous emission between the line of interest and additional Stokes lines thus turning the molecule into an ideal quantum emitter. At the mesoscopic level, experiments in the collective strong coupling regime of organic molecules with cavities or delocalized plasmonic modes have shown important light-induced modifications of material properties. Experimental and theoretical endeavors go into the direction of charge and energy transport enhancement [2-6], Förster resonance energy transfer (FRET) enhancement [7][8][9][10][11], modified chemical reactivity [12][13][14][15][16], polariton dynamics [17,18] etc. Oftentimes however, experiments rely on theoretical models developed for standard cavity quantum electrodynamics (cavity QED) [19][20][21] with two-level systems where one distinguishes between i) the Purcell regime, characterized by modifications of the spontaneous emission rates and ii) the strong coupling regime leading to the occurrence of hybrid light-matter states referred to as polaritons. Recent theoretical efforts aim at covering this gap by solving a generalized light-electronic-vibrations problem modeled as a Holstein-Tavis-Cummings Hamiltonian. Investigations aim at providing an understanding of the vibrationally induced cavity polariton asymmetry [18,22], vibrationally dressed polaritons [23], dark vibronic polaritons [24,25], developing a cavity Born-Oppenheimer theory [26,27] or deriving relevant simplified models for large scale numerics in the mesoscopic limit [28].We provide here an alternative path based on solving the Holstein-Tavis-Cummings dynamics at the level of operators rather than states. The basic model considers a molecular box (see Fig. 1(b)) comprised of an internal electronic transition coupled to any number of vibrational modes. Radiative decay and vibrational relaxation are included as stochastic source terms in a set of coupled standard quantum Langevin equations [29][30][31][32][33] for vibrational and polaron operators (similarly applied in optomechanical systems [34][35][36]). The method is nu...