  The nonlocal response of plasmonic materials and nanostructures is usually
described within a hydrodynamic approach which is based on the Euler-Drude
equation. In this work, we reconsider this approach within linear response
theory and employ Halevi's extension to this standard hydrodynamic model. After
discussing the impact of this improved model, which we term the Halevi model,
on the propagation of longitudinal volume modes, we accordingly extend the
Mie-Ruppin theory. Specifically, we derive the dispersion relation of
cylindrical surface plasmons. This reveals a nonlocal, collisional damping term
which is related to earlier phenomenological considerations of
limited-mean-free-path effects and influences both, peak width and amplitude of
corresponding resonances in the extinction spectrum. In addition, we transfer
the Halevi model into the time-domain thereby revealing a novel, diffusive
contribution to the current which shares certain similarities with
Cattaneo-type currents and analyze the resulting hybrid, diffusive-wave-like
motion. Further, we discuss the relation of the Halevi model to other
approaches commonly used in the literature. Finally, we demonstrate how to
implement the Halevi model into the Discontinuous-Galerkin Time-Domain
finite-element Maxwell solver and are able to identify an oscillatory
contribution to the diffusive current. The Halevi model thus captures a number
of relevant features beyond the standard hydrodynamic model. Contrary to other
extensions of the standard hydrodynamic model, its use in time-domain Maxwell
solvers is straightforward -- especially due its affinity to a class of
descriptions that allow for a clear distinction between bulk and surface
response. This is of particular importance for applications in nano-plasmonics
where nano-gap structures and other nano-scale features have to be modeled
efficiently and accurately.
