Second harmonic generation in metasurfaces with multipole resonant coupling

We report on the numerical demonstration of enhanced second harmonic generation (SHG) originating from collective resonances in plasmonic nanoparticle arrays. The nonlinear optical response of the metal nanoparticles is modeled by employing a hydrodynamic nonlinear Drude model implemented into Finite-Difference Time-Domain (FDTD) simulations, and effective polarizabilities of nanoparticle multipoles in the lattice are analytically calculated at the fundamental wavelength by using a coupled dipole-quadrupole approximation. Excitation of narrow collective resonances in nanoparticle arrays with electric quadrupole (EQ) and magnetic dipole (MD) resonant coupling leads to strong linear resonance enhancement. In this work, we analyze SHG in the vicinity of the lattice resonance corresponding to different nanoparticle multipoles and explore SHG efficiency by varying the lattice periods. Coupling of electric quadrupole and magnetic dipole in the nanoparticle lattice indicates symmetry breaking and the possibility of enhanced SHG under these conditions. By varying the structure parameters, we can change the strength of electric dipole (ED), EQ, and MD polarizabilities, which can be used to control the linewidth and magnitude of SHG emission in plasmonic lattices. Engineering of lattice resonances and associated magnetic dipole resonant excitations can be used for spectrally narrow nonlinear response as the SHG can be enhanced and controlled by higher multipole excitations and their lattice resonances. We show that both ED and EQ-MD lattice coupling contribute to SHG, but the presence of strong EQMD coupling is important for spectrally narrow SHG and, in our structure, excitation of narrow higher-order multipole lattice resonances results in five times enhancement.