LAUSR

The interaction between electromagnetic field and plasmonic nanostructures leads to both strong linear and nonlinear behaviors. In this paper, a time-domain hydrodynamic model for describing the motion of electrons in plasmonic nanostructures is presented, in which both surface and bulk nonlinearity are considered. A coupled Maxwell–hydrodynamic system capturing full-wave physics and free-electron dynamics is numerically solved with the parallel finite-difference time-domain method. The validation of the proposed method is presented by simulating a plasmonic metasurface. The linear response is compared with the Drude dispersion model and the nonlinear terahertz emission from a difference-frequency generation process is theoretically analyzed. The work is fundamentally important to design nonlinear plasmonic nanodevices, especially for efficient and broadband THz emitters.