  The unique structure of two-dimensional (2D) Dirac crystals, with electronic
bands linear in the proximity of the Brillouin-zone boundary and the Fermi
energy, creates anomalous situations where small Fermi-energy perturbations are
known to critically affect the electron-related lattice properties of the
system. The Fermi-surface nesting (FSN) conditions determining such effects via
electron-phonon interaction, require accurate estimates of the crystal's
response function $(\chi)$ as a function of the phonon wavevector q for any
values of temperature. Numerous analytical estimates of $\chi(q)$ for 2D Dirac
crystals beyond the Thomas-Fermi approximation have been so far carried out
only in terms of dielectric response function $\chi(q,\omega)$, for photon and
optical-phonon perturbations, due to relative ease of incorporating a
q-independent oscillation frequency in their calculation. However, models
accounting for Dirac-electron interaction with ever-existing acoustic phonons,
for which $\omega$ does depend on q and is therefore dispersive, are essential
to understand many critical crystal properties. The lack of such models has
often led to assume that the dielectric response function $\chi(q)$ in these
systems can be understood from free-electron behavior. Here, we show that,
different from free-electron systems, $\chi(q)$ calculated from acoustic
phonons in 2D Dirac crystals using the Lindhard model, exhibits a cuspidal
point at the FSN condition. Strong variability of $\frac{\partial\chi}{\partial
q}$ persists also at finite temperatures, while $\chi(q)$ may tend to infinity
in the dynamic case even where the speed of sound is small, albeit
nonnegligible, over the Dirac-electron Fermi velocity. The implications of our
findings for electron-acoustic phonon interaction and transport properties such
as the phonon line width derived from the phonon self energy will also be
discussed.
