Estimation of the Influence of a Noisy Environment on the Binary Decision Strategy in a Quantum Illumination Radar

A quantum illumination radar uses quantum entanglement to enhance photodetection sensitivity. The entanglement is quickly destroyed by the decoherence in an environment, although the sensitivity enhancement could survive thanks to quantum correlations beyond the entanglement. These quantum correlations are quantified by the quantum discord. Here, we use a toy model with an amplitude damping channel and Lloyd’s binary decision strategy to highlight the possible role of these correlations from the perspective of a quantum radar.


Introduction
Since the late 20th Century, radar technology has been used in many applications, especially for maritime and aeronautic purposes [1][2][3]. One of the most important subjects in radar technology concerns the detection of stealth targets in the context of background noise. In another way, the current development of quantum technologies provides new possibilities for remote detection, leading to the concept of quantum radar. This paper presents a "toy model" for quantum radar based on quantum entanglement between pairs of photons. Such a simple model does not aim to be realistic but rather provides pedagogical value concerning the potentiality created by the quantum radar.
The current development of quantum technologies for the transmission of information introduced the idea of "quantum radar", although this idea remained of little interest until Lloyd's article was published in 2008 [4]. In this article, Seth Lloyd showed that the quantum entanglement with pairs of photons can significantly improve the remote detection sensitivity in the optical frequency regime. This way of using entanglement for remote detection is called "quantum illumination" (QI). Since this article, interest in the field of quantum radar has grown. New theoretical and experimental research has been conducted on this subject [5][6][7][8][9][10][11][12]. Research around the quantum radar has moved from focusing on the individual photons to small bunches of photons [4,11]. In the same sense, research has moved from the optical frequency regime [4] to the microwave frequency regime [11][12][13], which is more suitable for radar applications but also more challenging. In this context, new technologies are currently being developed to make quantum illumination possible in the microwave regime. For instance, we can cite the Josephson junction, which enables the direct production of microwave-entangled photons at low temperature. There is also the coupling between an optical photon and a microwave photon [11]. Then, the Nitrogen-Vacancy centers (called NV centers) also permit the production of microwave entangled photons. Despite the great difficulties relating to the feasibility of such a quantum radar, this research field is highly active.
The quantum radar has the same purpose as the conventional radar, but the functioning relies on the principles of quantum mechanics.
We recall that a conventional radar works with the classical theory of electromagnetism based on Maxwell equations. To briefly summarize it, a radar is a device that sends an electromagnetic wave to detect a reflecting object that reflects a fraction of the incident wave to the radar. The radar scheme is characterized by the energy ratio called the "radar equation" between the received and the emitted wave. The received wave provides information on the detection and the position of an object (i.e., radar target) using signal processing.
In comparison, a quantum radar relies on quantum mechanics to work, but the current definition is not as clear. A quantum radar could be defined as a conventional radar that uses only a quantum electronic device to improve its sensitivity. Hence, if any quantum electronic device is found in the radar system, we could refer to it as a quantum radar. From another perspective, the use of a low number of photons instead of classical electromagnetic waves also means that we can refer to a quantum radar. Using only one of these aspects theoretically means that we can refer to a quantum radar. Given these arguments, the current definition of a quantum radar is still somewhat ambiguous and unclear, but the common thread between these definitions is the requirement of quantum mechanics [11,14].
Among the known types of quantum radar, the quantum illumination radar proposed by Lloyd [4] remains the most interesting because the enhancement of the sensitivity of the radar from the use of quantum entanglement is significant. Below, the quantum radar exclusively refers to the QI radar. Lloyd's quantum illumination is based on the use of pairs of entangled photons, whereby one photon is trapped in the radar system, while the second photon is emitted into space to be reflected by a target. The emitted photon is reflected by the target, and it comes back to the radar to be measured. In the radar, the detection strategy relies on a joint measurement of the pair of photons. This is a striking point that is nevertheless challenging in experimental terms, as specified in the current literature [11,12].
In Lloyd's article [4], the sensitivity enhancement induced by the entanglement could be maintained even if the entanglement is quickly lost during the propagation phase in the optical regime. That means some quantum correlations, i.e., some quantum information, could survive to the decoherence induced by the propagation environment which destroys the initial entanglement. The entanglement phenomenon underlines the presence of strong quantum correlations, but they are quickly lost due to the environment. The sensitivity resilience could be explained by quantum correlations beyond the entanglement [10]. Such correlations can be quantified by the quantum discord [15]. There is currently lack of studies examining the environmental influence (such as the atmosphere) on the quantum correlation evolution in the quantum illumination radar.
In this paper, we work with a toy model with Lloyd's quantum illumination scheme, introducing the damping effect of the propagation environment to follow the evolution of quantum information. The objective consists in studying a simple QI radar from the point of view of the quantum information theory introducing the quantum discord. Therefore, we study the quantum illumination radar for a pair of qubits entangled on thermal energy states since the detection process is photodetection. A generalized amplitude damping quantum channel is used to model the environmental influence acting on the thermal energy levels of one qubit alone. We use Lloyd's decision strategy to introduce a low average number n b of thermal photons by mode, representing the environmental thermal noise. The decision strategy is linked to the quantum channel through a damping probability p from the channel model [16]. For the sake of simplicity, we have limited ourselves to a low number of thermal modes d = 2 corresponding to qubit states widely used in quantum information theory [17,18]. The case d = 2 is a limiting case for small energy states, so our toy model does not aim to represent a realistic case.
This paper is divided into four sections. Section 2 explains the basis of the quantum illumination radar and the need to adopt a quantum formalism. The decision strategy is adapted to account for the environmental influence using an amplitude-damping quantum channel acting like a heat bath. In Section 3, we study the radar under the scope of quantum information theory to follow the evolution of information. Next, in Section 4, we introduce Lloyd's binary decision strategy according to the quantum channel used. We discuss the information evolution in the quantum channel linked to the decision theory. Section 5 provides the conclusion.

Description of a Quantum Illumination Radar
In this section, we introduce Lloyd's QI radar for pairs of entangled photons. We use a conventional radar scheme to explain the differences between the QI radar and a classical radar. The quantum formalism for the QI radar is introduced in Section 2.2. Then, the formalism used to define the amplitude damping channel is given in Section 2.3.

Description of the Radar System
The principle of a QI radar relies on the entanglement and on a joint measurement on the entangled system [4]. We start by considering a conventional radar scheme to highlight the differences between a quantum radar and a classical radar.
The QI radar scheme is built like the conventional radar scheme in a monostatic case. We have four main sections that permit us to define a quantum radar equation. In conventional radar theory, the radar scheme is constituted by the emitter (an emitting antenna), the propagation channel (the atmosphere), the target (a reflecting object), and the receiver (a receiving antenna) [1]. Following the same scheme, the quantum radar scheme is constituted by the emitter, which produces pairs of entangled photons in optical or microwave frequencies [11]. Next, the propagation channel describes the perturbative propagation medium, such as the atmosphere, acting on the emitted photon. The second photon is trapped in the radar, which is a non-perturbative medium. The radar target is an object that reflects the incident photon without absorbing it. Then, the receiver is a quantum sensor taking a joint measurement of the photon received and the trapped photon.
Both radar schemes are used to define the classical radar equation in Equation (1) and the quantum radar equation in Equation (2).
In classical radar theory, electromagnetic waves (EW) are used to detect distant objects. The energy used to produce the EW has to be significant, as the atmosphere induces an attenuation of the incident wave. This attenuation is due to scattering and absorption by the atmospheric molecules (H 2 , O 2 , etc.) [1]. Furthermore, an energy decay occurs because the EW propagates in all directions. Only a fraction of the incident wave is reflected by the target in all directions before it is received by the radar antenna. Therefore, the classical radar equation written in Equation (1) is an energy ratio between the received wave and the emitted wave.
In Equation (1), P t is the transmitted power over a time τ with the gain in the antenna G t . The factor 4πR 2 shows that the incident wave is irradiated in all directions of the space in a sphere of radius R. The factor 1/L 2 represents the wave attenuation due to the absorption and scattering over a round -trip. σ C is the Radar Cross Section (RCS) that indicates how an object reflects an incident wave. The RCS depends on the target geometry and on its nature [1]. A e = G rc λ 2 /(4π) is the effective area of the receiving antenna with a gain G rc , and it depends on the wavelength λ = c/f [3]. Here, the classical radar equation has the form of a signal-to-noise ratio (SNR).
In quantum illumination radar theory, we emit single photons, not EW. In the atmosphere, the emitted photon can be absorbed or scattered, in turn implying the loss of information, i.e., the loss of any useful signal. Hence, we add a probability factor P p for the photon not being absorbed or scattered by the propagation environment; the emitted energy and the received energy are the same. In addition, this photon energy E =hω is far lower than the energy of classical EW. Therefore, the quantum radar equation is not an energy ratio but a probability law to retrieve the emitted photon after its propagation and its reflection on the target. This probability law is written in Equation (2) following the form of Equation (1).
where P p is the probability of the photon not being absorbed and scattered by atmospheric molecules along the propagation. This probability is squared since the photon makes a round-trip between the radar and the target. σ Q is the Quantum Radar Cross Section (QRCS), which is the analog equivalent of the RCS, i.e., the way for an object to reflect the incident photon. Some research has already been conducted on the QRCS between a photon and reflecting targets with canonical geometries [19]. For photons with polarization and a non-polarized target, the interaction does not significantly alter the quantum state of the photon [20]. P a is the probability of detection of the quantum sensor. It depends on the effective area of the sensor and on its sensitivity threshold. For instance, detecting microwave photons is more challenging than detecting optical photons because their energy is much lower.
It should be noted that the QI radar equation is equal to the single-photon radar equation described in Section 4 since, in both radars, only one photon is emitted. The fundamental difference lies in the entanglement for the QI radar between photon A (ancilla) trapped in the radar and photon S (signal) that propagates. In the following, we describe in Figure 1 how a QI radar works using entangled photonic qubits.  Figure 1. The quantum illumination radar process has three steps. In Step 1, the entangled pair of photons is created inside the radar. In Step 2, we perform a separation. In 2a, the blue photon is trapped inside the radar, while in step 2b, the red photon propagates through a medium before being reflected by the target. The red photon returns to the receiver. In step 3, the pair of photons is gathered, and a joint measurement of the quantum state blue-red is taken.
In Figure 1, the QI radar principle for a pair of photonic qubits is described in three steps. These steps condition the quantum formalism used in Sections 2.2 and 2.3 and the decision theory used in Section 4.
Step 1 consists in creating the pair of entangled photonic qubits. Photon A (ancilla) is represented in blue and photon S (signal) is represented in red. Both are qubits, which means that we describe them with a state vector |ϕ i = A, S as a superposition of any two quantum states 0 and 1 : The state vector represents a quantum state ϕ of the quantum system i = A, S in a Hilbert space H i of dimension 2. The state vectors represent pure quantum states with a basis of vectors in H i . However, in practice, we usually prefer the density matrixρ = j,k p j,k |j, k j, k| on the eigenvectors |j, k with the associated probability p j,k . The density matrix notation is suitable to mixed states, and it is widely used in quantum information theory because a density operatorρ i represents all the available information about a quantum state on the system "i" [17] (see calculations in Section 3.1).

In
Step 2, the radar spatially separates the pair of photonic qubits. Photon A is kept in the radar. Photon S is emitted into space towards a reflecting object. In Step 2a, photon A is trapped in the radar using a mirror cavity without any loss. This is an ideal case because ensuring a lossless cavity is very difficult in practice. Another point is to prevent any perturbation in system A, particularly to prevent any measurement. In Step 2b, photon S travels through a propagation environment, and the quantum state of S is perturbed by the interaction with the environment. For instance, in the atmosphere, the perturbations are mainly due to atmospheric molecules (N 2 , O 2 , etc.). The perturbations undergone by photon S have an impact on the system AS since the system is initially entangled. After the propagation, photon S interacts with the reflecting object. We assume the photon is simply reflected by the object before it comes back to the radar. This approach is idealistic since in practice photon S can be absorbed or scattered during the propagation; furthermore, it could be also reflected anywhere by the object. Finally, photon S is caught by the receiver of the radar.
Both photons are gathered in the radar system in Step 3. The QI radar takes a joint measurement on system AS to identify its new quantum state. In reference to Lloyd's quantum illumination [4], the initial entanglement is lost due to the propagation in the environment, and therefore the quantum state is modified.
Step 2b introduces two significant sources of perturbation independent from any technological devices used in the radar. In this step, there is a perturbation due to the propagation environment and a perturbation due to the interaction with a distant object.
The interaction between the photon and the object is calculated with the QRCS. According to the references [19,20], the interaction effect on the quantum state of the photon is negligible because the interaction time is too fast. Therefore, considering the object as a simple reflecting plate, we suppose this interaction effect is negligible.
The propagation environment has a strong effect on the photon state since it induces the decoherence, which in turn destroys the initial entanglement. In Lloyd's article [4], the entanglement decay is immediate in the optical frequency regime.
In our toy model based on Lloyd's quantum illumination, we have decided not to focus on the object influence on the qubit state. Instead, we have focused our attention on the propagation environment and its influence on the qubit state. This influence on the qubit state S affects the state of the system AS due to the quantum correlations. A generalized amplitude damping channel is used in Section 2.3 to model the environmental influence, and the link with the decision theory (Section 4) is made with the damping parameter p used in the quantum channel. The channel permits one to compute the quantum information evolution as a function of p in Section 3.2.
Lloyd's quantum illumination is based on the photodetection of a target of reflectivity η with single-photon emissions in a bath of thermal photons. The thermal photons are assumed to have an average occupation probability n b = (1 − e −βhω )e −βhω over the thermal modes, where β = (k B T ) −1 with k B , the Boltzmann constant, and T is the temperature. According to Lloyd's article, the photodetector can distinguish d = W × t d modes during a detection event where W is the frequency bandwidth and t d is the temporal window for the detection. According to the channel used in Section 2.3, we have limited ourselves to the limiting case of d = 2 modes. The photodetector can detect at most one photon in a detection event that corresponds to dn b 1. This approximation is valid for optical frequency regimes, but this is not the case when we move on to microwave regimes. The binary decision strategy relies on photodetection, which is why we must work on thermal energy modes in the quantum damping channel to model the situation. The quantum channel depends on a damping parameter p, which is theoretically dependent on the environment but p does not represent the absorption/scattering probability P p of Equation (2). The decision strategy uses this parameter p to compose the probabilities of detection error and the signal-to-noise ratios in Section 4.
Thus, we must work with pairs of entangled thermal qubits in the quantum channel. The next section gives the formalism for the quantum states in the QI radar. Section 2.3 uses this formalism to define the amplitude-damping quantum channel.

Expression of the Quantum States in the QI Radar
Firstly, we discuss about the quantum entanglement and its consequences in the equations. Then, we introduce the quantum formalism for the photonic qubits and for the state of two entangled qubits.
The entanglement is a purely quantum phenomenon, and there is no classical counterpart [21]. A quantum system "1,2" in a quantum stateρ 1,2 constituted by two subsystems "1" and "2" in quantum statesρ 1 andρ 2 is entangled if we cannot describe it as a product state of its two subsystemsρ 1,2 =ρ 1 ⊗ρ 2 , where ⊗ is the Kronecker product. A product state is also called a "separable state", i.e., a state with no entanglement. Therefore, an entangled state induces strong quantum correlations between subsystems "1" and "2". This implies that the stateρ 1 of the subsystem "1" is conditioned by the stateρ 2 of the subsystem "2" and vice versa. The entanglement is the characteristic feature of quantum mechanics, and it represents the key resource for several information processes in quantum information theory [17,18]. Nevertheless, we have seen in Section 1 that some quantum correlations exist beyond entanglement. These correlations can be quantified by the quantum discord, which is different from the entanglement. The difference between the entanglement and the quantum discord is pointed out in Section 3.1. Now, we must define the formalism used for thermal qubits in Equation (3) on the thermal energy states ε i=1,2 , since Lloyd's quantum illumination is based on the photodetection.
Equation (3) is the reduced density operator on the subsystem S involved in the entangled stateρ AS, thermal in Equation (4) on the basis of the energy states {|ε 1 ε 1 , |ε 1 ε 2 , |ε 2 ε 1 , |ε 2 ε 2 } of two qubits.ρ This entangled state is designed for d = 2 modes, and it corresponds to the shape of entangled states used in Lloyd's article [4] : In the atmosphere, photon S is submitted to a thermalization by the environment that induces an evolution of its quantum state. An evolution ofρ S, thermal implies an evolution ofρ AS, thermal .
To model such an evolution, we use the article [16], which presents an experiment to study the photon thermalization. This evolution is described by a generalized amplitudedamping channel presented in Section 2.3 that we use to model the evolution ofρ AS, thermal . The evolution ofρ AS, thermal implies the evolution of its quantum information, as we will see in Sections 3 and 4.
The following section uses Equation (4) as the input state of the amplitude-damping channel. Note that for this QI radar, we do not consider the reflectivity of the target (η = 1), since we work only on the information evolution along the quantum channel modeling the environment.

The Amplitude Damping Channel for the QI Radar
The amplitude-damping quantum channel acts like a heat bath on the qubit S, which propagates in the environment, as seen in step 2b in Figure 1. More precisely, the quantum channel acts locally on the qubit S, but it does not act locally on the qubit A, which is kept in a quantum memory. Therefore, photon A trapped in the radar remains unperturbed by the quantum channel. Moreover, in the QI radar, no local measurement is performed on the subsystem A at any time, as mentioned in Section 2.1. The quantum channel makes the quantum state of Equation (4) evolve as a function of the parameter p described below.
The quantum channel used is taken from the reference [16]. We adapt it to the entangled state of Equation (4) acting only on the emitted photon S. We drop the subscript "thermal" to simplify the equations; then, we writeρ AS, thermal ≡ρ AS, in . The amplitude damping channel is defined in Equation (5).
where the Kraus operatorsK i acting locally on the subsystem S are defined in Equation (6) as functions of the Kraus operatorsΓ i of the original quantum channel in [16] : The Kraus operators verify the unitarity condition i=0..3K † iK i =Î AS . In Equation (6), there are two parameters that control the evolution ofρ AS, in . The first parameter ξ ∈ [0, 1/2] is the thermal population of the excited state |ε 2 . It controls the balance of population on qubit states. For the sake of simplicity, we assume photon S is fully mixed, so we take ξ = 1/2 for the computation in Section 3.2 corresponding to our qubit state defined in Equation (3). However, we keep the general symbol ξ in the equations.
The second parameter p = 1 − e −γt ∈ [0, 1] is the damping probability associated with the excited state |ε 2 as a function of the time t, where γ is the damping constant. The parameter γ depends on the environment and the frequency photon S. An optical photon is more impacted by the environment than a microwave photon, but setting a value for γ is not within the scope of this paper. γ is only taken as dependent on the environment. Therefore, the important parameter controlling the evolution of the stateρ AS, in is p ∈ [0, 1].
Using the quantum channel of Equation (5) on the density matrix in Equation (4) gives the output density matrixρ AS, out as a function of p and ξ written in Equation (7).
As explained in Section 1, the parameter p depends on the propagation environment. It is the main parameter that evolves in our toy model to represent the intensity of the action induced by the propagation environment on the quantum system AS. This density matrixρ AS, out is the starting point for the calculus of the quantum information evolution in Section 3.2. Parameter p also permits this quantum information evolution to be linked to the binary decision strategy developed in Section 4. To recap, the parameter has three roles to play. First, it describes the action intensity of the environment on the state described bŷ ρ AS, out . Second, it binds the model to a practical laboratory case [16]. Third, it binds the amplitude damping channel to Lloyd's binary decision strategy.
At this point, the amplitude damping channel action on the quantum state AS has been defined in Equation (7). We will compute the quantum information evolution in the QI radar in the next section.

Quantum Information in the QI Radar
The first section introduces quantum information theory using classical information theory. Using a particular quantum state, we compare the entanglement and the quantum discord to highlight its difference. In Section 3.2, we compute the quantum information evolution in the quantum channel.

Tools of the Quantum Information Theory
Quantum information theory is an extension of classical information theory for quantum mechanics [18]. The comparison between the mutual information definitions in both theories leads us to the definition of quantum discord. Secondly, we underline the difference between the discord and the entanglement by performing calculations for a Werner state.
In classical information theory, the information about a random variable X is defined by Shannon entropy H X (x) = − x∈X p x log(p x ). It quantifies the surprise, i.e., the information available for all the realizations x of probability p x for the random variable X. Considering two random variables X and Y , we define the mutual information I(X : Y ) in Equation (8) in two equivalent formulas, Equation (8a,b). The mutual information I(X : Y ) is the common information between these two variables X and Y .
In Equation (8), is the conditional Shannon entropy of X, with knowledge of the realization y ∈ Y , so it represents the information on X conditioned by the realizations y on Y . Then, we have the joint entropy H(X, is the joint information of the random variables X and Y . Now, we move to the definition of the mutual information in the quantum information theory. In quantum information theory, the analog of the Shannon entropy H(X) is the Von Neumann entropy S(ρ) = −Tr{ρ logρ} for a density matrixρ. S(ρ) quantifies the surprise, i.e., the information about a quantum state, since all information about this quantum state is encapsulated in the density matrixρ. The sum in classical information theory is replaced by the trace operation Tr{·} in quantum information theory. Similarly, the probability distributions are replaced by the density matrices that describe the quantum states as statistical mixtures of the eigenstates in a Hilbert space. As shown in Equation (8), the quantum mutual information can be defined in two ways in Equation (9); however, neither formula, Equation (9a) nor Equation (9b), is equivalent.
is the natural extension of the classical mutual information in Equation (8a), while Equation (9b) looks like the conditional information of Equation (8b). Compared to the classical case in Equation (8b), Equation (9b) introduces measurement theory, which is fundamental in quantum mechanics. Consequently, the latter definition reports how much information we can have on A if we measure S. Equation (9a) is called the "quantum mutual information", while Equation (9b) is called the "classical information" available in the quantum stateρ A conditioned by the measurement of the quantum stateρ S . To complete such measurements, we typically use Von Neumann projection operators {M (i) S . In Equation (9b), the quantum conditional entropy S(ρ A |ρ S ) defined in Equation (10) depends on the measurement operators.
where the probability on the stateρ A by a measurement on the state "i" is The quantum discord is then defined in Equation (11) as the difference between the quantum mutual information I(ρ AS ) in Equation (9a) and the classical information J (ρ AS ) in Equation (9b) that is accessible in subsystem A by means of a Von Neumann measurement taken of subsystem S.
This information quantity depends on the measurement operators {M (i) S } as it depends on J (ρ AS ). Typically, the Von Neumann operators are chosen to verify the minimum required for the quantum discord. Note as that the information J (ρ AS ) is not necessarily symmetrical, and the discord d(ρ AS ) has the same property. The logarithm functions used to calculate the discord on qubit states are in base 2.
After the quantum discord, we have to define the entanglement rate. Quantum mechanics only refers to the entanglement of a system, but quantum information theory enables us to quantify the degree of entanglement [17,21]. There are several ways to compute the entanglement rate that depend, for example, on the dimension of the quantum state. As we consider a pair of entangled qubits, we use the concurrence of Wootters to compute the entanglement rate using the entanglement of formation [22].
The entanglement rate E (ρ AS ) using the concurrence is calculated using a binary entropy , the concurrence of Wootters [23]. C(ρ AS ) is itself an entanglement measure such as C ∈ [0, 1]. If C = 1, the state is fully entangled and the entanglement vanishes for C = 0. We compute the concurrence C(ρ AS ) = max(0, which are the square-root eigenvalues in decreasing order of the spin-flip matrix R =ρ AS (σ y ⊗σ y )ρ * AS (σ y ⊗σ y ), whereσ y is a Pauli matrix. The symbol * denotes the complex conjugate.
At this point, we define the quantum discord and the entanglement rate. In the following, we highlight the difference between these two quantities using a Werner state (see [15] for details).
To show that both rates C(ρ AS ) and d(ρ AS ) are different, we take a Werner statê which is a mixture of a fully entangled state Ψ − Ψ − AS and a separable stateÎ AS /4 (= product state) [15]. The parameter z ∈ [0, 1] controls the combination of the two extreme states, so the Werner state is a partially entangled quantum state. In Figure 2, we plot the concurrence, the entanglement rate and the quantum discord as functions of parameter z. The three quantities decrease as z tends toward zero. We focus our attention on the rate E (ρ AS ) and d(ρ AS ). In Figure 2, the entanglement rate and the quantum discord are maximal for z = 1, which corresponds to a fully entangled state Ψ − Ψ − AS . E (ρ AS ) = 1 means the quantum state is fully entangled, and d(ρ AS ) = 1 means we have reached the maximum number of possible quantum correlations.
The entanglement rate and the discord are zero when z = 0, corresponding to the fully mixed stateÎ AS /4. E (ρ AS ) = 0 means that the state is no longer entangled, and we also have d(ρ AS ) = 0, so all quantum correlations have vanished. Now, looking at the intermediate z ∈]0, 1[, we can see that both rates take different values. The entanglement rate decay is faster than the discord rate as z tends towards zero. For z d = 1/3, we observe the entanglement rate E (ρ w ) = 0, the concurrence C(ρ w ) = 0 contrary to the quantum discord d(ρ w ) > 0. This means the stateρ w for z d = 1/3 is no longer entangled, but it still has quantum correlations. These correlations are not due to the entanglement as the state is separable; i.e., the state is not entangled. Figure 2 shows that entanglement and discord translate to two different kinds of quantum correlations, so there are quantum correlations beyond the entanglement.
It ought to be noted that we can have a quantum state that is not entangled with a non-zero discord, but the contrary is not possible. We cannot have an entangled state with a null discord. It is recalled that an entangled state cannot be written as a product state, called a separable state. Looking at Figure 2, when z d = 1/3, we can see that the state becomes separable, but it has a non-zero discord for z ∈ [0, 1/3]. Then, the entanglement is a quantum phenomenon, while the quantum discord is a quantum correlation measurement. The entanglement rate E (ρ w ) does not depend on projective measurements such as the discord d(ρ w ). For the latter, we have to be particularly attentive when selecting the measurement operators to verify Equation (11).

3DUDPHWHU
,QIRUPDWLRQELWV We have now defined the difference between the entanglement rate and the quantum discord. The next section uses the tools of information theory described here to compute the evolution of the quantum information along the amplitude damping quantum channel defined in Section 2.3.

Evolution of the Quantum Information in the QI Radar
In this section, we compute the entanglement rate and the quantum discord in the amplitude-damping quantum channel of Equation (5) as functions of the damping parameter p. We begin by calculating the entanglement rate before calculating the quantum discord.
The entanglement rate is calculated using the concurrence of Wootters C(ρ AS,out ), defined in Section 3.1. We compute the matrix R =ρ AS,out (σ y ⊗σ y )ρ * AS,out (σ y ⊗σ y ), and we extract its eigenvalues. In Equation (12), we write the square-roots of the eigenvalues √ λ i of the matrix R.
The eigenvalue √ λ 1 has a multiplicity of two, while √ λ 2 and √ λ 3 have a multiplicity of one. For each value of p, these eigenvalues are sorted in decreasing order to compute the concurrence C(ρ AS,out ) and next to evaluate the entanglement rate E (ρ AS,out ). Both rates are displayed in Figure 3, where we observe that the entanglement is lost before the parameter p reaches its limiting value 1. Indeed, we find C(ρ AS,out ) = 0 for p ≈ 0.83, so for p > 0.83, the system AS is no longer entangled. Therefore, the fully entangled state of Equation (7) becomes a separable state before the limiting value p = 1.

3DUDPHWHU ,QIRUPDWLRQELWV
( AS, out ) ( AS, out ) Now, we will evaluate the quantum discord d(ρ AS,out ) by calculating, in this order, the Von Neumann entropies S(ρ A,out ), S(ρ S,out ), S(ρ AS,out ), and S(ρ A,out |ρ S,out ) as functions of p. With these entropies, we compute the quantum mutual information I(ρ AS,out ) and the classical information J (ρ AS,out ) from Equation (9) to estimate the quantum discord using Equation (11).
Using Equation (7), we calculate the density matrices of both subsystems A and S using a partial trace operation :ρ A,out = Tr S {ρ AS,out } andρ S,out = Tr A {ρ AS,out }. Next, we calculate the Von Neumann entropies of both subsystems in Equation (13).
According to Equation (13a,b), the information in subsystems A and S is not affected by the action of the quantum channel. We have exactly the same amount of information for all values of parameter p ∈ [0, 1] for photon A and for photon S. Now, we apply the same calculation for the entropy S(ρ AS,out ) = −Tr{ρ AS,out logρ AS,out } on the entangled system AS. It should be noted that the logarithm used is in base 4, log(x) = ln(x)/ ln(4) as dim(ρ AS,out ) = 4.
The matrix logarithm log(ρ AS ) can be calculated analytically by diagonalizing the density matrixρ AS,out to obtain log(ρ AS,out ) = P log(D) P − 1, where P is the modal matrix and P −1 is its inverse. D is the diagonal matrix formed by the eigenvalues ofρ AS,out . The calculation is detailed in Appendix A within which we wrote out the Von Neumann entropy S(ρ AS,out ) = −H 1 − H 2 − H 5 − H 6 with coefficients {H i=1...6 }, also defined as functions of p in Appendix A. The evolution of S(ρ AS,out ) is plotted in Figure 4. We observe that the Von Neumann increases as p tends towards one as we lose the quantum correlations because of the quantum channel action. For p = 0, S(ρ AS,out ) = 0 corresponding to the fully entangled pure state and next the entropy increases because the uncertainty on the stateρ AS,out increases. Indeed, as p tends towards one, we progress to a fully mixed separable state.
Here, we compute the conditional entropy S(ρ A,out |ρ S,out ) S =Î A ⊗ |ε 2 ε 2 | S . We obtain the conditional entropy in Equation (14) that is plotted in Figure 4.
The conditional entropy increases as p tends towards one, and we observe its increase as the quantum stateρ AS,out loses its entanglement in Figure 3. The entanglement loss means we progressively lose the quantum correlations between the quantum stateρ A,out and the quantum stateρ S,out . Consequently, we increase the uncertainty between both states A and S when a projective measurement is completed forρ S,out .
Using the calculations of the entropies, we can now compute the quantum mutual information I(ρ AS,out ) and the classical information J (ρ AS,out ) as functions of the parameter p using Equation (9) to plot them in Figure 5. The quantum mutual information is maximal for the fully entangled state when p = 0, and it decreases to zero until p = 1. We observe the same evolution for the classical information J (ρ AS,out ). The quantum mutual information I(ρ AS,out ) represents the total information shared by both subsystems A and S. Its decreasing evolution in Figure 5 is in line with the entanglement decay observed in Figure 3. The entanglement decay means the destruction of the initial quantum information shared by both subsystems because the quantum correlations are vanishing. The classical information decay J (ρ AS,out ) shows that we lose the ability to extract information from stateρ A,out when projective measurements are taken in stateρ S,out . A comparison between Figures 3 and 5 shows that when the entanglement vanishes at p ≈ 0.83, the classical information J (ρ AS,out ) is close to zero but not the quantum mutual information I(ρ AS,out ).
Finally, we compute the quantum discord d(ρ AS,out ) as a function of parameter p in Figure 6, including the entanglement rates C(ρ AS,out ) and E (ρ AS,out ). We clearly see that the entanglement rate of decay is faster than the quantum discord decay as p tends towards one.
When p > 0.83, the entanglement rate vanishes but not the quantum discord. This means that some quantum correlations survived to the amplitude damping quantum channel in Equation (5) despite the stateρ AS,out | p>0.83 no longer being entangled, i.e., it is a separable state. When p tends towards one, the quantum discord vanishes as well, so all quantum correlations are lost and we obtain a maximally mixed separable state that is a product stateρ A,out ⊗ρ S,out = 1 4Î A ⊗Î S .

3DUDPHWHU
,QIRUPDWLRQELWV  Through calculations of the entanglement rate and the quantum discord in this section, we observed the presence of quantum correlations surviving slightly after the entanglement loss. In Lloyd's article [4], we keep the quantum advantage in terms of sensitivity even when the entanglement is lost. Reasoning with a quantum information approach allows us to understand that it is the existence of quantum information, i.e., quantum correlations, that produce the sensitivity enhancement. However, keep in mind that quantum correlations will certainly vanish in atmosphere after a finite propagation time, which leads to the question of the quantum radar ranging. Indeed, if entanglement is lost but not the quantum advantage, it is when the quantum correlations are totally lost that our sensitivity enhancement would collapse. Hence, even if the model presented is not realistic, it is worth linking the quantum channel action with the decision strategy presented in Section 4, since the entanglement rate and the quantum discord do not evolve at the same speed as function of the damping parameter p. This comparison between the amplitude damping channel and the binary decision strategy permit us to combine the propagation phase with the detection phase including thermal noise. The next section presents Lloyd's decision strategy, accounting for the quantum channel defined in Section 2.3. Section 4 shows the link between the quantum information evolution in the amplitude-damping quantum channel and Lloyd's binary decision strategy for the QI radar [4].

The Binary Decision Strategy for the QI Radar
In this section, we introduce Lloyd's binary decision strategy for the QI radar, establishing the link with the amplitude damping quantum channel of Section 2.3.
We start by explaining Lloyd's original decision strategy, and we then adapt this strategy according to the quantum channel action to calculate the signal-to-noise ratios (SNR) depending on the damping probability p.
In Lloyd's article [4], the binary decision strategy used relies on the discrimination of quantum states [24]. This discrimination of quantum states is enhanced by the entanglement in a QI radar compared to a single-photon radar. The binary decision can only provide the information on the absence or the presence of an object with a reflectivity 0 η 1 surrounded by a thermal noise. These hypotheses are, respectively, called hypothesis H 0 and hypothesis H 1 . Figure 7 provides an illustration of the situation depicted by both hypotheses. At the top of Figure 7, we have hypothesis H 1 , where a target can be detected. On the bottom, we have hypothesis H 0 , where only thermal noise photons are present. The construction of the approximations for the thermal quantum states associated with each hypothesis is possible for two reasons. The first reason is that the average number of thermal photons n b = (1 − e −βhω )e −βhω per thermal energy mode is very low : n b 1. This average number n b is calculated from Planck's law n(ω, T ) = (e βhω − 1) −1 . This approximation is suitable for the optical frequency regime but not for the microwave frequency regime. The second reason is that the photodetector can distinguish d = W × t d modes per detection event, as asserted in Section 2.1. Moreover, the photodetector can detect one photon at most per detection event, corresponding to a number of thermal photons detected dn b 1.
For Lloyd's QI radar, the detection works on a reflecting target with a reflectivity η, and the assumption was that the environment destroys the entanglement instantly, but the sensitivity enhancement is maintained. Atmosphere (N 2 , O 2   In hypothesis H 0 , only thermal photons can be detected, since no object is present.

Propagation in
In this paper, our toy model considers Lloyd's QI radar, but we only consider the influence of the propagation environment as we do not account for the target. Therefore, we consider the target as a perfectly reflecting object (η = 1) that reflects the incident photon S towards our quantum radar with certainty. Such a simplification is obviously unrealistic since the photon can be reflected anywhere in practice. We do assume that the photon comes back to the radar in the toy model. Thus, strictly speaking, we do not perform a true QI radar detection under such simplifications. Our objective is to study the evolution of quantum information in the damping channel with the signal-to-noise ratios (SNR) calculated with the binary decision strategy.
The assumptions for calculations are the same as those in the previous paragraph, but we limit ourselves to the case d = 2 modes, and hence the number of thermal photons seen by the photodetector is 2n b 1, which is obviously not realistic. Nevertheless, it permits us to work with thermal qubits in the decision theory and to use the amplitude damping channel acting as a heat bath on a qubit state as stated in [16]. The link with the channel is defined by the damping probability p = 1 − e −γt described in Section 2.3. This parameter p ∈ [0, 1] shows that the quantum channel action increases as p tends towards one. Next, we will take Lloyd's binary decision strategy and adapt it to our model.
In light of the reference article [4], we consider a single-photon radar and a quantum illumination radar. We begin by defining the thermal states on the single-photon radar for hypotheses H 0 and H 1 . Next, we use these thermal states to define the states on the pair of entangled qubits following the QI radar depicted in Figure 1 for both hypotheses H 0 and H 1 . In both radar scenarios, we compute the probability of detection error for single-shot measurements.
For the single-photon radar, according to Lloyd [4], the approximation of the thermal state found for hypothesis H 0 is given by Equation (15) where |vac is the vacuum state on the thermal energy modes and |ε i are the thermal modes populated by one photon. The photodetector cannot have a useful signal, since there is no object to reflect the photon. The photodetector can only see the vacuum or a thermal photon over the d = 2 modes.
In hypothesis H 1 , the object is present, and it can reflect the emitted photon. We obtain Equation (16), whereby we have a probability p of having only the thermal state of Equation (15) given the quantum channel influence, or we can retrieve photon S thermalized by the propagation environment.
where |ϕ S = 1/ √ 2(|ε 1 + |ε 2 ) is the qubit state of photon S. From this point, we have to define the probability of detection error for single-shot measurements. We use the reference [24] and Lloyd's article [4]. It consists in completing projective measurements on the positive and the negative parts of the operator (ρ 1 −ρ 0 ) in Equation (17).
The negative part corresponds to hypothesis H 0 , and the positive part corresponds to hypothesis H 1 . Therefore, the probability to detect one particular state is a conditional probability to obtain a positive/negative result, knowing the starting hypothesis H 1 or H 0 . The probabilities of detection error are conditional probabilities and are written in Table 1.
Following the QI radar scheme depicted in Figure 1, we take into account emitted photon S as well as ancilla photon A. We completed the same calculations for the pair of entangled qubits starting with hypothesis H 0 , whereby we obtained only a separable thermal state, since there is no object to detect. This quantum state is written in Equation (18) according to the reference [4].
In Equation (18), we can detect a thermal noise photon over the thermal energy modes |ε i corresponding to a product state such asÎ A ⊗Î S , while the alternative possibility is to detect only the vacuum |vac .
For hypothesis H 1 , we obtain Equation (19), whereby we have the probability p to obtain the thermal state of Equation (18) because of the quantum channel influence, or we can retrieve the entangled state |Ψ Ψ| AS .
In Equation (19), we have a decreasing probability (1 − p) to retrieve a useful signal as p tends towards one. The probabilities of detection error for single-shot measurements are calculated in Equation (20) with the operator (ρ (H 1 ) AS ) similarly to the singlephoton radar.
The positive part represents hypothesis H 1 and the negative part represents hypothesis H 0 . The probabilities of detection error are represented by the conditional probabilities to obtain a negative or a positive result given the starting hypothesis in Table 1. Table 1. Probabilities of error for single-shot measurements for the single-photon radar and the quantum illumination radar.

Radar Scheme Single-Photon Entangled Photons Outcomes no = (−) yes = (+) no = (−) yes = (+)
In Table 1, logically and in line with the observations of Lloyd, we observe that the probabilities of detection error are enhanced by the number d of modes involved in the pair of entangled qubits. Of course, our toy model gives only an enhancement by 1/2, which is the limiting case. Below, we calculate the SNR for the single-photon radar and for the QI radar in each hypothesis starting with the single-photon radar.
For the single-photon radar, the signal-to-noise ratios (SNRs) for hypotheses H 0 and H 1 consist in calculating the ratio of P e (+)/P e (−) to obtain the SNR written in Equation (21).
In Equation (21), all the SNR depend on the damping parameter p except for the SNR H 0 . Note when p = 1 − e −γt ,with a constant γ, tends towards one it corresponds to an infinite propagation, so we are sure to lose the emitted qubit or the quantum correlations in the QI radar. Looking at Equation (21) and recalling that n b 1, we see that the SNR H 0 depends exclusively on the thermal noise n b . We can associate the SNR H 0 with the probability of a false alarm to detect a signal when there is no object. As n b 1, we have SNR H 0 < 1, so the probability of false alarm is low. For SNR H 1 , we have pn b 1, since p ∈ [0, 1] and n b 1, so for a low value of p, we obtain SNR H 1 > 1. However, when p → 1, the SNR H 1 tends towards zero, so the quantum channel action makes the SNR collapse. For the SNR + , we compare the useful signal of hypothesis H 1 to the noise signal of hypothesis H 0 . As n b 1 and p → 1, we obtain a SNR + ≈ (1 − p)/p, so the SNR decreases quickly to zero because of the factor (1 − p). However, for small values of p, the SNR + is greater than one, because the quantum channel action is not strong enough.
To compare the detection efficiency between both radars, we look at the SNR + , since it takes into account both hypotheses. The quantum channel interacts with the emitted photon, producing the SNR decay until it vanishes for p = 1.
In the QI radar, the signal-to-noise ratios in Equation (22) are computed for hypotheses H 0 and H 1 by calculating the ratio of P e (+)/P e (−) as in the single-photon radar.
As shown by Seth Lloyd in their article [4], the SNR calculated from Table 1 benefits from the number of entangled modes. For instance, the SNR H 0 interpreted as the probability of false alarm is lower than in the SNR H 0 of the single-photon radar thanks to the factor of 1/2 originating from the entangled modes. For the SNR H 1 , as we have n b /2 1, SNR H 1 ≈ (1 − p)/p so the SNR H 1 is greater than one for small values of p. However, it vanishes as p tends towards one because of the thermalization induced by the quantum channel. For the SNR + , as n b /2 1, we also obtain SNR + ≈ (1 − p)/p, so the SNR is greater than one for low values of p but vanishes as p tends towards one. The SNR + shows the comparison between hypotheses H 0 and H 1 , and as for the single-photon radar, we lose the SNR + for p = 1. We plotted the SNR + of QI radar in Figure 8, which is linear with the damping probability p.

3DUDPHWHU
SNR QI SNR QI (+) Figure 8. The SNR + in the QI radar as a function of the parameter p ∈ [0, 1] with n b calculated for f = 5 × 10 14 Hz. The SNR has been normalized by its maximum because as the noise n b 1, the SNR is huge for small values of p.
The SNR H 1 and the SNR + in the single-photon radar and in the QI radar have the same behavior. The SNR is greater than one when p is low, but it finally vanishes when p tends toward one. The link between both radar situations is the amplitude damping channel acting on the emitted qubit, i.e., the qubit S in Figure 1. The quantum channel thermalizes the qubit S, which propagates, and the longer the propagation is, the greater the perturbation on the qubit state is. More precisely, the closer to one the parameter p is, the longer the propagation is. In both radars, the quantum channel perturbs the emitted photon, but, the most interesting observation is of the effect on the QI radar, as for p close to one, we have a greater probability to obtain a separable thermal stateÎ A /2 ⊗ρ 0 , as written in Equation (19).
We recall the different results in the damping channel model and in the binary decision strategy to clarify the link between both models.
The generalized amplitude-damping channel used from the article [16] permits us to describe the evolution of the entangled quantum stateρ AS,out . This channel models the decoherence of the maximally entangled quantum state as a function of the damping parameter p = 1 − e −γt , where γ depends on the atmosphere. Although the value γ is qualitative, we can safely say that the decoherence is greater in the optical domain compared to the microwave domain. In this paper, we stay in the optical domain to verify the assumptions in the binary decision strategy made in Section 4. At the input of the channel, we use the fully entangled state Ψ + Ψ + AS , and we obtain at the output a fully mixed separable stateÎ A /2 ⊗Î S /2. With the channel, we have the environmental action.
In Lloyd's binary decision strategy, we consider only the amount of noise for the photodetection to obtain the SNR according to the hypotheses H 0 and H 1 . In Equation (22), we have the SNR + , which compares hypothesis H 1 with hypothesis H 0 . This comparison between the hypotheses gives us the ratio between the useful signal and the noise from the environment. Looking at Figure 8 shows that the evolution of SNR + is linear. We have the maximum SNR + = 1 for p = 0, while the SNR + vanishes for p = 1. Please note that the calculated SNR + has been normalized because the thermal noise in the optical domain is very low. Furthermore, in Equation (20), the maximum of SNR + corresponds to the fully entangled state |Ψ Ψ| AS , while SNR = 0 corresponds to the fully mixed statê Lloyd's binary decision strategy has been elaborated with the damping probability p used in the channel (Equation (5)). Therefore, both models are linked by this quantity, which allows us to compare the states described byρ AS,out and |Ψ Ψ| AS .
Recall that an open quantum system is a system traveling into a perturbative environment that produces its evolution towards a fully mixed state [25]. Looking at Figure 6, we see that the stateρ AS,out is no longer entangled in the interval p ∈ [0.83, 1], but the quantum discord is different from zero. Hence, on this interval, the quantum stateρ AS,out is separable, but it still has quantum correlations. Now, as shown in Figure 8, the signal-to-noise ratio does not vanish: we have the 0 < SNR + < 1 on the interval p ∈ [0.83, 1]. The SNR + is not zero when the state is no longer entangled, but it still has some quantum correlations. This shows the role of quantum discord in the binary decision strategy for the QI radar.
The damping channel and the decision strategy are coupled by the parameter p, so the extremum states coincide. However, we do not have a perfect correspondence for the separable state, since in the damping channel (Equation (5)), we obtain a fully mixed stateÎ AS /4, and we obtain the fully mixed stateÎ A /2 ⊗ n b 2 i=1 |ε i ε i | for the decision strategy. The only difference comes from the noise n b introduced in the binary decision strategy, which does not appear explicitly in the amplitude-damping channel. In spite of this difference, the separable states obtained are similar.
Finally, when the discord becomes zero, this means that the quantum correlations vanish. We have the maximally mixed separable state in both models. The SNR + = 0 for this state because we cannot make a proper discrimination between the thermal qubit sent in the atmosphere and one thermal qubit from the thermal noise. We can see the vanishing discord as the operational limit scale of the QI radar.
Specifically, the parallel between the damping channel of Section 3 and the decision theory in Section 4 allows us to point out the role of discord in the noise resilience in the QI radar. This follows from the statement of the article [10]. Moreover, we used a link between both models that depends on the environment. This allows us to describe the environmental action on one hand and the signal-to-noise ratios depending on the thermal noise on another hand. In our toy model, we found that the discord is linked to the SNR calculated from the binary decision strategy. Consequently, it may be possible that the quantum discord is responsible of the quantum radar range from a quantum information point of view. Nevertheless, the presented toy model is unrealistic due to a lot of approximations, so we must discuss the limitations of such a model.
Firstly, we used a damping amplitude channel to model the propagation phase where we could use a master Equation [25]. However, the latter requires a lot of calculation based on an interaction Hamiltonian between the photon and the environment. This implies a good understanding of all interaction processes to define a suitable Hamiltonian for the entangled state. The damping amplitude channel is a sufficient tool of the quantum information theory to model the propagation phase because we consider the decoherence in an average process without detailing each phenomenon. The quantum channel approach gives an idea of the global evolution of the entangled stateρ AS,out in the atmosphere. A phenomenological approach would be better, but it is not in the scope of this paper.
Secondly, Lloyd's binary decision strategy relies on approximations on thermal states in the small thermal noise n b , so we are limited to the optical frequency regime for the QI radar. In addition, the link symbolized by p = 1 − e −γt between the amplitude damping quantum channel and the decision strategy is very simple. Moreover, it would be difficult to set a realistic value to γ to improve the toy model. We do not have an exact correspondence between the separable product states obtained with the amplitude damping channel and with the binary decision strategy. However, the common thread between these models is that both states are product states over two maximally mixed qubit systems; hence, there are no quantum correlations remaining. The correspondence is not as good as expected because of the approximations on thermal states and because the quantum channel used is quite simple. It acts like a heat bath on the emitted qubit without introducing an average thermal noise while the thermal noise is introduced by Lloyd's theory. Thirdly, the model is realized for qubit states on the thermal energy levels, so for a photodetector that can see d = 2 modes for an average thermal noise, n b 1. This allows us to link the model to the article [16], but this is an extremely limited case as the energy is very low due to the low number of thermal modes. Consequently, the model presented is not realistic, and in order to have a better approach, we should take a larger number of modes d. However, such a change also requires the modification of the amplitude damping channel used in Section 2.3. This is not as simple because the initial entangled state is more complex to define as for the quantum channel action on this state.

Conclusions
In this paper, we presented a toy model for a quantum illumination radar, focusing on the propagation channel that has a significant impact on radar detection. The QI radar uses pairs of entangled photons to perform a detection. We detailed the technical tools for a QI radar detection, and we introduced the quantum information theory to follow the evolution of quantum information. To highlight the environmental effect, we made a link between Lloyd's binary decision theory and an amplitude damping channel modeling a propagation phase through a perturbative environment such as the atmosphere. Binary decision theory relies on the discrimination of quantum states corresponding to detection hypotheses H 0 and H 1 . This decision theory has been adapted to the quantum channel used to perturb the quantum state of emitted photonic qubit S. The link between both models is performed using the damping parameter p acting on the thermal energy modes of the qubit propagating through the environment. The signal-to-noise ratios are calculated with the binary decision theory as functions of this damping probability p.
The calculations completed for the toy model suggest a link between the discord and sensitivity enhancement resilience in QI radar, as stated in [10]. In the amplitude damping channel based on the full thermalization of a photonic qubit, discord survived longer than entanglement. In the decision theory, the SNR + vanishes only when p = 1 corresponds to a vanishing discord. The disappearance of discord means the disappearance of quantum correlations inside the quantum system used for remote detection. The striking point in Lloyd's theory was the sensitivity enhancement resilience despite entanglement lost. This resilience can be explained by discord because in both models, the limiting quantum state is a product state of the formÎ AS /d, where d is the dimension of the quantum state. A product state is not entangled but, it could have quantum correlations in certain cases. Then, discord should partially explain the sensitivity resilience to decoherence produced by the propagation environment.
The toy model used is not a realistic model since it relies on many approximations, but it illustrates the role of quantum discord. In the quantum channel, we assumed the thermalization is fast and the only parameter is the damping probability p. In decision theory, we assumed a low thermal noise to make approximations for quantum states. Such an approximation on thermal noise is suitable for an optical frequency regime, but it does not work for a microwave frequency regime. In addition, we restricted ourselves to the limiting case of thermal qubits, which means we work with very low energy. To improve this model, we should work with photonic qudits, but that implies modifying the quantum channel used. Then, we have only considered single-photon emissions, which not very well adapted for radar detection in practice.
The results suggest further developing discord understanding for a QI radar. The current work takes part of a broader approach to the QI radar. It would be interesting to consider partially entangled states instead of fully entangled states. Moreover, a proper quantum radar would work with more modes and more photons to probe an object. To perform such a study, it becomes necessary to adopt the continuous variable formalism. Then, a quantum radar would work with microwave frequencies instead of optical frequencies, but it is theoretically and experimentally challenging because of the amount of noise in microwave frequencies. A possible way to overcome this difficulty is to increase the number of modes or the number of photons. At this stage of research, we have several ideas to explore. This paper showed that the quantum discord explains the noise resilience of quantum radar, but we can think about a way to use it as for entanglement. The number of studies on quantum discord in quantum information process increased several years ago, as has the number of experiments on quantum illumination. We are at the very beginning of the development of quantum radar theory. This paper also provides way to spread the idea of the quantum radar potentiality among the radar community. However, please note that we cannot prove with certainty that the QI radar will surpass conventional radar. We can at least think of using both to improve radar detection in particular situations. The eigenvalues found as functions of ξ ∈ [0, 1/2] and p ∈ [0, 1] are written out in Equation (A2) where ∆(ρ AS,out ) is defined in Equation (A3).