Preprint
Article

Linear Entropy in the Atom-Field Interaction at Finite Temperature

Altmetrics

Downloads

89

Views

29

Comments

0

This version is not peer-reviewed

Submitted:

06 October 2023

Posted:

07 October 2023

You are already at the latest version

Alerts
Abstract
We study the atom-field interaction at finite temperature and in the dispersive regime. We show that the master equation for this sytem may be solved with the use of superoperator techniques. We calculate the linear entropy in case the field is initially in a coherent state and the atom in a superposition of its ground and excited states.
Keywords: 
Subject: Physical Sciences  -   Atomic and Molecular Physics

1. Introduction

One of the mots important quantities to measure entanglement between two subsystems [1,2] is the von Neumann entropy [3]. It defines, among other things and very clearly, when two subsystems are entangled. In order to calculate the von Neumann entropy for one of the subsystems usually it is used the Araki-Lieb inequality [4,5], that allows, when initially the two sub-systems are found in pure states, to relate the two sub-systems entropies after they interact. To achieve this it is necessary to perform a partial trace on the system, such that we remain with only one subsystem described by its density matrix which is, in the most general case, a so-called "statistical mixture" state. It is more difficult the calculation of the entropy when one studies statistical mixtures [5,6] because it is a (logarithmic) function of the density matrix and the calculation may become extremely complicated [7,8].
One of the interactions where entanglement is usually studied is the resonant atom field interaction [9] where it is well known that the atomic inversion shows collapses and revivals of the Rabi oscillations [10,11,12]. Such revivals may be considered as a first measurement of the field as it gives light on the nature of the photon distribution of the cavity field: For instance, in the case that an initial squeezed state is considered, the atomic inversion shows so-called ringing revivals [10,11]. Entropy, together with the atomic inversion, may be used to obtain information about the non-classicality of a given state [6]. For example, if at an specific time, the atomic inversion is in the collapse region and the entropy of a field initially prepared in an squeezed state is close to zero it is well-known that a superposition of squeezed states may be generated [11].
However, if a superposition of coherent states, the so-called Schrödinger cat states, is considered as initial state, the revival of oscillations occurs sooner [12] than that of a single coherent state [13]. But, if a mixture of coherent states is initially considered, the revival of Rabi oscillations appear as in the case of a single coherent state. This indicates that the atomic inversion, together with the degree of mixedness of states may give information about the initial state used in a given interaction [5,14,15]. Such mixedness may be measured via the von Neumann entropy.
The von Neumann entropy is given by the expression
S = T r { ρ ln ρ } ,
where ρ , that defines the state of the system, is the density matrix and T r means the trace over the system’s degrees of freedom. Calculation of this quantity, however, may become a difficult task as mentioned before. A more convenient quantity to calculate because of its simplicity that however contains the same information as the von Neumann entropy is the linear entropy [16], defined as
S L = 1 T r { ρ 2 } .
The eigenbasis of the density matrix may be used to show that
T r { ρ 2 } = n ρ n 2 n ρ n = 1 ,
where the equality holds for pure states only. Therefore S L uniquely discriminates between pure and mixed states, just as the von Neumann entropy does.
In the atom-field interaction there are effects such as decoherence that play a strong role in the purity of the states of quantum systems [17]. Effects of environment strongly damage entanglement and the nonclassicality of states of light [18]. Entanglement may be even damaged by an environment where gravitational fluctuations are considered [19,20]. Moreover, studies to endure quantum coherence by using photon indistinguishability in noisy quantum networks have been performed [21]. When an environment is taken into account, however, complete information about an initial quantum state may be still recovered by using reconstruction techniques [22]. Such reconstruction technique are realized using quasiprobability distribution functions [23] that may be even generated in classical physics [24]. The damage becomes worst in the case that the environment is at finite temperature as the bath injects photons into the cavity destroying faster the possible coherences of the quantized field.
In this contribution we study this case and show how to calculate the linear entropy in the dispersive interaction between a two-level atom and a quantized field when we consider losses at temperature T > 0 . In the next Section we present the master equation that describes this system and show how to solve it, then in Section 3 we calculate the field linear entropy and do some numerical analysis in Section 4. Section 5 is left for conclusions.

2. Master Equation for a Real Cavity

The master equation for a two- level atom interacting with a quantized electromagnetic field in the dispersive regime, i.e., when the atom and the field stop exchhanging energy because they are sufficiently detuned, at finite temperature is given by [25,26,27] (we set = 1 )
d ρ ^ d t = i H ^ e f f , ρ ^ + k 1 2 a ^ ρ ^ a ^ a ^ a ^ ρ ^ ρ ^ a ^ a ^ + k 2 2 a ^ ρ ^ a ^ a ^ a ^ ρ ^ ρ ^ a ^ a ^ ,
where H ^ e f f = χ n ^ σ ^ z [25,26] is the so-called dispersive Hamiltonian, and k 1 = C 2 n ¯ , k 2 = C 2 n ¯ + 1 and n ¯ = 1 e v k k β T 1 . Here C is the decay constant and n ¯ is the average number of thermal photons. Also, the operators a ^ and a ^ are the annihilation and creation operators, respectively. By defining k 12 = k 1 + k 2 , the above master equation (4) may be rewritten as
d ρ ^ d t = i χ ρ ^ n ^ σ ^ z i χ n ^ σ ^ z ρ ^ + 2 k 1 a ^ ρ ^ a ^ + 2 k 2 a ^ ρ ^ a ^ k 12 ρ ^ n ^ k 12 n ^ ρ ^ 2 k 1 ρ ^ .
We may define the superoperators (see for instance [28])
J ^ 1 ρ ^ = 2 k 1 a ^ ρ ^ a ^ , J ^ 2 ρ ^ = 2 k 2 a ^ ρ ^ a ^ , R ^ ρ ^ = i χ ρ ^ n ^ σ ^ z i χ σ ^ z n ^ ρ ^ , S ^ ρ ^ = k 12 ρ ^ n ^ k 12 n ^ ρ ^ ,
to rewrite the master equation in a more compact form, namely
d ρ ^ d t = 2 k 1 ρ ^ + J ^ 1 ρ ^ + J ^ 2 ρ ^ + R ^ ρ ^ + S ^ ρ ^ .
This allows us to write a formal solution
ρ ^ t = e 2 k 1 t e J ^ 1 + J ^ 2 + R ^ + S ^ t ρ ^ 0 ,
that may be factorized in the form
ρ ^ t = e f 0 t e f 1 t J ^ 1 e f 2 t R ^ e f 3 t S ^ e f 4 t J ^ 2 ρ ^ 0 .
We may show that indeed (9) is a solution by deriving ρ ^ ( t ) with respect to time to obtain
d ρ ^ d t = d f 0 t d t + d f 1 t d t J ^ 1 ρ ^ t + d f 2 t d t e f 1 t J ^ 1 R ^ e f 1 t J ^ 1 ρ ^ t + d f 3 t d t e f 1 t J ^ 1 e f 2 t R ^ S ^ × e f 2 t R ^ e f 1 t J ^ 1 ρ ^ t + d f 4 t d t e f 1 t J ^ 1 e f 2 t R ^ e f 3 t S ^ J ^ 2 × e f 3 t S ^ e f 2 t R ^ e f 1 t J ^ 1 ρ ^ t .
We define the atomic superoperator L ^ ρ ^ = i χ σ ^ z ρ ^ i χ ρ ^ σ ^ z and, by using the relevant commutators given in Appendix A, and using that e w A ^ B ^ e w A ^ = B ^ + w A ^ , B ^ + w 2 2 ! A ^ , A ^ , B ^ + . . . we obtain the terms present in equation, (10) i.e.,
e f 1 t J ^ 1 R ^ e f 1 t J ^ 1 ρ ^ t = R ^ + f 1 t L ^ J ^ 1 ρ ^ t ,
e f 1 t J ^ 1 e f 2 t R ^ S ^ e f 2 t R ^ e f 1 t J ^ 1 = S ^ ρ ^ t + 2 k 12 f 1 t J ^ 1 ρ ^ t
and
e f 1 J ^ 1 e f 2 R ^ e f 3 S ^ J ^ 2 e f 3 S ^ e f 2 R ^ e f 1 J ^ 1 = [ e 2 k 12 f 3 + f 2 L ^ J ^ 2 + 4 k 1 k 2 k 12 f 1 e 2 k 12 f 3 + f 2 L ^ S ^ 4 k 1 k 2 f 1 e 2 k 12 f 3 + f 2 L ^ + 4 k 1 k 2 f 1 2 e 2 k 12 f 3 + f 2 L ^ J ^ 1 ] ρ ^ t .
By substituting these expressions in (10) and compare it with equation (7) we obtain the set of first order differential equations
d f 0 d t 4 k 1 k 2 f 1 e 2 k 12 f 3 + f 2 L ^ d f 4 d t = 2 k 1 , d f 1 d t + d f 2 d t f 1 L ^ + 2 k 12 d f 3 d t f 1 + 4 k 1 k 2 f 1 2 e 2 k 12 f 3 + f 2 L ^ d f 4 d t = 1 , d f 2 d t = 1 , d f 3 d t + 4 k 1 k 2 k 12 f 1 e 2 k 12 f 3 + f 2 L ^ d f 4 d t = 1 , e 2 k 12 f 3 + f 2 L ^ d f 4 d t = 1 .
The above system of differential equations depends on the superoperator L ^ , but, because there are no other superoperators in it, a simple solution may be found for the different functions
f 0 t , L ^ = L ^ + 2 k 12 4 k 1 2 t , + ln cos θ ( L ^ ) cos β ( L ^ ) 2 t + θ ( L ^ ) , f 1 t , L ^ = β ( L ^ ) 8 k 1 k 2 tan β ( L ^ ) 2 t + θ ( L ^ ) + L ^ + 2 k 12 8 k 1 k 2 , f 2 t = t , f 3 t , L ^ = L ^ 2 k 12 t + 1 k 12 ln cos β ( L ^ ) 2 t + θ ( L ^ ) cos θ ( L ^ ) f 4 t , L ^ = 2 cos 2 θ ( L ^ ) β ( L ^ ) tan β ( L ^ ) 2 t + θ ( L ^ ) + 2 L ^ + 2 k 12 cos 2 θ ( L ^ ) β ( L ^ ) 2 ,
with
β ( L ^ ) = 16 k 1 k 2 L ^ + 2 k 12 2
and
θ ( L ^ ) = arctan L ^ + 2 k 12 16 k 1 k 2 L ^ + 2 k 12 2
.

3. Linear entropy for the quantized field

Once we have calculated the total density matrix we may calculate the reduced field operator.
ρ ^ f t = T r A { ρ ^ t }
where ρ ^ 0 = ρ ^ f 0 ρ ^ A 0 is the initial density matrix for the total, atom-field, system. We consider an arbitrary initial field and the atom in a superposition of its excited and ground state, ρ ^ A 0 = ψ A 0 ψ A 0 , with ψ A 0 = 1 2 e + g therefore ρ ^ A 0 = 1 2 I ^ + e g + g e .
After some algebra we obtain
ρ ^ f t = r = 0 k = 0 d 1 t , r , k a ^ k e k 12 f 3 t , 0 n ^ × [ e i χ n ^ t a ^ r ρ ^ f 0 a ^ r e i χ n ^ t + e i χ n ^ t a ^ r ρ ^ f 0 a ^ r e i χ n ^ t ] × e k 12 f 3 t , 0 n ^ a ^ k ,
with
d 1 t , r , k = e f 0 t , 0 2 k 2 f 4 t , 0 r 2 r ! 2 k 1 f 1 t , 0 k k ! .
If we consider that ρ ^ f 0 = α α , with the coherent state given by [13] α = e α 2 2 n = 0 α n n ! n , we obtain
ρ ^ f t = r = 0 k = 0 n = 0 m = 0 2 d 1 t , r , k × α 2 r α n α m n ! m ! e α 2 k 12 f 3 t , 0 n + m × a ^ k n m a ^ k cos χ n m t
we need to calculate a ^ k n which is easy to see that is given by
a ^ k n = ( n + k ) ! n ! n + k
and substitituing d 1 t , r , k , us to write
ρ ^ f t = k = 0 n = 0 m = 0 e k 12 f 3 t , 0 n + m 2 k 1 f 1 t , 0 k × ( n + k ) ! ( m + k ) ! cos χ n m t k ! n ! m ! × α n + m e f 0 t , 0 + 2 k 2 f 4 t , 0 1 α 2 n + k m + k
setting j = n + k y l = m + k we find the reduced field operator
ρ ^ f t = e f 0 t , 0 + 2 k 2 f 4 t , 0 1 α 2 r = 0 g = r h = r × g ! h ! α e k 12 f 3 t , 0 g + h 2 r r ! g r ! h r ! × 2 k 1 f 1 t , 0 r cos χ g h t g h
Once we have calculated the reduced field operator we may calculate the field entropy S ^ f
S ^ f t = 1 T r ρ ^ f 2 t
Substitituing the found solution for the reduced field operator (21) in the above equation, after some algebra we obtain the field entropy
S ^ f t = 1 e f 0 t , 0 + 2 k 2 f 4 t , 0 1 α 2 2 × n = 0 k = 0 r = 0 j = 0 n ! j ! 2 k 1 f 1 t , 0 k + r n r ! j r ! k ! r ! × α e k 12 f 3 t , 0 2 j + n k r cos 2 χ j n t j k ! n k !

4. Numerical results

We plot the linear field entropy (23) with the help of numerical analysis. In the Figure 1, Figure 2 and Figure 3 we set the parameter k 1 = 0.01 and k 1 = 0.05 that correspond to a decay constant C = 0.08 and an average number of photons n ¯ = 0.25 . The linear entropy suffers oscillations that stabilize around a value of S f 0.33 . It may be seen that for smaller values of the coherent intensity, namely, α = 1 in Figure 1 such stabilization is faster. This is because of the fact that there are not many coherent photons, initially, inside the cavity and they are replaced by the thermal field as the cavity is at a finite temperature. Eventually, the cavity is filled by a thermal field
ρ T h = n = 0 n ¯ n ( n ¯ + 1 ) n + 1 | n n | ,
that has a linear entropy of [see equation (3)]
S T h = 1 n = 0 n ¯ n ( n ¯ + 1 ) n + 1 2 ,
that, for the parameters chosen in Figure 1, Figure 2 and Figure 3, S T h 0.33 .
The same occurs with the linear entropies plotted in Figure 4, Figure 5 and Figure 6. In these figures, we set k 1 = 0.005 and k 2 = 0.01 which deliver a decay constant C = 0.01 , smaller than the previous examples, which means that the emptying of the coherent filed and the filling of the thermal field are slower. On the other hand, the average number of photons for this choice of parameters is about n ¯ = 1 that produces a linear entropy for the thermal field of about S T h 0.66 , precisely the stabilization line shown in Figure 4, Figure 5 and Figure 6.

5. Conclusions

We have studied the atom-field interaction at finite temperature when the field is initially in a coherent state and in the dispersive regime. By using superoperator techniques, we manage to solve the master equation for this system and have calculated the linear entropy. We have shown that the linear entropy for the dispersive interaction subject to decay at T > 0 approaches the limit given by the thermal distribution of the bath.

Appendix A

Here we list the commutators between the different superoperators needed to solve equation (7),
S ^ , J ^ 1 ρ ^ = 2 k 12 J ^ 1 ρ ^ , R ^ , J ^ 2 ρ ^ = L ^ J ^ 2 ρ ^ , R ^ , J ^ 1 ρ ^ = L ^ J ^ 1 ρ ^ , S ^ , J ^ 2 ρ ^ = 2 k 12 J ^ 2 ρ ^ , J ^ 1 , J ^ 2 ρ ^ = 4 k 1 k 2 k 12 S ^ ρ ^ 4 k 1 k 2 ρ ^ , R ^ , L ^ ρ ^ = 0 , R ^ , S ^ ρ ^ = 0 , S ^ , L ^ ρ ^ = 0 , J ^ 1 , L ^ ρ ^ = 0 , J ^ 2 , L ^ ρ ^ = 0 .
It is worth to mention that the superoperators involved in the master equation, namely, J ^ 1 , J ^ 2 and S ^ commute in such a way that only give the same superoperators, allowing an algebraic solution.

References

  1. C. C. Gerry and P. L. Knight, Introductory Quantum Optics (Cambridge University Press, 2005).
  2. Abdel-Kahlek, K. Berrada, and H. Eleuch, Effect of the time-dependent coupling on a superconducting qubit-field system under decoherence: Entanglement and wehrl entropy, Ann. of Phys 361, 247–258. (2015). [CrossRef]
  3. J. von Neumann, Mathematical Foundations of Quantum Mechanics (Princeton University Press, 1955).
  4. H. Araki and E. H. Lieb, Entropy inequalities, Commun. Math. Phys. 18, 160–170. (1970).
  5. J. D. Phoenix and P. L. Knight, Fluctuations and entropy in models of quantum optical resonance, Annals of Physics 186, 381–407. (1988). [CrossRef]
  6. A. Zúñiga-Segundo, R. Juárez-Amaro, O. Aguilar-Loreto, and H. Moya-Cessa, Field’s entropy in the atom-field interaction: Statistical mixture of coherent states, Annals of Physics 379, 150–158. (2017). [CrossRef]
  7. B. M. Rodríguez-Lara and H. M. Moya-Cessa, Por qué y cómo encontramos funciones de matrices: entropia en mecánica cuántica, Revista Mexicana de Física E 51, 87–98 (2005).
  8. A. Motazedifard, S. A. Madani, and N. S. Vayaghan, Measurement of entropy and quantum coherence properties of two type-i entangled photonic qubits, Opt. Quant. Electr. 53, 378–401. (2021). [CrossRef]
  9. E. T. Jaynes and F. W. Cummings, Comparison of quantum and semiclassical radiation theories with application to the beam maser, Proc. IEEE. 51, 89–109. (1963).
  10. M. Satyanarayana, P. Rice, V. R., and H. Carmichael, Ringing revivals in the interaction of a two-level atom with squeezed light, J. Opt. Soc. Am. B 6, 228–237. (1989).
  11. H. Moya-Cessa and A. Vidiella-Barranco, Interaction of squeezed states of light with two-level atoms, J. of Mod. Optics 39, 2481–2499. (1992). [CrossRef]
  12. A. Vidiella-Barranco, H. Moya-Cessa, and V. Buzek, Interaction of superpositions of coherent states of light with two-level atoms, J. of Mod. Optics 39, 1441–1459. (1992). [CrossRef]
  13. R. J. Glauber, Coherent and incoherent states of the radiation field, Phys. Rev. 131, 2766 (1963). [CrossRef]
  14. S. M. Barnett, A. Beige, A. Ekert, B. M. Garraway, C. H. Keitel, V. Kendon, M. Lein, G. J. Milburn, H. Moya-Cessa, M. Murao, J. K. Pachos, G. M. Palma, E. Paspalakis, S. J. D. Phoenix, B. Piraux, M. B. Plenio, B. C. Sanders, J. Twanmley, A. Vidiella-Barranco, and M. S. Kim, Journeys from quantum optics to quantum technology, Progress in Quantum Electronics 54, 19–45. (2017). [CrossRef]
  15. J. Gea-Banacloche, Atom- and field-state evolution in the jaynes-cummings model for large initial fields, Phys. Rev. A 44, 5913–5931. (1991). [CrossRef]
  16. R. Juárez-Amaro, L. M. Rosales, J. A. Anaya-Contreras, Zúñiga-Segundo, and H. M. Moya-Cessa, Relation between the entropy and the linear entropy in the ion–laser interaction, Journal of Modern Optics 67, 805–810. (2020). [CrossRef]
  17. S. J. D. Phoenix, Wave-packet evolution in the damped oscillator, Phys. Rev. A 41, 5132–5138. (1990). [CrossRef]
  18. L. M. Arévalo-Aguilar and H. M. Moya-Cessa, Quantum bits and superposition of displaced fock states of the cavity field, Revista Mexicana de Física 48, 423 (2002). [CrossRef]
  19. H. Wu, X. Fan, and L. Chen, How gravitational fluctuations degrade the high dimensional spatial entanglement, Phys. Rev. D 106, 045023. (2022).
  20. H. Wu, X. Fan, and L. Chen, Twisting and entangling gravitons in high-dimensional orbital agular momentum states ia photon-graviton conversion, Phys. Rev. D 107, 125027 (2023). [CrossRef]
  21. A. Perez-Leija, D. Guzmán-Silva, L.-M. R. de J., M. Gräfe, M. Heinrich, H. M. Moya-Cessa, K. Busch, and A. Szameit, Endurance of quantum coherence due to particle indistinguishability in noisy quantum networks, npj Quantum Information 4, 45 (2018). [CrossRef]
  22. N. Yazdanpanah, M. K. Tavassoly, R. Juárez-Amaro, and H. M. Moya-Cessa, Reconstruction of quasiprobability distribution functions of the cavity field considering field and atomic decays, Optics Communications 400, 69–73. (2017). [CrossRef]
  23. H. Moya-Cessa, Entropy operator and associated wigner function, International Journal of Quantum Information 5, 149 (2023). [CrossRef]
  24. J. Moya-Cessa, L. Berriel-Valdos, and H. Moya-Cessa, Optical production of the husimi function of two gaussian functions, Applied Mathematics Infromation Sciences 2, 309–316 (2008).
  25. A. B. Klimov and L. L. Sanchez-Soto, Method of small rotations and effective hamiltonians in nonlinear quantum optics, Physical Review A 61, 063802 (2000). [CrossRef]
  26. S. M. Chumakov, A. B. Klimov, and C. Saavedra, Dispersive atomic evolution in a dissipative-driven cavity, Physical Review A 61, 033814 (2000). [CrossRef]
  27. A. B. Klimov, L. L. Sanchez-Soto, and J. Delgado, Mimicking a kerrlike medium in the dispersive regime of second-harmonic generation, Optics Communications 191, 419 (2001). [CrossRef]
  28. L. M. Arévalo-Aguilar and H. Moya-Cessa, Cavidad con pérdidas: una descripción usando superoperadores, Revista Mexicana de Física 42, 675 (1996).
Figure 1. Plot of the linear entropy, S f ( t ) , as a function of time. We set the parameters χ = 1 , k 1 = 0.01 , k 2 = 0.05 , and α = 1 .
Figure 1. Plot of the linear entropy, S f ( t ) , as a function of time. We set the parameters χ = 1 , k 1 = 0.01 , k 2 = 0.05 , and α = 1 .
Preprints 87136 g001
Figure 2. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 1 except the amplitude of the initial coherent field: α = 3 .
Figure 2. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 1 except the amplitude of the initial coherent field: α = 3 .
Preprints 87136 g002
Figure 3. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 1 except the amplitude of the initial coherent field: α = 5 .
Figure 3. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 1 except the amplitude of the initial coherent field: α = 5 .
Preprints 87136 g003
Figure 4. Plot of the linear entropy, S f ( t ) , as a function of time. We set the parameters χ = 1 , k 1 = 0.005 , k 2 = 0.01 , and α = 1 .
Figure 4. Plot of the linear entropy, S f ( t ) , as a function of time. We set the parameters χ = 1 , k 1 = 0.005 , k 2 = 0.01 , and α = 1 .
Preprints 87136 g004
Figure 5. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 4 except the amplitude of the initial coherent field: α = 3 .
Figure 5. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 4 except the amplitude of the initial coherent field: α = 3 .
Preprints 87136 g005
Figure 6. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 4 except the amplitude of the initial coherent field: α = 5 .
Figure 6. Plot of the linear entropy, S f ( t ) , as a function of time. We set the same parameters as in Figure 4 except the amplitude of the initial coherent field: α = 5 .
Preprints 87136 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated