We explore the collective density oscillations of a collection of charged massive Dirac particles, in one, two and three dimensions and their one dimensional superlattice. We calculate the long wavelength limit of the dynamical polarization function
analytically, and use the random phase approximation to obtain the plasmon dispersion. The density dependence of the long wavelength plasmon frequency in massive Dirac systems is found to be different compared to systems with parabolic, and gapless Dirac dispersion. We also calculate the long wavelength plasmon dispersion of a 1d metamaterial made from 1d and 2d massive Dirac plasma. Our analytical results will be useful for exploring the use of massive Dirac materials as electrostatically tunable plasmonic metamaterials and can be experimentally verified by infrared spectroscopy as in the case of graphene [L. Ju. et. al., Nat. Nanotechnol. 6, 630 (2011)].
