Principle and Applications of Multimode Strong Coupling Based on Surface Plasmons

In the past decade, strong coupling between light and matter has transitioned from a theoretical idea to an experimental reality. This represents a new field of quantum light–matter interaction, which makes the coupling strength comparable to the transition frequencies in the system. In addition, the achievement of multimode strong coupling has led to such applications as quantum information processing, lasers, and quantum sensors. This paper introduces the theoretical principle of multimode strong coupling based on surface plasmons and reviews the research related to the multimode interactions between light and matter. Perspectives on the future development of plasmonic multimode coupling are also discussed.


Introduction
Surface plasmon polaritons [1][2][3] (SPPs) are localized electromagnetic wave generated by the collective oscillation of electrons on metal surfaces, that are decayed exponentially in the direction perpendicular to the interface of both media [4]. It has become a popular research topic in the field of nanophotonics in recent years. Owing to their distinctive optical properties [5][6][7][8][9], SPPs can break through the optical diffraction limit and strengthen the near-field enhancement effect, providing an opportunity to explore the interaction between light and matter in micro/nano fields.
In the last few decades, owing to the potential applications of SPPs and excitons in sensors, light emitters and nano-optical devices, considerable efforts have been devoted to the coupling. SPPs represent a basic and interesting multimode coupling system. It is greatly affected by the structural parameters and can be adjusted independently. Strong coupling between SPPs and organic molecules have always been the subject of intense investigation [10] Essentially, SPPs can be considered as resonant-element-excited states [11][12][13], which are formed by the coherent coupling of quasi-free electron gas and the collective oscillation along the metal surface after the interaction between light and matter. Particularly, a new resonance, surface plasmon resonance (SPR), may be generated, whereas the resonance frequency in the metal SPPs is equal to the incident light [14][15][16]. It is worth noting that the frequency is affected by the size and morphology of the metal nanoparticles (NPs). On this basis, a variety of new composite structures can be produced by combining the metal structure with the material molecules. The study of the interaction between SPPs and molecules is of great significance not only in basic physics, but also in the research and development of various new nano-optoelectronic devices, as well as the development of quantum communication [17].
where ω A = ω A / ω F ≥ 0 is given in units of the frequency of the field, and g = g/ ω F is the (adimensional) coupling parameter. In addition, the operatorsâ andâ † are the one-mode annihilation and creation photon operators, respectively,Ĵ z ,Ĵ ± are the atomic relative population operator and atomic transition operators. Furthermore, the Dicke Hamiltonian model shows wide applications in many fields [58][59][60][61][62].
The famous Jaynes-Cummings (JC) model [63][64][65] is the simplest but nontrivial model for depicting the interaction of a two-level atom with an electromagnetic field at resonance, describing the radiation field and rotating wave approximation, additionally, it can be usually described as the most common soluble quantum-mechanical model of a single atom in an electromagnetic field, it represents the single-excitation sub-space of the full Jaynes-Cummings model: whereâ andâ † are also operators, which are the same as Equation (1), and ω pl represents the plasmon resonance frequency, while ω em is the emitter transition frequency, in addition, σ is the Pauli matrices for the two-level system. The decoherence and damping of the emitter and plasmon are neglected in this condition. The coupling strength between plasmon and emitter, g is as follows: where the transition dipole matrix element of emitter is → µ em and the local electric field generated by the plasmon at the position of the emitter is → ε 0 . For the electromagnetic fields which were confined into a cavity of mode volume V, the maximum g (for the emitter is at the maximum value of the field, whose dipole moment corresponds to the polarized direction of the field) is: Diagonalization of the Hamiltonian in Equation (2), the classical coupled harmonic oscillator model is often used to predict the upper polariton branches (UPB) and lower polariton branches (LPB) [66][67][68], it provides the mode frequencies of the polaritons: where the vacuum Rabi frequency Ω R is Ω R represents the frequency at which the energy oscillates between the plasmon and the emitter. Compared with the intrinsic linewidths of optical modes and excitons, Rabi splitting [69][70][71] is larger and more capable of maintaining a strong coupling state, which may lead to two independent modes: UPB and LPB, as shown in Figure 1a. While the two independent states are degenerated and owed equal lifetimes and effective masses, the strong coupling, which demonstrated anti-crossing and the appearance of two transitions of equal strength separated by Rabi splitting, the two new modes produced in the strong coupling state can be explained through mode hybridization, analogous to the molecular orbital theory [72,73]. In addition, the strong coupling regime illustrated by two coupled oscillators can be seen in Figure 1b [74], this system composed exciton A and plasmon. An important consequence of strong coupling is that the exciton should possess the plasmonic properties due to the coherent interactions, such as polarization and emission direction [75][76][77], and the exciton-plasmon coupling, including enhanced cooperative emission and enhanced exciton transport, in the excited state molecules and molecular assemblies interacting with optical modes localized in cavities and other nanostructures. A significant limitation of the single molecules strong-coupling is that the operation may be required under cryogenic temperature [78], and that the presence of nonlocality imposes, will limit strictly the practical applications [79].
where the vacuum Rabi frequency R Ω is R Ω represents the frequency at which the energy oscillates between the plasmon and the emitter. Compared with the intrinsic linewidths of optical modes and excitons, Rabi splitting [69][70][71] is larger and more capable of maintaining a strong coupling state, which may lead to two independent modes: UPB and LPB, as shown in Figure 1a. While the two independent states are degenerated and owed equal lifetimes and effective masses, the strong coupling, which demonstrated anti-crossing and the appearance of two transitions of equal strength separated by Rabi splitting, the two new modes produced in the strong coupling state can be explained through mode hybridization, analogous to the molecular orbital theory [72,73]. In addition, the strong coupling regime illustrated by two coupled oscillators can be seen in Figure 1b [74], this system composed exciton A and plasmon. An important consequence of strong coupling is that the exciton should possess the plasmonic properties due to the coherent interactions, such as polarization and emission direction [75][76][77], and the exciton-plasmon coupling, including enhanced cooperative emission and enhanced exciton transport, in the excited state molecules and molecular assemblies interacting with optical modes localized in cavities and other nanostructures. A significant limitation of the single molecules strong-coupling is that the operation may be required under cryogenic temperature [78], and that the presence of nonlocality imposes, will limit strictly the practical applications [79]. Furthermore, damping might be taken into consideration through the Heisenberg-Langevin method, in which both the plasmon and emitter are coupled to a library of harmonic-oscillator modes [80]. The Green's function approach [81] can be used to calculate the absorption spectra, while the Liouville approach [82] can be used to depict the system dynamics. Either the Liouville equations or the Heisenberg-Langevin [83,84] may both calculate the absorption and scattering cross-sections, through introducing the external electric field, ε , that acts on the plasmon:  Furthermore, damping might be taken into consideration through the Heisenberg-Langevin method, in which both the plasmon and emitter are coupled to a library of harmonicoscillator modes [80]. The Green's function approach [81] can be used to calculate the absorption spectra, while the Liouville approach [82] can be used to depict the system dynamics. Either the Liouville equations or the Heisenberg-Langevin [83,84] may both calculate the absorption and scattering cross-sections, through introducing the external electric field, ε, that acts on the plasmon: whereμ pl = d pl (â +â † ) is the plasmon dipole operator and d pl is the plasmon dipole moment. Owing to the interaction between the emitter and light being much weaker than that of the plasmonic nanostructure, the direct driving of the emitter by the external field can be ignored. Then, the optical response can be reached by solving the equations of motion in frequency or time domains, respectively. Several assumptions are involved in this method of analysis. By regarding the emitter as a two-level system, the complex energy level structure that exists in any quantum emitter is ignored. By regarding the plasmon field as a population of a single (quasi-) normal modes, additional modes that may overlap in frequency and any frequency dependence of the plasmon damping rate are ignored. Finally, it regards the emitter as a point dipole, not taking into account any effects due to magnetic field gradients in the spatial range Nanomaterials 2022, 12, 1242 5 of 25 of the emitter. Although the model is simplified, it can still process the experimental data quantitatively.
The solution of the quantum model requires large computational resources for accounting for all the potentially occupied plasmon substates. If the plasmon field can be approximated as the classical one, that is, significant acceleration can be obtained while the effect of quantum fluctuations in the field can be ignored, then in this condition, the Heisenberg-Langevin or Liouville equations may be reduced as the Maxwell-Bloch [85][86][87], which is semiclassical. Matthew's team [78] have verified that the predictions of the quantum model and the semiclassical model are almost identical in the case of parameters related to plasmon coupling, the Maxwell Bloch equation is as follows [88]: ..
where d em represents the emitter's transition dipole moment, γ pl is the plasmon decay rate, γ 1 is the emitter's energy decay rate, γ 2 is the emitter's dephasing rate, ρ em The pure dephasing of the plasmon can be ignored for the slowing energy damping. Differentiation of the first of Equation (9) with respect to time and substitution of the second equation into the first obtains ..
In the case of the linear response (scattering or absorption) of the coupled system, further approximations regarding the coherence between the emitter ground state and the excited state can also be ignored, and the emitter's population almost remains in the ground state; For example, ρ em 3 ≈ −1 and ρ em 2 1. Using µ em = d em ρ em 1 , we obtain: ..
Assuming the case where the dephasing of the emitter is much faster than energy decay, the total emitter decay rate is γ em ≈ γ 2 . This assumption is suitable for all experiments at or near room temperature. Equations (8) and (10) describe a pair of coupled classical harmonic oscillators. The displacement of each harmonic oscillator corresponds to the dipole moments of the plasmon and exciton, respectively. Generally, the electromagnetic near field couples via two oscillators. The plasmon oscillator is driven only by an external field. An intuitive and simple picture of this coupling exciton-emitter system is provided by the coupled oscillator model, and the oscillation of two coupled pendulums is easily solved, which can be found in reference [74,89].
In the absence of the external driving field, the coupled oscillator equation can be solved easily in the frequency domain to obtain the normal modal frequency of the system. If assuming the coupling strength, damping rate and detuning are smaller than resonance frequency (g, γ em , γ pl , ω em − ω pl ω pl ), and the normal-mode frequencies are: where now γ pl and γ em are the uncoupled SPPs and exciton decay rates, respectively. Comparison with Equations (5) and (6), shows that the results of the coupled-oscillator model are the same as the Jaynes-Cummings Hamiltonian under the constraints of zero damping and dephasing, if Ω R is real while the normal modes are non-degenerate. In the cases of ω em − ω pl ω pl , and γ em < γ pl (in all experiments that have been shown, both were fulfilled), this case reduces to Between over-damped and the under-damped cases, the transition can be expressed as a boundary between weak and strong coupling. However, strong coupling is considered to occur only when at least one complete Rabi oscillation occurs [90], when: Generally, many researchers use Equation (11) to judge whether the system is in the strong-coupling state or not. However, it is significant to realize that it does not represent a sharp threshold, when coupling plasmon-emitter systems cross this boundary, their behavior does not change qualitatively.
Moreover, the coupled oscillator equation can be easily generalized to N emitters independently coupled to a single plasmon mode, each with a different coupling constant g i . In the condition, the strong-coupling leads to where g i and where assuming the same linewidth, γ em , for all emitters. Up until now, the JC model proposed by Tavis and Cummings and its extension, the Tavis-Cummings (TC) model, are still the basis of studying the basic properties of quantum electrodynamics (QED) and understanding the existence of Rabi oscillations. To our knowledge, the SPPs-exciton strong coupling system is related to harmonic oscillators, and thus the definition of strong coupling is dependent on, among other things, the conventions of the particular field of physic involved; the system is in the strong coupling regime whenever the vacuum Rabi splitting (often called Rabi splitting) is experimentally observable [91].

Multiple Harmonic Oscillators
Furthermore, the multimode system is similar with multiple harmonic oscillators, considering the coupled three-harmonic-oscillator system as an example, as shown in Figure 2a. When coupling occurs, the frequencies of the oscillators will change, which is related to the coupling strength. Similarly, as shown in Figure 2b, when SPPs-double exciton coupling occurs, the plasmon energy level U pl , exciton A energy level U ea , and exciton B energy level U eb will change, and hybridization will occur, producing three hybridized plexciton states and double Rabi splitting. After the coupling, there are three unique solutions for U (U H , U M , and U L ), and thus the polariton energy dispersion comprises three energy branches. These are the aforementioned UPB and LPB, along with middle polariton branches (MPB). This change is evident in the scattering, absorption, and extinction spectra of the coupled system, with three peaks and two dips. When the plasmon resonance energy is coupled from low to high with excitons of constant energy, a double anticrossing image occurs, which is a typical characteristic of the multimode strong coupling between two excitons (excitons A and exciton B) and a plasmon, as illustrated in Figure 2c. comprises three energy branches. These are the aforementioned UPB and LPB, along with middle polariton branches (MPB). This change is evident in the scattering, absorption, and extinction spectra of the coupled system, with three peaks and two dips. When the plasmon resonance energy is coupled from low to high with excitons of constant energy, a double anticrossing image occurs, which is a typical characteristic of the multimode strong coupling between two excitons (excitons A and exciton B) and a plasmon, as illustrated in Figure 2c. To provide physical insights into plasmon-exciton multimode coupling, we further focus on the linear analysis of the scattering spectra, using the coupled oscillator model to fit the calculation results [93,94]. In the hybrid system, the coupling oscillators represent the plasmon, exciton A, and exciton B. The motion equations for the three oscillators are as follows [95,96]: x t x t x t g x t g x t F t γ ω To provide physical insights into plasmon-exciton multimode coupling, we further focus on the linear analysis of the scattering spectra, using the coupled oscillator model to fit the calculation results [93,94]. In the hybrid system, the coupling oscillators represent the plasmon, exciton A, and exciton B. The motion equations for the three oscillators are as follows [95,96]: where, x pl , x a and x b represent plasmon, exciton A oscillator, and exciton B oscillator, respectively. γ pl , γ a and γ b are the damping rate of plasmon, exciton A and exciton B, respectively, ω pl , ω a and ω b are the resonance frequencies of plasmon, exciton A and exciton B, respectively, g a is the coupling rate between the plasmon and exciton A, while g b is the coupling rate between the plasmon and exciton B, F pl , F a and F b represent the driving forces due to the external source. We assumed that both excitons A and B were both entirely driven by the plasmon oscillator; hence, we set F pl (t) = F pl e −iωt , where ω is the frequency of the electric field, F a (t) = 0, F b (t) = 0. Finally, x pl (t), x a (t) and x b (t) can be derived from Equations (16)- (18). In some specific hybrid nanostructures, the dimensions of the structures are small compared to the optical wavelength; thus, the scattering cross-section can be calculated in the quasi-limit [97]. In this limitation, the scattering cross-section is , where k = ωn/c is the wave vector of light. By substituting x pl (t) into (8π/3) · k 4 F pl x pl 2 and the use of the incident light energy E replaces the incident light frequency ω, the scattering cross-section can be obtained as follows: . The multimode coupling system can be explored by description of a three-coupled oscillator model composed of the optical modes [98,99], and the interaction Hamiltonian for the system can be expressed as [100]: are the undisturbed energies of the optical modes in excitons A and B, respectively, U a − i γ a 2 is the exciton energy, g b is the coupling constant between optical modes, and g a is the exciton-photon coupling potential. α pl , α b , and α a are the eigenvector components (Hopfield coefficients), the corresponding α pl 2 ,|a b | 2 , and |α a | 2 represent the weighting efficiencies and satisfy α pl . Thus, it is expected that there might exist three modes, namely the central mode under the unperturbed energy and the external two modes shifted by ± (g 2 a + g 2 b ) to a higher and lower energy. The eigenvectors Ψ of the UPB, MPB and LPB modes are [100][101][102]: where φ pl , φ b , and φ a are the unperturbed basis functions. With Equation (22), no component from the basis function φ pl is contained by the eigenvector of the central mode, of exciton A; for the central mode, the amplitude of the light field in the cavity containing this condition is zero. All weights of φ pl reside equally in the two external components. In comparison, if light is entering into the empty cavity, then the central mode will be bright (for the nonzero φ d component in Equation (22)). Surprisingly, the nature of the (3 × 3) matrix Hamiltonian of Equation (20) is shown, which has one basis function (φ pl ) that is coupled to the other two basic functions (φ b , φ a ). φ b and φ a by contrast are not coupled directly. The eigenvector Ψ MPB is the central mode. At the same time, these two external modes contain components from both φ pl and φ b , with thus observation in both directions.
This situation for the coupled two-cavity one-exciton case is similar to that found for the normal modes of three masses connected linearly by springs [103] or of linear triatomic molecules [104]. For which cases, there is one mode corresponding to the zero displacement of the central mass, and the two outer masses undergo tensile mode vibrations. This is consistent with Equation (22) where in the central mode the weight of the oscillator (φ pl ) connected with the other two is zero.

Scattering Spectra
For analyzing the coupling strength in the plasmonic system, which cannot be obtained directly from Fano's original spectrum, the coupled harmonic oscillator model served as a simulation of the scattering spectrum for describing the real physical situation, phenomeno-logically [95]. Thus, fitting of the calculated spectra into the coupled-oscillator model, which makes it possible to determine whether the system is in a strong-coupling regime or not, was performed for determining which phenomenon (hybridization or interference) is responsible for the spectra.
The coupled-oscillator equations can also be solved with a driving force to obtain the absorption and scattering cross sections. When the frequency of driving force is ω under the steady state, the solution is: The scattering cross-section is given by [95]: Similar conclusions can be obtained for absorption and extinction cross-sections [95]. From Equations (24) and (25), the calculated spectra indicate unique plasmon and exciton parameters and different coupling strengths.
In the previous literature [95,105], the following represent various coupled-oscillator equations with different coupling terms: ..
Generally, the coupling terms used in Equations (8) and (10) are preferred, since the equation for the leading choice of the coupling term is derived from a quantum mechanical model in the linear and classical limitations. However, if the coupling term in Equation (26) is used, the steady-state solution is different only in the last term of the denominator: This will give almost the same results as Equation (24), while g, γ em , γ pl ω; this has been the case for all reported experimental results [95,106,107].
Furthermore, a coupled harmonic oscillator model elucidates the nature of the interaction and reveals the coherent quantum superposition of the excitons, which are mediated by the plasmonic interaction. In addition, each component of the coupled system is described as a harmonic oscillator with its own resonance damping and frequency.
The theory discussed here is based on the characteristics of SPPs when coupled with excitons, which are confined to the metal surface at the nanometer scale, thus greatly compressing the spatial distribution of the electromagnetic field [108,109]. Furthermore, this spatial variation will introduce substantial theoretical complexities in the coupled oscillator models [77]. This forms a theoretical basis for exploring the coupling effect, clearly, the coupling and hybridization phenomena open an avenue for the enhancement of plasmonic resonances in nanostructures, which deserves further study. Two-dimensional transition metal dichalcogenides (TMDs) have recently attracted widespread attention owing to their unique electrical and optical properties. The TMDs display a strong enhancement of the Coulomb interaction [110]. The resulting bound electron-hole pairs, govern the electrical charge and optical transport properties, then affect the Coulomb interaction between the oscillator strength and exciton binding energy [111].

Recent Progress in Multimode
The exciton absorption of TMDs is located within the visible light band, showing broader application prospects and attracting the interest of researchers. TMDs have opened a new research challenge for strong coupling. Menon M. et al. [112] demonstrated the strong coupling between two-dimensional excitons and cavity photons in a monolayer MoS 2 -based microcavity with a Rabi splitting of 46 ± 3 meV for the first time. The same team [113] presented an approach for dynamically regulating the interaction between excitons in monolayer WS2 and microcavity photons at room temperature, and observed a Rabi splitting of 60 meV.
Furthermore, TMDs monolayers provide the opportunity to study multimode strong coupling with larger Rabi splitting, and temperature might be a key factor affecting the multiple coupling. Cuadra J. et al. [114] reported the strong interaction between localized surface plasmon resonance (LSPR) in silver nanoprisms and excitons and trions in monolayer tungsten disulfide (WS 2 ). The high density of the photonic states is shown in Figure 3a at the corners of the nanoprism, which overlaps with the WS 2 monolayer for efficient plasmon−exciton interaction, and the inset shows the SEM image of such a particle and a magnified view of a corresponding dark-field image. In addition, for T = 300 K, the Rabi splitting was 120 meV, whereas in a low temperature case, T = 6 K (Figure 3b, red, green, and blue solid lines represent the eigen energies extracted from the Hamiltonian analysis. Black solid lines indicate exciton, trion and plasmon resonances), three anticrossing bands, corresponding to UPB, MPB and LPB, respectively, with minimal splitting measured approximately 150 meV. The 30 meV increase can be explained by the detuned exciton and trion. Both plasmon and exciton resonances are narrowed upon cooling, and exciton line narrowing at low temperatures is consistent with previous results [115]. Figure 3c shows the Hopfield coefficients for the plasmon, exciton, and trion contributions to the UPB, MPB, and LPB states as a function of the plasmon resonance. Double Rabi splitting was observed in the dark-field scattering spectra using a plasmonic nanostructure interacting with two types of excitons (charged and neutral) in WS 2 with decreasing temperature. The degree of plasmon-exciton-trion coupling might be tailored by varying the temperature.
Moreover, similar exploration was conducted by other groups, and the size of the nanosystem was also related to the strength of the multimode coupling. Li B. et al. [116] demonstrated a large 300 meV Rabi splitting that was achieved under ambient conditions in a strong coupling regime by Ag-WS 2 heterostructure embedding. The system consisted of a 100 nm Ag mirror and a 30 nm MgF 2 spacer with the Ag nanodisks fabricated by e-beam lithography (EBL), as shown in Figure 4a. In the plasmon-exciton-cavity, the 99 nm Ag nanodisk was the best structure for the strongest interaction, with a 300 meV Rabi splitting among the three oscillators, which can be seen in Figure 4b. Red dots with error bars show energies obtained from the reflection spectrum. The horizontal black dashed lines represent the A-exciton and the microcavity resonant energy, respectively. The black slanted shortdashed line represents plasmon resonance mode. It depicts resulting dispersion curves of the Ag-WS 2 heterostructure that was embedded in the optical microcavity. They proposed the criteria of the strong coupling for the three oscillators, and the Hopfield coefficients of the three branches were calculated as shown in Figure 4c. Each branch consisted of part-plasmon, part-exciton, and part-cavity modes, demonstrating that the hybrid state was successfully generated by combining the Ag-WS 2 heterostructure with the optical microcavity. At a nanodisk diameter of 99 nm, the strongest interaction of 300 meV Rabi splitting was obtained. Thus, it is clearly evident that the size of the Ag nanodisk affects the strong coupling.
citon line narrowing at low temperatures is consistent with previous results [115]. Figure  3c shows the Hopfield coefficients for the plasmon, exciton, and trion contributions to the UPB, MPB, and LPB states as a function of the plasmon resonance. Double Rabi splitting was observed in the dark-field scattering spectra using a plasmonic nanostructure interacting with two types of excitons (charged and neutral) in WS2 with decreasing temperature. The degree of plasmon-exciton-trion coupling might be tailored by varying the temperature. Moreover, similar exploration was conducted by other groups, and the size of the nanosystem was also related to the strength of the multimode coupling. Li B. et al. [116] demonstrated a large 300 meV Rabi splitting that was achieved under ambient conditions in a strong coupling regime by Ag-WS2 heterostructure embedding. The system consisted  Figure 4a. In the plasmon-exciton-cavity, the 99 nm Ag nanodisk was the best structure for the strongest interaction, with a 300 meV Rabi splitting among the three oscillators, which can be seen in Figure 4b. Red dots with error bars show energies obtained from the reflection spectrum. The horizontal black dashed lines represent the A-exciton and the microcavity resonant energy, respectively. The black slanted short-dashed line represents plasmon resonance mode. It depicts resulting dispersion curves of the Ag-WS2 heterostructure that was embedded in the optical microcavity. They proposed the criteria of the strong coupling for the three oscillators, and the Hopfield coefficients of the three branches were calculated as shown in Figure 4c. Each branch consisted of part-plasmon, part-exciton, and part-cavity modes, demonstrating that the hybrid state was successfully generated by combining the Ag-WS2 heterostructure with the optical microcavity. At a nanodisk diameter of 99 nm, the strongest interaction of 300 meV Rabi splitting was obtained. Thus, it is clearly evident that the size of the Ag nanodisk affects the strong coupling. Jiang P. et al. [117] proposed a strong exciton-plasmon-exciton coupling system that consisted of a silver nanoprism separated from a monolayer WS2 by J-aggregates, as depicted in Figure 5a, which was performed by finite-difference time-domain (FDTD) simulations to obtain the optical response. Figure 5b depicts the scattering cross-section spectra of the hybrid nanostructure. The solid white lines indicate the three anticrossing bands corresponding to the UPB, MPB and LPB, and the solid white lines match well with the numerical simulation results; the splitting extracted between the UPB and MPB is 130 meV at the zero detuning while the splitting between the MPB and LPB is 170 meV. Since the damping losses were not considered, the calculated Rabi splittings are larger than the actual results slightly. Furthermore, Figure 5c shows the weighting efficiencies for the LSPR mode, exciton A of WS2, and J-aggregates exciton contributions to three polariton Jiang P. et al. [117] proposed a strong exciton-plasmon-exciton coupling system that consisted of a silver nanoprism separated from a monolayer WS 2 by J-aggregates, as depicted in Figure 5a, which was performed by finite-difference time-domain (FDTD) simulations to obtain the optical response. Figure 5b depicts the scattering cross-section spectra of the hybrid nanostructure. The solid white lines indicate the three anticrossing bands corresponding to the UPB, MPB and LPB, and the solid white lines match well with the numerical simulation results; the splitting extracted between the UPB and MPB is 130 meV at the zero detuning while the splitting between the MPB and LPB is 170 meV. Since the damping losses were not considered, the calculated Rabi splittings are larger than the actual results slightly. Furthermore, Figure 5c shows the weighting efficiencies for the LSPR mode, exciton A of WS 2 , and J-aggregates exciton contributions to three polariton branches as a function of the plasmon resonance. Overall, the research interest in TMDs mainly focuses on their distinctive optical properties: strong nonlinear optics response, effective valley-spin coupling, strong exciton binding energy, large spin-orbit interaction, and transition dipole momentum of the material. These properties of multimode strong coupling systems based on TMDs and their related characteristics are now sought after by researchers. Such a system has potential applications in optical modulators at the nanoscale and polaritonic devices based on ultrathin materials.

Coupling Related to Single-Dye Microcavity Nanostructure
Although most studies about strong coupling have been explored interactions between a single excitonic mode and a single photonic mode, some teams have focused on multimode mixed coupling, in which the excitonic mode is coupled to more than one photonic mode; however, these systems are complicated, and the quality of excitons are more or less affected.
Zhang K. et al. [118] presented a hybrid strong coupling between multiple photonic modes and excitons in an organic-dye-attached photonic quasicrystal. The electron beam evaporation method was used to deposit the dielectric multilayers, and TDBC dye was selected for the J-aggregates, acting as an organic semiconductor. Following the simulation parameters, SiO2∕Ta2O5 multilayers were deposited onto K9 glass substrate, and then the J-aggregates were directly spun onto the top surface, as schematically described in Figure 6a. The value of the coupling energy, namely the coupling energy between the photonic mode and the excitons, was found to be Overall, the research interest in TMDs mainly focuses on their distinctive optical properties: strong nonlinear optics response, effective valley-spin coupling, strong exciton binding energy, large spin-orbit interaction, and transition dipole momentum of the material. These properties of multimode strong coupling systems based on TMDs and their related characteristics are now sought after by researchers. Such a system has potential applications in optical modulators at the nanoscale and polaritonic devices based on ultrathin materials.

Coupling Related to Single-Dye Microcavity Nanostructure
Although most studies about strong coupling have been explored interactions between a single excitonic mode and a single photonic mode, some teams have focused on multimode mixed coupling, in which the excitonic mode is coupled to more than one photonic mode; however, these systems are complicated, and the quality of excitons are more or less affected.
Zhang K. et al. [118] presented a hybrid strong coupling between multiple photonic modes and excitons in an organic-dye-attached photonic quasicrystal. The electron beam evaporation method was used to deposit the dielectric multilayers, and TDBC dye was selected for the J-aggregates, acting as an organic semiconductor. Following the simulation parameters, SiO 2 /Ta 2 O 5 multilayers were deposited onto K9 glass substrate, and then the J-aggregates were directly spun onto the top surface, as schematically described in Figure 6a. The value of the coupling energy, namely the coupling energy between the photonic mode and the excitons, was found to be hΩ H = 67 meV, hΩ C = 93 meV and hΩ L = 85 meV. The six polariton bands of the calculated eigen energies in this system are described by the dashed lines in Figure 6b, and all these results fit well with the experimental data.  Figure 6b, and all these results fit well with the experimental data. Balasubrahmaniyam M. et al. [119] investigated and experimentally implemented the photonic analog of localization induced by ultra-strong interactions in a coupled threemode system. Furthermore, they demonstrated the ultra-strong coupling between a highly dispersive cavity plasmon mode and dimer excitons of J-aggregates Rhodamine B (RhB). A prism arrangement in a photon tunneling configuration was used for probing the dispersion, and two prisms (SF11, refractive index ~ 1.78) with 40 nm thick Ag films coated over their hypotenuse faces were brought close to each other to form the cavity, as illustrated in Figure 7a. Moreover, the position of the RhB in the cavity results in the formation of three antisymmetric plasmon (ASP) mode-exciton hybrid modes and splits the measured ASP transmission peak into three. Figure 7b shows the position of these peaks, where a permanent broad gap of 420 meV coupling is obtained, confirming the ultrastrong coupling between the exciton levels and the ASP mode.  Balasubrahmaniyam M. et al. [119] investigated and experimentally implemented the photonic analog of localization induced by ultra-strong interactions in a coupled three-mode system. Furthermore, they demonstrated the ultra-strong coupling between a highly dispersive cavity plasmon mode and dimer excitons of J-aggregates Rhodamine B (RhB). A prism arrangement in a photon tunneling configuration was used for probing the dispersion, and two prisms (SF11, refractive index~1.78) with 40 nm thick Ag films coated over their hypotenuse faces were brought close to each other to form the cavity, as illustrated in Figure 7a. Moreover, the position of the RhB in the cavity results in the formation of three antisymmetric plasmon (ASP) mode-exciton hybrid modes and splits the measured ASP transmission peak into three. Figure 7b shows the position of these peaks, where a permanent broad gap of 420 meV coupling is obtained, confirming the ultra-strong coupling between the exciton levels and the ASP mode.  Figure 6b, and all these results fit well with the experimental data. Balasubrahmaniyam M. et al. [119] investigated and experimentally implemented the photonic analog of localization induced by ultra-strong interactions in a coupled threemode system. Furthermore, they demonstrated the ultra-strong coupling between a highly dispersive cavity plasmon mode and dimer excitons of J-aggregates Rhodamine B (RhB). A prism arrangement in a photon tunneling configuration was used for probing the dispersion, and two prisms (SF11, refractive index ~ 1.78) with 40 nm thick Ag films coated over their hypotenuse faces were brought close to each other to form the cavity, as illustrated in Figure 7a. Moreover, the position of the RhB in the cavity results in the formation of three antisymmetric plasmon (ASP) mode-exciton hybrid modes and splits the measured ASP transmission peak into three. Figure 7b shows the position of these peaks, where a permanent broad gap of 420 meV coupling is obtained, confirming the ultrastrong coupling between the exciton levels and the ASP mode.  Hakala T. et al. [120] demonstrated the emission of three energy branches in the strong coupling between the surface plasmon polaritons and J-aggregates Rhodamine 6G (R6G) sandwich structure, with a glass substrate, a 45 nm silver film in the middle, and a resist layer on top, with double vacuum Rabi splitting energies up to 230 and 110 meV for the 200 mM sample R6G concentrations. It is the first time that in the strongly coupled SPPmolecule system, each of the three energy branches can be converted into photons, and the finding shows great foreground in the field of multimode energy transfer and hybridization.
In all, the one-dye microcavity nanostructure system can be used to provide multimode hybrid interactions, which may inspire related exploring on multimode light-matter interactions and achieving some potential applications, such as multimode sensors and spectroscopy.

Coupling Related to Two-Dye Microcavity Nanostructure
It is worth noting that current research leans towards more than two types of light components to form hybrid polaritons, one with strong dispersion and the other one nearly nondispersive, which can inspire related studies on hybrid light-matter interactions. To date, interactions between electronic states in organic states and plasmons have been undertaken intensively. In particular, when the emitter which is near the metal surface with the plasmon resonance, the spontaneous emission rate is modulated and the rate of energy transfer is changed [121,122]. However, most studies are conducted under the weak coupling region. On the other hand, pioneering studies have demonstrated that organic dye molecules with high oscillator strength form strong coupling region under plasmon excitation [123,124]. The coupling strength can be controlled through the plasmon energy and its line width, suggesting usability of a wide variety of organic molecules for achieving strong coupling.
Lidzey D. et al. [125] presented microcavities that can occur between the two cyanine dyes, with simultaneous strong coupling of the excitations of the individual dyes to a single cavity mode. The cyanine dye layers were spatially separated by a 100 nm-thick barrier layer of polystyrene to give an exciton separation of 140 nm. The cavity contained a dye concentration of approximately 2.3 × 10 20 cm −3 , where the 37 meV and 58 meV splitting branches were observed, whereas with a lower dye concentration of 1.2 × 10 20 cm −3 , the splitting branches were 18 meV and 44 meV. The exciton scattering rates were expected to be significantly enhanced because Frenkel excitons have large interaction cross-sections with molecular vibrations.
Melnikau D. et al. [126] introduced hybrid structures consisting of Au nanostars and Jaggregates of cyanine dyes, and Rabi splitting with an energy of up to 260 meV. In addition, the absorption spectrum of the complex hybrid system showed two pronounced dips at 590 and 642 nm (red curve), which corresponded to the maximum absorption wavelengths of the two different dyes (JC1 and S2165), and double Rabi splitting with the energies of 187 and 119 meV, as shown in Figure 8a. In addition, Figure 8b is the results from the model simulations [127] corroborated the experimental findings, when the position of the exciton resonance shifts to red or blue, which is relative to the maximum absorption of the nanostar, a unique asymmetric profile can be seen in the spectrum of the hybrid system.
Coles D. et al. [128] used strong coupling in an optical microcavity to mix the electronic transitions of two J-aggregated molecular dyes (TDBC and NK-2707) with a 200-nm-thick silver mirror, and used both non-resonant photoluminescence emission and photoluminescence excitation spectroscopy to show that hybrid polariton states act as an efficient and ultrafast energy-transfer pathway between the two exciton states. At normal incidence, the photon energy was 101 meV below the NK-2707 exciton (hereafter referred to ex1 for simplicity) and 250 meV below the TDBC exciton (hereafter ex2). A Rabi splitting energy of hΩ 1 = 73 meV was determined between the LPB and MPB, whereas there was a larger splitting of hΩ 2 = 155 meV between the MPB and UPB, consistent with the higher oscillator strength contributed from the TDBC J-aggregates. Coles D. et al. [128] used strong coupling in an optical microcavity to mix the electronic transitions of two J-aggregated molecular dyes (TDBC and NK-2707) with a 200nm-thick silver mirror, and used both non-resonant photoluminescence emission and photoluminescence excitation spectroscopy to show that hybrid polariton states act as an efficient and ultrafast energy-transfer pathway between the two exciton states. At normal incidence, the photon energy was 101 meV below the NK-2707 exciton (hereafter referred to ex1 for simplicity) and 250 meV below the TDBC exciton (hereafter ex2). A Rabi splitting energy of ℎΩ = 73 meV was determined between the LPB and MPB, whereas there was a larger splitting of ℎΩ = 155 meV between the MPB and UPB, consistent with the higher oscillator strength contributed from the TDBC J-aggregates.
It is worth noting that, apart from the interaction area length, the two-dyes microcavity nanostructure can increase the SPP-material interaction by creating a layer at the top of the interacting area, thereby decreasing the mode volume and preventing decay into the radiative mode, in which the upper exciton reservoir alone could be excited. This shows potential applications in integrated microcavity sensors and optical devices.

Multimode Coupling Related to Periodic Noble Metallic NPs-J-Aggregates Nanostructure
Organic semiconductor materials, such as J-aggregate [129][130][131] dyes, offer great potential because of their advantages in low-cost manufacturing and optical properties. Jaggregates support exciton states, which are electrically neutral electron/hole pairs generated by the absorption of photons. Excitonic states also exhibit strong nonlinear optical behavior that can be used to generate excitation sources for photon and transistor effects.
In addition, the periodic noble metallic nanoparticles (NPs) can be adjusted to modulate the plasmonic response of the nanostructure, therefore, the new periodic NPs monolayer nanostructure becomes a prime candidate for researchers exploring strong coupling phenomena.
Zhang K. et al. [132] demonstrated a nanostructured cavity with hybrid coupling among the molecular excitons, surface plasmon polaritons, and Fabry-Perot mode, where a J-aggregate-doped polyvinyl alcohol (PVA) layer was inserted between a silver grating and a thick silver film. They designed a nanostructure PVA with a thickness of 170 nm, inserted between a silver grating of 300 nm period and a 40 nm slit width, as can be seen in Figure 9a, the SPPs energy can be tuned through changing the period of the grating. In addition, for introduction of the excitons, the pure PVA was replaced with cyanine dye TDBC, as shown in Figure 9b. They observed three dispersive polariton bands with Rabi splitting energies of ℎΩ = 110.6 meV and ℎΩ = 122.5 meV, as shown in Figure 9c, indicating the strong coupling between the SPPs mode and the PB1 or PB2 mode (there are Figure 8. (a) Absorption spectra of gold nanostars (balck curve), pristine J-aggregates of JC1 (magenta curve) and S2165 (blue curve), and the hybrid structure (red curve). (b) Theoretical extinction spectra of gold nanostars (black curve) and their hybrid structure with J-aggregates (red curve). Adapted with permission from ref. [126]. Copyright 2013 IEEE.
It is worth noting that, apart from the interaction area length, the two-dyes microcavity nanostructure can increase the SPP-material interaction by creating a layer at the top of the interacting area, thereby decreasing the mode volume and preventing decay into the radiative mode, in which the upper exciton reservoir alone could be excited. This shows potential applications in integrated microcavity sensors and optical devices.

Multimode Coupling Related to Periodic Noble Metallic NPs-J-Aggregates Nanostructure
Organic semiconductor materials, such as J-aggregate [129][130][131] dyes, offer great potential because of their advantages in low-cost manufacturing and optical properties. J-aggregates support exciton states, which are electrically neutral electron/hole pairs generated by the absorption of photons. Excitonic states also exhibit strong nonlinear optical behavior that can be used to generate excitation sources for photon and transistor effects.
In addition, the periodic noble metallic nanoparticles (NPs) can be adjusted to modulate the plasmonic response of the nanostructure, therefore, the new periodic NPs monolayer nanostructure becomes a prime candidate for researchers exploring strong coupling phenomena.
Zhang K. et al. [132] demonstrated a nanostructured cavity with hybrid coupling among the molecular excitons, surface plasmon polaritons, and Fabry-Perot mode, where a J-aggregate-doped polyvinyl alcohol (PVA) layer was inserted between a silver grating and a thick silver film. They designed a nanostructure PVA with a thickness of 170 nm, inserted between a silver grating of 300 nm period and a 40 nm slit width, as can be seen in Figure 9a, the SPPs energy can be tuned through changing the period of the grating. In addition, for introduction of the excitons, the pure PVA was replaced with cyanine dye TDBC, as shown in Figure 9b. They observed three dispersive polariton bands with Rabi splitting energies of hΩ 1 = 110.6 meV and hΩ 2 = 122.5 meV, as shown in Figure 9c, indicating the strong coupling between the SPPs mode and the PB1 or PB2 mode (there are two unique solutions, thus two hybrid modes emerging). Due to the non-dispersion of the Fabry-Perot modes and the exciton, these mixed modes hardly disperse with the plane wave vector, as shown in Figure 9c, the two modes PB1 and PB2 were labeled. two unique solutions, thus two hybrid modes emerging). Due to the non-dispersion of the Fabry-Perot modes and the exciton, these mixed modes hardly disperse with the plane wave vector, as shown in Figure 9c, the two modes PB1 and PB2 were labeled. Wang H. et al. [133] constructed a hybrid system consisting of a 200 nm thick Au film patterned with a 2D square lattice gold nanohole array, and three different samples were prepared by spin-coating a uniform layer (300 nm thick) of sulforhodamine 101 (SR101) films with different concentrations, as can be seen in Figure 10a. Furthermore, by increasing the SR101 concentration, the Rabi splitting widened. The double Rabi splitting energies were 255 and 188 meV in steady-state transmission measurements, as illustrated in Figure 10b, and the dispersion was in excellent agreement with the typical signature of strong coupling. When the gold nanohole period is 380 nm, which matches the absorption of SR101 of the SPP resonance, the coupling strength maximizes, and when the plasmon mode periodicity is 350 nm, which matches the shoulder absorption of SR101, the double Rabi splitting is observable for the two periods. Yang H. et al. [134] investigated plasmon−exciton−plasmon (PEP) couplings in Ag−Jaggregates−Ag (AJA) nanostructures using the finite element method, which can be seen in Figure 11a. The grating period is 500 nm, which is not excessively large so as to superpose the lattice resonance with the LSPR. The three new hybridized plexciton modes were obtained by adjusting the geometry and the incident angle. As shown in Figure 11b, the Wang H. et al. [133] constructed a hybrid system consisting of a 200 nm thick Au film patterned with a 2D square lattice gold nanohole array, and three different samples were prepared by spin-coating a uniform layer (300 nm thick) of sulforhodamine 101 (SR101) films with different concentrations, as can be seen in Figure 10a. Furthermore, by increasing the SR101 concentration, the Rabi splitting widened. The double Rabi splitting energies were 255 and 188 meV in steady-state transmission measurements, as illustrated in Figure 10b, and the dispersion was in excellent agreement with the typical signature of strong coupling. When the gold nanohole period is 380 nm, which matches the absorption of SR101 of the SPP resonance, the coupling strength maximizes, and when the plasmon mode periodicity is 350 nm, which matches the shoulder absorption of SR101, the double Rabi splitting is observable for the two periods. two unique solutions, thus two hybrid modes emerging). Due to the non-dispersion of the Fabry-Perot modes and the exciton, these mixed modes hardly disperse with the plane wave vector, as shown in Figure 9c, the two modes PB1 and PB2 were labeled. Wang H. et al. [133] constructed a hybrid system consisting of a 200 nm thick Au film patterned with a 2D square lattice gold nanohole array, and three different samples were prepared by spin-coating a uniform layer (300 nm thick) of sulforhodamine 101 (SR101) films with different concentrations, as can be seen in Figure 10a. Furthermore, by increasing the SR101 concentration, the Rabi splitting widened. The double Rabi splitting energies were 255 and 188 meV in steady-state transmission measurements, as illustrated in Figure 10b, and the dispersion was in excellent agreement with the typical signature of strong coupling. When the gold nanohole period is 380 nm, which matches the absorption of SR101 of the SPP resonance, the coupling strength maximizes, and when the plasmon mode periodicity is 350 nm, which matches the shoulder absorption of SR101, the double Rabi splitting is observable for the two periods. Yang H. et al. [134] investigated plasmon−exciton−plasmon (PEP) couplings in Ag−Jaggregates−Ag (AJA) nanostructures using the finite element method, which can be seen in Figure 11a. The grating period is 500 nm, which is not excessively large so as to superpose the lattice resonance with the LSPR. The three new hybridized plexciton modes were obtained by adjusting the geometry and the incident angle. As shown in Figure 11b, the Yang H. et al. [134] investigated plasmon−exciton−plasmon (PEP) couplings in Ag−Jaggregates−Ag (AJA) nanostructures using the finite element method, which can be seen in Figure 11a. The grating period is 500 nm, which is not excessively large so as to superpose the lattice resonance with the LSPR. The three new hybridized plexciton modes were obtained by adjusting the geometry and the incident angle. As shown in Figure 11b  Li M. et al. [135] demonstrated a classic oscillator coupled model consisting of a 300 nm thick SiO2 substrate and a 200 nm thick silver film with 100 nm diameter periodic hole arrays. The molecular layer doped with polyvinyl alcohol (PVA) to enhance the localized SPP electric field was deposited on the 40 nm thick top layer and in the holes, as illustrated in Figure 12a. Different periods of theoretical curves are plotted and was in good agreement with simulated results, which can be seen in Figure 12b. In the lower polaritons band, the slight inconsistency is due to the dominating SPP modes in the hybrid system. Furthermore, the detailed formation of a strong coupling was proposed based on three states: First, the interaction between the SPP (1,0) mode with molecular excitons generates the first Rabi splitting, and then with the adjustment of the system structure, the interaction between SPP (1,1) mode and the excitons generates the second anticrossing effect. It is the double Rabi splitting that draws attention to exploring the enhanced strong coupling effects. A maximum effective coupling strength of 0.316 and Rabi splitting value of 663 meV.  Figure 12a. Different periods of theoretical curves are plotted and was in good agreement with simulated results, which can be seen in Figure 12b. In the lower polaritons band, the slight inconsistency is due to the dominating SPP modes in the hybrid system. Furthermore, the detailed formation of a strong coupling was proposed based on three states: First, the interaction between the SPP (1,0) mode with molecular excitons generates the first Rabi splitting, and then with the adjustment of the system structure, the interaction between SPP (1,1) mode and the excitons generates the second anticrossing effect. It is the double Rabi splitting that draws attention to exploring the enhanced strong coupling effects. A maximum effective coupling strength of 0.316 and Rabi splitting value of 663 meV. Evidently, both the selection of this type of J-aggregated molecular and grating period were dictated by the need to realize strong coupling, motivated by the exceptionally high oscillator strength and narrow resonances, even at room temperature. The corresponding work may inspire related studies on hybrid light-matter interactions and Evidently, both the selection of this type of J-aggregated molecular and grating period were dictated by the need to realize strong coupling, motivated by the exceptionally high oscillator strength and narrow resonances, even at room temperature. The corresponding work may inspire related studies on hybrid light-matter interactions and achieve potential applications in multimode lasers and optical micro-spectroscopy.

Multimode Coupling Related to Core-Shell Multicomponent Systems
Various interesting phenomena occur if the coupling between the cavity and the emitter is sufficiently large such that the characteristic interaction time exceeds all other decay channels. In this strong coupling regime, photon and matter excitations hybridize to form a new type of quasi particle. These polaritons are revealed by two new energy-shifted optical resonances. They possess both light-and matter-like properties, which make them suitable for a wide range of applications. Here, the controlled strong coupling of a single emitter to a cavity allows quantum-state mapping of a localized emitter qubit to a moving photon qubit. For this type of application, it is essential to achieve strong coupling of a single emitter.
Melnikau D. et al. [92] demonstrated a single hybrid structure with multiple spectral features, which was integrated with three different components, including core-shell Au@Ag nanorods (NRs), and two different dyes (5,5 ,6,6 -tetra-chloro-1,1 ,3,3 -tetraethylimidacarbocyanine iodide(TCI) and 2-[3-[1,1-dimethyl-3-(4-sulfobutyl)-1,3-dihydro-benzo[e]indol-2-ylidene]-propenyl]-1,1-dimethyl-3-(4-sulfobutyl)-1H-benzo[e]indolium hydroxide, (DBI)), as shown in Figure 13a, leading to a strong collective exciton-plasmon coupling. Figure 13b displays the theoretical results (solid lines) for the Au@Ag NRs with the two dyes fitted to the experimental data obtained for the hybrid plexciton states in the case of double Rabi splitting, and the estimated Rabi splitting values are 175 and 163 meV. All three components were involved in the strong interaction, and the total value of 338 meV reflects the onset of collective extended energy splitting between the states of the lower and upper plexcitons. Despite the fact that microcavities are considered to be an excellent alternative to plasmonic nanostructures for achieving strong coupling, the core-shell multicomponent systems show magneto-optical activity for hybridized plexciton states resulting from strong exciton-plasmon coupling in nanostructures, in which the plasmonic component induces magnetic properties in nonmagnetic organic fluorophores. This is significant for the development of new sensing systems based on magneto-optical activity, while expanding the portfolio of materials that can be used for optical information storage and processing.
Different multimode coupling with their applications is shown in Table 1.  Despite the fact that microcavities are considered to be an excellent alternative to plasmonic nanostructures for achieving strong coupling, the core-shell multicomponent systems show magneto-optical activity for hybridized plexciton states resulting from strong exciton-plasmon coupling in nanostructures, in which the plasmonic component induces magnetic properties in nonmagnetic organic fluorophores. This is significant for the development of new sensing systems based on magneto-optical activity, while expanding the portfolio of materials that can be used for optical information storage and processing.
Different multimode coupling with their applications is shown in Table 1. Multimode coupling related to core-shell multicomponent systems Nanoscale optical information storage and processing, etc.

Conclusions and Development Trends
In general, surface plasmons have become a popular topic in the field of photonics owing to their distinctive optical properties [136,137]. This is particularly evident on the upper surface in studying the coupling between surface plasmons and excitons. This paper reviewed the research progress in the field based on multimode coupling related to TMDs microcavity nanostructure, single-dye microcavity nanostructure, two-dyes microcavity nanostructure and core-shell multicomponent systems. This paper is introduced from a theoretical point of view. It has been shown that multimode coupling can not only reduce the size of quantum devices, but also has preliminarily proven the superiority to regulating the interaction between light and matter.
However, there remain numerous unknown problems regarding multimode coupling based on surface plasma nanostructures, which should be explored further in terms of the following aspects: (1) Strong coupling (the signature is a single Rabi splitting) can induce photon electricity and photon thermal processes, whereas multimode coupling is more complicated than general coupling, and may lead to a new series of energy conversion processes, such as photon-thermal-electricity and photon-electricity-thermal processes induced by surface plasmons. In addition, it provides a new route for further exploration of the physical mechanism of surface plasmons, which bears potential scientific applications for the development of new nano-optical devices.
(2) Strong coupling can be demonstrated including ultrafast tuning of strongly coupled metal-molecular aggregates via femtosecond pumping [138], dynamical modification of the polariton composition [139] and UV illumination of silver nanoparticle arrays [140]. However, there are few studies in this field of new physics, which concern potential phase changes, such as structure-induced phase changes and temperature-induced phase changes. Moreover, it may provide an exciting topic for future research.
(3) With the help of a database, we discovered that reports based on multimode splitting are limited, and that single Rabi splitting is the mainstream in current research fields. If the conditions are appropriately changed, future research hotspots would focus on whether a higher-mode process could occur.
(4) In practical applications, research into quantum information processing, higherorder processes, and nonlinear optics between surface plasmons and matter, are expected to be further discussed in terms of the physical essence of interactions.
(5) As an emerging discipline, machine learning applications in plasmonics have attracted considerable attention [141][142][143], such as both the forward prediction of far-field optical properties and the inverse prediction of on-demand dimensional parameters of nanoparticles [144], and the prediction of color of nano-structured surfaces through either the nanoparticle geometric parameters, or laser parameters [145]. However, there are few examples by using machine learning to study the multimode strong coupling, which may be one of the future focuses.
The multimode coupling of plasmons and excitons based on metal nanostructures have greatly promoted the development of strong coupling. In particular, it contains considerable potential for application in the design and improvement of organic photoelectric devices such as organic luminescence and organic batteries [146,147]. Therefore, surface plasma optics, produced by the combination of quantum optics and surface plasmon optics, may serve as a new development direction. This will not only increase the depth of research regarding the basic properties of surface plasmon excitons but also provide directions for solving current problems in quantum optics research.

Data Availability Statement:
The study did not report any data.