Universal freezing of quantum correlations within the geometric approach.

Quantum correlations in a composite system can be measured by resorting to a geometric approach, according to which the distance from the state of the system to a suitable set of classically correlated states is considered. Here we show that all distance functions, which respect natural assumptions of invariance under transposition, convexity, and contractivity under quantum channels, give rise to geometric quantifiers of quantum correlations which exhibit the peculiar freezing phenomenon, i.e., remain constant during the evolution of a paradigmatic class of states of two qubits each independently interacting with a non-dissipative decohering environment. Our results demonstrate from first principles that freezing of geometric quantum correlations is independent of the adopted distance and therefore universal. This finding paves the way to a deeper physical interpretation and future practical exploitation of the phenomenon for noisy quantum technologies.

In quantum mechanics, the mathematical description of a composite quantum system is based on both the superposition principle and the tensorial structure of the Hilbert space associated with it. The coexistence of these two principles makes the properties of generic states of a composite quantum system particularly weird, in the sense that they cannot be reproduced by any state of a classical system. Some of the most striking non-classical properties exhibited by quantum states can be collected under the name of quantum correlations12. Nowadays, there is universal consensus on the fact that the quantum correlations shared by two subsystems in a global pure state are entirely captured by entanglement and can be quantified by any valid entanglement measure13. On the other hand, it is also now quite clear that there exist non-entangled mixed states still manifesting some non-classical features, such as an unavoidable disturbance due to local measurements, which embodies the concept of quantum discord45. Therefore, in realistic open quantum systems, entanglement may represent only a portion, sometimes negligible, of the quantumness of correlations, while more general figures of merit to quantify quantum correlations are provided by suitable measures of discord-type correlations26. Despite an intense recent activity in investigating interpretation, quantification, and applications of discord and related quantifiers of quantum correlations2, these quantities still remain less understood than entanglement.

In order to unveil the most profound signatures of quantumness in composite systems, it is essential to identify mathematically rigorous and physically meaningful properties that differentiate the notion of discord-type quantum correlations from that of entanglement (and of classical correlations), and are manifested by any valid measure thereof. One such property common to generic discord-type correlation measures is, for instance, the absence of monogamy7. Besides the fundamental implications, this area of investigation has a technological motivation89, since quantum correlations beyond and even without entanglement have been shown to play a resource role for certain schemes of quantum computation101112, communication13141516, and metrology171819. Finding valuable and general traits of these quantum correlation resources, in particular for what concerns their dynamical preservation during the unavoidable interaction of a principal quantum systems with the surrounding environments, constitutes an important aim with a clear potential to lead to useful recipes for their practical exploitation.

Numerous works have in fact investigated the dynamics of general quantum correlations in open quantum systems undergoing various types of Markovian or non-Markovian evolution, as reviewed e.g. in Refs. 2,20,21. Although different measures of quantum correlations can exhibit distinct features and impose inequivalent orderings on the set of quantum states, it has emerged as a general trait that discord-type quantum correlations are more robust than entanglement against noise2223242526 (see also2728, for a critical assessment) and cannot generally vanish at a finite evolution time (due to the fact that zero-discord states are of null measure25, while entanglement can suffer so-called sudden death2930. However, a fascinating and nontrivial phenomenon of extreme robustness to noise exhibited by general quantum correlations deserves special attention, and is the subject of our investigation.

Namely, under local non-dissipative decoherence evolutions, it has been observed that a number of known discord-type measures all remain constant (‘frozen’) for a finite time interval in Markovian conditions63132, and for multiple intervals333435, or forever36, in non-Markovian conditions, when considering two non-interacting qubits initially in a specific class of Bell-diagonal states. This freezing phenomenon, not exhibited by any measure of entanglement, is quite appealing since it implies that every protocol relying on discord-type quantum correlations as a resource will run with a performance unaffected by noise in the specific dynamical conditions. Currently, the occurrence of freezing has been investigated by explicitly considering the evaluation of specific discord-type measures on a case by case basis63133 However, it is natural to ask whether this phenomenon is a mere mathematical accident due to the particular choices of quantum correlations quantifiers, or whether it must manifest independently of the adopted measure, thus having a universal character and promising to bear a deep physical meaning. This work addresses such an issue.

We prove that freezing occurs for any geometric measure of quantum correlations, whenever the distance defining the measure respects a minimal set of physical assumptions, namely dynamical contractivity under quantum channels, invariance under transposition, and convexity. The freezing phenomenon is therefore revealed as universal within the geometric approach to quantum correlations. Notice that our work differs from other complementary investigations of the freezing phenomenon3237. In particular, in a recent work37, the authors provide necessary and sufficient conditions for a general state to exhibit freezing under non-dissipative decoherence, according to some specific measure of discord. Here, instead, we focus on a specific class of initial states, and we identify the minimal set of conditions that any general distance-based measure of discord needs to satisfy in order to freeze. On some random family of initial states, it is certainly possible to see freezing according to one discord-type measure but not to another. What we prove here is that, for the specific class of Bell-diagonal states identified in Refs. 6, 31, all bona fide geometric quantifiers of quantum correlations (respecting the three physical assumptions mentioned above) undergo the same dynamics, featuring the freezing phenomenon. In proving the main result, we also introduce and characterise a global quantum control channel which can completely invert decoherence on a subset of Bell-diagonal two-qubit states, and can be of independent interest.

The paper is organised as follows. We first summarise the properties that any valid measure of quantum correlations is expected to hold, and provide general definitions for distance-based measures, distinguishing between those for entanglement and those for discord-type correlations. In the Results section, we illustrate the freezing phenomenon in geometric terms, by considering for convenience the specific case of the Bures distance-based measure of quantum correlations. We then present the main result, by proving that freezing must happen for any bona fide distance-based discord-type measure. Finally, we offer our conclusions in the Discussion section. Some technical bits are deferred to the Methods section.

Measures of quantum correlations

Here we recall the requirements that a valid measure of quantum correlations is expected to have, later focusing specifically on geometric (distance-based) definitions. For more details on the quantification of quantum correlations, the reader is referred to recent review articles12.

In this paper, we consider a state of a two-qubit system, with subsystems tagged and . As is well known138, entanglement quantifiers capture the degree of non-separability of the state of the global system; the corresponding distance-based measures of entanglement are calculated from the set of separable states39, i.e., states which can be written as convex combination of product states,

where is a probability distribution, while and are arbitrary ensembles of states for subsystem and , respectively.

Quantifiers of quantum correlations other than entanglement, the so-called discord-type measures, capture instead the minimal degree of disturbance on the state after local projective measurements on the system245. The projective measurements can be performed either on a subsystem only (which gives rise to one-way, asymmetric discord-type measures) or on both subsystems (which gives rise to two-way, symmetric discord-type measures). These two versions of discord-type measures have valuable operational meanings in different contexts2 and the corresponding distance-based measures are calculated, respectively, from the set of classical-quantum (CQ) and classical-classical (CC) states41. Explicitly, a CQ state (with respect to measurements on subsystem ) is a particular type of separable state, which can be written in the form

where is a probability distribution, denotes an orthonormal basis for subsystem A, and is an arbitrary ensemble of states for subsystem . Similarly, a CC state can be written in the form

where is a joint probability distribution, while and denote orthonormal bases for subsystem and , respectively. Clearly, the set of CC states is contained in the set of CQ states, which is a subset of the set of separable states.

In general, if not explicitly written, by “classical” states we hereafter mean CC states, i.e., states which are diagonal in a product basis, as defined by Eq. (3); these states correspond merely to the embedding of a bipartite probability distribution into the quantum formalism. We further specify that for Bell-diagonal states42, namely the specific class of two-qubit states considered in this work, the two notions of discord are completely equivalent, therefore our conclusions about the universality of the freezing will apply indifferently to both one-way and two-way geometric measures of discord-type correlations.

From a quantitative point of view, a valid entanglement measure, also known as an entanglement monotone3, is any real and nonnegative function on the set of states satisfying the following basic axioms143:

(E.i) if is a separable state as defined in Eq. (1);

(E.ii) is invariant under local unitaries, i.e. for any state and any local unitary operation acting on subsystem ();

(E.iii) is monotonically nonincreasing under local operations and classical communication (LOCC), i.e. for any state and any LOCC channel .

Furthermore, additional properties for an entanglement measure can include convexity,

(E.iv) is convex, i.e. , for any pair of states and any .

Notice, however, that while convexity is physically desirable (as it would mean that entanglement cannnot increase by mixing states), it is not an essential property, since there are valid entanglement monotones which are not convex44.

The theory of quantum correlations other than entanglement is not completely developed yet2,38, but we can nonetheless identify some desiderata for any quantifier thereof. A (two-way) discord-type measure is any real and nonnegative function on the set of states satisfying the following requirements:

(Q.i) if is a classical state as defined in Eq. (3);

(Q.ii) is invariant under local unitaries, i.e. for any state and any local unitary operation acting on subsystem ();

(Q.iii) is monotonically nonincreasing under local commutativity preserving quantum channels , i.e. for any state and any commutativity preserving map on subsystem , that is, when for arbitrary marginal states and ;

(Q.iv) reduces to an entanglement measure for pure states, i.e. for any pure state .

We introduce property (Q.iii) in analogy with property (E.iii) for entanglement. Namely, it is known that local commutativity preserving channels cannot create discord-type correlations, as they leave the set of classical states invariant45. We thus require that any valid measure of discord-type correlations should be monotonically nonincreasing under such channels. Notice, in particular, that for two qubits these channels include local unital channels22.

The above requirements need to be slightly modified if a one-way discord-type measure Q→, say with measurements on , is considered. Specifically, property (Q.i) becomes: Q→(ρ)=0 if ρ is a CQ state as defined in Eq. (2). Furthermore, a stricter monotonicity requirement supplements (Q.iii) for all valid one-way discord-type measures64647, namely

(Q.iii.bis) Q→ is monotonically nonincreasing under arbitrary local quantum channels on the unmeasured subsystem , that is, for any state ρ and any completely positive trace-preserving (CPTP) map on subsystem .

Properties (Q.ii) and (Q.iv) apply equally to two-way and one-way discord-type measures. The latter property just signifies that, in pure bipartite states, there is a unique kind of quantum correlations, arising in all but tensor product states. Even correlations stronger than entanglement, such as steering and nonlocality, just collapse back to non-separability in the case of pure states.

In order to investigate the freezing phenomenon, in this paper we resort to a geometric approach to define a very general class of valid measures of quantum correlations. According to such an approach, the entanglement and the discord-type correlations of a state can be quantified as the minimal distance from to the sets and of separable and classical states, respectively4048. In formulae,

where separable states are defined by Eq. (1), while classical states are defined respectively by Eqs. (3) and (2) depending on whether a two-way or one-way discord-type measure is considered. In these definitions, can denote in principle any suitable distance on the set of quantum states.

In order for the geometric measures and to respect the essential properties listed above, the distance needs to satisfy certain mathematical requirements49. Here, we identify a minimal set of three such requirements, that will be said to characterise as a bona fide distance. Given any states , , , and , these are:

(D.i) Contractivity under CPTP maps, i.e.

for any CPTP map ;

(D.ii) Invariance under transposition

(D.iii) Joint convexity, i.e.

for any .

Let us comment on the physical significance of these requirements.

On the one hand, the contractivity property (D.i) of , Eq. (6) is a fundamental requirement for a distance in quantum information theory49, and has two purposes. First, it makes a statistically relevant distance, due to the fact that non-invertible CPTP maps are the mathematical counterparts of noise and the latter cannot lead to any increase in the information related to the distinguishability of quantum states5051. Second, (D.i) makes the distance-based measures of entanglement and quantum correlations [Eqs. (4) and (5)] physically meaningful, by implying the essential properties (E.ii), (E.iii), (Q.ii), (Q.iii), and also (Q.iii.bis), where the latter applies to the corresponding distance-based one-way discord-type measure defined by choosing the set of CQ (rather than CC) states in Eq. (5). Notice that (D.i) implies in particular the standard property of invariance of the distance under unitary operations, since they can be seen as reversible CPTP maps; namely, , for any pair of states , , and any unitary .

On the other hand, the invariance of a distance under transposition, Eq. (7), is not typically discussed in the literature. However, transposition of an hermitian matrix, which amounts to complex conjugation in the computational basis, corresponds to a reflection in a -dimensional hyperplane. Property (D.ii) thus means that a distance on the set of quantum states is assumed to be invariant under reflections, which appears as a fairly natural requirement50. Notice that, together with (D.i), this property implies invariance of the distance under antiunitary operations. Any antiunitary matrix can be expressed as , where is a unitary matrix and denotes complex conjugation in the computational basis. We have then , for any pair of states , .

Finally, the joint convexity property (D.iii) of , Eq. (8), is also quite intuitive and it makes the corresponding distance-based entanglement measure convex, implying the desirable property (E.iv). Notice however that discord-type measures are, by contrast, neither convex nor concave52, as the set of classical states is not a convex set25.

There are a number of known distances which satisfy the three physical assumptions listed above50, and have been employed to define valid geometric measures of quantum correlations. Suitable examples for include in particular the relative entropy4053, the squared Bures distance65455, the squared Hellinger distance185657, and the trace (or Schatten one-norm) distance585960. Contrarily, the Hilbert-Schmidt distance does not respect the contractivity property (D.i), and earlier attempts to adopt it to define geometric measures of quantum correlations4861 have led to inconsistencies6263.

In this paper, we label a generic distance obeying properties (D.i), (D.ii), and (D.iii) as a bona fide one, and the associated distance-based quantities and as bona fide measures of entanglement and discord-type correlations respectively. Therefore, the main result of this paper will be a proof of the universality of the freezing for all geometric measures of quantum correlations constructed via bona fide distances as formalised in this Section.

Results

Freezing of quantum correlations measured by Bures distance

We now present the freezing phenomenon from a geometric perspective, by employing a particular bona fide measure of quantum correlations, that is the Bures distance-based measure6556465. We first recall all the basic ingredients for the complete description of the phenomenon.

The Bures distance between two states and is defined as

where

is the Uhlmann fidelity66.

The Bures distance arises from a specific case of a general family of Riemannian contractive metrics on the set of density matrices, characterised by Petz67 following the work by Morozova and Čencov68. It can be connected operationally to the success probability in ambiguous quantum state discrimination, and it has been successfully employed to define geometric measures of entanglement, quantum, classical, and total correlations6225564656970. Specifically, the Bures entanglement measure and the Bures discord-type measure are defined by choosing in Eqs. (4) and (5) the squared Bures distance , which is a bona fide one obeying properties (D.i), (D.ii), and (D.iii).

Bell-diagonal (BD) states, also referred to as T-states or two-qubit states with maximally mixed marginals42, are structurally simple states which nonetheless remain of high relevance to theoretical and experimental research in quantum information, as they include the well-known Bell and Werner states39 and can be employed as resources for operational tasks such as entanglement activation and distribution via discord-type correlations47717273. BD states are by definition diagonal in the basis of the four maximally entangled Bell states, and their Bloch representation in the computational basis is

where , is the identity matrix and are the Pauli matrices. Because is positive semidefinite, the vector characterising any BD state is constrained to lie within the tetrahedron with vertices , , and 424852; the vertices represent respectively the four Bell states

For convenience, in the following, we shall specify an arbitrary BD state by referring equivalently to its defining triple .

A non-dissipative quantum channel acting on a qubit induces decoherence with no excitation exchange between the qubit and its environment. We consider the evolution of two non-interacting qubits undergoing local identical non-dissipative decoherence channels. The action of any such channel on each single qubit is characterised by the following Lindblad operator631,

where is the decoherence rate, is the reduced state of subsystem () and represents the direction of the noise. Namely, the choice of respectively identifies decoherence in the Pauli basis for each qubit, and the corresponding channels are known in the quantum computing language as bit flip (), bit-phase flip (), and phase flip () channels3149. It is worth noting that one can easily derive the dynamics of the composite two-qubit system from the dynamics of the single qubits, since each of the two qubits is locally interacting only with its own environments (and not with the other qubit), so that they have independent dynamical evolutions21.

Equivalently, the evolution of a two-qubit state under local non-dissipative decoherence channels can be obtained in the operator-sum representation by the map

where the single-qubit Kraus operators are

with and being the direction of the noise as in Eq. (13). Interestingly, from Eq. (14) one can easily see that two non-interacting qubits initially in a BD state, undergoing local identical non-dissipative channels, maintain the BD structure for all time. More specifically, the triple characterizing the BD evolved state can be written as follows

where is the triple characterising the initial BD state .

For non-interacting qubits initially in a BD state , undergoing identical local non-dissipative channels, the freezing phenomenon for discord-type quantum correlations occurs if specific initial state conditions are satisfied. For convenience and without loss of generality, from now on we focus our analysis on the phase flip (pure dephasing) channel (), for which these initial conditions consist of the triples such that

The suitable initial conditions for bit flip () or bit-phase flip () channels can be obtained by setting analogous relations among the coefficients (Ref. 6). The BD states satisfying the constraint of Eq. (17) distribute within a two-dimensional surface inside the tetrahedron of all BD states, which is shown in Fig. 1 and will be referred to herein as (phase flip) freezing surface.

Figure 1

The phase flip freezing surface (meshed cyan) within the tetrahedron of all BD states (light yellow) represented in the space. The surface contains all and only the BD states with triple , and thus accommodates all the BD states respecting Eq. (17). Solid black lines represent the classical BD states, which lie on the axes. The dotted red lines represent the threshold points on the surface when , which occurs at the time defined in Eq. (18). For any state obeying the initial conditions of Eq. (17), we show that bona fide discord-type quantum correlations are frozen under local phase flip channels up to the time . As an example, the dashed blue line represents the dynamical trajectory of the initial BD state , which evolves under local phase flip channels moving towards the -axis with increasing time; the discord-type correlations are frozen in the initial segment of the trajectory up to the intersection with the red dotted line, and decay exponentially afterwards, as plotted in Fig. 2.

From Eq. (16) one can see that the time evolved state is a BD state characterised by the triple {c1(0)e-2γt,-c1(0)c1(0)e-2γt, c1(0)}, which means that it remains confined within the freezing surface at any time. An example of this dynamical trajectory is represented in Fig. 1 by the dashed blue line.

The above dynamical conditions give rise to a very peculiar evolution of the quantum correlations present in the time evolved state. Namely, defining a threshold time by

we find that the Bures distance-based measure of discord stays constant (freezes) for and then decays exponentially from onwards, as is shown in Fig. 2. This can be straightforwardly shown by exploiting the available closed formula for Bures discord-type correlations of BD states655. On the contrary, entanglement measured e.g. by undergoes a typical sudden death at a finite time29. We stress again that this behaviour of quantum correlations, here illustrated for , has been independently observed (on a case by case basis) for several valid discord-type measures in the aforementioned dynamical conditions6: this paper will provide a rigorous basis to establish its universality within the bona fide geometric approach.

Figure 2

The paradigmatic freezing of Bures distance-based quantum correlations for an initial BD state of the form . The solid blue line represents the time evolution of discord-type correlations and the dashed red line represents the time evolution of Bures distance-based entanglement.

The freezing phenomenon can be understood in geometric terms by looking at Fig. 1, which represents the phase flip freezing surface containing BD states of the form , thus containing all the states with initial conditions identified in Eq. (17). The solid black lines represent the axes in the space, which correspond to all the classical BD states. The dashed blue line represents a particular state evolution under local phase flip channels, leading to the freezing phenomenon. The intersection between the dotted red line and the dashed blue line represents the time such that , i.e., the threshold time (18) corresponding to the sudden change from freezing to decaying of quantum correlations.

For , i.e. when , one of the closest classical states to the evolved BD state is its Euclidean orthogonal projection onto the -axis, i.e. the BD classical state with triple . From Eq. (5) we know that the discord-type quantum correlations of for are given exactly by the (squared) Bures distance between the evolving state and this closest classical state . Interestingly, one can observe that this distance is constant for any , which indeed implies that the quantum correlations of are frozen for any , given arbitrary initial conditions on the freezing surface defined by (17).

For , i.e. when , one of the closest classical states to the evolved BD state is instead its Euclidean orthogonal projection onto the -axis, i.e. the BD classical state with triple , which is independent of time. Therefore the quantum correlations of the evolved state decrease for any , as the distance between the evolving state and the steady closest classical state decreases for any .

For two-qubit BD states of the form (17) undergoing local pure dephasing evolutions, freezing (up to a time ) of geometric quantum correlations measured by the Bures distance thus relies on the following two properties:

(F.i) (according to the Bures distance) one of the closest classical states to the evolved BD state is the classical BD state when ;

(F.ii) the (Bures) distance between the BD states with triples and does not depend on , that is

The main result of this paper will be to show that these two properties are satisfied by any contractive, transposition invariant, and convex distance, thus implying the freezing phenomenon for any bona fide distance-based measure of discord-type quantum correlations as defined above. Contrarily, we remark that the non-contractive Hilbert-Schmidt distance satisfies only the first property (F.i), whereas it does not manifest the kind of translational invariance expressed in Eq. (19), due to the fact that the trajectory of the evolved state is not parallel to the -axis according to the Euclidean geometry, as is shown in Fig. 1. As a result, the Hilbert-Schmidt geometric discord48, which is not a bona fide measure63, does not manifest freezing in the considered dynamical conditions, as previously observed6.

Freezing of quantum correlations for all bona fide geometric measures

We now prove that for any distance respecting the bona fide requirements (D.i), (D.ii), and (D.iii), the aforementioned freezing properties (F.i) and (F.ii) are satisfied.

Let us begin by proving that (D.i) (F.ii), i.e., that the very peculiar invariance on the phase flip freezing surface expressed in Eq. (19) follows from the contractivity of a distance .

Theorem 1. Any contractive distance satisfies the following translational invariance properties

and

where denotes any BD state on the phase flip freezing surface.

Proof. We first prove Eq. (20). Considering an evolution due to local complete pure dephasing channels acting on the two qubits, corresponding to the dynamical map as given in Eq. (14) for and , we get the inequality

where the first equality is due to the fact that

and the final inequality is due to contractivity of the distance .

We now introduce a global two-qubit rephasing channel with operator-sum representation

where

with the 8 Kraus operators satisfying , where are the four pure Bell states defined in Eq. (12). The global rephasing map , which is a physical CPTP map for all , transforms any two-qubit state into a BD state lying on the phase flip freezing surface, i.e., , where is a BD state with characteristic triple given by , with . Specifically, the action of on a BD state which already belongs to the phase flip freezing surface is

This map is therefore able to restore the lost coherence for any (even completely dephased) BD state on the freezing surface, thus effectively reverting their decoherence process.

We have then the inequality

where the first equality is due to the fact that

and the final inequality is again due to contractivity of the distance . By putting together the two opposite inequalities (22) and (27), we immediately get the invariance of Eq. (20) for any contractive distance.

To prove now the claim of Eq. (21), we introduce the unitary such that

By exploiting the invariance under unitaries of any contractive distance, and the just proven invariance expressed by Eq. (20), we finally have

that establishes Eq. (21), thus concluding the proof of the Theorem.

■

We now prove that (D.i), (D.ii), (D.iii) (F.i), that is, one of the closest classical states to a BD state of the form , with , is always the BD classical state of the form for any bona fide distance as defined above. We divide this result in two main steps, represented by Theorem 2 and Theorem 3. To bridge between the two Theorems we need four Lemmas which are formulated and proved in the Methods.

Let us begin by the following powerful result, which applies to all two-qubit BD states.

Theorem 2. According to any contractive, transposition invariant, and convex distance, one of the closest classical states to a BD state is always a BD classical state,

for some index and some coefficient .

Proof. For an arbitrary two-qubit state , described in the Bloch representation

by the triple , there exists another two-qubit state with the associated triple 66574. The specific triple for a two-qubit BD state as defined in Eq. (11) is and ; notice hence that for BD states.

Given any BD state , and any general two-qubit state , it holds that

where we have introduced an antiunitary operator which maps to . Explicitly, it acts in the computational basis as , where denotes complex conjugation in the standard basis, which amounts to transposition for quantum states. Eq. (33) thus follows by exploiting the invariance of under antiunitary operations which is provided by properties (D.i) and (D.ii).

By further using the convexity of the distance in the second argument,

which is automatically implied by its joint convexity in property (D.iii), one can also show that for any BD state it holds that

where is described in the Bloch representation by the triple . Explicitly,

where we have used in the first equality, convexity of in the inequality and in the final equality.

We will now consider the distance from to the (larger) set of two-qubit CQ states, and show that its minimum can be attained by a fully classical (CC) BD state, hence proving the main result of the theorem. Recall that any CQ two-qubit state is of the form

where , is an orthonormal basis for qubit and and are arbitrary states of qubit . Such a CQ state will have the associated triple , where

with being the Pauli matrices. For any state in this form, a second state

can be derived using the identities

where and are the traceless part of and . This state is manifestly CQ and it can be easily verified that it will have the associated triple with . From the inequality (35), we have in particular that for any BD state and any CQ state , so that in order to minimise it suffices to restrict ourselves to CQ states with associated triple .

Temporarily, we relax the restriction that is of unit length and consider the distance from the (even larger) set of states for which . This is a convex set and so, due to the convexity of the distance, any local minimum will be a global one. We can now use a trick analogous to the one used for Eq. (33), this time between with and and with and . We then see that

where we have introduced the unitary operator such that for any BD state , and we have exploited (D.i). A similar result holds when considering distances from states with either or , by using, respectively and . Also, a similar result holds by switching the vectors we consider, for any . From these observations we have that, if for some index , then the minimum distance is attained for , and viceversa.

We can then restrict our attention to states with and , where the index sets the nonzero vector element. From the previous results, we notice in fact that minimisation only needs to be performed over and as the distance can only decrease under any variation in any other single element. Furthermore, and appear only as a product in the density matrix, never on their own. This means that minimising over both is equivalent to setting and minimising only over , thus allowing us to reimpose the restriction that , thus coming back to analyse the distance from to CQ states. The remaining states over which the minimisation in the single parameter needs to be performed amount exactly to the set of BD classical states (aligned on the axes in Fig. 1), hence finding the minimum among these will return the global minimum for the distance from an arbitrary BD state to the set of two-qubit classical states, proving the claim.

■

What is left at this point is to perform the final minimisation, namely to determine exactly the values of the index and of the coefficient that entirely specify the closest classical state of Eq. (31) as a function of the coefficients defining any given BD state . When we restrict ourselves to BD states belonging to the phase flip freezing surface of Eq. (17), the solution is provided by Theorem 3, which makes use of the auxiliary results proven in the Methods.

Theorem 3. According to any convex and contractive distance, one of the closest classical states to a BD state of the form is

when , the BD classical state , i.e. the one with and in Eq. (31);

when , the BD classical state , i.e. the one with and in Eq. (31).

Proof. According to Theorem 2, one of the closest classical states to any BD state is a classical BD state. According to Lemmas A1 and A2, the closest classical BD state to a BD state of the form is either {c1,0,0} or {0,0,c3}. Finally, according to Lemmas 3 and 4, if then one of the closest classical BD states to is {c1,0,0}, whereas if then one of the closest classical BD states to is {0,0,c3}.

■

From a physical perspective, the most relevant implication of Theorem 1 and Theorem 3 is that the freezing phenomenon, as described earlier adopting the guiding example of the Bures distance, occurs in fact for any bona fide distance-based measure of quantum correlations whose underlying distance is invariant under transposition, convex and contractive.

Namely, the condition is equivalent to , so that from Theorem 3 we have that one of the closest classical states to the evolved BD state is for any . Therefore, referring to the definition of distance-based quantifier of discord-type quantum correlations given in Eq. (5), one has

and according to Theorem 1 this is constant and equal to for any .

On the other hand, the condition is equivalent to , so that from Theorem 3 we have that one of the closest classical states to the evolved BD state is the BD state for any . Therefore, in this case the geometric quantifier of discord-type quantum correlations is

that, due to the contractivity of , has to be monotonically nonincreasing for any , eventually decaying to zero.

Discussion

In this paper we have established from first principles the general character of an intriguing dynamical trait of quantum correlations other than entanglement, namely their freezing under given environmental and initial conditions. This phenomenon manifests for the class of Bell-diagonal states of two qubits, which often constitute the simplest yet highly relevant class of states for which one is able to analytically calculate measures of quantum, classical, and total correlations, for instance by distance-based (geometric) quantifiers. In particular, we have shown here that a specific class of Bell-diagonal states of two qubits, each undergoing local non-dissipative decoherence, manifests freezing of discord-type quantum correlations whenever the distance adopted to measure them is assumed to be invariant under transposition, dynamically contractive and convex. As these physical properties are instrumental to define valid distance-based measures of correlations, our result means that freezing of quantum correlations occurs independently of the adopted distance and is therefore universal within a bona fide geometric approach.

Frozen quantum correlations have been verified both theoretically6243132333775 and experimentally767778798081 by using specific measures of quantum correlations26, but until now it was an open problem whether all suitable discord-type quantifiers (including potentially new ones yet to be defined) would freeze in the same dynamical conditions. Our work rigorously contributes to the settling of this problem and provides elegant evidence strongly supporting the conclusion that freezing of quantum correlations is a natural physical phenomenon and not merely a mathematical accident. Notice that freezing in BD states, as described in this paper, has also been observed for some discord-type measures which do not manifestly enjoy a distance-based definition, such as the local quantum uncertainty618 and the interferometric power19. This leaves some room for further research aimed to prove the occurrence of freezing from only the basic properties (Q.i), (Q.ii), and (Q.iii) of quantum correlations, possibly without the need to invoke a geometric approach as considered in this work. Alternatively, our result might suggest that all measures of discord could possibly be recast into a geometric form via some bona fide distance, at least when restricted to BD states of two qubits (this is the case, for instance, for the conventional entropic measure of discord5, which becomes equivalent to the relative entropy-based discord31 for BD states); this would also be an interesting direction to explore, in a more mathematical context of information geometry.

We further remark that, although we have explicitly considered Markovian evolutions in our analysis, the freezing of quantum correlations also occurs in the presence of non-Markovian channels which can be described by a master equation with a memory kernel, as in the case of pure dephasing or decoherence under classical random external fields3334353682. Indeed, in these cases the dynamics of BD states can be formally written as in Eq. (13), but with replaced by a more general time-dependent rate . This can give rise to a dynamics with multiple intervals of constant discord333482, or discord frozen forever36 depending on the initial conditions. By our analysis, we conclude that those fascinating features, which might be observable e.g. in the dynamics of impurity atoms in Bose-Einstein condensates3683 are universal too and manifest when probed by any bona fide geometric discord-type measure .

Within this paper we introduced an intriguing global rephasing channel, which is able to reverse the effects of decoherence for certain two-qubit BD states. This physical CPTP channel may be of interest for applications other than proving the universality of the freezing phenomenon, for example quantum error correction49, where it is desirable to combat the effects of noise, typically manifesting via local bit flip, phase flip, or bit-phase flip channels. For suitable BD states, all these errors can be corrected by global maps such as the one in Eq. (24). We also note that the action of this channel resembles (but is different from) the physical situation of refocusing by dynamical decoupling control on qubits undergoing low-frequency pure dephasing8485. The further characterisation and experimental implementation of our global rephasing map for quantum information processing calls for an independent analysis which is beyond the scope of this paper.

From a fundamental perspective it is important to understand the deeper physical origin of frozen quantum correlations. There are reasons to reckon that the phenomenon is related to the complementary freezing of classical correlations. Typically, as observed so far using specific quantifiers, given particular dynamical and initial conditions as studied here, quantum correlations are initially frozen and classical correlations decay but, after a characteristic time , classical correlations freeze and quantum correlations decay316574. This has been linked to the finite-time emergence of the classical pointer basis within the fundamental theory of decoherence75788186. Nevertheless, classical correlations are still inconsistently defined in geometric approaches6587, and it remains unknown whether they exhibit freezing after for any bona fide distance. This is certainly an aspect deserving further investigation.

Very recently, some of us have shown that an even more fundamental property of quantum systems, namely coherence88 in a reference basis, can also remain frozen under local nondissipative decoherence channels for the same class of initial states as studied here89. Such a result holds more generally for a class of -qubit states with maximally mixed marginals (for any even ), which include and extend the two-qubit set discussed in this work. This suggests that multiqubit and multipartite quantum correlations can freeze as well under the same dynamical conditions3790, and the methods of this work can be readily employed to prove the universality of freezing within the geometric approach, in such a more general instance as well.

Our result has also an impact from an applicative point of view. The property of being unaffected by the noise for a given period of time makes quantum correlations other than entanglement important for emergent quantum technologies89. Despite numerous basic experimental investigations, this resilience has yet to be properly exploited as a resource for quantum enhanced protocols e.g. in communication, computation, sensing and metrology. The universality of the freezing phenomenon for geometric quantum correlations, in paradigmatic quantum states and dynamical evolutions as shown here, promises to motivate further research in this context.

Methods

Here we derive some technical results needed for the proof of Theorem 3.

Lemma A1. According to any contractive distance , it holds that:

among the BD classical states belonging to the -axis, the closest state to a BD state of the form is the orthogonal projection of onto the -axis, i.e. ;

among the BD classical states belonging to the -axis, the closest state to a BD state of the form is the orthogonal projection of onto the -axis, i.e. ;

Proof. Regarding point (1), we need to prove that for any

In fact

where the first equality is due to Theorem 1, which holds for any contractive distance, the second equality is due to the fact that

with representing complete local pure dephasing towards the -axis, and finally the inequality is due to contractivity of the distance . Regarding point (2) we need to prove that for any

In fact

where the first equality is due to Theorem 1, which holds for any contractive distance, the second equality is due to the fact that

with representing complete local pure dephasing towards the -axis, and finally the inequality is due to contractivity of the distance .

■

An analogous result does not hold for the classical BD states lying on the -axis. However, due to the following Lemma 2, we can discard the classical BD states on the -axis in order to find out the closest classical BD state to a BD state lying on the phase flip freezing surface.

Lemma A2. According to any contractive distance , it holds that:

the BD classical state on the -axis closest to the BD state is closer to than any classical BD state belonging to the -axis;

the BD classical state on the -axis closest to the BD state is closer to than any classical BD state belonging to the -axis.

Proof. Regarding point (1), we need to prove that for any

In fact

where the first equality is due to Theorem 1, which holds for any contractive distance, the second equality is due to the fact that

with representing complete local pure dephasing towards the -axis, and finally the inequality is due to contractivity of the distance .

Regarding point (2), we need to prove that

In fact

where the first equality is due to Theorem 1, the second equality is due to the fact that

with representing complete local pure dephasing towards the -axis, and finally the inequality is due to contractivity of the distance .

■

Lemma A3. According to any contractive distance , if then

Proof. Let us suppose that , then we have

where the first equality is due to the fact that

with and being unitaries, whereas the second equality is due to unitary invariance of any contractive distance .

■

Lemma A4. According to any contractive distance , for any the following holds:

Proof. Regarding Eq. (54), we have

where the equality is due to the fact that

with representing local pure dephasing towards the -axis until the time such that , whereas the inequality is due to contractivity.

Regarding Eq. (55) we have

where the equality is due to the fact that

with representing local pure dephasing towards the -axis until the time such that , whereas the inequality is due to contractivity.

■

Additional Information

How to cite this article: Cianciaruso, M. et al. Universal freezing of quantum correlations within the geometric approach. Sci. Rep. 5, 10177; doi: 10.1038/srep10177 (2015).