Acta Polytechnica doi:10.14311/AP.2014.54.0133 Acta Polytechnica 54(2):133–138, 2014 © Czech Technical University in Prague, 2014 available online at http://ojs.cvut.cz/ojs/index.php/ap STABILITY OF BOSE-EINSTEIN CONDENSATES IN A PT-SYMMETRIC DOUBLE-δ POTENTIAL CLOSE TO BRANCH POINTS Andreas Löhle, Holger Cartarius∗, Daniel Haag, Dennis Dast, Jörg Main, Günter Wunner Institut für Theoretische Physik 1, Universität Stuttgart, 70550 Stuttgart, Germany ∗ corresponding author: Holger.Cartarius@itp1.uni-stuttgart.de Abstract. A Bose-Einstein condensate trapped in a double-well potential, where atoms are incoupled to one side and extracted from the other, can in the mean-field limit be described by the nonlinear Gross-Pitaevskii equation (GPE) with a PT symmetric external potential. If the strength of the in- and outcoupling is increased two PT broken states bifurcate from the PT symmetric ground state. At this bifurcation point a stability change of the ground state is expected. However, it is observed that this stability change does not occur exactly at the bifurcation but at a slightly different strength of the in-/outcoupling effect. We investigate a Bose-Einstein condensate in a PT symmetric double-δ potential and calculate the stationary states. The ground state’s stability is analysed by means of the Bogoliubov-de Gennes equations and it is shown that the difference in the strength of the in-/outcoupling between the bifurcation and the stability change can be completely explained by the norm-dependency of the nonlinear term in the Gross-Pitaevskii equation. Keywords: Bose-Einstein condensates, PT symmetry, stability, Bogoliubov-de Gennes equations. 1. Introduction In quantum mechanics non-Hermitian Hamiltonians with imaginary potentials have become an impor- tant tool to describe systems with loss or gain effects [1]. Non-Hermitian PT symmetric Hamiltonians, i.e. Hamiltonians commuting with the combined action of the parity (P: x →−x, p →−p) and time reversal (T : x → x, p → −p, i → −i) operators, possess the interesting property that, in spite of the gain and loss, they can exhibit stationary states with real eigen- values [2]. When the strength of the gain and loss contributions is increased typically pairs of these real eigenvalues pass through an exceptional point, i.e. a branch point at which both the eigenvalues and the wave functions are identical, and become complex and complex conjugate. A promising candidate for the realisation of a PT symmetric quantum system are Bose-Einstein conden- sates. At sufficiently low temperatures and densities they can in the mean-field limit be described by the nonlinear Gross-Pitaevskii equation [3, 4]. If a conden- sate is trapped in a double-well potential it is possible to add atoms to one side of the double well and re- move atoms from the other. This leads to a gain or loss to the condensate’s probability amplitude. If the strength of both contributions is equal, the process can effectively be described by a complex external potential V (−x) = V ∗(x) rendering the Hamiltonian PT symmetric [5]. The experimental realisation of an open quantum system with PT symmetry will be an important step since the experimental verification has only been achieved in optics so far [6, 7]. The nonlinearity ∝ |ψ(x, t)|2 of the Gross-Pitaevskii equation introduces new interesting properties [8–17] such as PT broken states which appear for gain/loss contributions lower than those at which the corre- sponding PT symmetric states vanish. In optics the same effect can be observed for wave guides with a Kerr nonlinearity. They are described by an equa- tion that is mathematically equivalent to the Gross- Pitaevskii equation. It has been shown that the addi- tional features appearing in the presence of the non- linearity might be exploited to create uni-directional structures [18] or solitons [19–24]. Since PT symme- try in optics is extensively studied [5, 18–28] and has been experimentally realised [6, 7] these approaches seem to be very promising. From the mathematical point of view a new type of bifurcation appearing in the nonlinear PT sym- metric Gross-Pitaevskii equation is of special interest [8–17]. As in linear quantum systems, two PT sym- metric stationary states merge in an exceptional point if the strength of a parameter describing the gain and loss processes is increased. However, in contrast to its linear counterpart, the Gross-Pitaevskii equation possesses no PT broken states emerging at this excep- tional point. They already appear for lower strengths of the gain/loss parameter and bifurcate from one of the PT symmetric states in a pitchfork bifurcation. 133 http://dx.doi.org/10.14311/AP.2014.54.0133 http://ojs.cvut.cz/ojs/index.php/ap A. Löhle, H. Cartarius, D. Haag et al. Acta Polytechnica The new bifurcation point has been identified to be a third-order exceptional point [10]. For attractive nonlinearities one finds that the PT broken solutions bifurcate from the ground state. In this scenario the PT symmetric ground state is the only state which exists on both sides of the bifurcation and always possesses a real energy eigenvalue. The pitchfork bi- furcation is expected to entail a change of its stability. However, it is observed that this stability change does not occur exactly at the bifurcation but at a slightly different value of the gain/loss parameter. The discrepancy between the points of bifurcation and stability change seems to be surprising and does not appear in all similar systems. The mean-field limit of a two mode approximation with a Bose-Hubbard Hamiltonian [14–17] does not show this effect. The model has, however, two crucial differences from the treatment of Bose-Einstein condensates with the full Gross-Pitaevskii equation in Ref. [11]. The latter system contains a harmonic trap in which infinitely many stationary states can be found, whereas the nonlinear two-mode system exhibits only four states, viz. the two PT symmetric and the two PT broken states mentioned above. Furthermore the nonlinearity derived in [14, 15] is slightly different from that of the Gross-Pitaevskii equation. It has the form ∝ |ψ|2/||ψ||2, and hence does not depend on the norm of the wave function. Thus, there might be two reasons for the appearance of the discrepancy. It could have its origin in the existence of higher modes influencing the ground state’s stability or in the norm-dependency of the Gross-Pitaevskii nonlinearity. It is the purpose of this article to clarify this ques- tion. To do so we study a Bose-Einstein condensate in an idealised double-δ trap [8–10], a system of which al- ready its linear counterpart helped to understand basic properties of PT symmetric structures [29–33]. This system is described by the Gross-Pitaevskii equation, i.e. the contact interaction has the norm-dependent form ∝ |ψ(x, t)|2. However, it exhibits only four sta- tionary states of which two are PT symmetric and two are PT broken as in the two-mode model [14– 17]. Additionally, in a numerical study the structure of the nonlinearity can easily be changed such that the system’s mathematical properties can be brought into agreement with the mean-field limit of the Bose- Hubbard dimer. The article is organised as follows. We will intro- duce and solve the Gross-Pitaevskii equation of a Bose-Einstein condensate in a double-δ trap for an attractive atom-atom interaction in Section 2. Some properties of the stationary solutions which are im- portant for the following discussions are recapitulated. Then we will investigate the ground state’s stability in the vicinity of the bifurcation in Section 3. The Bogoliubov-de Gennes equations are solved for both types of the nonlinearity, and the origin of the discrep- ancy between the bifurcation and the stability change is discussed. Conclusions are drawn in Section 4. 2. Bose-Einstein condensates in the PT symmetric double-δ trap We assume the condensate to be trapped in an ide- alised trap of two delta functions, i.e. the potential has the shape [8] V (x) = −(1 − iγ) δ (x− b) − (1 + iγ) δ (x + b) , (1) where the units are chosen such that the real part due to the action of the δ functions has the value -1. It describes two symmetric infinitely thin potential wells at positions ±b. In the left well we describe an outflux of atoms by a negative imaginary contribution −iγ, and on the right side an influx of particles is described by a positive imaginary part +iγ of the same strength. This leads to the time-independent Gross-Pitaevskii equation in dimensionless form [8], [ − d2 dx2 − (1 − iγ) δ (x− b) − (1 + iγ) δ (x + b) + g ∣∣ψ(x)∣∣2]ψ(x) = −κ2ψ(x) (2) with the energy eigenvalue or chemical potential µ = −κ2. The parameter g is determined by the s-wave scattering length, which effectively describes the van der Waals interaction for low temperatures and densities. Physically it can be tuned via Feshbach resonances. Throughout this article we assume g to be negative, i.e. the atom-atom interaction is attractive. Solutions to the Gross-Pitaevskii equation are found with a numerical exact integration. The wave func- tion is integrated outward from x = 0 in positive and negative direction. To do so the initial values ψ(0),ψ′(0), and κ have to be chosen. The arbitrary global phase is exploited such that ψ(0) is chosen to be real. Together with ψ′(0) ∈ C and κ ∈ C we have five parameters which have to be set such that physically relevant wave functions are obtained. These have to be square-integrable and normalised in the nonlinear system. The five conditions ψ(∞) → 0, ψ(−∞) → 0, and ||ψ|| = 1 ensure that these conditions are fulfilled, and, together with the five initial parameters, define a five-dimensional root search problem, which is solved numerically. Figure 1 shows a typical example for all station- ary states found in the case g = −1 (red solid lines). Two states with purely real eigenvalues vanish for increasing γ in an exceptional point. Since κ is plot- ted and µ = −κ2 is the energy eigenvalue, the upper line corresponds to the ground state. Two complex and complex conjugate eigenvalues bifurcate from this ground state in a pitchfork bifurcation at a critical value γκ ≈ 0.3071. The same plots for g = −0.5 (blue dotted lines) and g = 0 (green dashed lines) demon- strate how the pitchfork bifurcation is introduced by the nonlinearity. 134 vol. 54 no. 2/2014 Stability of BECs in a PT-Symmetric Double-δ Potential 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.1 0.2 0.3 0.4 0.5 0.6 R e (κ ) γ (a) γ κ g = -1.0 g = -0.5 g = 0 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0 0.1 0.2 0.3 0.4 0.5 0.6 Im (κ ) γ (b) γ κ Figure 1. Real (a) and imaginary (b) parts of the stationary solutions found for the Gross-Pitaevskii equation (2) for the case g = −1 (red solid lines), which is compared with its linear counterpart g = 0 (green dashed lines) and the weaker interaction g = −0.5 (blue dotted lines). In the nonlinear case we observe two complex conjugate states bifurcating from the ground state in a pitchfork bifurcation. For g = −1 this occurs at γκ ≈ 0.3071. This eigenvalue structure is very generic for PT symmetric systems with a quadratic term in the Hamil- tonian. It is found for Bose-Einstein condensates in a true spatially extended double well in one and three dimensions [11, 12] which contain the nonlinearity ∝ |ψ|2 as well as for the mean-field limit of the Bose- Hubbard dimer [14–17] with a term ∝ |ψ|2/||ψ||2. The difference in the norm-dependency of the nonlinearity between the two systems does not lead to different eigenvalues, as has already been mentioned [17]. The difference appears, however, for the dynamics, where the norm plays a crucial role in the presence of gain and loss [13]. 3. Stability analysis of the ground state The linear stability is analysed with the Bogoliubov-de Gennes equations. They are derived under the assump- tion that a stationary state ψ0(x,t) is perturbed by a small fluctuation θ(x,t), i.e. ψ(x,t) = eiκ 2t [ ψ0(x) + θ(x,t) ] , (3a) where θ(x,t) = u(x)e−iωt + v∗(x)eiω ∗t. (3b) With this ansatz and a linearisation in the small quan- tities u and v one obtains from the Gross-Pitaevskii equation the coupled system of the Bogoliubov-de Gennes differential equations, d2 dx2 u(x) = [ −(1 + iγ)δ(x + b) − (1 − iγ)δ(x− b) + κ2 −ω + 2g ∣∣ψ0(x)∣∣2]u(x) + gψ0(x)2v(x), (4a) d2 dx2 v(x) = [ −(1 − iγ)δ(x + b) − (1 + iγ)δ(x− b) + (κ2)∗ + ω + 2g ∣∣ψ0(x)∣∣2]v(x) + gψ∗0 (x)2u(x). (4b) In Equation (3b) it can be seen that ω decides on the temporal evolution of the fluctuation. Real values of ω describe stable oscillations, whereas imaginary parts lead to a growth or decay of the fluctuation’s amplitude. Thus, ω measures the stability of the stationary solution ψ0 against small fluctuations. Numerically the Bogoliubov-de Gennes equations are solved with the same method as the station- ary states, i.e. the wave functions u and v are in- tegrated outward from x = 0. It can easily be seen that the Bogoliubov-de Gennes equations (4) are in- variant under the transformation u(x) → u(x)eiχ, v(x) → v(x)eiχ with a real phase χ. Similarly to the procedure for the integration of the stationary states this symmetry can be exploited. The remaining initial values with which the integration has to be started are u(0) ∈ C, Re(v(0)), u′(0) ∈ C, v′(0) ∈ C, ω ∈ C. In a nine-dimensional root search they have to be chosen such that the conditions u(±∞) → 0, v(±∞) → 0, and ∫ ∞ −∞ ∣∣u(x) + v∗(x)∣∣2 dx = 1 are fulfilled [11]. Two further symmetries can be exploited to re- duce the number of independent stability eigenval- ues that have to be calculated. The replacement (u,v,ω) → (v∗,u∗,−ω∗) leaves the ansatz (3b) invari- ant. Thus, if ω is a stability eigenvalue also −ω∗ is a valid solution. Furthermore for every eigenvalue ω there is also one solution with the eigenvalue ω∗ if the stationary state ψ0 is PT symmetric. This can be verified by applying the PT operator to the Bogoliubov-de Gennes equations (4). Due to these two symmetries it is sufficient to search for stability eigenvalues with Re(ω) > 0 and Im(ω) > 0. 3.1. Stability in the vicinity of the bifurcation The relevant question which has to be answered by our calculation is whether or not the discrepancy between the γ values of the pitchfork bifurcation and the sta- bility change appears for the double-δ potential. Thus we calculated the stability eigenvalue with Re(ω) > 0, Im(ω) > 0 for a range of γ around the bifurcation, which is shown in Figure 2 for g = −1. For increasing 135 A. Löhle, H. Cartarius, D. Haag et al. Acta Polytechnica 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 κ , ω γ γ κ Re(ω) Im(ω) Re(κ) Figure 2. Real (red solid line) and imaginary (green dashed line) part of the stability eigenvalue ω for the stationary ground state in the case g = −1. The stability change occurs at γω ≈ 0.3138. To illustrate the pitchfork bifurcation the real parts of κ of all stationary states are also shown (blue dotted lines). The value γκ ≈ 0.307 is marked by the black dashed- dotted line. Obviously there is a discrepancy between the two values. γ the eigenvalue ω switches from real to imaginary at γω ≈ 0.3138 marking the stability change. The pitchfork bifurcation is visible in the real parts of κ of all stationary states of the system. It is marked by the black dashed-dotted line. Obviously there is a discrepancy between γκ and γω. The difference ∆γ = γκ −γω (5) is ∆γ ≈−0.0067. The system does not possess any further stationary states besides those shown in Figure 1. Three of these four states are participating in the pitchfork bifurcation. Only the excited PT symmetric solution would be able to influence the dynamics of the ground state at a value γ 6= γω. However, it stays real for all parameters γ shown in Figure 2 and cannot cause any qualitatively different behaviour of the ground state’s dynamics. Thus, an influence of further states can be ruled out as the reason for the discrepancy in the double-well system of Refs. [11, 12]. The remaining difference between the Gross- Pitaevskii equation (2) and the two mode system of Refs. [14–17] is the norm-dependency of the non- linearity. Indeed, an influence of the norm is already present if the study done in Figure 2 is repeated for dif- ferent values of the nonlinearity parameter g. Figure 3 shows ∆γ as a function of g. A strong dependency is visible. Even the sign changes. For g / 0.4 the ground state becomes unstable at γω < γκ. For g → 0 the discrepancy vanishes as expected. 3.2. Norm-independent variant of the Gross-Pitaevskii equation An even clearer identification of ∆γ with the norm- dependency of the nonlinearity in the Gross-Pitaevskii -0.007 -0.006 -0.005 -0.004 -0.003 -0.002 -0.001 0 0.001 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 ∆ γ -g Figure 3. Difference between γκ and γω as defined in Equation (5) as a function of g. A strong dependency is clearly visible. equation (2) can be given with a small modification. The replacement g|ψ|2 → g|ψ|2∫ |ψ|2 dx (6) makes the Gross-Pitaevskii equation (2) norm- independent. Note that this is exactly the form of the mean-field limit of Refs. [14–17]. Since the stationary states are normalised to 1 they are not influenced by the replacement. However, it influences the dynam- ics and also the linear stability. The Bogoliubov-de Gennes equations have to be adapted. Assuming again a small perturbation of the form (3a) and (3b) a linearisation in u and v leads us to d2 dx2 u(x) = [ −(1 + iγ)δ(x + b) − (1 − iγ)δ(x− b) + κ2 −ω + 2g ∣∣ψ0(x)∣∣2]u(x) + gψ0(x)2v(x) −g ∣∣ψ0(x)∣∣2ψ0(x)S, (7a) d2 dx2 v(x) = [ −(1 − iγ)δ(x + b) − (1 + iγ)δ(x− b) + (κ2)∗ + ω + 2g ∣∣ψ0(x)∣∣2]v(x) + gψ∗0 (x)2u(x) −g ∣∣ψ0(x)∣∣2ψ∗0 (x)S (7b) with the integral S = ∫ [ v(x)ψ0(x) + u(x)ψ∗0 (x) ] dx. (7c) For a numerical solution of the modified Bogoliubov- de Gennes equations the value of S is included in the root search, i.e. for the integration of Equations (7a) and (7c) a value for S is guessed and subsequently compared with the result of Equation (7c). Since S is in general a complex value this increases the dimension of the root search to 11. An example for a typical result is shown in Figure 4. The discrepancy between the γ values at which the pitchfork bifurcation and the stability change occur vanishes. The two values agree, as is expected for a pitchfork bifurcation. This is true for all values of g. 136 vol. 54 no. 2/2014 Stability of BECs in a PT-Symmetric Double-δ Potential 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 κ , ω γ γ κ = γ ω Re(ω) Im(ω) Re(κ) Figure 4. The same as in Figure 2 but for the modi- fied Bogoliubov-de Gennes equations (7a)-(7c). Now the γ values at which the pitchfork bifurcation and the stability change occur are in perfect agreement. Thus, we conclude that the discrepancy appearing in Figure 2 for the Gross-Pitaevskii equation is solely a result of the norm-dependent nonlinearity ∝ |ψ|2 in the Hamiltonian and is not a consequence of the interaction with higher excited states. 3.3. Discussion The reason why the stability change does not occur exactly at the bifurcation can also be understood in- tuitively. Figure 1 indicates that the values γκ of the pitchfork bifurcation are not equal for all values of g. For g = 0 there is no pitchfork bifurcation. The two real eigenvalues κ vanish in a tangent bifurcation and two complex eigenvalues emerge. Only for nonvanish- ing g the pitchfork bifurcation exists and moves to smaller values of γ for increasing |g|. If now the norm of a wave function changes, this can also be understood as a variation of g, i.e. a wave function with a norm N = ||ψ||2 has the same effect as a wave function with the norm 1 and a modified value g̃ = Ng. The values γκ of the pitchfork bifurcation are obviously different for g and g̃. With this relation in mind it is not surprising that a fluctuation changing the norm of the wave function may cause a qualitative change of the condensate’s stability properties in the vicinity of γκ. 4. Conclusions We investigated the origin of the discrepancy between the value of the gain/loss parameters at which the ground state of a Bose-Einstein condensate in a double- well trap passes through a pitchfork bifurcation and at which its stability changes. In a naive expecta- tion these γ values should be identical. However, it is found that this is not exactly fulfilled. Since this discrepancy does not occur for a similar system, viz. the mean-field limit of a Bose-Hubbard dimer [14–17], we investigated the differences between the equations describing the two systems. It was found that the norm-dependency of the nonlinearity in the Hamiltonian ∝ |ψ|2 is unambiguously the origin of the discrepancy. It can be completely removed with the replacement |ψ|2 → |ψ|2/||ψ||2. An intuitive ex- planation can be given. Fluctuations which change the norm of a stationary state are able to shift the position γκ of the bifurcation. Since the dynamical properties of the wave func- tions are crucial for the experimental observability of a PT symmetric Bose-Einstein condensate it is important to know about all processes introducing possible instabilities. As has been shown in this article the stability relations are nontrivial close to branch points in condensate setups with gain and loss. The Gross-Pitaevskii has a norm-dependent nonlinearity, and therefore it should be clarified in future work which type of fluctuation influences the wave func- tion’s norm such that additional instabilities appear. In particular, it would be interesting to see, how the amplitude of a fluctuation is related to the size of the difference ∆γ. A deeper understanding of how this effect influences realistic setups generating the PT symmetric external potential [34], would also be of high value. Acknowledgements We thank Eva-Maria Graefe for stimulating discussions. References [1] N. Moiseyev. Non-Hermitian Quantum Mechanics. Cambridge University Press, Cambridge, 2011. [2] C. M. Bender, et al. Real spectra in non-Hermitian Hamiltonians having PT symmetry. Phys Rev Lett 80:5243, 1998. doi: 10.1103/PhysRevLett.80.5243 [3] E. P. Gross. Structure of a quantized vortex in Boson systems. Nuovo Cimento 20:454, 1961. doi: 10.1007/BF02731494 [4] L. P. Pitaevskii. Vortex lines in an imperfect Bose gas. Sov Phys JETP 13:451, 1961. [5] S. Klaiman, et al. Visualization of branch points in PT -symmetric waveguides. Phys Rev Lett 101:080402, 2008. doi: 10.1103/PhysRevLett.101.080402 [6] A. Guo, et al. Observation of PT -symmetry breaking in complex optical potentials. Phys Rev Lett 103:093902, 2009. doi: 10.1103/PhysRevLett.103.093902 [7] C. E. Rüter, et al. Observation of parity-time symmetry in optics. Nat Phys 6:192, 2010. doi: 10.1038/nphys1515 [8] H. Cartarius, et al. Model of a PT -symmetric Bose-Einstein condensate in a δ-function double-well potential. Phys Rev A 86:013612, 2012. doi: 10.1103/PhysRevA.86.013612 [9] H. Cartarius, et al. Nonlinear Schrödinger equation for a PT -symmetric delta-function double well. J Phys A 45:444008, 2012. doi: 10.1088/1751-8113/45/44/444008 [10] W. D. Heiss, et al. Spectral singularities in PT -symmetric Bose-Einstein condensates. J Phys A 46:275307, 2013. doi: 10.1088/1751-8113/46/27/275307 137 http://dx.doi.org/10.1103/PhysRevLett.80.5243 http://dx.doi.org/10.1007/BF02731494 http://dx.doi.org/10.1103/PhysRevLett.101.080402 http://dx.doi.org/10.1103/PhysRevLett.103.093902 http://dx.doi.org/10.1038/nphys1515 http://dx.doi.org/10.1103/PhysRevA.86.013612 http://dx.doi.org/10.1088/1751-8113/45/44/444008 http://dx.doi.org/10.1088/1751-8113/46/27/275307 A. Löhle, H. Cartarius, D. Haag et al. Acta Polytechnica [11] D. Dast, et al. A Bose-Einstein condensate in a PT symmetric double well. Fortschr Physik 61:124–139, 2013. doi: 10.1002/prop.201200080 [12] D. Dast, et al. Eigenvalue structure of a Bose-Einstein condensate in a PT -symmetric double well. J Phys A 46:375301, 2013. doi: 10.1088/1751-8113/46/37/375301 [13] H. Cartarius, et al. Stationary and dynamical solutions of the Gross-Pitaevskii equation for a Bose-Einstein condensate in a PT symmetric double well. Acta Polytechnica 53(3):259–267, 2013. [14] E. M. Graefe, et al. A non-Hermitian PT symmetric Bose-Hubbard model: eigenvalue rings from unfolding higher-order exceptional points. J Phys A 41:255206, 2008. doi: 10.1088/1751-8113/41/25/255206 [15] E. M. Graefe, et al. Mean-field dynamics of a non-Hermitian Bose-Hubbard dimer. Phys Rev Lett 101:150408, 2008. doi: 10.1103/PhysRevLett.101.150408 [16] E. M. Graefe, et al. Quantum-classical correspondence for a non-Hermitian Bose-Hubbard dimer. Phys Rev A 82:013629, 2010. doi: 10.1103/PhysRevA.82.013629 [17] E. M. Graefe. Stationary states of a PT symmetric two-mode Bose-Einstein condensate. J Phys A 45:444015, 2012. doi: 10.1088/1751-8113/45/44/444015 [18] H. Ramezani, et al. Unidirectional nonlinear PT -symmetric optical structures. Phys Rev A 82:043803, 2010. doi: 10.1103/PhysRevA.82.043803 [19] Z. H. Musslimani, et al. Optical solitons in PT periodic potentials. Phys Rev Lett 100:30402, 2008. doi: 10.1103/PhysRevLett.100.030402 [20] F. K. Abdullaev, et al. Dissipative periodic waves, solitons, and breathers of the nonlinear Schrödinger equation with complex potentials. Phys Rev E 82:056606, 2010. doi: 10.1103/PhysRevE.82.056606 [21] Y. V. Bludov, et al. Nonlinear patterns in Bose-Einstein condensates in dissipative optical lattices. Phys Rev A 81:013625, 2010. doi: 10.1103/PhysRevA.81.013625 [22] R. Driben, et al. Stability of solitons in parity-time-symmetric couplers. Opt Lett 36:4323–4325, 2011. doi: 10.1364/ol.36.004323 [23] F. K. Abdullaev, et al. Solitons in PT -symmetric nonlinear lattices. Phys Rev A 83:41805, 2011. doi: 10.1103/PhysRevA.83.041805 [24] Y. V. Bludov, et al. Stable dark solitons in PT -symmetric dual-core waveguides. Phys Rev A 87:013816, 2013. doi: 10.1103/physreva.87.013816 [25] K. G. Makris, et al. Beam dynamics in PT symmetric optical lattices. Phys Rev Lett 100:103904, 2008. doi: 10.1103/PhysRevLett.100.103904 [26] K. G. Makris, et al. PT -symmetric optical lattices. Phys Rev A 81:063807, 2010. doi: 10.1103/physreva.81.063807 [27] A. Ruschhaupt, et al. Physical realization of PT -symmetric potential scattering in a planar slab waveguide. J Phys A 38:L171, 2005. doi: 10.1088/0305-4470/38/9/L03 [28] R. El-Ganainy, et al. Theory of coupled optical PT -symmetric structures. Opt Lett 32:2632, 2007. doi: 10.1364/ol.32.002632 [29] V. Jakubský, et al. An explicitly solvable model of the spontaneous PT -symmetry breaking. Czech J Phys 55:1113–1116, 2005. doi: 10.1007/s10582-005-0115-x [30] A. Mostafazadeh. Delta-function potential with a complex coupling. J Phys A 39:13495, 2006. doi: 10.1088/0305-4470/39/43/008 [31] A. Mostafazadeh, et al. Spectral singularities, biorthonormal systems and a two-parameter family of complex point interactions. J Phys A 42:125303, 2009. doi: 10.1088/1751-8113/42/12/125303 [32] H. Mehri-Dehnavi, et al. Application of pseudo-Hermitian quantum mechanics to a complex scattering potential with point interactions. J Phys A 43:145301, 2010. doi: 10.1088/1751-8113/43/14/145301 [33] H. F. Jones. Interface between Hermitian and non- Hermitian Hamiltonians in a model calculation. Phys Rev D 78:065032, 2008. doi: 10.1103/physrevd.78.065032 [34] M. Kreibich, et al. Hermitian four-well potential as a realization of a PT -symmetric system. Phys Rev A 87:051601(R), 2013. doi: 10.1103/physreva.87.051601 138 http://dx.doi.org/10.1002/prop.201200080 http://dx.doi.org/10.1088/1751-8113/46/37/375301 http://dx.doi.org/10.1088/1751-8113/41/25/255206 http://dx.doi.org/10.1103/PhysRevLett.101.150408 http://dx.doi.org/10.1103/PhysRevA.82.013629 http://dx.doi.org/10.1088/1751-8113/45/44/444015 http://dx.doi.org/10.1103/PhysRevA.82.043803 http://dx.doi.org/10.1103/PhysRevLett.100.030402 http://dx.doi.org/10.1103/PhysRevE.82.056606 http://dx.doi.org/10.1103/PhysRevA.81.013625 http://dx.doi.org/10.1364/ol.36.004323 http://dx.doi.org/10.1103/PhysRevA.83.041805 http://dx.doi.org/10.1103/physreva.87.013816 http://dx.doi.org/10.1103/PhysRevLett.100.103904 http://dx.doi.org/10.1103/physreva.81.063807 http://dx.doi.org/10.1088/0305-4470/38/9/L03 http://dx.doi.org/10.1364/ol.32.002632 http://dx.doi.org/10.1007/s10582-005-0115-x http://dx.doi.org/10.1088/0305-4470/39/43/008 http://dx.doi.org/10.1088/1751-8113/42/12/125303 http://dx.doi.org/10.1088/1751-8113/43/14/145301 http://dx.doi.org/10.1103/physrevd.78.065032 http://dx.doi.org/10.1103/physreva.87.051601 Acta Polytechnica 54(2):133–138, 2014 1 Introduction 2 Bose-Einstein condensates in the PT symmetric double- trap 3 Stability analysis of the ground state 3.1 Stability in the vicinity of the bifurcation 3.2 Norm-independent variant of the Gross-Pitaevskii equation 3.3 Discussion 4 Conclusions Acknowledgements References