Anisotropy and asymptotic degeneracy of the physical-Hilbert-space inner-product metrics in an exactly solvable crypto-unitary quantum model

In quantum mechanics (formulated, say, in Schr\"{o}dinger picture) only the knowledge of a complete set of observables $\Lambda_j$ enables us to declare the related physical inner product (i.e., the Hilbert-space metric $\Theta$ such that $\Lambda_j^\dagger \Theta=\Theta\,\Lambda_j$, i.e., such that $\Theta=\Theta(\Lambda_j)$) unique. In many applications people simplify the model and consider just a single input observable (mostly an energy-representing Hamiltonian $\Lambda_1=H$) and pick up, out of all of the eligible metrics $\Theta=\Theta(H)$, just the simplest candidate (typically, in the case of the special self-adjoint input $H$ we virtually always work with trivial $\Theta=I$). As long as this forces us to admit only the self-adjoint forms of any other input observable $\Lambda_j$, the scope of the theory is, without any truly meaningful phenomenological reason, restricted. In our present paper we describe a strictly non-numerical $N$ by $N$ matrix model in which such a restriction is replaced by another, phenomenologically non-equivalent restriction in which $\Theta \neq I$ and in which the system reaches a collapse (i.e., a loss-of-bservability catastrophe) via unitary evolution.


Introduction
In many standard and routine applications of quantum theory the evolution in time is prescribed, in the so called Schrödinger representation [1], by Schrödinger equation where the state vector belongs to a physical Hilbert space of conventional textbooks, |ψ(τ ) ≻ ∈ H (T ) .The Hamiltonian is usually assumed time-independent and self-adjoint in H (T ) .Often, the theory is realized in the physical and, at the same time, userfriendly special space H (T ) = L 2 (R d ) of square-integrable coordinate-dependent functions ψ(x, τ ) =≺ x|ψ(τ ) ≻ in d dimensions.
Under these assumptions the evolution is unitary [2] and Eq. ( 1) is formally solvable, The practical construction of the wave functions usually proceeds via an approximate or exact diagonalization of h = h † [3].The description of the evolution remains equally routine for the time-dependent Hamiltonians h = h(τ ).One may also move from the primary Schrödinger representation to its equivalent Heisenberg-representation alternative.Via a suitable unitary operator one then obtains the Heisenberg-representation wave functions which are required not to vary with time [1].
A more challenging theoretical as well as conceptual scenario emerges when the Heisenbergrepresentation-inspired preconditioning of the wave-function ket-vector is chosen invertible but non-unitary [4], i.e., such that Then, product Θ(τ ) can be perceived as playing the role of a correct Hilbert-space metric in an "amended" physical Hilbert space H (A) such that the ket-vectors |ψ ∈ H (A) are simpler in comparison with their more conventional textbook avatars |ψ ≻ ∈ H (T ) .In most of the reviews of the state of art (cf., e.g., [5,6,7]) the authors are trying to cover the whole new terrain of the theory on a rather abstract level.At the same time, the simplification of the picture caused by the change of paradigm is usually mentioned just marginally, as an assumption or a tacit wish rather than as a rather difficult necessary condition of a consequent practical and constructive implementation of the formalism.
In our present paper we intend to pay more attention to the postulates of solvability and simplification.Indeed, these features of the quantum models of interest represent one of the not often emphasized keys to applicability of the whole non-unitary-preconditioning idea behind the Dyson-inspired mapping (3).
The presentation of our considerations and results will start in section 2 where we will summarize a few basic concepts forming the theory.In subsequent section 3 we will point out that in the context of rigorous mathematics the theory itself is still in the stage of formal development, characterized by the existence of a large number of open mathematical questions and challenges (cf., e.g., [8,9,10]).In this sense we decided to circumvent some of these challenges in the spirit of the words of warning in reviews [8,11].Thus, we will just consider a family of sufficiently innocent-looking benchmark models living in Hilbert spaces of an arbitrary finite dimension N.
One of the phenomenologically most relevant benefits of such a choice of models will be discussed in section 4. Its essence will be emphasized to lie in the possibility of using the well known mathematical ambiguity and flexibility of the inner-product metric Θ for the purposes of description of the quantum systems in an arbitrarily small vicinity of their singularities representing certain forms of a quantum catastrophe.
In section 5 we will modify the paradigm and extend the dynamical framework of Schrödinger representation which is inherently stationary.We will turn attention to a more explicit study of the time-dependent aspects of our class of benchmark (and exactly solvable) models.In this section we will emphasize that the dominant merit of our models lies in the closed-form availability of non-stationary metrics Θ = Θ(τ ).
The details of the construction are made explicit in sections 6 and 7 in which we present a mathematical core of our present message.A basic mathematical characteristic of our class of models will be shown to lie in the smoothness of the time-dependence of the innerproduct metrics Θ = Θ(τ ) and, first of all, in the existence of these operators for the times covering the whole interval connecting, in one extreme, the Hermitian quantum mechanics (characterized by the trivial and fully isotropic metric and reached, in our units, at τ = 0) with the other extreme of a "quantum catastrophic" alias "phase-transition" alias "fully degenerate" collapse of the system in the τ → 1 limit in which the inner product metric asymptotically and ultimately degenerates and ceases to exist.
A few concluding remarks will be added in sections 8 and 9.

An outline of theory
The conceptual consistency of the non-unitary Dyson's mapping (3) is based on the requirement of equivalence between the evaluations of the old and new inner product, The main reason why the non-unitarity Ω(τ ) = Ω † (τ ) in Eq. ( 3) is challenging is that the survival of the requirement of equivalence of physics in H (T ) and H (A) leads to the apparently counterintuitive definition (5) of the inner product in H (A) .Subsequently, it is fairly difficult to resist the temptation of introducing another, third, user-friendlier Hilbert space H (F ) in place of H (A) .In it, one re-accepts the manifestly unphysical but simplerto-use metric Θ (F ) = I which only has to be remembered as mathematically useful and preferable even though manifestly unphysical.
In the case of the simplest, manifestly time-independent non-unitary mappings Ω the trick (5) and transition to the "three-Hilbert-space" (THS) representation of a given quantum system proved particularly rewarding in applications, say, in nuclear physics [11,12].The isospectrality of the mappings of the Hamiltonians as induced by Eq. ( 3) has been used there to facilitate the practical variational estimates of the bound state energies of certain heavy nuclei.Later on, the THS formalism also found applications in relativistic quantum field theory.Emphasis has been redirected to the study of systems exhibiting the parity times timereversal symmetry alias PT −symmetry of the Hamiltonian (cf., e.g., the dedicated reviews [5,6] for an extensive information and detailed discussion).
Quickly, the growth of the scope of the theory has been also noticed and accepted in the other parts of physics like, say, experimental optics [13].Various unexpected consequences of the generalization h → H have been found inspiring, especially when the researchers managed to keep the evolution-generator PT −symmetric.This led to a new perception of Maxwell equations (in the so called paraxial approximation) and to the experiments using metamaterials with anomalous refraction indices [14,15,16]).
A purely quantum theoretical as well as phenomenological appeal of the THS approach reemerges when one opens the Pandora's box of time-dependent problems [17] - [27].First of all, it is necessary to imagine that the second Hilbert space becomes time-dependent in a way mediated and carried by the time-dependence of its metric Θ = Θ(τ ) [23,25].
Thus, one has H (A) = H (A) (τ ) and one must replace the time-evolution Schrödinger Eq. ( 1) valid in H (T ) by its manifestly non-Hermitian analogue [ 17,23].This version of evolution equation may still be considered and solved in the unphysical but user-friendlier Hilbert space H (F ) .In this case it is only necessary to keep in mind that the sophisticated non-Hermitian generator of evolution must be defined as composed of the observable Hamiltonian and of another operator called quantum Coriolis force [27].Marginally, it may be added that both of these components of the observable "instant energy" Hamiltonian H(τ ) have, in general, complex spectra [28].
Obviously, a consistent version of the formalism requires a cancelation of the non-Hermiticities as carried by G(τ ) and Σ(τ ).In a way explained in brief note [20] the latter cancelation is still absent in the systems with trivial G(τ ) = 0.In papers, [20,22], incidentally, such an option (simplifying Schrödingerian Eq. ( 7) and reminding us of Heisenberg representation) has been extended to cover also the non-vanishing but still simplified, viz., time-independent constant-operator Schrödingerian generators G(τ ) = G(0) = 0. Nevertheless, in the fully general version of the formalism one cannot rely upon similar simplifications.
In particular, the changes of the physical inner product in H (A) (τ ) need not be slow.Hence, the Coriolis-force operator Σ(τ ) treated as a difference between H(τ ) and G(τ ) need not be small, either.This might make the adiabatic approximation more or less useless [29,30].At the same time, we are often forced to assume the validity of the adiabatic approximation for practical purposes.This is the situation in which one needs a methodical guidance mediated, typically, by the exactly solvable examples.Via their deeper analysis one can identify the dynamical regimes in which Σ(τ ) can be kept small.
A family of models characterized by an explicit, closed-form knowledge of the relevant operators will be introduced and described in what follows, therefore.

Benchmark model
One of the most immediate consequences of relations ( 5) and ( 6) is that the self-adjointness of h(τ ) can equivalently be re-expressed as the metric-dependent quasi-Hermiticity [8] of An analogous relation is also required to be satisfied by any other candidate Λ(τ ) for an observable.
In the light of the Dieudonné's critical analysis [8] the authors of review [11] pointed out and strongly recommended that all of the eligible candidates for an observable (including, naturally, also the energy-representing Hamiltonian) should be, preferably, represented by operators which are, in the friendly Hilbert space H (F ) of mathematical preference, bounded.
In our present paper we followed the recommendation.

Bounded-operator Hamiltonians
In the context of prevailing model-building practice, the constraint of boundedness appeared over-restrictive [31,32,33].Indeed, multiple popular elementary forms of Hamiltonians have the form of a differential operator which is unbounded [34,35].In this light and in a way recalled at the beginning of section 8 below we decided to accept the constraint and to replace, in one of our related papers [36], the most common harmonic-oscillator ordinary-differential Hamiltonian by its truncated and shifted, N−dimensional diagonalmatrix equidistant-spectrum analogue In parallel, the methodical and pedagogical role of the most popular anharmonic-oscillatorlike Hermiticity-violating interactions [37,38] was transferred to the off-diagonal elements of certain real and, say, tridiagonal non-Hermitian N by N multiparametric matrices This enabled us to simplify the proofs of the required reality of the bound-state spectra.
For models (11) we managed to reduce these proofs to a virtually elementary spectralcontinuity or spectral-inertia argument, applicable at the not too large couplings g j at least.After an additional up-down symmetrization of the above matrix, i.e., after the choice of g N −1 = g 1 and g N −2 = g 2 , etc, we arrived at the final form of our benchmark toy-model Hamiltonian After such a choice of the class of models we also managed to parallel the phenomenologically relevant parity times time reversal symmetry of multiple common differential-operator toy-model Hamiltonians by a formally analogous PT −symmetry as imposed upon our matrices (12).This build-up of analogy merely required the specification of T as the transposition plus sign reversal.The parity-simulating indefinite square root of the unit matrix appeared then represented by an antidiagonal N by N matrix P with non-vanishing elements P j,N −j+1 = 1, j = 1, 2, . . ., N.
In our subsequent paper [39] we turned attention to one of the methodically most welcome features of the PT −symmetric and J = [N/2]−parametric benchmark Hamiltonian (12), viz., to the availability of the amazingly elementary geometric form of boundary ∂D of the J−dimensional compact domain D of the real parameters g j for which the spectrum of H (N ) (P T ) remains real.For our models (12) this boundary (or, in the language of physics, the horizon of the bound-state stability of the system [40]) has been shown to acquire, at any matrix dimension N, the same generic geometric form of surface of a smoothly deformed hypercube with protruded edges and vertices (cf.also [41] for more details).

Fall in instability
Initially, the family of our present toy models (12) has been developed with the purpose of having a tractable sample of a quantum analogue of a classical concept of an evolution singularity called, in the Thom's popular terminology [42,43], a "catastrophe".This aim of study has been made explicit in our paper [44].In place of our present variable τ measuring the time during the fall of the system into its degenerate singularity we used a different variable λ = 1 − τ 2 , in terms of which some of the formulae appeared simpler.
Indeed, we revealed that there exists a certain specific λ−parametrization of the couplings g j = g j (λ) in (12) such that • the two by two matrix H (2) (P T ) (λ) appears useful as a benchmark model of an energybifurcation scenario in which i.e., in which the spectrum is real iff λ ≥ λ 0 = 0 and in which the whole spectrum becomes completely degenerate iff λ = 0 while it finally gets purely imaginary iff λ < 0; • the three by three matrix H (P T ) (λ) has been found to serve as a benchmark model of a new, energy-trifurcation quantum catastrophe in which Again, the spectrum proved completely degenerate iff λ = 0. Up to the exceptional, λ−independent real level E [N/2] = 0 emerging at any odd N, the rest of the spectrum was again purely real or imaginary iff λ ≥ 0 or λ < 0, respectively; • the four by four matrix H (P T ) (λ) with spectrum found then an analogous interpretation of a benchmark quantum model admitting an energy-quadrifurcation.
Analogous quantum-catastrophic (QC) features have constructively been guaranteed to hold, for a special λ−dependence of model (12), at any integer N ≥ 2 (see more details in section 4 below).
In what follows these observations will inspire and enable us to simulate, at any N, the QC history starting at a conventional Hermitian N−level oscillator Hamiltonian at τ = 0.In the opposite, "asymptotic" extreme with τ → 1 the system will be found to collapse into a complete (i.e., N−tuple) energy-level degeneracy.

QC singularity
We will see below that during the whole evolution process our model is such that its nontrivial Hilbert-space metric Θ (N ) (τ ) can be calculated in closed form.This will enable us to attribute the phenomenon of the QC collapse to a manifest interplay between the ad hoc time-dependence of the Hamiltonian and the equally ad hoc time-dependence of the Hilbert-space physical inner-product metric.

Parametrization
One of the simplest forms of the above-mentioned λ−parametrizations of the couplings in Hamiltonians ( 12) is given by the formula which was proposed in paper [39], The merit of this parametrization is that in the whole interval of λ ∈ (0, 1) (or formally even of λ ∈ (0, ∞)) it guarantees the reality of the energy spectrum.The second merit of the λ−parametrization (16) lies in the fact that the boundary value of λ = 0 strictly separates the stable dynamical quantum regime (with λ > 0 yielding the real, "observable" N−plet of bound state energies) from the half-axis of λ < 0 (for which the system ceases to be observable).
In the algebraic terminology of paper [36] and of the older literature [45,46,47] one encounters the so called Kato's exceptional point (EP) of the N−th order at λ = 0.In paper [39] we restricted our attention to the models with small λs, therefore.We found that Eq. ( 16) may be further reparametrized in terms of another time variable t which has been found to measure a recovery from the QC singularity in a broader physical multiparametric domain D, The alternative time-parameter t ∈ (0, ∞) appeared suitable for our having Hamiltonian H (N ) (P T ) more comfortably tractable near its EP alias QC singularity.

Time-dependent metric
At small ts the evolution can be interpreted as the motion of system away from QC, towards a stable and less anisotropic dynamical regime.Parametrization (17) proves useful, first of all, at the very short times t ≪ 1 at which it effectively rescales and magnifies the interior of D in the vicinity of EP.Simultaneously, such an ad hoc change of scale did not lower the number of degrees of freedom -one could still work with as many as J alternative coupling constants G n ≥ 0.
Trivial selection of special couplings G n = 0 realizes a one-parametric, simplified but still instructive reduction of the picture of dynamics.After the hypothetical start of evolution at the λ = 0 singularity one ultimately reaches a manifestly Hermitian regime at λ = 1.In this sense, parameter λ = λ(t) ∈ (0, 1) measures the recovery.
At any fixed value of N and for any suitable Hamiltonian the physics varies with inner product (5).The reconstruction of all of the eligible metrics Θ may be found summarized in our dedicated work [48].It can be summarized as starting from Schrödinger equation where we replaced Hamiltonian H (N ) (P T ) by its Hermitian conjugate in H (F ) .The complete solution of the new equation opens the way towards the reconstruction of any metric from its spectral representation All of the parameters κ n > 0 are freely variable.This is an ambiguity which reflects the absence of an exhaustive information about the physics behind the quantum system in question (cf., e.g., Ref. [11] for explanation).In Refs.[48] and [49] we discussed the general recipe from a more formal point of view by which multiple suitable metrics Θ may always be assigned to a given Hamiltonian H via Eq.(19).We emphasized there that the construction is always ambiguous.Now, it is worth adding that the sufficiency of the solution of the auxiliary conjugate Schrödinger equation (18) has to be perceived as a another, serendipitious merit of the models living in finite dimensions N < ∞.

Two by two example
Let us pick up N = 2 and study formula (19) in more detail.Firstly, let us change the variables, t → r = r(t) = √ t > 0, yielding In terms of the pair of abbreviations u = √ 1 − r and v = √ 1 + r we may then calculate the maximal and minimal eigenvalues E + = r and E − = −r of (20) as well as the related respective real eigenvectors where the superscript T denotes transposition.
It remains for us to insert vectors ( 21) and ( 22) in the spectral expansion of the metric.Once we fix an inessential overall factor and denote κ + = sin α and κ − = cos α with 0 < α < π/2 we get the general N = 2 metric-operator matrix Its elementary form facilitates the direct determination of its eigenvalues, One easily verifies that the requirement of the necessary positivity of these eigenvalues is trivially satisfied at any square-root-time r = √ t such that 0 < r < 1.We may conclude that the standard probabilistic interpretation of our time-dependent N = 2 THS QC quantum model is determined not only by its one-parametric Hamiltonian (20) but also by the specification of the concrete value of variable α.Via Eq. ( 23) this choice selects one of the eligible inner products (5).This makes the Hilbert space of states fully defined and unique, with an asymmetry alias anisotropy of its geometry measurable simply by the difference 1 − r 2 sin 2 2α.

Anisotropy
Any energy-representing input Hamiltonian H(τ ) admits many alternative, nonequivalent predictions of the results of measurements [11].Formally, this is caused by the existence of multiple Hamiltonian-independent parameters as sampled by α in Eq. ( 23) at N = 2.The ambiguity must be removed via additional, physics-based constraints.
The situation becomes slightly different near the QC degeneracy because the EP-related confluence of the energy levels already implies that a consistent picture of reality and, in particular, the necessary regular inner-product metric ceases to be positive definite.
A universal remedy does not exist because a small increase/decrease of time may cause a large change of the metric in general.In what follows we will exclude such a formally admissible non-perturbative behavior of Θ(τ ) as unphysical.

Metrics as functions of time τ
In our older paper [44] we constructed the well-behaved, extrapolation-friendly metrics up to N = 3 or, in an implicit form (19), up to N = 5.The choice of parameters κ was dictated by a half-intuitive requirement of simplicity.Vague as such a recipe might have been, it found an independent support in a similarity of formulae at several Hilbert-space dimensions.The success of such a choice of parameters motivated also our present study.
Our present strategy will be based on a change of philosophy.In place of studying just a vicinity of QC using a small EP-unfolding time t (alias τ 1) we will search for an explicit description of the QC-emergence process in the whole interval of τ ∈ (0, 1).Thus, the time variable τ will replace the not too suitable time-like parameter λ(t) of Eq. ( 17)).Running in opposite direction, i.e., from the τ = 0 instant (at which our Hamiltonian is diagonal) to the QC limit of τ → 1 (in which our Hamiltonian becomes non-diagonalizable and merely Jordan-block representable).
The choice of the new "time of collapse" variable τ = √ 1 − λ will enable us to treat the fall of our system into its singularity as a process which started long before the catastrophe.In the final QC limit τ → 1 the time-dependent vectors |ψ n of Eq. ( 18) can be then perceived as getting mutually parallel.The operator (19) itself degenerates to a singular but particularly elementary matrix of rank one of course.
A disadvantage of the replacement of t or λ(t) by τ might be that the construction of the metric appeared easier at small λ.Fortunately, near the opposite extreme of τ = 0 the transition t → τ simplifies the Hamiltonians themselves, One is led to a replacement of expansion ( 19) by a less subtle, brute-force technique, rendered possible by the tridiagonal-matrix form of Hamiltonian.Such an idea proved productive, indeed.
Theorem 1 The metrics Θ (N ) (τ ) compatible with Hamiltonians (25) may be sought in the finite-sum form with sparse-matrix coefficients Proof follows from the observation that at any k = 1, 2, . . ., N, the set of all of the nonvanishing elements of matrix M (N ) (k) may be compressed and arranged into an auxiliary, k by (N − k + 1)−dimensional array with α 11 (k) = M (N ) 1k (k), etc.This reduces the proof to the inspection of the set of N 2 Dieudonnés linear algebraic compatibility relations written in the matrix form not all of which are independent (cf.Ref. [48] for details).

Onset of the process of degeneracy
From the point of view of potential applications of formula (26) it is important that at the beginning of the fall into instability (i.e., at the far-from-QC instant τ = 0) the Hamiltonians (25) will all coincide with the respective truncated and diagonal (i.e., Hermitian) harmonic-oscillator-like matrices.In this picture the spectrum of energies E (N ) n (τ ) will remain real but shrinking with the growth of the innovated time τ .In the limit τ → 1, i.e., at the very end of the fall of the system into QC singularity, the spectrum becomes completely degenerate, E (N ) n (1) = 0, n = 0, 1, , . . ., N − 1.In the latter limit the Hamiltonian (i.e., matrix H (N ) (1)) ceases to be diagonalizable and loses its standard physical tractability and interpretation.Vice versa, from the point of view of physics the description of the evolution as generated by H (N ) (τ ) will change at τ = 1, requiring an introduction of some new degrees of freedom beyond this instant, i.e., at τ > 1.The study of such a discontinuous switch to a new form of Hamiltonian at later times lies already beyond the scope of present paper: Interested readers might consult, e.g., a dedicated study [50].Beyond the framework of quantum theory, examples of such an EP-mediated phase transition may be found, e.g., in magnetohydrodynamics [51].
Temporarily, let us now return, in the context of interpretations, to the times of recovery t or λ(t).In these variables, the motion beyond the end-of-the-interval λ(t) = 1 appears much less exotic.Typically, the energies would not feel the change at all (cf., e.g., Eqs. ( 13), ( 14) or ( 15)).Still, the values of λ = λ (outer) > 1 remain mathematically less interesting because in the "outer" interval of parameters our Hamiltonian matrices become Hermitian.Thus, it is natural to require that in the extrapolated dynamical regime the metric remains constant and trivial, Θ (outer) ≡ I.In other words, we propose to match the non-Hermitian and Hermitian dynamical regimes strictly at λ(t) = 1, with One of the consequences of such a discontinuation of the model may be seen in the subsequent most natural reinterpretation of the point of matching λ(t) = 1: The process of the QC degeneracy is expected to proceed only at λ(t) < 1.The optimal QC-related metrics should be then required continuous just at the relevant times τ ∈ (0, 1), i.e., in particular, up to the very instant of the EP degeneracy.

Example: N = 2
A deeper phenomenological meaning of requirement (32) may be most immediately illustrated via the N = 2 model.Recalling the set of all metrics (23) we notice that they are numbered by the single optional real variable α.It is easy to see that each deviation of this parameter from its unique, anisotropy-minimizing and extrapolation-friendly value as deduced from formula (24) would necessarily violate constraint (32).
The existence of the closed formula for the spectrum enables us to see that the influence of α changes from very weak (in the QC vicinity, i.e., at t ≪ 1) to very strong (near the Hermitian dynamical regime where λ(t) = 1).The extrapolation-friendly choice of α = π/4 of papers [44] or [49] appears exceptional.Solely for this choice of the free parameter the difference between the eigenvalues of the metric will strictly vanish at λ(t) = 1.We may just repeat that for the value of α = π/4 the τ = 0 instant really carries the meaning of a Hermitian onset of the fall into QC singularity.
Once we preserve the latter exceptional choice of the parameter α = π/4 at all times τ ∈ (0, 1), the difference between the two eigenvalues of the N = 2 metric (measuring a Hilbert-space anisotropy) remains always minimized (cf.Eq. ( 24)).In this formulation, the selection of a "minimal Hilbert-space anisotropy" principle of paper [49] is certainly confirmed as optimal.

Eigenvalues of the metrics
In the spirit of the methodical project as outlined in paper [49] the uniqueness of the choice of α = π/4 at N = 2 should be, mutatis mutandis, extended and amended to apply at any N. Preliminarily, such an idea has been tested and found feasible in Ref. [44] where we recalled formula (19) and where we managed to evaluate, up to N = 5, the ketket-eigenvectors |ψ (N ) n in closed form.Now we intend to amend the recipe and to find and formulate a more general result.Our task may be separated into two subtasks.In the first one (to be dealt with in this section) the problem will be reconsidered at a few smallest dimensions N. We will reveal that a new and promising guide to extrapolations in N can and should be sought in a certain very regular sparse-matrix pattern emerging in the formulae for the metrics Θ (N ) (τ ) [cf. also Eq.Nr. (10) in [44] in this respect].
Secondly, in a genuine climax of our present constructive efforts (and in a way described in subsequent section 7) we will find that the latter observation opens the way towards a remarkably efficient study and closed-form evaluation of the eigenvalues θ (N ) n (τ ) of the metrics.Very well reflecting both the anisotropy and asymptotic degeneracy of the physical Hilbert space and, hence, carrying a perceivably more useful information about dynamics than the matrix elements of the metric-operator N by N matrices Θ (N ) (τ ) themselves.

N = 2, revisited
The explicit construction of metric Θ (N ) (τ ) via auxiliary Schrödinger Eq. ( 18) is not too easy even at N = 2, i.e., for our first nontrivial QC-related Hamiltonian matrix The efficiency of this construction remains comparable with the brute-force solution of Eq. ( 31) (cf.subsection 3 above).Nevertheless, it still makes sense to rederive metric Θ (2) , by the amended method, for the pedagogical purposes.We may start from the real-matrix ansatz with the subscripted vector κ containing two arbitrary positive components.Next we fix an overall multiplication constant by setting the determinant equal to one.This enables us to put b = sinh ν and choose ε = ±1 in a = ε cosh ν exp ̺ and d = ε cosh ν exp(−̺).Both of the new parameters ν and ̺ are assumed real.The metric must be positive so that we may only use ε = 1.Finally we check that the matrix constraint (31) degenerates to the single, time-reparametrization item Our conclusion is that for any given τ ∈ (0, 1) we may choose any real ̺ ∈ (0, ̺ max ) (note that this is the parameter which makes the main diagonal of the metric asymmetric).This choice enables us to evaluate ν = ν(τ, ̺) from the latter equation (this implies that at a fixed time, the value of ̺ max must be such that cosh ̺ max = 1/τ ).Summarizing, we may set α 11 (1) = cosh ν exp ̺, α 12 (1) = cosh ν exp(−̺) and α 11 (2) = sinh ν in Eq. ( 26) at N = 2.The resulting eigenvalues of the metric are both, by construction, positive.At the very start of the fall of the system into the catastrophe, i.e., at τ = 0 one has ̺ max (0) = ∞ so that there is no upper bound imposed upon ̺(0).Still, as long as one might like to have the trivial, isotropic initial value of Θ (2) (0) ∼ I (implying the special choice of ν(0) = 0 and ̺(0) = 0) the resulting metric becomes, up to the above-mentioned irrelevant overall multiplication factor, unique at τ = 0.
During the subsequent growth of τ < 1 the requirement of the minimization of the anisotropy leads to the rule ̺(τ ) = 0 (cf.Eq. ( 36)) so that the remaining variable ν < 0 may now be interpreted as another version of the time of the QC degeneracy which is just rescaled and, incidentally, inverted (cf.Eq. ( 35)).
Once we return to the standard variables we get our unique and minimally anisotropic metric in the virtually trivial form From this formula we may deduce the special, minimally anisotropic version of eigenvalues in the form compatible with their more strongly anisotropic generalization (24).

N = 3
Whenever one tries to move to the higher matrix dimensions N one encounters the technical problem of an increase of multitude of parameters.In the first nontrivial case with N = 3 let us first follow the N = 2 guidance (cf. the ultimate choice of ̺ = 0 in the preceding paragraph 6.1) and let us omit the discussion of the metrics with an asymmetric form of their main diagonal.
Once we also keep ignoring the other, irrelevant though still existing overall factor, we are, after some straightforward manipulations using Eq. ( 31), left with the last free parameter g in the metric Among its three readily obtainable eigenvalues the middle one (with an inverted-parabola dependence on τ ) remains positive for the parameters g < 1/τ 2 .
The change of sign of the remaining two eigenvalues takes place at the curves g = 1/τ 2 and g = 1/(2 − τ 2 ) in the g − τ plane.As a consequence, the correct and unique choice of the parameter is g = 1, yielding again the unique metric with the expected τ −dependence of the eigenvalues as given by Eq. ( 38).

N = 4
In the next step of our constructive considerations we are getting beyond the formulae derived in older papers.We succeeded because already the N = 3 formula (39) offered the hint.Thus, making use of the analogy and performing an extrapolation, it just proved sufficient to verify that the following tentative candidate for the metric obeys all the necessary and sufficient requirements.They include the validity of the Dieudonné's Eq. ( 31) as well as the feasibility of evaluation of the τ −dependent eigenvalues of the candidate for the metric.We immediately see that they behave as they should, We may now combine the results of preceding section with the contents of Theorem 1.Using an elementary insertion in Eq. ( 31) we may easily prove that in the expansions (26) of the metrics with minimal anisotropy the diagonal-matrix coefficients (27) may be defined, at all N, by the elementary formula Similarly, the closed formula is also available for the antidiagonal coefficients in M (N ) (N), Next, the bidiagonal matrix coefficients (28) may be defined, at all N, by the slightly less elementary general formula Due to the easily verified symmetry, the analogous formula exists for the coefficients in Up to now, unfortunately, we did not succeed in an extension of these observations to the tridiagonal sparse matrix coefficients (29), etc.Nevertheless, we believe that the task is not impossible.This belief seems supported by Theorem 1, i.e., by the reducibility of the N by N matrices M (N ) (k) with k = 3, 4, . . . to the respective auxiliary k by N − k + 1 arrays containing the non-vanishing matrix elements α jm (k) of M (N ) (k).
The first missing set of coefficients occurs at N = 5.Its values should be better rewritten in the compact form of array which we derived using the brute force construction based on Eq. (31).It again deserves the compact presentation as array The closed form of the latter result indicates that there might exist a not too complicated extrapolation recipe, with the help of which we would be able to determine the unique, minimally anisotropic metric at any dimension N.This belief seems further supported by the regularity and apparent extrapolation-friendliness of the next two sparse-matrix "missing" coefficients and .
They were obtained, with the assistance of the computerized symbolic manipulations, by the brute-force solution of the set of 49 linear algebraic Eqs.(31).In the above-described constructions of the N by N matrices of metric Θ (N ) we did not manage to find, unfortunately, any obvious general extrapolation tendency or pattern.

Eigenvalues at arbitrary N
For this reason, we turned attention from matrices to the perceivably simpler-to-display N−plets of their eigenvalues θ (N ) n (τ ).At the first few values of N we performed the bruteforce calculations.We met the ultimate success which may be given the following form of proposition.

Proposition 2
The time-dependent eigenvalues of Θ (N ) (τ ) may be written in the form of polynomials where the numerically evaluated values of the coefficients C (N ) nk may be found listed, up to N = 8, in the Pascal-like schemes of Tables 1 -4.3n in Eq. ( 44) After a cursory inspection of the latter four Tables 1 -4 one immediately finds, in all of them, a regularities resembling the well known Pascal triangle.The analogy is almost perfect.It enabled us to reveal the extrapolation pattern and to specify the recurrences for coefficients which appeared all solvable in terms of binomial coefficients.In other words, the closed-form eigenvalues (44) of the metrics have been obtained for any time τ and for any matrix dimension N.This is our present main result.
Theorem 3 The time-dependent eigenvalues of metrics Θ (N ) (τ ) of Eq. ( 26) are given by formula 4n in Eq. ( 44) where Proof is straightforward and proceeds by mathematical induction.

Unbounded differential-operator Hamiltonians
In the current literature the merits and applicability of the THS formalism are most often illustrated by the replacement of the most common harmonic or anharmonic oscillator by the Bender's and Boettcher's [52] family of non-Hermitian power-law-interaction models H (BB) (δ) = p 2 + i δ x 2+δ which are characterized just by the single real exponent δ ≥ 0.
In the historical perspective such a choice was surprising but well motivated by the needs of development of quantum field theory (cf., e.g., [53]) and/or of perturbation expansion methods [54] - [57].The use of nontrivial and, in general, manifestly Hamiltonian-dependent Hilbert-spacemetric operators may be perceived as an important innovation of the model-building in quantum theory.Nevertheless, whenever accepted as a sound theoretical tool in physics, its mathematical consistency must always be carefully reexamined.In this sense the proof of nonexistence of the metric operator for the most popular and phenomenologically highly relevant H (BB) (δ) [58] makes the study of this particular Hamiltonian far from being completed.Hence, the illustration purposes are, in the eyes of mathematicians, much better served by the bounded-operator Hamiltonians [11] as sampled in our paper.
This being said, it is still possible to conclude that in comparison with the conventional quantum theory using selfadjoint operators the question of interpretation of the Bender's and Boettcher's unbounded-Hamiltonian models, albeit still open [59], remains inspiring.Especially in the light of the abstract mathematical comments by Dieudonné [8] who pointed out that not only the necessary ad hoc specification but even the very proof of existence of the correct physical Hilbert space of states may be a highly nontrivial question.

Quantum catastrophe as an EP-related concept
The classical Thom's concept of a "catastrophe" [42] is based on the idea that in a certain dynamical regime an infinitesimal change of one or more relevant parameters leads to an abrupt change of the behavior of the system.In Ref. [44] we conjectured that a quantum analogue of such a singularity can be represented by the Kato's [45] EP singularity at which a multiplet of quantum bound-state energies merges and ceases to be observable.
In spite of the fact that one of the older reviews of the related theory (viz., Ref. [11]) already appeared as early as in 1992, an extension of its scope beyond the domain of nuclear physics was not too quick.Fortunately, multiple extensions already do exist at present, based on the innovations of the theory called PCT −symmetric [5] alias pseudo-or quasi-Hermitian [6] or Krein-space self-adjoint [60].In this context the construction of illustrative quantum catastrophes in [44] has been facilitated by the availability of Hamiltonians for which at least some of the Kato's EP singularities (which occurred, traditionally, just at the complex, "unphysical" values of parameters) became experimentally accessible, in principle at least [14,15].
In our present paper we decided to reanalyze, therefore, some of the latter models, with the emphasis put upon the open questions concerning a smooth transition from the classical to quantum dynamics.Our approach has been based on the use of the THS formalism, with the aim of a further amendment of our understanding of the possible definition and systematic description of a catastrophic evolution scenario in the language of quantum theory.

Closed formulae
Using a specific model we managed to cover several methodical topics including not only the ambiguity of the metric (i.e., of the necessary specification of one of the eligible physical Hilbert spaces) but also its descriptive aspects and appeal.Thus, we paid attention to the anisotropy of the alternative Hilbert-space geometries as well as to the constraints imposed by the necessary positive-definiteness of the mathematically correct inner products in these spaces.
We formulated several arguments in favor of our choice of the toy model.Firstly, it enabled us to employ the metodical idea of paper [49].Thus, we played with the free parameters in order to minimize the anisotropy of the physical Hilbert-space metric Θ.Secondly, in the context of study of the EP degeneracy as initiated in [44] we found it useful to invert the arrow of time.Thus, in place of the time t which starts at the EP singularity we used another time variable τ which runs in opposite direction.This enabled us to set the initial zero long before the fall of the system into its physical QC singularity.Thirdly, we found it productive to start our analysis from the systems with the smallest level-multiplicities N ≤ 4. Using the brute-force linear algebra methods we managed to construct the fully explicit matrices of the metrics which appeared (and were declared) optimal and unique.
Ultimately, the transparency and compact form of the results of the brute-force linearalgebraic calculations opened the way towards extrapolations.Their use (followed by the decisively facilitated formal proofs) finally enabled us to extend the validity of some of our previous empirical small−N observations to arbitrary Hilbert-space dimensions N.
On this background, the qualitative features of the QC process have been shown related to the explicit qunatitative properties of our manifestly time-dependent physical Hilbertspace metrics.Thus, in a climax of the story the picture of the N−tuple QC leveldegeneracy scenario was given the form supported by the closed formulae reflecting, via the eigenvalues of the metric, both the time-dependent anisotropy and asymptotic degeneracy of the system in question.

Summary
In our present toy model the evolution in time was explained as proceeding from an innocent-looking and safely Hermitian equidistant-energy-level onset prepared at an initial time τ = 0 up to an ultimate collapse realized via a complete, N−tuple EP degeneracy of the energy spectrum at the final QC time τ = 1.
Our paper offered a compact and consistent picture of the process in which the not quite expected exact solvability of our toy model enabled us to cover all times τ ∈ (0, 1).We were able to describe the quantum-evolution fall of the system in the level-degeneracy quantum catastrophe, and we were able to explain such a collapse as a consequence of an unlimited growth of the anisotropy of the underlying time-dependent Hilbert space H (A) .
We found it natural to characterize an optimal version of the latter process by a minimal spread of the set of eigenvalues θ (N ) n of the related physical inner-product metric Θ (N ) .We decided to make such a metric unique via a minimization of the latter anisotropyrepresenting spread, with an emphasis put upon the zero limit of the special measure ρ = max(θ m ) of the spread at the onset τ = 0 of the process.We managed to match our metric smoothly to both of its extremes, i.e., not only to the most common isotropic metric at τ = 0 but also to the asymptotically degenerate metric at τ = 1.In between these two extremes the operator (i.e., matrix) has been kept smooth, nontrivial and optimal during all τ ∈ (0, 1).From the point of view of phenomenology, we arrived at a benchmark quantum representation of the EP-related catastrophe in which the fall into the degeneracy appeared realized in finite time.
Our toy-model simulation of the catastrophe can be perceived as initiated by an arbitrary conventional unitary-evolution prehistory at τ < 0. According to the general principles of quantum theory the states of the system during its EP-related degeneracy at τ ∈ (0, 1) were assumed described differently, by a THS wave function ψ which evolves in time in a way which has thoroughly been described in Ref. [18].In our present paper we skipped most of the related technical details and we restricted our attention just to the description of the interplay of time-dependence between a pre-selected "non-Hermitian" benchmark Hamiltonians H (N ) (τ ) and one of the eligible "Hermitizing" metric operators Θ (N ) (τ ).
Our choice of the latter operator can be characterized as truly exceptional: In the context of mathematics its form has been shown to lead to a unique, minimally anisotropic geometry of the physical Hilbert space of states.In the context of physics we emphasized that a decisive merit of the time-dependence of our inner-product metric should be seen in the guarantee of existence of this metric up to an arbitrarily small vicinity of the ultimate catastrophic quantum collapse of the system.

Table 1 :
Pascal triangle for coefficients C

Table 2 :
Pascal-like triangle for coefficients C

Table 3 :
Pascal-like triangle for coefficients C

Table 4 :
Pascal-like triangle for coefficients C