“…In spite of the manifest non-Hermiticity of the P T −symmetric candidates H for Hamiltonians, these operators were shown eligible as generators of unitary evolution [12,14]. In this context, one of the basic methodical assumptions accepted in the current literature on BH models [19,20,[34][35][36][37] was that the infinite-dimensional matrix (5) as well as all of its separate submatrices (6) had to be complex symmetric, tridiagonal and P T −symmetric, with P equal to an antidiagonal unit matrix, and with symbol T representing an antilinear operation of Hermitian conjugation (i.e., transposition plus complex conjugation). This led to the conclusion (or rather conjecture) that in the EPN limit (i.e., at the instant of the loss of diagonalizability), the canonical representation of every N by N submatrix H (N) (γ (EP) ) can be given the form of the N by N Jordan matrix (16) with, due to P T −symmetry, η = 0.…”