Metamaterials, and their 2D counterpart metasurfaces, made of a number of artificial micro/nanostructures with subwavelength dimensions, have attracted increasing interest due to the extraordinary capability for exterior light manipulation. [1][2][3][4] The subwavelength thickness of metasurfaces along the propagating direction is helpful for the miniaturization of optical systems and improving the efficiency of photonic devices. [5,6] As for light manipulation, one important application is the optical field localization, which is traditionally done with Fabry-Perot cavities or some echo wall modes. However, their dimensions are usually comparable to the wavelength. If the optical field can be further localized in a subwavelength scale, or the optical field can be compressed in both the spatial and time domains, it would be very beneficial for nonlinear optical devices. [7] The investigation on metasurfaces with resonant building blocks has inspired rapid development for subwavelength optics, in which metasurfaces are used for novel wavefront engineering and extreme optical localization. [8] Various plasmonic and Mie metasurfaces have been proposed for enhancing light-matter interactions and localizing light fields in a subwavelength volume, which is usually accompanied by the generation of high-Q-factor resonances [9,10] ; for example, the high-Q "trapped mode" [11] with a sharp spectral response has been exploited for applications, including ultrasensitive sensing, [12] high-Q filtering, [13] and optical switching. [14] Recently, metasurfaces/metamaterials with strong toroidal excitations have been demonstrated as a novel route for realizing extreme optical field localization and higher-Q optical response. [15][16][17][18][19][20][21] Different from electric and magnetic excitations in traditional multipole expansions, toroidal excitations arise from the decomposition of the momentum tensors with currents flowing on the surface of the torus along its meridians. [22][23][24][25] The toroidal excitations taking advantage of weak free-space coupling properties are beneficial for achieving ultra-high-Q response in a metamaterial, which has inspired numerous applications, [24] such as metaswitches, [26,27] photodetectors, [28] metamodulators, [29] metasensors, [30][31][32] and nonlinear harmonic signal generators. [33] Their unique field localization configuration is also helpful for