“…We consider regular arrays from Si NPs with R = 65 nm, arranged in a 2D rectangular lattice with h y = 480 nm and h x = 580 nm. A direct comparison with full-field simulations ([31] Figure 1) ( [38] Figure 3) has confirmed a reliability of the coupled dipole approximation (1) for arrays with similar pitches and the same R. Under a normal incidence with E 0 = (E 0x , 0, 0) and H 0 = 0, H 0y , 0 , arrays with these geometrical parameters exhibit ED and MD CLRs ( [37] Figure 2b) at λ ≈ 490 nm and λ ≈ 586 nm, respectively, which is the result of a hybridization between ED (λ ≈ 450 nm) and MD (λ Figure 2), it is insightful to consider an incidence with only one θ x,y varied keeping the other θ y,x = 0. Following this approach, it is possible to study separately ED and MD CLRs, while, for any other oblique incidence with θ x = 0 and θ y = 0, one can expect the optical response to be a superposition of the studied examples.…”