AP07_2-3.vp 1 Introduction In physics there are usually one or more free parameters which specify some properties of the system. Aside from trivial parametric dependences that can be determined directly, e.g., by a change of units, variation in the values of some parameters can qualitatively affect the system’s behaviour. In quantum physics, where almost all physical properties are related to spectral properties of a chosen set of operators (in particular the Hamiltonian), qualitative changes of spectra due to varying parameters of the Hamiltonian are of promi- nent importance. ��-symmetry, which was introduced as a concept gene- ralising the usually demanded self-adjointness of the Hamiltonian, does not per se provide the reality of its spec- trum and consequently of the energies. Even if the spectrum is real it need not to be so for all values of the considered parameters. Therefore determining the regions in the para- metric space where the spectrum’s reality holds is a key step in finding the limits of the physical significance of any ��-sym- metric theory(1). Most attention has been concentrated on discrete spectra, for the following reasons. First, it is usually more convenient (at least for mathematically less careful and less rigorous phys- icists) to deal with isolated eigenvalues than with the continu- ous part of the spectrum, where all manipulations become in a way more mathematically intricate. Second, there are many systems whose spectrum is fully discrete and others where the essential spectrum is insensitive to parametric change. The points of eigenvalues’ complexification are thus the investi- gated limits of physical relevance. After examining many of the “classical” examples of ��-symmetric systems an appar- ently regular pattern was observed – a square-root singularity structure and a Jordan-block degeneracy. Though this is beyond the scope of this article, it can be useful for the reader to consider one more application of studying complexification and its properties. Complex (albeit not necessarily ��-symmetric) Hamiltonians do also arise with extending parameters of the Hermitian systems into the complex domain. The structure of the spectrum and espe- cially the distribution of the exceptional points (which are, in this case up to slight generalisation, the points of the spec- trum’s complexification) then influences the phase transitions and the chaotic behaviour of many systems. See e.g. [2] or [3] (and refs. therein). Outside of quantum mechanics, ��-symmetric Hamil- tonians are of use in the framework of magnetohydrody- namics. Here again the exceptional points and level crossings (we will discuss their close relatedness later) have a direct physical significance. For the connection between magneto- hydrodynamics and ��-symmetry, see [4, 5]. The concrete na- ture of the complexification is important, e.g., if one has to construct a perturbation theory around these singularities [6]. The purpose of this paper is to summarise some well- -known facts about complexification of energies and also to discuss some of the other (than usual) possible scenarios of complexification, especially for Dirac operators. In the first section we will make clear the terminology and compare par- ticular properties of self-adjoint and ��-symmetric parame- trically dependent Hamiltonians. In the second section, we will present the situation in relativistic models. 2 The basics of complexification In the following, we will investigate the spectra of opera- tors that depend on a single parameter. These operators will be called Hamiltonians, as is usual in the ��-symmetric con- text, even if the matrices used for illustration in this section are rather toy models than realistic generators of time evolu- tion. However, the physical background of what is discussed in the second section will be rather obvious. The investigated parameter will be denoted by c, the Hamiltonian H(c) and its eigenvalues E, unless stated other- wise. The point in the parametric space where complexi- fication occurs will be called the exceptional point (2) (EP) and denoted by cep; complexification means that some chosen subset (usually consisting of two eigenvalues) of the spectrum � �� H( )c is real for c c� ep and imaginary for c c� ep or vice versa. To avoid unnecessary complications we will examine the most common situation when the Hamiltonian is of the form H H V( )c c� �0 (1) Note that in spite of being quite restrictive within all imaginable parametric dependences, form (1) is useful in a broad class of physical applications. Obviously, one can get other than linear parametric dependences by means of reparametrisation. 50 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 47 No. 2–3/2007 How Do Energies Complexify? H. Bíla Some particular properties of the parametric dependence of eigenvalues with emphasis on their complexification are discussed. The non-diagonalisability of ��-symmetric matrix Hamiltonians in exceptional points is compared with level-crossing prohibition of Hermitian systems. For non-matrix Hamiltonians, the different way of complexification between Klein-Gordon and Dirac Hamiltonians is demonstrated. Keywords: ��-symmetry, exceptional points, Dirac equation All discussed Hamiltonians will be ��-symmetric. This means that there exists an anti-linear operator, conventionally written as a product of time reversal � and parity reflection �, (3) which commutes with the Hamiltonian: �� ��H H( ) ( )c c� (2) It is well known that condition (2) itself provides that the discrete spectrum consists of real eigenvalues and com- plex-conjugated pairs. Thanks to the simple parametric de- pendence of H it is reasonable to estimate that the eigenvalues are continuous functions of c.(4) From these two facts one can infer that at the exceptional point at least two eigenvalues must merge. Hence it could be useful to generalise the notion of EPs to those c where the spectrum is degenerated even if it is real (and non-degenerated) in the neighbourhood of the point. 2.1 Matrices and avoided crossings The simplest case of a ��-symmetric system is a paramet- rically dependent (finite-dimensional) matrix. Examples of these have been used many times to demonstrate various properties of ��-symmetric systems. For some infinite-di- mensional systems, non-self-adjointness is essential only in the subspaces spanned on the lowest-energy eigenvectors [7, 8] and matrix representations of the Hamiltonian on this subspace can be considered as a good approximation to the complete problem. Because for matrices of form (1) the eigenvalues are al- ways continuous functions of c and their total number is fixed, complexification can clearly happen only in two ways. In the first case, at the exceptional point there is an n-dimensional degenerate subspace with n independent eigenvectors. In the second case, there are only n k� independent eigenvectors with k � 0 and there is thus a Jordan-block structure. The important thing now is that the former case, let us call it diagonalisable EP, is much rarer than the latter (non-diago- nalisable). To see this, let us first examine the simplest 2×2 matrix case. A two-dimensional diagonalisable matrix with a degenerated spectrum is a multiple of the unit matrix and so if one assumes diagonalisable EP at cep, then H I( )c kep � and thus H I V0 � � �k c c( )ep . (3) The consequence is that H0 and V commute. This is something that is not satisfied in almost all physical situa- tions. However it must be emphasised that although the commutativity of H0 and V is a sufficient condition for the existence of diagonalisable EPs in any dimension of Hilbert space, it is a necessary condition for that only in the case of two dimensional systems. The stated argument cannot be easily generalised to more dimensions, and it needs a differ- ent approach to see the reasons why diagonalisable EPs are “prohibited”. Hermitian matrices are always diagonalisable and thus the only EP that can occur in a Hermitian parametrically depend- ent system is a diagonalisable one, which is usually referred to as a level crossing. That level crossings are in some way avoided for Hermitian systems is well known, but because the mecha- nism of this is exactly the same as the mechanism protecting against diagonalisable complexification in the ��-symmetric case, let us formulate the statement more precisely. For the purposes of this article, the following formulation of the “no-crossing” theorem will be sufficient: Let � be a set of all N-dimensional Hermitian matrices and S N: � � 2 be a mapping that maps any V � with entries vij to the vector S v v v v v v vN NN( ) ( , , , , , , , , )V � � � � �11 12 12 1 22 23� � . (4) Let H0 then be a fixed Hermitian operator and �C be a subset of � such that for all V �C the operator H V0 � c has a level crossing for at least one c � . Then S( )�C has zero mea- sure in �N 2 .(5) The message of the theorem is that parametrically de- pendent matrices of the form (1) which have a level crossing are extremely rare in the sense that they constitute a zero- -measure subset of all possible matrices of the considered type. A sketch of the proof can be outlined as follows: The characteristic polynomial �( ) det(E c E� � �H V I)0 (5) of a Hermitian matrix is real for all real c and E. It must have a multiple root in an exceptional point, which is otherwise stated as �( )cep � 0 , (6) �� �E c ep � 0 . (7) The dependence of � on c is polynomial and thus smooth and therefore also �� �c c ep � 0 ; (8) otherwise � would have complex roots at cep � � for some non-zero small �. Now, having a given H0, let us choose any cep and Eep and find any matrix � such that H V0 � c has an EP at c0 with the energy Eep. (6) One may ask what happens with the eigen- values if the matrix elements of V (or better the components xi of the vector x S( )V ) and the parameter c are slightly var- ied around the selected point. In the ( )N2 2� -dimensional space of variables xi, c and E the eigenvalues are located on the subset where �( , , )x c E , therefore one needs to solve the differential equation 0 � � � � �d d d d� � �� � �� � ( )x E E c cx . (9) From (6) and (7) it follows that if one starts at the EP the last two terms in (9) vanish and one gets 0 � �( )x � dx (10) This is an equation which defines an ( )N2 2� dimensional manifold in �N 2 . It is obvious now that the N2-dimensional vectors x for which S � 1( )x �C are restricted to lie on such manifolds clearly of zero measure since they have one dimen- sion less than the space of all Hermitian N×N matrices.(7) © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 51 Acta Polytechnica Vol. 47 No. 2–3/2007 What are the differences if one takes the ��-symmetric case instead of Hermitian one? The number of free real pa- rameters determining a ��-symmetric matrix is the same as for a Hermitian matrix and its characteristic polynomial is also real. The crucial difference is that since there is no guar- anteed reality of the spectrum, condition (8) does not hold. On the other hand since we know [9] that a diagonalisable ��-symmetric matrix with real eigenvalues can be trans- formed to a Hermitian one by a similarity transformation, and because similar matrices have the same characteristic polynomials, condition (8) is still valid for diagonalisable EPs of ��-symmetric matrices. However the argument does not apply for non-diagonalisable EPs, where an analogous simi- larity transformation does not exist. In this case, (9) can be solved by a function c(x) that determines the position of the selected EP under an arbitrary change of parameters xi . 2.2 Square-root dependence To illustrate the behaviour of eigenvalues on a concrete example, let us consider the simple two-dimensional matrices H � � � � � � � � � � �� c de de c iq iq 1 2 (11) with c d qi, , �. (8) The eigenvalues are given by 2 41 2 1 2 2 2E c c c c d� � � � �( ) . (12) The square root in (12) is a typical example of parametric dependence near the non-diagonalisable EP. Since (11) is a very simple (i.e. two-level) system the dependence is exactly the square root; for larger matrices of type (1) the square-root function is modified by some non-singular c-dependent fac- tor, but still lim c c E c � �� ep d d . (13) The essential fact is that in complex neighbourhood of the exceptional point the eigenvalues form a structure with two Riemann sheets typical for the square root – this is univer- sally true for matrices in all EPs where two levels cross. The scheme also holds for a large class of Schrödinger operators with ��-symmetric potentials. A square-root singularity was attested for example in the founding problem of ��-symme- try – the Hamiltonians � � � d d 2 2 2 x ix c( ) (14) of Bender, Boettcher and Meisinger [10]. Another example is the one-dimensional Laplacian at the interval ( , )0 d discussed in [11] with ��-symmetric boundary conditions(9) � � � �� �( ) ( ) ( )0 0b ic , � � �� �( ) ( ) ( )d b ic d (15) The system defined in such a way has several interesting properties: The total number of complex eigenvalues does not exceed one pair. The exceptional point usually has a Jor- dan-block degeneracy and the square-root-like behaviour of the energy, but the situation can be made slightly more com- plicated for example by keeping b c n2 2 1 2� � � fixed and varying c, which leads to unusual triple (still non-diagonalis- able) degeneracy in the EP. An even more extraordinary possibility is to put b � 0 and observe the dependence on c. In this setting the spectrum is real for all c and only one energy level is not constant; the exceptional points occur at k n� �, where n � and k E� . 52 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 47 No. 2–3/2007 Fig. 1: The eigenvalues of matrix (11) with c c c1 2� � and d e iq� � 1 (solid lines). The dashed lines represent an analogous Hermitian matrix where both off-diagonal ele- ments are �1. The graph illustrates the typical behaviour of energy levels of both Hermitian systems (avoided cross- ing at c � 0) and ��-symmetric systems (complexification at exceptional points c � �1 ). Fig. 2: The energy dependence of the Laplacian with boundary conditions (15). The left graph shows the level crossings for b � 0. The sit- uation on the right represents varying b for fixed value b c2 2 2 4� � . . Solid lines represent real eigenvalues while dashed and dotted lines are the real and imaginary parts of one of the complex eigenvalues respectively (the second one is conjugated). The values of the wave vector k E� are shown. Because actually no complexification occurs also (13) does not hold, but a Jordan-block structure is still maintained. As a summary to the first section it is worth noting that the merging of two (or ocassionally more) eigenstates and eigenvalues in the point of complexification, while being obligatory for matrices, is prevalent but not universal for Schrödinger operators. Complications come into play when the potential involves a singularity, typically a divergent term proportional to 1 x. In such situations the eigenfunction cor- responding to one of the real eigenvalues expected to merge at the EP ceases to be square integrable and thus in this point one sees a splitting of one real energy into two complex ener- gies (see for example [12]). We will also discuss this later with the radial Dirac equation. 3 Relativistic systems The equations of relativistic quantum mechanics provide a natural resource of systems with complexification phenom- ena. In contrast to non-relativistic ��-symmetry, here one does not need to add an “artificial” complex term and com- plexification is encountered within undoubtedly physically relevant problems. The point of complexification is usually regarded as the furthermost boundary of quantum-mechani- cally describable configuration. In the following subsections we will discuss in more detail three different relativistic models. 3.1 Dirac square well One of the simplest interactions that can be imagined is represented by a finite square well potential. The first dis- cussed relativistic system will be the spin 1 2 particle in such potential. The one-dimensional Dirac Hamiltonian in one suitable representation reads H � � � � � � � � �� � � �� m c x i i m c x d x x d � � � � ( , ) ( , ) ( ) ( ) 0 0 , (16) where c and d are the well’s depth and width, respectively, and is the mass of the particle, throughout the following deliber- ately put to be equal to one. The domain of the Hamiltonian (16) is � �DomH � � � �� �L L L L2 2 2 2( ) ( ) ( ) ( )� � � � . (17) The Hamiltonian has as usual the Dirac Hamiltonian con- tinuous spectrum � c � �� � � �( , ) ( , )1 1 which does not de- pend on c, d. Our interest is concentrated on the bound states. The search for stationary states consists of solving the prob- lem on the intervals II � ��( , )0 , I dII � ( , )0 and I dIII � �( , ). The normalisable solutions for energy E are � � � � � I I II II II III III � � � � � � � � K K K K x x x x e e e e , , . (18) where K refers to still unspecified complex constants and � � � � � � � � � � � � � � 1 1 ( ), , , , c E P E P P (19) the square roots are taken positive or positive imaginary (that is purely imaginary with positive imaginary parts). The matching conditions in 0 and d involve then only the wave-function, not the derivative, as H contains only the first derivative term. After eliminating all K s we get the equation ( )1 1� � �� ��e d (20) with � � � � � � P P . From the definitions above it can be seen that � is real positive or imaginary positive, but no real positive � can fulfill equation (20). It follows that E lies in the interval (�1, 1) and thus in the gap of the continuous spectrum, as might be expected. Now it would be useful to take a logarithm of both sides of (20) to get a real equation � � 2 1 1 1 1 2 12� � � � � � � � � �arctan ( )( ) ( )( ) ( ) mod E c E E c E d c E . (21) By investigating (21) we arrive at the following observa- tions about the dependence of the energies on parameters c and d: 1. Due to the symmetry of the system the change c c � yields E E � and thus one can restrict one’s attention to c � 0 without loss of generality. 2. Each energy decreases as c or d increases. 3. For the i-th bound state and fixed d there are critical val- ues c di min ( ) and c di max( ) such that the bound state exists only if c c ci i ( , ) min max . The energy of the bound state “sinks” into the continuous spectrum at these points: lim , lim . min max c c i c c i i i E E � � � 1 1 (22) 4. By analogy for fixed c there are d ci min ( ) and d ci max( ) and the bound state exists only if d d di i ( , ) min max . 5. d c i c ci min ( ) ( ) � � 2 2 � and d c i c ci max( ) ( ) ( ) � � � 2 1 2 � . When c � 2 the value d ci max( ) � � . The energies are numbered from 0 to �. We can formally consider the points with parameters c d ci, ( ) min or c d ci, ( ) max as exceptional points, but in these points no complexification takes place. The spin-1/2 particle in a square well thus behaves “correctly” and exhibits no complexification for strong potential; this can be placed in contrast to a scalar particle in an analogous potential de- scribed in the next subsection. 3.2 Klein-Gordon square well It is well known that a scalar particle is described by the Klein-Gordon equation. As for the involved potential we will consider the minimal coupling, and thus the equation for sta- tionary states has the form © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 53 Acta Polytechnica Vol. 47 No. 2–3/2007 (( ( )) ) ( )( , )E c x M xd x� � � �� � �0 2 2 2 0 , (23) where c is the well’s depth and d its width. It is worth mentioning that it is possible to describe system (23) by a Hamiltonian of form (1) (see [13]). The price paid is that one must introduce a two-component formalism, which leads to some ambiguities as long as the step from one to two components is not unique, and the resulting Hamiltonian is non-Hermitian(10). Since we are mainly interested in the spec- trum we will not follow this way. After some manipulations with matching conditions at 0 and d and logarithming as in the Dirac case, it is easy to obtain the secular equation for the eigenvalues, which reads arctan ( )( ) ( )( ) ( ) mod 1 1 1 1 2 1 2 2� � � � � � � � � E E c E c E d c E � . (24) Although (24) resembles (21), the solutions to it behave differently at larger values of c and d. In addition to levels analogous to those of the Dirac Hamiltonian, there is a sec- ond set of energies emerging from the lower continuum with increasing c or d that merge with the former in a traditionally shaped exceptional point (Fig. 3). This effect is enabled by the non-self-adjointness of the Hamiltonian, in contrast with the Dirac problem. 3.3 Relativistic coulomb Hamiltonian The last example is intended to show that also the Dirac equation can lead to complex energies. It is the well-known radial Coulomb Hamiltonian H � � � � � � � � � � � � � � � � � � c r m r r c r m r r � � � � , (25) where � �� 1 2 3, , ,� characterises the angular momentum and m is the mass, once more without loss of generality set to one in what follows. We will not discuss the algorithm of the solution since the reader is assumed to know it, and instead we write down the resulting formula E c nn � � � � � � �� � � �� � 1 1 2 � , (26) �� �2 2c . (27) Obviously c > � leads to complex values of En and all the levels with the same � complexify simultaneously. The energy dependence is clearly of square-root type, however only one real eigenvalue exists at c c� ep for each complex pair at c c� ep. To see how this is possible for the symmetric Hamil- tonian (25) one needs to examine it more closely. We can restrict ourselves to s-states (� � 1) because these complexify first and there is no qualitative difference for higher �. The Hamiltonian’s domain of definition trivially must be a subset of � �D L L L L� � �� � � �� �2 2 2 2( ) ( ) ( ) ( )� � � �H , (28) (the existence of H� for all � D is understood). If one has to establish the self-adjointness of H one has first to apply per partes integration on the integrals of wave function products and demand vanishing of boundary terms. The boundary at r � � causes no complications because �( )r 0 for r � trivially due to the square integrability. One must be more careful at r 0. The square integrability of � itself clearly does not forbid non-zero values at origin, however square integrability of H� can do better. To see whether there are functions from D with non-zero limit at r � 0 one may check the asymptotics of H�; around zero one can disregard the mass-proportional terms, and one has H � � � � � � � � 1 2 1 1 1 2 2 1 1 1 1 2 � � � � � � � � � � � � � � � � � � �� � � � cr r r cr � � � . (29) 54 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 47 No. 2–3/2007 Fig. 3: The lowest energy level plotted against depth c with fixed width d � 1 5. for relativistic square wells. The spin-1/2 case on the left hand side does not complexify while the spinless case on the right hand side does; in the latter graph the two levels that merge in the EP are shown. Fig. 4: Six lowest real energies of the Dirac Hamiltonian with Coulomb potential. The line styles distinguish between different values of �. This two-component function can be square-integrable 1. either if � �� r with � � 1 2 (and thus �( )0 0� ), because then both �� and � r are proportional to r � �1 and there- fore square integrable 2. or if the whole expression (29) vanishes. The second case leads to the differential equation for asymptotical behaviour of � near the origin, and the solution to it is � � �r (30) with the same � 0 as in (27) – � still being equal to 1. The so- lution proportional to r � is not interesting for being zero at r � 0 but the solution r � diverges. For � 1 2 it is not square integrable and can thus be ruled out, but if � 1 2 (it is c � 3 2 ) this solution lies in D. In other words, in this case there are functions from D which do not vanish at r � 0 and if we make D the definition domain of H the operator would not be symmetric because of the non-zero boundary term in per partes integration. After restricting the domain of H to functions which van- ish in the origin the operator is symmetric with point spec- trum (26), but still not self-adjoint, since the domain of the adjoint operator would be the entire set D, which is now for c � 3 2 different from Dom H. From the previous discussion it follows that there are two main differences between square-well and Coulomb interac- tion within the framework of the Dirac equation. The first, connected with the long range of Coulomb potential, is that the eigenvalues do not emerge from the upper continuum as the strength of the interaction increases (as they do in the case of a square well) but they all exist for any non-zero c. A similar difference is also present in a non-relativistic treatment of the problems. The second difference, more interesting for the purposes of this paper, the existence of complexification at a certain value of c, is a consequence of the 1 r singularity at the origin. This complexification is of an extraordinary type within the family of ��-symmetric models with only one real eigenvalue forming a complex pair(11). 4 Summary The purpose of this paper was two-fold: First to discuss the simple principles behind the complexification of eigenvalues and to show the connection between avoided level crossings in the Hermitian framework and avoided diagonalisable EPs in ��-symmetry. And second, to show that in less standard sit- uations (where standard means a matrix or a Schrödinger Hamiltonian), there are more ways in which the eigenvalues can cross the boundary between the physical world and the complex realm. It is by no means intended to state that the few examples presented here constitute in any sense a com- plete list of what can happen. Rather the aim of the article was to illustrate that some quite ordinary systems can behave con- trary to the usual “��-symmetric intuition”. Remarks (1) It may be useful to note that a systematic approach to ��-symmetry cannot be probably achieved without use of the Krein space concept [1]. When the spectrum of a Krein space operator is real, a transformation to a Hilbert space operator exists and the ��-symmetric system can be treated as in standard quantum mechanics. (2) Some authors distinguish between ”exceptional“ and ”di- abolic“ points regarding the type of complexification. Here we use the first term for all types. (3) Usually but not necessarily the standard realisations of � and � are used. A particularly useful non-standard choice is � = I, which identifies symmetric and ��-symmetric operators. (4) This need not to be valid for general H0 and V, but most physical systems do not spoil the assumption. In many situ- ations it can be proved through perturbation theory. (5) The mapping S is a bijection which creates N2-dimensional vector from N2 independent real parameters determining the Hermitian N×N matrix, i.e. the real diagonal matrix elements and real and imaginary parts of the entries above the diagonal. (6) Such matrices obviously exist. A particular example can be easily constructed in the diagonal basis of H0: Any matrix which is diagonal in this basis and has values v E cii ii� � �( ( ) )ep epH0 1 at least for two different � �i N 1� is good enough. (7) We suppose the measure in the space of matrices is gener- ated by an Euclidean norm A � � Aij2 or equivalent. On contrary there can be reasons for choosing measures con- centrated on less-dimensional subspaces – usually because of some symmetry prescription for the „randomly choosen“ Hamiltonian. Under such conditions the level crossings are not forbidden. (8) It is the most general ��-symmetric two-dimensional ma- trix with respect to � defined by Pauli matrix �3 and � be- ing simply complex conjugation. (9) Strictly speaking the system defined by (15) is not of type (1). But the linear parametrisation is still valid in a general- ised sense, for the Hamiltonian’s associated sequilinear form. (10)And also not manifestly Lorentz covariant. (11)To be more specific, the statement is true if the bound- ary condition �( )0 0� is given. Without it, demanding only square integrability of � and H�, one gets a clearly non-physical situation: if c � 3 2, then all real numbers become eigenvalues with eigenfunctions proportional to confluent hypergeometric function U multiplied by decay- ing exponential. References [1] Mostafazadeh, A.: Krein-Space Formulation of ��-Sym- metry, ���-Inner Products, and Pseudo-Hermiticity, Czech. J. Phys. Vol. 56 (2006), p. 919. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 55 Acta Polytechnica Vol. 47 No. 2–3/2007 [2] Heiss, W. D., Müller, M., Rotter, I.: Collectivity, Phase Transitions and Exceptional Points in Open Quantum Systems, Phys. Rev. Vol. E58 (1998), p. 2894. [3] Cejnar, P., Heinze, S., Dobeš, J.: Thermodynamic Anal- ogy for Quantum Phase Transitions at Zero Tempera- ture, Phys.Rev. Vol. C71 (2005), p. 011304. [4] Günther, U., Stefani, F., Gerbeth, G.: The MHD �2–Dy- namo, �2–Graded Pseudo-Hermiticity, Level Crossings and Exceptional Points of Branching Type. Czech. J. Phys. Vol. 54 (2004), p. 1075. [5] Günther, U., Stefani, F., Znojil M.: MHD �2–Dynamo, Squire Equation and ��-Symmetric Interpolation be- tween Square Well and Harmonic Oscillator, J. Math. Phys. Vol. 46 (2005), p. 063504. [6] Günther, U., Rotter, I., Samsonov, B. F.: Projective Hilbert Space Structures at Exceptional Points, arXiv:0704.1291. [7] Mostafazadeh, A., Batal, A.: Physical Aspects of Pseudo- -Hermitian and ��-Symmetric Quantum Mechanics, J. Phys. Vol. A37 (2004), p. 11645. [8] Quesne, C., Bagchi, B., Mallik, S., Bíla, H., Jakubský, V., Znojil, M.: ��-Supersymmetric Partner of a Short- -Range Square Well, Czech. J. Phys. Vol. 55 (2005), p. 1161. [9] Mostafazadeh, A.: Exact ��-Symmetry is Equivalent to Hermiticity, J.Phys. Vol. A36 (2003), p. 7081. [10] Bender, C. M., Boettcher, S., Meisinger, P.: ��-symmet- ric Quantum Mechanics, J. Math. Phys. Vol. 40 (1999), p. 2201. [11] Krejčiřík, D., Bíla, H., Znojil, M.: Closed Formula for the Metric in the Hilbert Space of a ��-Symmetric Model, J. Phys. Vol. A39 (2009). p. 10143. [12] Lévai, G.: Comparative Analysis of Real and ��-Sym- metric Scarf Potentials, Czech. J. Phys. Vol. 56 (2006) p. 953. [13] Feshbach, H., Villars, F.: Elementary Relativistic Wave Mechanics of Spin 0 and Spin 1/2 Particles, Rev. Mod. Phys. Vol. 30 (1958). Mgr. Hynek Bíla phone: +420 266 173 280 email: hynek.bila@ujf.cas.cz Nuclear Physics Institute Academy of Science of the Czech Republic 250 68 Řež, Czech Republic 56 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 47 No. 2–3/2007