The existence of an acoustic plasmon in extrinsic (doped or gated) monolayer graphene was found recently in an ab initio calculation with the frozen lattice [M. Pisarra et al., arXiv:1306.6273, 2013. By the fully dynamic density-functional perturbation theory approach, we demonstrate a strong coupling of the acoustic plasmonic mode to lattice vibrations. Thereby, the acoustic plasmon in graphene does not exist as an isolated excitation, but it is rather bound into a combined plasmonphonon mode. We show that the coupling provides a mechanism for the bidirectional energy exchange between the electronic and the ionic subsystems with fundamentally, as well as practically, important implications for the lattice cooling and heating by electrons in graphene.Known for its extraordinary properties and vast potential applications [1], graphene -a two-dimensional crystal comprised of a honeycomb lattice of carbon atomscontinues to receive much attention as it reveals new remarkable features [2][3][4][5][6][7]. For one of the recent findings, an acoustic plasmon (APl) (plasmon with linear wave-vector dispersion) has been predicted theoretically in an extrinsic free-standing monolayer graphene [7]. This finding is extraordinary considering that APl generation conventionally involves a surface state immersed in the bulk of a metal [8].Exhibiting linear wave-vector dispersion, acoustic APl persists down to low frequencies, where it can be expected to interact with phonon oscillations. The possibility of coupling these two types of elementary excitations motivates questions of fundamental physics as well as of potential applications. In this Letter we show that the APl -phonon coupling indeed occurs in the electron-doped graphene and it provides a mechanism for the bidirectional energy exchange between the electronic and ionic subsystems. The conventional treatment of lattice vibrations by frequency-independent densityfunctional perturbation theory (DFPT) [9] is inadequate for capturing the essentially dynamic nature of the coupled plasmon-phonon modes, and we therefore implement a fully dynamic approach treating the electron-hole, plasmon, and phonon elementary excitations on the equal footing [10].Our ab initio calculations for monolayer graphene employ the full-potential linear augmented plane-wave (FP-LAPW) code Elk [11]. The super-cell geometry is utilized with a separation of the layers in the z direction of 40 bohr, which effectively ensures the non-interaction between the layers. The local-density approximation to the exchange-correlation potential [12, 13] is used. * nazarov@gate.sinica.edu Acoustic-plasmon and phonons in graphene.-We start by reproducing the APl and phonon spectra of graphene without the coupling of the two excitations. In Fig. 1, left panel, the energy-loss function of graphene is plotted for a number of equidistant values of the wave-vector. The calculation with the carbon atoms fixed at their equilibrium positions has been used. The APl can be easily recognized by the linear dispersion of the peak ...