arxiv: v1 [math-ph] 15 Dec 2012

Transcription

1 arxiv: v [math-ph] 5 Dec 202 Damping and Pseudo-fermions F. Bagarello Dieetcam, Facoltà di Ingegneria, Università di Palermo, I-9028 Palermo, Italy fabio.bagarello@unipa.it home page: Abstract After a short abstract introduction on the time evolution driven by non self-adjoint hamiltonians, we show how the recently introduced concept of pseudo-fermion can be used in the description of damping in finite dimensional quantum systems, and we compare the results deduced adopting the Schrödinger and the Heisenberg representations.

2 I Introduction In a series of papers, []-[8], we have considered two operators a and b, with b a, acting on a Hilbert space H, and satisfying the commutation rule [a,b] =. A nice functional structure has been deduced under suitable assumptions, and some connections with physics, and in particular with what is usually called quasi-hermitian quantum mechanics and with the technique of intertwining operators, have been established. Following Trifonov, [9], we have called pseudo-bosons PB) the particle-like excitations associated to this structure. A similar analysis has also been carried out for what we have called nonlinear pseudo-bosons NLPB) in [0]-[2], and most of the original results have been recovered also in this more general situation, where the main ingredient is not the commutation rule between a and b but their raising and lowering properties, when applied to two fixed biorthogonal bases. The analytical treatment of both PB and NLPB turns out to be particularly difficult in the case where regularity is lost, that is, see below, when these bases are not Riesz bases. In this case, in fact, the intertwining operators appearing in the game whose square roots are metric operators in the sense of the literature on quasi-hermitian quantum mechanics, see [3, 4, 5] and references therein) turns out to be unbounded. For this reason, a large amount of mathematical care is required, and this makes the rigorous treatment of the physical system rather complicated, [2]. More recently, [6], we have introduced a similar extension for the canonical anticommutation relation CAR), following an original idea by Trifonov and his collaborators, [7], also briefly considered in [8] and in [9]. The idea is very similar to the previous one: we consider again two operators a and b, with b a, acting on a Hilbert space H, and satisfying the anticommutation rule {a,b} =, as well as a 2 = b 2 = 0. Of course, if b = a we go back to ordinary CAR. The functional structure that we deduce out of these operators is nice and friendly, since H is finite-dimensional. The related excitations are called pseudo-fermions, PF). One of the most important features of PF is that all the problems arising when dealing with PB, because of their in general) unbounded nature, are now absent. In this paper we focus on a particular aspect of PF, i.e. on their possible use in the analysis of damped quantum systems. More in details: in the description of damping in quantum optics certain finite dimensional, non self-adjoint, matrices are quite often used. These matrices are usually called effective hamiltonians. Solving the related Schrödinger equation, one sees that the wave-function of the system goes to zero for large t, at least if the parameters of the model are chosen properly. Already in [6] we have seen that one such a model, originally proposed in [20], can be easily written in terms of PF. However, we have not explored the damping 2

3 features of that system, because that was not our major interest. Here, on the other hand, we focus our attention exactly on this aspect, and we work both in the Schrödinger and in the Heisenberg representations, showing that analogous conclusions can be deduced. We also extend our analysis to a higher-dimensional Hilbert space, which we will relate to a two-dimensional family of PF. For both these examples we deduce damping. The paper is organized as follows: in the next section we briefly discuss some basic facts on PF and on the time evolution in the Heisenberg representation, when this is driven by some non self-adjoint hamiltonian. This, we believe, is useful since not all the readers are familiar with this kind of time-evolution, and surely less readers are familiar with pseudo-fermions. Section III is devoted to our first, two-dimensional, example, while in Section IV we consider a four-dimensional, example. For both these examples we compare the well known Schrödinger point of view with the modified Heisenberg evolution, and we show that, not surprisingly, the same conclusions can be deduced. Particularly relevant for us is the fact that, after a somehow standard transformation of the effective hamiltonian, we will naturally driven to introduce pseudo-fermions in the game. Section V contains an abstract generalization of the approach, while our conclusions are given in Section VI. II Pseudo-fermions and dynamics We begin this section by briefly reviewing the main definitions and results concerning PF in one dimension. The extension to higher dimensions will be discussed later on. The starting point is a modification of the CAR {c,c } = cc +c c =, {c,c} = {c,c } = 0, between two operators, c and c, acting on a two-dimensional Hilbert space H. The CAR are replaced here by the following rules: {a,b} =, {a,a} = 0, {b,b} = 0, 2.) where the interesting situation is when b a. These rules automatically imply that a non zero vector, ϕ 0, exists in H such that aϕ 0 = 0, and that a second non zero vector, Ψ 0, also exists in H such that b Ψ 0 = 0, [6]. Under these conditions it is possible to recover similar results as those for PB. In particular, we first introduce the following non zero vectors ϕ := bϕ 0, Ψ = a Ψ 0, 2.2) It might be worth stressing that we are not claiming here that PF are always relevant to explain damping. 3

4 as well as the non self-adjoint operators N = ba, N = N = a b. 2.3) We further introduce the self-adjoint operators S ϕ and S Ψ via their action on a generic f H: S ϕ f = ϕ n,f ϕ n, S Ψ f = Ψ n,f Ψ n. 2.4) n=0 n=0 Hence we get the following results, whose proofs are straightforward and will not be given here:. 2. aϕ = ϕ 0, b Ψ = Ψ ) Nϕ n = nϕ n, NΨ n = nψ n, 2.6) for n = 0,. 3. If the normalizations of ϕ 0 and Ψ 0 are chosen in such a way that ϕ 0,Ψ 0 =, then for k,n = 0,. ϕ k,ψ n = δ k,n, 2.7) 4. S ϕ and S Ψ are bounded, strictly positive, self-adjoint, and invertible. They satisfy S ϕ ϕ ϕ 2, S Ψ Ψ Ψ 2, 2.8) S ϕ Ψ n = ϕ n, S Ψ ϕ n = Ψ n, 2.9) for n = 0,, as well as S ϕ = S Ψ. Moreover, the following intertwining relations hold. S Ψ N = NS Ψ, S ϕ N = NS ϕ, 2.0) The above formulas show that i) N and N behave as fermionic number operators, having eigenvalues 0 and ; ii) their related eigenvectors are respectively the vectors of F ϕ = {ϕ 0,ϕ } and F Ψ = {Ψ 0,Ψ }; iii) a and b are lowering operators for F ϕ and F Ψ respectively; iv) b and a are rising operators for F ϕ and F Ψ respectively; v) the two sets F ϕ and F Ψ are 4

5 biorthonormal; vi) the well-behaved 2 operators S ϕ and S Ψ maps F ϕ in F Ψ and viceversa; vii) S ϕ and S Ψ intertwine between operators which are not self-adjoint, in the very same way as they do for PB. We refer to [6] for further remarks and consequences of these definitions. In particular, for instance, it is shown that F ϕ and F Ψ are automatically Riesz bases for H, and the relations between fermions and PF are discussed. II. The algebraic dynamics In [20] a non self-adjoint effective hamiltonian is considered in the description of the time evolution of the wave-function of a two-dimensional quantum system. The resulting dynamics shows a decay of this wave-function. This is just a single reference where such an approach is considered. Other references are, for instance, [7], [2] and, for more abstract considerations, [22]. The starting point is always the same: a Schrödinger equation i Ψt) = H eff Ψt), 2.) where H eff H eff. Actually, in [20] and [7], H eff is assumed to be pseudo-hermitian. We will not consider this extra assumption in the first part of this section. It is worth stressing here that formula 2.) should be considered as an assumption in itself, since in standard quantum mechanics the Schrödinger equation is usually if not always) assumed to hold for a self-adjoint hamiltonian. Of course, if H eff does not depend explicitly ontime, the formal solution of 2.) is Ψt) = e ihefft Ψ0). It is clear that, if H eff is a finite-dimensional matrix, this exponential can be defined using, for instance, a norm convergent series: e ihefft it) k :=. k! k=0 H k eff. On the other hand, if H eff is an unbounded operator, defining e ihefft is quite a delicate task, since we cannot even use the spectral theorem. We will not consider this aspect of the theory here, since we are concerned only with finite matrices. Solving 2.) means that we are adopting the Schrödinger representation. To move to the Heisenberg representation, we assume, as it is always done in ordinary quantum mechanics, that the time evolution of the mean values of the observables do not depend on the representation chosen. In other words, 2 i.e. self-adjoint, bounded, invertible, with bounded inverse 5

6 if X is an observable of our physical system, calling X eff t) its time evolution driven by H eff, we require that Ψt),XΨt) = Ψ0),X eff t)ψ0). Since this equality must hold for all possible choices of Ψ0), we deduce that X eff t) := e ih eff t X e ih eff t. 2.2) This isthetime evolution of X inthe Heisenberg representation. Notice that, if H eff = H eff =: H, we go back to the usual formula and to the unitary evolution, Xt) := e ih t X e ih t. Notice also that, while XY)t) = Xt)Yt), XY) eff t) X eff t)y eff t): the time evolution in 2.2) is not an automorphism of the set of the observables of the system. Another interesting feature of 2.2) is that this time evolution is stable under the adjoint: X eff t)) = X ) eff t). Also, it is clear that the Heisenberg equation of motion for X eff t) should be modified. For that, let us define the following effective commutator between two generic operators A and B: [A,B] eff := AB B A. 2.3) The effective commutator satisfies, for example, the following equality: given three generic operators A, B and C, then [AB,C] eff = A[B,C] eff +[A,C] eff B +AC C)B, which reduces to the standard formula when C = C, but not in general. Also, not surprisingly, the Jacobi identity ceases to be true for [.,.] eff. Still we have [A,B] eff = [A,B] eff. The differential equation for X eff t) is the following: i d dt X efft) = [X eff t),h eff ] eff = e ih eff t [X,H eff ] eff e ih eff t, 2.4) which returns the standard equation if H eff = H eff. This formula shows, among other things, that an integral of motion, Z, is not an operator commuting with H eff. What Z has to satisfy is clearly the following requirement: [Z,H eff ] eff = 0. In fact, if this is true, then Żefft) = 0. In [2] the authors have considered a slightly different ) point of view, ) splitting H eff into two parts, both self-adjoint: H eff = H 2 eff +H eff i H 2i eff H eff. In this way, rather than introducing an effective commutator, they deduce a differential equation for X eff t) in which both a commutator and an anti-commutator appear. 6

7 Something more can be said if H eff is crypto-hermitian, see [], i.e. if a positive selfadjoint, time-independent, operator Θ exists, bounded with bounded inverse, such that H eff = Θ H effθ. In this case we get the following: and the differential equation is Θ X eff t) = e ih efft Θ X ) e ih efft, d Θ X eff t) ) = ie ihefft [H eff,θ X]e ihefft. dt Let us remark that here, first of all, the commutator is the standard one. Also, at a first sight the time evolution for Θ X might appear as an automorphism. Nevertheless, it is not hard to check that this is not true, exactly because of the presence of Θ. In fact, Θ X eff t) is, in general, different from Θ X) eff t). A simple case when they coincide is when H eff is self-adjoint. The essential reason why formulas above look interesting is because only H eff appears, and not H eff. In particular, the differential equation shows that X efft) is constant when H eff commutes not with X but with Θ X. This is the formula which replaces, in this particular case, the one deduced above, [X,H eff ] eff = 0. Because of the properties of Θ, we know that Θ ±/2 can be defined, and h := Θ /2 H eff Θ /2 is a self-adjoint operator: h = h. Let us now introduce the following automorphism of BH), the set of the bounded operators on H: j Θ x) := Θ /2 xθ /2. This map is invertible, and it is clear that h = j Θ H eff). If we further define the following standard Heisenberg evolution, α t X) := e iht Xe iht, X BH), which has all the usual properties of the time evolution for systems driven by self-adjoint hamiltonians, we deduce that X eff t) = j Θ α t j Θ X)) ), 2.5) which has a first interesting consequence: this effective time evolution coincides with α t if this latter commutes with j Θ. Moreover, the interpretation is the following: to find X eff t), we have first to map X into j Θ X). Then we let this operator evolve using h, and finally we transform back this result, using j Θ ; therefore, the effective time evolution appears to be the composition of these three maps. 7

8 III An example from the literature The first example is taken by [7, 20], and its relation with PF was already partly considered in [6]. The starting point is the Schrödinger equation ) i Ψt) = H eff Ψt), with H eff = iγ a v 2 v iγ b, 3.) where γ a,γ b > 0 and v C. III. Schrödinger representation With a simple change of variable Φt) = e Γt Ψt), Γ = 2 γ a+γ b ), we deduce that i Φt) = HΦt), where H = iγ 2 +H eff = iγ v v iγ ), Φt) = Φ 0 t) Φ t) Here 2 is the two-by-two identity matrix, and γ = 2 γ a γ b ). The differential equations for Φ 0 t) and Φ t) are easily deduced. Let us introduce Ω := v 2 γ 2. Then { Φ0 t) = ΩΦ 0 t), Φ t) = ΩΦ t). ). 3.2) If Ω = 0 then the functions Φ 0 t) and Φ t) are linear in t, so that ) ) Ψt) = e Γt Φ 0 t) e γa+γ b) 2A t 0 +B 0 t) =, Φ t) e γa+γ b) 2A t +B t) where A 0,A,B 0 and B are fixed by the initial conditions. In particular, if Ψ0) = we find that A 0 = ϕ 0, A = ϕ, B 0 = γϕ 0 ivϕ and B = γϕ ivϕ 0. The norm of Ψt) is simply the norm in H = C 2. Therefore Ψt) 2 = e tγa+γ b) A 0 +B 0 t 2 + A +B t 2), ϕ 0 ϕ ), and this squared norm goes to zero when t. Let us consider now the case in which Ω > 0. In this case the solution can be written as Ψt) = e γa+γ A b) t 0 cos Ωt)+B 0 sin ) Ωt) 2 A cos Ωt)+B sin, Ωt) 8

9 where, again, A 0,A,B 0 and B are fixed by the initial conditions: A 0 = ϕ 0, A = ϕ, B 0 = Ω γϕ 0 ivϕ ) and B = Ω γϕ ivϕ 0 ). In this case Ψt) 2 is the product of a decaying exponential and an oscillating function. More explicitly, Ψt) 2 = e tγa+γ b) A 0 cos Ωt)+B 0 sin Ωt) 2 + A cos Ωt)+B sin ) Ωt) 2, which decays to zero independently of the particular values of ϕ 0 and ϕ. Let us finally consider the case in which Ω < 0. In this case the solution can be written as Ψt) = e γa+γ A b) t 0 exp Ω t)+b 0 exp ) Ω t) 2 A exp Ω t)+b exp. Ω t) Here, A 0,A,B 0 and B are related to the initial conditions as follows: ) A 0 = ϕ 0 + γϕ 0 ivϕ ), B 0 = 2 Ω 2 ) A = ϕ + γϕ ivϕ 0 ), B = 2 Ω 2 ϕ 0 + Ω γϕ 0 +ivϕ ) ϕ + Ω γϕ +ivϕ 0 ) In this case it is not completely clear that Ψt) goes to zero for increasing t. However, to get such a conclusion, it is enough to consider the asymptotic behavior of e γa+γ b)t e 2 Ω t, which is the only dangerous contribution appearing in Ψt) 2. The conclusion is that Ψt) 0, when t, if 2 Ω < γ a +γ b, which is always satisfied. The outcome of this analysis is, therefore, the following: the wave-function Ψt) solving the Schrödinger equation 3.) goes to zero independently of the value of Ω. ) ),. III.2 Heisenberg representation It is easy to see that Ω is related to the eigenvalues of H, which are λ ± := ± v 2 γ 2. Hence λ + = λ = 0 if Ω = 0, while λ ± := ± Ω when Ω > 0. Finally, λ ± := ±i Ω when Ω < 0. In all these cases the eigenstates of H can be written as iγ + ) ) Ω v η + =, η = iγ + ) Ω v In fact, Hη ± = λ ± η ±. Moreover, in general η + and η are not orthogonal: η +,η = 2γ γ i ) Ω, which is zero only if γ = 0 H = H ) or if γ = i Ω H = H ). Due to v 2 the relation between H and H eff we also find that H eff η ± = E ± η ±, E ± = i 2 γ a +γ b )± Ω. 3.3) 9 ).

10 In [6] we have analyzed, in connection with PF, the following non self-adjoint operator, ) H = iδ ω, 2 ω iδ which coincides with our H simply taking δ = 2γ and ω = 2v. As in [6], it is therefore possible to introduce two operators a and b, such that {a,b} =, a 2 = b 2 = 0, and H = Ω ba ). 2 Introducing also the number operators N = ba and N = N, and assuming, just to fix the ideas, that Ω > 0, we can write H = Ω N ) and H = Ω N ). Now, we want to 2 2 recover here the same damping we have found working in Schrödinger representation. For that, it is natural to consider the time evolution of the number operators N and N, which should be computed as in 2.2): N eff t) = e ih eff t N e ih eff t, N eff t) = e ih eff t Ne ih eff t = N eff t)). Due to the relation between H eff and H, H eff = H iγ, and to the above expressions for H and H, we find that, for each observable X, X eff t) = e 2Γt e iωnt Xe iωnt. 3.4) Therefore, recalling that the anticommutation rules imply that N 2 = N, N 2 = N, and so on, we find that N eff t) = e 2Γt Ne iωt +NN e iωt ) ), while N eff t) is simply its adjoint. Then, if we estimate the norm of N eff t), it is trivial to deduce that N eff t) 3e 2Γt, which goes to zero when t diverges. Hence, as expected, we recover damping also in Heisenberg picture, with no need of introducing extra assumptions on the parameters of the system. More connections with PF are discussed in [6]. IV A two-dimensional example We are going now to repeat the same analysis considering the following four-by-four effective hamiltonian: 0 αβω ω 2 ) 0 0 H eff = i ω ω 2 α+β)ω ω 2 ) 0 0 ω 2 αω 2 βω 2 αω 0 βω 2 β 2 ω 2 0 αω 2 βω where α β and ω > ω 2 > 0. It is evident that H eff H eff. We will show that, also for this hamiltonian, a decay can be deduced. 0,

11 IV. Schrödinger representation As before, we want to solve first the Schrödinger equation i Ψt) = H eff Ψt) and check whether damping is deduced, and in which conditions. Even if H eff has no direct physical interpretation, still it is interesting, since, with a simple transformation, we can transform the original Schrödinger equation into a new one, with a traceless hamiltonian H, which is also not selfadjoint. This is very close to what was done, for instance, in [7, 6, 20], and in the previous section. For that it is enough to proceed as in the previous example, defining Φt) = e Γt Ψt), H = iγ 4 +H eff, and Γ = α+β 2) ω ω 2 ), where 4 is the identity matrix in H = C 4. Then we get where H = i and i Φt) = HΦt), 2 α+β)ω ω 2 ) αβω ω 2 ) 0 0 ω ω 2 2 α+β)ω ω 2 ) 0 0 ω 2 αω 2 2 )ω +ω 2 ) 0 βω 2 β 2 ω )ω +ω 2 ) 4.) Φ 0 t) Φ t) Φt) = Φ 2 t). Φ 3 t) As in our first example, the matrix H is traceless. The related Schrödinger equation of motion produces, first of all, Φ 0 t) = ΩΦ 0 t), where Ω := 2 ω ω 2. Then Φ 0 t) = Ã0e Ωt + B 0 e Ωt, with Ã0 and B 0 fixed by the initial conditions. The second component of Φt), Φ t), can be deduced from Φ 0 t) since Φ t) = Moreover { αβω ω 2 ) ) 2 α+β)ω ω 2 )Φ 0 t) ) Φ 0 t). Φ 2 t) = A 2 e ω +ω 2 )t/2 +Ã2e Ωt + B 2 e Ωt, Φ 3 t) = A 3 e ω +ω 2 )t/2 +Ã3e Ωt + B 3 e Ωt. As for the asymptotic behavior of Ψt), we recall that ω > ω 2 > 0. Of course, to conclude that Ψt) = e Γt Φt) 0 when t, it is enough to consider the worse contribution, i.e. the one 4.2),

12 coming from Φ 3 t): in fact, if this goes to zero, then it is clear that all the other components go to zero as well, so that Ψt) 0 when t. We have, for t, { ) } { } t Ψ 3 t) exp 2 ω +ω 2 ) Γ t = exp αω 2 βω ). Recalling that α > β, Ψ 3 t) 0 if and only if αω 2 βω < 0, i.e. if α β < ω ω ) In other words, when condition 4.3) is satisfied, not only Ψ 3 t), but Ψt) itself decreases to zero for t diverging. IV.2 Heisenberg representation Let us now consider the following operators β 2 β 3 0 β a := β β 2 0 α α 2 β β α 2 +β) 0, b := a 2 := β β α) α β 2 α) β 2 β α) α α β 3 0 β α) α, b 2 := α) α 0 0 α α 2 +β α) β α) α) α) α α 2 β β β α 0 β 0 α 0 α 0 α αβ 0 0 They satisfy the following anticommutation rules: {a j,b k } = δ j,k 4, with {a j,a k } = {b j,b k } = 0, j,k =,2 3. Therefore they are two-dimensional pseudo-fermions. H can be written in terms of these operators as H = i ω b a +ω 2 b 2 a 2 ω +ω ). As in the previous example, the time evolution of a given operator X, can be written in a way which extends 3.4):,. X eff t) = e ih eff t X e ih efft = e 2Γ+ω +ω 2 )t e ω N +ω 2 N 2 )t Xe ω N +ω 2 N 2 )t, 3 Other examples of 4 4 matrices satisfying these rules can be found in [23] 2

13 where N j = b j a j and N j = N j, j =,2. Recalling that [N,N 2 ] = 0 and that Nj 2 = N j, we get N,eff t) = e 2Γ+ω 2)t N +N N e ωt ) ) 4 +N 2 +N 2 )e ω2t )+N 2 N 2 e ω2t ) ), and a similar expression can be deduced for N 2,eff t). Moreover, N j,eff t) is simply the adjoint of N j,eff t). Using N j and N j, j =,2, we get N,eff t) exp{ 2Γ ω ω 2 )t}, which goes to zero, when t, if and only if 2Γ > ω + ω 2, i.e., if condition 4.3) is satisfied: we recover exactly and not surprisingly) the same conclusion as in the Schrödinger representation. Remark: It might beinteresting to notice that the matrices a j andb j above can be related to the following, standard, two-dimensional fermion annihilation operators, A = 0 0 0, A = , as in [6]: a j = TA j T, b j = TA j T, j =,2, where T is the following non-singular matrix: 0 α β T = α β 0 α Here α 0 and α β. Therefore, this example can be considered as an explicit realization of the equivalence theorem stated in [6], concerning the relations between fermions and pseudofermions. V An abstract generalization We will now describe how the examples considered in the previous sections can be further extended. Let us consider two sets of pseudo-fermionic operators, a j and b j, with b j a j, j =,2,...,N, satisfying {a j,b k } = δ j,k, {a j,a k } = {b j,b k } = 0, 3

14 for all j and k, and let ϕ 0,0,...,0 and Ψ 0,0,...,0, with ϕ 0,0,...,0,Ψ 0,0,...,0 =, be the vacua of the a j s and b j s respectively: a j ϕ 0,0,...,0 = b j Ψ 0,0,...,0 = 0, for all j. Then, extending the procedure of Section II, the sets F ϕ and F Ψ of the functions ϕ n := ϕ n,n 2,...,n N = b n bn 2 2 bn N N ϕ 0,0,...,0, and Ψ n := Ψ n,n 2,...,n N = a ) n a 2) n2 a N )n N Ψ 0,0,...,0, n j = 0,, are biorthonormal bases of the 2 N -dimensional Hilbert space H N. Let now i Ψt) = H eff Ψt), with H eff H eff, be our original Schrödinger equation. As in Sections III and IV, we introduce a new wave function, Φt) = e Γt Ψt), where Γ is a real constant, to be fixed. The Schrödinger equation for Φt) is i Φt) = H N Φt), H N := H eff +iγ N, where N is the identity operator on H N. The examples discussed before show that, in some interesting situations, the hamiltonian H N has the following general form: H N = N Ω j N j 2 N Ω j N, j= j= where N j = b j a j is the j th pseudo-fermionic number operator, while the Ω j s are, in general, complex quantities. Let us now define the following quantity: T N := k,n Ψ n,h N ϕ k = k Ψ k,h N ϕ k, 5.) where the last equality is a consequence of the orthonormality of F ϕ and F Ψ, and of the fact that ϕ k is eigenstate of H N. T N appears a natural extension of the trace of an operator to the present settings, where a single orthonormal basis is replaced by two biorthonormal bases. From 5.) it is now easy to check that T N = 0 for all N =,2,... In the extended) Heisenberg representation, see Section II., the time evolution X eff t) = e ih eff t Xe ihefft of a given operator X, assume the following form: X eff t) = e t2γ+iω )+ +IΩ N )) N e itω j N j X N e itω ln l. 5.2) j= l= 4

15 Focusing now to the time behavior of the number operator N k, and recalling that N 2 j = N j and N 2 j = N j, for all j, we deduce that N k ) eff t) = e t2γ+iω )+ +IΩ N )) N j= N N +N j e itω j ))N k l= N +N l e itω l )) Damping can be found by estimating the norm of the various N k ) eff t). For that, we first recall that N j = N j, for all j. The simplest situation is when the complex frequencies are all real: in fact, if IΩ j ) = 0 for all j, we easily conclude that N k ) eff t) 3 2N e 2Γt, which goes to zero when t diverges for all possible values of positive Γ. Hence, in this case, damping is deduced. It is probably more interesting to consider the situation in which not all the Ω j s are strictly real. In particular, we consider here the opposite situation, i.e. the one in which IΩ j ) 0 for all j, since this is the most dangerous case. Now, the estimate above should be replaced by the following: N k ) eff t) e 2Γ+IΩ )+ +IΩ N ))t N ) 2+e tiω j ) 2, which, clearly, extends the previous one. In this case, damping is recovered if j= 5.3) Γ > 2 N IΩ j ), 5.4) j= while a similar and simpler) condition is recovered if only some of the Ω j s are not purely real. Notice that the inequality 5.4) becomes quite simple, Γ > 0, when IΩ j ) = 0, for all j. In this case, we recover our previous conclusion. If, inanalogywitht N, weintroducethegeneralizedtracealsoforh eff, recallingthatt N = 0 we find that T eff N := k Ψ k,h eff ϕ k = i2 N Γ. Therefore, in order for H eff to produce damping, the following inequality should be satisfied by the effective hamiltonian H eff : it eff N > 2N j= N IΩ j ). 5

16 We end this section with some remarks on the pseudo-fermionic structure related to this abstract system. The starting point is the pair of operators S ϕ and S Ψ which extends those of 2.4). They satisfy a multi-dimensional version of 2.0), N j S ϕ = S ϕ N j, N j S Ψ = S Ψ N j, for all j. Moreover, since the square root of both S ϕ and S Ψ can be defined, we can also introduce the self-adjoint operators n j = S /2 Ψ N jsϕ /2 = n j. Therefore, if the Ω j are real, we deduce that S Ψ H N = H N S Ψ and that h := S /2 Ψ H N S /2 ϕ = N Ω j n j 2 j= N Ω j N 5.5) j= is a self-adjoint operator. It is well known that the operator S /2 Ψ, as well as its inverse S/2 ϕ, can be used to define a different scalar product in H N. Then the conclusion is the following: a damping effect in H N, endowed with its natural scalar product <.,. >, can be obtained using certain non self-adjoint effective hamiltonians which, when considered again in H N, but endowed with a different scalar product, are self-adjoint. In fact, let us define a new scalar product.,. Ψ on H N as follows: f,g Ψ := S /2 Ψ f,s/2 Ψ g, for all f,g H N. Then, since the intertwining relation above implies that H N = S Ψ H N S Ψ, we get f,h N g Ψ = = S /2 Ψ f,s/2 Ψ H N g = S /2 Ψ H N f,s /2 Ψ g = H N f,g Ψ, S Ψ S Ψ H N S Ψf,g = for all f,g H N. Hence H N turns out to be self-adjoint with respect to.,. Ψ. Alternatively and more straightforwardly), formula 5.5) shows that the effective hamiltanians considered here and producing damping are similar but not unitarily equivalent) to self adjoint operators. A natural and quite interesting question to consider is therefore: is this a general requirement? We have no answer now, but this is part of our works in progress. VI Conclusions After some general remarks on the time evolution driven by a non self-adjoint hamiltonian, we have shown how this can be used in the analysis of the decay of finite dimensional systems 6

17 in connection with pseudo-fermions. In particular, in the examples discussed here, we have seen that the effective hamiltonians describing damping are related, in a very direct way, to pseudo-hermitian operators. We have also shown that the Heisenberg representation can be conveniently adopted in the analysis of the time evolutions of these systems. The next natural step of our research will be to check how much of this approach can be extended to infinite dimensional spaces and to try to answer to the following, rather general, question: can any effective damping be described in terms of some pseudo-hermitian hamiltonian? Acknowledgements The author acknowledges financial support by the MIUR. References [] F. Bagarello, Pseudo-bosons, Riesz bases and coherent states, J. Math. Phys., 50, DOI:0.063/ , ) 0pg) [2] F. Bagarello Construction of pseudo-bosons systems, J. Math. Phys., 5, doi:0.063/ , ) 0pg) [3] F. Bagarello Mathematical aspects of intertwining operators: the role of Riesz bases, J. Phys. A, doi:0.088/75-83/43/7/75203, 43, ) 2pp) [4] F. Bagarello, F. Calabrese Pseudo-bosons arising from Riesz bases, Bollettino del Dipartimento di Metodi e Modelli Matematici, 2, 5-26, 200) [5] F. Bagarello, Regular) pseudo-bosons versus bosons, J. Phys. A, 44, ) [6] F. Bagarello, Examples of Pseudo-bosons in quantum mechanics, Phys. Lett. A, 374, ) [7] S.T. Ali, F. Bagarello, J.-P. Gazeau, Modified Landau levels, damped harmonic oscillator and two-dimensional pseudo-bosons, J. Math. Phys., 5, ) [8] F. Bagarello, Pseudo-bosons, so far, Rep. Math. Phys., 68, No. 2, ) 7

18 [9] D.A. Trifonov, Pseudo-boson coherent and Fock states, arxiv: quant-ph/ , Differential Geometry, Complex Analysis and Mathematical Physics, eds. K. Sekigawa et al W. Scientific 2009), pp [0] F. Bagarello, Non linear pseudo-bosons, J. Math. Phys., J. Math. Phys., 52, 06352, 20) [] F. Bagarello, M. Znojil, Non linear pseudo-bosons versus hidden Hermiticity, J. Phys. A, , 20) [2] F. Bagarello, M. Znojil, Non linear pseudo-bosons versus hidden Hermiticity. II: The case of unbounded operators, J. Phys. A, 45, 53 doi:0.088/75-83/45//53, 202) [3] C. Bender, Making Sense of Non-Hermitian Hamiltonians, Rep. Progr. Phys., 70, ) [4] A. Mostafazadeh, Pseudo-Hermitian representation of Quantum Mechanics, Int. J. Geom. Methods Mod. Phys. 7, ) [5] M. Znojil, Three-Hilbert-space formulation of Quantum Mechanics, SIGMA 5, ) [6] F. Bagarello, Linear pseudo-fermions, J. Phys. A, in press [7] O. Cherbal, M. Drir, M. Maamache, D. A. Trifonov, Fermionic coherent states for pseudo- Hermitian two-level systems, J. Phys. A, 40, , 2007) [8] A. Mostafazadeh, Pseudo-Hermitian representation of Quantum Mechanics, J. Phys. A, , 2004) [9] C. M. Bender, S. P. Klevansky, PT-symmetric representations of fermionic algebras, Phys. Rev. A, 84, ) [20] Y. Ben-Aryeh, A. Mann, I. Yaakov, Rabi oscillations in a two-level atomic system with a pseudo-hermitian hamiltonian, J. Phys. A, , 2004) [2] E.-M. Graefe, M. Hoening, H. J. Korsch, Classical limit of non-hermitian quantum dynamics-a generalized canonical structure, J. Phys. A, 43, , 200) [22] F. Bagarello, M. Znojil, The dynamical problem for a non self-adjoint hamiltonian, Proceedings of IWOTA 20 July 26, 200, Berlin, Germany), Operator Theory: Advances and Applications, 22, 09-9, 202 8

19 [23] O. Cherbal, D. A. Trifonov, Extended PT- and CPT-symmetric representations of fermionic algebras, Phys. Rev. A 85, ) 9