Symmetry in the open-system dynamics of quantum correlations

We study the symmetry properties in the dynamics of quantum correlations for two-qubit systems in one-sided noisy channels, with respect to a switch in the location of noise from one qubit to the other. We consider four different channel types, namely depolarizing, amplitude damping, bit-flip, and bit-phase-flip channel, and identify the classes of initial states leading to symmetric decay of entanglement, non-locality and discord. Our results show that the symmetric decay of quantum correlations is not directly linked to the presence or absence of symmetry in the initial state, while it does depend on the type of correlation considered as well as on the type of noise. We prove that asymmetric decay can be used to infer, in certain cases, characteristic properties of the channel. We also show that the location of noise may lead to dramatic changes in the persistence of phenomena such as entanglement sudden death and time-invariant discord.


I. INTRODUCTION
Correlations of genuine quantum nature among the individual constituents of composite systems play a fundamental role in quantum physics.Entanglement is the paramount example of quantum correlations, considered by the founding fathers of quantum physics as the most bizarre aspect of this theory [1].Initially viewed as a mere philosophical subject, it gained popularity and importance with the development of quantum information theory, when it was recognized as a resource for several tasks such as teleportation, superdense coding, and quantum key distribution [2,3].
A specific type of quantum correlation possessed by some entangled states is associated to the concept of quantum non-locality, which implies that predictions of quantum mechanics cannot be simulated by a local hidden variable model [4][5][6].The presence of non-local correlations in bipartite quantum systems leads to violation of Bell-type inequalities, such as the Clauser-Horne-Shimony-Holt (CHSH) inequality [7].States violating Bell inequalities are crucial for certain quantum technologies such as secure quantum communication [8].
For mixed states, there exists a broader type of quantum correlations that does not occur in classical systems, namely quantum discord.It has been demonstrated that discordant states can perform more efficiently than their classical counterparts in certain applications [9].As a consequence, numerous different measures of discord have been introduced in the recent literature to characterize quantum correlations more general than entanglement [10][11][12][13][14][15][16].* hesaly@utu.fiLike all crucially quantum properties, entanglement, discord and non-locality are fragile and generally quickly disappear in presence of noise induced by the environment.Discord is clearly more robust than entanglement to the effects of noise, while entanglement may disappear after a finite time (sudden death of entanglement), discord decays asymptotically [17].Moreover, for certain types of local noise discord may remain constant in time, although the state of the system evolves (time-invariant discord) [18,19].
In this paper we investigate the dynamics of quantum correlations, namely discord, entanglement and nonlocality, for two-qubit systems subjected to various types of one-sided noisy channels.More specifically we focus on a class of states known as X states for which all quantum correlations can be calculated analytically, and we study their symmetry properties with respect to a change in the location of the noise from one qubit to the other.Asymmetric decay of entanglement was earlier studied briefly in [20].Changing the noise location corresponds, for example, to a situation in which Alice produces a state possessing certain quantum correlations and sends it to Bob, who measures the received state.The two parts of the quantum correlated state travel along different paths of equal length, named U (upper) and L (lower).An eavesdropper Eve may attack either one path or the other, her attack being modeled by introducing noise of different type (depolarizing, amplitude damping, bit-flip, and bit-phase-flip channel).
We are interested in understanding whether the dynamics of various quantum correlations is sensitive, and if so how much, to the location of Eve's attack, more precisely to whether she eavesdrops along the U or L path.Of course, we expect that the answer to this question will depend on both the initial quantum correlated state that Alice prepares, and the type of noise introduced by Eve.More precisely, we are interested in identifying the classes of initial states leading to symmetric behavior, with respect to noise acting on either U or L, for the different types of channels and for different types of quantum correlations.Moreover, we investigate how such classes of initial states change for different types of quantum correlations and if there are overlaps between the classes of initial states that lead to symmetric behaviour of quantum entanglement, non-locality and discord.
More specifically, our study gives an answer to the following interesting questions: 1) Is there any connection between symmetry properties of the initial state and symmetry properties of (i) the state dynamics and (ii) the dynamics of quantum correlations (discord, entanglement, non-locality) with respect to a switch in the noise between the channel?In other words, are certain symmetries in the initial state necessary to guarantee a symmetric dynamics of quantum correlations?
2) How sensitive are effects such as entanglement sudden death and time-invariant discord to the location of noise (along the U or L path)?Stated another way, are there situations for which entanglement sudden death or time-invariant discord occurs only when Eve attacks along the U (L) path but not when she attacks the L(U ) path, all other conditions being the same?
3) Assuming that we do not know which type of noise acts on either U or L, can Bob use the symmetric/asymmetric decay properties of quantum correlations to infer or characterize the type of noise, under minimal assumptions?
Our paper is structured as follows.In Sec.II we introduce the class of initial states considered and discuss the symmetry properties of the dynamics and of entanglement as measured by concurrence.In Sec.III we analyze some physical consequences of our findings such as the effect of noise location on entanglement sudden death, the use of asymmetric decay for channel discrimination, and the connection between entanglement decay and entropy.In Sec.IV we discuss the symmetry properties of non-locality and compare the classes of initial states leading to symmetric or asymmetric decay for entanglement and non-locality, respectively.In Section V we discuss the symmetry properties of a distance-based measure of quantum discord and the effect of the noise location on time-invariant discord.Finally in Sec.VI we summarize and present conclusions.
X states are a subclass of 2-qubit states appearing naturally in physical processes [21][22][23][24].We denote an arbi- trary X state by X, where The first condition is equivalent to tr where σ 1 , σ 2 , and σ 3 are the Pauli matrices.Here superscripts U and L refer to operators on Hilbert spaces H U and H L of qubits U and L, respectively.In the vector notation of basis B, the first qubit corresponds to U and the latter one to L. The conditions in Eq. ( 2) are satisfied for instance by all Bell-diagonal and Werner states, which form two important subclasses of X states.
Since entanglement is an essential feature of quantum mechanics, quantifying it has been an active field of research.Multiple entanglement measures have been defined, the most popular of which is concurrence [25].It has been shown [26] that, for all X states, concurrence can be obtained directly from the matrix elements as: Physical dynamics of a quantum state ρ is represented by completely positive and trace preserving linear maps called channels.A map Φ is a channel if and only if there exists such set of operators for all states ρ [27].The decomposition in Eq. ( 4) is called the Kraus decomposition and the operators K i are called the Kraus operators.
In this study we concentrate on one-sided channels, which means that the channel influences only one of the local qubit states at a time, see Fig. 1.This means that the Kraus operators acting on H U ⊗ H L are of the form K i := K i ⊗ I L , when the channel affects qubit U and of the form K i := I U ⊗K i when the channel affects qubit L.Here I U and I L are the identity operators of H U and H L , respectively.We denote the one-sided channels affecting qubits U and L by respectively.We say that the dynamics of a state ρ is symmetric if Φ U (ρ) = Φ L (ρ).Next we study the conditions leading to symmetric and asymmetric state dynamics and entanglement decay under the effects of different channels1 .

A. Depolarizing channel
Corresponding to Eq. ( 5) and ( 6), the dynamics of an arbitrary X state under one-sided depolarizing channels can be written as where p ∈ [0, 1] is the channel strength parameter, telling how strongly the channel influences states and we have denoted r jk := ρ jj − ρ kk , q := 1 − p.By comparing Eq. ( 7) and ( 8 A straightforward calculation shows that entanglement decays symmetrically, i.e.C(Φ U p (X)) = C(Φ L p (X)) ∀p ∈ [0, 1], if and only if ρ 33 = ρ 22 or ρ 11 = ρ 44 .The first condition is necessary for the symmetry of the initial state, as formulated in Eq. (2).Instead the second one, ρ 11 = ρ 44 , is not.
We conclude that symmetry of entanglement decay requires neither swap symmetry of the initial state nor symmetry in the dynamics of the state.

B. Amplitude damping channel
The dynamics of an arbitrary X state under one-sided amplitude damping channels can be written as Comparison of Eq. ( 11) and (12) shows, that the dynamics of a state is symmetric if and only if X = |11 11|, which is invariant.Since for this state ρ 22 = ρ 33 = ρ 23 = 0, symmetric state dynamics occurs only for a single state which is swap symmetric, unlike in the case of depolarizing channel.By (3), we get the concurrences of the output states as Now entanglement decay is symmetric if and only if ρ 33 = ρ 22 or ρ 44 = 0.The first condition is necessary for the symmetry of the initial state but the second one, ρ 44 = 0, is not related to it.
So again, an initially asymmetric state can lead to symmetric decay of entanglement.Also, initial states leading to asymmetric state dynamics can have symmetric decay for entanglement.

C. Entanglement decay in channel combinations
To generalize the analysis, we study the initial conditions for entanglement decay also in channel combinations.By combining local amplitude damping and depolarizing channels with equal channel strengths p, we can create new channels for the two-qubit system.In Fig. 2 we illustrate the most simple two channel combinations.We have solved the families of initial states leading to symmetric and asymmetric entanglement decay in each noise configuration (a)-(d) with the same analysis as used in Sec.II A-II B.
Trivially, entanglement decay is symmetric for all initial states in configuration (a) if Φ = Ψ.If in configuration (b) the channels are chosen as Φ = Ψ, the conditions for symmetric entanglement decay are the same as in the situation of channel Φ influencing just one of the qubits once.The same applies also for configurations (c) and (d): whenever Φ = Ψ = Ξ, the conditions for symmetric entanglement decay are the same as in the case of Φ affecting just one of the qubits once.This means, that in the sense of entanglement decay symmetry, adding identical copies of the same channel does not break the symmetry or create it.
On the other hand, if Φ = Ψ in (a), entanglement decay is symmetric if and only if ρ 22 = ρ 33 .Also, if two of the channels in (c) are different, entanglement decays symmetrically if and only if ρ 22 = ρ 33 .If in (b) and (d) two of the channels are different, entanglement decays symmetrically if and only if ρ 22 = ρ 33 or ρ 11 = ρ 44 = 0. First of these is just one of the swap symmetry conditions, leading trivially to symmetric entanglement decay but the second one is actually a condition which satisfies the non-trivial symmetry conditions of both depolarizing and amplitude damping channels, presented in Sec.II A-II B.

D. Bit-flip and bit-phase-flip channels
To avoid redundancy, we present here only the results for bit-flip and bit-phase-flip channels.For further details we refer the reader to appendices A-B.
We see that, for bit-flip channel, the dynamics of a state is symmetric if and only if ρ 11 = ρ 44 , ρ 22 = ρ 33 and ρ 14 , ρ 23 ∈ R. So, as in the case of amplitude damping channel, symmetric initial state is necessary but not sufficient condition for the symmetry of state dynamics.Due to the form of dynamics of coherences, it is not simple to solve analytically conditions for symmetric entanglement decay for the whole family of X states.By using a restrictive assumption, ρ 23 ∈ R or ρ 14 ∈ R, we can perform the analysis.For this subfamily of X states, entanglement decay is symmetric if and only if ρ 22 = ρ 33 or ρ 11 = ρ 44 .
Interestingly, the families of initial states leading to symmetric state dynamics in bit-flip and bit-phase-flip channels are identical.Also for bit-phase-flip channel it is difficult to solve analytically, when entanglement decay is symmetric for the whole family of X states.If we set the same restriction as used for bit-flip channel above, we see that the necessary and sufficient conditions for symmetric entanglement decay in bit-phase-flip channel are the same as for bit-flip channel.On the other hand, by setting ρ 23 , ρ 14 ∈ C\R in numerical tests, we could not find any initial states leading to symmetric entanglement decay for either of the channels.This serves as evidence for a claim that ρ 23 ∈ R or ρ 14 ∈ R is a necessary condition for an X state to have symmetric entanglement decay in these channels.
To conclude, as in Sec.II A-II B, also for bit-flip and bit-phase-flip channels an asymmetric initial state with asymmetric state dynamics can lead to symmetric entanglement decay.

III. OBSERVATIONS
In Sec.II A-II D we studied families of initial states leading to symmetric and asymmetric decay of entanglement under different local channels.Next we concentrate on analyzing the implications of the results by using depolarizing and amplitude damping channels as examples.For plotting purposes, we set the coherence terms in each initial state to be maximal:

A. Asymmetry in sudden death of entanglement
In Fig. 3 we present the concurrence of a state under amplitude damping channel.The plot shows that, for this choice of initial state, entanglement decay is sudden when amplitude damping noise affects qubit L, whereas entanglement decays asymptotically when the noise is acting on qubit U , instead.This means that, for this particular choice of initial state, the influence of the noise is significantly more harmful when qubit L is affected by the noise.

B. Entanglement decay as a resource
In this section we show that the asymmetry of entanglement decay can be used to gain information on one-sided channels.A similar protocol was introduced in Ref. [28], where behavior of quantum discord and negativity were used to discriminate between channels.
In Fig. 4(a)-(b) we present the concurrence in depolarizing channel for two initial states.By choosing ρ 44 = ρ 11 the noise is guaranteed to have identical influence independent of the location of the noise.The difference between the initial states used for Fig. 4(a) and 4(b) is that the values of ρ 22 and ρ 33 were swapped.The form of Eq. ( 9) and (10) shows that, decay of entanglement is invariant under swapping ρ 22 and ρ 33 , when ρ 44 = ρ 11 .
In Fig. 4(c) and 4(d) we present the plots of concurrence in amplitude damping channel for initial states used in Fig. 4(a) and 4(b), respecticely.We see that ), for both initial states X, when- ever p ∈ (0, 1).We notice that in Fig. 4(c) noise on qubit U has more harmful effect on the concurrence and in 4(d) the noise on qubit L is more harmful.The swapping of the majorization of concurrence in Fig. 4(c) and 4(d) can be seen directly from Eq. ( 13) and (14).In fact, swapping ρ 22 and ρ 33 just swaps the curves corresponding to C(Φ U p (X)) and C(Φ L p (X)).Plots (a)-(d) in Fig. 4 show that, by using these two initial states, we can gain information of the channel which affects the system if we can assume that the noise is either depolarizing or amplitude damping type.If we let the noise influence one of the qubits and make a tomographic measurement2 , there are two possibilities, when using the two initial states presented above: either the concurrence has the same value for both initial states or one initial state leads to higher value of concurrence than the other.In the first case we know for sure that the noise was caused by a depolarizing channel and in the latter case the noise must have been amplitude damping.This means that this pair of initial states can be used to distinguish the two channels.
On the other hand, this can be done without knowing the value of p. So, after determining which channel affected the state, we can also obtain the value of p for each channel by comparing the experimentally determined value of concurrence to the analytical solutions.Note that in the reasoning above we have not assumed anything about the location of the noise either.If we conclude that the noise was amplitude damping, we can compare the measured values of the concurrence for the two initial states.If the measured value of concurrence was smaller for the choice ρ 22 = 0, then we know that the noise was affecting qubit U .On the other hand if the value of the concurrence is smaller for the choice ρ 33 = 0, we know that the noise was influencing qubit L, instead.Same reasoning can be done also for the case when the noise is depolarizing type.Plots and the corresponding pair of initial states are presented in Fig. 5.The protocol is not restricted just to these two channels.By recalling the symmetry conditions of entanglement decay presented in Sec.II C-II D, one notices that we can now distinguish the set of channel combinations (a) and (c), the set of configurations (b) and (d), and individual amplitude damping, bit-flip and depolarizing channels.

C. Entanglement decay and entropy
In Ref. [20] Życzkowski & al. studied asymmetric entanglement decay from the point of view of classical and quantum subsystems.They characterized a subsystem as classical if its von Neumann entropy is smaller than the von Neumann entropy of the total system, and as quantum if it is not classical.It was shown, through an example state, that when one of the subsystems is classical and the other one is quantum, noise affecting the classical subsystem decreases entanglement faster than if it was influencing the quantum subsystem.
In Fig. 6 we present concurrence for another initial state as a function of channel strength parameter p in amplitude damping and depolarizing channels.Von Neumann entropy of the initial total system state is S(X) ≈ 0.40 and Von Neumann entropies of the re- Here the red and blue curves correspond to the cases with local noise on qubit U and L, respectively.We have chosen ρ11 = 0.9, ρ22 = 0, ρ33 = 0.08, and ρ44 = 0.02.
duced states of qubits U and L are S(tr L [X]) ≈ 0.14 and S(tr U [X]) ≈ 0.47, respectively.Now subsystem L is quantum but subsystem U is classical.By looking at Fig. 6 it is clear that the result of [20] does not hold generally, since in Fig. 6(a) entanglement decays faster when the noise affects qubit L but in Fig. 6(b) the situation is the opposite.This also shows that, whether local noise on qubit U results to faster or slower entanglement decay is not only a property of the initial state, but depends also on the channel.

IV. DECAY OF NON-LOCALITY
In their famous paper [29], Einstein, Podolsky and Rosen concluded, that some, hidden, variables should be added to quantum mechanics to restore locality to the theory.By assuming that the hidden variable theory is of local realistic nature, Clauser, Horne, Shimony and Holt derived the so-called CHSH inequality, which can be used to test the local hidden-variable theories [7].
The CHSH inequality can be written for a system in state ρ as F ≤ 2, where is the Bell function, â and â are some variables with values ±1 for qubit U and b and b are some variables with values ±1 for qubit L. Whenever F > 2, the locality assumption is violated, and we say that the state is non-local.CHSH inequality has been violated in experiments repeatedly, proving that the local hidden variable theories cannot be valid [30].
In [31] it was shown, that for 2-qubit states F = 2 √ u + ũ, where u and ũ are the two largest eigenvalues of U ρ = M T ρ M ρ , and M is a matrix defined by M i,j = tr[ρσ i ⊗ σ j ], where i, j ∈ {1, 2, 3}.For X states the eigenvalues of U ρ become Finally, we get F j = 2 max{F j 1 , F j 2 }, where F j 1 = u j 1 + u j 2 and F j 2 = u j 1 + u j 3 , and the superscript j ∈ {U, L} tells whether the channel influences qubit U or L. Next we study the behavior of F under local depolarizing and amplitude damping channels and their combinations.As in Sec.II, also here we are interested in whether the dynamics of F depends on the location of the noise.

A. Depolazing channel
For one-sided depolarizing channel, we see that This means that the effect of depolarizing noise on nonlocality of the state is independent of the location for all X states.Since certain initial states lead to asymmetric entanglement decay in depolarizing channel, this implies that non-locality and entanglement behave in different way in terms of location of the noise.

B. Amplitude damping channel
For one-sided amplitude damping channel, we see that Comparison of Eq. ( 17) and (19) shows that now the location of the noise makes a difference, unlike in the case of depolarizing channel.We see that The first condition is satisfied by all symmetric initial states, so, as in the case of concurrence, also symmetric decay of non-locality seems to be a direct consequence of symmetric initial state.On the other hand, the second condition does not require symmetric initial state.For all initial X states satisfying the second condition, the value of F is maximized with the choice  (a) Depolarizing channel.
state we get F = 2, which is not interesting in the context of decay of non-locality, since such state is initially local.
In Fig. 7 we present the behavior of F in amplitude damping channel.The plot illustrates the difference between the two noise locations: when the channel affects qubit U , the decay is twice as fast compared to the case of noise on qubit L. For comparison, we present also the plot of concurrence for the same initial state.In contrast to what happens to the Bell function, concurrence decays independently of the location of the noise.
We have gethered in Tab.I the families of initial states leading to interesting dynamics for concurrence and Bell function.The most interesting result is obtained with choices ρ 22 = ρ 33 , ρ 11 = ρ 44 = 0.For this family, depolarizing channel leads to asymmetric entanglement decay and symmetric decay of Bell function.Contrary to this, the result for amplitude damping channel is the opposite: symmetric entanglement decay and asymmetric decay of Bell function.This means that there is no hierarchy between the asymmetry of entanglement decay and decay of Bell function: asymmetry of one property does not imply or exclude the asymmetry of the other.

C. Channel combinations
To complete the study of Bell function decay under depolarizing and amplitude damping channels, we perform the analysis on combinations presented in Sec.II C. Trivially again F decays symmetrically in (a) if both channels are the same.Then again, if they are different, the decay is symmetric if and only if ρ 22 = ρ 33 .In fact, this is the same condition as for symmetric entanglement decay in this configuration.If in (b), (c), or (d) all channels are depolarizing (amplitude damping) type, the symmetry conditions for F decay are the same as for single depolarizing (amplitude damping) channel.So, as in the case of entanglement decay, also the symmetry of Bell function decay seems to be invariant under repetition of the same local noise.
On the other hand, if in combinations (b) -(d) there is at least one copy of each channel, F decays symmetrically if and only if ρ 22 = ρ 33 .We note, that this differs from the conditions of symmetric entanglement decay in this configuration.The only exception appears in configuration (c).If Φ and Ξ are amplitude damping channels and Ψ is depolarizing channel, F decays symmetrically for all initial X states.For the sake of example, let us assume a situation, in which symmetric decay of non-locality is desired.Configuration (a) can be divided into two cases in terms of symmetry conditions: Φ = Ψ, leading always to symmetric decay of F, and Φ = Ψ leading to symmetric decay of F if and only if ρ 22 = ρ 33 .In the latter case we can now achieve symmetry for all initial X states, by adding amplitude damping noise before depolarizing channel.It is worth noting, that there is something special about this configuration, since it is impossible to induce the symmetry by adding amplitude damping noise after the depolarizing channel or on the same side with the original amplitude damping channel.Also, if the original configuration has just a single amplitude damping channel, one can achieve symmetry by adding local amplitude damping and depolarizing noises on the qubits.
The analysis in Sec.II C shows, that such phenomenon does not occur for entanglement decay in simple combinations of depolarizing and amplitude damping channels: in the case described above, adding one amplitude damping or depolarizing channel in any possible location has no effect in the symmetry of entanglement decay.

V. DECAY OF TRACE DISTANCE DISCORD
Since all quantum correlations cannot be described by entanglement and non-locality, we conclude our study by considering the dynamics of a more general type of correlation, quantum discord.Due to the difficulty of computing and comparing the exact values of quantum discord, geometric measures have been developed.Geometric discord measures are based on the smallest distance between the given state ρ and the set of states with zero discord.A state ρ has zero discord if and only if it can be decomposed as where {|α j } j is a set of orthogonal vectors in H U and ρ L (j) are positive operators in H L .The choice of metric used to measure the distance determines the properties of geometric discord.A good choice for metric is trace distance D tr (ρ, ξ) = ||ρ−ξ|| tr , where ||ρ|| tr = tr[ ρ † ρ]/2 is the trace norm.With this norm, the trace distance discord can be defined as where the minimization is taken over the set all states ρ with zero discord.
In Ref. [32] it was shown, that for an arbitrary X state the trace distance discord can be calculated as Assuming ρ 23 , ρ 14 ∈ R simplifies Eq. ( 22) into A. Asymmetric discord dynamics In depolarizing, bit-flip and bit-phase-flip channels the parameters γ 1 , γ 2 and γ 3 evolve symmetrically with respect to noise location.This means that whenever the conditions of Eq. ( 25) or Eq. ( 26) are satisfied throughout the dynamics, trace distance discord behaves symmetrically.On the other hand, we see that in the case of Eq. ( 27), discord evolves symmetrically if and only if the parameter x has symmetric dynamics.This is equivalent to using initial state with x = 0.In the case of Eq. ( 27), this choice leads to D(X) = |γ 2 |/2.
Two exemplary cases of states evolving in subspaces of X states, defined by conditions in Eq. ( 25) - (27), are illustrated in Fig. 8.In Fig. 8(a) we see, that noise on qubit U makes the state enter or exit the green segment if and only if noise on qubit L does so.Contrary to this, the line between purple and yellow segments is crossed by noise on one qubit when noise on the other one does not make the state cross the line.Instead in Fig. 8(b), each segment boarder is crossed by noise on one qubit if and only if noise on the other one makes the state cross the line.This case is special also because decay of discord is symmetric in all sections, and thus through the whole dynamics.
The non-trivial cases occur when noise on one qubit causes the evolved state Φ p (X) to satisfy x 2 but noise on the other one leads to a state satisfying instead.This situation is illustrated in Fig. 8 (a) as the curved boarder between the purple and yellow segments.For depolarizing, bit-flip and bit-phase-flip channels this case never results to symmetric decay of discord.Thus we conclude, that symmetry of discord dynamics requires that both local channels map the state into the same segment of the state space at each value of p.For depolarizing and bit-phase-flip channels, initial state satisfying |γ 3 | ≥ |γ 1 | guarantees that the property is preserved throughout the dynamics.In the sense of Fig. 8 this means that any initial state inside the green segment never exits it during dynamics.Neither can any state outside the region enter it.As an interesting example, we study the maximally discordant mixed two-qubit states (MDM).In [24] the analytic form of MDM's was solved.They are all X states and it is easy to see that they satisfy |γ 3 | ≥ |γ 1 |, but their state dynamics is symmetric only in special cases.Thus all two-qubit MDM's have symmetric discord dynamics under depolarizing and bit-phase-flip channel.
For bit-flip channel the conditions are more restrictive.If  For bit-phase-flip channel, the set of states staying inside the region is more exclusive.Noise on qubit U preserves the property if and only if the initial state satisfies γ 2 = 0. Bit-phase-flip noise on qubit L keeps the state inside the region if and only if γ 2 = x = 0. Like in the case of the other channels, also here we see that initial states, satisfying γ 2 = 0 but violating x = 0, stay inside the region when noise affects qubit U but exit it when noise is applied on qubit L instead, resulting to asymmetric discord dynamics.
To conclude, we see that for all of the channels and all three segments there exist well-defined families of states which stay within the segment they started from.There is no preference in the direction of crossing the segment boundaries: each segment has initial states exiting it and initial states from other segments entering it during the evolution.The only exception is segment |γ 3 | ≥ |γ 1 | in depolarizing and bit-phase-flip channels: no initial state outside the segment can enter it nor does any state initially inside exit it.

B. Time-invariant quantum discord
Finally, we study the so-called time-invariant discord phenomenon, where the value of discord is not influenced by the channel.Local depolarizing noises never induce time-invariant discord.Instead bit-phase-flip and bit-flip channels do.Discord is invariant under local bit-phaseflip channel if and only if x = 0 and |γ 3 | < |γ 2 |.First of the conditions is preserved in the dynamics for all initial states satisfying it and the second one holds for the evolved state Φ U p (X) if and only if Φ L p (X) satisfies it.So we see that time-invariant discord occurs under bitphase-flip channel only if the discord dynamics is symmetric.
For a bit-flip channel, there are two disjoint families of states leading to time-invariant discord: initial states satisfying |γ 3 | ≥ |γ 1 | > 0, or |γ 3 | = 0 and ρ 23 = ρ 14 = 0.All initial states satisfying either of the conditions satisfies the condition through the whole dynamics.States in the first family have symmetric discord dynamics.Instead, for initial states in the latter family, discord dynamics is asymmetric.In fact the time-invariant discord occurs only when the local noise is applied on qubit L. If the noise affects qubit U instead, the value of discord goes to zero for all initial states of the family.

VI. SUMMARY AND DISCUSSION
We have studied the dynamics of concurrence, Bell function and trace distance discord under one-qubit channels and their combinations.The channels we considered were depolarizing, amplitude damping, bit-flip, and bit-phase-flip channels and simple combinations of depolarizing and amplitude damping channels.We saw that even though the input state, or even the dynamics of the state, is asymmetric, entanglement, non-locality and discord can decay symmetrically.We noticed that the families of asymmetric states leading to symmetric entanglement or Bell function decay are not the same for different channels or channel combinations.Thus, by measuring how concurrence or Bell function decays, one can deduce which noise was affecting the two-qubit system, which qubit it affected, and how large is the channel strength parameter p.
We also saw that, for some initial states, entanglement decay is sudden when the noise affects one qubit and asymptotic when the noise influences the other.Also, the same initial state can lead to sudden death of entanglement for one type of channel and asymptotic decay for another.For one-sided amplitude damping channel, initial states with symmetric entanglement decay can lead to asymmetric decay of Bell function and the opposite happens in depolarizing channel.This means that there is no natural hierarchy between the asymmetric decays of concurrence and non-locality.
Finally, we studied the dynamics of trace distance discord noticing that the total space of states can be divided into three disjoint regions.We saw that symmetric decay of discord requires that local noise affecting the system maps the initial state into the same region at the same value of parameter p.We characterized families of initial states leading to symmetric and asymmetric discord dynamics for each channel.We also characterized families of initial states leading to time-invariant discord and saw that for some of these families discord dynamics is symmetric but for one family time-invariant discord occurs only when noise affects qubit L.

FIG. 1 .
FIG. 1.Visualization of different noise locations.Alice sends a two-qubit system to Bob.Qubits U and L are transmitted through the upper and lower path, respectively.In (a) and (b) the local noise noise influences the qubit U and L, respectively.

FIG. 2 .
FIG. 2. Schematic illustrations (a), (b), (c), and (d) visualize four different geometric configurations for noise combinations.Here channels Φp, Ψp, and Ξp correspond to local amplitude damping or depolarizing noises with equal channel strength parameters p acting on qubits U and L.

FIG. 3 .
FIG.3.The red and blue curves correspond to situations with amplitude damping channel on qubit U and L, respectively.Influence of the noise on qubit L leads to sudden death of entanglement, but when the noise affects the qubit U , entanglement decays asymptotically.We chose the input state as ρ11 = 0.35, ρ22 = 0.4, ρ33 = 0.05 and ρ44 = 0.2.

FIG. 5 .
FIG. 5. Concurrence of an X state as a function of the channel strength parameter p under depolarizing [(a)-(b)] and amplitude damping [(c)-(d)] channel.Here the red and blue curves correspond to the cases with local noise on qubit U and L, respectively.In (a) and (c) we have chosen ρ11 = 0.2, ρ22 = 0.65, ρ33 = 0.15, and ρ44 = 0.In (b) and (d) we have chosen ρ11 = 0.2, ρ22 = 0.15, ρ33 = 0.65, and ρ44 = 0. Note that the situation is opposite to Fig. 4: now depolarizing channel causes asymmetric entanglement decay, and entanglement decays symmetrically in amplitude damping channel, instead.

FIG. 6 .
FIG.6.Concurrence of a state as a function of the channel strength parameter p for (a) amplitude damping and (b) depolarizing channel.Here the red and blue curves correspond to the cases with local noise on qubit U and L, respectively.We have chosen ρ11 = 0.9, ρ22 = 0, ρ33 = 0.08, and ρ44 = 0.02.

FIG. 7 .
FIG.7.Behavior of (a) Bell function and (b) concurrence as a function of channel parameter p.In each plot the red and blue lines correspond to amplitude damping channel affecting qubits U and L, respectively.Here we have chosen ρ11 = ρ44 = 0, ρ22 = 0.1875, and ρ33 = 0.8125.

FIG. 8 .
FIG.8.Visualization for possible trajectories of state dynamics (color online).The cases are divided between initial states satisfying (a) x = 0 and (b) x = 0.The space of states is represented by a clock-face and split into three colored segments.The channel strength p corresponds to the rotation angle θ(p) of the orange hand of the clock.The red and blue dashed circles represent possible trajectories of Φ U p (X) and Φ L p (X), respectively.As p increases, the clock hand rotates moving the evolved state along the dashed circle which corresponds to the noise location.The value of θ(0) is determined by the initial state and the behavior of θ(p) depends on the channel Φp and the initial state.

TABLE I .
Behavior of entanglement decay (ED) and decay of non-locality (ND) for different families of initial states.