Experimental Realization of Optimal Noise Estimation for a General Pauli Channel

Experimental Realization of Optimal Noise Estimation for a General Pauli Channel

A. Chiuri Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Roma, Italy Istituto Nazionale di Ottica (INO-CNR), Largo Enrico Fermi 6, I-50125 Firenze, Italy    V. Rosati Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Roma, Italy    G. Vallone Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Roma, Italy Department of Information Engineering, University of Padova, I-35131 Padova, Italy    S. Pádua Departamento de Física, Universidade Federal de Minas Gerais. Caixa Postal 702, Belo Horizonte, MG 30123-970, Brazil    H. Imai Dipartimento di Fisica “A. Volta” and INFN-Sezione di Pavia, via Bassi 6, 27100 Pavia, Italy    S. Giacomini Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Roma, Italy    C. Macchiavello Dipartimento di Fisica “A. Volta” and INFN-Sezione di Pavia, via Bassi 6, 27100 Pavia, Italy    P. Mataloni Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Roma, Italy Istituto Nazionale di Ottica (INO-CNR), Largo Enrico Fermi 6, I-50125 Firenze, Italy
August 26, 2019

We present the experimental realization of the optimal estimation protocol for a Pauli noisy channel. The method is based on the generation of 2-qubit Bell states and the introduction of quantum noise in a controlled way on one of the state subsystems. The efficiency of the optimal estimation, achieved by a Bell measurement, is shown to outperform quantum process tomography.

Figure 1: (a) Experimental setup. Photons A and B are spatially and temporally superimposed on a symmetric beam splitter (BS). The optical path delay allows us to vary the arrival time of the photons on the BS. Photons are collected by using an integrated system, composed by a GRIN lens (GL) and a single mode fiber, and then detected by single photon counters. The same setup allows us to perform the ancillary assisted quantum process tomography (AAQPT) after removing the BS. Quantum state tomography (QST) jam01pra on the output state is performed by using quarter-wave plates (QWPs), half-wave plates (HWPs), and polarizing beam splitters (PBSs). (b) Scheme of the implemented Pauli channel. , , represent the time intervals of , or activation. Both , , and the repetition time can be varied by a remote control.

Introduction.- Quantum noise is unavoidably present in any realistic implementation of quantum tasks, ranging from quantum communication protocols cira99pra to quantum information processing devices and quantum metrology giov04science ; giov06prl . The performance and the optimization of quantum tasks quite often depend on the level of noise which is present in the physical realization considered. It is therefore of great interest to develop experimental methods to estimate the level of noise in the system under examination as precisely as possible. Quantum process tomography (QPT) nielsen , which has already been implemented in various experimental realizations qpt-exp ; obri04prl , represents a well-known method to identify an unknown noise, but it lacks the notion of efficiency. In many realistic scenarios, however, some a priori information on the kind of noise is available and therefore the problem of measuring it is equivalent to estimate few noise parameters in the most efficient way. This is the context of quantum channel estimation qc-est , which is based on quantum state estimation theory, a merge of statistical estimation theory and quantum physics pioneered by Helstrom and Holevo Helstrom ; Holevo . The aim of this Letter is to provide the first genuine experimental application of quantum channel estimation theory. The experimental realization presented here is based on a quantum optical setup, but it opens new perspectives of applications to a great variety of physical scenarios and quantum technologies, from atomic to solid state systems.

A quantum channel estimation problem is generally formulated as follows. We need to estimate the true value of the -dimensional parameter for a given smooth parametric family of noisy quantum channels. The estimation scheme is twofold: first we prepare a quantum system described by the density operator as input to the channel , then we perform some quantum state measurement on the output state in order to estimate the channel parameters in the most efficient way. Thus, the problem is to seek an optimal input for the channel and an optimal measurement on the output state . The notion of optimality is here based on the minimization of the covariance matrix


where , and denote the estimated and the true values for the -th component of the noise parameter respectively, denotes the expectation with respect to the measurement procedure employed, and is a projective measurement operator.

In this work we present an experimental implementation of an optimal quantum channel estimation scheme for a qubit Pauli channel (PC). The action of such a family of channels on the density operator of a qubit can be described as nielsen


where is the identity operator, {} () are the three Pauli operators respectively, and {} represent the corresponding probabilities (). The family of Pauli channels represents a wide class of noise processes, that includes several physically relevant cases such as the depolarising channel, which will be considered in the following, the dephasing and the bit-flip channels.

The optimal channel estimation scheme is achieved as follows opt-qce . The optimal input state is represented by a Bell state for two qubits, for example the singlet state , where only one of the qubits is affected by the noisy channel while the other one is left untouched. The optimal measurement consists of a Bell measurement on the two qubits at the channel output, namely the projective measurement . The outcome probabilities then provide an optimal estimation of the channel parameters .

As mentioned above, this scheme is optimized by minimizing the covariance matrix of the estimation error (1). According to the quantum Cramér-Rao theorem Holevo , the minimum covariance matrix in this case is given by opt-qce


where is the quantum Fisher information matrix Holevo of a maximally entangled input state . We want to point out that this scheme is optimal for any number of input qubits. Actually, no additional entanglement among the input qubits and no collective measurements at the output can increase the efficiency of the present scheme opt-qce . Moreover, it can be also straightforwardly generalized to estimate any general noise process of the form (2), where the operators are replaced by any set of unitary operators such that Tr. The same scheme can also straightforwardly extended to estimate any generalized Pauli channel for quantum systems in arbitrary finite dimension opt-qce .

We will now present the experimental implementation of this optimal estimation scheme for a quantum optical setup, where the state of the two qubits is represented by polarization states of two photons and the action of the Pauli channel is introduced in a controlled way by employing liquid crystal retarders, as explained in the following. The method has been first applied to estimate a general Pauli channel, with independent values of the probabilities . Then it has been applied to a depolarizing channel (DC), namely the case of isotropic noise, with , where the parameter completely specifies the channel itself, and the minimum variance (Eq. (1) for the one-dimensional case) is given by . In this case the procedure simplifies and, in the following, we show that only two projective measurements, , are needed.

Figure 2: Experimental probabilities of measuring the four Bell states obtained for four different cases of anisotropic noise. The black boxes report the corresponding theoretical values. a) Identity: noiseless channel. b) : only one Pauli matrix, is acting on the state . c) Partially anisotropic DC: and operate for the same time interval, in fact the probabilities of measuring the states and are equal. d) Totally anisotropic DC: each Pauli operator operates for a different time interval.
Figure 3: Noise parameter estimation for the DC case. (a) Measured values of vs by implementing the projective measurements. Continuous red line corresponds to the theoretical behavior. Inset: experimental values of the standard deviations for the optimal protocol implemented by the projective measurements. They are obtained by propagating the Poissonian uncertainties. The solid line represents the expected theoretical behavior. b): Experimental probabilities associated to the experimental matrix vs . Values of are obtained by maximizing the fidelity between theoretical and experimental matrix . Error bars are calculated by considering the Poissonian uncertainty associated to the coincidence counts, and simulating different matrices of the process, obtaining, in this way, different values of . Inset: experimental values of the standard deviations for the AAQPT. These have been calculated by a simulation based on the poissonian uncertainty associated to the coincidence counts. The solid line represents the optimal bound. (c),(d) Experimental (left side) and theoretical (right side) matrices for and .

Experimental scheme.- Different techniques have been exploited to experimentally implement a PC acting on a single qubit state jinprl09 ; ricprl04 ; karjosab04 ; eise11pra . The optimal noise estimation protocol, proposed in this work, was implemented by the interferometric scheme shown in Fig.1a). Precisely, a two-photon entangled source cinpra04 generates the two-qubit singlet state , where two qubits are encoded in the polarization degree of freedom, with () referring to the horizontal (vertical) polarization of photons and . In our setup, the single qubit noisy channel is operating only on one of the two entangled particles (i.e. photon A). The general Pauli channel (PC) consists of a sequence of liquid crystal retarders (LC1 and LC2) in the path of photon . The LCs act as phase retarders, with the relative phase between the ordinary and extraordinary radiation components depending on the applied voltage . Precisely, and (Fig.1b) correspond to the case of LCs operating as half-wave plate (HWP) and as the identity operator, respectively. The LC1 and LC2 optical axes are set at and with respect to the V-polarization. Then, when the voltage is applied, the LC1 (LC2) acts as a () on the single qubit. We were able to switch between and in a controlled way and independently for both LC1 and LC2. The simultaneous application of on both LC1 and LC2 corresponds to the operation. We could also adjust the temporal delay between the intervals in which the voltage is applied to the two retarders. We define , , respectively as the activation time of the operators , or and is the period of the LCs activation cycle, as shown in Fig.1b).

Experimental implementation of the Pauli channel (anisotropic noise).- A general PC was generated by varying the four time intervals , , and . The intervals are related to the probabilities (), introduced in Eq. (2), by the following expression: . The probability of the identity operator is given by (with ).

To obtain the probabilities associated to the four projection operators , we measured the coincidence counts between the two outputs of the BS. In fact, these probabilities are related to the interference visibility measured by the interferometer in Fig.1a). The half-wave plate (HWP) and quarter-waveplate (QWP) of Fig.1a) were used to project the noisy state onto the four different Bell states.

Different configurations of the noisy channel were investigated by implementing the optimal noise protocol estimation for each configuration. A summary of four relevant experimental results, corresponding to different probabilities associated to the Bell states, are given in Fig.2.

In the measurements shown in Fig.2, case (a) correspond to a noiseless channel (identity transformation) while the cases b), c) and d), correspond to different complete noisy channels with (i.e. we set ). For each process, the first column shows the relative weights between the Pauli operators acting in the channel. From these values it is possible to calculate the theoretical ones. For instance, let us consider the process d) where the , and act respectively for , and . The expected values of are, for this process, , , and . The slight disagreement between the expected theoretical values and the experimental measured ones are mainly due to the finite rise and decay times of the electrical signal driving the LC devices.

We have implemented the protocol by using always the same input state and projecting it on the Bell basis. It is worth noting that this is totally equivalent to entering the PC with the four Bell states and to projecting them into the state.

Experimental implementation of the depolarizing channel (isotropic noise).- The condition corresponds to the depolarizing channel, with the three Pauli operators acting on the single qubit with the same probability . This parameter was changed by fixing the times and varying the period . The optimal protocol to estimate the value of was realized by using the Bell state , as mentioned above.

The DC was activated on photon . In this case the projective measurement , consisting of just two projectors, is sufficient to optimally estimate and has been performed for several noise degrees. For each level of noise, we estimated the channel parameter as where are the coincidences between the two outputs of the BS in interference condition and is the number of events in which the two photons are detected on the same BS output side. was estimated by knowing the amount of coincidences out of interference. The typical peak interference measured for the state as a function of the path delay is shown in Fig.4 of the Supplemental Material. In Fig.3a) we report the experimental values corresponding to the different values of . In the corresponding inset we show the errors evaluated by propagating the and Poissonian errors. They are in good agreement with the expected theoretical behavior.

Ancillary assisted quantum process tomography.- The experimental results, just discussed for the optimal estimation of the depolarizing channel, have been compared with the probability values of which can be obtained by exploiting the Ancillary assisted quantum process tomography (AAQPT) DL ; Alteprl03 ; mohs08pra ; karjosab04 . The action of a generic channel operating on a single qubit can be written as , where the matrix characterizes completely the process.

AAQPT is based on the following procedure: i) prepare a two-qubit maximally entangled state and reconstruct it by Quantum state tomography (QST) jam01pra ; ii) send one of the two entangled qubits through the channel ; iii) reconstruct the output two-qubit state by QST and obtain, in this way, the matrix from the two-qubit output density matrix. For a DC, the matrix is expressed as nielsen


We implemented the AAQPT algorithm by injecting the state into the DC and we reconstructed by QST the density matrices of the input and output states for several noise degrees [see Fig.1a)]. We obtained the experimental matrix for different values of and, for each value of , we found the parameter maximizing the fidelity between the experimental and the theoretical process matrices. The experimental results are shown in Fig.3b). Even in this case the theoretical behavior is fully satisfied. However, comparing these results with those obtained by the optimal protocol, we observe that the latter leads to the same results, but with a much lower number of measurements. In fact, in this case, only the two projections are needed while, to implement the AAQPT algorithm, 16 measurements are necessary. Moreover, by adopting our experimental setup we were able to demonstrate that the value of and the DC action do not depend on the input state. In fact the AAQPT was realized with all the four Bell states entering the DC, obtaining the same results of those shown in Fig. 3b).

It is worth noting that, even if the AAQPT gives a more complete information on the process compared to the implemented optimal protocol, the latter allows us to achieve a more accurate value of . The inset in Fig.3(a) shows that, for the optimal protocol, the measured standard deviation reaches the lower bound given by [i.e. the square root of Eq. (1) for the one-dimensional case divided by ], where is the dimension of the sample used to evaluate , thus demonstrating experimentally the attainability of the Cramér-Rao bound. We show in the inset in Fig.3(b) the standard deviations, well above the optimal bound, obtained with the AAQPT (see the Supplemental Material for details about the numerical estimation). The lower optimal bound represented by the black curve is below the experimental data, demonstrating that AAQPT is far away from the optimal estimation protocol presented in this work.

Conclusion.- An optimal protocol allowing the most efficient estimation of a noisy Pauli channel has been experimentally implemented in this work. The action of the noisy channel was introduced on one qubit of a maximally entangled pair in a controlled way. The efficiency of this method has been compared to the one achieved by quantum process tomography, demonstrating that the optimal protocol allows us to achieve the theoretical lower bound for the errors and to perform the estimate of the noisy channel with a lower number of measurements. This method can be profitably applied when some knowledge on the noise process is available and can be successfully implemented in quantum-enhanced technologies involving the management of decoherence.

This work was supported by EU-Projects CHISTERA-QUASAR and CORNER, by the CNR project ”Hilbert“, and by the FARI project 2010 of Sapienza Università di Roma. GV was partially supported by the Strategic-Research-Project QUINTET of the Department of Information Engineering, University of Padova and the Strategic-Research-Project QUANTUMFUTURE of the University of Padova. SP was supported by CNPq, FAPEMIG, INCT-quantum information and Italian-Brazilian Contract CNR-CNPq (Quantum information in a high dimension Hilbert space).


  • (1) J. I. Cirac et al., Phys. Rev. A 59, 4249 (1999).
  • (2) V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 1330 (2004).
  • (3) V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. Lett. 96, 010401 (2006).
  • (4) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge University Press, Cambridge, U.K., 2000).
  • (5) I. Bongioanni et al., Phys. Rev. A 82, 042307 (2010).
  • (6) J. L. O’Brien et al., Phys. Rev. Lett. 93, 080502 (2004).
  • (7) A. Fujiwara, Phys. Rev. A 63, 042304 (2001).
  • (8) C.W. Helstrom, Quantum Detection and Estimation Theory, (Academic Press: New York, 1976).
  • (9) A. Fujiwara and H. Imai, J. Phys. A 36, 8093 (2003).
  • (10) A. S. Holevo, Probabilistic and Statistical Aspects of Quantum Theory, (Amsterdam: North-Holland, 1982).
  • (11) C. Cinelli et al., Phys.Rev.A 70, 022321 (2004).
  • (12) Jin-Shi Xuet al., Phys.Rev.Lett. 103,240502 (2009).
  • (13) M. Ricci et al., Phys. Rev. Lett. 93, 170501 (2004).
  • (14) M. Karpinski, et al., JOSA B 25, 668 (2008).
  • (15) A. Shaham and H. S. Eisenberg, Phys.Rev.A 83, 022303 (2011).
  • (16) G. M. D’Ariano and P. Lo Presti, Phys. Rev. Lett. 86, 4195 (2001).
  • (17) M. Mohseni et al., Phys. Rev. A 77,032322 (2008).
  • (18) J. B. Altepeter, et al., Phys. Rev. Lett. 90, 193601 (2003).
  • (19) D. F. V. James et al., Phys.Rev.A 64 ,052312 (2001).

I Supplementary Information

Ancillary Assisted Quantum Process Tomography (AAQPT) standard deviations, shown in the inset in Fig.3b) of the main text, have been calculated by a simulation realized by MATHEMATICA 5.0. It was based on the following procedure:

  • a poissonian uncertainty was associated to the coincidence counts

  • hundred different matrices of the process were simulated

  • the fidelity between these matrices and the theoretical ones, reported in the main text, was calculated

  • the routine NMAXIMIZE allowed to find numerically the value of maximizing the Fidelity for each simulated process

  • this sample, composed of hundred values of , was used to calculate both average and standard deviation.

The experimental values of the inset in Fig.3b) correspond to a sample of 1600 coincidences per second. It can be evaluated that a number of coincidences per second, which is nearly two orders of magnitude larger, is needed to approach the standard deviation value obtained by the optimal protocol (i.e. the inset in Fig.3a) of the main text).

Figure 4: Dip (peak) of the measured coincidence counts as a function of the optical path delay for a state () entering the BS in absence of noise.
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
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

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 description