Elsevier

Annals of Physics

Volume 327, Issue 3, March 2012, Pages 893-913
Annals of Physics

N-site-lattice analogues of V(x)=ix3

https://doi.org/10.1016/j.aop.2011.12.009Get rights and content

Abstract

A discrete N-level alternative to the popular imaginary cubic oscillator is proposed and studied. As usual, the unitarity of evolution is guaranteed by the introduction of an ad hoc, Hamiltonian-dependent inner-product metric, the use of which defines the physical Hilbert space and renders the Hamiltonian (with real spectrum) observable. Due to the simplicity of our model of dynamics the construction of the set of eligible metrics is shown tractable by non-numerical means which combine the computer-assisted algebra with the extrapolation and/or perturbation techniques.

Highlights

► Elementary non-Hermitian quantum-lattice Hamiltonians H(a,z) with real spectra proposed. ► Exceptional points a=α calculated. ► All topologically nonequivalent stability-loss patterns shown reachable via “exponent” z. ► “Critical exponents” of PT-symmetry breakdowns assigned Fibonacci-related classification and calculated. ► Hermitization matrices Θ shown obtainable (and samples constructed) via symbolic manipulations and extrapolations in N.

Introduction

Many measured spectra of energies may be interpreted as excitations of a quasiparticle. The simplest fits of such a type employ the elementary one-dimensional differential Schrödinger equation d2dx2ψn(x)+V(x)ψ(x)=Enψn(x)ψ(±Λ)=0,Λ containing a real potential V(x). In more sophisticated models, potential V(x) may even be allowed complex, provided only that the spectrum itself remains real (cf. reviews [1], [2], [3], [4] for details).

Whenever V(x)V(x), the reality of the spectrum may be fragile and sensitive to perturbations [5]. A remarkable exception emerges with the robustly real spectra generated by many potentials with the property V(x)=V(x) called, in the literature, PT-symmetry [2] alias parity–pseudo-Hermiticity [3] alias Krein-space-Hermiticity [6].

In the early studies of this remarkable mathematical phenomenon the attention of the authors has mainly been restricted to the imaginary cubic oscillator example V(x)=V(IC)(x)=ix3 modified, possibly, by some other, asymptotically subdominant terms (cf., e.g., papers by Caliceti et al. [7], Alvarez [8] or Bender et al. [9], [10] and several further authors [11]). In what follows, we shall mainly feel inspired by this choice of V(x) as well.

The essence of our present message will lie in the recommendation of a drastic simplification of the necessary mathematics. In a way explained in Section 2 this will be achieved by means of the replacement of the differential Schrödinger equation (1) by its discrete, difference-equation analogue. A few simplest illustrations will be then added in Section 3 where the discretization of the coordinate will be shown to facilitate the study of the reality (i.e., in principle, observability) of the spectrum.

The desirable flexibility of our present discrete simulation of the dynamical energy spectra will be shown achieved by a supplementary one-parametric deformation of the potential resembling slightly the influential proposal of deformation ix3(ix)3+δ by Bender and Boettcher [10]. The details will be described in Section 4. The subsequent Section 5 will then recall the known theory and explain some of its details via the simplest possible example with N=2. In Sections 6 The first nontrivial, 7 The, 8 Extrapolations and models with we shall finally present some applications of this theory to the models with N=4,N=6 and general N=2K8, respectively. Section 9 is summary.

Section snippets

Discrete Schrödinger equations

The spectrum of many Krein-space-Hermitian Hamiltonians HH has been found robustly real and bounded below [10], [12], [13]. In the spirit of the general theory as outlined, first, by Scholtz et al. [14], one can conclude that the apparent non-Hermiticity is “false” and that it may be reinterpreted as a mere consequence of an inappropriate choice of the Hilbert-space representation L2(R)H(F). Hence, the abstract remedy is straightforward and lies in an interpretation-mediating transition to a

The series of one-parametric Hamiltonians

Let us pick up the potential V(xj)=ixj3 and insert it in Eq. (6). Once we abbreviate a=h5/8 we obtain the following one-parametric sequence of Hamiltonians H(2)(a)=(2ia112+ia),H(3)(a)=(28ia10121012+8ia),H(4)(a)=(227ia10012ia10012+ia10012+27ia),. The first element H(2)(a) of this series is particularly elementary. For this reason it has also been chosen and studied in the methodical part of Ref. [15] (cf. its section II. B).

For the sake of brevity we shall restrict our attention

The series of two-parametric Hamiltonians

In the examples of the preceding section we encounter another suspicious feature even at any fixed and finite dimension N=2K. Indeed, the critical left and right points a=±α(2K) of the loss of the reality of the spectrum seem to be exclusively connected with the confluence of the first and second lowest excited states ε2 and ε3 or, symmetrically, of their equally “privileged” mirror partners εN2 and εN1.

In what follows we shall explain this apparent privilege as an artifact caused by the too

Interpretation

Usually [25] people decide to work in a fixed, specific Hilbert space H(S) and treat a given quantum system as physical if and only if its time evolution is unitary. Naturally, this is the strategy which does not change if the space H(S) proves endowed with a general, nontrivial metric Θ=Θ>0 (i.e., with its inner product defined in terms of this metric, cf. Eq. (2) above). One still has to guarantee that every candidate for an operator of observable proves self-adjoint with respect to this

The complete set of pseudometrics

The general four-parametric Hermitian N=4 candidates for the metric (= “pseudometrics”) may be obtained, most easily, from Eq. (10) again. These solutions appear symmetric with respect to their second diagonal. They may be written in the closed four-parametric form of matrix Θ[H(4)(a,z)](k,m,r,h)=[kmikwWZm+ikwrhi(kw+ra)WWh+i(kw+ra)rmikwZWm+ikwk] where we introduced function w=w(z,a)=3za (this quantity must be positive and, for z>0, larger than a) and where we abbreviated W=W(k,m,r)=w2k+r

The N=6 model

For the not too large matrix dimensions N=2K the technique of the construction of the general pseudometrics Θ via the solution of the linear algebraic equation Eq. (10) remains feasible and straightforward. With the growth of N the only difficulty emerges in connection with the printed presentation of the resulting multiparametric set of pseudometrics. For this reason it makes sense to find a sufficiently representative set of a few key parameters. For each such choice, moreover, it becomes

Extrapolations and models with N8

At any dimension N=2K the dynamical contents of our toy models is controlled by the Hamiltonian (which varies with the “kinematical” parameter a(α(N),α(N)), “dynamical” parameter z(,) and real spectral locus {εj(a,z)}) and by the Hermitian, positive-definite metric Θ(N) specified by an N-plet of parameters {k,m,}D.

Naturally, practical implementations of such a recipe will require also a determination of the whole positivity domain D, or of its suitable non-empty subdomain at least. This

Summary

Our first tests of the idea of discretization proved disappointing. We observed that the growing-dimension series (7) of the simple-minded one-parametric descendants of the popular differential imaginary cubic oscillator do not offer a sufficiently rich variability of the coupling-dependence of the energy spectra. Fortunately, the merely slightly more sophisticated and re-scaled two-parametric choice (8) of the model has been shown to offer a flexibility of the spectral loci which covers a

Acknowledgments

This work is supported by the GAČR grant No. P203/11/1433, the MŠMT “Doppler Institute” project No. LC06002 and the Institutional Research Plan AV0Z10480505.

References (32)

  • H. Langer et al.

    Czech. J. Phys.

    (2004)
  • M. Znojil

    SIGMA

    (2008)
  • M. Znojil

    J. Phys. A

    (2010)
  • M. Znojil

    Phys. Rev. D

    (2009)

    Phys. Rev. D

    (2009)
  • P. Dorey et al.

    J. Phys. A

    (2007)
  • C.M. Bender

    Rep. Progr. Phys.

    (2007)
  • A. Mostafazadeh

    Int. J. Geom. Methods Mod. Phys.

    (2010)
  • M. Znojil

    SIGMA

    (2009)
  • M. Znojil

    J. Math. Phys.

    (2004)
    M. Znojil

    Phys. Lett. B

    (2007)
    Y.N. Joglekar et al.

    Phys. Rev. A

    (2010)
    A. Cavaglia et al.

    J. Phys. A

    (2011)
    H. Schomerus

    Phys. Rev. A

    (2011)
    Y.N. Joglekar et al.
  • E. Caliceti et al.

    Comm. Math. Phys.

    (1980)
  • G. Alvarez

    J. Phys. A: Math. Gen.

    (1995)
  • C.M. Bender et al.

    Phys. Rev. D

    (1997)
  • C.M. Bender et al.

    Phys. Rev. Lett.

    (1998)
  • E. Delabaere et al.

    J. Phys. A: Math. Gen.

    (2000)
    M. Znojil

    Phys. Lett. A

    (2010)
    J. Zinn-Justin et al.

    J. Phys. A

    (2010)
  • P. Dorey et al.

    J. Phys. A: Math. Gen.

    (2001)
  • G. Lévai et al.

    J. Phys. A: Math. Gen.

    (2000)
    G. Lévai et al.

    J. Phys. A: Math. Gen.

    (2003)
  • Cited by (11)

    View all citing articles on Scopus
    View full text