Plasmonic evolution maps for planar metamaterials

1. INTRODUCTION

Over the past decade, planar metamaterials that act as two-dimensional single-layer counterparts to fully three-dimensional metamaterials have drawn intense attention [1–4]. As compared to three-dimensional metamaterials, planar metamaterials are considerably simpler in terms of fabrication by top-down techniques. In addition to a negative refractive index, which is the most remarkable property of metamaterials [5], many other attractive phenomena for condensed-state and quantum-classical analogies have been reported in planar metamaterials. These phenomena include the Fano resonance [6,7], the electromagnetically induced transparency [8–10], near-zero refractive index [11,12], and the topological edge states [13–15]. From the microwave range to the terahertz, infrared, and optical ranges, a variety of planar metamaterials have been designed and applied in perfect absorber [16], sensing [17,18], circular dichroism [19–21], wave plates [22–24], coherent modulation [25,26], electrical modulators [27–34], metalenses [35–38], and optical rulers [39].

Accompanying this fast development, the underlying theoretical methods always play a fundamental role in designing planar metamaterials with particular phenomena and applications. Considering that the unit cell of a planar metamaterial is usually made of a plasmonic atom or molecule, it becomes critical to understand the mode’s origin and features during the analysis. Some classical electromagnetic analysis theories and methods, such as the Mie theory [40], the coupled-mode theory [9,41], the multipole-expansion method [42,43], and the inductive-capacitive (LC) equivalent-circuit model [44], have been widely applied to predict the resonant modes of simple planar metamaterials. For complex structures, the electromagnetic simulation tools, such as the finite-element method [18–21,27–33], the finite-difference time-domain method [9,11,16,17,35–38], and the moment method [45], have become more universal for verifying a mode’s wavelength and features. In addition, Rechberger et al. and Prodan et al. have established particular plasmonic hybridization models (PHMs) to well describe the mode coupling and energy states in complex plasmonic molecules [46–49].

In many studies, mode analysis directly based on the final field/charge/current distribution is simple and efficient for us to identify a mode’s type and order in simple plasmonic atoms, such as nanodisks [11], nanorods [35], and L-shaped nanostructures [24]. However, for more complex plasmonic atoms, such as square rings [50], split-ring resonators [27–33], H-shaped nanostructures [51], and S-shaped nanostructures [23], it is difficult to deeply understand a mode’s origin because the final field/charge/current distribution does not reveal the complicated formation process of a plasmonic mode. Based on the PHM, we can efficiently and clearly understand the coupling mechanism of low-order modes in plasmonic molecules [47–49,52–55]. However, such a model will become less efficient or even invalid for high-order modes in a plasmonic molecule constructed by complex plasmonic atoms.

Considering the above issues, in this paper, we would like to propose a new analysis method based on natural evolution theory in order to understand the modes’ origin and coupling mechanism in planar metamaterials constructed by complex plasmonic atoms and molecules. We assume that a complex plasmonic atom/molecule initially grows from one basic element, and introduce several structural and mode evolution operators to construct evolutional maps that can well describe the modes’ formation process. We will also clearly reveal the dependency of a structure’s symmetry on the formation of a plasmonic mode based on our proposed mode evolution method (MEM).

2. RESULTS

A. Planar Metamaterial Samples

We consider two typical complex planar metamaterials in our study, i.e., the complementary split-ring resonator (CSRR) and the H-shaped air slit. As shown in Fig. 1(a), the CSRR contains a square air ring with one or more splits on the upper/bottom sides of the air ring. The side length and the slit width of the square air ring are defined as $l$ and $w$, respectively. The width of split is denoted as $g$, and $\mathrm{\Delta }x$ represents the horizontal displacement of the split from the center. In the experiment, six CSRR planar metamaterial samples with different values for the split number and symmetry were fabricated. A 50-nm thick gold layer was first deposited on the glass substrate by thermal evaporation (Oerlikon Univex 250). The patterns of CSRR planar metamaterials were then transferred onto the gold layer by using focused ion beam milling (FEI Helios NanoLab 650). It can be seen from the scanning electron microscope (SEM) images that all of the samples have high quality, with good homogeneity and consistency [Fig. 1(b)]. The structural parameters of these samples are approximately $l=550\mathrm{nm}$, $w=50\mathrm{nm}$, and $g=80\mathrm{nm}$. The period of each array is 625 nm in both the $x$ and $y$ directions. Sample I without any split was fabricated for comparison purposes. Sample II and Sample III have one split at the upper center ($\mathrm{\Delta }x=0$) and the upper right side ($\mathrm{\Delta }x=150\mathrm{nm}$), respectively. Sample IV, Sample V, and Sample VI have two splits at the upper and bottom sides of the air ring with different values of $\mathrm{\Delta }x$. As far as the mirror and center inversion symmetries are concerned, Sample I and Sample IV satisfy both the $x/y$ mirror and center inversion symmetries, while Sample II and Sample VI satisfy only the $y$ mirror symmetry and the center inversion symmetry, respectively. Sample III and Sample V show neither mirror nor center inversion symmetry.

 

Fig. 1. Schematic diagram and SEM images of the (a), (b) CSRR and (c), (d) H-shaped complementary planar metamaterials. From Sample I to Sample VI, different numbers of splits are introduced and arranged in designed symmetry/asymmetry. $\mathrm{\Delta }x$ represents the horizontal displacement of the split from the center. The symmetry/asymmetry of the H-shaped air slit is modulated by introducing a vertical displacement $\mathrm{\Delta }y$ of the horizontal air slit shift along the $y$ direction. The scale bars in the figure are all 500 nm.

Download Full Size | PPT Slide | PDF

For the H-shaped air slit, the length and the width of the vertical air slits are defined as $l$ and $w$, and the length and width of the horizontal air slit are defined as $g$ and $h$, respectively [Fig. 1(c)]. A vertical displacement $\mathrm{\Delta }y$ for the horizontal air slit is introduced along the vertical direction. As indicated by Fig. 1(d), we fabricated three H-shaped air slit planar metamaterial samples with different values of $\mathrm{\Delta }y$, i.e., Sample ${\mathrm{I}}^{\prime }$ ($\mathrm{\Delta }y=0$), Sample ${\mathrm{II}}^{\prime }$ ($\mathrm{\Delta }y=20\mathrm{nm}$), and Sample ${\mathrm{III}}^{\prime }$ ($\mathrm{\Delta }y=40\mathrm{nm}$). The period of each array is 305 nm in both the $x$ and $y$ directions, and the other structural parameters are approximately $l=200\mathrm{nm}$, $w=50\mathrm{nm}$, $g=100\mathrm{nm}$, and $h=40\mathrm{nm}$. Sample ${\mathrm{I}}^{\prime }$ satisfies the $x/y$ mirror symmetry and center inversion symmetry, while Sample ${\mathrm{II}}^{\prime }$ and Sample ${\mathrm{III}}^{\prime }$ have only $y$ mirror symmetry.

B. Reflection Spectra and Conventional Mode Analysis

The reflection spectra of the samples were measured with a visible-near-infrared microspectrophotometer (Jasco MSV-5200). A commercial software package, Lumerical FDTD Solutions, was used to simulate the reflection spectra of samples. In our simulation, the plane wave source was used to excite the eigenmodes of planar metamaterials. The periodic boundary condition was used along the $x$ and $y$ directions, and the perfectly matched layer (PML) absorbing boundary condition was used along both the $z$ directions. The grid size was set as 2 nm for the $x$, $y$, and $z$ directions. The complex permittivity of the gold was fitted to the Johnson and Christy data.

Figure 2 shows the experimental and simulated reflection spectra of the CSRR samples when the light is normally incident to the metamaterials with an electric field polarized along the $y$ direction. The experimental spectra match well with the simulation spectra in terms of the resonance number and location. It is found that the resonance number has a close relationship with the symmetry of the CSRRs. Specifically speaking, Sample I, which has the highest symmetry, shows two main resonances in the near-infrared band [Fig. 2(a)]. One resonance is located at 1057 nm (I-1). Another resonance located at around 3200 nm is not shown here due to the spectrum limit in the experiment. When one split is introduced without displacement ($\mathrm{\Delta }x=0$), the horizontal symmetry of the square air ring is broken, and we can observe three resonances (II-1 at 999 nm, II-2 at 1295 nm, and II-3 at 2165 nm) in the spectra of Sample II [Fig. 2(b)]. For Sample III, with a positive displacement of 150 nm, both the horizontal and vertical symmetries of the square air ring are broken, and there is one more minor resonance (III-4 at 1221 nm) in the spectra [Fig. 2(c)]. In comparison, when two splits without displacement are introduced in Sample IV, both the horizontal and vertical symmetries are recovered, and the resonance number is reduced to two again [Fig. 2(d)]. Similar to Sample III, Sample V also shows four resonances due to the completely broken horizontal and vertical symmetries [Fig. 2(e)]. The last sample, Sample VI, only has a center inversion symmetry, but its reflection spectra are similar to those of Sample IV, which also satisfies the center inversion symmetry [Fig. 2(f)]. Actually, more symmetry breaking of CSRRs usually corresponds to more resonances because high-order modes prefer to stay in the low-symmetry structures [56–58]. It should be noted that when the plane wave is polarized along the $x$ direction, the resonances are also sensitive to the symmetries of the six CSRR samples. The relation between resonance and structural symmetry is similar to the $y$-polarization case.

 

Fig. 2. Experimental and simulated reflection spectra of (a) Sample I, (b) Sample II, (c) Sample III, (d) Sample IV, (e) Sample V, and (f) Sample VI. The solid curve in each experimental spectrum serves as a guide for easier comparison.

Download Full Size | PPT Slide | PDF

One conventional method for the mode analysis of SRRs/CSRRs is simulating the final electric field and charge distributions of resonant modes. From Figs. 3(a), 3(d), 3(g), and 3(j), it is clear to see that the electric-field amplitude distribution of the representative modes in Fig. 2 retains the structural symmetry of the samples well. The higher-order mode contains more subregions with concentrated electric fields. Based on the corresponding charge distribution, as presented in Figs. 3(b), 3(e), 3(h), and 3(k), II-5 and III-5 show a clear dipole-like charge distribution, and the higher-order modes show more pairs of dipole-like charges. However, it is difficult to understand the modes’ origin of these higher-order modes based on the final electric field and charge distributions.

 

Fig. 3. Conventional mode analysis for Samples I–VI. (a), (d), (g), (j) Electric field distributions; (b), (e), (h), (k) charge distributions; (c), (f), (i), (l) current distributions. The big arrows serve as a guide to show the current directions at the metal/air surface. The symmetrically distributed currents in the plasmonic atom samples (I–III) are highlighted by green color, and the currents in the right atom of the plasmonic molecule samples (IV–VI) are highlighted by blue color.

Download Full Size | PPT Slide | PDF

Another widely used method is using an LC equivalent-circuit model to predict the resonant modes of SRRs/CSRRs [44]. In particular, the resonant frequency of the fundamental LC mode is written as $f=(2\pi \sqrt{L_{\mathrm{eff}}{C}_{\mathrm{eff}}}{)}^{-1}$, where $L_{\mathrm{eff}}$ and $C_{\mathrm{eff}}$ are the effective inductance and capacitance, respectively. According to the Babinet principle [59], CSRRs work as dual SRRs. In our study, light is normally incident to the samples, and the electric moments associated with the induced currents are perpendicular to the plane of the CSRRs. Therefore, the interaction of the induced electric moments with the electric field of the incident wave is negligible, while the circulating surface current distribution is generally used to identify the possible LC resonance in the CSRRs [60]. As shown in Figs. 3(c), 3(f), 3(i), and 3(l), II-5 and III-5 are proposed to originate from the fundamental LC resonance because the light can drive a couple of circulating surface currents in the split of each CSRR unit cell. The current flows become more complex for the higher-order modes, and it is difficult to identify these modes based on the LC model.

Figures 4(a) and 4(b) show the experimental and simulated reflection spectra for three H-shaped air slit planar metamaterial samples when the plane wave is polarized along the $y$ direction. Sample ${\mathrm{I}}^{\prime }$ is a symmetrical structure with only one resonance ${\mathrm{I}}^{\prime }-1$ located at 1446 nm. Sample ${\mathrm{II}}^{\prime }$ shows a new resonance ${\mathrm{II}}^{\prime }-2$ at 889 nm when its vertical symmetry is broken. As vertical symmetry is further broken in Sample ${\mathrm{III}}^{\prime }$, the additional resonance ${\mathrm{III}}^{\prime }-2$ at 837 nm is enhanced with a slight blueshift. The corresponding electric-field and charge distributions of resonant modes are given in Figs. 4(c) and 4(d). It is clear to see that the electric field of mode ${\mathrm{I}}^{\prime }-1$ is concentrated in the horizontal air slit and mode ${\mathrm{I}}^{\prime }-1$ contains three pairs of dipole-like charges. For ${\mathrm{II}}^{\prime }-1$, in addition to the concentrated electric field in the horizontal air slit, we can also observe “hot spots” at the lower corners that connect the horizontal and vertical air slits. For ${\mathrm{II}}^{\prime }-2$, only hot spots are found near the upper corners. Both ${\mathrm{II}}^{\prime }-1$ and ${\mathrm{II}}^{\prime }-2$ contain three pairs of dipole-like charges with different charge distributions. More distinct hot spots can be observed from ${\mathrm{III}}^{\prime }-1$ and ${\mathrm{III}}^{\prime }-2$ when the vertical symmetry is further broken. It is difficult to understand the origin of these modes based on the final electric field and charge distributions. From the surface current distributions in Fig. 4(e), the modes $\mathrm{I}-1$, ${\mathrm{I}}^{\prime }-1$, and ${\mathrm{II}}^{\prime }-1$ are proposed to originate from the fundamental LC resonance because the light can drive a couple of circulating surface currents in the gap of two vertical slits. This LC resonance is less influenced by the vertical displacement of the horizontal air slit because the total length of the circulating surface currents is not changed. The higher-order modes, however, are also difficult to be identified based on the LC model. When the plane wave is polarized along the $x$ direction, the resonances are almost not influenced by the vertical displacement of the horizontal air slit.

 

Fig. 4. Conventional mode analysis for three H-shaped complementary planar metamaterials. (a), (b) Experimental and simulated reflection spectra; (c) electric field distributions; (d) charge distributions; (e) current distributions. The symmetrically distributed currents are highlighted by green color.

Download Full Size | PPT Slide | PDF

C. Mode Analysis Based on PEMs

The above mode analysis based on conventional methods shows obvious limitations for deeply understanding the modes’ origin in complex planar metamaterials. Here we propose a new method called the MEM with a core analysis tool, i.e., plasmonic evolution maps (PEMs) in order to give a qualitative explanation of the modes’ origin. The PEMs contain three elements, i.e., structural operator, mode operator, and characteristic spectra/charge distribution. In the following contents, we will illustrate how the PEMs work for different plasmonic atom samples (Samples I–III and Samples ${\mathrm{I}}^{\prime }–{\mathrm{III}}^{\prime }$) and plasmonic molecule samples (Samples IV–VI).

Figures 5(a)–5(d) present the PEMs for Sample I, which gradually grows from a horizontal air slit (primary state) into the square air ring (final state). The structural operators during the evolutional process include growth and crossover, while the related structural operators include expansion, preservation, and mutation. From Fig. 5(a), we can see that the primary horizontal air slit shows only one dipole-like resonance at 1612 nm. When the L-shaped air slit is constructed via the structural operator growth, besides the dipole-like resonance at 3205 nm, it can support two additional higher-order plasmonic modes due to the increased air slit’s length, i.e., a quadrupole-like charge at 1716 nm and a hexapole-like charge at 1107 nm [Fig. 5(b)]. It should be noted that these higher-order modes are not real electric quadrupoles or hexapoles. They are actually the higher-order surface plasmon polariton (SPP) modes that satisfy the standing-wave formula $\lambda =2n_{\mathrm{eff}}L/N$, where $L$ is the total length of the air slit, $n_{\mathrm{eff}}$ is the effective refractive index of the SPP waveguide, and $N$ is the order of eigenmode. For the L-shaped air slit with a length of 1100 nm, the calculated eigenmodes are located at 3410 nm ($N=1$, $n_{\mathrm{eff}}=1.55$), 1705 nm ($N=2$, $n_{\mathrm{eff}}=1.55$), and 1173 nm ($N=3$, $n_{\mathrm{eff}}=1.6$), which are close to the simulation results. As a result, a mode operator expansion is defined to describe the mode evolution during the growth of the structure. Under expansion operation, the quantity of the mode is continuously increasing and the mode’s order clearly depends on how many pairs of electric dipoles it has.

 

Fig. 5. Mode analysis based on PEMs for (a)–(d) Sample I and (e)–(h) Sample II. D, dipole-like; Q, quadrupole-like; H, hexapole-like; O, octopole-like charge distributions; SCD/SCQ/SCH/SCO, symmetrically coupled dipole-like/quadrupole-like/hexapole-like/octopole-like charge distributions; ASSCD/ASSCQ/ASSCH, antisymmetrically coupled SCD/SCQ/SCH modes.

Download Full Size | PPT Slide | PDF

When the L-shaped air slit further connects with a vertical air slit on the right side, the formed U-shaped air slit shows 10 pole-like, hexapole-like, and dipole-like charge distributions at 1052, 1794, and 4688 nm, respectively [Fig. 5(c)]. This L-U formation does not belong to a growth process because the number of dipolar pairs of three modes is not continuous. In fact, the 10 pole-like and hexapole-like charge distributions are symmetrically coupled hexapole-like (SCH) and symmetrically coupled quadrupole-like (SCQ) charge distributions considering that the U-shaped air slit possesses $y$ mirror symmetry, and the driving electric field shows symmetric polarization along the $x$ direction. The real mode order for SCH and SCQ is $N=3$ and $N=5$, respectively. In other words, the SCH/SCQ in the U-shaped air slit (child structure) can be regarded as the overlapping (crossover) result of the two symmetric hexapole-like/quadrupole-like charge distributions in the left and right L-shaped air slits (parent structures). As a result, a mode operator preservation is defined to describe the mode evolution during the crossover of the structure. The preservation operation means that the child structure retains the same or similar genes (charge distributions) of the parent structures and that its charge distribution is the internal overlapping result of the parent structures after crossover. The preservation operation may narrow the linewidth of the eigenmode (such as SCQ).

The final square air ring can be regarded as the second-generation child structure of two parent U-shaped air slits after crossover [Fig. 5(d)]. Considering that the square air ring possesses $x$ mirror symmetry and the driving electric field shows antisymmetric polarization along the $y$ direction, resonances 1 and 2 are actually the antisymmetrically coupled SCH (ASSCH) and antisymmetrically coupled SCD (ASSCD) modes, respectively. Under the preservation operation, the ASSCH mode shows an enhanced line shape compared to the SCH mode due to the constructive coupling along the vertical air slits, while the ASSCD mode (two parallel short dipoles) shows a great blueshift with respect to the SCD mode (a long dipole) due to the destructive coupling along the vertical air slits. In particular, there is no antisymmetrically coupled SCQ (ASSCQ) in the spectrum due to the destructive coupling between two antisymmetric SCQ modes along the vertical air slits. This missing mode phenomenon behaves like a mutation in the evolution process. Therefore, a mode operator mutation is defined to describe the unpreserved charge distributions during the crossover.

Figures 5(e)–5(h) present the PEMs for Sample II. As compared to the square air ring, this CSRR with one split shows two main different points during the mode evolution. One difference is that the L-U formation still belongs to the growth process because the mode’s quantity is further increased and the mode’s order is still continuous from Figs. 5(f) and 5(g). Another difference is that the mutation operation is absent during the crossover process of two asymmetrical U-shaped air slits arranged in a $y$ mirror symmetry [Fig. 5(h)]. Due to completely broken symmetry, the mode evolution of Sample III always belongs to the growth process from the primary state to the final state.

Figure 6 presents the PEMs for Sample ${\mathrm{II}}^{\prime }$, which can be divided into two back-to-back U-shaped air slits. For the upper U-shaped air slit, a dipole-like charge evolutional process (D-SCD) can be observed [Figs. 6(a)–6(c)], while the lower U-shaped air slit corresponds to both dipole (${\mathrm{D}}^{\prime }-{\mathrm{SCD}}^{\prime }$) and quadrupole-like charge (${\mathrm{Q}}^{\prime }-{\mathrm{SCQ}}^{\prime }$) evolutional processes [Figs. 6(d)–6(f)]. The formations of SCD, ${\mathrm{SCD}}^{\prime }$, and ${\mathrm{SCQ}}^{\prime }$ modes in Sample ${\mathrm{II}}^{\prime }$ are similar to those in Figs. 5(a)–5(c). As shown in Fig. 6(g), the final resonances 1 and 2 of the H-shaped air slit are the crossover results of the upper and lower U-shaped air slits. Based on the charge distribution in Fig. 6(h), resonances 1 and 2 are the combination of $\mathrm{SCD}+{\mathrm{SCD}}^{\prime }$ and $\mathrm{SCD}+{\mathrm{SCQ}}^{\prime }$ via the preservation operation, respectively. The PEMs for Sample ${\mathrm{I}}^{\prime }$ and Sample ${\mathrm{III}}^{\prime }$ are similar to those of Sample ${\mathrm{II}}^{\prime }$.

 

Fig. 6. Mode analysis based on PEMs for Sample ${\mathrm{II}}^{\prime }$.

Download Full Size | PPT Slide | PDF

For Sample IV, Sample V, and Sample VI, which are plasmonic molecule samples, we define a structural operator (group) and two mode operators (preservation and learning) to describe the mode evolution in the plasmonic molecules. Preservation means that the charge distribution of each plasmonic atom keeps unchanged when they are combined in a group. Learning means that the charge distribution of a plasmonic atom in the group is affected by neighboring atoms. As shown in Fig. 7(a), for Sample V, resonance 1 is the preservation result of modes L1 and R1, while resonances 2 and 3 are the learning results from modes L2 and R2 by comparing with the charge distributions at the square and circle positions. As a comparison, the resonances of Sample VI, which contains the same two plasmonic atoms arranged in center inversion symmetry, are all due to the preservation operation during the group process [Fig. 7(b)]. The PEMs for Sample IV, which contains the same two plasmonic atoms arranged in mirror symmetry, are similar to Sample VI. This is because the atoms in Sample IV and Sample VI have the same order eigenmodes due to a fixed total length, and the mode preservation during the crossover and group processes does not change the number and resonant position of eigenmodes.

 

Fig. 7. Mode analysis based on PEMs for (a) Sample V and (b) Sample VI. The square and circle in (a) correspond to the positions of modes 2 and 3, respectively.

Download Full Size | PPT Slide | PDF

Based on these PEMs for plasmonic atom and molecule samples, it is found that the use of structural and mode operators strongly depends on the structural symmetry of samples. The physical mechanism behind the mode operators includes the near- or far-field excitation of plasmonic modes as well as the intrinsic or external coupling between plasmonic modes. The PEMs and the related mechanism construct the MEM for mode analysis in the plasmonic atoms and molecules, as summarized in Table 1.

Table 1. Summarization of MEM for Plasmonic Atoms and Molecules

View Table

3. DISCUSSION AND CONCLUSION

For a particular plasmonic atom or molecule, how to construct suitable PEMs is very important for the application of the MEM. General guidance for this issue is given as follows. First, the PEMs of a plasmonic molecule are fixed via the group of constituted plasmonic atoms. Second, the PEMs of a constituted plasmonic atom start from the growth process based on a primary state. If the intermediate state of the plasmonic atom shows a structural symmetry or the continuity of the mode’s order is broken, the growth process should be terminated and replaced by the crossover process. Third, the choice of primary and intermediate states is not fixed, and thus a constituted plasmonic atom may have different PEMs. For instance, the PEMs of Sample ${\mathrm{II}}^{\prime }$ can be constructed following the sequence of either L-U-H or L-T-H. Both of them can explain the final resonant modes well with different coupling mechanisms. Based on the above guidance, we have successfully verified the PEMs for other typical planar metamaterials, including the S-shaped planar metamaterial [23], the electric CSRR [29], and the hybrid CSRR [61]. In addition to the complementary planar metamaterials, the PEMs can also be applied for the normal type. For example, we have constructed the PEMs for the normal SRRs related to Sample I and Sample II and found that the SRRs show a one-to-one correspondence of plasmonic resonances to those of the CSRRs in Fig. 5. Therefore, the PEMs also satisfy the Babinet principle [59].

In principle, the MEM can be applied to arbitrary plasmonic atoms and molecules. For a plasmonic atom, the MEM works as a kind of ab initio theory applied in nanophotonics and plasmonics. For a plasmonic molecule, the MEM becomes a supplement to the widely used PHM, which contains one structural operator (“hybridization”) and two mode operators (“bonding” and “anti-bonding”), as well as a core analysis tool (energy-state maps). The PHM is mainly used among lower orders of plasmonic modes such as dipole–dipole and dipole–quadruple coupling [47–49,52–55], while there is no limitation to the mode’s order based on the MEM. A combination of the MEM and PHM is able to provide comprehensive mode analysis for complex metamaterials.

In addition to mode analysis, some other potential applications are expected based on the MEM. For instance, we can freely tailor the bright (superradiant) and dark (subradiant) modes according to the symmetry-dependent PEMs of different plasmonic atoms and molecules, which is helpful in further studies on the coherent Mie scattering, Fano resonance, and electromagnetically induced transparency. We can also design a genetic metasurface based on the MEM. A conventional metasurface usually deals with information based on the Berry phase change of plasmonic atoms [35–38]. In comparison, a genetic metasurface will deal with information by writing or reading the evolutional relationships of internal plasmonic atoms.