Classical electronic and molecular dynamics simulation for optical response of metal system

An extended molecular dynamics simulation that incorporates classical free electron dynamics in the framework of the force-field model has been developed to enable us to describe the optical response of metal materials under the visible light electric field. In the simulation, dynamical atomic point charges follow equations of motion of classical free electrons that include Coulomb interactions with the oscillating field and surrounding atomic sites and collision effects from nearby electrons and ions. This scheme allows us to simulate an interacting system of metals with molecules using an ordinary polarizable force-field and preserves energy conservation in the case without applying an external electric field. As the first applications, we show that the presented simulation accurately reproduces (i) the classical image potential in a metal-charge interaction system and (ii) the dielectric function of bulk metal. We also demonstrate (iii) calculations of absorption spectra of metal nano-particles with and without a water solvent at room temperature, showing reasonable red-shift by the solvent effect, and (iv) plasmon resonant excitation of the metal nano-particle in solution under the visible light pulse and succeeding energy relaxation of the absorbed light energy from electrons to atoms on the metal and to the water solvent. Our attempt thus opens the possibility to expand the force-field based molecular dynamics simulation to an alternative tool for optical-related fields.