Nonclassical Effects Based on Husimi Distributions in Two Open Cavities Linked by an Optical Waveguide

Nonclassical effects are investigated in a system formed by two quantum wells, each of which is inside an open cavity. The cavities are spatially separated, linked by a fiber, and filled with a linear optical medium. Based on Husimi distributions (HDs) and Wehrl entropy, we explore the effects of the physical parameters on the generation and the robustness of the mixedness and HD information in the phase space. The generated quantum coherence and the HD information depend crucially on the cavity-exciton and fiber cavity couplings as well as on the optical medium density. The HD information and purity are lost due to the dissipation. This loss may be inhibited by increasing the optical susceptibility as well as the couplings of the exciton-cavity and the fiber-cavity. These parameters control the regularity, amplitudes, and frequencies of the generated mixedness.

Towards understanding the dynamics of the nonclassical effects, the Husimi quasiprobability distributions are utilized as the preferred tool to investigate the phase space information and the mixedness. Phase-space Husimi distributions (HDs) [33] are usually used to investigate quantum effects. They are always positive distributions [34], and they are crucial tools for exploring the

Physical Model
Our considered system is formed by two open cavities (A and B). These cavities are linked by a waveguide and filled by a linear optical medium. Each cavity also contains a quantum well placed in a position corresponding to the maximum of the electromagnetic field [47][48][49]. We focus on the case where: (1) the driving lasers are in incident pumping; and (2) the laser, the cavities, the quantum wells, and the fiber modes have the same frequency (resonance). Additionally, we concentrate on the case of the short fiber interaction, which can be realized in usual experimental situations [39,50]. The condition of the short fiber limit is defined by 2lν 2πc 1, where l designates the fiber length, c is the light speed in the fiber, andν is the cavity dissipation. The interaction Hamiltonian is given by: The last two terms describe the interaction between the two exciton-photon (EP) systems and the fiber mode.f is the annihilation operator for the optical fiber mode, and λ f designates the coupling constants between the optical fiber and the microcavity modes.Ĥ k Int−EP is the interaction Hamiltonian of the i-EP system,Ĥ whereâ i (â † i ) designate the annihilation (creation) operators of the cavities. χ k = χ (k = A, B) present the linear optical susceptibility of the k-optical medium inside the cavities.b k is the annihilation excitonic operator. The second term of the Hamiltonian represents the k-EP interaction with the coupling constants λ k (k = A, B). ε k designates the amplitudes of the k-pumps laser.
The master equation of the system is (we takeh = 1) [51]: The γ k and κ k represent the normalized excitonic spontaneous emission and the cavity dissipation rates, respectively.
For weak pumping and weak dissipations, the off-diagonal terms 2â i ρâ † i and 2b i ρb † i can be neglected in Equation (3) [52,53]. Therefore, the density matrix can then be factorized, and Equation (3) becomes: where the non-Hermitian operator H E has the following expression: From Equation (4), the dynamics of the system density matrix,ρ(t), are governed by the wave function |ψ(t) as: The wave function |ψ(t) satisfies: The model of Equation (1) describes two exciton-photon systems interacting with a fiber field; therefore, the total number of the excitations is identified by (f †f +â † iâ i +b † ib i ). The Hilbert space of the wave function of the system is assumed under the the following conditions: (1) the number of excitations is restricted to 3; (2) only single photon processes are considered. Therefore, the two-exciton basis space is: and the wave function |ψ(t) can be represented as: where |β i = |Cav A , Cav B , Fib designate the cavities-fiber states, and are given by: |β 1 = a 1 |000 + a 9 |010 + a 13 |100 + a 17 |110 Using Equation (7), the dynamics of the time-dependent coefficients a n (t)(n = 1 − 18) are governed by the following differential equations: whereγ i = γ i 2 . λ k are the coupling rates between excitons and photons, while ε k are the laser pump amplitudes. In the numerical calculations, we take γ i = κ i = κ. The phase space information and the mixture will be studied via the Husimi distributions. In the two-exciton basis states {| i }(i = 1 − 4), the reduced density matrix of the two-exciton system is represented as: Then the reduced density matrices of the k-exciton (k = A, B) are defined by After tracing the states of the cavity and fiber fields and the B-exciton, we quantify the Husimi distribution and Wehrl entropy of the A-exciton to explore the phase space information and the mixedness of the A-exciton. In the numerical calculations, the time t represents a unitless time normalized to τ c , where τ c is the round trip time for a photon inside the cavity. We also normalize the parameters of the system to 1/τ c as: The the Husimi distribution and Wehrl entropy are investigated under the effects of the mentioned physical parameters after applying the time normalization for τ c of the round-trip time of photons in the cavity. For our investigation, the values of these parameters were chosen so that they are close to that of a typical experiment [54][55][56]. The choice of λ i between 1 and 0.004 gives us a margin for the Rabi frequency between 20 THz and 80 GHz. The choice of the dissipations γ i = κ i = 0.025 gives a lifetime of excitons of 0.34ps and a lifetime of photons of 0.17ps (equivalent to a reflection rate of the sample of 95%).

Phase Space Information of the Husimi Distribution
The HD is a tool to investigate the phase space information for a spin system. Here, we take two 1 2 -spins (two-exciton system). The phase space information is determined by the angles θ and φ [37]. The loss of information means that the HD is invariant in the phase space, and thus it is an inductor of the mixture.
Using the basis: {| i }, the two-exciton Bloch coherent states can be written as [57]: Consequently, the HD of the two-exciton system ρ E s (t) is: and its partial HD of the k-exciton (say k = A) is given by: where σ A i (i = x, y, z) represent the the expectation values of the Pauli matrices for the A-exciton associated with its reduced density matrix, where ρ ij (t) = i |ρ E s (t)| j are the elements of the two-exciton state ρ E s (t). Figures 1 and 2 show the dependence of the HD of the A-exciton H(θ, φ) = H k (θ, φ) on the evolution parameter, the external environment, the cavity-exciton, and the fiber-cavity coupling rates as well as on the optical susceptibility.    H(θ, φ, t) of the A-exciton when the total system is initially in the state |110 (a 17 = 1) where θ ∈ [0, 4π] and φ ∈ [0, 2π]. The values of the evolution parameter t were chosen based on the results of the Wehrl entropy mixedness. The A-exciton is in a pure state at t p = 0, 2.25π, ..., and it is in a mixed state at several times, one of them being t m = 2.6101π (see Figure 3a). To show the effect of the cavity-exciton interaction on the initial Husimi distribution H(θ, φ, t p ), the Husimi distribution is displayed in Figure 1a. The Husimi distribution of the initial A-exciton state ρ A (0) = |1 A 1 A | has a regular oscillatory behavior. We note that the 3D-Husimi distribution H(θ, φ) depends only on the angle θ, and that it has a regular 2π-oscillatory surface. The maxima of the HD are at θ = 2nπ, (n = 1, 2, ...) whereas its minima are at θ = (2n + 1)π, (n = 0, 1, 2, ...). In this case, the A-exciton HD oscillate between its extreme values 1+cos θ 4π . For θ = 0 we have H min = 0, and for θ = π we get H max = 1 2π 0.16. The counter-plot of the HD appears as color lines in the θ − φ plane, while parallel lines of the φ-axis illustrate the dependence of the maxima and minima of the HD on the phase space angular variables. Figure 1b, at t m = 2.6101π, shows the effect of the [cavity-exciton]-fiber-[cavity-exciton] interactions, with the strong couplings λ A = λ B = λ f = 1, on the A-exciton HD that leads to a decreasing of the amplitudes of the HD, where the maxima and minima of the HD tend to the stationary HD surface H(θ, φ) = 1 2 H max 0.08 (that appears clear with a large value of damping). The dependence of the HD on the θ-angle can be controlled by the evolution parameter t. Figure 1c, shows the Husimi distribution for an optical susceptibility χ = 2 with the same taken strong interaction couplings in Figure 1a. In this case, we observe that the oscillation frequency of HD, due to the unitary interactions, is at the same values as the case χ = 0. While the oscillation amplitudes of HD increase. The linear optical susceptibility leads to an increase of the dependence of the HD on the θ-angle, and the maxima and minima of the HD return to their initial values approximately. We deduce that when the cavities are filled by a linear optical medium, the phase space information is more robust against the interactions. Figure 1d shows the effect of the dissipation cavities on the maxima and minima of HD for the case where the same strong interaction couplings taken in Figure 1a. The amplitudes of the HD oscillations are affected, where the maxima and minima of HD remain at the same stationary value. The HD oscillations vanish after incasing the dissipations. The increase of the dissipation leads to a stationary HD surface of a 0.8 value (see Figure 1d); therefore, it leads to an erasing of the phase space information.
Based on the property that HD depends only on the angle θ, we can use the partial 2D-Husimi distribution H(θ) = H(θ, π) for the fixed value φ = π to analyze the phase space information under the effects of the external environment, the cavity-exciton, and the fiber-cavity coupling rates as well as the optical susceptibility; see Figure 2.
In Figure 2a, we display the dependence of the 2D-Husimi distribution H(θ) on the angle θ when the coupling strength of the fiber-cavity interaction is reduced to λ f = 0.2. We observe that the 2D-Husimi distribution H(θ) has regular oscillations with a 2π-period around the stationary value H(∞) = 0.08. For the case χ = 0 and λ f = 0.2 and the maxima of the HD are at θ = (2n + 1)π, (n = 0, 1, 2, ...) whereas its minima are at θ = 2nπ, (n = 1, 2, ...); this is an opposite behavior to the case where the strong fiber-cavity interaction coupling λ f = 1.
The dashed curve shows the effect of the optical susceptibility χ = 2 for the case λ f = 0.2. We note the same observations of the case λ f = 1. Dash-dot and doted curves show the effect of the cavity dissipation. The observations of the effects of the external environment and the optical susceptibility, for λ f = 1, are confirmed for λ f = 0.2. Figure 2b shows the sensitivity of the HD oscillations to the exciton-cavity coupling constants. If they are reduced to λ A = λ B = 0.2, the amplitudes of the HD oscillations present notable changes. By comparing the cases (λ A , λ B , λ f ) = (0.2, 0.2, 1) and (λ A , λ B , λ f ) = (1, 1, 0.2), we find that the amplitudes of the HD oscillations increase due to enhancement of the optical susceptibility. We find that the weak exciton-cavity couplings lead to an increasing of the dependence of the HD on the θ-angle, and the phase space information is more robust against the interactions.
Finally, we can deduce that the HD oscillations may be controlled by the exciton-cavity and fiber-cavity coupling constants as well as by the susceptibility of the optical medium.

Wehrl Entropy
One of the most important applications of the Husimi distribution in quantum information is the Wehrl entropy (WE) [35]. It is used to quantify the mixture (purity loss/mixedness) for open systems in the phase space determined by the angels θ ∈ [0, π] and φ ∈ [0, 2π] [58]. The partial Wehrl density of the k-exciton D k (θ k , φ k , t), is defined by: from this expression we can deduce the partial WEs of the A-exciton as: If the two excitons start in the lower state | ↓ A ↓ B , then the initial value of the partial WE is calculated as [57]: The WE function S(t) [57] verifies, The initial value of S(0) is for a pure state. The maximal value S max = ln(4π) indicates that the state is totally mixed. The Wehrl entropy is a good quantifier that provides important information about the transition processes (mixedness) between a pure state and a steady mixed state. In the real quantum devices, the decoherence is unavoidable, which destroys the pure state. In other words, the Wehrl entropy mixedness is the measure of the environment effect on the system.
For the case of one-qubit systems, there is a relation between the Wehrl entropy and the other entropies as von Neumann entropy and linear entropy [59]. It was proven that they exhibit similar behaviors.
In Figure 3, the WE function S(t) is plotted for different values of λ i (A, B, f ), χ and κ. The solid curve of Figure 3 shows the effect of the strong [cavity-exciton]-fiber-[cavity-exciton] interaction couplings λ A = λ B = λ f = 1 in the absence of the optical medium and the dissipation. The Wehrl entropy S(t) shows that the interactions lead to a generation of the mixedness of the A-exciton with a quasi-regular oscillatory behavior. In each period, the Wehrl entropy reaches its maximum value S max (t) = ln(4π), and remains there for a short interval. In these stability intervals, the A-exciton is in a maximally mixed state. The generated maximally mixed state as an initial state can be useful to realize quantum information processes. The Wehrl entropy returns instantaneously to its initial value at the end of each period (this means that the A-exciton returns to its initial pure state |1 A ).   The dashed curve of Figure 3a shows how the optical susceptibility affects the dynamical behavior of the generated A-exciton mixedness. The amplitudes and the period of the oscillations of the Wehrl entropy are decreased by an increasing of the optical susceptibility, i.e., the generated mixedness is reduced due to the increase of the optical susceptibility.
The dash-dot and doted curves show the effect of the external environment on the generated A-exciton mixedness. The increase of the cavity dissipations leads to the increase of the minima of the Wehrl entropy. The WE function S(t) reaches stationary behavior quickly, i.e., it leads to an acceleration of the A-exciton, causing it to be in a steady mixed state for a long time. The generated maximally mixed state can be used as an initial state to realize quantum computing processes [60,61] and qubit-channel metrology [62].
The effect of the strength of the fiber-cavity interaction on the generated mixedness is shown in Figure 3b. The Wehrl entropy is displayed when the coupling strength of the fiber-cavity interaction is reduced to λ f = 0.2. We observe that the WE function S(t) has a regular oscillatory behavior with a π period. The period, the amplitude, and the stability intervals of the generated mixedness are less than the case where λ f = 1. The dashed curve of Figure 3b shows that the Wehrl entropy is improved by increasing the optical susceptibility, unlike the one in the case of the strong fiber-cavity interaction coupling. In the absence of the dissipation (see the dash-dot and doted curves of Figure 3b), we note that the decrease of the fiber-cavity coupling leads to delaying the WE from reaching the stationary value S(∞) ≈ ln(4π) and the steady mixed state.
Our physical interpretations of the HD-information and the Wehrl entropy-mixedness in qubit-cavity systems linked by a waveguide mode contribute to updating more potential applications for distributed quantum information processing that are based on realizing a multi-target-qubit unconventional geometric phase gate in a multi-cavity system [40], as well as entanglement stabilization protocols between two superconducting qubits that are coupled to distant cavities [41] and the architecture of the cavity-based quantum network [42]. It is found that certain coupling parameters of the exciton-cavity and the fiber-cavity interactions control the effects of the optical linear medium and the environment, which are major obstacles to the applications of the quantum information technology.

Conclusions
In this investigation, we have studied the Husimi distributions and Wehrl entropy in a system formed by two open cavities linked by a waveguide. Each cavity contains one exciton and is filled by an optical linear medium. Based on Husimi distributions and Wehrl entropy, we have explored the effects of the optical susceptibility, the cavity dissipation, the coupling of the exciton-cavity and the fiber-cavity interactions. These parameters lead to notable changes in the dynamical behavior of the HD information and the generated mixedness. The generated HD information and the qubit purity are lost completely due to the excitonic spontaneous emission and the cavity dissipations. The mixedness and the HD information may be inhibited by increasing the optical susceptibility and the couplings of the exciton-cavity and the fiber-cavity. Therefore, we can deduce that the HD-information and mixedness depend crucially on the cavity-exciton and fiber cavity interactions as well as on the optical medium density. These physical parameters control the regularity, amplitudes, and frequencies of the generated mixedness. Since the studied system is very compact and can be realized experimentally, we expect that it would be very useful in the design of quantum components involved in quantum information processes.