Symmetric and asymmetric solitons supported by a PT-symmetric ...

Vol. 26, No. 6 | 19 Mar 2018 | OPTICS EXPRESS 6949

Symmetric and asymmetric solitons supported by a PT-symmetric potential with saturable nonlinearity: bifurcation, stability and dynamics P ENGFEI L I , 1,* C HAOQING DAI , 2 R UJIANG L I , 3

AND

YAQIN G AO 4

1 Department

of Physics, Taiyuan Normal University, Taiyuan, 030031, China of Sciences, Zhejiang Agriculture and Forestry University, Lin’an 311300, China 3 College of Information Science and Electronic Engineering, Zhejiang University, Hangzhou, 310027, China 4 Department of Physics, Taiyuan University of Science and Technology, Taiyuan 030024, China 2 School

* [email protected]

Abstract: The symmetry breaking bifurcation of solitons in an optical waveguide with focusing saturable nonlinearity and parity-time (PT)-symmetric complex-valued external potentials is investigated. As the soliton power increases, it is found that the branches of asymmetric solitons split off from the base branches of PT-symmetric fundamental soliton. The bifurcation diagrams, consisting essentially of the propagation constants of optical solitons, indicate that symmetric fundamental and multipole solitons, as well as asymmetric solitons can exist. The stabilities and the dynamics characteristics of solitons are comprehensively investigated. We find the different instability scenarios of the symmetric solitons, but the symmetry breaking bifurcation is caused only by the onset of instability of the symmetric fundamental solitons. This result is further confirmed by the numerical examples with the different saturable nonlinearity parameters. In particular, we find that the soliton power and the stability of soliton at the bifurcation points are significantly changed by varying the strength of the saturable nonlinearities. These results provide additional way to control symmetry breaking bifurcations in PT-symmetric optical waveguide. © 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement OCIS codes: (190.0190) Nonlinear optics; (190.6135) Spatial solitons; (350.5500) Propagation.

References and links 1. C. M. Bender and S. Boettcher, “Real spectra in non-Hermitian Hamiltonians having P T symmetry,” Phys. Rev. Lett. 80(24), 5243–5246 (1998). 2. C. M. Bender, D. C. Brody, and H. F. Jones, “Complex extension of quantum mechanics,” Phys. Rev. Lett. 89(27), 270401 (2002). 3. C. M. Bender, S. Boettcher, and P. N. Meisinger, “PT-symmetric quantum mechanics,” J. Math. Phys. 40(5), 2201–2229 (1999). 4. R. El-Ganainy, K. G. Makris, D. N. Christodoulides, and Z. H. Musslimani, “Optical P T-symmetric structures,” Opt. Lett. 32(17), 2632–2634 (2007). 5. S. Klaiman, U. Gunther, and N. Moiseyev, “Visualization of branch points in P T-symmetric waveguides,” Phys. Rev. Lett. 101(8), 080402 (2008). 6. V. S. Bagnato, D. J. Frantzeskakis, P. G. Kevrekidis, B. A. Malomed, and D. Mihalache, “Bose-Einstein condensation: twenty years after,” Rom. Rep. Phys. 67(1), 5–50 (2015). 7. A. Lupu, H. Benisty, and A. Degiron, “Switching using PT symmetry in plasmonic systems: positive role of the losses,” Opt. Express, 21(18) 021651–21668 (2013). 8. H. Alaeian, and J. A. Dionne, “Non-Hermitian nanophotonic and plasmonic waveguides,” Phys. Rev. B 89(7), 075136 (2014). 9. Y. Choi, J. K. Hong, J. H. Cho, K. G. Lee, J. W. Yoon, and S. H. Song, “Parity-time-symmetry breaking in double-slab surface-plasmon-polariton waveguides,” Opt. Express, 23(9) 11783–11789 (2015). 10. H. Alaeian, and J. A. Dionne, “Parity-time-symmetric plasmonic metamaterials,” Phys. Rev. A 89(3), 033829 (2014). 11. N. Lazarides, and G. P. Tsironis, “Gain-driven discrete breathers in P T-symmetric nonlinear metamaterials,” Phys. Rev. Lett. 110(5), 053901 (2013). 12. N. M. Chtchelkatchev, A. A. Golubov, T. I. Baturina, and V. M. Vinokur, “Stimulation of the fluctuation superconductivity by P T symmetry,” Phys. Rev. Lett. 109(15), 150405 (2012).

#319669 Journal © 2018

https://doi.org/10.1364/OE.26.006949 Received 15 Jan 2018; revised 21 Feb 2018; accepted 22 Feb 2018; published 7 Mar 2018


13. A. Ruschhaupt, F. Delgado, and J. G. Muga, “Physical realization of P T-symmetric potential scattering in a planar slab waveguide,” J. Phys. A: Math. Gen. 38(9), L171–L176 (2005). 14. A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. Volatier-Ravat, V. Aimez, G. A. Siviloglou, and D. N. Christodoulides, “Observation of P T-symmetry breaking in complex optical potentials,” Phys. Rev. Lett. 103(9), 093902 (2009). 15. C. E. Rüter, K. G. Makris, R. El-Ganainy, D. N. Christodoulides, M. Segev, and D. Kip, “Observation of parity-time symmetry in optics,” Nat. Phys. 6(3), 192–195 (2010). 16. C. M. Huang, F. W. Ye, Y. V. Kartashov, B. A. Malomed, and X. F. Chen, “P T symmetry in optics beyond the paraxial approximation,” Opt. Lett. 39(18), 5443–5446 (2014). 17. Z. H. Musslimani, K. G. Makris, R. El-Ganainy, and D. N. Christodoulides, “Optical solitons in P T periodic potentials,” Phys. Rev. Lett. 100(3), 030402 (2008). 18. F. K. Abdullaev, Y. V. Kartashov, V. V. Konotop, and D. A. Zezyulin, “Solitons in P T-symmetric nonlinear lattices,” Phys. Rev. A 83(4), 041805(R) (2011). 19. X. Zhu, H. Wang, L. Zheng, H. Li, and Y. J. He, “Gap solitons in parity-time complex periodic optical lattices with the real part of superlattices,” Opt. Lett. 36(14), 2680–2682 (2011). 20. H. G. Li, Z. W. Shi, X. Jiang, and X. Zhu, “Gray solitons in parity-time symmetric potentials,” Opt. Lett. 36(16), 3290–3292 (2011). 21. S. M. Hu, X. Ma, D. Lu, Z. Yang, Y. Zheng, and W. Hu, “Solitons supported by complex P T-symmetric Gaussian potentials,” Phys. Rev. A, 84(4), 043818 (2011). 22. B. Midya, and R. Roychoudhury, “Nonlinear localized modes in P T-symmetric Rosen-Morse potential wells,” Phys. Rev. A 87(4), 045803 (2013). 23. S. Nixon, L. J. Ge, and J. K. Yang, “Stability analysis for solitons in P T-symmetric optical lattices,” Phys. Rev. A 85(2), 023822 (2012). 24. D. A. Zezyulin, and V. V. Konotop, “Nonlinear modes in the harmonic P T-symmetric potential,” Phys. Rev. A 85, 043840 (2012). 25. V. Achilleos, P. G. Kevrekidis, D. J. Frantzeskakis, and R. Carretero-González, “Dark solitons and vortices in P T-symmetric nonlinear media: From spontaneous symmetry breaking to nonlinear P T phase transitions,” Phys. Rev. A 86(1), 013808 (2012). 26. M. A. Miri, A. B. Aceves, T. Kottos, V. Kovanis, and D. N. Christodoulides, “Bragg solitons in nonlinear P Tsymmetric periodic potentials,” Phys. Rev. A 86(3), 033801 (2012). 27. P. F. Li, B. Liu, L. Li, and D. Mihalache, “Nonlinear parity-time-symmetry breaking in optical waveguides with complex Gaussian-type potentials,” Rom. J. Phys. 61(3–4), 577–594 (2016). 28. U. A. Khawaja, S. M. Al-Marzoug, H. Bahlouli, and Y. S. Kivshar, “Unidirectional soliton flows in P T-symmetric potentials,” Phys. Rev. A 88(2), 023830 (2013). 29. P. F. Li, L. Li, and B. A. Malomed, “Multisoliton Newton’s cradles and supersolitons in regular and parity-timesymmetric nonlinear couplers,” Phys. Rev. E 89(6), 062926 (2014). 30. X. F. Zhu, Y. G. Peng, and D. G. Zhao, “Anisotropic reflection oscillation in periodic multilayer structures of parity-time symmetry,” Opt. Express, 22(15), 18401–18411 (2014). 31. C. Q. Dai, and X. G. Wang, “Light bullet in parity-time symmetric potential,” Nonlinear Dyn. 77(4), 1133–1139 (2014). 32. Y. Chen, and Z. Y. Yan, “Multi-dimensional stable fundamental solitons and excitations in PT-symmetric harmonicGaussian potentials with unbounded gain-and-loss distributions,” Commun. Nonlinear. Sci. Numer. Simulat. 57, 34–46 (2018). 33. J. K. Yang, “Symmetry breaking of solitons in one-dimensional parity-time-symmetric optical potentials,” Opt. Lett. 39(19), 5547–5550 (2014). 34. J. K. Yang, “Can parity-time-symmetric potentials support families of non-parity-time-symmetric solitons,” Stud. Appl. Math. 132(4), 332–353 (2014). 35. J. K. Yang, “Symmetry breaking of solitons in two-dimensional complex potentials,” Phys. Rev. E 91(2), 023201 (2015). 36. S. M. Hu, and W. Hu, “Defect solitons in saturable nonlinearity media with parity-time symmetric optical lattices,” Physica B. 429, 28–32 (2013). 37. P. Cao, X. Zhu, Y. J. He, and H. G. Li, “Gap solitons supported by parity-time-symmetric optical lattices with defocusing saturable nonlinearity,” Opt. Commun, 316, 190–197 (2014). 38. Y. V. Kartashov, V. A. Vysloukh, and L. Torner, “Soliton trains in photonic lattices,” Opt. Express. 12(13), 2831–2837 (2004). 39. M. Matuszewski, B. A. Malomed, and M. Trippenbach, “Spontaneous symmetry breaking of solitons trapped in a double-channel potential,” Phys. Rev. A 75(6), 063621 (2007). 40. A. Sacchetti, “Universal critical power for nonlinear Schrödinger equations with symmetric double well potential,” Phys. Rev. Lett. 103(19), 194101 (2009). 41. J. K. Yang, “Classification of solitary wave bifurcations in generalized nonlinear Schrödinger equations,” Stud. Appl. Math. 129(2), 133–162 (2012). 42. J. K. Yang, and T. I. Lakoba, “Universally-Convergent Squared-Operator Iteration Methods for Solitary Waves in General Nonlinear Wave Equations,” Stud. Appl. Math. 118(2), 153–197 (2007).


43. Z. G. Chen, M. Segev, T. H. Coskun, D. N. Christodoulides, and Y. S. Kivshar, “Coupled photorefractive spatial-soliton pairs,” J. Opt. Soc. Am. B. 14(11), 3066–3077 (1997). 44. J. W. Fleischer, T. Carmon, and M. Segev, “Observation of discrete solitons in optically induced real time waveguide arrays,” Phys. Rev. Lett. 90(2), 023902 (2003).

1.

Introduction

Parity-time (PT) symmetric systems originate from an observation in non-Hermitian quantum mechanics and they have attracted a growing interest in recent years. A simple one-dimensional PT-symmetric system can be described by the Schrödinger equation with a complex potential, where the real part of the potential is an even function of the coordinates and the imaginary part is an odd function. Although PT-symmetric systems are non-Hermitian, they possess all-real eigenvalue spectra, as demonstrated in seminal works by Bender and Boettcher [1–3]. The analysis was initially motivated by the development of quantum theories, then it soon spread out to the various PT-symmetric physical systems including optics and photonics [4, 5], Bose-Einstein condensates [6], plasmonic waveguides and meta-materials in which losses appear due to the metallic absorption [7–11], and superconductivity [12], etc. In optical waveguides, the propagation of light can often be described theoretically by Schrödinger-type equation under the paraxial approximation, where the refractive indices of optical waveguides can be described by real-valued potential functions and the imaginary parts of complex potentials stand for the effects of losses and gains in the optical waveguides. A judicious balance of gain and loss constitutes a PT-symmetric optical waveguide, then the propagation of light can be governed by a Schrödinger equation with a PT-symmetric Hamiltonian. PT-symmetric optical system, constructed by two coupled optical waveguides with equal gain and loss, was firstly investigated theoretically and experimentally in the linear regime [4, 13–15]. Beyond the paraxial approximation, rigorous theory based on Maxwell’s equations was also presented [16]. When nonlinear effects are involved in the PT-symmetric optical waveguides, the interplay between PT-symmetry and nonlinear effect will give rise to novel phenomena, such as PT-symmetric stationary nonlinear modes [17–27] and novel nonlinear dynamic properties [28–30]. Furthermore, theoretical studies were extended to the PT-symmetric spatiotemporal localizations. Spatiotemporal solitons (light bullets) were found in PT-symmetric potentials with cubic nonlinearity and inhomogeneous diffraction [31], as well as an inhomogeneous nonlinear medium with unbounded gain and loss distributions [32]. Recently, non-PT-symmetric solitons (asymmetric solitons) in PT-symmetric potentials have aroused great attention. Different from the aforementioned solitons, they possess asymmetric profiles due to the symmetry breaking of solitons. It has been shown that this symmetry breaking cannot occur in generic PT-symmetric potentials. But for a special class of PT-symmetric potentials, a family of stable PT-symmetry breaking solitons with real eigenvalue spectra have been found [33]. For this type of potentials, a precondition of the existence of non-PTsymmetric solitons is miraculously satisfied [34, 35], where the model describing the symmetry breaking bifurcations of PT-symmetric solitons is mainly focused on cubic nonlinearity. In theoretical works, the model has been considered in PT-symmetric optical waveguides with saturable nonlinearities [36, 37], however, only the properties of the symmetric solitons have been investigated. In previous papers, the symmetry breaking of PT-symmetric solitons was decided by the shapes or the parameters of the PT-symmetric potentials. The main aim of this paper is to discuss how the effect of the saturable nonlinearity can vary the symmetry breaking bifurcation, which could provide an additional way to transform the soliton from a symmetric to an asymmetric profile in PT-symmetric system. In the following sections, we will carry out a detailed investigation of symmetric and asymmetric solitons in a PT-symmetric potential with focusing saturable nonlinearity and study the key features of two types of solitons in one dimensional Schrödinger equation. Especially, we will show that the points of the symmetry


breaking are to be significantly changed by adjusting the strength of the saturable nonlinearity and present the dependence of the soliton power at the point of symmetry breaking bifurcation on the modulation strength of the potential. This may provide another way to control the symmetry breaking bifurcations in PT-symmetric optical waveguides by changing the strength of the saturable nonlinearity. The paper is organized as follows. In Sect. 2, the governing model is introduced. In Sect. 3, we show that asymmetric solitons, as well as the fundamental and multipole symmetric solitons can exist in one dimensional Schrödinger equation with the focusing saturable nonlinearity. The bifurcation phenomena and the dependence between the propagation constants and the input power are discussed. Numerical results on the dependence of the soliton power at the point of symmetry breaking bifurcation on the modulation strength of the potential are also presented. In Sect. 4, the stability of symmetric and asymmetric solitons with different parameters of saturable nonlinearity are analyzed systematically, and the corresponding nonlinear evolution dynamics are studied by performing direct numerical simulations. Finally, we will drawn our conclusions in Sect. 5. 2.

The governing model

In paraxial regime, optical wave propagation in a planar graded-index waveguide with focusing saturable nonlinearity is governed by the following one dimensional normalized nonlinear Schrödinger equation ∂ψ ∂ 2 ψ |ψ| 2 ψ + + U ψ + σ i = 0, (1) (ξ) ∂ζ ∂ξ 2 1 + S |ψ| 2 where ψ (ζ, ξ) is the normalized optical field envelope function, U (ξ) ≡ V (ξ) + iW (ξ) is the normalized complex potential, σ = 1 represents the focusing nonlinearity, and S is the saturable parameter standing for the degree of saturable nonlinearity [38]. When the parameter of focusing nonlinearity σ = 0, Eq. (1) describes a linear PT-symmetric system, where it has purely real eigenvalues below an exceptional point and a pair of complex conjugated eigenvalues above the point. Only for the purely real eigenvalues, the Hamiltonian commutes with PT operator. When the parameter of focusing nonlinearity σ = 1, we focus on the solutions with physical relevance, i.e., the solutions are stationary, and the corresponding eigenvalues are all real. To obtain optical solitary wave solutions, we assume the stationary solutions in the form of ψ (ζ, ξ) = φ (ξ) eiβζ , where φ (ξ) ≡ φ R + iφ I is a complex function with the real part φ R and imaginary part φ I , and β is the propagation constant. Substitution into Eq. (1) yields d2 φ |φ| 2 φ + U φ + σ − βφ (ξ) = 0. (ξ) (ξ) dξ 2 1 + S |φ| 2

(2)

In this paper, we consider a special class of PT-symmetric potential U (ξ) = g 2 (ξ) + i

dg (ξ) , dξ

where g (ξ) is a real and even function. We take g (ξ) in the form of ξ + ξ0 ξ − ξ0 g (ξ) = W0 sech + sech , χ0 χ0

(3)

(4)

where W0 represents the modulation strength of the potential in Eq. (3), ξ0 and χ0 are related to the separation and width of PT-symmetric potential, respectively. Then the real and imaginary parts of the PT-symmetric potential can be constructed as 2 ξ + ξ0 ξ − ξ0 2 V (ξ) = W0 sech + sech , (5) χ0 χ0


W0 ξ + ξ0 ξ + ξ0 ξ − ξ0 ξ − ξ0 W (ξ) = − sech tanh + sech tanh . χ0 χ0 χ0 χ0 χ0

(6)

1.8

3

(a)

V( W(

2

g(

)

(b)

) )

1.2 SS1

1

SS2 SS3

0.6

AS

0

-1 -15

0.6

0.0 -10

-5

0

5

10

15

(c)

0.5

0.8

1.0

1.5 P

2.0

2.5

(d)

R

0.4

3.0

R

I

I

0.4

0.2

0.0

-15

0.0

-10

-5

0

5

10

15

-15

-10

-5

0

5

10

15

Fig. 1. Diagram of bifurcation and the corresponding symmetric and asymmetric solitons. (a) Profiles of the real and even function g(ξ) (black dashed curve), real (blue solid curve) and imaginary (thin red solid curve) components of the potential. (b) The propagation constants versus the soliton power, where the symmetric fundamental solitons (blue solid curve labeled with SS1), symmetric dipole solitons for the excited state (black dash-dot curve labeled with SS2), symmetric tripole solitons for the excited state (thin gray solid curve labeled with SS3), as well as asymmetric solitons (red dotted curve labeled with AS) are shown. (c) and (d) Distributions of fundamental symmetric and asymmetric solitons with the soliton power P = 1.5, where the blue solid and dashed red curves represent the real part and imaginary part, respectively. Here, the parameters are ξ0 = 2, χ0 = 1, and W0 = 1.5.

3.

Solitons and symmetry breaking bifurcations

In this section, we aim to search for the symmetric and asymmetric solitons, and illustrate the characteristics of the bifurcations of symmetry breaking in the special class of PT-symmetric potential in Eq. (3). First, we consider the focusing nonlinearity with the saturable parameter S = 1. As an example, the real and even function g (ξ) is chosen in the form of a double-hump hyperbolic secant function, and the PT-symmetric potential in Eqs. (5) and (6) is in the form of a Scarff-II potential with double-hump profiles. We consider the PT-symmetric potential with parameters ξ0 = 2, χ0 = 1, and W0 = 1.5. The profiles of the real function g (ξ), real and imaginary parts of the PT-symmetric potential are shown in Fig. 1(a). The propagation constants and the corresponding solitons can be obtained by solving numerically the stationary Eq. (2) with the power-conserving squared-operator iteration method [42]. Numerical results indicate that there exist two families of solitons, i.e., PT-symmetric solitons (symmetric solitons) and non-PT-symmetric solitons (asymmetric solitons). The diagrams of bifurcation, consisting essentially of the propagation constants ∫ +∞ of solitons, show that there only exist symmetric solitons at a low soliton power (P = −∞ |ψ| 2 dξ). The solitons with symmetric profiles are found, including the symmetric fundamental solitons (denote with SS1), the symmetric dipole and tripole solitons for the excited states (denote with SS2 and SS3). With the increase of soliton power, a new branch splits off from the base branch of PT-symmetric fundamental solitons at the point of P = 0.7. Beyond this point, the symmetry of the fundamental


PT-symmetric solitons are broken with the increasing of soliton power, but the symmetry of the multipole soliton is still preserved. The propagation constants of symmetric and asymmetric solitons (denote with AS) form a pitchfork bifurcations shown in Fig. 1(b), which is also called as symmetry-breaking bifurcations in the models with real potential [39–41]. We present the profiles of symmetric and asymmetric solitons with soliton power P = 1.5 in Fig. 1(c) and Fig. 1(d), respectively. It is a natural result for the symmetric fundamental soliton that the complex-valued function possesses symmetric real profile and antisymmetric imaginary profile, as shown in Fig. 1(c). For asymmetric soliton, the profiles of the real and imaginary parts are all asymmetric, as shown in Fig. 1(d).

Fig. 2. Profiles of solitons for different soliton power. (a) Symmetric fundamental solitons, (b) symmetric dipole solitons for the excited state, (c) symmetric tripole solitons, (d) asymmetric solitons. Here, the parameters are the same as those in Fig. 1.

Figure 2 shows the symmetric and asymmetric solitons with the soliton power varying from 0.1 to 3. Fig. 2(a) presents the profiles of the fundamental solitons. The profiles of the symmetric dipole and tripole solitons for the excited states are shown in Fig. 2(b) and Fig. 2(c), respectively. It is worth to note that the profiles of the symmetric fundamental solitons and the symmetric dipole solitons for the excited state are all double humps, but the profiles of the real and imaginary parts have an opposite symmetry. The profiles of asymmetric solitons are presented in Fig. 2(d), where the results show that asymmetric solitons can exist ranging from 0.7 to 3 for the parameter of soliton power. With the increasing of soliton power, the profiles of asymmetric solitons are gradually concentrated into one of the humps of the PT-symmetric potential. To further illustrate the effect of the focusing saturable nonlinearity on the symmetry breaking of solitons, the diagrams of bifurcation are depicted with the different saturable nonlinear parameters. In Fig. 3(a), the diagram of bifurcation is presented with a weak saturable nonlinearity S = 0.1. The fundamental solitons are symmetric when the soliton power is below 0.5 (shown with blue solid curve), and the asymmetric solitons begin to emerge when the soliton power exceeds 0.5 (shown with red dotted curve). The point of symmetry breaking bifurcation appears at P = 0.5 as shown by the black vertical dashed line. The profiles of solitons, corresponding to the blue solid curve for P 6 0.5 and red dotted curve in Fig. 3(a), present a transformation of the fundamental solitons from symmetric to asymmetric distributions, as shown in Fig. 3(d). Diagram of bifurcation with moderate saturable nonlinearity (S = 0.5) is depicted in Fig. 3(b), where the point of symmetry breaking bifurcation appears at P = 0.6. For comparison, diagram


Fig. 3. Diagrams of bifurcation and profiles of solitons for different saturable nonlinear parameters. (a), (b) and (c) are diagrams of bifurcations with saturable nonlinear parameters S = 0.1, S = 0.5, and S = 1, respectively, where the vertical dashed lines show the soliton powers at bifurcation points. (d), (e) and (f) depict that the profiles transform from the symmetric fundamental solitons into asymmetric solitons due to symmetry breaking. The vertical plane shows the soliton power at bifurcation points, where the saturable nonlinear parameters are S = 0.1, S = 0.5, and S = 1, respectively. Here, the other parameters are the same as those in Fig. 1.

of bifurcation with saturable nonlinear parameter S = 1 is also plotted in Fig. 3(c) with the point of symmetry breaking bifurcation P = 0.7. The corresponding solitons of a transformation from symmetric to asymmetric distributions are also shown in Fig. 3(e) and Fig. 3(f), respectively. These numerical results indicate that the soliton power at the bifurcation points can be changed by the different values of the saturable parameter. Therefore, the the symmetry breaking bifurcation may be controlled by the different strengths of the saturable nonlinearities instead of the shape and parameter of the PT-symmetric potential. It is also relevant to address the dependence of the symmetry breaking bifurcation on the modulation strength of the potential W0 with different saturable nonlinear parameters. In Fig. 4, the top row shows the symmetry breaking bifurcations with modulation strength of the potential W0 = 2, where Figs. 4(a), 4(b) and 4(c) present the diagrams of bifurcation with S = 0.1, 0.5 and 1, respectively. The bottom row presents the symmetry breaking bifurcations with modulation strength of the potential W0 = 2, where Figs. 4(d), 4(e) and 4(f) present the diagrams of bifurcation with S = 0.1, 0.5 and 1, respectively. The numerical results demonstrate that the soliton power at the point of symmetry breaking bifurcation increases with the decrease of the modulation strength of the potential W0 , which indicates that the symmetry of soliton needs a much higher soliton power to break for the weak modulation strength of the potential.


4

3.4

3.6

(a) 3.6

(c)

(b)

W0=2

W0=2

W0=2

S=0.1

S=0.5

S=1

3.2

β

3

3.2

2.8

2.8

1

2

3

1

2

3

1 1.3

2

3

2

3

(f)

W0=1

W0=1

S=0.5

S=0.1

1

1

0.9

(e)

(d) W0=1

β

2.6

S=1

0.8 0.7

0.7 0.6

1

2

3

1

P

2

P

3

0.5

1

P

Fig. 4. Diagrams of bifurcation depending on the different modulation strength of the potential W0 with different saturable nonlinear parameters. (a), (b) and (c) depict the diagrams of bifurcation with modulation strength of the potential W0 = 2, and the focusing saturable nonlinearity parameters S = 0.1, 0.5 and 1, respectively. The red dotted curves and blue solid curves represent the propagation constants of symmetric and asymmetric solitons. (d), (e) and (f) show the diagrams of bifurcation with modulation strength of the potential W0 = 1, the focusing saturable nonlinearity parameters S = 0.1, 0.5 and 1, respectively. Here, the other parameters of potential are the same as those in Fig. 1.

4.

Linear stability analysis and dynamics

In order to examine the stability of these solitons, it is necessary to explore the perturbed solutions for a known solution φ (ξ) with the form h i ∗ ψ (ζ, ξ) = eiβζ φ (ξ) + u (ξ) e δζ + v ∗ (ξ) e δ ζ , (7) where φ (ξ) is the stationary solution with a real propagation constant β, u (ξ) and v (ξ) are small perturbations with |u|, |v| |φ|. Substituting Eq. (7) into Eq. (1) and keeping only the linear terms, we obtain the following linear eigenvalue problem L11 L12 u u i =δ . (8) L21 L22 v v Here, L11 , L12 , L21 and L22 are expressed as the following forms L11 =

d2 2σ |φ| 2 σS |φ| 4 − + U − β + 2 , dξ 2 1 + S |φ| 2 1 + S |φ| 2 L12 =

2σφ2 1 + S |φ| 2

−

σS |φ| 2 φ2 2 , 2 1 + S |φ|

(9)

(10)


σS |φ| 2 φ∗2 2σφ∗2 L21 = , 2 − 1 + S |φ| 2 1 + S |φ| 2 L22 = −

(11)

2σ |φ| 2 d2 σS |φ| 4 ∗ − U + β − + 2 , dξ 2 1 + S |φ| 2 1 + S |φ| 2

(12)

where δ is a complex eigenvalue of the linear problem for Eq. (8). The positive real part of complex eigenvalue can be used to measure the instability growth rate of the perturbation. If δ contains a real part, it is indicated that the soliton φ (ξ) is linearly unstable, otherwise, φ (ξ) is linearly stable.

Max(δR )

0.4 0.3

(a) SS1 SS2 SS3 AS

0.2 0.1 0.5

Max(δR )

0.3

1.5 P

2

2.5

(b)

3

SS1 SS2 SS3 AS

0.2 0.1 0.5

Max(δR )

1

0.2

1

1.5 P

2

2.5

(c)

3

SS1 SS2 SS3 AS

0.1

0.5

1

1.5 P

2

2.5

3

Fig. 5. Eigenvalues of linear stability analysis for the different parameters of saturable nonlinearity. (a), (b) and (c) depict the largest real part of δ for the symmetric and asymmetric solitons with focusing saturable nonlinearity parameters S = 0.1, 0.5 and 1, respectively. Here, the parameters are the same as in Fig. 3.

Figure 5 presents the dependence of the largest real part Max(δR ) of the eigenvalues of linear stability analysis for the symmetric and asymmetric solitons on the soliton power P with different saturable nonlinearity parameters S. The largest real part of the eigenvalues of linear stability analysis for the symmetric fundamental solitons, the symmetric dipole and tripole solitons for the excited states, as well as the asymmetric solitons are plotted with blue solid curves, black short dot-dashed curves, gray dashed curves and red curves with circles for S = 0.1, S = 0.5 and S = 1 in Figs. 5(a), 5(b) and 5(c), respectively. For the saturable parameter S = 0.1, all of the symmetric solitons are linearly stable when the soliton power is lower than 0.5, then the value of


Max(δR ) begins to increase, which indicates that the symmetric fundamental solitons turn into a unstable region ranging from P = 0.5 to 3. It is in this region that the asymmetric solitons can coexist with the symmetric solitons, as shown in Fig. 5(a). The symmetric dipole solitons for the excited states are linearly stable in the region of P = 0.5 to 3, but the symmetric tripole solitons possess a weak unstable region ranging from P = 1.6 to 3. Asymmetric solitons become unstable when the soliton power increases upto P = 2.5. The results of the linear stability analysis for S = 0.5 and S = 1 are shown in Figs. 5(b) and 5(c), where the symmetric fundamental solitons turn into a unstable region at the points of P = 0.6 and P = 0.7, respectively. We also found that the stable regions for symmetric tripole solitons and asymmetric solitons shrink with the increasing of the saturable nonlinearity parameter. By comparing Fig. 3 with Fig. 5, it is found the asymmetric solitons always exist in that regions where the symmetric fundamental solitons begin to turn into an unstable region. This indicates that the above mentioned a special class of PT-symmetric potential is only one of the necessary conditions for the symmetry-breaking bifurcations, and the symmetry breaking has been induced by the instability of the symmetric fundamental solitons. Furthermore, the points of symmetry breaking bifurcation are changed by the different values of the saturable parameter in Fig. 3, but essentially, these bifurcation points are caused by the onset of instability of the symmetric fundamental solitons with different saturable parameters, as shown those in Fig. 5. The results of linear stability analysis for solitons are further confirmed by employing direct numerical simulations of Eq. (1). Numerical results are summarized in Fig. 6 and Fig. 7 with the saturable parameter S = 1. The eigenvalue spectra of Eq. (8) for symmetric and asymmetric solitons are shown in the left column in Fig. 6 with the soliton power P = 1, where the eigenvalue spectra of linear stability analysis for the symmetric solitons are displayed in Figs. 6(a), 6(b) and 6(c). It is found that the real component (δR ) and imaginary component (δI ) of the eigenvalues of linear stability analysis appear in pairs for the symmetric solitons, so these eigenvalue spectra are quartet symmetry. The eigenvalues of linear stability analysis also appear in pairs for the asymmetric soliton, however, the real components of these eigenvalues are asymmetric, as depicted in Fig. 6(d). From the eigenvalue spectra of linear stability analysis, we can found that the symmetric fundamental soliton is strongly linearly unstable in Fig. 6(a). The corresponding evolution plot shows that the optical field intensities oscillate periodically between two humps of the PT-symmetric potential after a short distance of stationary evolution, as shown in Fig. 6(e). The symmetric solitons for the excited states, as well as the asymmetric soliton, are weakly and linearly unstable. To confirm that, the initial profiles are perturbed by a 5% random noise. The results of the direct numerical simulations indicate that these solitons can propagate stably at the initial stage, then the amplitudes of soliton appear small oscillations due to the weak instability, and these small oscillations are continually amplified with the increasing of propagation distance which eventually leads to the broken-down profiles, as presented in Figs. 6(f), 6(g) and 6(h), respectively. The eigenvalue spectra and the dynamics evolutions of symmetric and asymmetric solitons are also presented in Fig. 7 with different soliton power. By comparing Fig. 7(a) with Fig. 6(a), we can see that the symmetric fundamental soliton have a larger instability growth rate for the soliton power P = 2.5. The evolution plot shows that the symmetric fundamental soliton has eventually broken down after oscillating periodically in Fig. 7(e). Similar to Fig. 6(b), the symmetric dipole soliton for the excited state is still weakly and linearly unstable for the soliton power P = 2.5. The eigenvalue spectrum of linear stability analysis and the corresponding evolution perturbed by a 5% random noise are shown in Figs. 7(b) and 7(f), respectively. The stability of the symmetric tripole soliton and the asymmetric soliton decrease when the soliton power P = 2.5, as shown in Figs. 7(c) and 7(d). The evolution plots present these unstable propagations, even though the initial profiles haven’t been perturbed by the random noise, as shown in Figs. 7(g) and 7(h). Finally, we would like to note that, the control of symmetry breaking bifurcations can be


Fig. 6. The eigenvalue spectra of linear stability analysis for the solitons with saturable parameter S = 1 and the corresponding evolution plots. (a), (b) and (c) show the eigenvalue spectra of linear stability analysis for SS1, SS2 and SS3, (d) shows the eigenvalue spectrum of AS with the soliton power P = 1, respectively. (e), (f), (g), and (h) The corresponding evolution plots of the field intensity. Here the other parameters are the same as in Fig. 1.

realized experimentally by considering the waveguide fabrication with photorefractive crystals, where the PT-symmetric Scarff-II potential with double-hump profiles in Eqs. (5) and (6) can be realized by Ti in-diffusion, gain and loss can be realized with Fe-doped [15]. The strength of the saturable nonlinearity can be dynamically adjusted with an applied voltage in the biased photorefractive crystals [43, 44]. For the controllable symmetry breaking bifurcations of solitons, it includes two steps: generating a symmetric or an asymmetric soliton and adjusting the point of symmetry breaking bifurcation. Then soliton can be transformed from a symmetric to an asymmetric profile by varying the strength of the saturable nonlinearity, and vice versa. 5.

Conclusions

In conclusion, we have performed a systematic study for the symmetry breaking of solitons in PT-symmetric optical waveguides with focusing saturable nonlinearity. Our analysis has shown that the branches of asymmetric soliton bifurcated from the base branches of PT-symmetric fundamental soliton with the increasing of the soliton power, causing by the onset of instability of the symmetric fundamental solitons. The stabilities of the symmetric and asymmetric soliton solutions have been comprehensively analyzed by employing linear stability analysis. By comparing numerical examples of the


Fig. 7. The eigenvalue spectra of linear stability analysis for the solitons with saturable parameter S = 1 and the corresponding evolution plots. (a), (b) and (c) show the eigenvalue spectra of linear stability analysis for SS1, SS2 and SS3, (d) shows the eigenvalue spectrum of AS with the soliton power P = 2.5. (e), (f), (g), and (h) The corresponding evolution plots of the field intensity. Here the other parameters are the same as those in Fig. 1.

different saturable nonlinearity parameters, we found that the above-mentioned a special class of PT-symmetric potentials is a necessary condition for the symmetry breaking bifurcations. Under such a condition, the symmetry breaking of the solitons have been essentially caused by the instabilities of the fundamental solitons. Especially, these numerical results of the linear stability analysis have also indicated that the soliton power at the bifurcation point can be changed by the different saturable nonlinearity parameters. The different instability scenarios of solitons, including the eigenvalue spectra of linear stability analysis and the corresponding evolution plots, have been revealed by using direct numerical simulations. The results have shown that the stability of the symmetric and asymmetric solitons decrease in PT-symmetric optical waveguides with the focusing saturable nonlinearity when the soliton power increases. Thereby, a further study with these effects is of great interest for the controllable symmetry breaking bifurcations of solitons in PT-symmetric optical waveguides with saturable nonlinearity. Extension of this symmetry breaking in saturable nonlinear media to higher spatial dimensions is also an important direction.


Funding Doctoral Scientific Research Foundation of Taiyuan Normal University (Grant No. I170144); Zhejiang Provincial Natural Science Foundation of China (Grant No. LY17F050011); Open Fund of IPOC (BUPT); Doctoral Scientific Research Foundation of Taiyuan University of Science and Technology (Grant No. 20152043); National Natural Science Foundation of China (Grant No. 11747063). Disclosures The authors declare that there are no conflicts of interest related to this article.