Supersymmetry and exceptional points

A conceptual bridge is provided between SUSY and the three-Hilbert-space upgrade of quantum theory a.k.a. ${\cal PT}-$symmetric or quasi-Hermitian. In particular, a natural theoretical link is found between SUSY and the presence of Kato's exceptional points (EPs), both being related to the phenomenon of degeneracy of energy levels. Regularized spiked harmonic oscillator is recalled for illustration.


Introduction
Ten years ago people perceived the concepts of supersymmetry and of exceptional points as the two remote if not fully separate subjects of research. For example, during the international workshop "Supersymmetric Quantum Mechanics and Spectral Design" (July 18 -30, 2010, Benasque, Spain, organized by A. Andrianov et al [1], with proceedings [2]), just one talk was devoted to the role of exceptional points in supersymmetric quantum systems. In parallel, no speaker mentioned supersymmetry during the workshop "Exceptional points in physics" (November 2 -5, 2010, NITheP, Stellenbosch, South Africa, organized by W. D. Heiss [3]). Also the same year, the participants of the 9th continuation of the series "Pseudo-Hermitian Hamiltonians in Quantum Physics" (June 21 -24, 2010, Zhejiang University, Hangzhou, China, organized by J.-D. Wu [4], with proceedings [5]) did not yet identify their field of interest as an emergent bridge between the two concepts, from either the mathematical and physical points of view.
At present, the situation is different. The use of the mathematical notion of exceptional points (EPs, [6]) makes an impression of being, in multiple branches of physics, ubiquitous [7,8]. The same growth of popularity became also characteristic for the concept of hidden Hermiticity alias pseudo-Hermiticity [9] or, better, quasi-Hermiticity [10,11], most widely known under the physicsemphasizing nicknames of "PT −symmetry" [12] or "gain-loss balance" [13]. In comparison, the concept of supersymmetry (SUSY, [14]) seems, undeservedly, overshadowed by these developments.
In our present paper we intend to demonstrate that the concept of SUSY witnesses a revival, in particular, due to its remarkable overlap with both of the EP-and PT −related areas of research. The new conceptual bridges are emerging between these notions in a way which we intend to describe in what follows. On the SUSY side, in essence, the core of the mutual influence lies in the emergence of the new, unconventional realizations of the SUSY algebra in quantum mechanics. One witnesses there a shift of attention from Hermitian to non-Hermitian Hamiltonians. On the EP side, in parallel, precisely such a shift makes the EP singularities accessible. Indeed, they are suddenly found as lying on boundaries of the relevant parametric domains admitting the standard probabilistic interpretation of the systems.
Our detailed scrutiny of this territory will proceed step by step. In section 2 an introductory review of the formalism will be provided summarizing the basic features of the PT −emphasizing quantum mechanics of closed systems. Subsequently, in section 3 we will turn attention to a few most relevant features of the role of EPs in such a theory or theories. In section 4 we will finally combine the preceding material in a climax outlining the bridges between SUSY and EPs. We will show, in essence, that the recent transfer of the applications of SUSY to quasi-Hermitian quantum systems extends the class of eligible models while still finding its natural limitations of applicability precisely at the same natural EP boundaries.
A few complementary aspects of our message will be also discussed in section 5. We will emphasize there that there exists a close relationship between the EP-related properties of the SUSY and non-SUSY quantum systems. All of these systems share the limitations of their unitarity and observability (i.e., of their probabilistic interpretation) strictly along the mathematical parameter-range boundaries formed by the EP singularities.

The concept of hiddenly Hermitian Hamiltonians
In the present paper our perception of quantum physics will be restricted to the theory in which a non-Hermitian but hiddenly Hermitian Hamiltonian H generates a strictly unitary evolution of quantum system in question (see the resolution of an apparent paradox as given in the early review of this approach to quantum mechanics in [11]).

Non-Hermitian physics: closed versus open systems
In a concise introductory comment let us point out that the quantum systems under our present consideration are usually called "closed systems". Carefully, one has to separate the study of these models (with real energies) from the descriptions of other, manifestly non-unitary quantum systems (with complex energies) called "open". Unfortunately, very similar language is often used in the literature on both the closed and open quantum dynamics. In a word of warning against possible misunderstandings let us recall, for illustration, the existence of potentially misleading overlaps in terminology. Often, people indiscriminately speak about "non-Hermitian quantum mechanics" having in mind either the resonances in unstable, open quantum systems [15,16] or the safely stable bound states in a closed, unitary setup [12]. Needless to add that the specific trademark logo of PT −symmetry born in the context of quantum mechanics [17,18] and quantum field theory [19] became quickly popular also in the non-quantum world [13]. In optics, for example, the idea proved particularly successful [20]. Many research teams using Maxwell equations revealed that the propagation of light acquires remarkable properties when studied in a non-conservative medium characterized by PT −symmetry redefined, ad hoc, as a balance between loss and gain. Thirdly, the same name "PT −symmetry" extended its implementations also to non-linear systems and equations, with the deepest relevance and applicability even in nanotechnologies [8].
Out of such a broad area of physics, our present attention will only be paid to the quantum closed systems.

Triplets of Hilbert spaces: A key to the unitarity paradox
The famous Stone's theorem [21] implies that in a preselected Hilbert space (denoted, say, by a dedicated symbol K), the Schrödinger equation cannot describe a unitary (i.e., norm-preserving) evolution of a state |ψ(t) whenever the Hamiltonian happens to be non-Hermitian, H = H † . Fortunately, in a way dating back to Dyson's studies of ferromagnetism [22] it became clear that in many realistic quantum models the non-unitarity of the evolution in K may be declared, under certain conditions, irrelevant. Indeed, in many practical applications, it may make sense to treat H as a mere auxiliary isospectral avatar of a "true and correct" Hamiltonian which is "acceptable" (i.e., self-adjoint, h = h † ) in another, true and correct, physical Hilbert space of states (to be denoted here by another dedicated symbol L). In other words, the Dyson-attributed preconditioning (2) changes the picture while making the theory compatible with the Stone's theorem as well as with the conventional-textbook closed-system unitary-evolution requirements. Dyson's numerous followers (cf., e.g., review [11]) were pragmatic practitioners who usually started their considerations from the principle of correspondence. Having specified a realistic (typically, a fermionic atomic-nucleus-model) Hamiltonian h acting in L they needed to evaluate its low-lying spectrum. In such a (typically, variational) setting they revealed that the calculations may be perceivably simplified, via preconditioning (2), after an educated ad hoc guess of a suitable invertible Hermiticity-breaking map Ω −1 : L → K. This map represented, typically, correlations between fermions living in L in the language of certain effective bosons living in K. In practice, such a strategy really accelerated the calculations yielding a perceivably simplified picture of the underlying physical reality.
The Dyson-recommended preconditioning h → H is, almost without exceptions, motivated by the simplification of the solution of Schrödinger Eq. (1). Naturally, the price to pay is not entirely negligible because the non-Hermiticity of H (which must still be assumed diagonalizable [9]) implies that its right and left eigenvectors are not just mutual Hermitian conjugates. For this reason (see also [23]) it makes sense to complement the bra-vector conjugates n| of the conventional right eigenkets |n by the left-eigenbras of H (denoted as eigenbrabras n|) with, naturally, different conjugates |n = |n (that's why we amended the notation).
For the same reason the full and explicit description of a pure time-dependent state of the systems requires not only the knowledge of the ket-state solutions |ψ(t) of the evolution equation (1) but also the complementary, non-Hermiticity-related knowledge of the brabra (or ketket) solutions of the following complementary evolution equation in K, The two Schrödinger equations are interrelated by Eq. (2) and by the fact that h = h † in L. Thus, in K we can speak about quasi-Hermiticty of H [10] and write In the stationary cases this observation makes the second Schrödinger equation redundant because we may set, simply, |ψ(t) = Θ |ψ(t) . Moreover, after we introduce such a biorthonormal basis {|n , |n } in K which diagonalizes H, we have relations We become able to define the third Hilbert space H when changing the inner product (•, •) between any two elements ψ a and ψ b in K. Thus, once we define we may conclude that there are no reasons for calling H "non-Hermitian". It is manifestly non-Hermitian in K but this Hilbert space is unphysical [23]. A precise and detailed mathematical formulation of this conclusion was given in survey [24]: Lemma 1 [24] The observable predictions concerning a closed quantum system and using the space-Hamiltonian doublets (H, H) or (L, h) are equivalent.
In the theories of closed quantum systems based on equivalences (H, H) ≡ (L, h) the terminology is still unsettled. The underlying innovative formulation of quantum mechanics is often called quasi-Hermitian [11] or, more or less equivalently, PT −symmetric [12] alias three-Hilbert-space (3HS, [23]) alias crypto-Hermitian [25] alias pseudo-Hermitian [9]. Mathematicians, in contrast, seem to speak just about non-self-adjoint operators without adjectives [26]. In what follows, the acronym 3HS will be mostly employed.

Non-uniqueness of Hilbert space H
In practice people usually make use of just one of three alternative consequences of the equivalence in Lemma 1. The first one is trivial: whenever the evaluation of predictions remains feasible before preconditioning (2) (i.e., in the initial representation (L, h)), there is no reason to leave the comfortable constructive single-space textbook recipes. The selection of one of the other two options, i.e., the description of reality in (H, H)−representation is usually motivated either by the prohibitively complicated aspects of Hamiltonian h, or by a technical unfriendliness of Hilbert space L. The first version of the "non-Hermitian", (H, H)−based option was made popular, many years ago, by Dyson [22]. Its key technical ingredient may be seen in accessibility of an initial knowledge of mapping Ω. Still, even with this knowledge, the not too elementary transition from the traditional physical Hilbert space L to its ad hoc auxiliary partner K (in which the Hamiltonian becomes non-Hermitian) must really very strongly be motivated by the associated emergent technical simplifications [11].
In our present paper we intend to prefer the second "non-Hermitian" option which, in essence, means that one wants to work in a maximally user-friendly Hilbert space K from the very beginning. The analysis then starts from a presumably physical diagonalizable Hamiltonian candidate H and from its diagonalization leading to formulae (5).
It is necessary to add that when we decided to set, in our preceding considerations, |ψ(t) = Θ |ψ(t) , we simplified our analysis only at an expense of ignoring the fact that during the construction of the basis we did not take into considerations that either the biorthonormality and completeness relations are both invariant under a simultaneous rescaling of |n → |n /̺ n and |n → |n × ̺ n . As we already emphasized in [27], this is a well known fact [11] which implies that the reconstruction of the physical Hilbert space H (i.e., the formulation of the physical contents of the theory) is in fact ambiguous and, whenever one wishes to suppress this ambiguity, the recipe requires more information.
The "missing" information may find its origin in physics (one can introduce some other observables [28], cf. also reviews [11,12]). Some additional mathematical requirements may also help (in [29], for example, we proposed to minimize the anisotropy of the geometry in the physical Hilbert space). Still, in both contexts the essence of the problem is that the innocent-looking rescaling redefines the metric. Thus, just one of its special cases was used in Eq. (6). The general Hamiltonian-compatible form of the physical metric is in fact non-unique and multiparametric, This implies that also the inner product (•, •) and space H become variable, ambiguous and ̺−dependent.
For practical purposes, there exist multiple methods of removal of such an unpleasant nonuniqueness of the theory. Their samples may be found, e.g., in [11,12] or in [29].

The concept of exceptional points
One of the main distinguishing features of the parameter-dependent non-Hermitian Hamiltonians H(λ) is that even in the absence of any symmetry (including, in a very prominent place, supersymmetry) one can still encounter a degeneracy of the energy levels in their spectra [30]. According to Kato [6], such a degeneracy proves to be of a fundamental importance in perturbation theory. For this reason he conjectured to call the underlying values of parameters exceptional points (EPs). In a way depending on the mathematical context these values λ (EP ) were characterized either by their analytic-function background (i.e., by their connection with the branch point singularities, see the next paragraph), or by the role played in spectral theory. The latter area of mathematical applications also led to the Kato's best known definition of the EPs as the points at which the algebraic multiplicity of an energy eigenvalue becomes different from its geometric multiplicity (see p. 62 in loc. cit.).
For the practical users of the idea it is sufficient to know that at the instant of the EP degeneracy of the energies one also observes a parallelization of some of the eigenvectors. Still, in the context of physics it took time before people imagined that the originally not too emphasized theoretical as well as experimental accessibility of these degeneracies could open a broad new area of many new EP-related phenomena [3, 7, 31].

The birth of the concept in perturbation theory
From the strictly historical point of view the current interest in exceptional points grew from several independent sources. All of them are, directly or indirectly, related to the abstract mathematical concept of an analytic function F(λ). One should remember that F(λ) is, in general, multivalued at its (complex) variable λ so that, strictly speaking, it is not a function but rather a collection of several single-valued functions which are represented, up to a suitable radius of convergence λ max , by the Taylor series. For a quick clarification of this concept it is sufficient to recall the square-root function S(λ) which is, in the complex plane of λ, two-valued. Its two branches S ± (λ) = ± √ λ are already single-valued functions which coincide at λ = 0. They may be defined, locally, by their respective Taylor series T ± (λ). Thus, e.g., at a small complex ξ in λ = λ(ξ) = 1/2 + ξ these two Taylor series are infinite, and they both have the radius of convergence equal to ξ max = 1/2. In Kato's book [6] the Taylor series represent the roots z n = z n (λ) of secular equations P(z) = det[H(λ) − z I] = 0 where the operator H(λ) (or, in simpler cases, an N by N matrix H (N ) (λ)) plays the role of quantum Hamiltonian. The explicit Taylor-series constructions of its boundstate-energy roots z n (λ) = E n (λ) = E n (0) + λ E ′ n (0) + . . . are then deduced and shown to be of paramount importance in practice. What is essential is the knowledge of the radius of convergence λ max . A key to the answer is that the collection of roots E n (λ) often admits a unifying analyticfunction reinterpretation E(λ). Its separate single-valued sheets represent energy levels E n (λ). The radius of convergence λ max acquires an elementary meaning of the distance of the origin from the nearest EP singularity.

EPs in two by two effective Hamiltonians
In recent review paper [32] the authors emphasized that the study of certain complicated integrable models (IM) in 1+1 dimensions may often profit from the so called ODE/IM correspondence. The acronym ODE abbreviates here the ordinary differential equations, and the benefits are not only bidirectional but also deeply relevant in the present context. This was demonstrated, earlier, in [33] where the same correspondence formed a background of the famous proof of reality of the boundstate energy spectra of certain ordinary differential three-parametric quantum Hamiltonians H which were found non-Hermitian in K = L 2 (S) (here the symbol S stands for a suitable complex contour).
In a brief addendum [34] the authors of the proof paid more attention to the structure of the EP boundaries ∂D of the domain of parameters at which the spectrum remains real. They pointed out that in a small EP vicinity the passage of parameters through the boundary can be visualized via an effective two-by-two Hamiltonian This was finally the simplification which helped them to clarify some aspects of the reality proof in [33] including, in particular, the fact that in their model one can also detect a certain form of hidden supersymmetry. All of these observations contributed to the motivation of our present study. First of all we imagined that a replacement of a complicated Hamiltonian by its "effective" two-by-two-matrix simulation offers a useful and universal mathematical tool of a qualitative description of the spectrum near its EP singularity. Secondly, the efficiency of the use of the effective N by N Hamiltonians seems to survive even a transition to the more complicated EP singularities, i.e., far beyond the class of their comparatively elementary ODE realizations [35]. Thirdly, even if one stays just with the most conventional ODE models, the effective-matrix simplifications may be expected to open a viable way towards a semi-explicit construction of the physical Hilbert space H, i.e., of the last -and often omitted -step of the theory.
In the special case of ordinary differential Schrödinger equations, virtually all of the other available methods of such a construction of H and Θ seem to be prohibitively difficult [36]. In the approximation by the two by two matrix simulation (8), the latter last step remains elementary. Indeed, once we insert matrix H (2) (δ) of Eq. (8) together with a real and symmetric general matrix ansatz for the corresponding Hilbert-space metric Θ (2) in the hidden Hermiticity condition (4) we obtain, with obvious proof, an exhaustive answer at any positive δ = η 2 and positive η, Lemma 2 The necessary positivity of metric (9) is guaranteed whenever its free real parameter b is such that b 2 < 1.

EPs in crypto-Hermitian models with local potentials
In the physics-oriented literature, by far the most popular class of illustrative "non-Hermitian" Hamiltonians H with real spectra is being obtained via a judicious modification of some of the standard self-adjoint models of textbooks.

Regularized harmonic oscillator
In the above-cited papers [32,33,34] the authors studied a specific three-parametric family of non-Hermitian one-dimensional Schrödinger ODEs describing bound states in a local potential and having the standard eigenvalue-problem form H(M, A, ℓ) ψ = E ψ where and where x moved along a suitable ad hoc complex contour. Using sophisticated mathematical methods these authors contributed to the clarification of the not quite expected fact that inside a suitable domain of parameter D (with, as we mentioned above, EP boundary ∂D) the bound-state spectrum {E n } proved real and discrete even though the potential itself was not real.
These observations were preceded by our letter [37] in which we showed that in the oneparametric special case of H(1, 0, ℓ) (which appeared to be exactly solvable), many results become obtainable via entirely elementary methods. Also, the exact diagonalizability of the special case with Hamiltonian H(1, 0, ℓ) or, equivalently, the exact solvability of the following Schrödinger ODE bound state problem with α = ℓ + 1/2 ≥ 0, will prove useful, in what follows, for illustration of several nontrivial SUSY-EP correspondences.
In a preparatory step let us now briefly summarize some of the most relevant features of the bound-state solutions of Eq. (10). Firstly, in terms of the confluent hypergeometric special functions their general form is known, Secondly, these functions belong to K = L 2 (−∞, ∞) if and only if the infinite series terminates. This yields the compact formula for the wave functions They are numbered by the so called quasi-parity quantum number Q = ±1 and by another, common quantum number n = 0, 1, . . .. The power and exponential factors are accompanied here by the classical Laguerre polynomials, What is most important (cf. [38]) is that the operator Q with quasi-parity eigenvalues Q has later been reinterpreted as the operator of charge C which defines, together with conventional parity P, the most popular Hilbert-space special metric Θ = CP [9,12]. As a consequence, the values of energies are all real whenever α ≥ 0, E = E Q,n = 2 − 2Qα + 4n , Q = ±1 , n = 0, 1, 2, . . . .
Nevertheless, a warning must follow: The validity of the strictly mathematical result (12) of Ref. [37] does not imply that all of these real numbers can really be treated as bound state energies. The reasons are rather subtle, made explicit only long after the publication of loc. cit. Their essence lies in the emergence of EPs and in the ambiguity of metrics. The effect is most easily clarified via Fig. 1. In this picture the leftmost, most obvious EP degeneracy is easily seen to occur at α = α (EP ) 0 = 0. This is an instant of complexification, very often called spontaneous breakdown of PT −symmetry [12]. This is a very special EP mathematical singularity. In [39] it has been given an alternative physical meaning and name of quantum phase transition of the first kind. Schematically, the underlying mechanism of the complexification is well explained by the two by two matrix model (8) of Ref. [34]. In the limit δ → 0 this matrix degenerates to the manifestly non-diagonalizable (i.e., unobservable and phenomenologically unacceptable) Jordan-block matrix J (2) (0). It is less easily seen that all of the subsequent integers α = α (EP ) k = k = 1, 2, . . . are also exceptional points (cf. [37] for the proof). In [40] we explained their specific nature, and we proposed to call these unavoided level crossings without complexification the quantum phase transitions of the second kind. In the light of the Laguerre-polynomial solvability of the present model the quantum phase transition degeneracies of the second kind can even be given the form of exact identities L

Regularized three-body Calogero model
In [41] we re-interpreted Schrödinger Eq. (10) as a two-body special case of an integrable A−body Calogero model of the rational PT −symmetric A n series [42,43]. Although the discussion in loc. cit. did not involve the question of EPs, we are now re-attracting attention to this family of models because their exact solvability again opens the way of our understanding the presence, localization and properties of the EPs in a less elementary setting. For the sake of brevity let us only recall here the first nontrivial three-body version of the HamiltonianH Interested readers may find all of the technical details and formulae in [41]. Here we will only recall a few most relevant features of the model and add a few EP-related remarks.
which defines the spectrum. From our present point of view such a formula demonstrates, immediately, that the EP singularites may be expected to occur, in these integrable models, at any number of particles A. At the same time, an increase in the complexity of the spectrum at A = 3 (made explicit by the occurrence of an additional quantum number k) indicates that not all of the unavoided level crossings must necessarily be the EP singularities. Some of them just reflect the presence of a conventional (e.g., rotational [41]) symmetry in the corresponding partial differential Schrödinger equation. A conversion of formula (13) into a picture (sampled by Fig. 2 here) just reconfirms that in contrast to the model of the preceding subsection, some of the energy levels become degenerate due to the classical symmetries of interaction rather than due to the emergence of an EP singularity. At the higher quantum numbers this phenomenon really does occur in a way mentioned, in Fig. 2, in the right column of comments. Unfortunately, we do not have enough space here for a sufficiently exhaustive discussion of the differences or analogies between the qualitatively slightly different dynamical scenarios with A = 2 and A = 3 or A ≥ 3. In fact, one of our main reasons is that after one turns attention to the supersymmetric extensions of these models, the number of open questions starts growing rather quickly even during the study of the perceivably simpler illustrative calogerian model with A = 2.

The concept of hiddenly Hermitian supersymmetry
Even though the SUSY-based arrangement of existing elementary particles into supermultiplets [44] did not work in particle physics [45], the idea itself proved extremely successful after its transfer to quantum mechanics. It helped to classify a family of exactly solvable models [46,47]. In the words of the latter reference, "in a supersymmetric theory commutators as well as anticommutators The degeneracies and coupling-dependence of the spectrum of the non-Hermitian threeparticle Calogero model of Ref. [41]. Levels E (n,k) = 5 + 4n + 6k ± 3 √ 1 + 2g are numbered by a pair of integers n, k = 0, 1, . . .. All of the visible exceptional points g (EP ) are marked by vertical lines.
appear in the algebra of symmetry generators". For our present purposes we shall use the explicit forms of these generators forming the specific Lie superalgebra sl(1|1) (see Eq. (18) below). The same or similar ideas were also actively used during the birth of quasi-Hermitian theories [34,48,49,50]. Surprisingly enough, during or after all of these promising developments, attention remained restricted to the spectral problem, i.e., to the solution of the non-Hermitian Schrödinger equations in K. Virtually no results were obtained in the direction of construction of the corresponding physical Hilbert space of states H, i.e., of the amended operator Θ = I of the correct physical Hilbert-space metric.
The latter gap in the literature and, in particular, in the SUSY-related implementations of the 3HS theory is in fact not too surprising because the task is known to be truly difficult [9,36]. At the same time, the knowledge of the geometry of H becomes crucial whenever the metric Θ deviates from the unit operator "too much" (see, e.g., the thorough analyses of this topic in [51,52]). For this reason, at least the use of qualitative considerations and approximate methods is needed when one decides to turn attention to the dynamical regime close to an EP extreme.

Harmonic oscillator example
Let us now return to our quasi-Hermitian harmonic-oscillator model (10) and let us amend slightly the notation: In terms of the complex coordinate we shall consider the regularized quasi-Hermitian harmonic-oscillator Hamiltonians and their two shifted versions such that α = |γ| and β = |γ + 1|. In the light of the results of paragraph 3.3.1 we may conclude that all of the integer values of γ are exceptional points, {γ (EP ) } = Z. These values leave both of our sub-Hamiltonians non-diagonalizable so that they have to be excluded from our consideration as manifestly unphysical. Under this constraint, in terms of the so called superpotential we may define operators A = ∂ x + W and B = −∂ x + W and form the set of generators Interested readers may find an exhaustive description of details in [53] [cf. also the present Eq. (12) for energies].
In this setup our present intention is to pay attention to the occurrence and role of the EPs which were not properly treated in paper [53]. We will fill the gap in what follows. For such a purpose we will only need to recall, in explicit form, the existing results concerning the spectrum of the γ−dependent super-Hamiltonian G = G(γ). Its parameter γ ∈ R taken from superpotential (17) is in fact the key quantity which determines the behavior of the (super)energies. In the light of this fact let us now split the real line of γs into three subintervals.
As we indicated above, all of the integer values of γ have to be excluded from our considerations because at these values we have γ = γ (EP ) ∈ Z so that our quantum model loses any acceptable probabilistic interpretation [9]. At the half-integer values of γ = −K − 1/2 the equidistance now characterizes the non-anomalous part of the doubly degenerate spectrum, i.e., the spectrum without the lowermost (2K − 1)−plet of levels (counting their SUSY-related degeneracy).

Corollary 6
In any preceding crypto-Hermitian harmonic-oscillator realization of SUSY algebra (18) the range of the admissible parameters γ always remains restricted to a finite interval of unit length, with the EP boundaries given by the two neighboring integers.

Regularized Hermitian limit
In 1984 Jevicki and Rodrigues [54] pointed out that formalism of supersymmetric quantum mechanics breaks down for many phenomenologically interesting potentials with singularities (i.e., e.g., with the repulsive centrifugal-like barriers). In [53] we showed that what could help is a complex-shift regularization, and we asked what happens with the SUSY structures in the limit c → 0 + . Our considerations were illustrated by the regularized harmonic oscillator of Eq. (10). Thanks to the exact solvability of this model we were able to summarize our observations in Table  Nr. 3 of Ref. [53] and in the appended spectrum-explaining comments.
From the present point of view the latter results deserve a comment. The point is that the centrifugal singularity remains regularized at any non-vanishing shift-parameter c > 0. In loc. cit. we were interested also in the behavior of the system in the limit c → 0 + . In this limit the regularization is removed so that the norm of at least some of the c > 0 solutions diverges. In loc. cit. we reviewed the situation and we managed to show that some of the solutions still remain acceptable by obeying the emergent standard boundary conditions in the origin. The resulting, drastically restricted family of normalizable c = 0 bound states is displayed in Fig. 3.

SUSY and local potentials
Whenever the local potentials are being involved (which seems to be one of the characteristic aspects of many successful applications of various SUSY models), it will be absolutely necessary to return to the older [10] as well as newer [51,55] mathematical criticism of the very foundations of the 3HS formalism. Indeed, also in this more formal direction of research many important questions remain open at present [56].
In the field of mathematics, several interrelated challenges also survive concerning the constructions of the Hamiltonian-dependent Hilbert-space metrics Θ = Θ(H) [29]. The abstract freedom of our choice between alternative formulations of quantum mechanics (sampled by the pair (H, H) and (L, h) as mentioned in Lemma 1) is of a very restricted practical use. In applications, such a parameter Figure 3: Degeneracies of the c → 0 spectra of the quasi-Hermitian and supersymmetric harmonic oscillator (10). In the conventional regular γ = −1/2 special case the supersymmetric Hamiltonian becomes self-adjoint. Only four exceptional points γ ∈ {−2, −1, 0, 1} remain phenomenologically relevant since the sets of bound states become incomplete out of the innermost admissible interval of γ ∈ (−1, 0). choice is usually more or less unique, determined by the criteria of feasibility of the constructions and calculations. This is a pragmatic attitude which is being accepted in virtually all of the applied hidden-Hermiticity settings. A decisive factor leading to the choice of a specific representation (H, H) is usually found in the user-friendly nature of the Schrödinger equations when solved in K. In this sense the field of SUSY quantum mechanics seems particularly suitable and promising fur future reseaarch, especially due to its deep algebraic as well as analytic background.

Analyticity and EPs
In the conventional quantum physics using self-adjoint Hamilonians the EP singularities are often perceived as interesting and, perhaps, mathematically useful but not too phenomenologically important complex-plane points of coincidence of an n = m pair of the analytically continued bound state energies, E n (λ (EP ) ) = E m (λ (EP ) ) [35,57]. In a historical perspective this encouraged the study of operators of observables (i.e., most often, Hamiltonians) which were analytic functions of their parameters, H(λ) = H(0) + λ H ′ (0) + . . .. Such an assumption proved to be a natural source of overlaps between perturbation-theory mathematics and the PT-symmetry-related physics [58,59]. Typically, this helped to clarify the connections between the time-independent Schrödinger equations H(λ)ψ n (λ) = E n (λ)ψ n (λ) and their solutions based on the power-series Rayleigh-Schrödinger ansatzs ψ n (λ) = ψ n (0) + λ ψ ′ n (0) + . . . and E n (λ) = E n (0) + λ E (1) n + . . . . In a way discussed by Kato [6], at least a subset of all of the (real or complex) energy eigenvalues E n (λ) may be then interpreted as evaluations of a single analytic, multisheeted function E(λ).
Kato also emphasized that a nice explicit illustration of the related mathematics may be then provided by certain ad hoc finite-dimensional N by N matrix toy-model Hamiltonians H (N ) (λ). In the present context of the SUSY -EP correspondence this could inspire a further study of various quasi-Hermitian SUSY Hamiltonians and, in particular, their unphysical EP limits (plus their small perturbations -see, e.g., a few related remarks in [31]) in their canonical block-diagonal representations. For example, in a small vicinity of unperturbed α = α where the blocks would be just the Jordan-block matrices J (2) (η) obtained form the matrix of Eq. (8) in its non-diagonalizable extreme, J (2) (η) = lim δ→0 H (2) (δ) + η I (2) . In the future applications involving some less elementary and, in particular, some supersymmetric models one may expect to obtain various less schematic and mathematically more interesting EP-related matrix structures. They could serve as approximations of realistic Hamiltonians along the lines studied, e.g., in [52].

Outlook
In our paper we skipped all outlines of physics behind the 3HS theory, and also of the multiple early contacts of this theory with SUSY (an extensive information on this topic may be found, e.g., in [60]). In any forthcoming phenomenology-oriented analysis more attention should certainly be paid, therefore, to the interlaced SUSY and EP aspects of the constructions of the physical Hilbert spaces of states H as well as to the role played by some other, non-Hamiltonian operators of the hiddenly Hermitian observables.
During the future developments of the field of the SUSY-EP contacts one should certainly return, last but not least, also to the realistic-physics inspiration provided by the Dyson-inspired studies of non-SUSY interacting boson Hamiltonians H as well as by the recently revealed Dysonunrelated hiddenly Hermitian but non-SUSY parallel developments involving, for example, the clusters-coupling Hamiltonians H [61].