Asymmetric quantum hypothesis testing with Gaussian states

Asymmetric quantum hypothesis testing with Gaussian states

Gaetana Spedalieri Department of Computer Science, University of York, York YO10 5GH, United Kingdom    Samuel L. Braunstein Department of Computer Science, University of York, York YO10 5GH, United Kingdom
Abstract

We consider the asymmetric formulation of quantum hypothesis testing, where two quantum hypotheses have different associated costs. In this problem, the aim is to minimize the probability of false negatives and the optimal performance is provided by the quantum Hoeffding bound. After a brief review of these notions, we show how this bound can be simplified for pure states. We then provide a general recipe for its computation in the case of multimode Gaussian states, also showing its connection with other easier-to-compute lower bounds. In particular, we provide analytical formulae and numerical results for important classes of one- and two-mode Gaussian states.

pacs:
03.67.-a, 89.70.Cf, 03.67.Hk, 03.65.Ta, 02.10.Ud

I Introduction

Quantum hypothesis testing (QHT) is a fundamental topic in quantum information theory Wilde (); Nielsen (), playing a non-trivial role in protocols of quantum communication and quantum cryptography QCprot (); QCprot2 (). The typical formulation of QHT is given in terms of quantum state discrimination QHT (); QHT2 (); RMP (); Helstrom (), where a certain number of generally non-orthogonal quantum states (the quantum hypotheses) have to be discriminated by means of a quantum measurement. In particular, the simplest scenario regards the statistical discrimination between two non-orthogonal quantum states, corresponding to the ‘null’ and the ‘alternative’ quantum hypotheses, occurring with some a priori probabilities. In symmetric testing, these hypotheses have the same cost QHT (); QHT2 (); Helstrom () and the goal is to minimize the mean error probability of confusing them by suitably optimizing the quantum measurement.

For such a basic problem, we know closed analytical formulae identifying both the minimum error probability, given by the Helstrom bound Helstrom (), and the optimal quantum detection, expressed in terms of the Helstrom matrix Helstrom (). Furthermore, we can also use an easier-to-compute bound which becomes tight in asymptotic conditions. This is the recently-introduced quantum Chernoff bound QCB (), for which we know simple formulae in the case of multi-mode Gaussian states QCB2 (), (i.e., those states with Gaussian Wigner function RMP ()).

In this paper, we consider asymmetric QHT, where two quantum hypotheses have different associated costs QHT (); QHT2 (); Helstrom (). In this approach, we aim to minimize the probability that the alternative hypothesis is confused for the null hypothesis, an error which is known as ‘false negative’. This minimization has to be done by suitably constraining the probability of another possible error, known as a ‘false positive’, where the null hypothesis is confused for the alternative hypothesis. This is clearly the best approach for instance in medical-type testing, where the null hypothesis typically represents absence of a disease, while the alternative corresponds to the presence of a disease.

Asymmetric QHT is typically formulated as a multi-copy discrimination problem, where a large number of copies of the two possible states are prepared and subjected to a collective quantum measurement. From this point of view, the aim is to maximize the error-exponent describing the exponential decay of the false negatives, while placing a reasonable constraint on the false positives. For this calculation, we can rely on two mathematical tools. The first is the quantum relative entropy RMP () between the two states, while the other is the recently-introduced quantum Hoeffding bound (QHB) QHB (), which performs the optimization of the error-exponent while providing a better control on the false positives.

In this work, we start by giving some basic notions on asymmetric QHT and briefly reviewing the QHB, also showing how its computation simply reduces to the quantum fidelity Fidelity () in the presence of pure states. Then, we provide a general recipe for computing this bound in the case of multimode Gaussian states, for which it can be expressed in terms of their first- and second-order statistical moments. In the general multimode case, we derive a relation between the QHB and other easier-to-compute bounds, which are based on well-known mathematical inequalities. Finally, we derive analytical formulas and numerical results for the most important classes of one-mode and two-mode Gaussian states.

By developing the theory of asymmetric QHT for Gaussian states, our work could be useful in tasks and protocols involving Gaussian quantum information RMP (), including technological applications of quantum channel discrimination (e.g., quantum illumination Qill (); Qill2 () or quantum reading Qread (); Qread2 (); Qread3 (); Qread4 ()) where we are interested in increasing our ability to accept one specific quantum hypothesis.

Ii Brief review of asymmetric testing

ii.1 Basic formulation

In binary QHT we consider a quantum system which is prepared in some unknown quantum state , which can be  or . For instance we can imagine one party, say Alice, who prepares such a system. This system is then passed to Bob, who does not know which choice Alice has made. Thus, Bob must decide between the following two hypotheses

(1)
(2)

In order to discriminate between these two hypotheses, i.e., distinguish between the two states, Bob applies a quantum measurement, generally described by a positive operator valued measure (POVM). Without loss of generality, Bob can always reduce his measurement to be a dichotomic POVM with  Helstrom (). The outcome , with POVM operator , is associated to the null hypotheses , while the other outcome , with POVM operator , is associated with the alternative hypothesis .

Since the two quantum states  and are generally non-orthogonal, there is a non-zero error probability to confuse the two hypotheses. We can identify two different types of error: Type-I and type-II errors, with associated conditional error probabilities. By definition, the type-I error, also known as a ‘false-positive’, is where Bob accepts the alternative hypothesis when the null hypothesis  holds. We have a corresponding error probability expressed by

(3)

Then, the type-II error or ‘false-negative’ is where Bob accepts the null hypothesis when the true hypothesis is the alternative. This error occurs with conditional probability

(4)

Note that we can introduce other probabilities, but they are fully determined by and . For instance, we may also consider the ‘specificity’ or ‘true-negativity’ of the test which is the success probability of identifying the null hypothesis, i.e., which is simply given by . Similarly, we may also consider the ‘sensitivity’ or ‘true-positivity’ of the test which is the success probability of identifying the alternative hypothesis, i.e., .

The costs associated with the two types of error can be very different especially in the medical and histological settings. For instance, in a medical test,  is typically associated with no illness, while  with the presence of the disease. It is therefore clear that we would like to have tests where the false-negative probability (or rate)  is the lowest possible, so that ill patients are not diagnosed as healthy. For this reason, in a medical setting, hypothesis testing is almost always asymmetric, meaning that we aim to minimize one of the two conditional error probabilities.

ii.2 Multi-copy formulation

In general we can formulate the problem of QHT as an -copy discrimination problem QHT (); QHT2 (). This means that Alice has quantum systems which are prepared in two possible multi-copy states

(5)

These systems are passed to Bob who performs a collective measurement on them. As before, this general POVM can be chosen to be dichotomic with .

The error probabilities now depend on the number of copies . In particular, the probability of false positives is given by

(6)

and the probability of false negatives is

(7)

In the limit of a large number of copies , these probabilities go to zero exponentially, i.e., we have

(8)

where the coefficients

(9)
(10)

are called the ‘error-exponents’ or ‘rate limits’ QHB ().

Bob’s aim is to maximize the error exponent , so that the error probability of false negatives has the fastest exponential decay to zero. This must be done while controlling the rate of false positives. Here a well known result is the ‘quantum Stein lemma’ QHB () which connects  with the quantum relative entropy between the single-copy states  and . For a large number of copies , there is a dichotomic POVM such that the error probability of the false positives is bounded

(11)

and the error probability of false negatives goes to zero with error-exponent

(12)

More powerfully, we may use the notion of the QHB QHB (). For , there is a dichotomic POVM such that the error-exponent of false positives is lower-bounded by a positive parameter

(13)

and the error-exponent of false negatives satisfies

(14)

where is the QHB defined by

(15)

where

(16)

is the ‘s-overlap’ between the single-copy states and . Note that the quantum Hoeffding bound enforces a stronger constraint on false-positives, since these are bounded at the level of the error-exponent and not at the level of the error probability as happens for the quantum relative entropy bound.

Iii Asymmetric testing with pure states

Asymmetric testing becomes very simple when one of the states (or both) is pure. In this case, we can in fact relate the QHB to the quantum fidelity between the two states.

Let us start by considering the case where only one of the states is pure, e.g., . We can write JPAgae ()

(17)

where is the fidelity between and . Eq. (17) implies . By using the latter inequality in Eq. (15), we derive the fidelity-bound

(18)

This bound can be further simplified by explicitly performing the maximization with regard to the parameter . After a simple calculation we find

(19)

which depends on the comparison between the parameter and the fidelity of the two states.

More specifically, in the discrimination of two pure states, we find that the previous fidelity-bound becomes tight

(20)

In fact, for pure states and , and for any , we can write

(21)

Therefore we can replace in the QHB of Eq. (15), which implies Eq. (20NoteGG ().

Iv Asymmetric testing with Gaussian states

iv.1 Basics of bosonic systems and Gaussian states

A bosonic system of modes is a quantum system described by a tensor product Hilbert space and a vector of quadrature operators BraR1 (); BraR2 ()

(22)

These operators satisfy the vectorial commutation relations SPEDcomm ()

(23)

where is the symplectic form, defined as

(24)

Correspondingly, a real matrix is called ‘symplectic’ when it preserves by congruence, i.e., .

By definition, we say that a bosonic state is ‘Gaussian’ when its phase-space Wigner representation is Gaussian RMP (). In such a case, we can completely describe the state by means of its first- and second-order statistical moments. These are the mean value or displacement vector , and the covariance matrix (CM) with generic element

(25)

where denotes the anticommutator. The CM is a real symmetric matrix, which must satisfy the uncertainty principle RMP ()

(26)

An important tool in the manipulation of Gaussian states is Williamson’s theorem RMP (): For any CM , there is a symplectic matrix such that

(27)

where

(28)

The matrix  is the ‘Williamson form’ of , and the set is the ‘symplectic spectrum’ of . According to the uncertainty principle, each symplectic eigenvalue must satisfy the condition , with for all if and only if the Gaussian state is pure.

iv.2 Computation of the quantum Hoeffding bound

Our goal is to find a general recipe for the calculation of the QHB for Gaussian states. We start from the general formula in Eq. (15) involving the logarithm of the -overlap defined in Eq. (16). Given two -mode Gaussian states, and , we can write an explicit Gaussian formula for the -overlap in terms of their statistical moments (, ) and (, ). This is given by QCB2 (); JPAgae ()

(29)

where is the difference between the mean values, while and depends on the CMs and . In particular, introducing the two real functions

(30)
(31)

we can write the formulas

(32)

and

(33)

where and are the symplectic spectra of the two states, with and being the symplectic matrices which diagonalize the two CMs according to Williamson’s theorem, i.e.,

(34)

Substituting Eq. (29) into Eq. (15), corresponds to explicitly computing the logarithmic term , yielding

(35)

In particular for zero-mean Gaussian states we have and the previous expression simplifies to

(36)

iv.3 Other computable bounds

Note that computing the -overlap and its logarithmic form  could be difficult due to the presence of the symplectic matrices, and , in the term in Eq. (33). A possible solution is to compute an upper bound, known as the ‘Minkowski bound’, which is based on the Minkowski determinant inequality Bathia () and depends only on the two symplectic spectra QCB2 (). Specifically, we have , where

(37)

and

(38)

Another easy-to-compute upper bound is the ‘Young bound’ , which is based on Young’s inequality Young () and satisfies

(39)

where QCB2 ()

(40)

and

(41)

Taking the negative logarithm of Eq. (39), we can write the following inequality for the QHB

(42)

where

(43)
(44)

In the specific case where one of the two Gaussian states is pure, we can compute their fidelity and apply the upper bound given in Eqs. (18) and (19), which becomes tight when both states are pure [see Eq. (20)]. In particular, for two multimode Gaussian states and , we can easily write their fidelity in terms of the statistical moments JPAgae ()

(45)

where . As a result, we can use Eq. (19) with

(46)

V Discrimination of one-mode Gaussian states

In this section, we examine the case of one-mode Gaussian states. This means we fix in the previous formulas of Sec. IV, with matrices becoming , vectors becoming 2-dimensional, and symplectic spectra reducing to a single eigenvalue. For instance, the -overlap can be more simply computed using the expressions

(47)
(48)

In particular, here we shall derive the analytic formulas for the QHB for two important classes: Coherent states (in Sec. V.1) and thermal states (in Sec. V.2).

v.1 Asymmetric testing of coherent amplitudes

The expression of the QHB is greatly simplified in the case of one-mode coherent states and . Since both states are pure, the QHB is equal to the fidelity bound in Eq. (19), i.e., . Therefore, it is sufficient to compute the fidelity between the two coherent states, which is given by

(49)

so that , and we can write

(50)

Assuming that we impose a good control on the rate of false positives (so that ), then the error-exponent for the false negatives is simply given by . More explicitly, this corresponds to an asymptotic error rate

(51)

Note that, if we have poor control on the rate of false positives, i.e., , then the QHB is infinite. This means that the probability of false negatives goes to zero super-exponentially, i.e., more quickly than any decreasing exponential function.

v.2 Asymmetric testing of thermal noise

In this section we derive the QHB for one-mode thermal states and , with variances equal to and , respectively (in our notation, , where is the mean number of thermal photons). These Gaussian states have zero mean () and CMs in the Williamson form and (so that ). Thus, we can write

(52)

and derive

(53)

This is the -overlap to be used in the QHB of Eq. (15).

Given two arbitrary and , the maximization in Eq. (15) can be done numerically. The results are shown in Fig. 1 for thermal states with variances up to vacuum units (equivalent to mean thermal photon). From the figure we can see an asymmetry with respect to the bisector which is a consequence of the asymmetric nature of the hypothesis test. The bottom-right part of the figure is related to the minimum probability of confusing a nearly-vacuum state () with a thermal state having one average photon (). By contrast, the top-left part of the figure is related to the probability of confusing a thermal state having one average photon () with a nearly-vacuum state (). These probabilities are clearly different.

Figure 1: (Color online). We plot the QHB associated with the discrimination of two thermal states: as null hypothesis, and as alternative hypothesis. We consider low thermal variances and we have set for the false positives.

We are able to derive a simple analytical result when we compare a thermal state with the vacuum state. Let us start by considering the vacuum state to be the null hypothesis () while the thermal state is the alternative hypothesis (). In this specific case, we find

(54)

and we get

(55)

Since is a constant, the maximization of over corresponds to minimizing the function , whose minimum occurs at . As a result, we have

Since , we can write the QHB in terms of the mean number of thermal photons, i.e.,

(56)

This is the optimal error exponent for the asymptotic probability of false negatives, i.e., of confusing a thermal state with the vacuum state.

Let us now consider the thermal state to be the null hypothesis () while the vacuum state is the alternative hypothesis (). In this case, we derive

(57)

which leads to the following expression for the QHB

(58)

This is related to the minimum probability of confusing the vacuum state with a thermal state. Note that this is very different from Eq. (56).

Vi Discrimination of two-mode Gaussian states

In this section we consider two important classes of two-mode Gaussian states. The first is the class of Einstein-Podolsky-Rosen (EPR) states, also known as two-mode squeezed vacuum states. The second (broader) class is that of two-mode squeezed thermal (ST) states, for which the computation of the QHB is numerical.

vi.1 Asymmetric testing of EPR correlations

The expression of the QHB in the case of EPR states is easy to derive. Since EPR states are pure, the QHB is given by of Eq. (19). As a result, we need only to compute the fidelity between the two states.

An EPR state has zero mean and CM

(59)

with , is the identity matrix and

(60)

Given two EPR states with parameters and , their fidelity is computed via Eq. (45), yielding

(61)

where . After simple algebra, we find

(62)

to be used in Eq. (19).

vi.2 Squeezed thermal states

In this section we consider symmetric ST states , which are Gaussian states with zero mean and CM

(63)

where and  Bfide (); Bfide2 () (in particular, without loss of generality, we can assume ). These are called symmetric because they are invariant under permutation of the two modes asy ().

Note that, for , we have no correlations, and the ST state is a tensor-product of thermal states, i.e., . For the correlations are maximal, and the ST state becomes an EPR state, i.e., . Finally, for , we have maximal separable correlations. In other words, is the separable ST state with the strongest correlations (e.g., highest discord).

The symplectic decomposition of a symmetric ST state is known. From the CM of Eq. (63), one can check that the symplectic spectrum is degenerate and given by the single eigenvalue

(64)

The symplectic matrix which diagonalizes in Williamson form is given by

(65)

where

(66)

As a result, the s-overlap between two symmetric ST states, and , can be computed using the simplified formulas

(67)
(68)

where () is the degenerate eigenvalue of (), computed according to Eq. (64), and () is the corresponding diagonalizing symplectic matrix, computed according to Eqs. (65) and (66).

Let us start with simple cases involving the asymmetric testing of correlations with specific ST states. First we consider the asymmetric discrimination between the uncorrelated thermal state as null hypothesis and the correlated (but separable) ST state as alternative hypothesis. A false negative corresponds to concluding that there are no correlations where they are actually present Other (). It is straightforward to derive their degenerate symplectic eigenvalues which are simply and . Then, we have , while can be easily computed from Eqs. (65) and (66). By substituting these into Eqs. (67) and (68), we can compute the s-overlap and therefore the QHB via Eq. (15). The results are plotted in Fig. 2, for values of thermal variance up to (i.e., from zero to mean photon) and small values of the parameter , bounding the rate of false-positives. As expected, the QHB improves for decreasing and increasing .

Figure 2: (Color online). Asymmetric discrimination between the thermal state and the ST state with maximal separable correlations. We plot the QHB as a function of the thermal variance and the false-positive parameter . As we can see the QHB improves for lower and for higher .

Now let us consider the asymmetric discrimination between and the EPR state , i.e., the most correlated and entangled ST state Other (). Thanks to the simple symplectic decomposition of the EPR state (), we can further simplify the previous Eqs. (67)-(68) and write

(69)

with being given by Eq. (59). As before, we compute the QHB which is plotted in Fig. 3, for and . As expected the QHB improves for decreasing and increasing . Note a discontinuity identifying two regions, one where the QHB is finite, and the other where it is infinite (white region in the figure).

Figure 3: (Color online). Asymmetric discrimination between the thermal state and the EPR state . We plot the QHB as a function of the thermal variance and the false-positive parameter . The QHB improves for lower and for higher . In particular, there is a threshold value after which the QHB becomes infinite (white region).

In fact, by expanding the term in Eq. (15) for , that we find

(70)

where

(71)

For values of and such that , we find that the term diverges at the border, making the QHB infinite. For a given , this happens when

(72)

Finally, we consider the most general scenario in the asymmetric testing of correlations with ST states. In fact, we consider two generic ST states, and , with the same thermal noise but differing amounts of correlation. For this computation, we use Eqs. (64)-(66) with or , to be replaced in Eqs. (67)-(68), therefore deriving the s-overlap and the QHB. At small thermal variance () and for the numerical value , we plot the QHB as a function of the correlation parameters and . As we can see from Fig. 4, the QHB is not symmetric with respect to the bisector (where it is zero) and increases away from this line.

Figure 4: (Color online). Asymmetric discrimination between two ST states with the same thermal variance () but different correlations and . Setting , we plot the QHB as a function of and . We can see that the QHB increases orthogonally to the bisector .

Vii Conclusion

In this work we have considered the problem of asymmetric quantum hypothesis testing by adopting the recently-developed tool of the quantum Hoeffding bound (QHB). After a brief review of these notions, we have shown how the QHB can be simplified in some cases (pure states) and estimated using other easier-to-compute bounds based on simple algebraic inequalities.

In particular, we have applied the theory of asymmetric testing to multimode Gaussian states, providing a general recipe for the computation of the QHB in the Gaussian setting. Using this recipe, we have found analytic formulas and shown numerical results for important classes of one-mode and two-mode Gaussian states. In particular, we have studied the behavior of the QHB in the low energy regime, i.e., considering Gaussian states with a small average number of photons.

Our results could be exploited in protocols of quantum information with continuous variables. In particular, they could be useful for reformulating Gaussian schemes of quantum state discrimination and quantum channel discrimination in such a way as to give more importance to one of the quantum hypotheses. This asymmetric approach could be the most suitable in the development of quantum technology for medical applications.

Acknowledgments

This work has been supported by EPSRC (EP/J00796X/1). G.S. has been supported by an EPSRC DTA grant. The authors thank C. Ottaviani and S. Pirandola for enlightening discussions.

References

  • (1) M. M. Wilde, Quantum Information Theory (Cambridge University Press, Cambridge, 2013).
  • (2) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press 2000).
  • (3) C. Silberhorn, T. C. Ralph, N. Lutkenhaus, and G. Leuchs, Phys. Rev. Lett. 89, 167901 (2002).
  • (4) A. M. Lance, T. Symul, V. Sharma, C. Weedbrook, T. C. Ralph, and P. K. Lam, Phys. Rev. Lett. 95, 180503 (2005).
  • (5) C. Weedbrook, S. Pirandola, R. Garcia-Patron, N. J. Cerf, T. C. Ralph, J. H. Shapiro, and S. Lloyd, Rev. Mod. Phys. 84, 621 (2012).
  • (6) C. W. Helstrom, Quantum Detection and Estimation Theory, Mathematics in Science and Engineering, Vol. 123 (Academic Press, New York, 1976).
  • (7) A. Chefles, Contemp. Phys. 41, 401 (2000).
  • (8) S. M. Barnett and S. Croke, Advances in Optics and Photonics 1, 238-278 (2009).
  • (9) K. M. R. Audenaert, J. Calsamiglia, L. Masanes, R. Munoz-Tapia, A. Acın, E. Bagan, and F. Verstraete, Phys. Rev. Lett. 98, 160501 (2007).
  • (10) S. Pirandola, and S. Lloyd, Phys. Rev. A 78, 012331 (2008).
  • (11) K. M. R. Audenaert, M. Nussbaum, A. Szkola, and F. Verstraete, Commun. Math. Phys. 279, 251 (2008).
  • (12) R. Jozsa, J. Mod. Opt. 41, 2315 (1994).
  • (13) S.-H. Tan, B. I. Erkmen, V. Giovannetti, S. Guha, S. Lloyd, L. Maccone, S. Pirandola, and J. H. Shapiro, Phys. Rev. Lett. 101, 253601 (2008).
  • (14) S. Lloyd, Science 321, 1463 (2008).
  • (15) S. Pirandola, Phys. Rev. Lett. 106, 090504 (2011).
  • (16) S. Pirandola, C. Lupo, V. Giovannetti, S. Mancini, and S. L. Braunstein, New J. Phys. 13, 113012 (2011).
  • (17) G. Spedalieri, C. Lupo, S. Mancini, S. L. Braunstein, and S. Pirandola, Phys. Rev. A 86, 012315 (2012).
  • (18) C. Lupo, S. Pirandola, V. Giovannetti, and S. Mancini, Phys. Rev. A 87, 062310 (2013).
  • (19) G. Spedalieri, C. Weedbrook, and S. Pirandola, J. Phys. A: Math. Theor. 46, 025304 (2013).
  • (20) We can write , where can be neglected and
  • (21) S. L. Braunstein, and P. van Loock, Rev. Mod. Phys. 77, 513 (2005).
  • (22) S. L. Braunstein, and A. K. Pati, Quantum Information with Continuous Variables (Kluwer Academic, Dordrecht, 2003).
  • (23) More generally, for any two vectorial operators and , we can express their commutation relations in the compact form .
  • (24) R. Bhatia, Matrix Analysis (Springer-Verlag, New York, 1997).
  • (25) W. H. Young, Proc. R. Soc. London, Ser. A 87, 331 (1912).
  • (26) S. Pirandola, A. Serafini, and S. Lloyd, Phys. Rev. A 79, 052327 (2009).
  • (27) S. Pirandola, New J. Phys. 15, 113046 (2013).
  • (28) Extension to asymmetric ST states is only technical.
  • (29) For brevity we do not consider the other case where the ST state is the null hypothesis and the thermal state is the alternative hypothesis. This case is included in the our final analysis for generic ST states.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
267155
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description