###### Abstract

We analyze nonlinear effects in optically driven arrays of nonlinear metallic nanoparticles. We demonstrate that such plasmonic systems are characterized by a bistable response, and they can support the propagation of dissipative switching waves (or plasmonic kinks) connecting the states with different polarization. We study numerically the properties of such plasmonic kinks which are characterized by a subwavelength extent and a tunable velocity.

Subwavelength plasmonic kinks in arrays of metallic nanoparticles

Roman E. Noskov, Pavel A. Belov, and Yuri S. Kivshar

nanometa@gmail.com

OCIS codes: (190.1450) Bistability; (190.6135) Spatial solitons; (250.5403) Plasmonics; (250.6715) Switching.

## References and links

- [1] M.I. Brongersma and P.G. Krik, Eds. Surface Plasmon Nanophotonics (Spinger, Dordrecht, 2007), 268 pp.
- [2] S. Maier, Plasmonics: Fundamentals and Applications (Springer, New York, 2007), 219 pp.
- [3] D.K. Gramotnev and S.I. Bozhevolnyi, “Plasmonics beyond the diffraction limit,” Nat. Photonics 4, 83-91 (2010).
- [4] J. Takahara, S. Yamagishi, H. Taki, A. Moromoto, and T. Kobayashi, “Guiding of a one-dimensional optical beam with nanometer diameter,” Opt. Lett. 22, 475-477 (1997).
- [5] K. Li, M. Stockman, and D. Bergman, “Self-similar chain of metal nanospheres as an efficient nanolens,” Phys. Rev. Lett. 91, 227402 (2003).
- [6] N.C. Panoiu and R.M. Osgood, “Subwavelength nonlinear plasmonic nanowire,” Nano Lett. 4, 2427-2430 (2004).
- [7] W. Fan, S. Zhang, N.C. Panoiu, A. Abdenour, S. Krishna, R.M. Osgood, K.J. Malloy, and S.R.J. Brueck, “Second harmonic generation from a nanopatterned isotropic nonlinear material,” Nano Lett. 6, 1027-1030 (2006).
- [8] J.A.H. van Nieuwstadt, M. Sandke, S. Enoch, and L. Kuipers, “Strong modification of the nonlinear optical response of metallic subwavelength hole arrays,” Phys. Rev. Lett. 97, 146102 (2006).
- [9] G.A. Wurtz, R. Pollard, and A.V. Zayats, “Optical bistability in nonlinear surface-plasmon polaritonic crystals,” Phys. Rev. Lett. 97, 057402 (2006).
- [10] A.A. Zharov, R.E. Noskov, and M.V. Tsarev, “Plasmon-induced terahertz radiation generation due to symmetry breaking in a nonlinear metallic nanodimer,” J. Appl. Phys. 106, 073104 (2009).
- [11] R.E. Noskov, A.A. Zharov, and M.V. Tsarev, “Generation of widely tunable continuous-wave terahertz radiation using a two-dimensional lattice of nonlinear metallic nanodimers,” Phys. Rev. B 82, 073404 (2010).
- [12] Y. Liu, G. Bartal, D.A. Genov, and X. Zhang, “Subwavelength discrete solitons in nonlinear metamaterials,” Phys. Rev. Lett. 99, 153901 (2007).
- [13] A. Marini, D.V. Skryabin, B. Malomed, “Stable spatial plasmon solitons in a dielectric-metal-dielectric geometry with gain and loss,” Opt. Express 19, 6616 (2011).
- [14] A. Marini, A.V. Gorbach, and D.V. Skryabin, “Coupled-mode approach to surface plasmon polaritons in nonlinear periodic structures,” Opt. Lett. 35, 3532 (2010).
- [15] F. Ye, D. Mihalache, B. Hu, and N.C. Panoiu, “Subwavelength plasmonic lattice solitons in arrays of metallic nanowires,” Phys. Rev. Lett. 104, 106802 (2010).
- [16] K.M. Leung, “Optical bistability in the scattering and absorption of light from nonlinear microparticles,” Phys. Rev. A 33, 2461 (1986).
- [17] S. Yong and D. Stroud, “Surface-plasmon dispersion relations in chains of metallic nanoparticles: An exact quasistatic calculation,” Phys. Rev. B 69, 125418 (2004).
- [18] P.B. Johnson and R.W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 9, 4370 (1972).
- [19] E.D. Palik, Handbook of Optical Constants of Solids (Academic, Orlando, 1985).
- [20] V.P. Drachev, A.K. Buin, H. Nakotte, and V.M. Shalaev, “Size dependent for conduction electrons in Ag nanoparticles,” Nano Lett. 4, 1535 (2004).
- [21] M.J. Weber, Handbook of Optical Materials (CRC Press, 2003).
- [22] I.V. Shadrivov, A.A. Zharov, N.A. Zharova, and Y.S. Kivshar, “Nonlinear magnetoinductive waves and domain walls in composite metamaterials,” Photonics Nanostruct. Fundam. Appl. 4, 69 (2006).
- [23] N.N. Rosanov, N.V. Vysotina, A.N. Shatsev, I.V. Shadrivov, and Y.S. Kivshar, “Hysteresis of switching waves and dissipative solitons in nonlinear magnetic metamaterials,” JETP Lett. 93, 743 (2011).
- [24] O.M. Braun and Yu.S. Kivshar, The Frenkel-Kontorova Model: Concepts, Methods, and Applications (Springer-Heidelberg, 2004), 498 pp.

## Appendix A Introduction

When the size of conventional optical devices is reduced to nanoscale, the spatial confinement of light becomes inherently limited by diffraction. However, metal-dielectric structures allow to achieve a subwavelength light confinement with surface plasmon polaritons (SPPs), or plasmons [1, 2, 3]. In particular, by using plasmons excited in an array of resonantly coupled metal nanoparticles [4, 5], one can spatially confine and manipulate optical energy over distances much smaller than the wavelength.

Because of the strong enhancement of the field induced by the excitation of SPPs and increased optical nonlinearity, surface plasmons can be employed for the realization of a variety of nonlinear optical effects. In particular, several nonlinear optical processes have been demonstrated in plasmonic nanostructures, e.g., optical limiting and self-phase modulation in arrays of structured nanoparticles [6], second-harmonic generation and bistability in nanostructured metal films [7, 8, 9] and terahertz radiation generation in nonlinear metallic nanodimers [10, 11]. In addition, strong geometric confinement can boost efficiency of nonlinear optical effects, including the existence of subwavelength solitons in metal-dielectric multilayers [12, 13, 14] and arrays of metal nanowires [15]. The solitons supported by such media result from a balance between tunneling of surface plasmon modes and nonlinear self-trapping.

Here we study nonlinear effects in optically-driven arrays of nonlinear metal nanoparticles. In addition to the subwavelength plasmonic solitons analyzed earlier in other plasmonic systems [12, 13, 14, 15], we demonstrate that an array of nonlinear metal nanoparticles can support the propagation of dissipative switching waves (or plasmonic kinks) connecting the states with different polarization. These kinks become possible due to a bistable nonlinear response of each nanoparticle [16], and they have a subwavelength extent and a finite (or zero) velocity.

The paper is organized as follows. In Sec. 2, we describe our nonlinear model and discuss its realistic parameters. Section 3 summarizes our results for the study of the basic properties of plasmonic kinks and their generation. Finally, Sec. 4 concludes the paper.

## Appendix B Model and basic equations

We consider a linear array of identical spherical metallic (silver) nanoparticles embedded into a SiO host medium with permittivity , and is exited by the optical field (see Fig. 1). We assume that the particle radii and the center-to-center distance are nm and nm, respectively. Ratio satisfies the condition , so that we can employ the point dipole approximation [17]. In the optical spectral range, a linear part of silver dielectric constant can be written in a generalized Drude form , where , eV, eV [18] (hereinafter we accept time dependence); whereas dispersion of SiO can be neglected since for wavelengths 350 - 450 nm [19]. Nonlinear dielectric constant of silver is , where is the local field inside -th particle. We keep only cubic susceptibility due to spherical symmetry of particles. According to the model suggested in Ref. [20] and confirmed in experiment, 10 nm radii Ag spheres possess a remarkably high and purely real cubic susceptibility esu, in comparing to which the cubic nonlinearity of SiO is negligibly weak ( esu [21]).

To study nonlinear dynamics of the particle polarizations, , we employ the model based on the dispersion relation method. This technique was firstly applied to investigation of light scattering by nonlinear metal nanodimers [10, 11] and nonlinear effects in magnetic metamaterials [22, 23], and here we extend it to the case of an array of metal nanoparticle. According to this model, temporal dynamical response of metal nanoparticle arrays excited by the optical field with the frequency lying close to the frequency of the surface plasmon resonance of an individual particle, , can be described by the following system of coupled equations for slowly varying amplitudes of the particle dipole moments:

(1) |

where

and

are in charge of dipole-dipole interactions between particles, and are dimensionless slowly varying amplitudes of the particle dipole moments and external electric field, respectively, the indices ’’ and ’’ stand for the transverse and longitudinal components with respect to the array direction, , , describes both thermal and radiation losses of particles, , and . Equations (1) describe temporal nonlinear dynamics of an array of metal nanoparticles driven by arbitrary external optical field with the frequency . We stress that model takes into account all particle interactions through the dipole fields, and it can be applied both to finite and infinite arrays.

## Appendix C Bistable response and plasmonic kinks

In general, there are two fundamental nonlinear phenomena inherent in our system: modulational instability and bistability. Here we discuss the letter phenomenon. When the spatially distributed system demonstrates bistable response, it can support simultaneously two different stable states forming switching waves or kinks. However, to observe kinks and study its properties, one should satisfy one more condition - the stability criterium, i.e. provide a range of the external field intensities (lying inside the bistability region) which is free from modulational instability. Our findings show that for an array of nanoparticles this condition is fulfilled for the transversal excitation only. Therefore, we, next, consider just this case.

To begin with, we excite the infinite array by an homogenous electric field with the transversal polarization: . In this case, all particle dipole moments remain the same, , and the stationary solution of Eq. (1) can be written as follows

(2) |

where

Transition from to has been made via the replacement and taking into account symmetry structure of the series. When

the polarization becomes a three-valued function of leading to bistability (see Fig. 2(a)).

To analyze switching waves, we perform numerical simulations of Eq. (1) for a finite array (with 200 nanoparticles) at initial conditions corresponding to different branches of the homogeneous stationary solution (2) over different parts of the array. In Fig. 2(b) the profile of the stationary kink is presented. However, our calculations show that kinks can move as well. We have investigated impact of the field intensity on temporal kink dynamics. Characteristic results are summarized in Figs. 3(a,c). To simplify analysis, we fixed and set as a step-like function of time, as shown in Figs. 3(b,d). Constant results in the constant kink velocity, . That is why a step-like time dependency of allows us to control overall bistability zone.

Figure 3(a) demonstrates dynamical behavior of the system when increases from the lower to the upper threshold of the bistability region, as shown in Fig. 3(b). In this case the kink state changes from standing to moving, and for lying close to the middle of the bistability region the second kink is appeared at the left edge of the array, moving towards the first one. This is the straight consequence of cooperation between edge effects and modulational instability, which, in the end, crushes both kinks.

To avoid appearance of modulational instability, we also carry out the simulation when decreases from the middle to the lower threshold of the bistability region. The result is presented in Fig. 3(c). Here the kink velocity changes from the finite value to zero, and modulational instability growth is not observed.

Figure 4 shows the kink velocity as a function of the external field intensity . We notice that there exists a finite range of the field intensities where the kink does not move, i.e. . Generally speaking, this is a common property of nonlinear discrete systems [24] where the effective periodic potential created by the lattice requires a finite value of the applied external force to start kink’s motion. The existence of the kinks with zero velocity suggest a straightforward way to construct spatially localized modes by taking two standing kink of the opposite polarity. The resulting localized modes form the continuous families, and they resemble the similar localized modes discussed earlier for metal-dielectric multilayers [12, 13, 14] and arrays of metal nanowires [15]. On the other hand, we may expect generating such subwavelength plasmon solitons by considering interaction of two slowly moving kinks of the opposite polarity.

Here, we demonstrate this latter approach. We study numerically the interaction of two kinks of the opposite polarities, and select the value of the external field to realize the kinks with the near-zero velocity [see Fig. 5(a)]. We find that the kinks start moving slowly towards each other, then stop and form a stable dissipative plasmon soliton, as shown in Fig. 5(b). The generated plasmonic mode is localized only on 6 particles, and its width is about 150 nm. Because corresponds to the radiation wavelength nm, this plasmon soliton is truly subwavelength because it extends for .

It remains open if the kinks motion may demonstrate bistability, as was recently revealed for a similar nonlinear system with local coupling [23]. This will be a subject for future studies.

## Appendix D Conclusions

We have demonstrated that an array of interacting nonlinear metal nanoparticles can support a novel type of nonlinear localized modes, plasmonic kinks, which describe switching waves connecting two different states of polarization of the metal particles. We have generated various types of such kinks, analyzed numerically their velocity as a function of the external field intensity, and revealed that the array of metal nanoparticles can support deeply subwavelength dissipative plasmon solitons which in our case can be regarded as a bound state of two kinks of the opposite polarity.

## Acknowledgement

The authors acknowledge a support from Australian Research Council and Russian Ministry of Education and Science, as well as fruitful discussions with A.A. Zharov and N.N. Rosanov.