A new time-dependent density-functional method for molecular plasmonics: Formalism, implementation, and the Au144(SH)60 case study

We describe the implementation and application of a recently developed time-dependent density-functional theory (TDDFT) algorithm based on the complex dynamical polarizability to calculate the photoabsorption spectrum of large metal clusters, with specific attention to the field of molecular plasmonics. The linear response TDDFT equations are solved in the space of the density fitting functions, so the problem is recast as an inhomogeneous system of linear equations whose resolution needs a numerical effort comparable to that of a SCF procedure. The construction of the matrix representation of the dielectric susceptibility is very efficient and is based on the discretization of the excitation energy, so such matrix is easily obtained at each photon energy value as a linear combination of constant matrix and energy-dependent coefficients. The code is interfaced to the Amsterdam Density Functional (ADF) program and is fully parallelized with standard message passing interface. Finally, an illustrative application of the method to the photoabsorption of the Au144(SH)60 cluster is presented.