Next Article in Journal
Effect of Twisting Phases on Linear–Circular Polarization and Spin–Orbital Angular Momentum Conversions in Tightly Focused Vector and Scalar Beams
Next Article in Special Issue
A Tunable Optoelectronic Oscillator with Phase-to-Amplitude Modulation Transformation via an Acetylene Reference Cell
Previous Article in Journal
Implementing Predictive Models in Artificial Intelligence through OCT Biomarkers for Age-Related Macular Degeneration
Previous Article in Special Issue
Flexible, Stretchable, Tunable, and Switchable DFB Laser
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Statistical Mixture of Kaleidoscope States Interacting with a Two-Level Atom: Entropy and Purification

by
Jorge A. Anaya-Contreras
1,†,
Arturo Zúñiga-Segundo
1,† and
Héctor M. Moya-Cessa
2,*,†
1
ESFM Departamento de Física, Edificio 9 Unidad Profesional Adolfo López Mateos, Instituto Politécnico Nacional, México 07738, Mexico
2
Instituto Nacional de Astrofísica, Óptica y Electrónica, Calle Luis Enrique Erro No. 1, Sta. María Tonantzintla, Puebla 72840, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Photonics 2023, 10(2), 150; https://doi.org/10.3390/photonics10020150
Submission received: 19 October 2022 / Revised: 4 January 2023 / Accepted: 30 January 2023 / Published: 31 January 2023
(This article belongs to the Special Issue Lasers and Dynamic of Systems)

Abstract

:
We investigate some of the fundamental features of the interaction of a mixture of coherent states, namely, a Kaleidoscope states mixture, with two-level atoms in the Jaynes–Cummings model framework. We begin our analysis by calculating the von Neumann entropy of the field, which is determined with the help of the virtual atom method. The oscillations appearing in the entropy indicate a state of purity greater than the initial state, i.e., a purification of the initial state due to a transfer of coherence from the atom to the field. In this oscillatory region, we obtain a negative Wigner function that hints at a (noisy) multiple Schrödinger cat.

1. Introduction

The relation between classical and quantum phenomena has attracted attention over the years, becoming a very discussed problem in quantum mechanics. The main problem is that what causes, in the macroscopic world, the quantum interference of superposition states and the n-qubit entanglement states may not be observed [1,2,3].
Schrödinger cat [4,5,6] states, or superpositions of coherent states, have attracted the attention of researchers due to their fundamental features. Of particular interest is the case of the superposition of two (or more) coherent states [7,8], where, because of quantum interference, their properties are very different from the properties of the constituent coherent states, as well as from the incoherent superposition or statistical mixtures of such states. For example, the superposition exhibits sub-Poissonian photon statistics, higher-order squeezing, and oscillations in the photon number distribution [8]. These properties clearly differentiate the state of the superposition and statistical mixture of two coherent states [7]. Because superpositions of macroscopically distinguishable states (or Schrödinger cat-like states) may be produced by using coherent states, the problem is important for the quantum theory of measurement [9]. Several schemes have already been proposed to produce a superposition of coherent states, for instance the non-linear interaction of the field in a coherent state with a Kerr-like medium can produce their superposition [4]. Another possible way would be through the interaction between quantized fields, initially prepared in coherent states, with two-level atoms [7,8] or ion laser interactions [9].
On the other hand, entropy [9,10,11] is one of the main tools for measuring entanglement [2,12]. Precisely, one of the main tasks in the present manuscript is to calculate the entropy of the field for a mixture of Kaleidoscope states or a statistical mixture of coherent states, in the Jaynes–Cummings framework, which we will do with the aid of the Araki–Lieb inequality [13]. In principle, it seems impossible to use the Araki–Lieb inequality to calculate the field entropy because the field is in a statistical mixture of states [14]. However, via purification of the mixed density matrix of the quantized field [15], we will be able to use that inequality in order to calculate the field von Neumann entropy, either for the atom or the field.
In the next section, we define the Kaleidoscope states that we will consider as the initial state in their interaction with a two-level atom. Section 3 deals precisely with this interaction and, there, we calculate the field entropies and their Wigner functions [16], where we show that the statistical mixture of coherent states, i.e., the Kaleidoscope states mixture, goes to a minimum at half the revival time, independently of its number of components. Finally, in Section 4, we summarize our conclusions.

2. Kaleidoscope States

Kaleidoscope states are a particular superposition of n coherent states, and are defined in reference [17], as:
λ 1 α | ψ 1 α λ 2 α | ψ 2 α λ 3 α | ψ 3 α λ n α | ψ n α = 1 n 1 1 1 1 1 ω * ω * 2 ω * ( n 1 ) 1 ω * 2 ω * 4 ω * 2 ( n 1 ) 1 ω * ( n 1 ) ω * 2 ( n 1 ) ω * ( n 1 ) ( n 1 ) | α | ω α | ω 2 α | ω ( n 1 ) α ,
where λ k α are normalization constants and ω = exp i 2 π n with n, k integers which satisfy 1 k n .
We can observe that the Vandermonde matrix that transforms the vectors in Equation (1) is proportional to the so-called discrete Fourier transform, also known as the quantum Fourier transform [18].
On the other hand, an initial statistical mixture of n coherent states may be written as
ρ ^ F ( 0 ) = 1 n k = 1 n | ω ( k 1 ) α ω ( k 1 ) α | .
This density matrix may be diagonalized by the virtual atom method [15], where the important issue is to establish a connection between this density and the virtual ( V ) density operator ρ ^ V = | ψ V ψ V | , which is calculated from the pure state
| ψ V = λ 1 α | ψ 1 α | A 1 + λ 2 α | ψ 2 α | A 2 + λ 3 α | ψ 3 α | A 3 + + λ n α | ψ n α | A n ,
where { | A k } is the virtual atom basis and | ψ k α are the Kaleidoscope states, whenever k = 1 , . . , n . After tracing the density operator ρ ^ V = | ψ V ψ V | , over the virtual atom states { | A k } , we obtain the block diagonal density matrix
ρ ^ F ( 0 ) = k = 1 n λ k α | ψ k α ψ k α | ,
where
λ k α = e | α | 2 n m = 1 n ω * m ( k 1 ) exp ω m | α | 2 .
Similarly, by tracing the density ρ ^ V = | ψ V ψ V | over the field basis, we obtain the virtual atom (VA) density matrix
ρ ^ V A = 1 n 1 α | ω α * α | ω 2 α * α | ω ( n 1 ) α * α | ω α 1 α | ω α * α | ω ( n 2 ) α * α | ω 2 α α | ω α 1 α | ω ( n 3 ) α * α | ω ( n 1 ) α α | ω ( n 2 ) α α | ω ( n 3 ) α 1 ,
and, taking into account that α | ω k α = α | ω ( n k ) α * and α | ω k α = exp ( | α | 2 ) exp ( ω k | α | 2 ) , the Equation (6) may be rewritten as:
ρ ^ V A = exp ( | α | 2 ) n exp ( | α | 2 ) exp ( ω * | α | 2 ) exp ( ω * 2 | α | 2 ) exp ( ω * ( n 1 ) | α | 2 ) exp ( ω * ( n 1 ) | α | 2 ) exp ( | α | 2 ) exp ( ω * | α | 2 ) exp ( ω * ( n 2 ) | α | 2 ) exp ( ω * ( n 2 ) | α | 2 ) exp ( ω * ( n 1 ) | α | 2 ) exp ( | α | 2 ) exp ( ω * ( n 3 ) | α | 2 ) exp ( ω * | α | 2 ) exp ( ω * 2 | α | 2 ) exp ( ω * 3 | α | 2 ) exp ( | α | 2 ) ,
where
| A k 1 n 1 ω ( k 1 ) ω 2 ( k 1 ) ω 3 ( k 1 ) ω ( n 1 ) ( k 1 ) ,
are eigenvectors of ρ ^ V A , whose eigenvalues are given by (5). For details of the above results, we refer to the Appendixes Appendix A and Appendix B.

3. Interaction of Kaleidoscope States with a Two-Level Atom

The interaction between a quantized field and a two-level atom (under rotating wave approximation) is given by the Jaynes–Cummings interaction Hamiltonian [19] (for simplicity we have set = 1 ),
H ^ I = λ a ^ σ + a ^ σ + ,
where we have considered on-resonance conditions (equal field and atomic transition frequencies). In the above equation, λ is the coupling constant, a ^ and a ^ and the annihilation and creation operators, respectively, and σ + and σ are the raising and lowering Pauli operators, respectively. The evolution operator, U ^ = exp ( i H ^ I t ) , in the 2 × 2 basis is given by
U ^ = cos λ t a ^ a ^ i V ^ sin λ t a ^ a ^ i V ^ sin λ t a ^ a ^ cos λ t a ^ a ^ ,
with V ^ and V ^ the London phase operators [20]. For details of the above and next results, we refer to the appendices.
If we consider the atom initially prepared in the state | e and the field in a statistical mixture of n coherent states, i.e., the system is initially prepared in ρ ^ ( 0 ) = ρ ^ F ( 0 ) | e e | , with ρ ^ F ( 0 ) defined by (2), then the time evolved density matrix is given by
ρ ^ = k = 1 n | C k C k | | C n S 1 | + k = 2 n | C k 1 S k | | S 1 C n | + k = 2 n | S k C k 1 | k = 1 n | S k S k | ,
with
| C k = e | α | 2 2 m = 0 + α n m + k ( n m + k ) ! cos λ t n m + k + 1 | n m + k , ( k = 1 , 2 , n 1 ) ,
| C n = e | α | 2 2 m = 0 + α n m ( n m ) ! cos λ t n m + 1 | n m ,
| S k = i e | α | 2 2 m = 0 + α n m + k 1 ( n m + k 1 ) ! sin λ t n m + k | n m + k , ( k = 1 , 2 , n ) .
In the study of the interaction between two subsystems, namely atom (A) and field (F), the subsystem entropies S A and S F , and the total entropy S A F play an essential role and they obey the Araki–Lieb inequality [13],
| S A S F | S A F S A + S F .
Particularly, when the states of the subsystems are prepared in pure states, we have that S A F = 0 , and the entropies of the subsystems will be equal S A = S F . This situation is precisely our case because the states of the field (3) and atom are pure states, and we can apply the virtual atom method [15] in order to find the entropies.
Therefore, tracing over the field states, we obtain the reduced block diagonal virtual atom density matrix
ρ ^ V A = C 1 | C 1 C 1 | S 1 * 0 0 C 1 | S 1 S 1 | S 1 0 0 0 0 C n | C n C n | S n * 0 0 C n | S n S n | S n ,
whose eigenvalues are:
λ k ± = 1 2 C k | C k + S k | S k ± 1 2 C k | C k S k | S k 2 + 4 | C k | S k | 2 ,
for k = 1 , 2 , 3 , 4 , , n , such that the field entropy can be written in the form
S F = k = 1 n λ k + ln λ k + k = 1 n λ k ln λ k .
Finally, tracing the density (11) over the atomic states, we obtain the reduced density matrix for the field as
ρ ^ F = k = 1 n | C k C k | + k = 1 n | S k S k | .
As two coherent states are sufficiently apart when α 2 , they may be considered orthogonal. They may be considered orthogonal as α | α = exp ( 2 | α | 2 ) and, as α becomes larger, the exponential approaches zero. On the other hand, the Kaleidoscope-states are orthogonal for any α > 0 , and ρ ^ V A will be a block diagonalizable density as was shown in Equation (16).
The field entropy is plotted as a function of the scaled time λ t in Figure 1, for Kaleidoscope-states with different values of statistical mixture of coherent states, namely, n = 1 , 2 , 4 , 8 and 16, and α = 6.0 . Their initial values are calculated as ln ( n ) and are indicated in the figure by dotted lines. It may be seen that the entropies have similar behavior: each one possesses a global minimum of about λ t 19 ; for all values of n 2 the field entropy is below the dotted line corresponding to their initial values, making it clear that a purification process takes place. This may be explained as follows: the photon distribution for a mixture of coherent states is a Poissonian distribution, just as the one for a coherent state. This means that the collapses and revivals occur exactly at the same time for each of the components, namely t R / 2 = π | α | / λ , with t R the so-called revival time. This phenomenon happens because each of the components of the mixture does not interfere with any other components as the state considered is not a superposition of states but rather a statistical mixture. Once the interaction begins, each component is divided into two counter-rotating components (in phase space, see Figure 2) that produce the purest state precisely at half the revival time [21]. In this region and in each case, the field becomes purer than its initial state and oscillations appear for n 2 . Such purification occurs because of a transfer of coherence from the atom to the field [22].
Here, some questions arise: what will the form of the above-mentioned purer states for different values of a statistical mixture of coherent states n about λ t 19 be? What will the form of the state where the field entropy reaches its maximum value for different values of a statistical mixture of coherent states n be? In order to try to answer the above questions, we calculated the field Wigner function for the reduced density matrix field defined in Equation (19). The Wigner function may be written as [16]
W F ( β ) = 2 π Tr F ρ ^ F D ^ ( β ) ( 1 ) n ^ D ^ ( β ) , = 2 π k = 1 n C k | ρ ^ F D ^ ( β ) ( 1 ) n ^ D ^ ( β ) | C k + 2 π k = 1 n S k | ρ ^ F D ^ ( β ) ( 1 ) n ^ D ^ ( β ) | S k ,
where D ^ ( β ) = exp ( β a ^ β * a ^ ) is the displacement operator.
In Figure 2a–c, we show the field Wigner function at λ t = 0 , and we clearly see the two, four, and eight peaks corresponding to each coherent state for n = 2 , 4 and 8, respectively. Additionally, in Figure 2d–f, we show the field Wigner function corresponding to the Kaleidoscope-State with n = 2 , 4 and 8 respectively when the field would become a purer state at time λ t 19.15 . These Wigner functions resemble a Schrödinger cat state of 2 , 4 and 8 components, where we note the characteristic interference structure. We clearly see the formation of quantum interferences halfway between the n humps. As we mentioned above, once the interaction starts, each of the components is divided into two counter-rotating terms and, at half the revival time, the humps recombine with the traveling contributions from other component states; that, however, is not a coherent process, as each one of the components interferes only with itself. At that time, the Wigner function becomes negative, which is a clear signature of a non-classical state. This fact, plus the purification of the field shown in Figure 1, hints that multiple Schrödinger cats are being generated. The frequency of the interference structure increases as the separation distance α increases [23]. For example, setting α = 4 and n = 2 , the entropy has a similar behavior as in Figure 1 but, as we can see in Figure 3a, now its minimum is around λ t 12.5 , and its corresponding field Wigner function has an interference structure with a lower frequency, as is shown in Figure 4c. Finally, when λ t 19 , the initial mixture of Kaleidoscope States (as we showed in Figure 2a–c) gains purity as was suggested by the entropy behavior, and the negativity of the field Wigner function are an indicator of the non-classical properties of the state [23] at λ t 19.15 .
On the other hand, for the case of the maximum value of the field entropy λ t 2.0 as it is shown in Figure 3a, we see that each peak of the initial coherent state in Figure 4a, split into two counter-rotating terms, where each one of the components interferes only with itself and spread on the circle of radius α = 4 , and the Wigner function, may have negative values, as is shown in Figure 4b, which is an indicator of the non-classical properties of the state [23].
For the sake of completeness, the atom entropy is plotted in Figure 3b as a function of the scaled time λ t with n = 2 and α = 4.0 . This entropy quickly saturates to the value of ln ( 2 ) , and it may be seen that the entropy has similar behavior for other values of n and α .
The above results can be compared with previous ones in Figure 1 for n = 16 and α = 6 . In that case, it was shown that the field entropy has its minimum around λ t 19.15 . In Figure 4d, we show the field Wigner function at λ t = 0 , and we clearly see the sixteen peaks corresponding to each of n = 16 coherent states. For the case of the maximum value of the field entropy ( λ t 1.15 ), we see that each peak of the initial coherent state splits into two counter-rotating terms, where each component interferes only with itself, and it the formation of the quantum interference structure on the circle of radius α = 6 may also be seen. In this case, the Wigner function has no negative values. Finally, when the time goes to λ t 19.15 , the peaks interfere and the formation of the quantum interference structures may be seen halfway between the 16 humps, which is a clear signature of a non-classical state.

4. Conclusions

We have shown that the Jaynes–Cummings interaction with an initial Kaleidoscope mixture of coherent states may be modeled by the virtual Hamiltonian method by extending the atomic Hilbert space such that a virtual pure state may be associated as an initial wavefunction. In particular, we have seen that the purification procedure takes us from a mixed field density matrix to a pure wave function that involves a virtual 2 n -level atom, as can be seen in the 2 n terms in Equation (19). Finally, we should mention that the effects presented in the field entropy for the initial field state given by a statistical mixture of the constituent states are reflected in the appearance of oscillations that give rise to negative Wigner functions resembling multiple Schrödinger’s cats.

Author Contributions

Conceptualization, J.A.A.-C.; investigation, all authors; methodology, all authors; software, all authors; validation, all authors; formal analysis, all authors; writing—original draft preparation, all authors; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Initial Field Density

Following the virtual atom method to the density matrix (2) we propose:
| ψ V = | α n | a 1 + | ω α n | a 2 + | ω 2 α n | a 3 + + | ω ( n 1 ) α n | a n ,
where for simplicity we represented the ket in the Hilbert space by the standard basis of n dimensional Euclidean space as:
| a 1 = 1 0 0 0 , | a 2 = 0 1 0 0 , | a 3 = 0 0 1 0 , , | a n = 0 0 0 1 ,
(we do here as it is done when dealing with Pauli-spin operators: one may go to the 2 × 2 matrix form or work with the well known σ operators). From Equation (1) and considering orthogonal properties of Discrete Fourier Transform (DFT), one has
| α / n | ω α / n | ω 2 α / n | ω ( n 1 ) α / n = 1 n 1 1 1 1 1 ω ω 2 ω ( n 1 ) 1 ω 2 ω 4 ω 2 ( n 1 ) 1 ω ( n 1 ) ω 2 ( n 1 ) ω ( n 1 ) ( n 1 ) λ 1 α | ψ 1 α λ 2 α | ψ 2 α λ 3 α | ψ 3 α λ n α | ψ n α .
Thus, according to the above formula (A3), we can rewrite Equation (A1) as:
| ψ V = 1 n λ 1 α | ψ 1 α + λ 2 α | ψ 2 α + + λ n α | ψ n α | a 1 + 1 n λ 1 α | ψ 1 α + ω λ 2 α | ψ 2 α + + ω ( n 1 ) λ n α | ψ n α | a 2 + 1 n λ 1 α | ψ 1 α + ω 2 λ 2 α | ψ 2 α + + ω 2 ( n 1 ) λ n α | ψ n α | a 3 + 1 n λ 1 α | ψ 1 α + ω ( n 1 ) λ 2 α | ψ 2 α + + ω ( n 1 ) ( n 1 ) λ n α | ψ n α | a n = λ 1 α | ψ 1 α 1 n | a 1 + | a 2 + | a 3 + + | a n + λ 2 α | ψ 2 α 1 n | a 1 + ω | a 2 + ω 2 | a 3 + + ω ( n 1 ) | a n + λ n α | ψ n α 1 n | a 1 + ω ( n 1 ) | a 2 + ω 2 ( n 1 ) | a 3 + + ω ( n 1 ) ( n 1 ) | a n = λ 1 α | ψ 1 α | A 1 + λ 2 α | ψ 2 α | A 2 + + λ n α | ψ n α | A n ,
where
| A 1 | A 2 | A 3 | A n = 1 n 1 1 1 1 1 ω ω 2 ω ( n 1 ) 1 ω 2 ω 4 ω 2 ( n 1 ) 1 ω ( n 1 ) ω 2 ( n 1 ) ω ( n 1 ) ( n 1 ) | a 1 | a 2 | a 3 | a n ,
with
| A k = 1 ω ( k 1 ) ω 2 ( k 1 ) ω ( n 1 ) ( k 1 ) = a 1 | A k a 2 | A k a 3 | A k a n | A k ,
then these new virtual atom states { | A i i } are orthogonal kets because the discrete Fourier transform (DFT) is an orthogonal matrix too.
Using the above results it can be shown that the density matrix associated with | ψ V defined by Equation (A1) will be
ρ ^ V = 1 n i = 1 n j = 1 n | ω ( i 1 ) α | a i a j | ω ( j 1 ) α | = i = 1 n j = 1 n λ i α λ j α * | ψ i α | A i A j | ψ j α | ,
where we have used relation (A4). Taking the trace of the density matrix (A7) over the field states we get the virtual atom density operator
ρ ^ V A = 1 n i = 1 n j = 1 n ω ( j 1 ) α | ω ( i 1 ) α | a i a j | = i = 1 n j = 1 n λ i α λ j α * ψ j α | ψ i α | A i A j | .
But in reference [17], it was remarked that Kaleidoscope states are orthogonal, then
ρ ^ V A = 1 n i = 1 n j = 1 n ω ( j 1 ) α | ω ( i 1 ) α | a i a j | = i = 1 n λ i α | A i A i | .
We will now do an analysis of the eigenvalue problem for ρ ^ V A , which we may start writing
ρ ^ V A | A k = 1 n i = 1 n j = 1 n ω ( j 1 ) α | ω ( i 1 ) α | a i a j | A k = λ k α | A k , = e | α | 2 n i = 1 n j = 1 n exp ω * ( j i ) | α | 2 1 n ω ( j 1 ) ( k 1 ) | a i = λ k α | A k , = b 1 b 2 b 3 b n . ,
where
b l = e | α | 2 n j = 1 n exp ω * ( j l ) | α | 2 1 n ω ( j 1 ) ( k 1 ) , = e | α | 2 n 1 n ω ( l 1 ) ( k 1 ) m = 1 l n l ω m ( k 1 ) exp ω * m | α | 2 , = e | α | 2 n a l | A k m = l n + l 1 ω * m ( k 1 ) exp ω m | α | 2 ,
with 1 l n . When l = 1 , we have
b 1 = e | α | 2 n a 1 | A k m = 1 n ω * m ( k 1 ) exp ω m | α | 2 = λ k α a 1 | A k ,
as a 1 | A k 0 , then we have eigenvalue
λ k α = e | α | 2 n m = 1 n ω * m ( k 1 ) exp ω m | α | 2 .
In a similar fashion, for l > 1 we have
b l = e | α | 2 n a l | A k m = l n + l 1 ω * m ( k 1 ) exp ω m | α | 2 , = e | α | 2 n a l | A k m = 1 n + l 1 ω * m ( k 1 ) exp ω m | α | 2 e | α | 2 n a l | A k m = 1 l 1 ω * m ( k 1 ) exp ω m | α | 2 , = e | α | 2 n a l | A k m = 1 n ω * m ( k 1 ) exp ω m | α | 2 e | α | 2 n a l | A k m = 1 l 1 ω * m ( k 1 ) exp ω m | α | 2 + e | α | 2 n a l | A k m = n + 1 n + l 1 ω * m ( k 1 ) exp ω m | α | 2 , = e | α | 2 n a l | A k m = 1 n ω * m ( k 1 ) exp ω m | α | 2 e | α | 2 n a l | A k m = 1 l 1 ω * m ( k 1 ) exp ω m | α | 2 + e | α | 2 n a l | A k m = 1 l 1 ω * ( m + n ) ( k 1 ) exp ω m + n | α | 2 .
Because ω n = exp i 2 π n n = 1 , we can cancel the last two sums in Equation (A14), and we have that
b l = e | α | 2 n a l | A k m = 1 n ω * m ( k l ) exp ω m | α | 2 = λ k α a l | A k ,
as a l | A k 0 , then we have eigenvalue
λ k α = e | α | 2 n m = 1 n ω * m ( k l ) exp ω m | α | 2 .
With the help of Equations (A12) and (A15), the eigenvalue Equation (A10) can be written as,
ρ ^ V A | A k = λ k α a 1 | A k a 1 | A k a n | A k a n | A k = λ k α | A k .
Now, taking the trace of the density matrix (A7) over field states we get the field density:
ρ ^ F = 1 n i = 1 n | ω ( i 1 ) α ω ( i 1 ) α | = i = 1 n λ i α | ψ i α ψ i α | .

Appendix B. Jaynes-Cummings Dynamics

The Jaynes-Cummings time evolution operator is given by
U ^ ( t ) = cos λ t a ^ a ^ i V ^ sin λ t a ^ a ^ i V ^ sin λ t a ^ a ^ cos λ t a ^ a ^ ,
where V ^ and V ^ are the London phase operators [20]. If we consider the initial density matrix ρ ^ ( 0 ) = ρ ^ F ( 0 ) | e e | , where ρ ^ F ( 0 ) was defined in Equation (4), we obtain:
ρ ^ ( t ) = U ^ ( t ) ρ ^ ( 0 ) U ^ ( t ) = U ^ ( t ) ρ ^ F ( 0 ) 0 0 0 U ^ ( t ) , = k = 1 n | C k C k | | C n S 1 | + k = 2 n | C k 1 S k | | S 1 C n | + k = 2 n | S k C k 1 | k = 1 n | S k S k | ,
where,
| C k = λ ( k + 1 ) α cos λ t a ^ a ^ | ψ ( k + 1 ) α , ( k = 1 , 2 , 3 , , n 1 ) , | C n = λ 1 α cos λ t a ^ a ^ | ψ 1 α , | S k = i λ k α V ^ sin λ t a ^ a ^ | ψ k α , ( k = 1 , 2 , 3 , , n ) .
Nevertheless, Kaleidoscope states can be rewritten as:
| ψ k α = e | α | 2 / 2 λ k α m = 0 + α n m + k 1 ( n m + k 1 ) ! a ^ ( n m + k 1 ) | 0 , = e | α | 2 / 2 λ k α m = 0 + α n m + k 1 ( n m + k 1 ) ! | n m + k 1 .
With the above relation, we can write
| C k = e | α | 2 2 m = 0 + α n m + k ( n m + k ) ! cos λ t n m + k + 1 | n m + k , ( k = 1 , 2 , 3 , , n 1 ) , | C n = e | α | 2 2 m = 0 + α n m ( n m ) ! cos λ t n m + 1 | n m | S k = i e | α | 2 2 m = 0 + α n m + k 1 ( n m + k 1 ) ! sin λ t n m + k | n m + k , ( k = 1 , 2 , 3 , , n ) .
From the time evolution density matrix (A20), we can obtain the reduced block diagonal atom density matrix by tracing over the field states, so that
ρ ^ A = C 1 | C 1 C 1 | S 1 * 0 0 C 1 | S 1 S 1 | S 1 0 0 0 0 C n | C n C n | S n * 0 0 C n | S n S n | S n ,
whose characteristic polynomial is:
det ρ ^ A λ I ^ = k = 1 n det C k | C k λ C k | S k * S k | C k S k | S k λ = 0 ,
and we get the eigenvalues
λ k ± = 1 2 C k | C k + S k | S k ± 1 2 C k | C k S k | S k 2 + 4 | C k | S k | 2 ,
for k = 1 , 2 , 3 , , n . By recalling the properties of the trace of the block diagonal matrix we calculate the atomic von Neumann entropy defined as
S A = Tr ρ ^ A ln ρ ^ A , = k = 1 n C k | C k ln k = 1 n C k | C k k = 1 n S k | S k ln k = 1 n S k | S k , = k = 1 n λ k + ln λ k + k = 1 n λ k ln λ k .
Finally, by tracing the density (A20) over the atomic states, we obtain the reduced density matrix for the field as
ρ ^ F = k = 1 n | C k C k | + k = 1 n | S k S k | .

References

  1. Man, Z.X.; Zhang, Y.J.; Su, F.; Xia, Y.J. Entanglement dynamics of multiqubit system in Markovian and non-Markovian reservoirs. Eur. Phys. J. D 2010, 58, 147–151. [Google Scholar] [CrossRef]
  2. Nourmandipour, A.; Tavassoly, M.K.; Rafiee, M. Dynamics and protection of entanglement in n-qubit systems within Markovian and non-Markovian environments. Phys. Rev. A 2016, 93, 022327. [Google Scholar] [CrossRef] [Green Version]
  3. Doll, R.; Wubs, M.; Hänggi, P.; Kohler, S. Incomplete pure dephasing of N-qubit entangled W states. Phys. Rev. B 2007, 76, 045317. [Google Scholar] [CrossRef] [Green Version]
  4. Yurke, B.; Stoler, D. Generating quantum mechanical superpositions of macroscopically distinguishable states via amplitude dispersion. Phys. Rev. Lett. 1986, 57, 13–16. [Google Scholar] [CrossRef] [PubMed]
  5. Brune, M.; Haroche, S.; Raimond, J.M.; Davidovich, L.; Zagury, N. Manipulation of photons in a cavity by dispersive atom-field coupling: Quantum-nondemolition measurements and generation of “Schrödinger cat” states. Phys. Rev. A 1992, 45, 5193–5214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Davidovich, L. From quantum to classical: Schrödinger cats, entanglement, and decoherence. Phys. Scr. 2016, 91, 063013. [Google Scholar] [CrossRef] [Green Version]
  7. Moya-Cessa, H.; Vidiella-Barranco, A. On the Interaction of Two-level Atoms with Superpositions of Coherent States of Light. J. Mod. Opt. 1995, 42, 1547–1552. [Google Scholar] [CrossRef]
  8. Vidiella-Barranco, A.; Moya-Cessa, H.; Bužek, V. Interaction of Superpositions of Coherent States of Light with Two-level Atoms. J. Mod. Opt. 1992, 39, 1441–1459. [Google Scholar] [CrossRef]
  9. Barnett, S.; Phoenix, S. Entropy as a measure of quantum optical correlation. Phys. Rev. A 1989, 40, 2404–2409. [Google Scholar] [CrossRef]
  10. Phoenix, S. Wave-packet evolution in the damped oscillator. Phys. Rev. A 1995, 41, 5132–5138. [Google Scholar] [CrossRef]
  11. Barnett, S.M.; Beige, A.; Ekert, A.; Garraway, B.M.; Keitel, C.H.; Kendon, V.; Lein, M.; Milburn, G.J.; Moya-Cessa, H.M.; Murao, M.; et al. Journeys from quantum optics to quantum technology. Prog. Quantum Electron. 2017, 54, 19–45. [Google Scholar] [CrossRef]
  12. Golkar, G.; Tavassoly, M. Atomic motion and dipole–dipole effects on the stability of atom–atom entanglement in Markovian/non-Markovian reservoir. Mod. Phys. Lett. A 2019, 34, 1950077. [Google Scholar] [CrossRef]
  13. Araki, H.; Lieb, E.H. Entropy inequalities. Commun. Math. Phys. 1970, 18, 160–170. [Google Scholar] [CrossRef]
  14. Zúñiga-Segundo, A.; Juárez-Amaro, R.; Aguilar-Loreto, O.; Moya-Cessa, H.M. Field’s entropy in the atom–field interaction: Statistical mixture of coherent states. Ann. Phys. 2017, 379, 150–158. [Google Scholar] [CrossRef]
  15. Anaya-Contreras, J.A.; Zúñiga-Segundo, A.; Moya-Cessa, H.M. Purification of the atom-field interaction Hamiltonian. Phys. Open 2019, 1, 100007. [Google Scholar] [CrossRef]
  16. Moya-Cessa, H.; Knight, P.L. Series representation of quantum-field quasiprobabilities. Phys. Rev. A 1993, 48, 2479–2481. [Google Scholar] [CrossRef]
  17. Koçak, A.; Pashaev, O.K. Special functions with mod N Symmetry Kaleidosc. Quantum Coherent States. J. Phys. Conf. Ser. 2019, 1194, 012059. [Google Scholar] [CrossRef]
  18. Perez-Leija, A.; Andrade-Morales, L.A.; Soto-Eguibar, F.; Szameit, A.; Moya-Cessa, H.M. The Pegg–Barnett phase operator and the discrete Fourier transform. Phys. Scr. 2016, 91, 043008. [Google Scholar] [CrossRef] [Green Version]
  19. Jaynes, E.; Cummings, F. Comparison of quantum and semiclassical radiation theories with application to the beam maser. Proc. IEEE 1963, 51, 89–109. [Google Scholar] [CrossRef] [Green Version]
  20. Susskind, L.; Glogower, J. Quantum mechanical phase and time operator. Phys. Phys. Fiz. 1964, 1, 49–61. [Google Scholar] [CrossRef] [Green Version]
  21. Gea-Banacloche, J. Atom- and field-state evolution in the Jaynes-Cummings model for large initial fields. Phys. Rev. A 1991, 44, 5913–5931. [Google Scholar] [CrossRef] [Green Version]
  22. Arancibia-Bulnes, C.A.; Moya-Cessa, H.; Sánchez-Mondragón, J.J. Purifying a thermal field in a lossless micromaser. Phys. Rev. A 1995, 51, 5032–5034. [Google Scholar] [CrossRef]
  23. Kenfack, A.; yczkowski, K. Negativity of the Wigner function as an indicator of non-classicality. J. Opt. B Quantum Semiclassical Opt. 2004, 6, 396–404. [Google Scholar] [CrossRef]
Figure 1. The evolution of the field entropy of Kaleidoscope-States as a function of the scaled time λ t and different values of statistical mixture of coherent states n = 1 , 2 , 4 , 8 and 16, with α = 6.0 . Their initial values are calculated as ln ( n ) and are indicated in the figure by dotted lines.
Figure 1. The evolution of the field entropy of Kaleidoscope-States as a function of the scaled time λ t and different values of statistical mixture of coherent states n = 1 , 2 , 4 , 8 and 16, with α = 6.0 . Their initial values are calculated as ln ( n ) and are indicated in the figure by dotted lines.
Photonics 10 00150 g001
Figure 2. Wigner function for Kaleidoscope-State with α = 6.0 for several values of the time λ t and n. For λ t = 0 : (a) n = 2 , (b) n = 4 and (c) n = 8 ; for λ t 19.15 : (d) n = 2 , (e) n = 4 and (f) n = 8 .
Figure 2. Wigner function for Kaleidoscope-State with α = 6.0 for several values of the time λ t and n. For λ t = 0 : (a) n = 2 , (b) n = 4 and (c) n = 8 ; for λ t 19.15 : (d) n = 2 , (e) n = 4 and (f) n = 8 .
Photonics 10 00150 g002
Figure 3. The evolution of the: (a) field entropy and (b) and atom entropy; for a Kaleidoscope-State as a function of the scaled time λ t with n = 2 and α = 4.0 . The initial value of field entropy is calculated as ln ( 2 ) and is indicated in the figure by a dotted line.
Figure 3. The evolution of the: (a) field entropy and (b) and atom entropy; for a Kaleidoscope-State as a function of the scaled time λ t with n = 2 and α = 4.0 . The initial value of field entropy is calculated as ln ( 2 ) and is indicated in the figure by a dotted line.
Photonics 10 00150 g003
Figure 4. Field Wigner function for Kaleidoscope-State for several values of the time λ t and n. For n = 2 (a) λ t = 0 , (b) λ t 2.0 and (c) λ t 12.5 . And for n = 16 (d) λ t = 0 , (e) λ t 1.15 and (f) λ t 19.15 .
Figure 4. Field Wigner function for Kaleidoscope-State for several values of the time λ t and n. For n = 2 (a) λ t = 0 , (b) λ t 2.0 and (c) λ t 12.5 . And for n = 16 (d) λ t = 0 , (e) λ t 1.15 and (f) λ t 19.15 .
Photonics 10 00150 g004
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.

Share and Cite

MDPI and ACS Style

Anaya-Contreras, J.A.; Zúñiga-Segundo, A.; Moya-Cessa, H.M. Statistical Mixture of Kaleidoscope States Interacting with a Two-Level Atom: Entropy and Purification. Photonics 2023, 10, 150. https://doi.org/10.3390/photonics10020150

AMA Style

Anaya-Contreras JA, Zúñiga-Segundo A, Moya-Cessa HM. Statistical Mixture of Kaleidoscope States Interacting with a Two-Level Atom: Entropy and Purification. Photonics. 2023; 10(2):150. https://doi.org/10.3390/photonics10020150

Chicago/Turabian Style

Anaya-Contreras, Jorge A., Arturo Zúñiga-Segundo, and Héctor M. Moya-Cessa. 2023. "Statistical Mixture of Kaleidoscope States Interacting with a Two-Level Atom: Entropy and Purification" Photonics 10, no. 2: 150. https://doi.org/10.3390/photonics10020150

APA Style

Anaya-Contreras, J. A., Zúñiga-Segundo, A., & Moya-Cessa, H. M. (2023). Statistical Mixture of Kaleidoscope States Interacting with a Two-Level Atom: Entropy and Purification. Photonics, 10(2), 150. https://doi.org/10.3390/photonics10020150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop