Abstract
A novel fully atomistic multiscale classical approach to model the optical response of solvated real-size plasmonic nanoparticles (NPs) is presented. The model is based on the coupling of the Frequency Dependent Fluctuating Charges and Fluctuating Dipoles (ωFQFμ), specifically designed to describe plasmonic substrates, and the polarizable Fluctuating Charges (FQ) classical force field to model the solvating environment. The resulting ωFQFμ/FQ approach accounts for the interactions between the radiation and the NP, as well as with the surrounding solvent molecules, by incorporating mutual interactions between the plasmonic substrate and solvent. ωFQFμ/FQ is validated against reference TD-DFTB/FQ calculations, demonstrating remarkable accuracy, particularly in reproducing plasmon resonance frequency shifts for structures below the quantum-size limit. The flexibility and reliability of the approach are also demonstrated by simulating the optical response of homogeneous and bimetallic NPs dissolved in pure solvents and solvent mixtures.