arxiv: v1 [quant-ph] 25 Sep 2015

Transcription

1 Neutral kaons as an open quantum system in a second quantization approach Kordian Andrzej Smoliński arxiv: v quant-ph] 5 Sep 05 Department of Theoretical Physics, Faculty of Physics Applied Computer Sciences, University of odz, ul. Pomorska 49/53, ódź, Pol (Dated: September 4, 05) Abstract We have shown that it is possible to formulate the consistent probability-preserving description of the CP-symmetry-violating evolution of a system of decaying particles. This has been done within the framework of quantum mechanics of open systems. This approach allows the description of both the exponential decay flavour oscillations. We have solved explicitly the Kossakowski-indblad master equation for a system of particles with violated CP symmetry found the evolution of any observable bilinear in creation annihilation operators. The choice of a concrete observable can be done by the proper choice of initial conditions for the system of differential equations. We have calculated the evolution as well as mean values of the observables most interesting from the physical point of view, we have found their lowest order difference with the CP-preserved values. PACS numbers: Ta, Yz,.30.Er, 4.40.Df K.A.Smolinski@merlin.phys.uni.lodz.pl

2 I. INTRODUCTION In recent years Bell inequalities ] were tested in systems of correlated neutral K, 3] or B mesons 4]. However, the quantum mechanical analysis of such systems encounters an important difficulty, namely, the irreversibility of time-evolution of unstable particles. Only the complete system consisting of the decaying particle as well as the decay products undergoes unitary evolution, which is actually described by quantum field theory. But in Einstein-Podolsky-Rosen (EPR) correlation experiments 5], it is more useful to consider the decaying particles only, neglecting the evolution of decay products. In the usual approach one introduces a non-hermitian Hamiltonian, as it was done in the classical works of Weisskopf Wigner 6]. However, the non-hermitian Hamiltonian does not provid an unambiguous way of calculating the probability of finding the system consisting of a few such particles in a given state after the measurement (such an approach causes also other ambiguities, see e.g. 7] for a discussion). Indeed, using this approach, we must reinterpret quantum mechanical results with the help of probability theory, especially if we would like to analyze correlation experiments with unstable particles, namely, neutral K or B mesons 8]. From the quantum mechanical point of view, this means that unit trace positivity of the density matrix for the system under consideration must be preserved 9]. Fortunately, it is possible to resolve the above-mentioned problems using the open quantum systems theory 0]. The idea that unstable particles can be treated as open quantum systems was proposed first by R. Alicki ] was developed by various authors in different contexts (see e.g., 3]). Recently, this approach was successfully applied to the systems of particles with flavour oscillations (like in the case of neutral kaons or B mesons) 4 6], it has been used successfully in the description of EPR correlations evolution of entanglement in K 0 K 0 system 7, 8]. Inthis theory theevolution of a quantum system is described by the master equation 9], which can be treated as the replacement for either von Neumann or Heisenberg equations, for the picture of quantum mechanics in use. Here, we follow the approach presented in 6], where systems with an arbitrary number of particles were described with the use of the second quantization formalism, which seems to be the most natural language for a system with varying number of particles. The paper is organized as follows. In Sect. II we introduce the notation conventions used through the paper. In Sect. III we introduce the master equation for a system of neutral kaons with violated CP symmetry we find the evolution of the observables in the Heisenberg picture,, in Sect. IV, we analyze the evolution of the total number of particles strangeness observable as well as the numbers of each of neutral kaon flavor. Finally, we conclude the paper in Sec. V. II. PREIMINARIA We start with introducing some notations, conventions definitions used throughout the paper, as well as give experimental values of some parameters important in the description of the systems of K 0 K 0. We define one-particle states K 0 K 0 as K 0 = a 0, (a) K 0 = b 0, (b)

3 where a b fulfills the usual canonical commutation relations a,a ] = b,b ] =, a,b] = a,b ] = 0. (a) (b) States K 0 K 0 are orthonormal, i.e. K 0 K 0 = K 0 K 0 =, (3a) K 0 K 0 = 0, (3b) they are eigenstates of the strangeness operator S = a a b b: S K 0 = K 0, S K 0 = K 0. (4a) (4b) For our convenience we will use the following short-h notation for multiparticle states: #K 0 = n,#k 0 = n n, n. (5) Therefore n, n = n! n! ( a ) n( b ) n 0 (6) S n, n = (n n) n, n. (7) However, the states K 0 K 0 are not eigenstates of time evolution. Actually, under time evolution the well behaving states are K 0 S K0, defined as follows K 0 S = p K 0 +q K 0, (8a) K 0 = p K 0 q K 0, (8b) with mean life times τ S = s τ = s, respectively 0]. Moreover, if p q then CP symmetry is violated. Indeed, p q = + ǫ ǫ (9) with p + q =, p q = A, (0a) (0b) A = Re(ǫ)/(+ ǫ ) is a measure of violation of CP symmetry (hereafter we use the so-called Wu-Yang phase convention ] in which ǫ = ǫ). Experimentally obtained values are A = 0.33% ǫ = ]. Notice also that the basis (8) is not orthogonal since K 0 S K0 = A. () 3

4 III. MASTER EQUATION The aim of this section is to find the evolution of the multiparticle system of K 0 K 0. We know that the evolution of a single neutral kaon can be described as an open system obeying the Kossakowski-indblad master equation 5]. In this case the one-particle Hamiltonian for this master equation was written in the base (8) reads H = A { ms K 0 S K0 S +m K 0 K0 A ( m i 4 Γ) K 0 S K0 +( m+ i 4 Γ) K 0 K0 S ]}, (a) while indbladians are of the form = Γ A S A Γ + m ( 0 K 0 Γ S A 0 K 0 ), (b) ( Γ ) ( = A Γ i m Γ A 0 K 0 Γ A Γ i m ) ] 0 K 0 Γ S, (c) K ( + )/, where m = (m S +m )/ = MeV/c is the K 0 mean mass, m = m m S = h/s 0], Γ = (Γ S +Γ )/, Γ = Γ S Γ. The state of a single particle system under the evolution given by the master equation with operators () is in general a mixed state of a single particle vacuum (see 5]). After the projection of this state on the single-particle sector one gets the state obtained by the non-trace-preserving equations, considered usually in the literature (see e.g. 3]). A. Master equation in second quantization formalism Now, let us go to the description of a multiparticle system. We achieve this by the replacement of intertwining operators with the appropriate combinations of annihilation creation operators, as was done in 6]. In the base () this Hamiltonian these indbladians take the form H = m ( a a+b b ) pq ( m+ i A A Γ ) a b qp ( m i A A Γ ) b a (3a) ( = Γ S A Γ + m p a+ q Γ +A ( Γ ) Γ i m +A a q Γ = p +A K = Γ( a a+b b ) pq A ) b A A ( Γ ia m ) a b qp A, (3b) ( Γ A Γ i m Γ )b, (3c) ( Γ+iA m ) b a et us consider a single particle state described by a density matrix of the form (3d) ρ = p K 0 K 0 +p K 0 K 0 +w K 0 K 0 +w K 0 K 0 +( p p ) 0 0 (4) 4

5 with 0 p, 0 p, 0 p +p, w p p. The evolution of this state in the Schrödinger picture is governed by the master equation t ρ(t) = ih,ρ(t)]+{k,ρ(t)}+ i i ρ(t) i. (5) By a straightforward calculation, using relation (8) between the bases { K 0 S, K0 } { K 0, K 0 }, one can check that the choices of Eqs. () (3) for operators appearing in Eq. (5) lead to the same evolution equation for ρ(t). Moreover, the evolution of the two particle system obtained from Eqs. (3) (5) lead to the evolution of the two particle state obtained by taking a symmetrized (the particles are indistinguishable) tensor product of single particle evolution (see 7]). Now, consider the Kossakowski-indblad master equation for the evolution of observable Ω(t) in the Heisenberg picture t Ω(t) = Ω(t)], (6) where Ω(t)] = ih,ω(t)]+ i { } i,ω(t)] i + i Ω(t), i], (7) we choose H i (i =,) in the form (3). We assume that Ω(t) can be written as a bilinear form in annihilation creation operators Ω(t) = ω aa (t)a a+ω ab (t)a b+ω ba (t)b a+ω bb (t)b b, (8) with the condition ω ba (t) = ωab (t) to guarantee that Ω(t) is Hermitian. B. General solution Here we want to find the evolution of any observable Ω(t) which is bilinear in annihilation creation operators. To achieve this aim, we begin with the following observation: a a] = Γa a pq A ( Γ i m) a b qp A ( Γ+i m) b a, a b] = Γa b qp +A ( Γ i m) a a qp A ( Γ+i m) b b, b a] = Γb a pq +A ( Γ+i m) a a pq A ( Γ i m) b b, (9a) (9b) (9c) b b] = Γb b pq +A ( Γ+i m) a b qp +A ( Γ i m) b a. (9d) We see that if Ω(t) is of the form (8) at a given moment of time t 0, it must preserve this form for all the time t t 0. If we take into account the linearity of Ω(t) linear independence of operators a a, 5

6 a b, b a, b b, we get the following system of first-order differential equations for ω: ω aa (t) = Γω aa (t) qp +A ( Γ i m) ω ab (t) pq +A ( Γ+i m) ω ba (t), (0a) ω ab (t) = pq A ( Γ i m) ω aa (t) Γω ab (t) pq +A ( Γ+i m) ω bb (t), (0b) ω ba (t) = qp A ( Γ+i m) ω aa (t) Γω ba (t) pq A ( Γ i m) ω bb (t), (0c) ω bb (t) = qp A ( Γ+i m) ω ab (t) pq A ( Γ i m) ω ba (t) Γω bb (t). (0d) Using straightforward methods we can find that the general solution of (0) is of the form ω aa (t) = e Γt { ( Γt) +cos( mt) ] ω aa (0) q ( p Γt) isin( mt) ] ω ab (0) A p ( +A q Γt) +isin( mt) ] ω ba (0) + A ( +A Γt) cos( mt) ] } ω bb (0), (a) ω ab (t) = e Γt { p ( q Γt) isin( mt) ] ω aa (0) ( p + ( Γt) +cos( mt) ] ω ab (0)+ A +A A p +A q ) ( q Γt) cos( mt) ] ω ba (0) ( Γt) +isin( mt) ] } ω bb (0), (b) { ω ba (t) = e Γt +A q ( A p Γt) +isin( mt) ] ω aa (0) + +A ( ) q ( A p Γt) cos( mt) ] ω ab (0)+ ( Γt) +cos( mt) ] ω ba (0) q ( p Γt) isin( mt) ] } ω bb (0), (c) { ω bb (t) = e Γt +A ( A Γt) cos( mt) ] ω aa (0) +A q ( A p Γt) +isin( mt) ] ω ab (0) p ( q Γt) isin( mt) ] ω ba (0) + ( Γt) +cos( mt) ] ω bb (0) }. (d) 6

7 The Hermicity of Ω(t) is preserved (i.e. the condition ω ab (0) = ω ba (0) implies ω ab(t) = ω ba (t) for all the time t 0), since from Eqs. (9) (0) we have ( ) p = +A q q A p. () Now let us observe that the different choices of initial conditions for ω s allow us to get the time evolution of different physically interesting observables of the form (8). IV. EVOUTION OF THE OBSERVABES AND THEIR AVERAGES Now we are prepared to find the evolution of some physically interesting observables, as well as the time dependence of their mean values. We begin with the total number of particles strangeness observables, next we analyze the numbers of each neutral kaon flavor. A. Total number of particles strangeness For the total number of particles observable N, we have N(0) = a a+b b, so ω aa (0) = ω bb (0) = ω ab (0) = ω ba (0) = 0, finally we get { ( N(t) = e Γt +A Γt) +A cos( mt) ] a a p ( +A q Γt) ia sin( mt) ] a b q ( A p Γt) +ia sin( mt) ] b a + ( A Γt) A cos( mt) ] } b b. (3) The expectation value of N(t) in the state n, n is N(t) = e Γt { ( A Γt) A cos( mt)] (n+ n) ( A Γt) cos( mt) ] (n n) } (4) is depicted in Fig.. If we take A = 0 in Eq. (4), we obtain the case with preserved CP symmetry, for which we recover the usual result for mean number of particles ( N(t) CP = e Γ S t +e ) Γ t (n+ n). (5) The leading part of the difference between CP-violated CP-preserved values is N(t) N(t) CP = A ( e Γ S t +e ) Γ t e Γt cos( mt) ] (n n)+o(a ), (6) is shown in Fig., too. 7

8 FIG. (Color online). Mean number of particles vs time for n + n =,...,5 (bottom-up) the difference between N(t) for CP-violated CP-preserved case vs time for n+ n = 4 n n = 0,...,4 (top-down). Forstrangeness S we have S(0) = a a b b, so ω aa (0) = ω bb (0) =, ω ab (0) = ω ba (0) = 0, S(t) = e Γt { +A A ( Γt) +cos( mt) ] a a p A ( +A q Γt) isin( mt) ] a b A q p A ( Γt) +isin( mt) ] b a + A A ( Γt) cos( mt) ] b b }. (7) The expectation value of S(t) in the state n, n is therefore S(t) = e Γt { cos( mt) A A ( Γt)] (n n) ( +A Γt) cos( mt) ] (n+ n) } (8) is shown in Fig.. As previously, when A = 0 we get the CP-preserved value 8

9 FIG. (Color online). Mean strangeness vs time for n+ n = 4 n n = 0,...,4 (bottom-up) with the strangeness oscillation phenomenon explicitly visible when n n 0 the difference between S(t) for CP-violated CP-preserved cases vs time for n+ n =,...,5 n n = 0 (bottom-up). S(t) CP = e Γt cos( mt)(n n), (9) so the leading term of the difference between CP-violated CP-preserved values is S(t) S(t) CP = A ( e Γ S t +e Γ t ) e Γt cos( mt) ] (n+ n)+o(a ) (30) is shown also in Fig.. 9

10 B. Number of K 0 K 0 For the numbers of K 0 K 0 we have N K 0(0) = a a N K 0(0) = b b, so N K 0(t) = e Γt { ( Γt) +cos( mt) ] a a ( Γt) isin( mt) ] a b p q +A q ( A p Γt) +isin( mt) ] b a + +A ( A Γt) cos( mt) ] } b b (3a) N K 0(t) = e Γt Their mean values are N K 0(t) = e Γt N K 0(t) = e Γt { A +A ( Γt) cos( mt) ] a a A p +A q q p ( Γt) +isin( mt) ] a b ( Γt) isin( mt) ] b a + ( Γt) +cos( mt) ] } b b. (3b) { ( Γt) +cos( mt) ] n+ +A ( A Γt) cos( mt) ] n } (3a) { A ( +A Γt) cos( mt) ] n+ ( Γt) +cos( mt) ] n }, (3b) respectively. We show these time evolution in the Figs Notice that after a suitable period of time (approx. 8ns) different initial states can give the same mean number of K 0 or K 0. If we put A = 0, we get again CP-preserved values N K 0(t) CP = e Γ St +e Γ t n+ n +e Γt n n cos( mt) (33a) N K 0(t) CP = e ΓSt +e Γ tn+ n e Γt n n cos( mt). (33b) The leading term of the differences for the CP-violated CP-preserved vaues are N K 0(t) N K 0(t) CP = A ( e Γ S t +e ) Γ t e Γt cos( mt) ] n+o(a ), (34a) N K 0(t) N K 0(t) CP = A ( e Γ S t +e ) Γ t e Γt cos( mt) ] n+o(a ). (34b) 0

11 FIG. 3 (Color online). Mean number of K 0 vs time for n =,,3 (bundles bottom-up) n = 0,, (bottom-up inside bundles). C. Number of K 0 S K0 et us return to the short long living states of neutral kaon, K 0 S K0, respectively. If we define another two pairs of annihilation creation operators, such that K 0 S = c S 0, (35a) K 0 = c 0, (35b) we can find, with use of Eq. (8), that c S = p a+q b, c = p a q b. (36a) (36b) The operators c S c fulfill almost the usual canonical commutation relations with the following exception: c S,c ] = c,c S ] = A. (37) This reflects the fact that K 0 S K0 are not orthogonal, see Eq. ().

12 FIG. 4 (Color online). Mean number of K 0 vs time for n =,,3 (bottom-up inside bundles) n = 0,, (bundles bottom-up). Now we define multiparticle states for K 0 S K0 as #K 0 S = n n S = (c S )n n! 0, #K 0 = n n = (c )n n! 0. (38a) (38b) These states can be expressed in terms of the sates n,k as follows n S = n = n (n ) p n k q k n k,k, k n (n ) ( ) k p n k q k n k,k k k=0 k=0 (39a) (39b)

13 For the number of K 0 S, N K 0(0) = c S S c S = +A a a+ A p q a b+ +A q p b a+ A b b, so N K 0 S (t) = e Γt { ] e Γt +A +A e Γt +A cos( mt) a a + ] p e Γt A +A q e Γt +ia sin( mt) a b + ] q e Γt A A p e Γt ia sin( mt) b a + A e Γt +A e Γt A cos( mt) ] } b b. (40) Similarly, for the number of K 0, N K 0(0) = c c = +A a a A p q a b +A q p b a+ A b b, N K 0 (t) = e Γt { ] e Γt +A +A e Γt +A cos( mt) a a ] p e Γt A +A q e Γt ia sin( mt) a b ] q e Γt A A p e Γt +ia sin( mt) b a + A e Γt +A e Γt A cos( mt) ] } b b. (4) Now we are prepared to find the expectation values of N K 0 S (t) N K 0 (t) in the states n S n, respectively. First, observe that Thus a b n S = q p a a n S, a b n = q p a a n, (4a) b a n S = p q b b n S, b a n = p q b b n. (4b) ) ] N K 0 S (t) n S = (e e Γt A Γt A )N(0) A ei mt (e Γt e i mt S(0) n S (43a) ) ] N K 0 (t) n = (e e Γt A Γt A )N(0) A e i mt (e Γt e i mt S(0) n. (43b) Of course n S n are eigenvectors of N(0): N(0) n S = n n S, N(0) n = n n, (44a) (44b) 3

14 but they are not eigenvectors of S(0) since n (n ) S(0) n S = n n S k p n k q k n k,k, k k=0 n (n ) S(0) n = n n ( ) k k p n k q k n k,k. k k=0 (45a) (45b) Notice that in view of Eq. (45) the states n S n are not eigenvectors of N K 0 S N K 0 respectively, as one can naïvely expect. However, the expectation values of S(0) in the states n S n are well established can be calculated as because n S S(0) n S = A n, (46a) n S(0) n = A n, (46b) n ( ) n k p (n k) q k = n q ( p + q ) n = n q. (47) k k=0 In conclusion, we get that the mean number of K 0 S K0 in the states n S n, respectively, evolves in time according to the Geiger-Nutall law n S N K 0 S (t) n S = ne Γ St, (48a) n N K 0 (t) n = ne Γ t. (48b) The states n S n are not orthogonal, thus we expect nonzero expectation values of N K 0 S (t) N K 0 in the states n n S, respectively. Indeed, N K 0 S (t) n = e Γt A (e i mt A e ) ] A )S(0) A Γt (e i mt e Γt N(0) n (49a) thus, N K 0 (t) n S = e Γt A A (e i mt A e ) ] )S(0) A Γt (e i mt e Γt N(0) n S, (49b) n N K 0 S (t) n = na e Γ t, (50a) n S N K 0 (t) n S = na e Γ St. (50b) So, finally we see that we can detect a fractionof order A of the other flavor in multiparticle short- or long-lived neutral kaon states. V. CONCUSIONS We have shown that it is possible to formulate the consistent probability-preserving description of the CP-symmetry-violating evolution of a system containing any number of 4

15 decaying particles. This has been done within the framework of quantum mechanics of open systems based on the approach developed in 5, 6]. To achieve this aim we have considered master equations built up from creation annihilation operators which generate dynamical semigroups that can describe the exponential decay flavour oscillations for a system of many particles. It should be noted that this dynamical semigroup was used for the description of the entire system of unstable particles, not only for decoherence effects, as it was done in, e.g., 9, 4]. We have used the fact that the decay of the particle can be regarded as a Markov process (the probability of decay of a particle is constant, so it does not depend on its history) we have solved explicitly the Kossakowski-indblad master equation for a system of particles with violated CP symmetry. To avoid working with infinite numbers of density matrix elements we performed calculations in the Heisenberg picture of quantum mechanics. This allowed us to deal with a low-dimensional system of linear differential equations for evolution of any observable bilinear in creation annihilation operators. The choice of a concrete observable is then reduced to the proper choice of initial conditions for the system of differential equations. To show explicitly the effectiveness of the introduced approach we have calculated the evolution of the observables which are most interesting from the physical point of view. We have also found the evolution of mean values of these observables as well as their lowest order difference with the CP-preserved values. It seems to us that the presented analysis of time evolution of neutral kaons could be applied to the description of EPR states. Of course one must construct such a state by labeling annihilation creation operators either by a discrete index, as it is done in quantum optics, or by a continuous parameter, as it is done in quantum field theory (beware that in the latter case annihilation creation operators are actually operator valued distributions). Moreover, we think that our approach would be interesting helpful also in analysis of mechanical systems electric circuits. This follows from the fact that there is a strong analogy with nonreversible classical mechanics quantum systems with CP violation ]. Finally, we would like to point out that all the presented results are also valid for neutral B mesons after appropriate change of notation values of physical quantities. This follows from the fact that B mesons evolve according to the same scheme as kaons. ACKNOWEDGMENTS This work is supported by the Polish National Science Centre under Contract No. 04/5/B/ST/007 by the University of odz. ] J. S. Bell, Physics (NY), 95 (964); J. F. Clauser, M. A. Horne, A. Shimony, R. A. Holt, Phys. Rev. ett. 3, 880 (969); J. F. Clauser A. Shimony, Rep. Prog. Phys. 4, 88 (978). ] F. Uchiyama, Phys. ett. A 3, 95 (997); A. Bramon M. Nowakowski, Phys. Rev. ett. 83, (999); R. Foadi F. Selleri, Phys. Rev. A 6, 006 (999); R. H. Dalitz G. Garbarino, Nucl. Phys. B 606, 483 (00); R. Bertlmann, W. Grimus, B. C. Hiesmayr, Phys. ett. A 89, (00); R. A. Bertlmann B. C. Hiesmayr, Phys. Rev. 5

16 A 63, 06 (00); A. Bramon G. Garbarino, Phys. Rev. ett. 89, 6040 (00); M. Genovese, Phys. Rev. A 69, 003 (004). 3] R. A. Bertlmann, K. Durstberger, B. C. Hiesmayr, Phys. Rev. A 68, 0 (003). 4] A. Go, J. Mod. Opt. 5, 99 (004); A. Go et al., Phys. Rev. ett. 99, 380 (007); A. Bramon, R. Escribano, G. Garbarino, J. Mod. Opt. 5, 68 (005). 5] A. Einstein, B. Podolsky, N. Rosen, Phys. Rev. 47, 777 (935); D. Bohm, Quantum Theory (Prentice-Hall, Englewood Cliffs, 95). 6] V. Weisskopf E. Wigner, Z. Phys. 63, 54 (930); 65, 8 (930). 7] T. Durt, Int. J. Mod. Phys. B 7, (03). 8] J. Ellis, J.. opez, N. E. Mavromatos, D. V. Nanopoulos, Phys. Rev. D 53, 3846 (996); J. Bernabéu, N. E. Mavromatos, P. Villanueva-Pérez, Phys. ett. B 74, 69 (03). 9] F. Benatti R. Floreanini, Phys. ett. B 389, 00 (996); Nucl. Phys. B 488, 335 (997); Phys. ett. B 40, 337 (997); Nucl. Phys. B 5, 550 (998); Phys. Rev. D 57, R33 (998); F. Benatti, in Irreversibility Causality Semigroups Rigged Hilbert Spaces, ecture Notes in Physics, Vol. 504, edited by A. Bohm, H.-D. Doebner, P. Kielanowski (Springer, Berlin Heidelberg, 998) pp. 4 47; F. Benatti, R. Floreanini, R. Romano, Nucl. Phys. B 60, 54 (00). 0] R. Alicki M. Fannes, Quantum Dynamical Systems (Oxford University Press, 00); H.-P. Breuer F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 00). ] R. Alicki, Rep. Math. Phys., (977); 4, 7 (978); in Quantum Dynamical Semigroups Applications, ecture Notes in Physics, Vol. 77, edited by R. Alicki K. endi (Springer, Berlin Heidelberg, 007) pp. 46. ] A. Weron, A. K. Rajagopal, K. Weron, Phys. Rev. A 3, 736 (985). 3] Y.-X. iu, Ş. K. Özdemir, A. Miranowicz, N. Imoto, Phys. Rev. A 70, (004). 4] R. A. Bertlmann, W. Grimus, B. C. Hiesmayr, Phys. Rev. A 73, 0540 (006). 5] P. Caban, J. Rembieliński, K. A. Smoliński, Z. Walczak, Phys. Rev. A 7, 0306 (005). 6] K. A. Smoliński, Open Syst. Inf. Dyn., (04). 7] P. Caban et al., Phys. ett. A 357, 6 (006). 8] P. Caban, J. Rembieliński, K. A. Smoliński, Z. Walczak, Phys. ett. A 363, 389 (007). 9] A. Kossakowski, Rep. Math. Phys. 3, 47 (97); V. Gorini, A. Kossakowski, E. C. G. Sudarshan, J. Math. Phys. 7, 8 (976); G. indblad, Commun. Math. Phys. 48, 9 (976). 0] K. A. Olive et al., Chin. Phys. C 38, (04). ] T. T. Wu C. N. Yang, Phys. Rev. ett. 3, 380 (964). ] J.. Rosner, Am. J. Phys. 64, 98 (996); D. Cocolicchio. Telesca, Found. Phys. ett., 3 (998); V. A. Kostelecký Á. Roberts, Phys. Rev. D 63, (00); M. Caruso, H. Fanchiotti, C. A. Garcia Canal, Ann. Phys. (N.Y.) 36, 77 (0). 6