Heisenberg limited quantum metrology under the effect of dephasing

Heisenberg limited quantum metrology under the effect of dephasing

Yuichiro Matsuzaki matsuzaki.yuichiro@lab.ntt.co.jp NTT Basic Research Laboratories, NTT Corporation, 3-1 Morinosato-Wakamiya, Atsugi, Kanagawa, 243-0198, Japan. NTT Theoretical Quantum Physics Center, NTT Corporation, 3-1 Morinosato-Wakamiya, Atsugi, Kanagawa 243-0198, Japan.    Simon Benjamin Department of Materials, University of Oxford, OX1 3PH, United Kingdom    Shojun Nakayama National Institute of Informatics, 2-1-2 Hitotsubashi, Chiyoda-ku, Tokyo 101-8430, Japan.    Shiro Saito NTT Basic Research Laboratories, NTT Corporation, 3-1 Morinosato-Wakamiya, Atsugi, Kanagawa, 243-0198, Japan.    William J. Munro NTT Basic Research Laboratories, NTT Corporation, 3-1 Morinosato-Wakamiya, Atsugi, Kanagawa, 243-0198, Japan. NTT Theoretical Quantum Physics Center, NTT Corporation, 3-1 Morinosato-Wakamiya, Atsugi, Kanagawa 243-0198, Japan. National Institute of Informatics, 2-1-2 Hitotsubashi, Chiyoda-ku, Tokyo 101-8430, Japan.

Quantum sensors have the potential to outperform their classical counterparts. For classical sensing, the uncertainty of the estimation of the target fields decreases with measurement time . On the other hand, by using quantum resources, one can decrease the uncertainty with time to - the Heisenberg limit. However as quantum states are susceptible to dephasing, it has been not clear whether one can achieve the Heisenberg limit for a measurement time scaling longer than the coherence time of those states. Here, we propose a scheme that estimates the amplitude of globally applied fields with Heisenberg limited sensitivity in principle for an arbitrary time scale under the effect of dephasing. We use one way quantum computing based teleportation between qubits to prevent the correlation between the quantum state and its local environment from increasing, and have shown that such a teleportation protocol can suppress the local dephasing while the information from the target fields keeps growing. Our methods has the potential to realize a quantum sensor with a sensitivity far beyond that of any classical sensors.


It is well known that two-level systems are attractive candidates to realize ultrasensitive sensors as the frequency of the qubit can be shifted by a coupling to target fields. Such a frequency shift induces a relative phase between the qubits basis states which can be simply measured in a Ramsey type experiment. This method has been used to measure magnetic fields, electric fields, temperature, and even gravity waves Budker and Romalis (2007); Balasubramanian and et al (2008); Maze and et al (2008); Degen et al. (2016); Takano et al. (2016). With the typical classical sensor measurement device (including SQUID’s Simon (1999), Hall sensors Chang et al. (1992), and force sensors Poggio and Degen (2010)), the uncertainty in the estimation of the target fields decreases only slowly as with measurement time . This scaling is known as the shot noise limit Huelga et al. (1997). Quantum and especially qubit-based sensors can in principle decrease the uncertainty with time to - the Heisenberg limit Said et al. (2011); Higgins et al. (2007); Waldherr et al. (2012). This can be achieved using the coherence properties of the qubit and several experimental demonstrations using quantum feedback have shown sub-shot noise sensitivity measuring the amplitude of target fields Said et al. (2011); Higgins et al. (2007); Waldherr et al. (2012). However as quantum states are fragile to decoherence, it has generally been considered that such sub-shot noise scaling can only be realized if the measurement time is much smaller than the coherence time of the qubit Demkowicz-Dobrzanski et al. (2012); Said et al. (2011). Recently several approaches using quantum error correction Gottesman (2009) and dynamical decoupling Viola et al. (1999); Taylor et al. (2008); De Lange et al. (2011) have been proposed to circumvent this limitation. Using quantum error correction, one can measure an amplitude of the target field with the Heisenberg limited sensitivity under the effect of specific decoherence such as bit flip errors Dür et al. (2014); Herrera-Martí et al. (2015); Arrad et al. (2014); Kessler et al. (2014); Cohen et al. (2016); Unden et al. (2016); Sisi and et al (2017), while dynamical decoupling makes it possible to estimate the frequency of time-oscillating fields with sensitivity beyond the shot-noise limit on time scale longer than the coherence time Schmitt and et al (2017); Boss et al. (2017). However, there is no currently known metrological scheme to achieve Heisenberg limited sensitivity when measuring the amplitude of the target fields with dephasing.

In this letter we propose a scheme to realize a Heisenberg-limited sensing of the amplitude of target fields under the effect of dephasing. We use a similar concept to the quantum Zeno effect (QZE) Misra and Sudarshan (1977); Itano et al. (1990); Facchi et al. (2001). For shorter time scales than the correlation time of the environment , the interaction with the environment induces a quadratic decay rate that is much slower than an exponential decay Nakazato et al. (1996). Hence frequent measurements to reset the correlation with the environment can keep the states in the initial quadratic decay region, which suppresses the decoherence Misra and Sudarshan (1977); Itano et al. (1990); Facchi et al. (2001). However, if we naively apply the QZE in quantum metrology, the frequent measurements freezes all dynamics so that the quantum states cannot acquire any information from the target fields. Instead, we use quantum teleportation (QT) based on the concepts from one-way quantum computation Raussendorf and Briegel (2001); Barjaktarevic et al. (2005); Silva et al. (2007); Olmschenk et al. (2009); Baur et al. (2012) to reset the correlation between the system and environment Averin et al. (2016). If we transfer the quantum states to a new site, we can prevent the correlation between the system and environment in the previous site from increasing, and the quantum state are then only affected by a slow quadratic decay due to the local environment in the new site. This noise suppression has been proposed and demonstrated by superconducting qubits Averin et al. (2016). The crucial idea in this paper is to use this one-qubit teleportation-based noise suppression for quantum metrology. Interestingly, although the QT protocol eliminates the deteriorate effect due to the dephasing from the local environment, we can accumulate the phase information from the global target fields during this protocol. We have shown that, as long as nearly perfect QT is available, we can achieve the Heisenberg limit sensing with dephasing. Moreover, we have found that, even when QT is moderately noisy, the sensitivity of our protocol can be much better than that of the standard Ramsey measurement.

Our system and environment can be described by a Hamiltonian of the form Hornberger (2009) where () denotes the system (environmental) Hamiltonian while denotes an interaction between the system and environment. Here is the usual pauli Z operator of the -th qubit with frequency , while and denote the environmental operator at that -th site. () denotes an identity operator for the system (environment). Transforming to an interaction picture, we have where . The initial state is given as where we further assume is in thermal equilibrium such that and our noise is non-biased such that for all . Since the initial state is separable, we consider only the first site by tracing out the other sites. Solving Schrodinger equation, we obtain

using a second order perturbation expansion in where denotes a time Hornberger (2009). Tracing out the environment, we obtain

where we define the correlation function of the environment as . If we are interested in a time scale much shorter than the correlation time of the environment, we can approximate the correlation function as . For most of the solid state systems, the correlation time is much longer than the coherence time of the qubit De Lange et al. (2010); Yoshihara et al. (2006); Kakuyanagi et al. (2007); Kondo et al. (2016), and so this condition is readily satisfied for many systems. In this case where denotes an error rate for and denotes a unitary operator at a site . Since the error rate has a quadratic form against the time , the decoherence effect is negligible for short time scales , which has been discussed in the field of the QZE Misra and Sudarshan (1977); Itano et al. (1990); Facchi et al. (2001). On the other hand, if we consider longer time scales with the same environment, error accumulation will destroy the quantum coherence of the qubit.

Let us now describe the noise suppression technique using QT. It begins with allowing free evolution of the qubit for a time (where is the total time and is the number of time QT to be performed). After this, QT sends to a different site . The quantum states then starts interacting with a new local environment described by a density matrix . The error rate will be suppressed due to the quadratic decay Averin et al. (2016). Performing QT times (each time to a fresh qubit) yields at site . For large this approaches the pure state meaning the QT approach can suppress the local dephasing.

Our general approach described above has used a perturbative analysis typically valid only for a short time scale. However we need to examine the dynamics of our system for arbitrary time scales, and so we will consider a more specific noise model given by during the evolution for a time (with giving the dephasing rate). This model is consistent with the general short time scale noise model described above by choosing . Typical dephasing models Palma et al. (1996); De Lange et al. (2010); Kakuyanagi et al. (2007); Yoshihara et al. (2006) show this behavior if the correlation time of the environment is much longer than the dephasing time.

Figure 1: Schematic illustration of qubits in a ring structure to measure globally applied fields at the Heisenberg limit. Half of the qubits contain an information of the target fields as a probe while the remaining half are used as an ancilla’s for one qubit teleportation. By performing a controlled phase gate between the probe qubit and ancillary qubit, single qubit measurements on the probe qubit with feedforward operations depending on the measurement results, one can teleport a quantum state from the original site to the right neighboring site Browne and Briegel (2006); Raussendorf and Briegel (2001). After the teleportation, the measured qubit becomes the new ancilla which we initialize into . By implementing the teleportation frequently, dephasing from the local environment can be suppressed while quantum states keep accumulating a phase information from the target global fields.

Now let us turn our focus to using the QT scheme to enable quantum metrology with Heisenberg-limited sensitivity. Consider the situation in which the qubit frequency is shifted depending on the amplitude of the target fields, and so a measurement of the qubit’s frequency shift allows us to infer the amplitude of the target field. Such a qubit frequency shift is estimated from a relative phase between quantum states. The key idea is to use the QT in a ring arrangement where a qubit has a tunable interaction with it’s nearest neighbor qubits, as illustrated in Fig. 1. Half of the qubits are used to probe the target fields while the remaining qubits are used as ancilla for QT (probe qubits are located between two ancillary qubits). The QT is enacted by performing a control-phase gate between a probe qubit and an ancilla qubit, followed by a measurement on the probe qubit (and single qubit corrections depending on the measurement result). This QT approach has been widely used in one-way quantum computation Raussendorf and Briegel (2001); Browne and Briegel (2006).

Our scheme to measure the amplitude of target fields is as follows: First, we prepare a probe state of located at the site (). Second we then let the state evolve for a time and then teleport (using QT) the state of the probe qubit to the next site using the ancillary qubit (we assume that our gate operations are much faster than ). Third we then in step 3 repeat the second step times while in the fourth step we let this state evolves for a time , and readout the state by measuring . Finally, we repeat these steps times during the measurement time where is the repetition number. When is even, the density matrix before the readout is described as where . Here, can be interpreted as an QT improved coherence time. If is odd, we obtain the same density matrix for the probe qubits at a site . In the situation that the target field to be sensed is weak () we can estimate the uncertainty in our estimator as


where . Setting we obtain for large . The scaling shows we can achieve the Heisenberg limit.

Let us now analyze how imperfect QT affects the performance of our sensing scheme. In our previous ideal QT analysis, the uncertainty of the estimation monotonically decreases as the number of the QT increases. However in realistic situations, errors caused during the QT operations will limit our achievable sensitivity and there will be an optimal number of QT’s that we can perform. We will consider a simple error model, a depolarization channel that makes the state completely mixed with a probability after the QT operation. In such a case, after the imperfect QT is performed on say the th site, the state evolves to . The uncertainty of the QT based metrology estimation for our weak target field is then . For , we obtain when and . This allows us to achieve Heisenberg limited sensitivity.

In quantum metrology one typically considers the scaling law in the limit of long , which we will now discuss. We can minimize this uncertainty by setting as long as is satisfied. In such a case we obtain which for gives the standard Ramsey uncertainty Huelga et al. (1997) where we replace with (because the standard Ramsey scheme can utilize every qubit to probe the target fields without ancillary qubits). For . we can treat as a continuous variable, and we can analytically minimize the uncertainty as for where we choose . In this case, we have a constant factor improvement over the standard Ramsey scheme for a longer . Actually, as long as , our scheme is better than that standard Ramsey scheme (). Also, if we have , we obtain . So our sensor has an advantage with finite errors by the QT.

There are of course other sources of decoherence that can not be suppressed by the QT protocol. For instance, If our quantum systems are affected by high-frequency noise with a short correlation time, the decay is not quadratic in nature but more exponential like. Energy relaxation in a high temperature environment is known to induce such noise Gardiner and Zoller (2004). Consider the situation in which an initial state evolves under the effect of both low-frequency dephasing and high-frequency noise for a time . In such a case where denotes the decay rate associated with the high-frequency noise ( gives the same noise model we used previously). It is then straightforward to calculate the uncertainty of the estimation under the effect of this noise with imperfect QT as . Choosing (where ), we minimize this with respect to time as

A numerical minimization of the uncertainty with as can be done. In Fig. 2, we plot versus and where is the uncertainty for the standard Ramsey scheme. Our plots shows that our scheme has a better performance than the standard Ramsey scheme for and .

Figure 2: Performance of our teleportation based scheme against the standard Ramsey scheme. We plot against and . Our scheme has a better performance than the standard Ramsey scheme for and .

A natural question we have not addressed so far is whether entanglement improves our teleportation-based sensing. We could for instance create a GHZ state composed of qubits. For a given qubits, we create GHZ states with this size, and the number of the GHZ states is . By letting the GHZ states evolves with low-frequency dephasing for a time where we ignore the high-frequency decoherence for simplicity, we have

for where denotes the dephasing rate for a single qubit. To readout the GHZ states, we measure a projection operator defined by where . We consider an imperfect QT as follows. If we teleportate a state of from sites to sites, we obtain a state of where denotes the error rate on a single qubit, denotes the ideal state (that we could obtain by a perfect QT), and denotes a decohered state. Here, we assume that any error on a single qubit makes the probe state into (The corresponds to our previous case where we used separable states). It is straightforward in this GHZ entangled situation to estimate our sensitivity as

where and . By setting , we reproduce the results discussed in Matsuzaki et al. (2011); Chin et al. (2012) as an entanglement based sensor with low-frequency dephasing. We can minimize this uncertainty with to obtain . For the ideal perfect QT (), we achieve the Heisenberg limit by choosing and . However, for and , with and , we can minimize the uncertainty as . This is a constant factor improvement over our scheme using separable states. This is consistent with the fact arbitrary small noise can sometimes make the entanglement sensor almost equivalent with separable sensors Huelga et al. (1997).

In conclusion, we have proposed in this letter a scheme to achieve Heisenberg limited quantum sensing of the amplitudes of globally applied fields. We have found that frequent implementations of quantum teleportation provide a suitable circumstance for sensing where the dephasing is suppressed while the information from the target fields is continuously accumulated. If perfect quantum teleportation is available, we can achieve the Heisenberg limit. Moreover, even when quantum teleportation is moderately noisy, our protocol still shows quantum enhancement over the standard Ramsey scheme.

This work was supported by JSPS KAKENHI Grant No. 15K17732 and partly supported by MEXT KAKENHI Grant No. 15H05870. S.C.B. acknowledges support from the EPSRC NQIT Hub, Project No. EP/M013243/1.


  • Budker and Romalis (2007) D. Budker and M. Romalis, Nature Physics 3, 227 (2007).
  • Balasubramanian and et al (2008) G. Balasubramanian and et al, Nature 455, 648 (2008).
  • Maze and et al (2008) J. Maze and et al, Nature 455, 644 (2008), ISSN 0028-0836.
  • Degen et al. (2016) C. Degen, F. Reinhard, and P. Cappellaro, arXiv preprint arXiv:1611.02427 (2016).
  • Takano et al. (2016) T. Takano, M. Takamoto, I. Ushijima, N. Ohmae, T. Akatsuka, A. Yamaguchi, Y. Kuroishi, H. Munekane, B. Miyahara, and H. Katori, Nature Photonics (2016).
  • Simon (1999) J. Simon, Advances in Physics 48, 449 (1999).
  • Chang et al. (1992) A. Chang, H. Hallen, L. Harriott, H. Hess, H. Kao, J. Kwo, R. Miller, R. Wolfe, J. Van der Ziel, and T. Chang, Applied physics letters 61, 1974 (1992).
  • Poggio and Degen (2010) M. Poggio and C. Degen, Nanotechnology 21, 342001 (2010).
  • Huelga et al. (1997) S. Huelga, C. Macchiavello, T. Pellizzari, A. Ekert, M. Plenio, and J. Cirac, Phys. Rev. Lett. 79, 3865 (1997).
  • Said et al. (2011) R. Said, D. Berry, and J. Twamley, Physical Review B 83, 125410 (2011).
  • Higgins et al. (2007) B. L. Higgins, D. W. Berry, S. D. Bartlett, H. M. Wiseman, and G. J. Pryde, Nature 450, 393 (2007).
  • Waldherr et al. (2012) G. Waldherr, J. Beck, P. Neumann, R. Said, M. Nitsche, M. Markham, D. Twitchen, J. Twamley, F. Jelezko, and J. Wrachtrup, Nature nanotechnology 7, 105 (2012).
  • Demkowicz-Dobrzanski et al. (2012) R. Demkowicz-Dobrzanski, J. Kolodynski, and M. Guta, Nature Communications 3 (2012).
  • Gottesman (2009) D. Gottesman, Quantum Information Science and Its Contributions to Mathematics 68, 13 (2009).
  • Viola et al. (1999) L. Viola, E. Knill, and S. Lloyd, Physical Review Letters 82, 2417 (1999).
  • Taylor et al. (2008) J. Taylor, P. Cappellaro, L. Childress, L. Jiang, D. Budker, P. Hemmer, A. Yacoby, R. Walsworth, and M. Lukin, Nature Physics 4, 810 (2008).
  • De Lange et al. (2011) G. De Lange, D. Ristè, V. Dobrovitski, and R. Hanson, Physical review letters 106, 080802 (2011).
  • Dür et al. (2014) W. Dür, M. Skotiniotis, F. Froewis, and B. Kraus, Phys. Rev. Lett. 112, 080801 (2014).
  • Herrera-Martí et al. (2015) D. A. Herrera-Martí, T. Gefen, D. Aharonov, N. Katz, and A. Retzker, Phys. Rev. Lett. 115, 200501 (2015).
  • Arrad et al. (2014) G. Arrad, Y. Vinkler, D. Aharonov, and A. Retzker, Phys. Rev. Lett. 112, 150801 (2014).
  • Kessler et al. (2014) E. M. Kessler, I. Lovchinsky, A. O. Sushkov, and M. D. Lukin, Phys. Rev. Lett. 112, 150802 (2014).
  • Cohen et al. (2016) L. Cohen, Y. Pilnyak, D. Istrati, A. Retzker, and H. Eisenberg, Phys. Rev. A 94, 012324 (2016).
  • Unden et al. (2016) T. Unden, P. Balasubramanian, D. Louzon, Y. Vinkler, M. B. Plenio, M. Markham, D. Twitchen, A. Stacey, I. Lovchinsky, A. O. Sushkov, et al., Phys. Rev. Lett. 116, 230502 (2016).
  • Sisi and et al (2017) Z. Sisi and et al, arXiv preprint arXiv:1706.02445 (2017).
  • Schmitt and et al (2017) S. Schmitt and et al, Science 356, 832 (2017).
  • Boss et al. (2017) J. Boss, K. Cujia, J. Zopes, and C. Degen, Science 356, 837 (2017).
  • Misra and Sudarshan (1977) B. Misra and E. G. Sudarshan, Journal of Mathematical Physics 18, 756 (1977).
  • Itano et al. (1990) W. Itano, D. Heinzen, J. Bollinger, and D. Wineland, Phys. Rev. A 41, 2295 (1990).
  • Facchi et al. (2001) P. Facchi, H. Nakazato, and S. Pascazio, Phys. Rev. Lett. 86, 2699 (2001).
  • Nakazato et al. (1996) H. Nakazato, M. Namiki, and S. Pascazio, Int. J. Mod. B 10, 247 (1996).
  • Raussendorf and Briegel (2001) R. Raussendorf and H. Briegel, Phys. Rev. Lett. 86, 5188 (2001).
  • Barjaktarevic et al. (2005) J. Barjaktarevic, R. McKenzie, J. Links, and G. Milburn, Physical review letters 95, 230501 (2005).
  • Silva et al. (2007) M. Silva, V. Danos, E. Kashefi, and H. Ollivier, New Journal of Physics 9, 192 (2007).
  • Olmschenk et al. (2009) S. Olmschenk, D. Matsukevich, P. Maunz, D. Hayes, L.-M. Duan, and C. Monroe, Science 323, 486 (2009).
  • Baur et al. (2012) M. Baur, A. Fedorov, L. Steffen, S. Filipp, M. Da Silva, and A. Wallraff, Phys. Rev. Lett. 108, 040502 (2012).
  • Averin et al. (2016) D. Averin, K. Xu, Y. Zhong, C. Song, H. Wang, and S. Han, Phys. Rev. Lett. 116, 010501 (2016).
  • Hornberger (2009) K. Hornberger, in Entanglement and Decoherence (Springer, 2009), pp. 221–276.
  • De Lange et al. (2010) G. De Lange, Z. Wang, D. Riste, V. Dobrovitski, and R. Hanson, Science 330, 60 (2010).
  • Yoshihara et al. (2006) F. Yoshihara, K. Harrabi, A. Niskanen, and Y. Nakamura, Phys. Rev. Lett. 97, 167001 (2006).
  • Kakuyanagi et al. (2007) K. Kakuyanagi, T. Meno, S. Saito, H. Nakano, K. Semba, H. Takayanagi, F. Deppe, and A. Shnirman, Phys. Rev. Lett. 98, 047004 (2007).
  • Kondo et al. (2016) Y. Kondo, Y. Matsuzaki, K. Matsushima, and J. G. Filgueiras, New Journal of Physics 18, 013033 (2016).
  • Palma et al. (1996) G. M. Palma, K. A. Suominen, and A. K. Ekert, Proc. R. Soc. London. Ser.A 452, 567 (1996).
  • Browne and Briegel (2006) D. E. Browne and H. J. Briegel (2006), quant-ph/0603226.
  • Gardiner and Zoller (2004) C. W. Gardiner and P. Zoller, Quantum Noise (Springer, Berlin, 2004).
  • Matsuzaki et al. (2011) Y. Matsuzaki, S. Benjamin, and J. Fitzsimons, Phys. Rev. A 84, 012103 (2011).
  • Chin et al. (2012) A. W. Chin, S. F. Huelga, and M. B. Plenio, Phys. Rev. Lett. 109, 233601 (2012).
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