refractive index dielectric (e.g., silicon (Si), germanium and gallium phosphide, etc.) nanostructures have emerged. [1-4] The low loss feature [5-8] together with the unique optical responses associated with the electric and magnetic multipole resonances provides an additional degree of freedom for light manipulation. [1-3,9-11] Recently, in addition to the electric and magnetic multipole moments, toroidal multipole moments, which are derived by multipolar decomposition of induced currents in dielectric resonators, have been attracting significant attention. [12] The lowest order toroidal multipole, i.e., the toroidal dipole (TD), is a current loop along a meridian of a torus in a dielectric resonator. The TD mode has the far-field radiation pattern identical to that of the electric dipole (ED) mode, and thus the destructive interference of the two modes produces a non-radiating state. [12,13] The state is called "anapole." The radiation-less anapole state can be realized in a simple dielectric nanodisk (ND) structure when the diameter-to-height ratio is higher than ≈5. [12] At the anapole state, electromagnetic fields are confined inside a ND, [14] which are utilized for the enhancement of the nonlinear optical responses, [15-19] the emission rate by the Purcell effect, [20] the Raman scattering, [21] and the photocatalytic activity. [22] TD moments exist not only in an isolated dielectric nanostructure, but also in coupled dielectric nanostructures like oligomers [23,24] and two dimensional arrays. [25-28] Basharin et al. suggested the existence of TD resonances in cluster arrays of subwavelength high-index dielectric cylinders in the terahertz frequency region. [25] Liu et al. theoretically demonstrated a high Q-factor response and resultant strong field confinement by TD resonances in square arrays of silicon split nanodisks. [26] Luo et al. proposed a metasurface consisting of two silicon split-ring resonators and theoretically demonstrated a high-Q TD resonance in the near-IR region. [27] These researches indicate that compared to TD resonances of isolated nanostructure, those of coupled nanostructures have higher degree of freedom in designing TD mode-related optical responses. However, experimental studies on TD responses of coupled dielectric nanostructures are still limited [29,30] and the application for the control of light-matter interactions has never been reported. In this work, we propose a two-dimensional hexagonal array of very thin (20-50 nm in thicknesses) circular Si NDs as the metasurface supporting coupled TD modes. We first calculate