Rabi oscillations and stimulated mode conversion on the subwavelength scale

We study stimulated mode conversion and dynamics of Rabi-like oscillations of weights of guided modes in deeply subwavelength guiding structures, whose dielectric permittivity changes periodically in the direction of light propagation. We show that despite strong localization of the fields of eigenmodes on the scales below the wavelength of light, even weak longitudinal modulation couples modes of selected parity and causes periodic energy exchange between them, thereby opening the way for controllable transformation of the internal structure of subwavelength beams. The effect is reminiscent of Rabi oscillations in multilevel quantum systems subjected to the action of periodic external fields. By using rigorous numerical solution of the full set of the Maxwell's equations, we show that the effect takes place not only in purely dielectric, but also in metallic-dielectric structures, despite the energy dissipation inherent to the plasmonic waveguides. The stimulated conversion of subwavelength light modes is possible in both linear and nonlinear regimes. ©2015 Optical Society of America OCIS codes: (050.6624) Subwavelength structures; (260.2710) Inhomogeneous optical media References and links 1. I. I. Rabi, “On the process of space quantization,” Phys. Rev. 49(4), 324–328 (1936). 2. X. G. Zhao, G. A. Georgakis, and Q. Niu, “Rabi oscillations between Bloch bands,” Phys. Rev. B Condens. Matter 54(8), R5235–R5238 (1996). 3. M. C. Fischer, K. W. Madison, Q. Niu, and M. G. Raizen, “Observation of Rabi oscillations between Bloch bands in an optical potential,” Phys. Rev. A 58(4), R2648–R2651 (1998). 4. J. N. Winn, S. Fan, J. D. Joannopoulos, and E. P. Ippen, “Interband transitions in photonic crystals,” Phys. Rev. B 59(3), 1551–1554 (1999). 5. R. W. Robinett, “Quantum wave packet revivals,” Phys. Rep. 392(1-2), 1–119 (2004). 6. D. Marcuse, Theory of Dielectric Optical Waveguides (Academic, 1991). 7. K. O. Hill, B. Malo, K. A. Vineberg, F. Bilodeau, D. C. Johnson, and I. Skinner, “Efficient mode conversion in telecommunication fiber using externally written gratings,” Electron. Lett. 26(16), 1270–1272 (1990). 8. K. S. Lee and T. Erdogan, “Fiber mode coupling in transmissive and reflective tilted fiber gratings,” Appl. Opt. 39(9), 1394–1404 (2000). 9. Y. V. Kartashov, V. A. Vysloukh, and L. Torner, “Resonant mode oscillations in modulated waveguiding structures,” Phys. Rev. Lett. 99(23), 233903 (2007). 10. K. G. Makris, D. N. Christodoulides, O. Peleg, M. Segev, and D. Kip, “Optical transitions and Rabi oscillations in waveguide arrays,” Opt. Express 16(14), 10309–10314 (2008). 11. K. Shandarova, C. E. Rüter, D. Kip, K. G. Makris, D. N. Christodoulides, O. Peleg, and M. Segev, “Experimental observation of Rabi oscillations in photonic lattices,” Phys. Rev. Lett. 102(12), 123905 (2009). 12. C. D. Poole, C. D. Townsend, and K. T. Nelson, “Helical-grating two-mode fiber spatial-mode coupler,” J. Lightwave Technol. 9(5), 598–604 (1991). 13. D. McGloin, N. B. Simpson, and M. J. Padgett, “Transfer of orbital angular momentum from a stressed fiberoptic waveguide to a light beam,” Appl. Opt. 37(3), 469–472 (1998). 14. C. N. Alexeyev and M. A. Yavorsky, “Generation and conversion of optical vortices in long-period helical core optical fibers,” Phys. Rev. A 78(4), 043828 (2008). 15. C. N. Alexeyev, T. A. Fadeyeva, B. P. Lapin, and M. A. Yavorsky, “Generation of optical vortices in layered helical waveguides,” Phys. Rev. A 83(6), 063820 (2011). 16. Y. V. Kartashov, V. A. Vysloukh, and L. Torner, “Dynamics of topological light states in spiraling structures,” Opt. Lett. 38(17), 3414–3417 (2013). #231816 $15.00 USD Received 2 Jan 2015; revised 18 Feb 2015; accepted 18 Feb 2015; published 3 Mar 2015 © 2015 OSA 9 Mar 2015 | Vol. 23, No. 5 | DOI:10.1364/OE.23.006731 | OPTICS EXPRESS 6731 17. G. K. L. Wong, M. S. Kang, H. W. Lee, F. Biancalana, C. Conti, T. Weiss, and P. St. J. Russell, “Excitation of orbital angular momentum resonances in helically twisted photonic crystal fiber,” Science 337(6093), 446–449 (2012). 18. V. S. Shchesnovich and S. Chávez-Cerda, “Bragg-resonance-induced Rabi oscillations in photonic lattices,” Opt. Lett. 32(13), 1920–1922 (2007). 19. B. Terhalle, A. S. Desyatnikov, D. N. Neshev, W. Krolikowski, C. Denz, and Y. S. Kivshar, “Dynamic diffraction and interband transitions in two-dimensional photonic lattices,” Phys. Rev. Lett. 106(8), 083902 (2011). 20. J. Takahara, S. Yamagishi, H. Taki, A. Morimoto, and T. Kobayashi, “Guiding of a one-dimensional optical beam with nanometer diameter,” Opt. Lett. 22(7), 475–477 (1997). 21. Q. Xu, V. R. Almeida, R. R. Panepucci, and M. Lipson, “Experimental demonstration of guiding and confining light in nanometer-size low-refractive-index material,” Opt. Lett. 29(14), 1626–1628 (2004). 22. S. Maier, Plasmonics: Fundamentals and Applications (Springer, 2007). 23. F. Ye, D. Mihalache, B. Hu, and N. C. Panoiu, “Subwavelength plasmonic lattice solitons in arrays of metallic nanowires,” Phys. Rev. Lett. 104(10), 106802 (2010). 24. C. O. M. S. O. L. Multiphysics, www.comsol.com


Introduction
Rabi oscillations are the periodic transitions (or revivals) between two stationary states of a quantum system driven by the periodic external field.This is purely resonant effect, whose efficiency strongly depends on the frequency detuning between the external field and own oscillation frequency of the quantum system.Since their prediction by I. Rabi in his seminal paper [1], such revivals and related resonant phenomena have been studied in a variety of atomic, optical, and condensed matter systems [2][3][4][5].Especially fruitful was the extension of the concept of Rabi oscillations to various weakly guiding paraxial optical structures where light propagation is described by the famous Schrödinger equation that also governs the evolution of the wavefunction in quantum-mechanical systems [6][7][8][9][10][11][12][13][14][15][16].In such paraxial structures, Rabi oscillations manifest themselves as the periodic energy exchange between linear guided modes of the multimode structure, stimulated by weak periodic modulation of the refractive index in the direction of light propagation that plays the role similar to that of the external driving field for multilevel quantum system.Besides purely fundamental aspects, in optical systems the interest to Rabi oscillations is explained by the potential rich practical applications of this effect for the controllable shaping of laser radiation.Rabi oscillations were observed in dynamical long-periodic gratings created in optical fibers [7,8], they were studied in various shallow multimode waveguides [9] and modulated periodic lattices [10,11].Especially interesting are realizations of Rabi oscillations in the two-dimensional systems, where longitudinal refractive index modulations may couple light modes with different topological charges [12][13][14][15][16][17].Notice, that sometimes equations, similar to those describing dynamics of Rabi oscillations, may be encountered in periodic guiding structures without any longitudinal variations [18,19].
So far, the dynamics of optical Rabi oscillations was analyzed mostly in the paraxial multimode guiding structures, where all characteristic scales (beam and waveguide widths) substantially exceed the wavelength.In this case the coupled-mode approach applied to the Schrödinger equation, governing light evolution, yields equations for mode amplitudes analogous to those for populations of levels in a driven quantum system.On the other hand, last decade has witnessed rapid advances in nanotechnologies that already allow fabrication of optical waveguides with subwavelength dimensions [20,21].Considerable waveguide depths can be achieved (strong guiding regime), so that eigenmodes of such structures also become subwavelength.Especially fruitful is the combination of dielectrics and metals allowing to confine light due to the excitation of surface plasmon polaritons (SPPs) [22].In such structures one has to resort to solution of the full set of Maxwell's equations, since the paraxial approximation fails to describe light evolution.A general question arises -are Rabi oscillations possible in the modulated subwavelength structures and to which extent this effect can be used for shaping of narrow light beams?
In this paper we address Rabi oscillations in subwavelength dielectric and metal-dielectric waveguides and show that this effect persists at the subwavelength scales and can be used to control the modal structure of the output beams.Longitudinal modulation of the dielectric permittivity in the waveguide region results in a highly selective excitation of modes of proper parity, while efficiency of this process is controlled by the detuning of modulation frequency from the resonant value.Efficient mode conversion can be achieved at very short propagation distances 20 m m  .

The model and coupled mode theory
We consider the propagation of a TM-polarized light beam along the z -axis of a subwavelength waveguide with longitudinally modulated dielectric permittivity.The evolution of nonzero components of the electric x z ( ,0, ) and magnetic y (0, , 0) H = H fields for selected polarization state is governed by the reduced system of Maxwell's equations: , where we have excluded the longitudinal component ¶ ¶ of the electric field for convenience; 0 e and 0 m are the vacuum permittivity and permeability; and w is the frequency of light.The function ( , )  x z e describes the relative permittivity of the guiding structure.In the case of dielectric waveguide we set where bg e is the relative background permittivity.The function ( 2) describes the profile of Gaussian waveguide with width a and depth p , which oscillates along the propagation direction with the amplitude d and frequency z w .For convenience we denote dynamically varying part of dielectric permittivity as In the absence of longitudinal modulation ( 0) d = , all eigenmodes propagate independently and the power concentrated within each mode remains unchanged.Longitudinal permittivity modulation inside the waveguide couples different modes and stimulates energy exchange between them.Before we go to the numerical simulations of Eq. ( 1) and ( 2), it is instructive to consider the dynamics of this process using coupled-mode approach [6,23].In this approach we use Lorentz reciprocity formula .We substitute fields in such form into reciprocity formula, integrate it with respect to transverse coordinate x , and, taking into account the orthogonality condition z ( ) , arrive to the coupled-mode equations for amplitudes ( ) where z -dependent coefficients ( ) nm z  are given by: Taking into account variation of lin sin( ) z z de w  along the propagation distance, one can see that energy exchange between modes n and m will be most efficient if the resonance condition 0 m n z b b w - = is satisfied, provided that the modes have the same symmetry, so that corresponding integral in the numerator of (4) is nonzero.

Numerical results and discussion
Rabi oscillations occur in multimode guiding structures.We search numerically for eigenmodes of the unmodulated ( 0) d = subwavelength Gaussian waveguide in the form { ( , ), ( , )} . The components of vectors ( ), ( ) x x E H describe the transverse field distributions, while n b are the propagation constants of the modes with indexes 1,2,... n = .We consider an illustrative situation when the waveguide supports three guided modes.This occurs for the following parameters: bg 2.25 e = , 2 p= , and 250 nm a = . The width of the waveguide a is substantially smaller than the wavelength 632.8 nm l = . The profiles of the transverse x ( ) E and longitudinal z ( ) E electric field components in all three guided modes are shown in Figs.1(a) and 1(b), respectively.Modes are confined on the scale comparable with the wavelength l .The longitudinal field component is comparable in amplitude with the transverse field one.The symmetry of Maxwell's equations dictates that z E is antisymmetric for symmetric x E and that the number of nodes in z E always exceeds by one the number of nodes in x E .In the following, when we mention the parity of the mode, we refer to the shape of the x E component.The propagation constants of the eigenmodes depicted in Figs.1(a) and 1(b) are given by is the vacuum wavenumber.From the coupled mode equation [Eq.( 3)], due to nonzero coefficient 13  , the coupling is expected to occur for the first and third modes with equal parity at .Notice, that efficient coupling between first and second modes is impossible for symmetric permittivity perturbations ( , ) ( , )  x z x z e e -= , since corresponding coupling coefficient 12  vanishes (coupling of these modes requires perturbations that break the symmetry of e , such as transverse oscillations of the waveguide center).  1) with finite-element method [24].In Fig. 2(a) only the first mode was provided at the input, while in Fig. 2  The efficiency of mode conversion can be quantified using instantaneous weights of the modes n is accompanied by the growth of the weight of the third mode 3 n and vice versa, that indicates on resonant coupling of these two modes.At the same time, second mode is not excited and its weight 2 n remains nearly zero.While modal field structure is in fact recovered after one period of oscillations, radiation and backward reflections, unavoidable for modulation depths d and periods 2 / z p w considered here, slightly reduce the output quantity 1 2 3 n n n + + .This circumstance is not accounted for by coupled-mode approach and can be captured only by solution of the full set of Maxwell's equations.Radiative losses rapidly decrease with decrease of d .Stimulated mode conversion is the resonant effect.While in resonance nearly all power is transferred between modes, the efficiency of conversion drops down with increase of the detuning of modulation frequency from the resonant one.In Fig. 4(a) we show maximal weight of the third mode max For small d values resonant excitation of higher-order modes (in the case when structure supports several of them) can be made highly selective due to small resonance width.Another important quantity that characterizes the process considered here is the conversion length, or period of Rabi  oscillations c z , defined here as a distance at which most of the power returns back to the input eigenmode at z r w w = .In accordance with coupled-mode approach, where coupling distance evolves , one observes in Fig. 4(c) rapid decrease of c z with increase of the permittivity modulation depth d .One should stress that in this structure the complete transformation of the shape of subwavelength light spots can be achieved at propagation distances as short as 20 m  m , provided that the permittivity modulation depth d is large enough.
So far, Rabi oscillations were considered only in dielectric guiding structures.The question arises whether similar effect can be observed in metal-dielectric waveguides, where the presence of metal allows to easily achieve light confinement at the subwavelength scales.We thus consider metal-dielectric-metal waveguide, whose permittivity is given by Eq. ( 2 ) for 300 nm x > . The shapes of eigenmodes of such a waveguide strongly differ from those in purely dielectric structures and the modes are more confined, but one still can see that the first and third modes are symmetric, while second mode is antisymmetric [Figs.1(c) and 1(d)].We found that despite losses in metal, the resonant longitudinal permittivity modulation in the dielectric region still results in efficient coupling between first and third modes in such a structure, as shown in Figs.2(c) and 2(d).
Finally, Rabi oscillations are possible in the presence of focusing nonlinearity in the dielectric.We account for nonlinear response of the material by assuming that permittivity (2) includes an additional term  in (4)], and also shifts the exact resonance frequency to a higher value [Fig.4(d)].The latter fact can be explained taking into account that nonlinearity acts differently on propagation constants of eigenmodes, leading to higher growth of 1 b in comparison with 3 b , thereby increasing the difference 1 3 b b -.The nonlinear resonance curve becomes slightly asymmetric.

Conclusion
In this paper, we have demonstrated that the Rabi oscillations are possible in various deeply subwavelength guiding structures (with widths down to 100 nm ) including subwavelength dielectric waveguides and metal-dielectric-metal guiding structures, where the dynamics is governed by Maxwell's equations.Rabi oscillation for subwavelength modes may occur at very short propagation distances (about tens of microns), without considerable backward reflections, in both linear and nonlinear regimes.This enables controllable shaping of deeply subwavelength lights spots at extremely short distances.
in the modulated waveguide can be represented as a superposition eigenmodes with z -dependent amplitudes ( ) n a z .A similar representation can be used for u u { , } E H but with const n a =

Fig. 1 .
Fig. 1.Transverse (a) and longitudinal (b) electric field distributions in three linear eigenmodes of Gaussian waveguide with 2 p = and 250 nm a = .Transverse (c) and longitudinal (b) electric field distributions in three eigenmodes supported by the dielectric slab with 600 nm a = and 2 p = surrounded by metal with Re 20 m e = -.The efficient coupling of the first and third modes of the subwavelength waveguide due to resonant ( ) z r w w = longitudinal modulation is illustrated in Figs.2(a) and 2(b).The results are obtained by direct solution of Maxwell's Eqs.(1) with finite-element method[24].In Fig.2(a)only the first mode was provided at the input, while in Fig.2(b) only the third mode was provided.One can observe gradual transformation of the initial bell-shaped field distribution into three-hump distribution, characteristic for the third guided mode.In agreement with Eqs.(3), after accumulation of the power in the third mode, the reverse process starts and power flows (b) only the third mode was provided.One can observe gradual transformation of the initial bell-shaped field distribution into three-hump distribution, characteristic for the third guided mode.In agreement with Eqs.(3), after accumulation of the power in the third mode, the reverse process starts and power flows back into the first mode.This results in periodic Rabi oscillations having for 0.1 d = the period as short as 86 m m , irrespectively of the input mode index n .

Fig. 2 .
Fig. 2. (a),(b) Dynamic of resonant nonparaxial mode conversion at 0.1 d = in dielectric structure at 0 0.389 z r k w w = » , when only first or third mode is provided at the input.The propagation distance and transverse window size are 130 m m and 4 m m , respectively.(c),(d) Dynamics of mode conversion at 0.1 d = and 0 0.232 z r k w w = » in metal-dielectric-metal structure.The propagation distance and transverse window size are 75 m m and 2 m m , respectively.Only 2 xEdistribution is shown.
coefficients determining modal composition of the transverse component of the electric field x ( , ) E x z at the distance z , while x y ( ), ( ) n n E x H x characterize field distributions in eigenmodes.The dependencies ( ) n z n corresponding to dynamics from Figs. 2(a) and 2(b) are shown in Figs.3(a) and 3(b).One can see that decrease in 1

3 n
acquired upon propagation as a function of the modulation frequency in the case, when only first mode is launched into structure at 0 z = .The resonance is symmetric in z w -its width depends on the depth of the longitudinal permittivity modulation.The width of the resonance z dw , defined at the level d almost linearly [Fig.4(b)].

Fig. 4 .
Fig. 4. (a) Maximal weight of the third mode versus modulation frequency at 0.1 d = in the linear medium.Resonance width (b) and effective coupling distance at z r w w = (c) versus refractive index modulation depth d in the linear medium.(d) Maximal weight of the third mode versus modulation frequency in focusing medium for peak nonlinear contribution to the refractive index nl 0.02 n = and 0.1 d = .
n > is the nonlinear coefficient.Solution of the Maxwell's Eqs.(1) reveals that weak focusing nonlinearity reduces the conversion length c z [this is a direct consequence of the increase of the effective waveguide depth leading to growth of the overlap integral 13