Abstract: We study Dirac plasmons and their damping in spatially separated N-layer graphene structures at finite doping and temperatures. The plasmon spectrum consists of one optical excitation with square-root dispersion and N – 1 acoustical excitations with linear dispersion, which are undamped at zero temperature and finite doping within a triangular energy region outside the electron-hole continuum. In the long-wavelength limit the energy and weight of the optical plasmon modes increase, respectively, as the square root and linearly with N in agreement with recent experimental findings. The energy and weight of the upper-lying acoustical branches also increase with N. This increase is strongest for the uppermost acoustical mode, and we find that its energy can exceed at some value of momentum the plasmon energy in an individual graphene sheet. Meanwhile, the energy of the low-lying acoustical branches decreases weakly with N as compared with the single acoustical mode in double-layer graphene structures. Our numerical calculations provide a detailed understanding of the overall behavior of the wave-vector dependence of the optical and acoustical multilayer plasmon modes and show how their dispersion and damping are modified as a function of temperature, interlayer spacing, and inlayer carrier density in (un)balanced graphene multilayer structures. DOI: 10.1103/PhysRevB.87.085401