# Bright Solitons in aPT-Symmetric Chain of Dimers

Authors: Omar B. Kirikchi; Alhaji A. Bachtiar; Hadi Susanto
Journal: Advances in Mathematical Physics
 (2016)

---

## Abstract

We study the existence and stability of fundamental bright discrete solitons in a parity-time- (PT-) symmetric coupler composed by a chain of dimers that is modelled by linearly coupled discrete nonlinear Schrödinger equations with gain and loss terms. We use a perturbation theory for small coupling between the lattices to perform the analysis, which is then confirmed by numerical calculations. Such analysis is based on the concept of the so-called anticontinuum limit approach. We consider the fundamental onsite and intersite bright solitons. Each solution has symmetric and antisymmetric configurations between the arms. The stability of the solutions is then determined by solving the corresponding eigenvalue problem. We obtain that both symmetric and antisymmetric onsite mode can be stable for small coupling, in contrast to the reported continuum limit where the antisymmetric solutions are always unstable. The instability is either due to the internal modes crossing the origin or the appearance of a quartet of complex eigenvalues. In general, the gain-loss term can be considered parasitic as it reduces the stability region of the onsite solitons. Additionally, we analyse the dynamic behaviour of the onsite and intersite solitons when unstable, where typically it is either in the form of travelling solitons or soliton blow-ups.


---

## Body

## 1. Introduction

A system of equations isPT-symmetric if it is invariant with respect to combined parity (P) and time-reversal (T) transformations. The symmetry is interesting as it forms a particular class of non-Hermitian Hamiltonians in quantum mechanics [1], which may have a real spectrum up to a critical value of the complex potential parameter, above which the system is in the “broken PT symmetry” phase [2–4].The most basic configuration havingPT symmetry is a dimer, that is, a system of two coupled oscillators where one of the oscillators has damping losses and the other one gains energy from external sources. Considerably, dimers are also the most important PT systems as the concept of PT symmetry was first realised experimentally on dimers consisting of two coupled optical waveguides [5, 6] (see also the review [7] for PT symmetry in optical applications). The experiments have been rapidly followed by many other observations of PT symmetry in different branches of physics, from mechanical to electrical analogues (see the review [8]).When nonlinearity is present in aPT system, one may have nontrivial behaviours that do not exist in the linear case, such as the presence of blow-up dynamics in the parameter region of the unbroken phase in the linear counterpart [9–11]. When nonlinear dimers are put in arrays where elements with gain and loss are linearly coupled to the elements of the same type belonging to adjacent dimers, one can also obtain a distinctive feature in the form of the existence of solutions localized in space as continuous families of their energy parameter [12]. The system therefore has two arms with each arm described by a discrete nonlinear Schrödinger equation with gain or loss. Here, we study the nonlinear localised solutions, which loosely we also refer to as bright discrete solitons, and their stability analytically and numerically.In the continuous limit, the coupled equations without gain-loss have been studied in [13–16], where it has been shown that the system admits symmetric, antisymmetric, and asymmetric solitons between the arms. Unstable asymmetric solutions bifurcate from the symmetric ones through a subcritical symmetry breaking bifurcation, which then become stable after a tangent (saddle-center) bifurcation. When one adds a gain and loss term in each arm, one obtains PT-symmetric couplers, which have been considered in [17–22]. In the presence of the linear-gain and loss terms, asymmetric solitons cease to exist, while antisymmetric solitons are always unstable [20], even though those with small amplitudes can live long due to weak underlying instability [17]. Symmetric solitons can be stable in a similar fashion to those in the system without gain-loss [20].The stability of bright discrete solitons inPT-symmetric couplers was discussed in [12] using variational methods, where it was shown that symmetric onsite solutions can be stable and there is a critical solution amplitude above which the PT symmetry is broken. The case when the polarity of the PT-symmetric dimers is staggered along the chain is considered in [23]. The same equations without gain and loss were considered in [24] where the symmetric soliton loses its stability through the symmetry-breaking bifurcation at a finite value of the energy, similarly to that in the continuous counterpart [13–16]. Recently, a similar PT chain of dimers with a slightly different nonlinearity was derived [25] to describe coupled chains of parametrically driven pendula as a mechanical analogue of PT-symmetric systems [26]. The stability of bright discrete solitons was established through the applications of the Hamiltonian energy and an index theorem. The nonlinear long-time stability of the discrete solitons was also established using the Lyapunov method in the asymptotic limit of a weak coupling between the pendula [27].In this work, we determine the eigenvalues of discrete solitons inPT-symmetric couplers analytically using asymptotic expansions. The computation is based on the so-called method of weak coupling or anticontinuum limit. The application of the method in the study of discrete solitons was formulated rigorously in [28] for conservative systems. It was then applied to PT-symmetric networks in [29, 30]. However, no explicit expression of the asymptotic series of the eigenvalues for the stability of discrete solitons has been presented before. Here, in addition to the asymptotic limit of weak coupling between the dimers, we also propose to consider expansions in the coefficient of the gain-loss terms. In this case, explicit computations of the asymptotic series of the eigenvalues become possible.The manuscript is outlined as follows. In Section2, we present the mathematical model. In Section 3, we use perturbation theory for small coupling to analyse the existence of fundamental localised solutions. Such analysis is based on the concept of the so-called anticontinuum limit approach. The stability of the solitons is then considered analytically in Section 4 by solving a corresponding eigenvalue problem. In this section, in addition to small coupling, the expansion is also performed under the assumption of small coefficient of the gain-loss term due to the nonsimple expression of the eigenvectors of the linearised operator. The findings obtained from the analytical calculations are then compared with the numerical counterparts in Section 5. We also produce stability regions for all the fundamental solitons numerically. In this section, we present the typical dynamics of solitons in the unstable parameter ranges by direct numerical integrations of the governing equation. We present the conclusion in Section 6.


## 2. Mathematical Model

The governing equations describingPT-symmetric chains of dimers are of the form [12](1)u˙n=iσun2un+iϵΔ2un+γun+ivn,v˙n=iσvn2vn+iϵΔ2vn-γvn+iun.The derivative with respect to the evolution variable (i.e., the propagation distance, if we consider their application in fiber optics) is denoted by the overdot, un=un(t) and vn=vn(t) are complex-valued wave function at site n∈Z, ϵ>0 is the constant coefficient of the horizontal linear coupling (coupling constant between two adjacent sites), Δ2un=(un+1-2un+un-1) and Δ2vn=(vn+1-2vn+vn-1) are the discrete Laplacian term in one spatial dimension, the gain and loss acting on complex variables un and vn are represented by the positive coefficient γ; that is, γ>0. The nonlinearity coefficient is denoted by σ, which can be scaled to +1 without loss of generality due to the case of focusing nonlinearity that we consider. Bright discrete soliton solutions satisfy the localisation conditions un,vn→0 as n→±∞.The focusing system has static localised solutions that can be obtained from substituting(2)un=Aneiωt,vn=Bneiωtinto (1) to yield the equations(3)ωAn=An2An+ϵAn+1-2An+An-1-iγAn+Bn,ωBn=Bn2Bn+ϵBn+1-2Bn+Bn-1+iγBn+An,where An and Bn are complex-valued and the propagation constant ω∈R.


## 3. Solutions of Weakly Coupled Equations

In the uncoupled limit, that is, whenϵ=0, chain (1) becomes the equations for the dimer. The static equation (3) has been analysed in detail in [29, 31], where it was shown that there is a relation between ω and γ above which there is no time-independent solution to (3) (see also the analysis below). When ϵ is nonzero, but small enough, the existence of solutions emanating from the uncoupled limit can be shown using the Implicit Function Theorem. The existence analysis of [25] can be adopted here despite the slightly different nonlinearity as the Jacobian of our system when uncoupled shares a rather similar invertible structure (see also [29, 30] that have the same nonlinearity in the governing equations but different small coupling terms). However, below we will not state the theorem and instead derive the asymptotic series of the solutions.Using perturbation expansion, solutions of the coupler (3) for small coupling constant ϵ can be expressed analytically as(4)An=An0+ϵAn1+⋯,Bn=Bn0+ϵBn1+⋯.By substituting the above expansions into (3) and collecting the terms in successive powers of ϵ, one obtains at O(1) and O(ϵ), respectively, the equations(5)An0=Bn0ω-Bn0Bn∗0-iγ,Bn0=An0ω-An0An∗0+iγ,(6)An1=Bn1ω-2Bn0Bn∗0-iγ-Bn02Bn∗1-Δ2Bn0,Bn1=An1ω-2An0An∗0+iγ-An02An∗1-Δ2An0.It is well-known that there are two natural fundamental solutions representing bright discrete solitons that may exist for anyϵ, from the anticontinuum to the continuum limit, that is, an intersite (two-excited-site) and onsite (one-excited-site) bright discrete mode. Here, we will limit our study to these two fundamental modes.


### 3.1. Intersite Soliton

In the uncoupled limit, the mode structuresAn(0) and Bn(0) for the intersite soliton are of the form(7)An0=a^0eiϕan=0,1,0otherwise,Bn0=b^0eiϕbn=0,1,0otherwise,with [31](8)a^0=b^0=ω∓1-γ2,sin⁡ϕb-ϕa=γ,which is an exact solution of (5). Note that (8) will have no real solution when |γ|>1. This is the broken region of PT-symmetry. The parameter ϕa can be taken as 0, due to the gauge phase invariance of the governing equation (1) and henceforth ϕb=arcsin⁡γ and π-arcsin⁡γ. The former phase corresponds to the so-called symmetric configuration between the arms, while the latter is called antisymmetric one. Herein, we also refer to the symmetric and antisymmetric soliton as soliton I and II, respectively. Equation (8) informs us that ω>1-γ2>0 and ω>-1-γ2 are the necessary conditions for solitons I and II, respectively.For the first-order correction due to the weak coupling, writingAn(1)=a~n,1eiϕa and Bn(1)=b~n,1eiϕb and substituting them into (6) will yield(9)
a
~
n
,
1
=
b
~
n
,
1
=
1
2
a
^
0
n
=
0,1
,
1
a
^
0
n
=
-
1,2
,
0
otherwise
.Equations (4), (7), (8), and (9) are the asymptotic expansion of the intersite solitons. One can continue the same calculation to obtain higher order corrections. Here, we limit ourselves to the first-order correction only, which is sufficient to determine the leading order behaviour of the eigenvalues later.


### 3.2. Onsite Soliton

For the onsite soliton, that is, a one-excited-site discrete mode, one can perform the same computations to obtain the mode structure of the form(10)An0=a^0eiϕan=0,0otherwise,Bn0=b^0eiϕbn=0,0otherwise,with (8). After writing An(1)=a~n,1eiϕa and Bn(1)=b~n,1eiϕb, the first order correction from (6) is given by(11)a~n,1=b~n,1=1a^0n=0,±1,0otherwise.


## 3.1. Intersite Soliton

In the uncoupled limit, the mode structuresAn(0) and Bn(0) for the intersite soliton are of the form(7)An0=a^0eiϕan=0,1,0otherwise,Bn0=b^0eiϕbn=0,1,0otherwise,with [31](8)a^0=b^0=ω∓1-γ2,sin⁡ϕb-ϕa=γ,which is an exact solution of (5). Note that (8) will have no real solution when |γ|>1. This is the broken region of PT-symmetry. The parameter ϕa can be taken as 0, due to the gauge phase invariance of the governing equation (1) and henceforth ϕb=arcsin⁡γ and π-arcsin⁡γ. The former phase corresponds to the so-called symmetric configuration between the arms, while the latter is called antisymmetric one. Herein, we also refer to the symmetric and antisymmetric soliton as soliton I and II, respectively. Equation (8) informs us that ω>1-γ2>0 and ω>-1-γ2 are the necessary conditions for solitons I and II, respectively.For the first-order correction due to the weak coupling, writingAn(1)=a~n,1eiϕa and Bn(1)=b~n,1eiϕb and substituting them into (6) will yield(9)
a
~
n
,
1
=
b
~
n
,
1
=
1
2
a
^
0
n
=
0,1
,
1
a
^
0
n
=
-
1,2
,
0
otherwise
.Equations (4), (7), (8), and (9) are the asymptotic expansion of the intersite solitons. One can continue the same calculation to obtain higher order corrections. Here, we limit ourselves to the first-order correction only, which is sufficient to determine the leading order behaviour of the eigenvalues later.


## 3.2. Onsite Soliton

For the onsite soliton, that is, a one-excited-site discrete mode, one can perform the same computations to obtain the mode structure of the form(10)An0=a^0eiϕan=0,0otherwise,Bn0=b^0eiϕbn=0,0otherwise,with (8). After writing An(1)=a~n,1eiϕa and Bn(1)=b~n,1eiϕb, the first order correction from (6) is given by(11)a~n,1=b~n,1=1a^0n=0,±1,0otherwise.


## 4. Stability Analysis

After we find discrete solitons, their linear stability is then determined by solving a corresponding linear eigenvalue problem. To do so, we introduce the linearisation ansatzun=(An+ϵ~(Kn+iLn)eλt)eiωt and vn=(Bn+ϵ~(Pn+iQn)eλt)eiωt,|ϵ~|≪1, and substitute this into (1) to obtain the linearised equations at O(ϵ~):(12)λKn=-An2-ωLn-ϵLn+1-2Ln+Ln-1+γKn-Qn,λLn=3An2-ωKn+ϵKn+1-2Kn+Kn-1+γLn+Pn,λPn=-RBn2+3IBn2-ωQn-ϵQn+1-2Qn+Qn-1-2RBnIBn+γPn-Ln,λQn=3RBn2+IBn2-ωPn+ϵPn+1-2Pn+Pn-1+2RBnIBn-γQn+Kn,which have to be solved for the eigenvalue λ and the corresponding eigenvector [{Kn},{Ln},{Pn},{Qn}]T. As the stability matrix of the eigenvalue problem (12) is real valued, λ¯ and -λ are also eigenvalues with corresponding eigenvectors [{Kn¯},{Ln¯},{Pn¯},{Qn¯}]T and [{Kn},{-Ln},{Pn},{-Qn}]T with γ→-γ, respectively. Therefore, we can conclude that the solution un is (linearly) stable only when R(λ)=0 for all eigenvalues λ.


### 4.1. Continuous Spectrum

The spectrum of (12) will consist of continuous spectrum and discrete spectrum (eigenvalue). To investigate the former, we consider the limit n→±∞, introduce the plane-wave ansatz Kn=k^eikn, Ln=l^eikn, Pn=p^eikn, and Qn=q^eikn, k∈R, and substitute the ansatz into (12) to obtain(13)λk^l^p^q^=γξ0-1-ξγ100-1-γξ10-ξ-γk^l^p^q^,where ξ=ω-2ϵ(cos⁡k-1). The equation can be solved analytically to yield the dispersion relation(14)λ2=4ϵωcos⁡k-1-4ϵ2cos⁡k-12-ω2-1+γ2±4ϵcos⁡k-1-2ω1-γ2.The continuous spectrum is therefore given by λ∈±[λ1-,λ2-] and λ∈±[λ1+,λ2+] with the spectrum boundaries(15)λ1±=i1-γ2+ω2±2ω1-γ2,λ2±=i1-γ2+8ϵω+16ϵ2+ω2+21-γ2±ω-4ϵ,obtained from (14) by setting k=0 and k=π in the equation.


### 4.2. Discrete Spectrum

Following the weak-coupling analysis as in Section3, we will as well use similar asymptotic expansions to solve the eigenvalue problem (12) analytically; that is, we write(16)□=□0+ϵ□1+ϵ□2+⋯,with □=λ,Kn,Ln,Pn,Qn. We then substitute the expansions into the eigenvalue problem (12).At orderO(1), one will obtain the stability equation for the dimer, which has been discussed for a general value of γ in [31]. The expression of the eigenvalues is simple, but the expression of the corresponding eigenvectors is not, which makes the result of [31] rather impractical to use. Therefore, here we limit ourselves to the case of small |γ| and expand (16) further as(17)□j=□j,0+γ□j,1+γ2□j,2+⋯,j=0,1,2,…. Hence, we have two small parameters, that is, ϵ and γ, that are independent of each other. For the sake of presentation, the detailed calculations are shown in the Appendix. Here we will only cite the final results.


#### 4.2.1. Intersite Soliton I

The intersite soliton I (i.e., the symmetric intersite soliton) has three pairs of eigenvalues for smallϵ and γ. One pair bifurcate from the zero eigenvalue. They are asymptotically given by(18)λ=ϵ2ω-1+γ22ω-1+⋯+Oϵ,λ=2ω-2-γ2ω-42ω-2+⋯+ϵω-2-γ2ω4ω-2+⋯+Oϵ3/2,2ω-2-γ2ω-42ω-2+⋯+ϵ1ω-2-γ2ω4ω-23/2+⋯+Oϵ3/2.


#### 4.2.2. Intersite Soliton II

The intersite soliton II, that is, the intersite soliton that is antisymmetric between the arms, has three pairs of eigenvalues given by(19)λ=ϵ2ω+1-γ22ω+1+⋯+Oϵ,λ=i2ω+2-γ2ω+42ω+2+⋯-iϵω+2+γ23ω4+35ω3+136ω2+208ω+1088ω+2ω3+6ω2+12ω+8+⋯+Oϵ3/2,i2ω+2-γ2ω+42ω+2+⋯+iϵ1ω+2+γ2ω4+21ω3+104ω2+184ω+1088ω+2ω3+6ω2+12ω+8+⋯+Oϵ3/2.


#### 4.2.3. Onsite Soliton I

The onsite soliton has only one eigenvalue for smallϵ given asymptotically by(20)λ=2ω-2-γ2ω-42ω-2+⋯+ϵ2ω-2-γ2ω2ω-23/2+⋯+⋯.


### 4.3. Onsite Soliton II

As for the second type of the onsite soliton, we have(21)λ=i2ω+2-γ2ω+4ω+2+⋯+2iϵ1ω+2-γ2ωω+23/2+⋯+⋯.


## 4.1. Continuous Spectrum

The spectrum of (12) will consist of continuous spectrum and discrete spectrum (eigenvalue). To investigate the former, we consider the limit n→±∞, introduce the plane-wave ansatz Kn=k^eikn, Ln=l^eikn, Pn=p^eikn, and Qn=q^eikn, k∈R, and substitute the ansatz into (12) to obtain(13)λk^l^p^q^=γξ0-1-ξγ100-1-γξ10-ξ-γk^l^p^q^,where ξ=ω-2ϵ(cos⁡k-1). The equation can be solved analytically to yield the dispersion relation(14)λ2=4ϵωcos⁡k-1-4ϵ2cos⁡k-12-ω2-1+γ2±4ϵcos⁡k-1-2ω1-γ2.The continuous spectrum is therefore given by λ∈±[λ1-,λ2-] and λ∈±[λ1+,λ2+] with the spectrum boundaries(15)λ1±=i1-γ2+ω2±2ω1-γ2,λ2±=i1-γ2+8ϵω+16ϵ2+ω2+21-γ2±ω-4ϵ,obtained from (14) by setting k=0 and k=π in the equation.


## 4.2. Discrete Spectrum

Following the weak-coupling analysis as in Section3, we will as well use similar asymptotic expansions to solve the eigenvalue problem (12) analytically; that is, we write(16)□=□0+ϵ□1+ϵ□2+⋯,with □=λ,Kn,Ln,Pn,Qn. We then substitute the expansions into the eigenvalue problem (12).At orderO(1), one will obtain the stability equation for the dimer, which has been discussed for a general value of γ in [31]. The expression of the eigenvalues is simple, but the expression of the corresponding eigenvectors is not, which makes the result of [31] rather impractical to use. Therefore, here we limit ourselves to the case of small |γ| and expand (16) further as(17)□j=□j,0+γ□j,1+γ2□j,2+⋯,j=0,1,2,…. Hence, we have two small parameters, that is, ϵ and γ, that are independent of each other. For the sake of presentation, the detailed calculations are shown in the Appendix. Here we will only cite the final results.


### 4.2.1. Intersite Soliton I

The intersite soliton I (i.e., the symmetric intersite soliton) has three pairs of eigenvalues for smallϵ and γ. One pair bifurcate from the zero eigenvalue. They are asymptotically given by(18)λ=ϵ2ω-1+γ22ω-1+⋯+Oϵ,λ=2ω-2-γ2ω-42ω-2+⋯+ϵω-2-γ2ω4ω-2+⋯+Oϵ3/2,2ω-2-γ2ω-42ω-2+⋯+ϵ1ω-2-γ2ω4ω-23/2+⋯+Oϵ3/2.


### 4.2.2. Intersite Soliton II

The intersite soliton II, that is, the intersite soliton that is antisymmetric between the arms, has three pairs of eigenvalues given by(19)λ=ϵ2ω+1-γ22ω+1+⋯+Oϵ,λ=i2ω+2-γ2ω+42ω+2+⋯-iϵω+2+γ23ω4+35ω3+136ω2+208ω+1088ω+2ω3+6ω2+12ω+8+⋯+Oϵ3/2,i2ω+2-γ2ω+42ω+2+⋯+iϵ1ω+2+γ2ω4+21ω3+104ω2+184ω+1088ω+2ω3+6ω2+12ω+8+⋯+Oϵ3/2.


### 4.2.3. Onsite Soliton I

The onsite soliton has only one eigenvalue for smallϵ given asymptotically by(20)λ=2ω-2-γ2ω-42ω-2+⋯+ϵ2ω-2-γ2ω2ω-23/2+⋯+⋯.


## 4.2.1. Intersite Soliton I

The intersite soliton I (i.e., the symmetric intersite soliton) has three pairs of eigenvalues for smallϵ and γ. One pair bifurcate from the zero eigenvalue. They are asymptotically given by(18)λ=ϵ2ω-1+γ22ω-1+⋯+Oϵ,λ=2ω-2-γ2ω-42ω-2+⋯+ϵω-2-γ2ω4ω-2+⋯+Oϵ3/2,2ω-2-γ2ω-42ω-2+⋯+ϵ1ω-2-γ2ω4ω-23/2+⋯+Oϵ3/2.


## 4.2.2. Intersite Soliton II

The intersite soliton II, that is, the intersite soliton that is antisymmetric between the arms, has three pairs of eigenvalues given by(19)λ=ϵ2ω+1-γ22ω+1+⋯+Oϵ,λ=i2ω+2-γ2ω+42ω+2+⋯-iϵω+2+γ23ω4+35ω3+136ω2+208ω+1088ω+2ω3+6ω2+12ω+8+⋯+Oϵ3/2,i2ω+2-γ2ω+42ω+2+⋯+iϵ1ω+2+γ2ω4+21ω3+104ω2+184ω+1088ω+2ω3+6ω2+12ω+8+⋯+Oϵ3/2.


## 4.2.3. Onsite Soliton I

The onsite soliton has only one eigenvalue for smallϵ given asymptotically by(20)λ=2ω-2-γ2ω-42ω-2+⋯+ϵ2ω-2-γ2ω2ω-23/2+⋯+⋯.


## 4.3. Onsite Soliton II

As for the second type of the onsite soliton, we have(21)λ=i2ω+2-γ2ω+4ω+2+⋯+2iϵ1ω+2-γ2ωω+23/2+⋯+⋯.


## 5. Numerical Results

We have solved the steady-state equation (3) numerically using a Newton-Raphson method and analysed the stability of the numerical solution by solving the eigenvalue problem (12). Here we will compare the analytical calculations obtained above with the numerical results.First, we consider the discrete intersite soliton I. We show in Figure1 the spectrum of the soliton as a function of the coupling constant ϵ for ω=1.2 and γ=0,0.5. On the real axis, one can observe that there is only one unstable eigenvalue that bifurcates from the origin. As the coupling increases, the bifurcating eigenvalue enters the origin again when ϵ→∞. Hence, in that limit we obtain a stable soliton I (i.e., a stable symmetric soliton). The dynamics of the nonzero eigenvalues as a function of the coupling constant is shown in the right panels of the figure, where one can see that the eigenvalues are on the imaginary axis and simply enter the continuous spectrum as ϵ increases.Figure 1
Eigenvalues of intersite soliton I withω=1.2, γ=0 (a, b) and 0.5 (c, d). Dots are from the numerics and solid lines are the asymptotic approximations in Section 4.2.1. The collection of dots forming black regions in the right column corresponds to the continuous spectrum.
(a) (b) (c) (d)In Figure2, we plot the eigenvalues for ω large enough. Here, in the uncoupled limit, all the three pairs of eigenvalues are on the real axis. As the coupling increases, two pairs go back toward the origin, while one pair remains on the real axis (not shown here). In the continuum limit ϵ→∞, we therefore obtain an unstable soliton I (i.e., an unstable symmetric soliton).Figure 2
The same as in Figure1 with γ=0 (a) and γ=0.5 (b), but for ω=5.2. In this case, all the eigenvalues are real.
(a) (b)In both figures, we also plot the approximate eigenvalues in solid (blue) curves, where good agreement is obtained for smallϵ.From numerical computations, we conjecture that if in the limitϵ→0 all the nonzero eigenvalues λ satisfy λ2>λ1-2 (see (15)), then we will obtain unstable soliton I in the continuum limit ϵ→∞. However, when in the anticontinuum limit ϵ→0 all the nonzero eigenvalues λ satisfy λ1+2<λ2<λ2-2, we may either obtain a stable or an unstable soliton I in the continuum limit.Next, we consider intersite solitons II (i.e., antisymmetric intersite solitons). Shown in Figure3 is the spectrum of the discrete solitons for two values of ω. In both cases, there is an eigenvalue bifurcating from the origin. For the smaller value of ω (Figures 3(a)–3(c)), we have the condition that all the nonzero eigenvalues λ satisfy λ2<λ2-2 in the anticontinuum limit ϵ→0. The collision between the eigenvalues and the continuous spectrum as the coupling increases creates complex eigenvalues. In the second case using larger ω (Figures 3(d)–3(f)), the nonzero eigenvalues λ satisfy λ2>λ1-2 when ϵ=0. Even though not seen in the figure, the collision between one of the nonzero eigenvalues and the continuous spectrum also creates a pair of complex eigenvalues. Additionally, in the continuum limit both values of ω as well as the other values of the parameter that we computed for this type of discrete solitons yield unstable solutions.Figure 3
The spectra of intersite soliton II withω=5.2 (a, b, c) and 8.2 (d, e, f) and γ=0.5. The most left panels are the spectra in the complex plane for ϵ=1. (b, c, e, f) The eigenvalues as a function the coupling constant. Solid blue curves are the asymptotic approximations.
(a) (b) (c) (d) (e) (f)We also study onsite solitons. Shown in Figures4 and 5 is the stability of discrete solitons types I and II, respectively.Figure 4
(a, c) The spectrum of onsite soliton I in the complex plane forϵ=0.2. (b, d) The eigenvalue as a function of the coupling and its approximation from Section 4.2.3. (a–d) are for ω=1.2 and 5.2, respectively. Here, γ=0.1.
(a) (b) (c) (d)Figure 5
The same as in Figure4, but for onsite soliton II with ω=5.2 (a, b) and 8.2 (c, d). (a, c) are with ϵ=1.
(a) (b) (c) (d)Figure4(a) shows that for (ω-1-γ2) small enough we will obtain stable discrete solitons. For coupling constant ϵ small, we indeed show it through our analysis depicted as the blue solid line. Numerically, we obtain that this soliton is also stable in the continuum limit ϵ→∞. However, when ω is large enough compared to 1-γ2, even though initially in the uncoupled limit the nonzero eigenvalue λ satisfies λ2<λ2-2, one may obtain an exponential instability (i.e., instability due to a real eigenvalue). Figure 4(c) shows the case when the discrete soliton is already unstable even in the uncoupled limit due to the nonzero eigenvalue that is already real-valued.Figure5 shows that the antisymmetric solitons are generally unstable due to a quartet of complex eigenvalues, as shown in Figures 5(a) and 5(c). When the coupling is increased further, there will be an eigenvalue bifurcating from ±λ1- that will move towards the origin and later becomes a pair of real eigenvalues. These solitons are also unstable in the continuum limit.Unlike intersite discrete solitons that are always unstable, onsite discrete solitons may be stable. In Figure6, we present the (in)stability region of the two types of discrete solitons in the (ϵ,ω)-plane for three values of the gain-loss parameter γ. Discrete solitons are unstable above the curves in panel (a) and between the curves in panel (b). Indeed as we mentioned before, for soliton I there is a critical ω that depends on γ below which the soliton is stable in the continuum limit, while soliton II is always unstable in that limit. Another difference between the two figures is that the stability curves in Figure 6(a) generally correspond to an eigenvalue crossing the origin that becomes real-valued, while the curves in the other panel are due to the appearance of a quartet of complex eigenvalues. In general, we obtain that the gain-loss term can be parasitic as it reduces the stability region of the discrete solitons.Figure 6
The stability region of the onsite soliton type I (a) and II (b) in the(ϵ,ω)-plane for several values of γ. The solutions are unstable above the curves in panel (a) and between the curves in panel (b).
(a) (b)Finally, we present in Figure7 the time dynamics of some of the unstable solutions shown in the previous figures. What we obtain is that typically there are two kinds of dynamics, that is, in the form of travelling discrete solitons or solution blow-ups. The first type was the typical dynamics of the intersite soliton I. The second dynamics is typical for the other types of unstable discrete solitons.Figure 7
The typical dynamics of the instability of the discrete solitons in the previous figures. Here,γ=0.5 and ϵ=1. Depicted in (a–h) are |un|2 and |vn|2, respectively. (a–h) The dynamics of intersite soliton I with ω=1.2, intersite soliton II with ω=5.2, and onsite solitons I and II both with ω=5.2.
(a) (b) (c) (d) (e) (f) (g) (h)


## 6. Conclusion

We have presented a systematic method to determine the stability of discrete solitons in aPT-symmetric coupler by computing the eigenvalues of the corresponding linear eigenvalue problem using asymptotic expansions. We have compared the analytical results that we obtained with numerical computations, where good agreement is obtained. From the numerics, we have also established the mechanism of instability as well as the stability region of the discrete solitons. The application of the method in higher dimensional PT-symmetric couplers (see, e.g., [32]) is a natural extension of the problem that is addressed for future work. Additionally, we also address the computation of eigenvalues of discrete solitons in the neighbourhood of broken PT-symmetry as future investigations.



---

*Source: 9514230-2016-12-04.xml*