Three-dimensional Spatiotemporal Accessible Solitons in a PT-symmetric Potential
Three-dimensional Spatiotemporal Accessible Solitons in a PT-symmetric Potential
Journal of the Optical Society of Korea. 2012. Dec, 16(4): 425-431
Copyright ©2012, Optical Society of Korea
  • Received : August 17, 2012
  • Accepted : October 11, 2012
  • Published : December 25, 2012
Export by style
Cited by
About the Authors
Wei-Ping, Zhong
Department of Electronic and Information Engineering, Shunde Polytechnic, Guangdong Province, Shunde 528300, China
Milivoj R., Belić
Texas A&M University at Qatar, P.O. Box 23874 Doha, Qatar
Tingwen, Huang
Texas A&M University at Qatar, P.O. Box 23874 Doha, Qatar

Utilizing the three-dimensional Snyder–Mitchell model with a PT -symmetric potential, we study the influence of PT symmetry on beam propagation in strongly nonlocal nonlinear media. The complex Coulomb potential is used as the PT -symmetric potential. A localized spatiotemporal accessible soliton solution of the model is obtained. Specific values of the modulation depth for different soliton parameters are discussed. Our results reveal that in these media the localized solitons can exist in various shapes, such as single-layer and multi-layer disk-shaped structures, as well as vortex-ring and necklace patterns.
Recently, the study of systems exhibiting parity-time ( PT ) symmetry has drawn a great deal of attention. The underlying idea is to extend canonical quantum mechanics by introducing a class of non-Hermitian Hamiltonians which exhibit entirely real eigenvalue spectra below a certain phase-transition point [1] , although the potentials in these Hamiltonians are complex-valued. A necessary condition for the Hamiltonian to be PT -symmetric is that its potential V ( x ), being complex, is subject to the spatial-symmetry constraint V ( x ) = V * ( -x ). The complex PT -symmetric potentials can be realized in the most straightforward way in optics, by combining the spatial modulation of the refractive index with properly placed gain and loss [2] . This possibility has excited extensive theoretical [3 , 4] and experimental [5] research. Pioneering theoretical works [2 , 3 , 4] stimulated recent experimental studies that eventually resulted in the observation of the PT symmetry breaking in both active [5] and passive [6] optically coupled systems. This will probably enable manufacturing of integrated PT photonic devices with extraordinary capabilities, such as double-refraction or energy flow tailoring. A new direction in nonlinear optical research, concerning PT optical lattices [7] and the related PT -based solitons, can be envisaged. The existence and propagation dynamics of the one-dimensional (1D) optical solitons in a PT -symmetric linear periodic potential have been examined in detail in Ref. [8] . Hence, further exploration of general properties of solitons in multidimensional PT -symmetric potentials is warranted.
An intriguing feature of PT -symmetric potentials is the spontaneous breakdown of PT symmetry above a threshold level of the strength of the imaginary part of the potential. Above that level the eigenfunctions of the Hamiltonian cease to be the eigenfunctions of the PT operator, even though the PT symmetry remains in force. The complex Coulomb potential, which will be used as the PT potential in this paper, was among the early potentials studied in the PT -symmetric setting [9] . It turned out, however, that in 1D it cannot be treated on the real x -axis, but on some trajectories in the complex x plane [10 , 11] . In the multidimensional case it is treatable in the quantum mechanical sense, but as a part of a more complicated PT -symmetric potential involving parabolic and quartic potentials [11] . We treat it as a nonlinear optical problem, which may not make much sense quantum mechanically, but represents a viable system optically.
Optical spatial solitons, which are self-trapped optical beams that exist by virtue of the balance between diffraction and nonlinearity, have lately been extensively studied in nonlocal nonlinear media [12 , 13] . It has been found that the nonlocality can prevent the collapse of self-focusing beams in media with cubic nonlinearity [14] , suppress azimuthal instabilities of vortex solitons [15] , and stabilize Laguerre soliton clusters, azimuthons, and multipole solitons [12] .
The evolution of optical beams in nonlocal nonlinear media is governed by the nonlocal nonlinear Schrödinger (NNS) equation [12 , 16] . Of particular importance is the case, referred to as the strongly nonlocal case, in which the characteristic length of nonlocality is much larger than the beam width. In 1997, Snyder and Mitchell [16] simplified the NNS equation in the case of strong nonlocality to a linear model, called the Snyder-Mitchell model [13] . Subsequently, Assanto et al demonstrated theoretically [17] and experimentally [18] that nematic liquid crystals are some of the strongly nonlocal nonlinear (NN) media. The Snyder-Mitchell model can support solitons with new properties, the so-called “accessible solitons” [16] , which are described by the solutions of a linear differential equation, tantamount to the high-dimensional quantum harmonic oscillator [12 , 13] . Owing to this feature, strong nonlocality exerts a stabilizing influence on the dynamics of solutions, as the solutions of linear systems cannot be unstable or chaotic.
Light bullets (LBs), or optical spatiotemporal solitons [19] , in which both the diffraction and group-velocity dispersion are balanced by the nonlinearity, are challenging subjects in multidimensional nonlinear optics [20] . In addition to their fundamental significance as particle-like waves, light bullets can find applications in long and short-distance communications, all-optical switching, and digital computing, among others [21] . In this paper we demonstrate that a class of new 3D light bullets originating in a PT -symmetric potential can be supported by the strongly nonlocal nonlinear media. We display rather unusual properties of these 3D light bullets. However, the main aim of the present work is to perform a detailed study of the 3D NNS equation in the strongly NN media in the presence of a PT -symmetric potential and in a region where optical solitons can exist. Thus, we concentrate on obtaining the localized solutions of the 3D strongly NNS equation, without regarding them a priori as stationary eigenfunctions of the corresponding Hamiltonian operator, with the corresponding real or complex eigenvalues.
The rest of the article is organized as follows. In Section II, the 3D Snyder-Mitchell model with a PT -symmetric potential is introduced, and the localized soliton solutions are constructed. The properties of 3D localized accessible solitons with the PT -symmetric potential are explored. We also illustrate and discuss some examples of the exact solutions obtained in Section III. The article is concluded in Section IV.
We begin the analysis from the scaled (3+1)D spatiotemporal nonlinear Schrödinger equation [13 , 22 , 23]
Lager Image
which governs the propagation of a slowly varying field envelope u along the propagation coordinate z in a nonlinear nonlocal optical medium. Here ∇ 2 is the full 3D spatiotemporal Laplacian,
Lager Image
is the radial coordinate, z is the retarded time in the reference frame moving with the pulse, and
Lager Image
represents the nonlocal nonlinearity induced by the optical beam intensity I =| u | 2 . W denotes the external potential. We assume N to be of the form
Lager Image
Lager Image
is the normalized symmetric real response function of the medium whose characteristic length determines the degree of nonlocality. The setting chosen in this paper is reminiscent of the X wave generation geometry [24] ; the wave equations are similar and there exist natural linear and nonlinear regimes of wave packet dynamics. The major difference is that we consider a nonlocal medium with anomalous dispersion.
In the case of strong nonlocality, Eq. (1) can be simplified to the 3D Snyder-Mitchell model in spherical coordinates [13 , 22 , 25] :
Lager Image
where s is a parameter proportional to the beam power. To analyze PT -symmetry relevant for Eq. (2), we choose the complex Coulomb potential
Lager Image
where μ (≠0) is a real constant. Obviously, when μ = 0, Eq. (2) is simplified to the general 2D NNS equation in strongly NN media [25] . The second term in Eq. (2) represents diffraction, the third term originates from the optical nonlinearity, and the fourth term is the external potential function. It is noted that Eq. (2) is rather generic; one may assign different physical interpretations to essentially the same type of equation by choosing different physical systems and variables. For example, in a typical quantum mechanical setting one may understand Eq. (2) as the scaled Schrödinger equation for the wave function of a particle moving in the potential
Lager Image
An often analyzed similar PT version involves both quadratic and quartic anharmonic oscillator terms
Lager Image
[26] . Hence, it is the parameter μ that decides whether Eq. (2) is PT -symmetric or not. It is easy to see that V * ( -r ) = V ( r ). Thus, Eq. (2) is a system with the PT-symmetric potential.
We treat equation (2) in spherical coordinates, by the method of separation of variables. Defining the complex optical field as u ( z , r , θ , ϕ ) = F ( z , r ) Y ( θ , ϕ ), with separated angular variables, the separation yields the following two equations:
Lager Image
Lager Image
where l is a non-negative integer. Equation (3A) has the spherical harmonics as the solution,
Lager Image
Lager Image
are the associated Legendre polynomials with
Lager Image
and ϕ is the azimuthal angle. The parameter 0≤ q ≤1 determines the modulation depth of the beam intensity. The parameter m is a real nonnegative integer, called the topological charge.
Now, we consider the solution of Eq. (3B). Following Refs. [12 , 25] , we define F ( z , r ) = A ( z , r ) e iB (z,r) , where A ( z , r ) and B ( z , r ) are the real functions of z and r . With this variable change and after a little algebra, we transform Eq. (3B) into two coupled equations for A and B :
Lager Image
Lager Image
These equations can be treated by the self-similar method [25] . To treat Eqs. (5), the amplitude A ( z , r ) and the phase B ( z , r ) of the beam are further defined as [25] :
Lager Image
Lager Image
where Ω( z , r ) is the self-similarity variable and w ( z ) is the pulse width. As a consequence of these definitions, Eq. (5B) yields:
Lager Image
Lager Image
Note that the parameter b remains invariant on propagation. Furthermore, it should be noted that Eqs. (7) and (8) are universally applicable to all types of self-similar pulses. Here and in what follows, the symbols containing subscript “0” are used to represent the initial values of the corresponding parameters, at distance z = 0.
By means of Eqs. (6), (7) and (8), a nonlinear differential equation for F is derived from Eq. (5A),
Lager Image
where cz is the derivative of c with respect to z . In order to solve Eq. (9), we introduce another variable transformation
Lager Image
from Eq. (9) one obtains:
Lager Image
To simplify Eq. (10), yet another variable transformation is introduced, Ω 2 = R ; thus, we find:
Lager Image
In the end, to obtain tractable solutions, we restrict their generality by choosing special forms of some parameters:
Lager Image
Lager Image
where n is a nonnegative integer and w 0 (=constant) is the initial beam width of the pulse. Using Eq. (12B) we find that Eq. (11) becomes:
Lager Image
The solution of Eq. (13) can be written in terms of the confluent hypergeometric functions. Thus, when w = w 0 , an exact single accessible light bullet solution to Eq. (2) with a PT -symmetric potential can be written as:
Lager Image
where Ylm ( θ , ϕ ) are the spherical harmonics and 1 F 1 is the confluent hypergeometric function of the first kind. It is straightforward to see that | u ( z , r , θ , ϕ )| vanishes at r → ∞, i.e., Eq. (14) represents a localized solitary solution. Arbitrariness in the choice of the soliton parameters n , m and l included in the above solution (3) implies that the beam field u ( z , r , θ , ϕ ) may possess a rich structure. For μ = 0 the solution (14) goes to the solution (21) in Ref. [25] , apart from a constant factor. It should be noted that the solution (14) differs from the solution in the absence of the PT potential, i.e. when μ=0, by a complex factor exp [ i ( μr / w 0 - μ2z / 2 )]. Thus, the influence of the complex Coulomb potential is to modulate the accessible solitons that exist in the absence of it. The intensity distributions | unlm | 2 remain the same in both cases.
To better understand the 3D soliton dynamics, we introduce some special types of localized solutions for the
Lager Image
Intensity distributions of the LBs with the spherical structures for m = 0 and n = 0,2,4 from left to right, respectively.
Lager Image
LB distributions with disk and ring-shaped profiles for m = 0; the parameter are: (a) ; 2,l = 4; (b) n = 1, l = 4; (c) n = 0, l = 6.
optical field expressed by Eq. (14) via suitable selections of the nonnegative integer parameters ( n , m , l ). We focus attention on the distributions of the optical intensity I = | u | 2 . In the following examples, we further fix the beam width w 0 = 1.
First, we address the simplest case m = 0. Then the LB intensity does not depend on the modulation depth q . Because the parameters n and l are arbitrary, various light bullet structures can be obtained. If the parameters n and I are chosen as zero, from Eq. (14) we find that 1 F 1 (0,3/2, r 2 ) = 1; the beam is then called the fundamental light bullet and forms a single-layer sphere, see Fig. 1 (a). Keeping l = 0 and increasing the parameter n , multilayered structures are obtained; a typical example is presented in Fig. 1 (b) for n = 2. Similarly, we can construct a higher-order fundamental LB for larger n , i.e., for n = 4 the intensity distribution is exhibited in Fig. 1 (c). In general, there exist n +1 layered spheres for such a soliton.
For n ≠0 and l ≠0, the higher-order LBs can be excited. Figure 2 (a) displays the results for n = 2, l = 4, which features three coaxial rings in the mid-plane; there are five layers stacked along the τ -axis. Figure 2 (b) displays similar structures for n = 1, l = 4, with two rings in the mid-plane, and two disks above and below the central rings. Finally, Fig. 2 (c), corresponding to n = 0, l = 6, shows a single middle ring and, once again, ring and disk-shaped objects along the τ -axis.
For q = 1 in Eq. (14), we obtain the vortex ring beam for a nonnegative integer l and a positive integer m . As an interesting case we pick the parameters m = l (≠0). An example of such a vortex ring LB is shown in Fig. 3 . The soliton parameters are: (a) m = l = 4, n = 0; (b) m = l = 3, n = 1; (c) m = l = 2, n = 2. The number of rings in the horizontal direction is determined by n . Changing the modulation depth from q = 1 to 0 < q < 1, we find that the LBs modulate azimuthally, see Fig. 4 . It is noted that the outer modulation is more distinct than the inner one. Actually, the formation of a vortex ring beam is the result of the periodic azimuthal modulation functions cos( ) and sin( ).
For m > 0 ( m ≠ l) an integer, and in the limit q = 1,
Lager Image
Intensity profiles of the vortex ring LBs for q = 1 and m = l. (a) m = l = 4, n = 0; (b) m = l =3, n = 1; (c) m = l =2, n = 2.
Lager Image
Intensity distribution of LBs from Fig. 3. Setup is the same as in Fig. 3, except for q = 0.95.
Lager Image
Vortex-ring solitons, for q = 1 and m ≠ l. Parameters (n,l,m) have the following values: (a) (0,4,1); (b) (1,4,2); (c) (2,4,2). The figure layout is as in Fig. 3.
the multilayered vortex ring LBs along the vertical τ -axis are found. In Fig. 5 we depict some properties of the vortex ring LBs. It is seen that for the same m , the larger the parameter n , the larger the soliton radius in the horizontal plane. The optical intensity is zero at the τ -axis, which is the location of the topological defect.
For q = 0 and l = m in Eq. (14), we obtain single-layer and multi-layer necklace beams for a positive integer m . A typical example of such a necklace is shown in Fig. 6 for l = m = 4, along with the axisymmetric distribution.
Self-trapped localized structures with a large number of azimuthal petals and multi-layered necklaces may exhibit a strong effective stabilization in strongly NN media [12 , 13] . Figure 7 displays the intensity distribution of multi-layer
Lager Image
Single- and multi-layer necklace LBs in the horizontal plane for l = m = 4 and n = 0,1,2 from left to right. The setup is as in Fig. 3, except for q = 0.
Lager Image
Structures of multi-layer necklace solitons in the vertical direction. The setup is the same as in Fig. 6, except for l ≠ m.
necklace solitons in the vertical direction, which exhibit similar patterns. These examples are obtained for positive integers ( n , l , m ) in Eq. (14). The parameters are: (a) (0,4,2); (b) (1,4,3); (c) (2,4,3). In these solutions, the necklace structure is still formed, due to the periodic azimuthal modulation. Note that these solitons form multi-layers, with the outer ones more strongly modulated than the inner counterparts.
Interesting structures are seen in Figs. 6 and 7 . We find that the larger the parameter m , the larger the necklace radius. It is seen that the soliton distributions change regularly with the azimuthal angle. The number of beads in each layer is determined by m , and the number of layers is determined by n . These solitons contain 2 m ( n +1) necklaces and form n +1 necklace layers in the horizontal direction. The number of necklace layers in the vertical direction is determined by l .
We have introduced a class of self-trapped LB solutions of the NNS equation in the strongly NN media with a PT -symmetric potential in the form of a complex Coulomb potential. We were not much concerned with the details of the eigenvalue spectrum, but with the solution of the NNS equation with a specific PT -symmetric potential. Analytical accessible soliton solutions are obtained with the help of the self-similar method for solving such evolution partial differential equations. They are given in terms of the confluent hypergeometric function of the first kind and spherical harmonics. We find that in addition to the fundamental LBs, these solutions may come in the form of 3D single-layer and multi-layer disk-shaped, vortex ring and necklace LBs.
This work was supported by the National Natural Science Foundation of China under Grant No. 61275001 and by the Natural Science Foundation of Guangdong Province, China, under Grant No. 1015283001000000. The work at the Texas A&M University at Qatar is supported by the NPRP 09-462-1-074 project of the Qatar National Research Fund.
Bender C. M. , Boettcher S. (1998) “Real spectra in non- Hermitian Hamiltonians having PT symmetry” Phys. Rev. Lett. 80 5243 - 5246    DOI : 10.1103/PhysRevLett.80.5243
El-Ganainy R. , Makris K. G. , Christodoulides D. N. , Musslimani Z. H. (2007) “Theory of coupled optical PT-symmetric structures” Opt. Lett. 32 2632 - 2634    DOI : 10.1364/OL.32.002632
Christodoulides D. N. , Lederer F. , Silberberg Y. (2003) “Discretizing light behavior in linear and nonlinear waveguide lattices” Nature (London) 424 817 - 823    DOI : 10.1038/nature01936
Abdullaev F. Kh. , Kartashov Y. V. , Konotop V. V. , Zezyulin D. A. (2011) “Solitons in PT-symmetric nonlinear lattices” Phys. Rev. A 83 041805-1 - 041805-4    DOI : 10.1103/PhysRevA.83.041805
Ruter C. E. , Makris K. G. , EI-Ganainy R. , Christodoulides D. N. , Segev M. , Kip D. (2010) “Observation of parity-time symmetry in optics” Nat. Phys. 6 192 - 195    DOI : 10.1038/nphys1515
Guo A. , Salamo G. J. , Duchesne D. , Morandotti R. , Volatier-Ravat M. , Aimez V. , Siviloglou G. A. , Christodoulides D. N. (2009) “Observation of PT-symmetry breaking in complex optical potentials” Phys. Rev. Lett. 103 093902-1 - 093902-4    DOI : 10.1103/PhysRevLett.103.093902
Makris K. G. , El-Ganainy R. , Christodoulides D. N. , Musslimani Z. H. (2010) “PT-symmetric optical lattices” Phys. Rev. A 81 063807-1 - 063807-10    DOI : 10.1103/PhysRevA.81.063807
Musslimani Z. H. , Makris K. G. , El-Ganainy R. , Christodoulides D. N. (2008) “Optical solitons in PT periodic potentials” Phys. Rev. Lett. 100 030402-1 - 030402-4    DOI : 10.1103/PhysRevLett.100.030402
Znojil M. , Levai G. (2000) “The Coulomb-harmonic oscillator correspondence in PT symmetric quantum mechanics” Phys. Lett. A 271 327 -    DOI : 10.1016/S0375-9601(00)00400-X
Levai G. (2009) “Spontaneous breakdown of PT symmetry in the complex Coulomb potential” Pramana 73 329 - 335    DOI : 10.1007/s12043-009-0125-5
Lévai G. , Siegl P. , Znojil M. (2009) “Scattering in the PTsymmetric Coulomb potential” J. Phys. A 42 295201 -    DOI : 10.1088/1751-8113/42/29/295201
Zhong W. P. , Yi L. (2007) “Two-dimensional Laguerre-Gaussian soliton family in strongly nonlocal nonlinear media” Phys. Rev. A 75 061801-1 - 061801-4    DOI : 10.1103/PhysRevA.75.061801
Zhong W. P. , Yi L. , Xie R. H. , Belic M. , Chen G. (2008) “Robust three-dimensional spatial soliton clusters in strongly nonlocal media” J. Phys. B: At. Mol. Opt. Phys. 41 025402 -    DOI : 10.1088/0953-4075/41/2/025402
Bang O. , Krolikowski W. , Wyller J. , Rasmussen J. J. (2002) “Collapse arrest and soliton stabilization in nonlocal nonlinear media” Phys. Rev. E 66 046619-1 - 046619-5    DOI : 10.1103/PhysRevE.66.046619
Buccoliero D. , Desyatnikov A. S. , Krolikowski W. , Kivshar Y. S. (2008) “Spiraling multivortex solitons in nonlocal nonlinear media” Opt. Lett. 33 198 - 200    DOI : 10.1364/OL.33.000198
Snyder A. , Mitchell J. (1997) “Accessible solitons” Science 276 1538 - 1541    DOI : 10.1126/science.276.5318.1538
Conti C. , Peccianti M. , Assanto G. (2003) “Route to nonlocality and observation of accessible solitons” Phys. Rev. Lett. 91 073901-1 - 073901-4    DOI : 10.1103/PhysRevLett.91.073901
Conti C. , Peccianti M. , Assanto G. (2004) “Observation of optical spatial solitons in a highly nonlocal medium” Phys. Rev. Lett. 92 113902-1 - 113902-4    DOI : 10.1103/PhysRevLett.92.113902
Silberberg Y. (1990) “Collapse of optical pulses” Opt. Lett. 15 1282 - 1284    DOI : 10.1364/OL.15.001282
Malomed B. A. , Mihalache D. , Wise F. , Torner L. (2005) “Spatiotemporal optical solitons” J. Opt. B: Quantum Semiclassical Opt. 7 R53 - R72    DOI : 10.1088/1464-4266/7/5/R02
Abdullaev F. K. , Konotop V. V. 2004 Nonlinear Waves: Classical and Quantum Aspects Kluwer Academic Publishers Dordrecht, Netherlands
Zhong W. P. , Belić M. , Xie R. , Huang T. , Lu Y. (2010) “Three-dimensional spatiotemporal solitary waves in strongly nonlocal media” Opt. Commun. 283 5213 - 5217    DOI : 10.1016/j.optcom.2010.08.004
Zhong W. P. , Belić, M. , Assanto G. , Malomed B. A. , Huang T. (2011) “Light bullets in the spatiotemporal nonlinear Schrödinger equation with a variable negative diffraction coefficient” Phys. Rev. A 84 043801-1 - 043801-8    DOI : 10.1103/PhysRevA.84.043801
Faccio D. , Averchi A. , Couairon A. , Kolesik M. , Moloney J. V. , Dubietis A. , Tamosauskas G. , Polesana P. , Piskarskas A. , Di Trapani P. (2007) “Spatio-temporal reshaping and X wave dynamics in optical filaments” Opt. Express 15 13077 - 13095    DOI : 10.1364/OE.15.013077
Zhong W. P. , Belic M. (2009) “Three-dimensional optical vortex and necklace solitons in highly nonlocal nonlinear media” Phys. Rev. A 79 023804-1 - 023804-6    DOI : 10.1103/PhysRevA.79.023804
Znojil M. (2000) “Quasi-exactly solvable quartic potentials with centrifugal and Coulombic terms” arXiv:math-ph/0002036v2    DOI : 10.1088/0305-4470/33/22/320