Nonclassical States for Non-Hermitian Hamiltonians with the Oscillator Spectrum

We show that the standard techniques that are utilized to study the classical like properties of the pure states for Hermitian systems can be adjusted to investigate the classicality of pure states for non-Hermitian systems. The method is applied to the states of complex-valued potentials that are generated by Darboux transformations and can model both non-PT-symmetric and PT-symmetric oscillators exhibiting real spectra.


Introduction
In quantum mechanics one finds two important classes of states of the radiation field: Fock and coherent states. The former, introduced by Fock at the dawn of quantum theory, contain a precise number of photons and produce average fields equal to zero. With exception of the vacuum, the Wigner distribution [1] exhibits negative values in some regions of the phase-space when it is evaluated with the Fock states. Thus, the properties of the quantum states of radiation fields that are occupied by a finite number of photons are far from the Maxwell theory. On the other hand, the coherent states were formally introduced in quantum optics by Glauber [2], although the first antecedents can be traced back to the Schrödinger papers on quantization, see e.g. [3]. They produce average fields different from zero as well as nonnegative symmetrical Wigner distributions. Although the coherent states are constructed as superpositions of Fock states, it is remarkable that their properties are very close to the Maxwell theory. This is because such states satisfy the notion of full coherence introduced by Glauber, while the Fock states (for n = 0) lack second (and higher) order coherence, and thus they are nonclassical. The set of nonclassical states includes also squeezed states [4], even and odd coherent states (also called Schrödinger cats) [5], binomial states [6,7], photon-added coherent states [8], etc. The difference between the classical and nonclassical properties of a given quantum state is strongly linked to the notion of entanglement [9], which is a fundamental concept required in the development of quantum computation and quantum information [10]. Noticeably, while investigating the quantum properties of a beam splitter, it has been found that the "entanglement of the output states is strongly related to the nonclassicality of the input fields" [11]. That is, the nonclassicality of an arbitrary (input) state can be tested by detecting entanglement in the output of a quantum beam splitter [11,12].
In the present work we study the nonclassical properties of the states of non-Hermitian Hamiltonians with real spectrum that are generated by Darboux transformations [13,14]. In general, the systems associated with non-Hermitian Hamiltonians are subject of investigation in many branches of contemporary physics [15]. The applications include the study of unstable (decaying) systems [16], light propagating in materials with complex refractive index [17], P T -symmetric [18,19] and P T C-symmetric [20] interactions, multiphoton transition processes [21], diverse measurement techniques [22], and coherent states [3,23,24], among others. In the case of the Darboux deformed non-Hermitian systems, the eigenfunctions obey a series of oscillation theorems [25] that permit their study as if they were associated to Hermitian systems. Indeed, the representation space of such non-Hermitian systems can be equipped with a bi-orthogonal structure that provides complete sets of orthonormal states [13]. We shall focus on non-Hermitian oscillators since their bases of eigenfunctions can be used to construct optimized binomial states [26] and generalized coherent states [24], which in turn give rise to a wide diversity of pure states as particular cases. One of the main results reported in this paper is to show that the techniques used to study classical like properties for Hermitian systems can be adjusted to investigate the classicality of the states of Darboux deformed non-Hermitian systems.
The outline of the paper is as follows: In Section 2 we revisit the main properties of the non-Hermitian oscillators we are interested in. The main differences between the conventional orthogonality and the bi-orthogonality are discussed in detail. In Section 3 we introduce the different bi-orthogonal superpositions of states that are to be analyzed, these include the optimized binomial states and diverse forms of generalized coherent states. Section 4 deals with the nonclassicality of the states defined in Section 3. Final comments are given in Section 5. We have added three appendices with relevant information that is used throughout the manuscript concerning operator algebras (Appendix A), and criteria of nonclassicality (Appendix B).

Non-Hermitian oscillators
The solution to the eigenvalue problem of the dimensionless (mathematical) oscillator H osc = − d 2 dx 2 + x 2 is given by the discrete eigenvalues E osc n = 2n + 1 and normalized eigenfunctions ϕ n (x) = e −x 2 /2 2 n n! √ π H n (x), n = 0, 1, . . . , where H n (x) stands for the Hermite polynomials [27]. Using the linearly independent solutions of the eigenvalue equation H osc u = −u and a set of real numbers {a, b, c, λ} fulfilling 4ac − b 2 = 4λ 2 , it can be shown that is a real-valued function which is free of zeros in R, and defines the dimensionless potential as a Darboux transformation of V osc (x) = x 2 [13]. Here Erf(x) defines the error function.
is a complexvalued function. Nevertheless, it may be shown that the imaginary part of such potential is continuous in R and satisfies the condition of zero total area [25]: The condition (4) implies a balanced interplay between gain and loss of probability that does not depend on any other symmetry of either ImV λ (x) or ReV λ (x). For instance, besides fulfilling (4), the profile of the real and imaginary parts of the oscillator depicted in Fig. 1(a) is not particularly special (this has been obtained from (3) with a = 5, b = 8, c = 5 and λ = 3). Looking for concrete symmetries we can take b = 0 to get even and odd forms of ImV λ (x) and ReV λ (x) respectively. The result is a non-Hermitian oscillator which is invariant under P T -transformations [18]. This case is illustrated in Fig. 1(b) for a = c = λ = 1, and b = 0. On the other hand, for b = ±2 √ ac, with a and c positive numbers, the parameter λ is equal to zero and the expression (3) decouples into the pair of real-valued functions where the condition γ > √ π/2 must be satisfied to avoid singularities. Fig. 1(c) shows V − λ=0 (x; γ) with a = 1 and c = 0.88. The concomitant V + λ=0 (x; γ) corresponds to the specular reflection of V − λ=0 (x; γ). The Hermitian oscillators (5) were first found by Abraham (a) non P T -symmetric (b) P T -symmetric and (c) Hermitian. In (a) and (b) the imaginary part (dashedred) satisfies the condition of zero total area (4). The zones of gain and loss of probability are shadowed in orange and red respectively. In all the cases the curve in blue corresponds to the real part of V λ (x) and the dotted-black curve represents the shifted oscillator V osc (x) − 2, the latter is obtained from (3) as a particular case.
and Moses through the systematic use of the Gelfand-Levitan equation [28], and then recovered by Mielnik as an application of his generalized factorization method [29].
Additionally, the expressions (5) represent a family of systems that converge to the conventional oscillator in the limit |γ| → ∞. Indeed, it may be verified that The oscillator V osc (x) − 2 is thus a member of the family of complex-valued potentials V λ (x) introduced in (3). This is depicted in the plots of Fig. 1 as a reference.
• Fundamental solutions. Another remarkable profile of the complex-valued oscillators (3) is that the functions are normalizable solutions of the related eigenvalue equation H λ ψ n+1 (x) = E n+1 ψ n+1 (x) for E n+1 = 2n + 1. The additional normalizable function also solves the eigenvalue equation defined by H λ , but it is not derivable from the set ϕ n (x) and belongs to the energy E 0 = −1. The spectrum of the potential V λ (x) is therefore composed of the equidistant energies E n = 2n − 1, n ≥ 0.
• Bi-orthogonality. Although the functions ψ n (x) are normalizable, they form a peculiar set since ψ 0 (x) is orthogonal to all the ψ n+1 (x) while the latter are not mutually orthogonal. As a result, in contrast with the Hermitian case, the norm of any superposition of states ψ n+1 (x) depends not only on the modulus of the related coefficients but also on the phase shift between them. For instance, the norm of satisfies ||Ψ|| 2 = |θ| 2 + |µ| 2 + 2|θµ| ξ n+1,m+1 cos(δ + γ n+1,m+1 ), δ = arg(θ * µ), where z * is the complex-conjugate of z ∈ C. The numbers ξ n+1,m+1 and γ n+1,m+1 are respectively the modulus and argument of the product between ψ n+1 (x) and ψ m+1 (x). Thus, depending on the phase-shift δ, the non-orthogonality of the set ψ n+1 (x) produces the oscillations of ||Ψ||. The complexity of such a dependence increases with the number of elements in the superposition.
One can face the above difficulties by introducing a bi-orthogonal system (see relevant information in Ref. [24]), formed by the eigenfunctions ψ n (x) of H λ and those of its Hermitian-conjugate H † λ ≡ H λ , written as ψ m (x). The main point is that the bi-product ψ m , ψ n = R ψ * m (x)ψ n (x)dx is equal to zero if n = m, and serves to define the bi-norm Therefore, besides the conventional normalization ψ n (x)/||ψ n ||, we have at hand the binormalization ψ n (x)/||ψ n || B , with ||ψ n+1 || B = 2(n + 1) for n ≥ 0. The bi-norm ||ψ 0 || B of the ground state depends on the set {a, b, c, λ} [24]. The real and imaginary parts of ψ n (x), as well as the probability density |ψ n (x)| 2 , behave qualitatively equal in both normalizations but their bi-normalized values are usually larger than those obtained with the conventional normalization. Such a difference is reduced as n increases. This property is illustrated in Fig. 2 for the first three eigenfunctions ψ n (x), and the corresponding probability densities |ψ n (x)| 2 , of the complex-valued oscillator depicted in Fig. 1(a).
The bi-orthogonal approach avoids the interference produced by the non-orthogonality. For instance, if the sates in (9) are substituted by their bi-normalized versions, we obtain the bi-orthogonal superposition The bi-norm of the latter state does not depend on the phase-shift so that Ψ B (x) is uniquely bi-normalized to 1 by the constant (|θ| 2 +|µ| 2 ) −1 . An additional property of the bi-orthogonal superpositions is that the values of their real and imaginary parts, as well as the values of the corresponding probability densities, are shorter than those obtained from the conventional approach. We show this property in Fig. 3 for the superpositions (9) and (11) with n = 0, m = 1, |θ| = |µ| = 1/ √ 2, and three different values of the phase-shift δ. Formally, one may say that Ψ(x) and Ψ B (x) represent two different superpositions of the states ψ n+1 (x) and ψ m+1 (x). In the sequel we shall take full advantage of the mathematical simplifications offered by the bi-orthogonal approach.
• Operator algebras. It can be shown that there exist at least two different algebras of operators associated with the eigenstates of the complex-valued oscillator V λ (x) [30].  Fig. 1(a). The eigenfunctions are depicted in the Argand-Wessel representation with Re ψ n (x) and Im ψ n (x) in the horizontal and vertical axes, respectively. They have been constructed from (7) and (8)  They are generated by two different pairs of ladder operators (see Appendix A for details). The first pair, A and A + , together with the Hamiltonian H λ , satisfy the quadratic polynomial (Heisenberg) algebra The second pair of ladder operators, denoted by C w and C + w , together with the Hamiltonian H λ , and an additional operator I w , satisfy the distorted (Heisenberg) algebra where w is a non-negative parameter that defines the 'distortion' suffered by the oscillator algebra when one substitutes the operators C w , C + w , and H λ for the conventional boson operators,â,â † , andn, respectively. The algebras (13) and (14) will serve to the analysis of classicality of the states associated with the complex-valued oscillators (3).

Bi-orthogonal superpositions
We shall study the properties of diverse superpositions of states ψ n (x) by using (i) Bi-normalization to obtain regular probability densities.  (ii) Bi-orthogonality to avoid the interference associated with non-orthogonality.
The space of states of H λ , denoted H λ , consists of all the bi-orthogonal superpositions such that In turn, H λ denotes the space spanned by the states ψ n (x). The series decomposition of φ(x), which is concomitant to φ(x) in (16), is obtained from (15) by changing ψ k (x) → ψ k (x). To simplify the notation, hereafter we use no label to denote the bi-orthogonal properties of the system {ψ n (x), ψ m (x)}.
The approach we are going to use identifies a threefold partnership between H λ , H λ , and the conventional oscillator H osc , where the eigenvalues play a main role. That is, the triad {H λ , H osc , H λ } represents a system for which the energies E n+1 = E osc n = 2n + 1 are threefold degenerate while the energy E 0 = −1 is only twofold degenerate. In this sense, our model is a generalization of the supersymmetric approach [31] since H λ and H λ are treated as two different faces of the same system, which can be studied in much the same way as in the Hermitian approaches [24].
We are interested in two different classes of superpositions. The optimized superpositions are such that their nonclassicality can be manipulated by tuning the number of elements as well as the state of lowest energy in the sum. The generalized coherent states are constructed in terms of the algebras underlying the complex-valued oscillators (3) and form an over-complete set in H λ . Next, we analyze these classes in detail.

Optimized binomial states
Let us consider a superposition of N + 1 adjacent states where the non-negative integer r ≥ 0 determines the state ψ r (x) of the lowest energy that is included in the packet, and the coefficients c b k give rise to the binomial distribution That is, the sum (17) is bi-normalized and includes only the eigenstates belonging to the energies E r , E r+1 , . . . , E r+N . The distribution (18) means that the probability of finding the system in the state ψ k+r (x) is weighted by the probability p of having success k times in n trials. The superpositions (17) are called optimized binomial states [26] since (for r = 0) they converge to the conventional binomial superposition at the oscillator limit.
To investigate the properties of φ b (x) let us take p ≈ 0. The modulus |c b 0 | of the first coefficient in the expansion is close to 1 and the other coefficients are almost zero. Then, the state ψ r (r) is very influential in the behavior of the entire packet since it is weighted by c b 0 in the superposition. The situation changes if p ≈ 1 since in this case the modulus |c b N | of the last coefficient in the sum is close to 1, while the moduli of the other coefficients are almost zero. In such case, ψ N +r (x) is the state of major influence in the packet.
Another important feature of the optimized binomial states is that their mean energy E b depends on the number r which labels the lowest of the involved energies That is, E b is bounded from below by r, and from above by N and r. Fixing E b and r in (19) we can write Then, from (18) we arrive at the Poisson distribution: In this form, as a limit version of the optimized binomial state φ b (x), we introduce the optimized Poisson state In the quantum oscillator limit (6), the superposition (23) acquires the form which coincides with the conventional coherent state for r = 0, see Eq. (B-25) of Appendix B.

Generalized coherent states
The ground state ψ 0 (x) is annihilated by both pairs of ladder operators A, A + , and C w , C + w , see Appendix A. In any case we can construct coherent states as eigenstates of the respective annihilator operator. We write The super-index γ stands for natural (N ), distorted (w), natural displaced (N d ), and distorted displaced (w d ). Such nomenclature refers to the form in which the coefficients c (γ) k have been selected, as it is explained below. For simplicity, we use Dirac notation. If the superposition (25) is an eigenvector of the natural annihilation operator A with complex eigenvalue z, then it is a natural coherent state [24]. We write 1 Notice that also the ground state |ψ 0 is eigenvector of A, but its complex eigenvalue is equal to zero. Moreover, as |z = 0 = |ψ 1 , we see that the eigenvalue z = 0 is twice degenerate. The straightforward calculation shows that the vectors |ψ 0 and |z minimize the uncertainty relation (A-9). In this sense |ψ 0 and |z represent two different types of minimal uncertainty states.
Similarly, the distorted coherent states [24] are defined as eigenvectors of the distorted annihilation operator C w with complex eigenvalue z. In this case we write with (w) k = w(w − 1) · · · (w − k + 1) being the Pochhammer symbol [27]. The eigenvalue z = 0 is also twice degenerate because |ψ 0 and |z = 0, w = |ψ 1 are annihilated by C w . As |ψ 0 and |z, w minimize the uncertainty relation (A-16), they represent two different types of minimal uncertainty states.
Another interesting superposition (25) is obtained by demanding that the vector |φ (γ) be the displaced version of a fiducial state. The states |ψ 0 and |ψ 1 are invariant under the action of the operator exponentiations e −z * A and e −z * Cw , because they are annihilated by both operators A and C w . The latter means that |ψ 0 and |ψ 1 can be used as fiducial states. Considering the operators one has the 'displaced' states |φ (N d ) = D(z)|ψ 1 and |φ (w d ) = D w (z)|ψ 1 . Remark that the ground state |ψ 0 is invariant under the action of D(z) and D w (z), so it is also a displaced state.
In the rest of this work we shall omit the description of the states |φ (N d ) since their properties are qualitatively similar to those of the (distorted) displaced coherent states: In contrast to the previous superpositions, the states |z, w d do not minimize the uncertainty relation (A- 16) for arbitrary values of z, but only for z = 0. The latter means that |ψ 0 and |ψ 1 are the only displaced coherent states of minimal uncertainty.

Nonclassical states for non-Hermitian oscillators
The P -representation introduced by Glauber [32] and Sudarshan [33] defines a limit for the classical description of radiation fields. If the P -function is not accurately interpretable as a probability distribution then the field "will have no classical analog" [2]. This notion of nonclassicality applies immediately to the Fock states |n for n = 0 because their P -functions are as singular as the derivatives of the delta function δ (2) (z) = δ(Re(z))δ(Im(z)), and are negative in some regions of the complex z-plane. The latter means that the light beams represented by any of the Fock states |n = 0 cannot be described in terms of the electrodynamics introduced by Maxwell. In turn, the Glauber states [2], usually denoted |α , are classical because their P -function is precisely δ (2) (z−α). The vacuum |0 is also classical as it is a coherent state with complex eigenvalue equal to zero. In general, the states for which the quadrature variances obey the inequalities (∆x) 2 ≥ 1 state using the Wigner distribution [1], which is regular and always exists. Nonclassical states have Wigner functions which either are negative in some regions of the phase-space or are squeezed in one of the variables of the phase-space.
In modern days, the nonclassicality of a state can be tested by another excellent method, which follows from the use of the quantum beam splitter. It creates entangled states at its output ports while at least one of its inputs are fed with a nonclassical state. The device is used frequently both in theory and experiment, since it not only tests the nonclassicality of the state, but also quantifies the degree of nonclassicality in an efficient way [11,12]. There are some well-known nonclassical states in the literature, which include squeezed states [4], even and odd coherent states (also called Schrödinger cats) [5], binomial states [6,7], photon-added coherent states [8], etc. Apart from the nonclassical states that arise from the harmonic oscillator, there exist many other such states emerging from different generalizations of coherent states; see, for instance [36]. Other interesting nonclassical states have also been shown to originate from some mathematical frameworks, in particular, from the noncommutative systems reported in [37][38][39]. The striking feature of such systems is that they give rise to well-defined nonclassical states although the corresponding models are non-Hermitian.
Quite recently it has been shown that the appropriate generalizations of the oscillator algebra permit the construction of nonlinear coherent states that satisfy a closure relation which is expressed uniquely in terms of the Meijer G-function [40]. This property automatically defines the delta distribution as the corresponding P -representation. However, in the same work, it is also shown that does not exist a classical analog for such states since they lack second-order coherence and exhibit antibunching. Thus, although their P -representation is a delta function, the nonlinear coherent states studied in [40] are not classical since they are not fully coherent in the sense established by Glauber [2]. Similar results have been obtained for the para-Bose oscillator [41,42]. The algebras (13) and (14) that we have used to construct the coherent states of Section 3.2 satisfy the requirements established in [40]. Considering also that the bi-orthogonality permits to operate the non-Hermitian oscillators as in the conventional Hermitian case, it may be shown that the corresponding P -representation is proportional to the delta distribution. Nevertheless, the latter means the existence of a classical analog only if the states of Section 3.2 are fully coherent.
Instead of investigating whether our coherent states lack second (or higher) order coherence, we shall use some of the techniques described in Appendix B to study their classicalness. Namely, we are going to study the corresponding variances, Mandel parameter, Wigner distribution and purity.

Nonclassical optimized binomial states
Let us investigate the properties of the distorted quadratures (A-15) in terms of the optimized binomial states (17)- (18). The variances (∆X w ) 2 and (∆P w ) 2 may be expressed in the form (B-17) of Appendix B: where we have used (A-16) and In the above expressions Ω k,r = 1 − δ k,r δ r,0 , and Ω k,r = Ω k,r − δ k,1 δ r,0 − δ k,0 δ r,1 . Hereafter we make θ k = 0 in c k = |c b k |e iθ k . It may be shown that there is always a subset of the mean energies (19)- (20) for which (∆X w ) 2 < 1 4 I w b . Thus, it is always possible to find a mean energy E b producing the squeezing of the distorted quadrature X w . The latter is illustrated in Fig. 4 for r = 2, N = 10 and two different values of w. In turn, the squeezing of P w occurs for r = 0 only, and it is not as strong as the one suffered by X w since (∆P w ) 2 is just slightly larger than 1 4 I w b . In general, the optimized binomial states (17) are of minimal uncertainty (i.e., they are classical states) only at the lowest mean energy for either (i) r = 0 and any value of w, or (ii) r = 2 and w = 0.
In the quantum oscillator limit (6) one has Clearly, we recover the conventional binomial state for r = 0. This state is squeezed in the quadraturex, for N = 10 and p = 0.5, as shown in the Wigner distribution of Fig. 5(a).
If r ≥ 1, besides the squeezing ofx, the Wigner distribution exhibits negative values in different zones of the phase-space, see Figs. 5(b) and 5(c). The latter results confirm that the optimized binomial states (34) are stronger nonclassical than the conventional ones. To get more insights about the optimized binomial states (34) we have depicted the related Mandel parameter Q, as a function of p, in Fig. 6. For p = 1 the superposition (17) is reduced to the state ϕ N +r (x), so that all the curves shown in Fig. 6 converge to −1 as p → 1. In turn, for r = 0, all the curves take the value −1 at p = 0 since only the Fock state ϕ r =0 (x) is included in the sum. Thus, for r = 0 and any p, the Mandel parameter Q in this case is always negative. Then, the optimized binomial states (34) are sub-Poissonian for r = 0, and thus do not admit any classical description. Let us apply also the beam-splitter technique (see Appendix B) to the states we are dealing with. We may study the purity (linear entropy) S L in terms of either the transmission coefficient T of the beam-splitter or the probability p that defines the superposition (34). In the former case it is convenient to fix N and E b to parameterize the purity with r and p. Then one finds that the purity is equal to zero for T = 0 and T = 1, see Fig. 7. However, it is important to remark that such results give no special information since they correspond to complete reflectance (i.e., the beam-splitter is indeed a perfect mirror) and complete transparency (no beam-splitter), so the output is respectively the vacuum and the state which was injected in input. In general, for other values of T , the purity reaches its maximum at T = √ 0.5, i.e., when the testing beam-splitter is 50/50. To analyze S L in terms of the probability p we must fix N and T . Then, with exception of the purity for r = 0, we find S L > 0 for any value of p (see Fig. 8). The latter confirms, once again, that the optimized binomial states φ b (x) are nonclassical. Concerning the case r = 0, the purity is equal to zero at p = 0 only. Thus, the superposition (34) only includes the vacuum state ϕ 0 (x), which is clearly classical.

Nonclassical optimized Poisson states
It may be shown that the optimized Poisson states (23) are minimum uncertainty states for concrete values of the parameters like either r = 1, ω = 1, or r = 2, ω = 0. In any case, they become states of minimum uncertainty as |z| → ∞. For other values of the parameters either (∆X w ) 2 or (∆P w ) 2 is squeezed, see Fig. 9. In the quantum oscillator limit (6), the superposition (23) acquires the form of Eq. (24). As indicated above, for r = 0 this state is reduced to the conventional coherent state, so that the Mandel parameter is equal to zero, as expected. For other values of r, the optimized Poisson states (24) are sub-Poissonian. In fact, in Fig. 10 we see that the Mandel parameter becomes zero as |z| → ∞, and it is equal to −1 for r = 0 at z = 0 (i.e. at the lowest mean energy E b = 2r − 1). In other words, for r = 0 and finite values of z, the optimized Poisson states φ osc P (x) are nonclassical. The curves correspond to r = 0 (solid-red), r = 1 (dotted-blue) and r = 2 (dashed-black).
For completeness, we show the Wigner function of φ osc P (x) in Fig. 11. Clearly, for r = 0 one has a conventional (classical) coherent state. Other values of r give rise to both, the squeezing ofx and the negativity of the Wigner function in some regions of the phase-space. • Photon added coherent states. The expression (24) resembles the structure of the photon-added coherent states for the conventional harmonic oscillator [8], written in Dirac notation as In both cases the superpositions are infinite and start with the eigenstate |ϕ r . The main difference lies on the coefficients, for if we make z = α, then c (pa) k = (k + 1) r Figure 12: (Color online) The purities of the photon-added coherent state (24) and the Poisson state (35) are respectively depicted in red and black. In (a) the mean energy is fixed, E b = 9, and the purities are plotted in terms of the transmission coefficient T with r = 2 (thick-solid for (24) and solid for (35)) and r = 4 (shot-dashed for (24) and long-dashed for (35)). In (b) the transmission coefficient is fixed, T = √ 0.5, and the purities are shown in terms of z with r = 4 (shot-dashed for (24) and long-dashed for (35)) and r = 6 (thick-solid for (24) and solid for (35)). For comparison, in both figures the disk and triangle mark the purity evaluated at the same parameters.
If r = 0 we see that, up to the normalization constant e −|α| 2 /2 , the coefficients are the same. For other values of r the superpositions (24) and (35) are different in general.
Let us compare |α (pa) with |φ osc P in order to identify which one is stronger nonclassical. Using the beam-splitter technique, we find that the nonclassicality of these states depends on different parameters. For instance, in Fig. 12(a) we have fixed the value of the mean energy as E b = 9 and, after adjusting the parameters z, α, and r, the purities S L (φ osc P ) and S L (α (pa) ) have been depicted in terms of the transmission coefficient T . As a global property, we see that the purity is maximum at T = √ 0.5 in all the cases. Besides, we can appreciate that S L (φ osc P ) < S L (α (pa) ). Thus, for the parameters used in the figure the state |α (pa) is stronger nonclassical than |φ osc P . The situation changes if we fix T , and adjust the other parameters to plot the purity S L in terms of |z|. Then S L (φ osc P ) > S L (α (pa) ) holds for small values of |z| and |α|, and S L (φ osc P ) < S L (α (pa) ) is valid for large values of such parameters. The latter is illustrated in Fig. 12(b) for T = √ 0.5 and |z| = |α|.

Nonclassical natural coherent states
By construction, the natural coherent states {|ψ 0 , |z } minimize the uncertainty relation (A-9) of Appendix A to investigate their properties in the oscillator limit (6), we first calculate the purity S L using Eq. (B-23) of Appendix B with r = 1, c k = c (N ) k , and N → +∞ in the Fock basis |ϕ k . The result is depicted in Fig. 13(a). As we can see, this function is always different from zero and decreases as |z| → ∞ for any value of the transmission coefficient T . That is, the states |z osc are nonclassical. The latter is confirmed by noticing that the Mandel parameter Q is always negative, as it is shown in Fig. 13(b).
On the other hand, the variances of the physical quadraturesx andp can be expressed In Fig. 14, we show the regions of the complex z-plane where the squeezing of eitherx (grey zones) orp (mesh, dashed blue zones) occurs. Notice that the squeezing is maximum at the phase values φ = 0, π, and φ = π 2 , 3π 2 , respectively. The white zones are the regions of nosqueezing, they overrun four distinguishable areas defined along the phase values φ = kπ 4 , k = 1, 3, 5, 7, and a circle of radius |z| ≈ 5 that is centered at the origin. Comparing the Figs. 13 and 14, we see that the purity and the Mandel parameter give information that is complementary to that obtained from the variances of the quadratures. Namely, S L and Q indicate strong nonclassicality for the states |z osc that satisfy the condition |z| 5, where no squeezing is expected for neitherx norp. dashed blue) zones correspond to the values of z ∈ C for whichx (p) is squeezed. In turn, the grey and blue zones are centered at φ = 0, π and φ = π 2 , 3π 2 , respectively.
In Fig. 15 we show the density plots of the Wigner function of |z osc for different values of z. In particular, for z = 0 the natural coherent state |z osc is nonclassical because it coincides with the first excited state of the oscillator |ϕ 1 , as shown in Fig. 15(a). If |z| > 5 and φ = 0 the quadraturex is squeezed; this is illustrated in Fig. 15(b) with z = 30. In turn, for z = i30 one has the opposite result,p is squeezed; see Fig. 15(c).

Even and odd natural coherent states
We can use the natural coherent states (26) to construct additional nonclassical states. Following [5] we introduce the vectors and call them even and odd natural cat states, respectively. The related Mandel parameter Q is always negative, see Fig. 16(a), so that |z ± are nonclassical. The variances of the natural quadratures (A-8) can be expressed in the form (B-17) of Appendix B, with A = X N , B = P N , and In this case the quadratures are squeezed in very localized regions that are distributed along the real and imaginary axes of the complex z-plane. This is illustrated in Fig. 16(b), where we can appreciate a 'discretization' of the eigenvalue z as follows. The quadrature X N is maximally squeezed along the real axis, in intervals |z   1 | defining the green zones that are closest to the origin in Fig. 16(b). For the even natural cats we write j = e, with |z (e) 1 | defining the blue zones that are closest to the origin in the same figure. The intersection of the above intervals is empty, and they are such that the points |z k+1 | < · · · . The phases constraining the green and blue zones are defined in intervals φ k > 0, that are shorter as k increases. A similar description holds for the squeezing of P N along the imaginary axis.
In the oscillator limit (6), the natural cat states |z osc ± preserve their nonclassicality. In Fig. 17 the Wigner function becomes negative in diverse zones of the complex z-plane. Such regions are less evident for |z| → ∞.  The red (even) and black (odd) regions correspond to the values of z ∈ C for which P N is squeezed. Blue (even) and green (odd) regions correspond to the squeezing of X N .

Nonclassical distorted coherent states
By construction, the distorted coherent states {|ψ 0 , |z, w } are minimal uncertainty states with respect to the uncertainty relation (A-16). In the oscillator limit, the Mandel parameter Q depends on the parameter of distortion. In Fig. 18(a), we see that the states |z osc , w are sub-Poissonian for any value of |z| whenever w = 1. For arbitrary values of w > 1, they are nonclassical only in the interval (0, |z w |). Here the value of |z w | is defined by w: it is larger if the value of w is increased.
The variances of the physical quadraturesx andp can be expressed in the form given by Eq. (B-17) of Appendix B, with A =x, B =p, and where For w ≈ 1, the quadraturesx andp are squeezed along the real and imaginary axes, respectively. However, for larger values of w, they are squeezed along the imaginary and real axes. Thus, for w >> 1 the squeezing is rotated by π 2 , as shown in Fig. 18(b).

Even and odd distorted coherent states
For the even and odd distorted cats the variances of the distorted quadratures (A-15) can be expressed in the form given by Eq. (B-17) of Appendix B, with A = X w , B = P w , and In Fig. 19, we show the regions of the complex z-plane where X w and P w are squeezed for the even distorted cat |z + , w with the indicated value of w. Depending on the phase φ, the squeezing is present in zones defined by the interval 0 < |z| 7 and aligned with either the real or the imaginary axes. Such zones are of maximum width at |z| ≈ 3.5 and are slenderer as either |z| → 0 or |z| → ∞.
In the oscillator limit (6), the variances ofx andp are calculated using the Eq. (B-17)  of Appendix B, with and f 2 (|z| 2 , w) given in (41). The quadraturesx andp are squeezed in narrow zones of the complex z-plane that are defined along the real and imaginary axis, respectively. As in the previous cases, there is no squeezing for small values of |z|. Moreover, the squeezing regions are slenderer as |z| → ∞. To illustrate the phenomenon, in Fig. 20 we show the behavior of the U -functions that define the variances (∆x) 2 and (∆p) 2 of the even distorted cat |z osc + , w . As U   On the other hand, the straightforward calculation shows that neither X w nor P w are squeezed for the odd distorted cats |z − , w . The same holds for the physical quadratureŝ x andp, evaluated in the oscillator limit state |z osc − , w .
The squeezing ofx andp for |z osc , w d is quite similar to that obtained for the distorted coherent states of the previous section.
On the other hand, for the even and odd displaced coherent states |z ± , w d = 1 ± 1 F 1 (w,1,−|z| 2 ) , and h 2 (|z| 2 , w) given in (46). From Fig. 23 we see that, along the real axis, the function U is negative. The conclusion is that only for w = 1 we obtain the squeezing of either X w or P w . Remarkably, for this value of the parameter of distortion, the distorted and the displaced coherent states coincide.

Conclusions
We have studied the nonclassical properties of some pure states of the non-Hermitian Hamiltonians with real spectrum that are generated by Darboux (supersymmetric) transformations. Our interest has been focused on non-Hermitian oscillators since their eigenfunctions form bi-orthogonal bases of states in a straightforward form [24]. Depending on the parameters that define the Darboux transformation, these oscillators may be either parity-time-symmetric or non-parity-time-symmetric [13]. The pure states studied in this work are constructed as superpositions of the related bi-orthogonal basis, they include the so-called optimized Binomial and optimized Poisson states [26,30], together with the even and odd versions of the generalized coherent states that are associated with the dynamical algebras of the non-Hermitian oscillators [24,30].
One of the main results reported in this paper is to show that the techniques used to study classical like properties for Hermitian systems can be adjusted to investigate the classicality of the states of non-Hermitian systems. The latter is relevant considering that, although the non-Hermitian systems are widely used in contemporary physics [16][17][18][19][20][21][22], the classicality of the related quantum states is rarely studied. The main difficulty is that non-Hermiticity of a given operator implies that the orthonormal properties of the related eigenfunctions are not granted a priori. Since most of the approaches implemented to analyze classicality assume orthonormality of states, it is not clear how to translate such methods to the non-Hermitian case in general. Such a difficulty can be overpassed by considering a bi-orthogonal basis for the space of states of the non-Hermitian system, just as we have shown throughout this work. It is expected that our approach can be applied in the study of resonances [16] (where the non-Hermiticity of the Hamiltonian defining the interaction of a scatterer with a projectile leads to the temporal trapping of the later by the former) and in the propagation of light in materials with complex refractive index [17] (where the imaginary part of the refractive index may model absorption or amplification of electromagnetic signals).
Concerning the non-Hermitian oscillators we have shown that the corresponding biorthogonal bases converge to the Fock basis in the appropriate limit (i.e., after canceling the imaginary part of the potential and properly selecting the parameters), so that the bi-orthogonal superpositions of such states are reduced to pure states of the conventional harmonic oscillator. In particular, we have found that the generalized coherent states of the non-Hermitian oscillators are not minimum uncertainty states for the physical position and momentum quadratures in the oscillator limit. Nevertheless, they exhibit some interesting nonclassical behavior. For instance, depending on the module and complex-phase of the coherence parameter α, the physical quadratures can be squeezed for the class of generalized coherent states that we have called "natural" (see Section 3.2). The squeezing is clearly exhibited in the corresponding Wigner distribution and revealed by the Mandel parameter for which it is shown that the photon number distribution is sub-Poissonian. The other class of generalized coherent states, called "distorted", have similar quadrature squeezing properties. However, the Mandel parameter reveals a different behavior since the photon distribution becomes sub-Poissonian for some values of |α|, but it converges to a Poissonian distribution at the limit |α| → ∞. The latter implies a nonclassical-toclassical transition for the distorted coherent states (see [8] for the results on the matter for photon-added coherent states). Such a transition can be also manipulated with the parameter w that defines the related dynamical algebra.
Our method can be extended to construct and study intelligent states associated with the dynamical algebras of the non-Hermitian oscillators. The model involves finite difference equations, which deserve special treatment. Work in this direction is in progress. ∆A∆B = 1 2 | C |. If U 1 = U 2 = 0 we have two different cases (i) U 1 and U 2 are both positive, then ∆A > ∆B and we say that B is squeezed (ii) U 1 and U 2 are both negative, then ∆A < ∆B and we say that A is squeezed. This criterion is used along the paper to analyze the inequalities (A-9) and (A-16) in their respective state spaces.
On the other hand, it is well known that the Mandel parameter [47] dictates a sub-Poissonian, Poissonian and super-Poissonian photon number distribution for −1 ≤ Q < 0, Q = 0 and Q > 0, respectively. Nonclassicality corresponds to the sub-Poissonian distributions, which is associated with the squeezing of the photon number 0 ≤ (∆N ) 2 < N .
• Beam-splitter technique. The action of a beam splitter on a given state |in = |Φ 1 ⊗ |Φ 2 is that it produces non-separable outputs |out in general [11]. If |out is entangled, then the signal |Φ 1 is nonclassical, even if the ancilla |Φ 2 is a classical state. The latter criterion is used with |Φ 2 = |ϕ 0 , which is classical, and |Φ 1 being any of the bi-orthogonal superpositions at the oscillator limit (6). Thus, we may write where the super-label "osc" means that the oscillator limit (6) has been applied. The action of the beam-splitter is represented by the unitary operator [11] BS = exp θ 2 a † 1 ⊗ a 2 e iφ − a 1 ⊗ a † 2 e −iφ , (B-20) where a k and a † k are the ladder operators acting on the input states |Φ osc k , with k = 1, 2. Up to a phase, the reflection and transmission coefficients of the beam-splitter are respectively given by R = sin θ 2 and T = cos θ 2 , with θ ∈ [0, π). Using |in osc , the straightforward calculation shows that the output state is of the form where Γ k,p = e −i(φ−π)p R p T k (k + p)! k!p! c k+p−r . Classical states satisfy the separability condition S L = 0 while the maximal entanglement is obtained for S L = 1. Then, the nonclassicality is associated with 0 < S L ≤ 1.