Nonlocal Hydrodynamic Model with Viscosive Damping and Generalized Drude-Lorentz Term

The response of plasmonic metal particles to an electromagnetic wave produces significant features at the nanoscale level. Different properties of the internal composition of a metal, such as its ionic background and the free electron gas, begin to manifest more prominently. As the dimensions of the nanostructures decrease, the classical local theory gradually becomes inadequate. Therefore, Maxwell's equations need to be supplemented with a relationship determining the dynamics of current density which is the essence of nonlocal plasmonic models. In this field of physics, the standard (linearized) hydrodynamic model (HDM) has been widely adopted with great success, serving as the basis for a variety of simulation methods. However, ongoing efforts are also being made to expand and refine it. Recently, the GNOR (general nonlocal optical response) modification of the HDM has been used, with the intention of incorporating the influence of electron gas diffusion. Clearly, from the classical description of fluid dynamics, a close relationship between viscosive damping and diffusion arises. This offers a relevant motivation for introducing the GNOR modification in an alternative manner. The standard HDM and its existing GNOR modification also do not include the influence of interband electron transitions in the conduction band and other phenomena that are part of many refining modifications of the Drude-Lorentz and other models of metal permittivity. In this article, we present a modified version of GNOR-HDM that incorporates the viscosive damping of the electron gas and a generalized Drude-Lorentz term. In the selected simulations, we also introduce Landau damping, which corrects the magnitude of the standard damping constant of the electron gas based on the size of the nanoparticle. We have chosen a spherical particle as a suitable object for testing and comparing HD models and their modifications because it allows the calculation of precise analytical solutions for the interactions and, simultaneously, it is a relatively easily fabricated nanostructure in practice. Our contribution also includes our own analytical method for solving the HDM interaction of a plane wave with a spherical particle. This method forms the core of calculations of the characteristic quantities, such as the extinction cross-sections and the corresponding components of electric fields and current densities.