Interaction of high order solitons with external dispersive waves

Interaction of high order solitons with external dispersive waves

I. Oreshnikov Department of Physics, St. Petersburg State University, 198504, Ulyanovskaya st. 1, Peterhof, St. Petersburg, Russian Federation    R. Driben ITMO University 197101, Kronverksky pr. 49, St. Petersburg, Russian Federation Department of Physics and CeOPP, University of Paderborn, Warburger Str. 100, D-33098 Paderborn, Germany
July 14, 2019

We have presented theoretical and numerical studies on interactions of dispersive waves with second order solitons. Disperive wave with considerable intensity resonantly colliding with the second order solitons can lead to acceleration/decceleration of the later with subsequent central frequency shifts, but still well preserve the oscillating structure of the 2-solitons. The 2-soliton creates an effective periodical refractive index profile and thus reflected and transmitted dispersive waves generate new spectral bands with sharp peak structures. When resonant dispersive waves interact with two second order solitons bounding them, multiple scattering with multiple frequency conversion of the radiation occurs. Thus 2-solitons and radiation trapped in between them can produce effective solitonic cavity with ”flat” or ”concave mirrors” depending on the intensity of the input.

Optical solitons attract attention of researches for decades due to the profound interest of their fundamental physics and technological applications book (). In addition to the fundamental solitons, integrable models and physical media described by nearly-integrable equations RMP () support high order -solitons, with , which are oscillating pulses periodically restoring their shape at distances that are multiples of the fundamental soliton period Satsu (). Experimental realization of 2- and 3-order optical solitons dates back to 1983 Stolen (). Initial narrowing of higher-order solitons was observed for soliton’s orders up to Mollenauer (). Higher-order solitons were also demonstrated in the cavity of a mode-locked dye laser Salin (). Strongly oscillating higher-order solitons find natural applications for the pulse compression Chan (); Li () and frequency conversion Lee (), while the breakup of high order solitons initiates the extremely important process of the supercontinuum generation SC ().

Propagation of optical pulse close to zero dispersion wavelength makes important consideration of high order dispersion terms and in particular the third order dispersion (TOD) term. Second order solitons (2-solitons) were found to be practically robust under the influence of moderate TOD Wai (), preserving their periodically oscillating nature. Very recently resonant radiation of 2-solitons was investigated DribenYulin () revealing fascinating structure of radiation band consisting of multiple frequency peaks. While emitting the mentioned above radiation 2-soliton remained its fundamental features for many periods. Interestingly, oscillating dissipative solitons featuring similar radiation behavior were reported too KudlinskiA (); KudlinskiB ().

In is relevant to mention that extensive studies were performed considering scattering of weak DWs on fundamental solitons, with soliton’s properties remaining invariant yulin (); pre (); Efimov (); Efimov2 (); Skryabinoverview (); Conforti (). On the other hand, it was also demonstrated that strong resonant collision of DWs with solitons can accelerate or decelerate solitons. These accelerations manifest themselves in frequency domain by upshifts or downshifts of central frequencies of the solitons DribenMitschke (); Demircan (); Demircan2 (); Tartara (). Generation of new bands of frequencies due to the scattering of DWs from solitons can be viewed as an alternative technique for generation of broad and coherent supercontinuum Demircan2 () without the high order soliton fission.

Taking one step further a very interesting scenario with DWs interacting with two co-propagating solitons Resonator (); Resonator2 (); Kudlinski (). DWs temporally located in between the solitons interact with one soliton after the other bouncing of them. This process can repeat itself many times resulting in creation of solitonic cavity behaving like two mirrors for the DWs. In case of DWs with considerable intensity, the two solitons get attracted one towards the other, resulting in collision or even fusion Fusion (). This phenomenon was proven to be responsible for the appearance of multiple soliton knot patterns DribenMitschke (); Resonator2 (); NC () during the complex supercontinuum generation process. Such a solitonic cavity was very recently successfully realized in experimental conditions Kudlinski (). Also a ”convex mirror” cavity with mutually repulsing via DW two dark solitons was demonstrated darksolitons ().

In the present work we aim to study interactions of 2-solitons with external dispersive waves (DWs), demonstrating a possibility to manipulate propagation trajectories of the 2-solitons, to create broad continuum with rich spectral properties and also to create solitonic cavities consisting of two 2-solitons and DWs bouncing between them.

Dynamics of optical pulse in the vicinity of zero dispersion wavelength (ZDW) is governed by nonlinear Schrödinger equation with the inclusion of TOD.


where and designate second and third order dispersion terms, while represents the nonlinear parameter. Considering standard fiber with and launching a soliton at central wavelength of , we will work with the following values of dispersion parameters: and . The nonlinear coefficient is taken .

The injected light consist of two pulses launched with different frequencies namely soliton and DW

The 2-soliton has the initial form

with temporal width of and the peak power of corresponding fundamental soliton . DW is given by

Here represent the amplitude of DW, the temporal delay of the DW with the respect to the soliton, is the temporal width of dispersive wave, and is the angular frequency deviation of the DW from the soliton given by .

The result of collision of a relatively wide DW () centered at and having a very small amplitude of with the 2-soliton is presented in Fig. 1(a, b). Fig. 1(a) demonstrates dynamics in temporal domain and it reveals that the DW is partially scattered and partially transmitted through the 2-soliton that creates an effective periodic refractive index potential. In general one can archive much higher percentage of scattering by optimizing the incident wavelength of the DW within the resonance band, but here we aim to inspect also the transmitted part of the DW. Looking at the spectral domain representation of the dynamics in Fig. 1(b) we can see several new frequency bands arising in a course of evolution. The leftmost band below consisting of multiple peaks is the resonant or Cherenkov radiation Akh () emitted form the 2-soliton. Its central spectral position as well as the separation between the frequency peaks is described by DribenYulin ():

Here is frequency detuning of the resonant wave from the soliton frequency, is the soliton velocity, is a period of the oscillations of 2-soliton and N is an integer number.

The region of transmitted and reflected radiation at the output of the fiber is shown in Fig. 1(c). The resonance frequencies for the radiation generated by FWM mixing between incident DW and 2-soliton were derived by using the perturbation theory similar to the analysis for the fundamental soliton in pre (). In contrast with the fundamental soliton case, the resonant radiation becomes polychromatic, with frequency detunings for incident and scattered fields being tied by the following relation


We solved the (2) graphically in order to obtain the spectral positions of the transmitted and reflected peaks and presented it in Fig. 1(d). For example the intersection of the blue curve with the horizontal magenta line represents the scattered and transmitted DWs of the zero order , that corresponds to the strongest reflected peak at about and strongest transmitted peak at . The intersection of the blue curve with the horizontal green line stand for and so on. Comparison between these results and those obtained by direct numerical simulations (Fig. 1(c)) is visualized by vertical dashed lines. We can find an excellent agreement between the two approaches with a discrepancy of only few nanometers.

Figure 1: (Color online) Scattering of a wide DW on 2-soliton. Dynamics in (a) temporal domain with shown instead of the conventional intensity in order to observe better less intense regions. (b) Spectral domain dynamics presented in logarithmic scale. (c) The spectral density at the end output of the fiber with the region of reflected and transmitted DW zoomed (from to ). (d) is the graphical solution to (2) predicting the position of the resonance frequencies. Blue curve is the left hand side of the equation; horizontal lines correspond to the right hand’s side for starting from the bottom.

Using more narrow incident DWs with the temporal width equal to that of the soliton () we have produced systematic simulations to provide the dependence of the scattered radiation central wavelength vs. the incident central wavelength as shown in Fig. 2. Solid curve is the plot of the analytical relation that is obtained from graphical solution of (2) for the fundamental harmonic . Shifting and thus modifying the place of collision with the respect to soliton’s shape at the moment of the collision did not significantly influence the results.

Figure 2: (Color online) Central wavelength of the scattered radiation as a function of incident. Solid line is the curve predicted by (2) for fundamental () harmonic. Square markers are found by numerical modeling.

Now we will inspect what happens when DW with more considerable intensity collides with the 2-soliton. In this case the velocity of the soliton will be modified and we can no longer rely on analytical approximations, but to study only results of numerical simulations. We collide DW with the 2-soliton as in the previous chapter, but this time we increase its input intensity to , that is still several times smaller that the input peak intensity of the 2-soliton participating in the collision and we will take the width of the DW equal to that of the soliton (). A typical result of collision of strong DW with the 2-soliton still resisting to the splitting is presented in the Fig. 3. The soliton’s trajectory is considerably shifted, consequently it central wavelength is downshifted from to . We can also observe in Fig. 3 (lower panel) generation of significant new light spectral bands in region between the soliton’s central wavelength and DW. Increasing the intensity of the incident DW much higher we will be observing the splitting of the 2-soliton.

Figure 3: (Color online) Interaction of a strong DW with 2-soliton. Dynamics in (a) temporal domain showing . (b) Spectral domain dynamics representation in logarithmic scale.

Even more interesting interaction scenario occurs when DW bounces off two co-propagating 2-solitons. We consider an input consisting of two well separated 2-solitons and DW launched in between them such as

If the wavelength of DW is located well inside the resonance band we will be able to observe its multiple reflections from the 2-solitons If the intensity of the DW is not very high the trajectories and the spectral characteristics of the two solitons are almost not modified and reflections occur pretty periodically to certain wavelength regions. For example the DW initially launched at with is being reflected several times to the region around and then back to its original spectral location (Fig. 4(a, b)). Interestingly the spectrum of the two solitons manifests multiple periodically alternating light and dark regions spanned for this particular example from to . Such an illumination spectrum that can be further controlled by the parameters of the participating solitons can be suggested for photonics applications requiring particular spectral design.

Figure 4: (Color online) Solitonic cavity created by two 2-solitons and DW bouncing between them. (a) Temporal domain representation with shown. (b) Spectral domain representation with output spectrum profile shown above the panel.

As mentioned above 2-solitons feature high degree of robustness when colliding with DWs and thus we can build a closed cavity consisting of two 2-solitons acting like concave mirrors without observing their splitting. An example of such cavity is demonstrated in Fig. 5(a, b). Here we used the same type of input light but increased an intensity of the incident DW to leaving the rest of the parameters unchanged. In this example the central wavelength of the left and right 2-solitons get upshifted and downshifted respectively. Also the periodicity of illuminated regions such as in Fig. 4 gets strongly distorted (Fig. 5(b)).

Figure 5: (Color online) Solitonic ”concave mirror” cavity created by two 2-solitons and stronger DW bouncing between them. (a) Temporal domain representation with shown. (b) Spectral domain representation.

To conclude, scattering of DWs on high order solitons was presented. The structure of the 2-soliton provides an effective periodical potential leading to polychromatic scattering and transmission for the incident DWs. In spectral domain it manifest itself by the generation of frequency comb like novel bands. Position of the peaks of the transmitted and reflected bands is predicted by excellently overlapping analytical and numerical approaches. Interaction of DWs with two 2-solitons bounding them and creation of solitonic cavities was also demonstrated. The presented nonlinear dynamics besides the fundamental interest suggest prospective photonic devices schemes for frequency conversion and broadband light generation with pre-designed spectrum.

Funding Information

R.D. gratefully acknowledges the support by the Russian Federation Grant 074-U01 through ITMO Early Career Fellowship scheme.


  • (1) Kivshar Y. S. and Agrawal G. P., Optical Solitons: From Fibers to Photonic Crystals (Academic Press, San Diego, 2003).
  • (2) Y. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 763-915, (1989).
  • (3) J. Satsuma and N. Yajima, Suppl. Progr. Theor. Phys. 55, 284–306 (1974).
  • (4) R. H. Stolen, L. F. Mollenauer, W. J. Tomlinson, Opt. Lett. 8, 186 (1983).
  • (5) L. F. Mollenauer, R. H. Stolen, J. P. Gordon, W. J. Tomlinson, Opt. Lett 8, 289 (1983).
  • (6) F. Salin, P. Grangier, G. Roger, A. Brun, Phys. Rev. Lett. 56, 1132 (1986).
  • (7) K. C. Chan and M. S. F. Liu, IEEE J. Quantum Electron. 31, 2226-2235 (1995)
  • (8) Qian Li, J.N. Kutz, P.K.A. Wai, J. Opt. Soc. Am. B, 27 2180–2189 (2010)
  • (9) K. S. Lee, J. A. Buck, J. Opt. Soc. Amer. B, Opt. Phys. 20, 514-519 (2003)
  • (10) J. M. Dudley, G. Genty, and S. Coen, Rev. Mod. Phys. 78, 1135-1184 (2006).
  • (11) P. K. A. Wai, C. R. Menyuk, Y. C. Lee, and H. H. Chen, Opt. Lett. 11, 464-466 (1986).
  • (12) A. Bendahmane, F. Braud., M. Conforti, B. Barviau, A. Mussot, and A. Kudlinski, Optica 1, 243-249, (2014).
  • (13) M. Conforti, S. Trillo, A. Mussot, and A. Kudlinski, Sci. Rep. 5 9433 (2015).
  • (14) R. Driben, A.V. Yulin, and A. Efimov, Opt. Express 23 (15), 19112-19117 (2015).
  • (15) A.V. Yulin, D.V. Skryabin and P. St. J. Russell, Optics Lett. 29, 2411-2413, (2004).
  • (16) D.V. Skryabin and A.V. Yulin, Phys. Rev. E 72, 016619 (2005).
  • (17) A. Efimov, A. V. Yulin, D. V. Skryabin, J. C. Knight, N. Joly, F. G. Omenetto, A. J. Taylor, and P. Russell, Phys. Rev. Lett. 95, 213902 (2005)
  • (18) A. Efimov, A. J. Taylor, A. V. Yulin, D. V. Skryabin, and J. C. Knight, Opt. Lett. 31(11), 1624–1626 (2006).
  • (19) D. V. Skryabin and A. V. Gorbach, Rev. Mod. Phys. 82, 1287–1299 (2010).
  • (20) A Bendahmane, A Mussot, M Conforti, A Kudlinski, Optics Express 23, 16595 (2015)
  • (21) R. Driben, F. Mitschke, and N. Zhavoronkov Optics Express, 18, 25993 (2010).
  • (22) A. Demircan, Sh. Amiranashvili, and G. Steinmeyer, Phys. Rev. Lett. 106, 163901 (2011).
  • (23) A. Demircan, S. Amiranashvili, C. Brée, and G. Steinmeyer, Phys. Rev. Lett.110(23), 233901 (2013).
  • (24) L. Tartara JOSA B 32, 395 (2015).
  • (25) A. V. Yulin, R. Driben, B. A. Malomed, and D. V. Skryabin, Opt. Exp. 21, 14481-14486 (2013).
  • (26) R. Driben, A.V. Yulin, A. Efimov and B.A Malomed, Opt. Express 21 (16), 19091–19096 (2013).
  • (27) S. F. Wang, A. Mussot, M. Conforti, X. L. Zeng, and A. Kudlinski, Optics Letters 40, 3320 (2015).
  • (28) R. Driben and I.V. Babushkin, Optics Letters 37, 24, 5157–5159 (2012).
  • (29) R. Driben, B. A. Malomed, A.V.Yulin, and D.V. Skryabin, Physical Review A 87, 063808 (2013).
  • (30) I. Oreshnikov, R. Driben, and A.V.Yulin,
  • (31) N. Akhmediev and M. Karlsson, Phys. Rev. A 51, 2602 (1995).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description