Using a combination of numerical and analytical techniques we demonstrate that a metal stripe surrounded by the active and passive dielectrics supports propagation of stable spatial surface-plasmon solitons. Our analytical methods include the multiple scale reduction of the Maxwell’s equations to the coupled Ginzburg-Landau system, and the soliton perturbation theory developed in the framework of the latter.
© 2011 Optical Society of America
Surface plasmon polaritons (SPPs) provide one of the favored approaches to realization of on-chip photonic devices and are a well established tool in sensing applications. While SPPs are exponentially localized in the direction perpendicular to the metal-dielectric interface by the natural boundary conditions, one should take a special care about suppression of their in-plane diffraction . An interesting alternative to various geometrical methods providing lateral confinement of SPPs is to use the concept of spatial solitons, where diffraction is suppressed by the nonlinearity induced focusing, see, e.g., [2, 3]. Characteristic loss length of SPPs is typically order of few tens of microns. In order to talk about diffraction compensation over such distances one should work with the beams having diffraction length shorter than the loss length. Such strong diffraction requires very large powers to achieve the soliton effect. One can structure the metal surface and thus weaken the diffraction and hence power requirements, and work with the discrete SPP solitons [4, 5].
On the other hand, the spatial soliton concept can be extended by complementing the diffraction vs nonlinearity balance with the gain vs loss balance. Thus completely solving the problem of the soliton decay due to linear absorption. The paradigm model in this context is the complex Ginzburg-Landau (GL) equation [6, 7]. Cubic GL equation has been derived for the SPPs at the boundary with an active dielectric . However, SPP solitons reported in the above work demonstrate substantial instabilities developing, both, in the soliton core and background . The nature of these instabilities is the same as the one known in other physical contexts with the supercritical transition between the zero and non-zero spatially homogeneous solutions . Making this transition subcritical is known to favor stabilization of solitons. This can be achieved, e.g., if the quintic nonlinearity is accounted for [6, 8, 9]. However, the higher order nonlinearities usually become important only for sufficiently high intensities. An alternative method to suppress the soliton instabilities is to drain the energy generated by the gain in the cubic GL equation using the coupling of the latter to the second purely lossy GL equation [10–12].
Below we consider a thin metal stripe, where SPPs at the two interfaces are weakly coupled. Dielectric at one of the interfaces is assumed active and at the other is passive. Using the method of multiple scales applied to the Maxwell equations we derive the system of two coupled GL equations for the SPP amplitudes and report stable spatial SPP solitons. Note, that the amplification of SPPs without considering transverse effects has been previously studied theoretically [13, 14] and experimentally using rhodamine 6G, quantum dots and erbium [15–17].
2. Reduction of the Maxwell’s equations to the coupled Ginzburg-Landau system
We start our analysis from the time independent Maxwell equations ∇⃗ × ∇⃗ × E⃗ = ɛE⃗. The spatial coordinates are normalized to k = λ/(2π), where λ is the vacuum wavelength. x is the coordinate perpendicular to the metal-dielectric interfaces, z is the propagation coordinate, while the in-plane diffraction of SPPs happens along y direction, see Fig. 1(a). The metal stripe has the width w along the x-axis and is infinite along y and z. ɛmetal = ɛ″m + iɛ″m is the dielectric constant inside the metal (ɛ′m < 0, ɛ″m > 0). Inside the active layer at the top of the metal stripe we have . Here χ3 = χ′3 + iχ″3, χ″3 > 0 and ɛ″a < 0 corresponding to the nonlinear absorption and linear gain, respectively. Our model of the active dielectric corresponds to the two-level atom medium in the regime far from saturation . ɛ in the passive dielectric layer is . We have expressed the real parts of the linear dielectric constants of the top layer and of the substrate as , because it is more straightforward to assume that in the zero approximation of our perturbation theory the SPPs at both interfaces are identical, so that the α terms appear only in the higher orders. This assumption helps to avoid unnecessary complexity of the analytical calculations.
Assuming that ɛ″m, α, ɛ″a and |χ3||E⃗|2 are of the order s (where s is the small dummy parameter) and that the width w of the metal stripe is much greater than the skin depth , we have the following expressions for the SPPs at the two interfaces:18]. Proceeding to the order s3/2 one obtains a system of equations for the correction δE⃗ = (δEx,δEz)T (|δE⃗| ∼ O(s3/2)): L̂δE⃗ =(N̂ – L̂)(ψpe⃗p + ψae⃗a) – b⃗, where and
The propagation equations for the amplitudes ψp and ψa can be derived by projecting the left and right hand side of the equation for δE⃗ onto e⃗p,a and assuming that the interaction of the two SPPs through the overlap of their tails becomes important only in the s3/2 order. Note, that the calculation of the projection integrals and use of the boundary conditions for all the components of the electric and magnetic fields reveal that the left hand side yields a nontrivial contribution due to discontinuities of Ex and finite values of Ez at the interfaces . The terms containing ϕp,a partly cancel out and partly can be expressed via ψp,a. The resulting system is the two coupled complex GL equations:
3. Subcritical transition to the spatially homogeneous SPPs
If coupling between the two interfaces is neglected, κ = 0, then the trivial solution ψa = 0 loses its linear stability at g = l, in favor of |ψa|2 = (g – l)/(fϒ″). The latter obviously exists only for g > l, so that the transition between the trivial and nontrivial states is supercritical. Thus any SPP soliton existing on the zero background for g > l suffers from the instability of its tails .
The transition between the zero and non-zero SPPs is however more involved for κ ≠ 0. The trivial solution ψa = ψp = 0 becomes unstable for g > gth, where the latter solves the cubic equationFig. 1(b) shows the range of g and Δβ values, where the non-zero solution coexist with the zero one. One can see, that the range of the suitable gain is maximal for Δβ = 0. An example of the subcritical dependence of |R|2 vs g is shown in Fig. 1(c). From this analysis we conjecture that in the subcritical case the soliton solutions on the zero background, if they exist, will have stable tails. In our numerical examples we assume ɛd = 1.8, the metal is silver, and λ = 530nm, so that ɛ′m = −15, ɛ″m = 0.4. Under these condition the dimensionless gain values used in Figs. 1 and 2 are scaled back to the physical units by multiplying them with ∼ 12μm−1. Thus to see the subcritical transition the gain should reach values boosting the intensity by the factor e after ∼ 20μm of the propagation distance. For comparison, the gain length demonstrated with dyes  was close to this value, while in the quantum dot case  it was one order of magnitude larger. The intensity plot shown in Fig. 1(c) corresponds to the nonlinear response of the material and pump powers such that the nonlinearity induced correction to the refractive index (given by fϒ|R|2) has a realistic value ∼ 10−3.
4. Spatial solitons
To analyse the existence of the soliton solutions we develop an approximate analytical theory and use a numerical approach based on the Newton-Raphson iterations. Differentiating the power integral and using Eqs. (2) we find the following balance equationEq. (4) to derive an evolution equation for this parameter and determine its stationary values. For the single cubic-quintic GL, the soliton of the ideal nonlinear Schrodinger equation has been successfully used in this approach [8, 9]. Neglecting all the dissipative terms l = g = ϒ″ = 0 in Eqs. (2), and assuming Δβ = 0, still leaves us with the system, where the solitons can not be found in a closed analytical form. Therefore we need to make further simplifying assumption that κ is a small parameter and |ψa| ≫ |ψp|. In this case Eq. (4), we find the equation for η:
Stationary (dη/dz = 0) values of η obey the cubic equation for η2, which roots predict the soliton families. Within the range of the subcritical behavior of the homogeneous solution we have found that there can be up to two positive solutions for η2 corresponding to the two different solitons. Maximum of |ψa| corresponding to these two branches is plotted as a function of the gain parameter g, in the Fig. 2(a) with the full lines. The soliton profile corresponding to the large amplitude branch is shown in Fig. 2(b) with the full line.
Although the soliton profiles are predicted by the perturbation theory quite accurately, they are not exact. Therefore, the use of numerical methods is necessary. We seek solitons in the form ψa,p = fa,p(y)eiμz, approximate with the 2nd order finite differences, and apply the Newton-Raphson iterations to the resulting set of the algebraic equations. Parameter μ is treated as one of the unknowns, while the perturbation theory provides an initial guess for the iterative procedure. Note that, since μ is an unknown, the system formally contains N equations and N + 1 unknowns. However, we reduce the number of unknowns by one, fixing Imψa to zero at y = 0, which always can be achieved by an appropriate selection of the arbitrary common phase of ψp and ψa. Numerically found branches of the solitons and an example of the soliton profile are shown in Fig. 2(a) and 2(b), respectively, alongside with the ones found using our analytical approach. Scaling back into physical units one finds that the soliton width is order of 1μm.
Very good quantitative agreement between our numerical and analytical approaches is evident. Note, however, that the perturbation theory predicts that the branch of the small amplitude solitons exists not only in the subcritical range (g < gth), but also for g > gth. While the former prediction is well confirmed by the numerical results, the latter one is not. The large amplitude solitons exist on both sides from gth, while their tails remain stable only for g < gth. We have also solved Eqs. (2) directly using a split-step method and have found that the large amplitude branch of the solitons remains stable over and well beyond practical propagation distances providing g < gth, while the small amplitude ones are unstable. The instability of the latter is not unexpected, since the small amplitude solitons separate the attraction basins of the stable zero solution and of the stable large-amplitude soliton.
In this work we demonstrate that the ubiquitous dielectric-metal-dielectric geometry can support propagation of stable spatial plasmon solitons. The stability is achieved by means of the coupling of the actively pumped SPP at one interface to the passive (lossy) one at the other. The lossy SPP absorbs the perturbations that would otherwise destabilize the zero background of the soliton. We have used the multiple scale approach to demonstrate the reduction of the Maxwell’s equations to the analytically and numerically tractable system of the coupled Ginzburg-Landau equations. We have found the gain interval of the coexistence of the stable zero amplitude solution with the stable spatially homogeneous SPP. Applying perturbative and numerical methods we have proved that within this interval there exist stable spatial plasmon solitons, such that the in-plane diffraction and absorption are both balanced through the interplay of the linear gain and complex nonlinearity, which was assumed coming from the dielectric. From our analysis it follows that, unlike the dielectric-metal-dielectric case, the metal-dielectric-metal geometry is not favorable for the existence of this type of solitons, while future studies in the multilayered structures can lead to interesting results.
We acknowledge useful comments from A. V. Gorbach and W. J. Firth.
References and links
1. D. K. Gramotnev and S. I. Bozhevolnyi, “Plasmonics beyond the diffraction limit,” Nat. Photonics 4, 83–91 (2010). [CrossRef]
6. N. N. Ahmediev and A. Ankiewicz, Solitons: Nonlinear Pulses and Beams (Chapman and Hall, 2007).
7. A. Marini and D. V. Skryabin, “Ginzburg-landau equation bound to the metal-dielectric interface and transverse nonlinear optics with amplified plasmon polaritons,” Phys. Rev. A 81, 033850 (2010). [CrossRef]
8. B. A. Malomed, “Evolution of nonsoliton and quasiclassical wavetrains in nonlinear Schrdinger and Korteweg -de Vries equations with dissipative perturbations,” Physica D 29, 155–172 (1987). [CrossRef]
10. B. A. Malomed and H. G Winful. “Stable solitons in two-component active systems,” Phys. Rev. E 53, 5365 (1996). [CrossRef]
11. J. Atai and B. A. Malomed, “Stability and interactions of solitons in two-component active systems,” Phys. Rev. E 54, 4371 (1996). [CrossRef]
12. W. J. Firth and P. V. Paulau, “Soliton lasers stabilized by coupling to a resonant linear system,” Eur. Phys. J. D 59, 13–21 (2010). [CrossRef]
13. D. J. Bergman and M. I. Stockman“Surface plasmon amplification by stimulated emission of radiation: Quantum generation of coherent surface plasmons in nanosystems,” Phys. Rev. Lett. 90, 027402 (2003). [CrossRef] [PubMed]
15. M. A. Noginov, V. A. Podolskiy, G. Zhu, M. Mayy, M. Bahoura, J. A. Adegoke, B. A. Ritzo, and K. Reynolds, “Compensation of loss in propagating surface plasmon polariton by gain in adjacent dielectric medium,” Opt. Express 16, 1385 (2008). [CrossRef] [PubMed]
16. P. M. Bolger, W. Dickson, A. V. Krasavin, L. Liebscher, S. G. Hickey, D. V. Skryabin, and A. V. Zayats, “Amplified spontaneous emission of surface plasmon polaritons and limitations on the increase of their propagation length,” Opt. Lett. 35, 1197–1199 (2010). [CrossRef] [PubMed]
18. D. V. Skryabin, A. Gorbach, and A. Marini, “Surface induced nonlinearity enhancement of TM-modes in planar subwavelength waveguides,” J. Opt. Soc. Am. B 28, 109–114 (2011). [CrossRef]