Acta Polytechnica doi:10.14311/AP.2014.54.0116 Acta Polytechnica 54(2):116–121, 2014 © Czech Technical University in Prague, 2014 available online at http://ojs.cvut.cz/ojs/index.php/ap A BOSE-EINSTEIN CONDENSATE WITH PT-SYMMETRIC DOUBLE-DELTA FUNCTION LOSS AND GAIN IN A HARMONIC TRAP: A TEST OF RIGOROUS ESTIMATES Daniel Haag, Holger Cartarius, Günter Wunner∗ Institut für Theoretische Physik 1, Universität Stuttgart, 70550 Stuttgart, Germany ∗ corresponding author: wunner@itp1.uni-stuttgart.de Abstract. We consider the linear and nonlinear Schrödinger equation for a Bose-Einstein condensate in a harmonic trap with PT -symmetric double-delta function loss and gain terms. We verify that the conditions for the applicability of a recent proposition by Mityagin and Siegl on singular perturbations of harmonic oscillator type self-adjoint operators are fulfilled. In both the linear and nonlinear case we calculate numerically the shifts of the unperturbed levels with quantum numbers n of up to 89 in dependence on the strength of the non-Hermiticity and compare with rigorous estimates derived by those authors. We confirm that the predicted 1/n1/2 estimate provides a valid upper bound on the shrink rate of the numerical eigenvalues. Moreover, we find that a more recent estimate of log(n)/n3/2 is in excellent agreement with the numerical results. With nonlinearity the shrink rates are found to be smaller than without nonlinearity, and the rigorous estimates, derived only for the linear case, are no longer applicable. Keywords: PT symmetry, Bose-Einstein condensates, perturbed harmonic oscillator. 1. Introduction Bose-Einstein condensates with PT -symmetric loss and gain have been proposed [1] as a first experimental realisation of PT symmetry in a real quantum system. The idea is to confine the condensate in a double well potential, and to create a PT -symmetric situation by coherently injecting atoms into one well and removing them from the other. A particular difficulty in the theoretical treatment is that, on account of the s- wave scattering of the atoms, the Gross-Pitaevskii equation, which describes the condensates, contains a term g|ψ|2, and is -shell-escapebthus nonlinear in the wanted wave function. In a series of papers both for realistic set-ups as well as for delta-function models of the double wells we have shown [2–7] that the nonlinearity introduces new features in the evolution of the eigenvalue spectrum as the non-Hermiticity is increased but yet PT symmetry of the wave function is preserved if both nonlinearity and non-Hermiticity are not too strong. Bose-Einstein condensates are usually trapped in harmonic potentials produced by counterpropagating laser beams. Therefore condensates with an addi- tional PT -symmetric double well potential can be regarded as a “perturbation” of the harmonic oscilla- tor. Aducci and Mityagin [8] and Siegl and Mityagin [9] have recently analysed perturbations of harmonic- like operators from a mathematical point of view. The second paper in particular also allows for singular perturbations, such as delta-functions. These authors have proved that the eigenvalues of the perturbed operator eventually become simple and the root system forms a Riesz basis. Their results are valid if the following criterion is fulfilled: ∀m,n ∈ N ∣∣〈ψm|B̂|ψn〉∣∣ ≤ M mαnα , α > 0, M > 0, (1) where B̂ is the perturbation operator, M is a constant that depends on the type of the perturbation, and the ψm are the harmonic oscillator eigenfunctions. A further prediction is that the shrink rate of the disks into which the eigenvalues can be shifted is proportional to 1/n2α, with α the exponent appearing in (1). It is the purpose of this paper to test the esti- mates numerically for the example of a Bose-Einstein condensate in a double well potential confined by a harmonic trap. To this end we consider the model of two PT -symmetric delta-function wells, since in this case simple analytical estimates can be obtained, whereas in the case of the realistic double well dis- cussed in Refs. [4, 5] complicated estimates in terms of hypergeometric functions result. It should be mentioned that of course because of their simplicity delta functions have been widely used in the literature in the context of PT symmetry. Spec- tral properties of scattering and bound states in PT - symmetric double- and multiple-delta function poten- tials have been investigated e.g. in Refs. [2, 3, 11–18]. In all these papers no external potential was present, in addition to the delta potentials. A paper in which PT -symmetric point interactions were studied embedded in an external potential is 116 http://dx.doi.org/10.14311/AP.2014.54.0116 http://ojs.cvut.cz/ojs/index.php/ap vol. 54 no. 2/2014 BEC with PT-Symmetric Double-Delta Loss and Gain: Test of Estimates that by Jakubský and Znojil [19], who positioned the delta functions in an infinitely high square well and analysed the spectrum in dependence on the position of the delta functions within the well. Their work was extended by Krejčiřík and Siegl [20, 21] who replaced the delta functions by PT -symmetric Robin boundary conditions at the edges of the square well. We note that the spectrum of a harmonic oscillator perturbed by two identical real-valued point interactions has been analysed by Fassari and Rinaldi [22]. However, to our knowledge the situation of two PT -symmetric delta functions in an external harmonic potential has not yet been investigated. While in the square well the delta functions can be placed only within the well, the harmonic oscillator potential has the advantage that the delta functions can in principle be shifted to any position on the real axis. 2. Bose-Einstein condensate in a PT -symmetric harmonic trap At low temperatures and densities Bose-Einstein con- densates are well described by the Gross-Pitaevskii equation [23, 24] ( − d2 dx2 + V (x) + g|ψ|2 ) ψ = µψ. (2) Here ψ denotes the condensate wave function, the eigenvalue µ is the chemical potential, and V (x) is the trapping potential to confine the condensate. The nonlinear term in (2) arises from the s-wave scatter- ing interaction of the atoms; g is a measure for the strength of this interaction. We consider a harmonic trapping potential and model a PT -symmetric double well with equilibrated loss and gain by imaginary delta functions. Thus the Hamiltonian that we consider here is given by Ĥ = − d2 dx2 +x2+iγ ( δ(x−b)−δ(x+b) ) +g ∣∣ψ(x)∣∣2. (3) Here ±b denotes the position of the imaginary deltas and γ the strength of the non-Hermiticity. We will later consider the effects of the nonlinearity on the spectrum, but for the time being we assume that the nonlinearity is negligible in order to be in a position to compare with the predictions of Mityagin and Siegl [9]. The eigenvalues of the unperturbed spectrum are given by µn = 2n+ 1 (n = 0, 1, 2, . . . ). Figure 1 shows the unperturbed spectrum together with the wave functions of the lowest five states. The dashed vertical lines designate different positions at which the delta functions are placed. From the figure it is already obvious that only such states will be significantly affected by the perturbation which are within the classically allowed region at the positions of the delta functions. By contrast, states for which the delta functions lie in the classically forbidden (exponentially decaying) regime will not be affected. This means that 0 2 4 6 8 10 −3 −2 −1 0 1 2 3 -0.5 0 0.5 -0.5 0 0.5 -0.5 0 0.5 -0.5 0 0.5 -0.5 0 0.5 V ψ x Figure 1. Unperturbed spectrum of the harmonic oscillator and the lowest five wave functions. Vertical dashed lines designate different positions of the PT - symmetric delta-function perturbations. as the delta functions are shifted further and further out an increasing number of low-lying eigenvalues will not be changed by the perturbation. It is easy to show that the Mityagin-Siegl crite- rion (1) is fulfilled for the imaginary delta function perturbation in (3) since∣∣〈ψm|B̂|ψn〉∣∣ = |γ|∣∣ψm(b)ψn(b) −ψm(−b)ψn(−b)∣∣. (4) The oscillator eigenfunctions have either even or odd parity, therefore |〈ψm|B̂|ψn〉| = 0 if m and n are both even or odd, while for m even and n odd, and vice versa, we have∣∣〈ψm|B̂|ψn〉∣∣ = 2|γ|∣∣ψm(b)∣∣∣∣ψn(b)∣∣ ≤ 2|γ|C̃m−1/4n−1/4. (5) In the last line we have exploited an inequality given by Mityagin and Siegl [9] which is valid for 2(2n+ 1) ≥ b2. The prediction then is that eigenvalues can only move to a distance M, where M is a constant, uniform for all eigenvalues, and in particular, that there exists an n0 such that for n ≥ n0 all eigenvalues stay in disjoint disks with shrinking radii, with the shrink rate being bounded from above by Cn−1/2. In fact, one of the authors of Ref. [9] has pointed out [10] that for the case of two delta potentials the shrink rate can be estimated even more precisely to behave as log(n)/n3/2. Note that no such statements can be made from their theorems for the case that the nonlinearity is also included as a perturbation. 3. Eigenvalue spectra The (real and complex) eigenvalues µ of the Hamilto- nian (3) and its eigenstates ψ are obtained, in both the linear and the nonlinear case, by integrating the wave functions outward from x = 0, and varying the initial values of Re ψ(0), ψ′(0) ∈ C, and µ ∈ C to find square-integrable normalised solutions (the arbitrary global phase is exploited by choosing Im ψ(0) = 0). Figure 2 shows for the lowest 30 levels the real and imaginary parts of the eigenvalues of (3) (for g = 0) as 117 D. Haag, H. Cartarius, G. Wunner Acta Polytechnica 0 10 20 30 40 50 60 70 0 5 10 15 20 25 30 µ r γ −10 −5 0 5 10 0 5 10 15 20 25 30 µ i γ Figure 2. Real parts (top) and imaginary parts (bottom) of the eigenvalues of (3) (with g = 0) for b = 0.2, evolving from the unperturbed levels with n = 0, . . . , 29, in dependence on the strength γ of the non-Hermiticity. functions of the strength γ of the non-Hermiticity, for a position of the delta functions at b = 0.2, close to the centre of the oscillator. One recognises that successive pairs of eigenvalues coalesce at branch points, from where onwards they turn into complex conjugate pairs. The branch points are shifted to larger values of γ as one goes up in the spectrum. What is surprising is that both the real and the imaginary part of the complex eigenvalues emerging from the branch point of the eighth and ninth excited state experience huge shifts, but eventually saturate, like the other levels. For the ninth excited state, the real parts remain approximately constant beyond the branch points, and the imaginary parts quickly tend to zero. In Figure 3 we move the delta functions further away from the centre of the harmonic oscillator to the classical turning points b = 1 and √ 7 of the unperturbed levels with n = 0 and n = 3, respectively. At the turning points the unperturbed wave functions enter into the classically forbidden, i.e. exponentially decreasing region. It is therefore no surprise that for b = 1 the ground state no longer “feels” the delta functions, and no longer unites with the first excited state at a branch point. Rather its eigenvalue remains real for any strength of the non-Hermiticity. For b = √ 7 it is the states with n = 0, 1, 2, 3 which exhibit this property. We note that similar behaviour was found by Jakub- ský and Znojil [19] in their PT -symmetric square well 0 10 20 30 40 50 60 70 0 5 10 15 20 25 30 µ r γ 0 10 20 30 40 50 60 70 0 5 10 15 20 25 30 µ r γ Figure 3. Real parts of the eigenvalues of (3) (with g = 0), evolving from the unperturbed levels with n = 0, . . . , 29, with the delta functions placed into the classical turning point of the unperturbed ground state (b = 1, top) and the third excited state (b = √ 7, bottom). model. In their terminology, energy levels which coa- lesce at a branch point and turn complex are called “fragile”, while energy levels whose eigenvalues remain real for any strength of the perturbation are called “robust”. The latter also include states where the delta functions happen to be at or close to a node of the wave function, and therefore remain unaffected by the perturbation. Examples of this can be seen in Figure 3. In Figure 2 the ground state coalesces with the first excited state, while in Figure 3, at b = 1 (top panel), it has become a single real level for any value of γ, and the first excited state coalesces with the second excited one. The question arises how the transition between the different coalescence behaviour occurs. This is illustrated in Figure 4, where the real and imaginary parts of the eigenvalues emerging from the ground state and the first two excited levels are shown as functions of γ, for three positions of the delta functions around b ≈ 0.9. It is evident that at b = 0.897 the ground and first excited state still coalesce, giving rise to a pair of complex conjugate eigenvalues which “collides” with the real eigenvalue of the second excited level. The latter then turns into another pair of complex conjugate eigenvalues. The behaviour is similar for b = 0.915 but here the pair of complex eigenvalues resulting from the merger of the ground and first excited state disappears by splitting into 118 vol. 54 no. 2/2014 BEC with PT-Symmetric Double-Delta Loss and Gain: Test of Estimates 1 2 3 4 5 0 1 2 3 4 5 µ r γ b=0.897 b=0.915 b=0.925 −2 −1.5 −1 −0.5 0 0.5 1 1.5 2 0 1 2 3 4 5 µ i γ b=0.897 b=0.915 b=0.925 Figure 4. Change in the coalescence and bifurcation behaviour of the three lowest eigenstates as a function of γ at three different positions of the delta functions in the vicinity of b = 0.9. (Top: real parts, bottom: imaginary parts of the eigenvalues.) two real eigenvalues the lower of which remains real for any γ, while the higher after a small interval of γ coalesces with the real eigenvalue of the second excited state and a new pair of complex eigenvalues is born. Finally, at b = 0.925 the transition has occured, the ground state has become a single real level, and the first two excited states come together. We note again that similar behaviour was found by Krejčiřík and Siegl [20] in their studies of eigenvalues in a square well with Robin boundary conditions (cf. Figure 6 in [20]). We now proceed to results for nonvanishing non- linearity. Figure 5 shows the spectrum for a value of g = 5 with the delta functions placed at b = 1. The overall behaviour is similar (many states coalesce at branch points), but there are significant differences. Like in our previous studies of PT -symmetric Hamil- tonians with nonlinearity [2–6, 25], pairs of complex conjugate eigenvalues (imaginary parts not shown) appear before the branch points are reached. Again there exist “robust” levels whose eigenvalues remain real for any value of γ. 4. Eigenvalue shifts For vanishing nonlinearity, the upper part of Figure 6 shows, for γ = 1 and different positions of the delta functions, the shifts of the eigenvalues as a function of the quantum number n in comparison with the predicted shrink rates of n−1/2 of [9] and the improved 0 10 20 30 40 50 60 0 1 2 3 4 5 6 7 8 9 µ r γ Figure 5. Real parts of the eigenvalues of (3) with fi- nite nonlinearity g = 5, evolving from the unperturbed levels with n = 0, . . . , 29, for b = 1 in dependence on the strength of the non-Hermiticity. Blue lines denote purely real eigenvalues, whereas real parts of complex eigenvalues (single lines) are shown in red. 10−4 10−3 10−2 10−1 100 1 10 100 ∆ µ n ∝ n−0.5 ∝ log(n)/n3/2 b = 1 b = 2 b = 3 b = 5 0.1 1 1 10 100 ∆ µ n ∝ n−0.5 ∝ log(n)/n3/2 b = 1 b = 2 b = 3 b = 5 Figure 6. Shift of the eigenvalues for γ = 1 and different positions of the delta functions in compar- ison with the mathematically estimated n−1/2 and log(n)/n3/2 dependence. Top: vanishing nonlinearity, g = 0, bottom: g = 2. estimate of log(n)/n3/2 [10]. It can be recognized that the n−1/2 dependence is indeed an upper bound on the shrink rate, and, moreover, that the improved estimate is in excellent agreement with the numerical data, irrespective of the position of the deltas! The picture changes when the nonlinearity is switched on. The bottom part of Figure 6 shows eigenvalue shifts for g = 2, and different positions of the imaginary delta potentials. From the comparison of the scales of the vertical axes in Figure 6 it can 119 D. Haag, H. Cartarius, G. Wunner Acta Polytechnica 10−2 100 102 1 10 ∆ µ n γ = 5 γ = 10 γ = 15 γ = 20 Figure 7. Shift of the real and complex eigenvalues for growing values of γ for b = 0.2, and g = 0. be seen that with nonlinearity the shifts even for the highest states are still bigger than 0.1, while they have already dropped well below 10−2 in the linear case. Furthermore, the shrink rate is found to be slower (approximately proportional to n−0.37) than predicted by both rigorous mathematical estimates for the linear case. We therefore find significant differences in the shrink rates with and without nonlinearity. For small γ, all eigenvalues are still real. The ques- tion therefore suggests itself how the situation changes when eigenvalues are involved that are shifted into the complex plane, which happens for increasing γ. Figure 7 shows for the eigenvalue spectra with b = 0.2, vanishing nonlinearity and different values of γ the distances |∆µ| of the real or complex eigenvalues from their original values (cf. also Figure 2). The outlier at n = 8 and 9 is caused by the giant change of the real and imaginary parts of the complex eigenvalues beyond the branch point of the corresponding states observed already in Figure 2. A similar outlier occurs around n = 75, and a monotonous decrease of the eigenvalue shifts in our calculations only sets in for n ≈ 80. Finally, in Figure 8 we investigate the effect of the nonlinearity on the eigenvalue shifts for γ = 1.5 and b = 0.5 and b = 1. We find that the eigenvalue shifts oscillate around straight lines with slopes -0.37, ir- respective of the strength of the nonlinearity. The amplitude of the oscillations, however, and their num- bers depend on the position of the delta functions. Again the slope of the estimate proportional to n−1/2, also shown in the Figure, is steeper than the actual slopes found in the numerical results. 5. Conclusions We have carried out a numerical analysis of a PT - symmetric double delta perturbation of the harmonic oscillator. We have also considered the case were in addition a Gross-Pitaevskii nonlinearity proportional to |ψ|2 is present. With the latter, the system can be considered as a model of a Bose-Einstein condensate in a double well with loss and gain of atoms. 10−1 100 1 10 ∆ µ n g = 1 g = 2 g = 3 g = 4 10−1 100 1 10 ∆ µ n ∝ n−0.37 ∝ n−0.5 Figure 8. Shift of the eigenvalues N = 2n + 1 for four different values of the nonlinearity at γ = 1.5 for b = 0.5 (top) and b = 1 (bottom). We have checked that the Mityagin-Siegl criterion for the perturbed eigenfunctions to form a Riesz basis is fulfilled, and compared rigorous mathematical es- timates for the shrink rate of the eigenvalue shifts in dependence on the harmonic oscillator quantum num- ber n for various strengths of the non-Hermiticity and the nonlinearity. We have verified that in the linear case the mathematical prediction for the shrink rates proportional to 1/n1/2 is a valid estimate, and that the improved estimate proportional to log(n)/n3/2 is in excellent agreement with the behaviour of the shrink rates found in the numerical results. By contrast, with nonlinearity we find slopes of approximately −0.37, less steep than both mathematical estimates. Evi- dently, to derive estimates also for the nonlinear case remains a mathematical challenge. A peculiarity that is found is the occurrence of outliers in the eigenvalue spectra which appear beyond branch points with unusually large real and imaginary parts of their complex conjugate eigenvalues. They are the reason why for growing strength of the non- Hermiticity the asymptotic shrink rate behaviour is attained only for high values of n. The nature of these outliers and their mathematical importance should certainly be clarified in future studies. Acknowledgements We are grateful to Petr Siegl for helpful explanations and discussions. We are also grateful to Boris Mityagin for very helpful remarks, and for communicating his improved estimate for the two-delta potentials. Furthermore we 120 vol. 54 no. 2/2014 BEC with PT-Symmetric Double-Delta Loss and Gain: Test of Estimates thank two anonymous referees for valuable comments. In particular, one referee points out that a deeper analytic insight in the bottom of the spectra could probably be obtained using the strategies described in Ref. [22] for the self-adjoint case, if adapted to the PT -symmetric perturbation. This is certainly also a useful suggestion for future work. References [1] S. Klaiman, et al. Visualization of branch points in PT -symmetric waveguides. Phys Rev Lett 101:080402, 2008. doi: 10.1103/PhysRevLett.101.080402 [2] 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 [3] 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 [4] 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. [5] 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 [6] 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 [7] 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 [8] J. Adduci, et al. Eigensystem of an L2-perturbed harmonic oscillator is an unconditional basis. Cent Eur J Math 10:569, 2012. doi: 10.2478/s11533-011-0139-3 [9] B. Mityagin, et al. Root system of singular perturbations of harmonic oscillator type operators, 2013. ArXiv:1307.6245. [10] B. Mityagin, private communication, January 2014. [11] H. F. Jones. The energy spectrum of complex periodic potentials of the Kronig-Penney type. Phys Lett A 262:242, 1999. doi: 10.1016/S0375-9601(99)00672-6 [12] Z. Ahmed. Energy band structure due to a complex, periodic, PT -invariant potential. Phys Lett A 286:231, 2001. [13] E. Demiralp. Bound states of n-dimensional harmonic oscillator decorated with dirac delta functions. J Phys A 38:4783, 2005. doi: 10.1088/0305-4470/38/22/003 [14] E. Demiralp. Properties of a pseudo-Hermitian Hamiltonian for harmonic oscillator decorated with dirac delta interactions. Czech J Phys 55:1081, 2005 doi: 10.1007/s10582-005-0110-2 [15] A. Mostafazadeh. Delta-function potential with a complex coupling. J Phys A 39:13495, 2006. doi: 10.1088/0305-4470/39/43/008 [16] 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/12503 [17] 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 [18] H. Uncu, et al. Bose-Einstein condensate in a harmonic trap with an eccentric dimple potential. Las Phys 18:331, 2008. [19] 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 [20] D. Krejčiřík, et al. PT -symmetric models in curved manifolds. J Phys A 43:485204, 2010. doi: 10.1088/1751-8113/43/48/485204 [21] P. Siegl. PT -symmetric square well-perturbations and the existence of metric operator. Int J Theor Phys textbf50:991, 2011. doi: 10.1007/s10773-010-0593-x [22] S. Fassari, et al. On the spectrum of the Schrödinger equation of the one-dimensional harmonic oscillator perturbed by two identical attractive point interactions. Rep Math Phys 69:353, 2012. [23] E. P. Gross. Structure of a quantized vortex in Boson systems. Nuovo Cimento 20:454, 1961. [24] L. P. Pitaevskii. Vortex lines in an imperfect Bose gas. Sov Phys JETP 13:451, 1961. [25] 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 121 http://dx.doi.org/10.1103/PhysRevLett.101.080402 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.1002/prop.201200080 http://dx.doi.org/10.1088/1751-8113/46/37/375301 http://dx.doi.org/10.1103/physreva.87.051601 http://dx.doi.org/10.2478/s11533-011-0139-3 http://dx.doi.org/10.1016/S0375-9601(99)00672-6 http://dx.doi.org/10.1088/0305-4470/38/22/003 http://dx.doi.org/10.1007/s10582-005-0110-2 http://dx.doi.org/10.1088/0305-4470/39/43/008 http://dx.doi.org/10.1088/1751-8113/42/12/12503 http://dx.doi.org/10.1088/1751-8113/43/14/145301 http://dx.doi.org/10.1007/s10582-005-0115-x http://dx.doi.org/10.1088/1751-8113/43/48/485204 http://dx.doi.org/10.1007/s10773-010-0593-x http://dx.doi.org/10.1088/1751-8113/46/27/275307 Acta Polytechnica 54(2):116–121, 2014 1 Introduction 2 Bose-Einstein condensate in a PT-symmetric harmonic trap 3 Eigenvalue spectra 4 Eigenvalue shifts 5 Conclusions Acknowledgements References