Overcoming Computational Bottlenecks in Quantum Hydrodynamics: A Volume-Based Integral Formalism

TL;DR

This paper introduces a volume-based integral formalism for quantum hydrodynamics that significantly improves computational efficiency in modeling the optical response of metallic nanoparticles, enabling advanced material responses to be simulated more easily.

Contribution

It develops a volume integral equation method for the Self-Consistent Hydrodynamic Drude Model, enhancing computational efficiency and adaptability for complex quantum plasmonic systems.

Findings

VIE method achieves similar performance to DE methods with less computational cost.

Mesoscopic material-response functions can be extracted without microscopic calculations.

The approach is adaptable to various material models and geometries.

Abstract

Mesoscopic models of the optical response of metals have emerged as fundamental building blocks in quantum plasmonics, in principle overcoming the computational bottlenecks of ab initio techniques by implementing aspects of the atomistic description of the metal in otherwise classical calculations. Nonetheless, even these approaches are eventually hindered by demanding computations due to sophisticated material response. Here, this issue is addressed for the advanced Self-Consistent Hydrodynamic Drude Model (SC-HDM), which captures both nonlocal electron dynamics and electron spill-out, through a Volume Integral Equation (VIE) method. Adopting an IE-based method shifts perspective from the commonly employed Differential Equation (DE)-based ones, demonstrating significant computational efficiency. The VIE approach is a valuable methodological scaffold: It addresses SC-HDM and simpler models, but can also be adapted to more advanced ones. For spherical nanoparticles (NPs), using the inherent symmetries, similar performance for three increasingly complicated material models is achieved, breaking the taboo that increased sophistication in material response requires taxing simulations. Mesoscopic material-response functions can be readily extracted from the VIE implementation, thus circumventing the need for lengthy microscopic calculations. This method opens a new way of modeling quantum hydrodynamic NPs and will serve as essential benchmarking tool for recipes addressing more complicated geometries.