Plasmonic systems, such as metal nanoparticles, are widely used in different areas of application, going from biology to photovoltaics. The modeling of the optical response of such systems is of fundamental importance to analyze their behavior and to design new systems with required properties. When the characteristic sizes/distances reach a few nanometers, nonlocal and spill-out effects become relevant and conventional classical electrodynamics models are no more appropriate. Methods based on the Time-Dependent Density Functional Theory (TD-DFT) represent the current reference for the description of quantum effects. However, TD-DFT is based on knowledge of all occupied orbitals, whose calculation is computationally prohibitive to model large plasmonic systems of interest for applications. On the other hand, methods based on the orbital-free (OF) formulation of TD-DFT can scale linearly with the …