“…The extinction spectra of the pristine and the relaxed matrix are calculated using the DDA method, 39,40 as in our previous work. 41 Recently, this method has been proved suitable for both dielectric and plasmonic systems, 34,35,[42][43][44] and extended into multipole analysis for the large particle array. [45][46][47] Before introducing the calculation processes of array extinction spectra, here is a brief introduction of the DDA method used here, and its accuracy is evaluated.…”