Analogue of double-$\Lambda$-type atomic medium and vector dromions in a plasmonic metamaterial

We consider an array of the meta-atom consisting of two cut-wires and a split-ring resonator interacting with an electromagnetic field with two polarization components. We prove that such metamaterial system can be taken as a classical analogue of an atomic medium with a double- $\Lambda$-type four-level configuration coupled with four laser fields, exhibits an effect of plasmon induced transparency (PIT), and displays a similar behavior of atomic four-wave mixing (FWM). We demonstrate that when nonlinear varactors are mounted onto the gaps of the split-ring resonators the system can acquire giant second- and third-order Kerr nonlinearities via the PIT and a longwave-shortwave interaction. We also demonstrate that the system supports high-dimensional vector plasmonic dromions [i.e. (2+1)-dimensional plasmonic solitons with two polarization components, each of which has a coupling between a longwave and a shortwave], which have very low generation power and are robust during propagation. Our work gives not only a plasmonic analogue of the FWM in coherent atomic systems but also provides the possibility for obtaining new type of nonlinear polaritons in plasmonic metamaterials.


I. INTRODUCTION
Electromagnetically induced transparency (EIT), a very intriguing phenomenon occurring in atomic gases, has been intensively investigated due to its interesting physical properties and promising practical applications. The basic mechanism of EIT is the existence of a destructive quantum interference effect between two pathways of atomic transitions induced by a control laser field, through which the absorption of a probe laser field can be largely cancelled [1].
Especially, the plasmonic analogue of EIT, called plasmon-induced transparency (PIT) [9][10][11], has become a very important platform for exploring EIT-like physical properties of plasmonic polaritons and for designing new types of metematerials.
Similar to EIT, PIT is resulted from a destructive interference effect between wideband bright and narrowband dark modes in artificial atoms (called meta-atoms). A typical character of PIT is the opening of a deep transparency window within broadband absorption spectrum, together with a steep dispersion and greatly reduced group velocity of plasmonic polaritons. PIT metamaterials can work in different frequency regions (including micro [10] and terahertz [11,14,17] waves, infrared and visible radiations [9,12,16]), and may be used to design novel, chip-scale plasmonic devices (including highly sensitive sensors [13,14], optical buffers [15,17], and ultrafast optical switches [17], etc.) in which the radiation damping can be significantly eliminated, very intriguing for practical applications.
In this article, we give a positive answer for the above question. The metamaterial we consider is assumed to be an array of meta-atoms [see Fig. 2(a)], i.e. the unit cells consisting of two cut-wires (CWs) and a split-ring resonator (SRR) [see Fig. 2(b)], interacting with an electromagnetic (EM) field with two polarization components. We show that such plasmonic metamaterial system may be taken as a classical analogue of an atomic medium with a double-Λ-type four-level configuration coupled with four (two probe and two control) laser fields [see Fig. 1(a) ], exhibits an effect of PIT and displays a similar behavior of atomic four-wave mixing (FWM).
Based on this classical analogue, we further show that, if nonlinear varactors are mounted onto the gaps of the SRRs, the system can acquire giant second-and third-order Kerr nonlinearities many orders of magnitude larger than conventional nonlinear optical media.
Using a method of multiple scales, we derive coupled envelope equations, which include dispersion, diffraction, and the Kerr nonlinearities and govern the evolution of the two polarization components of the EM field. We demonstrate that the system supports a new type of nonlinear plasmonic polaritons, i.e. high-dimensional vector plasmonic dromions, which are (2+1)-dimensional plasmonic solitons with two polarization components. Each polarization component has a coupling between a longwave and a shortwave, which have very low generation power and are robust during propagation. The results presented here not only gives a close metamaterial analogue of the EIT and FWM in multi-level atomic systems, useful to illustrate and find novel interference and nonlinear properties in solid systems, but also provides a way to obtain new type of plasmonic polaritons via suitable design of plasmonic metamaterials.
The main body of the article is arranged as follows. In Sec II, we give a simple introduction of the four-level atomic model allowing EIT and FWM, describe the metamaterial model, and show the similarity between the two models. The propagation of linear plasmonic polaritons in the metamaterial is discussed in detail. In Sec. III, we derive the coupled nonlinear envelope equations and present the vector plasmonic dromion solutions when the nonlinear varactors are mounted onto the gaps of the SRRs. Lastly, in Sec IV we give a discussion and a summary of our work. Details of some calculating results are given in five appendixes.
, which has always a large absorption peak at ω = 0 for arbitrary Ω c1 and Ω c2 .

II. EIT-BASED ATOMIC FWM AND ITS METAMATERIAL ANALOGUE
A. EIT-based FWM in a double-Λ-type four-level atomic system For a detailed comparison with the metamaterial model presented in the next subsection, we first give a brief introduction on a lifetime-broadened atomic gas with a double-Λ-type four-level configuration, shown in Fig. 1(a). In this system, two weak probe laser fields with central angular frequencies ω p1 and ω p2 and wavevectors k p1 and k p2 drive respectively the transitions |1 ↔ |3 and |1 ↔ |4 , and two strong control laser fields with central angular frequencies ω c1 and ω c2 and wavevectors k c1 and k c2 drive respectively the transitions |2 ↔ |3 and |2 ↔ |4 . The total electric fields in this system is given , where e jn and E jn (j = p, c; n = 1, 2) are respectively the unit vector denoting the polarization direction and the envelope of the corresponding laser field. Note that for simplicity all the laser fields are assumed to be injected in the same (i.e. z) direction (which is also useful to suppress Doppler effect). Under electric-dipole approximation and rotating-wave approximation (RWA), the Hamiltonian of the system in interaction picture readŝ )/ are respectively one-, two-, and three-photon detunings, with E l the eigenenergy of the atomic state |l (l = 1, 2, 3, 4); Ω p1 = (e p1 · p 31 )E p1 / , Ω p2 = (e p2 · p 41 )E p2 / , Ω c1 = (e c1 · p 32 )E c1 / , and Ω c2 = (e c2 · p 42 )E c2 / are respectively the half Rabi frequencies of the probe and the control laser fields, with p jl the electric dipole moment related to the transition |j ↔ |l . The Hamiltonian (1) allows three bright states and one dark state [50]. The dark state reads which is a superposition of only the two lower states |1 and |2 and has a zero eigenvalue.
The dynamics of the atoms is governed by the optical Bloch equation When the two weak probe fields are applied, the ground state |1 is not depleted much.
In this case, the Bloch equation reduces to with d j1 = ∆ j + iγ j1 with γ j1 = Γ 1j /2 (j = 2, 3, 4). Equations (3a)-(3c) describe the dynamics of three coupled harmonic oscillators [51], where σ 31 and σ 41 are bright oscillators due to their direct coupling to the probe fields Ω p1 and Ω p2 , but σ 21 is a dark oscillator because it has no direct coupling to any of the two probe fields.
The dynamics of the probe fields is governed by the Maxwell equation Here the polarization intensity is given by P = N 0 [σ 31 e i(k p1 z−ω p1 t) + σ 41 e i(k p2 z−ω p2 t) + c.c.], with N 0 the atomic density. Under a slowly-varying envelope approximation (SVEA), the Maxwell equation reduces to with κ 13 = N 0 |e p1 · p 13 | 2 ω p1 /(2 ε 0 c) and κ 14 = N 0 |e p2 · p 14 | 2 ω p2 /(2 ε 0 c). For simplicity, we assume the two control fields are strong enough and thus have no depletion during the evolution of the probe fields; additionally, the diffraction effect is negligible, which is valid for the probe fields having large transverse size.
We see that a transparency window is opened in the profile of Im(K + a ) near ω = 0; the transparency window becomes larger when the control fields are increased. The opening of the transparency window (called EIT transparency window) is due to the EIT effect contributed by the control fields. The blue solid cure in the figure is Im(K − a ) as a function of ω, which however has always a large absorption peak near ω = 0 irrespective of the value of the control fields. Below, for convenience we shall call the normal mode with the linear dispersion relation K + a (K − a ) as EIT-mode (non-EIT-mode). The double-Λ-type four-level system can be used to describe a resonant FWM process in atomic systems [1,[29][30][31][32]. The first laser field (i.e. the control field tuned to the |2 ↔ |3 transition with the half Rabi frequency Ω c1 ) and the second laser field (i.e. the probe field tuned to the |1 ↔ |3 transition with the half Rabi frequency Ω p1 ) can adiabatically establish a large atomic coherence of the Raman transition, described by the off-diagonal density matrix element σ 21 . The third laser field, i.e. the control field tuned to the |2 ↔ |4 transition with the half Rabi frequency Ω c2 , can mix with the coherence σ 21 to generate a fourth field with the half Rabi frequency Ω p2 resonant with the |1 ↔ |4 transition. For details, see Refs. [1,[29][30][31][32] and references therein. We assume that an incident gigahertz radiation E = e x E x +e y E y [with E j = E j0 e −iωpt +c.c. The dependence on the excitation condition for the absorption spectrum can be briefly explained as follows. A sole CW in the meta-atoms is function as an optical dipole antenna and thus serves as a bright (or radiative) oscillator, which can be directly excited by the incident radiation. The surface current in an excited SRR can be clockwise or anticlockwise direction, indicating that there is no direct electric dipole coupling with the incident radiation and hence the SRR serves as a dark or trapped oscillator with long dephasing time [55].  Based on the solution given by Eq. (B1), we deduce that, in the case of ω 3 = ω p , γ 3 = 0 and g 1 = g 2 , Eq. (6) allows a "dark state" (i.e. the state where both the bright oscillators are not excited, i.e. q 10 = q 20 = 0) exists, if This "dark state" condition is equivalent to the one obtained in the four-level double-Λ-type atomic system, given by Eq. (2). Obviously, the PIT-mode shown in Fig. 2(c) corresponds to the case κ 2 = −κ 1 , where the minus symbol can be understood as a π-phase difference resulting in a cooperative coupling effect, which is assumed in all the numerical calculations carried out below.
The equation of motion of the EM wave is governed by the Maxwell equation where 1(2) , with N m the unit-cell density, e the unit charge, and χ (1) D the optical susceptibility of the hosting material. We assume the distance between the meta-atoms is large so that the interaction between them can be neglected.
Assuming the central frequency of the incident radiation ω p is near the natural frequencies of the Lorentz oscillators described by Eq. (6) [55], a resonant interaction occurs between the incident radiation and these oscillators. To deal with the propagation problem of the plasmonic polaritons in the system analytically, we assume E j (r, t) = E j (z, t)e i(kpz−ωpt) + c.c.
are slowly-varying envelopes and ∆ l = ω p − ω l is a small detuning. With this ansatz and under RWA, Eq. (6) is simplified into the reduced Lorentz equation with d j = ∆ j + γ j /2. We see that the reduced Lorentz equation (9) describing the unit cell has the same form as the optical Bloch equation (3) describing the four-level double-Λ atom.
Consequently, each unit cell in the metamaterial is analogous to a four-level double-Λ-type atom in the atomic gas presented in the last subsection. That is to say, the unit cell is indeed a meta-atom, where the bright-oscillator excitation in the CW A (CW B) driven by E x (E y ) is equivalent to the dipole-allowed transition |1 ↔ |3 (|1 ↔ |4 ) driven by the probe field Ω p1 (Ω p2 ), and the dark-oscillator excitation in the SRR is equivalent to the dipole-forbidden transition |1 → |2 in the four-level double-Λ-type atom. We also see that the coupling between the CW A (CW B) and the SRR, described by κ 1 (κ 2 ), is equivalent to the control field Ω c1 (Ω c2 ) driven the atomic transition |2 ↔ |3 (|2 ↔ |4 ).
Under SVEA, the Maxwell equation in the metamaterial reads It is easy to get the linear dispersion relation where , with κ f j = κ j /(2ω p ) and g f j = g j /(2ω p ) (j = 1, 2). As expected, the metamaterial system allows two normal modes with the linear dispersion relation respectively given by K + m and K − m . In fact, K + m (K − m ) is a PIT-mode (non-PIT-mode) of the system, as explained below.
The character of the above two normal modes can be clearly illustrated by plotting K + When plotting the figure, the system parameters are taken from Appendix B, and additional parameters are chosen by κ 0 = 10 10 kg/(cm · s 2 · C) [56] and ∆ j = 0 (j = 1, 2, 3). We see that Im(K + m ) displays a transparency window (called PIT transparency window) near ω = 0, analogous to the EIT transparency window in Im(K + a ) of the four-level double-Λ-type atomic system [red dashed line and green dashed-dot line in Fig. 1(b)]. The steep slope of Re(K + m ) indicates a normal dispersion and a slow group velocity of the plasmonic polariton. As the coupling strength between the CWs and the SRR gets larger (i.e. the separations d x and d y is reduced), the PIT transparency window becomes wider and deeper, and the slope of Re(K + m ) gets flatter. The opening of the PIT transparency window is attributed to the destructive interference between the two bright oscillators and the dark oscillator through cooperative near-field coupling. Shown in Fig. 3(b) is the imaginary (red solid line) and the real (blue dashed line) of K − m , which is nearly independent on the coupling constant κ 1 (κ 2 = −κ 1 ). We see that Im(K − m ) has a single, large absorption peak and Re(K − m ) has an abnormal dispersion near ω = 0, analogous to Im(K − a ) of the double-Λ-type atomic system [blue solid line in Fig. 1(b)].

C. Propagation of linear plasmonic polaritons via an analogous FWM process of atomic system
As indicated above, the meta-atoms in the present metamaterial system are analogous to the four-level atoms with the double-Λ-type configuration, and hence an analogous resonant FWM phenomenon for the plasmonic polaritons is possible. That is to say, if initially only one polarization-component of the EM wave (e.g. x-component) is injected into the metamaterail, a new polarization-component (e.g. y-component) will be generated through two equivalent control fields (i.e. the couplings between the SRR and CWs, described by κ 1 and κ 2 ). To illustrate this, we present the solution of the ML equations (9) and (10) which can be obtained by using Fourier transform [32,57].
and F ± 0 is the initial amplitude of the normal mode K ± m determined by given excitation condition. We assume initially only the x-component of the EM field in input to the system, i.e. the initial condition for the EM field is given by E x (0, t) = 0, E y (0, t) = 0. By Eq. (12) we have For simplicity, we consider the adiabatic regime where the power series of K ± m and G ± on ω converge rapidly. By taking K ± m = K ± 0 + ω/V ± g + O(ω 2 ) and G ± = G ± 0 + O(ω), we readily obtain where τ ± = t − z/V ± g , with V ± g ≡ (∂K ± m /∂ω) −1 | ω=0 being the group-velocity of the normal mode K ± m . The conversion efficiency of the FWM is given by η where L is the medium length. For the case κ 2 = −κ 1 , one has Im(K − 0 ) ≫ Im(K + 0 ), which means that the K − m mode decays away rapidly during propagation and hence can be safely neglected. Then Eq. (14) is simplified as We see that the x-and y-polarization components of the EM wave have matched group velocity V + g . The expression of the FWM conversion efficiency reduces into Shown in Fig. 4 is the FWM conversion efficiency η as a function of the dimensionless optical depth (κ 0 g f 1 /γ 1 )L for ∆ 1 = ∆ 2 = 0 (blue dashed line) and for ∆ 1 = ∆ 2 = 5γ 1 (red solid line). When plotting this figure, we have set ∆ 3 = 0 and γ 3 ≈ 0 in order for a better analogue to the atomic system. The influence of γ 3 can be effectively reduced by introducing a gain element into the gaps of the SRRs (see the discussion in Sec. IV). From the figure, we see that for the case of exact resonance (i.e. ∆ 1 = ∆ 2 = 0), the FWM efficiency η increases and rapidly saturates to 25% when the dimensionless optical depth (κ 0 g f 1 /γ 1 )L ≈ 5 (i.e. L ≈ 0.9 cm), indicating a unidirectional energy transmission from E x to E y . For the case of far-off resonance (i.e. ∆ 1 = ∆ 2 = 5γ 1 ), the FWM efficiency displays a damped oscillation in the interval 0 < (κ 0 g f 1 /γ 1 )L < 250, indicating a back-and-forth energy exchange between E x and E y ; eventually the efficiency reach to the steady-state value 25% when (κ 0 g f 1 /γ 1 )L ≥ 300 (see the inset). Interestingly, the value of the FWM conversion efficiency may reach to η ≈ 76% at (κ 0 g f 1 /γ 1 )L ≈ 15 (i.e. L ≈ 3 cm).

III. VECTOR PLASMONIC DROMIONS IN THE PIT METAMATERIAL
Note that when deriving Eq. (10), the diffraction effect has been neglected, which is invalid for the plasmonic polaritons with small transverse size or long propagation distance; furthermore, because of the highly resonant (and hence dispersive) character inherent in the PIT metamaterial, the linear plasmonic polaritons obtained above inevitably undergo significant distortion during propagation. Hence it is necessary to seek the possibility to obtain a robust propagation of the plasmonic polaritons in the PIT metamaterial. One way to solve this problem is to make the PIT system work in a nonlinear propagation regime.
In recent years, nonlinear metamaterials have attracted much attention due to their potential applications (see Ref. [58] and references therein). One suitable way to design a nonlinear PIT metamaterial in microwave and lower THz ranges is to use nonlinear insertions onto the meta-atoms [59,60]. Here, as suggested in Refs. [24,59,60], we assume the nonlinear insertion in the PIT metamaterial are varactor diodes, which are mounted onto the gaps of the SRRs [60] [see Fig. 2(b)].

A. Nonlinear envelope equations
Since the introduction of the nonlinear element onto the SRRs, Eq. (6c) should be replaced by [60] where α and β are nonlinearity coefficients, described in Appendix E.
Due to the quadratic and cubic nonlinearities in Eq. (17), the input EM field (with only a fundamental wave) will generate longwave (rectification), and second harmonic compo- with θ p = (2k p + ∆k)z − 2ω p t and ∆k a detuning in wavenumber. The oscillations of the Lorentz oscillators in the meta-atoms have the form q j (r, t) = q dj (r, t) Substituting these expressions into Eqs. (6a), (6b), (8), and (17), and adopting RWA and SVEA, we obtain a series of equations for the motion of q µj and E µl , listed in Appendix C.
We solve the equations for q µj and E µl by using the standard method of multiple scales [61]. Take the asymptotic expansion q f j = ǫq (2) f l + · · · , and E dl = ǫ 2 E (2) dl + · · · (here ǫ is a dimensionless small parameter characterizing the amplitude of the incident EM field), and assume all quantities on the right sides of the asymptotic expansion as functions of the multiscale variables [61] z l = ǫ l z (l = 0, 1, 2) and t l = ǫ l t (l = 0, 1). Substituting the expansion into the equations for q µj and E µj and comparing powers of ǫ, we obtain a chain of linear but inhomogeneous equations which can be solved order by order.

The leading order [i.e. O(ǫ)] solution reads E
(1) f x = F + e iθ + and E (1) f y = G + F + e iθ + , where θ + = K + m z 0 − ωt 0 and F + is a slowly-varying envelope function to be determined in higherorder approximations. The expression of G + is given in Sec. II C. Here we consider only the PIT (i.e. K + m ) mode because the non-PIT (i.e. K − m ) mode decays rapidly during propagation, as indicated in the last section. The solution for q where K + 2 ≡ [∂ 2 K + m /∂ω 2 ]| ω=0 describes the group-velocity dispersion of the fundamental wave; R 0 is a coefficient characterizing the coupling between the fundamental and long waves, with the expression given in Appendix D; χ here V + p ≡ [n + m (0)/c] −1 is the phase-velocity of the longwave, defined by n + m (0) = n + m | ωp=0,ω=0 with n + m (ω; ω p ) = c[k p (ω p ) + K + m (ω; ω p )]/(ω p + ω). It is easy to obtain with X j = ω j ω 3 − κ 2 j . Under the PIT condition [i.e. (κ j /2ω p ) 2 ≫ γ j γ 3 /4, j = 1, 2], the real parts of the nonlinearity susceptibilities χ obtained by neglecting the diffraction term in Eq. (20) and then plugging the derived expression for Q + into Eq. (18). The second term in Eq. (22) is due to the longwave-shortwave interaction, where one can clearly see that when the group-velocity of the shortwave is nearly equal to the phase-velocity of the longwaves (i.e. V + g ≈ V + p ), the effective third-order nonlinear susceptibility χ (3) eff can be further enhanced, which can be realized by adjusting the coupling parameters κ 1 and κ 2 . Compared with the conventional PIT-based metamaterials [23,24] with one bright and one dark oscillators, the present FWM-based metamaterial consisting double bright and one dark oscillators allows for a stronger bright-dark coupling, resulting in an enhanced Kerr nonlinearity. As a result, the coupling strength for the longwave-shortwave interaction in the present FWM-based metamaterial is √ 2 time larger than that in the conventional PIT-based metamaterials.

B. Vector plasmonic dromions
Equations (18) and (20) show that the self-interaction of the shortwave F + can stimulate the generation of the longwave Q + [Eq. (20)], and at the same time the longwave Q + has a back-action to the shortwave F + [Eq. (18)]. Such equations admit solutions describing the excitation of high-dimensional, two-component (vector) nonlinear plasmonic polaritons in the PIT metamaterial. To demonstrate this, we convert them into the dimensionless form Here R x (R y ) is the typical radius of the incident EM field in the x (y) direction; τ 0 is the typical pulse duration of the probe field; U 0 (V 0 ) is the typical amplitude of the longwave (shortwave) envelope; L disp = −τ 2 0 /Re(K + 2 ), L diff ≡ n D ω p R 2 x /c and L nln = (2n D c)/[ω p U 2 0 Re(χ + )] are, respectively, the typical dispersion length, the typical diffraction length, and the typical nonlinearity length. Note that in obtaining Eq. (23), we have neglected the small imaginary parts of χ (3) + and K + 2 , which is reasonable under the PIT condition as discussed above.
In favor of the formation of plasmonic dromions, we take the following two assumptions.
First, we shall assume R x ≪ R y and thus g d0 ≪ 1 so that the original (3+1)-dimensional nonlinear problem can be reduced into a (2+1)-dimensional one. Second, we assume the contribution of the dispersion, diffraction and the nonlinearity effects are of the same level, which can be achieved by taking L diff = L disp = L nln and thus we obtain τ 0 = R x −ω p n D Re(K + 2 )/c and U 0 = [c/(ω p R x )] 2/Re[χ with v 1 = v + 2|u| 2 , where, for convenience, we have performed a 45-degree rotation of The DS-I equation (24) can be exactly solved via the Hirota's bilinear method [63] and various dromion solutions can be obtained.

Shown in
In doing so, we first integrate Eq. (24b), yielding where the nonzero boundary conditions V 11 | τ 1 →−∞ and V 12 | ξ 1 →−∞ are given by [64] In the numerical simulation, the integral Eq. Within the forth-order approximation, the explicit expression for the EM field in the metamaterial takes the form When u and v are taken as the dromion solution given above, we obtain a vector plasmonic dromion since the EM field (28) has two polarization components, with each component a plasmonic dromion. Note that, different from the result in the scalar model considered before [24], the polarization of the EM field obtained here can be actively selected by adjusting the separation between the CWs and SRR [i.e. d x and d y in Fig. 2(b) and hence the coupling constants κ 1 and κ 2 ], which can be served as a polarization selector for practical applications [65,66].
The threshold of the power density of the vector plasmonic dromion given above can be estimated by using Poynting vector. Based on the above system parameters, the average power of the vector plasmonic dromion is estimated as P = 6.1 mW.
We see that due to the resonant character of the PIT effect in the system, extremely low generation power is required for generating the vector plasmonic dromion.

IV. DISCUSSION AND SUMMARY
It should be mentioned that in writing the dark-state condition (7), the damping coefficient γ 3 of the dark oscillator in the meta-atoms is assumed to be small. However, due to the Ohmic loss inherent in the metal that construct the metamaterial, by our numerical calculation the numerical value of γ 3 is about 0.18 GHz, which, though smaller than damping coefficients of the CWs (γ 1 = γ 2 = 2.1 GHz), is still large and has inevitably detrimental impact on the PIT quality. In order to improve the performance of the PIT, one can suppress γ 3 by introducing a gain element into the SRR of the meta-atoms. One possible way is the use of tunneling diodes that have negative resistance and hence may provide gain to the PIT-based metamaterial [67,68]. Such method has been recognized to be useful for suppressing and even cancelling γ 3 , particularly in microwave and THz regimes.
In conclusion, in this article we have considered a plasmonic metamaterial interacting with an EM field with two polarization components. We have proved that such metamaterial can be taken as a classical analogue of an atomic gas with a double-Λ-type four-level configuration coupled with four laser fields, displays an PIT effect and an equivalent process of atomic FWM. We have shown that, when the nonlinear varactors are mounted onto the gaps of the SRRs, the metamaterial system can acquire giant second-and third-order Kerr nonlinearities via the PIT and the longwave-shortwave interaction. We have also shown that the system supports high-dimensional vector plasmonic dromions, which have very low generation power and are robust during propagation. Our work not only contributes a plasmonic analogue of atomic EIT and FWM but also provides a way for generating novel plasmonic polaritons, and hence opens a new avenue on the exploration of PIT effect in metamaterials.