Citation:

Improving the precision of optical metrology by detecting fewer photons with biased weak measurement


  • Light: Science & Applications  10, Article number: (2021)
More Information
  • Corresponding author:
    Geng Chen (chengeng@ustc.edu.cn)ChuanFeng Li (cfli@ustc.edu.cn)
  • Received: 26 November 2020
    Revised: 15 April 2021
    Accepted: 21 April 2021
    Published online: 17 May 2021

doi: https://doi.org/10.1038/s41377-021-00543-4

  • In optical metrological protocols to measure physical quantities, it is, in principle, always beneficial to increase photon number n to improve measurement precision. However, practical constraints prevent the arbitrary increase of n due to the imperfections of a practical detector, especially when the detector response is dominated by the saturation effect. In this work, we show that a modified weak measurement protocol, namely, biased weak measurement significantly improves the precision of optical metrology in the presence of saturation effect. This method detects an ultra-small fraction of photons while maintains a considerable amount of metrological information. The biased pre-coupling leads to an additional reduction of photons in the post-selection and generates an extinction point in the spectrum distribution, which is extremely sensitive to the estimated parameter and difficult to be saturated. Therefore, the Fisher information can be persistently enhanced by increasing the photon number. In our magnetic-sensing experiment, biased weak measurement achieves precision approximately one order of magnitude better than those of previously used methods. The proposed method can be applied in various optical measurement schemes to remarkably mitigate the detector saturation effect with low-cost apparatuses.
  • 加载中
  • [1] Helstrom, C. W. Quantum Detection and Estimation Theory (Academic, 1976).
    [2] Holevo, A. S. Probabilistic and Statistical Aspects of Quantum Theory (North-Holland, 1982).
    [3] Yuan, H. D. & Fung, C. H. F. Quantum parameter estimation with general dynamics. npj Quantum Inf. 3, 14 (2017). doi: 10.1038/s41534-017-0014-6
    [4] Wineland, D. J. et al. Spin squeezing and reduced quantum noise in spectroscopy. Phys. Rev. A 46, R6797 (1992). doi: 10.1103/PhysRevA.46.R6797
    [5] Caves, C. M. Quantum-mechanical noise in an interferometer. Phys. Rev. D 23, 1693 (1981). doi: 10.1103/PhysRevD.23.1693
    [6] Lee, H., Kok, P. & Dowling, J. P. A quantum Rosetta stone for interferometry. J. Mod. Opt. 49, 2325–2338 (2002). doi: 10.1080/0950034021000011536
    [7] Braunstein, S. L. Quantum limits on precision measurements of phase. Phys. Rev. Lett. 69, 3598–3601 (1992). doi: 10.1103/PhysRevLett.69.3598
    [8] Giovannetti, V., Lloyd, S. & Maccone, L. Quantum-enhanced measurements: beating the standard quantum limit. Science 306, 1330–1336 (2004). doi: 10.1126/science.1104149
    [9] Tan, S. H. et al. Quantum illumination with Gaussian states. Phys. Rev. Lett. 101, 253601 (2008). doi: 10.1103/PhysRevLett.101.253601
    [10] Pirandola, S. Quantum reading of a classical digital memory. Phys. Rev. Lett. 106, 090504 (2011). doi: 10.1103/PhysRevLett.106.090504
    [11] Nair, R. & Yen, B. J. Optimal quantum states for image sensing in loss. Phys. Rev. Lett. 107, 193602 (2011). doi: 10.1103/PhysRevLett.107.193602
    [12] Bollinger, J. J. et al. Optimal frequency measurements with maximally correlated states. Phys. Rev. A 54, R4649 (1996). doi: 10.1103/PhysRevA.54.R4649
    [13] Walther, P. et al. De Broglie wavelength of a non-local four-photon state. Nature 429, 158–161 (2004). doi: 10.1038/nature02552
    [14] Afek, I., Ambar, O. & Silberberg, Y. High-NOON states by mixing quantum and classical light. Science 328, 879–881 (2010). doi: 10.1126/science.1188172
    [15] Goda, K. et al. A quantum-enhanced prototype gravitational-wave detector. Nat. Phys. 4, 472–476 (2008). doi: 10.1038/nphys920
    [16] Grangier, P. et al. Squeezed-light-enhanced polarization interferometer. Phys. Rev. Lett. 59, 2153–2156 (1987). doi: 10.1103/PhysRevLett.59.2153
    [17] Xiao, M., Wu, L. A. & Kimble, H. J. Precision measurement beyond the shot-noise limit. Phys. Rev. Lett. 59, 278–281 (1987). doi: 10.1103/PhysRevLett.59.278
    [18] Treps, N. et al. A quantum laser pointer. Science 301, 940–943 (2003). doi: 10.1126/science.1086489
    [19] Giovannetti, V., Lloyd, S. & Maccone, L. Advances in quantum metrology. Nat. Photonics 5, 222–229 (2011). doi: 10.1038/nphoton.2011.35
    [20] Schnabel, R. Squeezed states of light and their applications in laser interferometers. Phys. Rep. 684, 1–51 (2017). doi: 10.1016/j.physrep.2017.04.001
    [21] Aharonov, Y., Albert, D. Z. & Vaidman, L. How the result of a measurement of a component of the spin of a Spin-1/2 particle can turn out to be 100. Phys. Rev. Lett. 60, 1351–1354 (1988). doi: 10.1103/PhysRevLett.60.1351
    [22] Aharonov, Y. & Vaidman, L. in Time in Quantum Mechanics (eds Muga, J. G., Mayato, R. S. & Egusquiza, I. L. ) 369–412 (Springer, 2002).
    [23] Hallaji, M. et al. Weak-value amplification of the nonlinear effect of a single photon. Nat. Phys. 13, 540–544 (2017). doi: 10.1038/nphys4040
    [24] Hosten, O. & Kwiat, P. Observation of the spin hall effect of light via weak measurements. Science 319, 787–790 (2008). doi: 10.1126/science.1152697
    [25] Dixon, P. B. et al. Ultrasensitive beam deflection measurement via interferometric weak value amplification. Phys. Rev. Lett. 102, 173601 (2009). doi: 10.1103/PhysRevLett.102.173601
    [26] Brunner, N. & Simon, C. Measuring small longitudinal phase shifts: weak measurements or standard interferometry? Phys. Rev. Lett. 105, 010405 (2010). doi: 10.1103/PhysRevLett.105.010405
    [27] Xu, X. Y. et al. Phase estimation with weak measurement using a white light source. Phys. Rev. Lett. 111, 033604 (2013). doi: 10.1103/PhysRevLett.111.033604
    [28] Ferrie, C. & Combes, J. Weak value amplification is suboptimal for estimation and detection. Phys. Rev. Lett. 112, 040406 (2014). doi: 10.1103/PhysRevLett.112.040406
    [29] Knee, G. C. & Gauger, E. M. When amplification with weak values fails to suppress technical noise. Phys. Rev. X 4, 011032 (2014). doi: 10.1103/PhysRevX.4.011032
    [30] Dressel, J. et al. Colloquium: understanding quantum weak values: basics and applications. Rev. Mod. Phys. 86, 307–316 (2014). doi: 10.1103/RevModPhys.86.307
    [31] Zhang, L. J., Datta, A. & Walmsley, I. A. Precision metrology using weak measurements. Phys. Rev. Lett. 114, 210801 (2015). doi: 10.1103/PhysRevLett.114.210801
    [32] Alves, G. B. et al. Weak-value amplification as an optimal metrological protocol. Phys. Rev. A 91, 062107 (2015). doi: 10.1103/PhysRevA.91.062107
    [33] Vaidman, L. Weak value controversy. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 375, 20160395 (2017). doi: 10.1098/rsta.2016.0395
    [34] Xu, L. et al. Approaching quantum-limited metrology with imperfect detectors by using weak-value amplification. Phys. Rev. Lett. 125, 080501 (2020). doi: 10.1103/PhysRevLett.125.080501
    [35] Harris, J., Boyd, R. W. & Lundeen, J. S. Weak value amplification can outperform conventional measurement in the presence of detector saturation. Phys. Rev. Lett. 118, 070802 (2017). doi: 10.1103/PhysRevLett.118.070802
    [36] Zhang, Z. H. et al. Ultrasensitive biased weak measurement for longitudinal phase estimation. Phys. Rev. A 94, 053843 (2016). doi: 10.1103/PhysRevA.94.053843
    [37] Li, D. M. et al. Optical rotation based chirality detection of enantiomers via weak measurement in frequency domain. Appl. Phys. Lett. 112, 213701 (2018). doi: 10.1063/1.5019816
    [38] Li, D. M. et al. A chiral sensor based on weak measurement for the determination of Proline enantiomers in diverse measuring circumstances. Biosens. Bioelectron. 110, 103–109 (2018). doi: 10.1016/j.bios.2018.03.033
    [39] Jaynes, E. T. in Probability Theory: The Logic of Science (ed. Bretthorst, G. L. ) (Cambridge University Press, 2003).
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(5)

Article Metrics

Article views(129) PDF downloads(54) Citation(0) Citation counts are provided from Web of Science. The counts may vary by service, and are reliant on the availability of their data.

Improving the precision of optical metrology by detecting fewer photons with biased weak measurement

  • 1. CAS Key Laboratory of Quantum Information, University of Science and Technology of China, 230026 Hefei, China
  • 2. CAS Center For Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, 230026 Hefei, Anhui, China
  • 3. National Laboratory of Solid State Microstructures and College of Engineering and Applied Sciences, Nanjing University, Nanjing, China
  • 4. Collaborative Innovation Center of Advanced Microstructures, Nanjing University, 210093 Nanjing, China
  • Corresponding author:

    Geng Chen, chengeng@ustc.edu.cn

    Chuan-Feng Li, cfli@ustc.edu.cn

doi: https://doi.org/10.1038/s41377-021-00543-4

Abstract: In optical metrological protocols to measure physical quantities, it is, in principle, always beneficial to increase photon number n to improve measurement precision. However, practical constraints prevent the arbitrary increase of n due to the imperfections of a practical detector, especially when the detector response is dominated by the saturation effect. In this work, we show that a modified weak measurement protocol, namely, biased weak measurement significantly improves the precision of optical metrology in the presence of saturation effect. This method detects an ultra-small fraction of photons while maintains a considerable amount of metrological information. The biased pre-coupling leads to an additional reduction of photons in the post-selection and generates an extinction point in the spectrum distribution, which is extremely sensitive to the estimated parameter and difficult to be saturated. Therefore, the Fisher information can be persistently enhanced by increasing the photon number. In our magnetic-sensing experiment, biased weak measurement achieves precision approximately one order of magnitude better than those of previously used methods. The proposed method can be applied in various optical measurement schemes to remarkably mitigate the detector saturation effect with low-cost apparatuses.

Introduction
  • Scientific communities pursue higher precision in the measurement of various quantities. Quantum metrology can potentially surpass classical protocols by exploiting quantum re- sources1-11, e.g., NOON states and squeezed states12-18. However, these quantum resources are intricate to prepare and control with currently available techniques19, 20. Another solution is to directly increase the copies of meter state, e.g., the number of photons n in the measurement of optical phase with an interferometer. In this case, the signal-to-noise ratio is proportional to $\sqrt n$. A main constraint of this method is the detector saturation effect (DSE) which occurs in various measurement scenarios and eventually damages the precision. Therefore, how to alleviate DSE and further enhance the precision is an interesting problem worth investigating. From a practical point of view, this requires a small fraction of photons being detected while maintaining almost the same metrological information as that contained in all incident photons. This requirement seems to be paradoxical because discarding photons inevitably leads to loss of information in general.

    Standard weak measurement (SWM) is an innovative method to determine small physical quantities that are impractical to measure using conventional measurement (CM)21-25. Especially, for longitudinal optical phase measurement, the method of measuring the spectrum shift of a broad-band light beam in SWM has been verified to be preferable than CM both in theory26 and experiment27. Although some theoretical papers argue that SWM is suboptimal28-30, it is widely acknowledged that the optimal precision of SWM is comparable with that of CM, even though only a small fraction of photons are post- selected for detection31, 32. In other words, SWM simultaneously amplifies profile shift and reduces the average photon number received by the detector. Based on this point, Vaidman conjectured that SWM can effectively alleviate DSE and outperform CM with incident photon number above the saturation threshold of detectors33. Recently, several theoretical and experimental studies have confirmed this advantage of SWM and proved that SWM offers an improved precision compared to CM in the presence of DSE34.

    In principle, two premises endow us with the ability to mitigate DSE and attain better precision; namely, fewer photons are detected and these photons contain more metrological information than the discarded ones. SWM satisfies these two premises by increasing the strength of post-selection, and eventually, an improved precision can be acquired35. In an ideal case, a stronger post-selection in SWM necessarily results in a larger factor of amplification; however, in practice, this simple post-selection cannot be arbitrarily strong, and SWM can only alleviate DSE to a limited extent34. By contrast, biased weak measurement (BWM) employs an additional reduction of photons in the post-selection by introducing a pre-coupling, and the remaining photons have shown to be extremely sensitive to the estimated parameter both in theory and experiments36-38. In this work, we demonstrate that BWM can effectively mitigate DSE and obtain much higher precision than both CM and SWM. Briefly speaking, BWM is impervious to DSE because an extinction point appears in the spectrum distribution; therefore, the number of detected photons is greatly reduced. What is more, these photons are much more sensitive to the estimated parameter than those detected in CM and SWM. The advantage of BWM is rigorously cast in terms of FI, and the results demonstrate that the Fisher information (FI) of BWM can grow persistently with increasing n. By contrast, the accessible FI of CM and SWM is much less since DSE dominates the detector response for a much lower incident photon number n. The advantages of BWM are experimentally demonstrated through the sensing of a static magnetic field, where the highest precision of BWM exceeds that of SWM by nearly one order of magnitude, and this superiority is further contrasted compared to CM. We believe that the proposed method can shed light on various measurement scenarios suffering from DSE.

Results
  • In the following discussion, we take $\hbar = 1$. Without loss of generality, we consider a scheme to measure a small optical delay ${{\rm{\tau }}} = \frac{\phi }{{cp_0}}$, which introduces an additional optical phase ϕ between two orthogonal polarization components $|0 > $ and $|1 > $ for the photons with momentum and speed. The photon momentum is $p_0 = \frac{{2\pi }}{{\lambda _0}} = \frac{{\omega _0}}{c}$, where λ00) denotes the central wavelength (frequency) of incident light. Theoretically, the coupling strength k = cτ can be estimated by the interaction between the system (which is initialized to $|\varphi _i > = \frac{1}{{\sqrt 2 }}(|0 > + |1 > )$) and the meter (which is initialized to ${\int} {{{{\mathrm{d}}p}}|{{\rm{\psi }}}} \left({{p}} \right) > $ and is assumed to have a Gaussian profile with mean value p0 and variance $({\Delta}{\mathrm{p}})^2$) where the Hamiltonian is ${{H}} = {{k}}\delta ({{t}} - t_0)\hat A\hat P$, in which $\hat A = \left| {0 > < 0} \right| - |1 > < 1|$ is the system operator, and $\hat P$ is the momentum operator of the photon.

    The CM, SWM, and BWM schemes to measure τ are diagrammed in Fig. 1. For CM, the system-meter coupling can be described by the unitary transformation ${{U}} = e^{ - ik{\hat {A}}{\hat {P}}}$ and the final joint state is given as follows:

    Fig. 1  Diagram of CM, SWM, and BWM schemes.

    The procedures to implement CM, SWM, and BWM are shown in a-c, respectively. All these three methods start from a system initialization, and all involve a coupling between the system and meter. In CM, a projective measurement is made on the system and the coupling strength can be estimated from the change in photon counting for each pixel. In SWM, a post-selection is applied to amplify the shift of the mean value of observable $\hat P$. This post-selection makes SWM more robust to DSE compared with CM; however, the detector array eventually saturates for these two methods when too many photons are received. In BWM, an additional pre-coupling is introduced before the coupling; consequently, the post-selection discards more photons than SWM and leads to an extinction point for the distribution of p, which enables the detector array to work below the threshold for a much larger photon number. Furthermore, the position of this extinction point is extremely sensitive to the coupling strength. Therefore, BWM is more robust to DSE and eventually attains better precision than those of CM and SWM
    $$ |\psi\, > \,_{\mathrm{joint}} = {\int} {{\mathrm{d}}p} \left( {e^{ipk}\left| {0\, > + e^{ - ipk}} \right|1 > } \right)|\psi \left( p \right) > $$ (1)

    where p is the eigenvalue of ${\hat{ P}}$. Then the system is projected on a certain basis that leads to an appreciable selection probability, e.g., $\frac{1}{{\sqrt 2 }}(\left| {0 > - i} \right|1 > )$, which leads to an unnormalized redistribution of p to be

    $$ D(p)_{\mathrm{CM}} = {\mathrm{sin}}^2\left( {\frac{\pi }{4} + pk} \right)\left| {\, < \,\psi \left( p \right)} \right|\psi \left( p \right)\, > \,|^2 $$ (2)

    and the shift of the mean value of p when ${{k}}p_0 \ll 1$ is calculated as

    $$ \delta p_{\mathrm{CM}} = \frac{{2k({\mathrm{{\Delta}}}p)^2{\mathrm{cos}}(kp_0)}}{{\sin \left( {2kp_0} \right) + e^{2k^2({\mathrm{{\Delta}}}p)^2}}} \simeq 2k\left( {{\mathrm{{\Delta}}}p} \right)^2 $$ (3)

    For SWM, a normal post-selection into $|{{\rm{\varphi }}}\, > \, _f = \frac{1}{{\sqrt 2 }}(e^{ - i{\, \it\, {\epsilon }}}\left| {0 > - e^{i{\it{\epsilon }}}} \right|1 > )$ is made on the system, and the distribution of p in this post-selected meter state is given as:

    $$ D(p)_{\mathrm{SWM}} = {\mathrm{sin}}^2\left( {pk + {\it{\epsilon }}} \right)|\, < \,\psi \left( p \right)|\psi \left( p \right)\, > \,|^2 $$ (4)

    when τ, ϵ ≪ 1, the value of k can be estimated through the shift of the mean value of p, which can be calculated as follows:

    $$ \delta p_{\mathrm{SWM}} = 2k({\mathrm{{\Delta}}}p)^2{\mathrm{cot}}{\it{\epsilon }} \simeq \frac{{2k\left( {{\mathrm{{\Delta}}}p} \right)^2}}{{\it{\epsilon }}} $$ (5)

    which is amplified by a factor of ${\mathrm{cot}}\epsilon$ compared to the shift in CM. The price for this amplification is post-selecting the photons with probability ${{O}}({\it{\epsilon }}^2)$.

    As shown in Fig. 1c for the BWM procedures, here, the main difference from SWM is an additional step to bias the meter before the coupling that encodes the parameter. Specifically, a predetermined delay $\frac{{\beta }}{{{c}}}$ is introduced between the two components of the system observable with β satisfying ${{p}}_0\beta + {\it{\epsilon }} = m\pi$ (m is an integer), and the corresponding distribution of p of the post-selected meter state is given as:

    $$ D(p)_{\mathrm{BWM}} = {\mathrm{sin}}^2(p\left( {k + \beta } \right) + {\it{\epsilon }})| < \psi \left( p \right)|\psi \left( p \right) > |^2 $$ (6)

    It is evident that when k = 0, an extinction point appears for p = p0, as shown in Fig. 1b. It has been suggested that the position of this extinction point is extremely sensitive to k31; i.e., even very small k yields a perceptible shift of this extinction point. The mean value shift of p for m = 0 in BWM is calculated as follows:

    $$ \delta p_{\mathrm{SWM}} \simeq \frac{{2k(p_0)^2}}{{\it{\epsilon }}} $$ (7)

    Since p0 is usually larger than its uncertainty Δp by at least one order of magnitude for a visible laser beam, and according to Eqs. (5) and (7), the mean value shift in BWM scheme is much larger than that in SWM. Correspondingly, this pre-coupling leads to an additional reduction of photons in the post-selection, which cannot be achieved by decreasing ϵ, and the post-selection probability is ${\mathrm{O}}(({\mathrm{{\Delta}}}p{\it{\epsilon }}/p_0)^2)$.

    The CM method is equivalent to an interferometer, in which two outcomes are obtained through a projective measurement, and approximately half of the photons are detected for each outcome26. Normally, τ can be estimated by simply summing up the photon number change over all the components of p of each outcome. As opposed to CM, SWM postselects a small fraction of photons, and the spectrum shift of SWM is amplified by the weak value ${{A}}_\omega = \frac{{ < \varphi _f\left| A \right|\varphi _i > }}{{ < \varphi _f|\varphi _i > }} = i{\mathrm{cot}}{\it{\epsilon }}$.

    Compared to SWM, BWM applies a pre-coupling procedure, which introduces an extinction effect and gives a further amplified mean value shift of p. From the above discussions, it can be concluded that among these three schemes, BWM acquires the largest meter shift and detects the fewest photons. As a result, as shown in Fig. 1, when a sufficiently large number of incident photons induces a flattening distribution on the detector array for CM and SWM schemes, the response of the detectors in BWM is maintained in the dynamic range of the detectors. Therefore, one can expect that BWM is more robust against DSE and will eventually outperform CM and SWM, and we give firm evidence for this advantage with both numerical calculation and experimental demonstration.

  • In this subsection, the advantage of BWM to mitigate DSE is verified by calculating the FI in a specific measurement scenario, which provides a lower bound for the uncertainty of the estimation of a parameter39. The value of FI is calculated by summing up the FI obtained for each component of p. Consider the specific experiment scheme shown in Fig. 2, which is proposed to sense the magnetic induction strength B. The Hamiltonian ${{H}} = {{k}}\delta \left({{{t}} - t_0} \right)\hat A\hat P$ couples the system and the meter with strength ${{k}} = {{VBl}}/p_0$, where V and l are the Verdet constant and length of the Faraday crystal, respectively, ${\hat{\mathrm A}} = |R > < R| - |L > < L|$ is the system operator, $|{{R}}\, > $ (|$|{{L}} > $) is the right(left) circularly polarized component of light and $\hat P$ is the momentum operator of photons. Therefore, B can be determined by studying the distribution of p (as suggested in Eqs. (4) and (6)).

    Fig. 2  Experimental setup to sense static magnetic field with CM, SWM, and BWM schemes.

    A pulsed laser centering at λ0 = 796 nm is horizontally polarized ($|{{H}} > $) after passing through a polarized beam splitter and then rotated to be 45° diagonal polarized ($|{{H}} + {{V}} > $) by a half-wave plate (HWP). The phase compensation plate introduces a biased phase between the $|{{H}} > $ and $|{{V}} > $ components, which are transformed to $|{{R}} > $ and $|{{L}} > $ after passing through HWP and the quarter-wave plate (QWP). A Faraday crystal (FC) placed between two electric coils introduces a relative phase between $|{{R}} > $ and $|{{L}} > $ that is proportional to the magnetic induction strength B. Both projective measurement in CM and post-selecting in SWM and BWM are implemented using a polarizer. By dispersing the beam with a grating, both the change of photon numbers and spectrum redistribution are recorded by a CMOS, and B can be determined from the corresponding spectrum shift

    Note that ${\mathrm{{\Delta}}}p \ll p_0$ and the light propagates along a single direction; thus, we can measure the distribution of p using a spectrometer based on the relation ${{p}} = \frac{{2{\pi}}}{{\lambda }}$. To be specific, the light is dispersed on the grating, and the photons with momentum pj are received by the jth pixel of the detector array, which is a complementary metal-oxide semiconductor (CMOS) in our experiment. The distribution of p is recorded as a frame by reading the number of excited electrons on each pixel.

    The value of FI can be calculated by summing up the FI of all the pixels on the CMOS, and the FI of the jth pixel can be obtained from the probability of exciting kj electrons. When the total number of incident photon is n, kj can be calculated as follows:

    $$ {{P}}\left( {k_j{\mathrm{|}}B} \right) = \mathop {\sum }\limits_{N_j} R_s\left( {k_j{\mathrm{|}}N_j} \right)P(N_j|\bar n_j\left( B \right),\sigma _j(B)) $$ (8)

    here, $P(N_j|\bar n_j\left(B \right), \sigma _j\left(B \right))$ is the Gaussian distribution with average photon number $\bar n_j\left(B \right) = n\mathop {\int }\nolimits_j^{} {\mathrm{d}}pD(p)_{{\mathrm{SWM}}({\mathrm{BWM}})}$ and standard deviation $\sigma _j(B)$ determined by analyzing the recorded frames for estimation. Note that $D(p)_{{\mathrm{SWM}}({\mathrm{BWM}})}$ is related to B through ${{k}} = {{VBl}}/p_0$; thus, the average photon number on each pixel is determined by B. $R_s\left({k_j{\mathrm{|}}N_j} \right)$ is the probability of generating kj electrons when the jth pixel receives exactly Nj photons, and the concrete expression gives a quantitative description of the response model of CMOS (see "Materials and methods" for details).

    The calculated FI against n is shown in Fig. 3a, in which we set m = 5 and B = 0.028T to be consistent with those applied in experiment. As can be seen, the FI of CM firstly reaches its maximum when n is ~5 × 106 since the DSE dominates the response for some of the pixels. At this stage, each pixel responds in the dynamic range for SWM and BWM, and the elicited FI grows with increasing n. When n increases to 108, SWM loses its advantage because DSE begins to undermine the performance of SWM, and the FI decreases gradually. When n exceeds 109, nearly all the pixels saturate in SWM, and the distribution carries negligible information about B. Consequently, the elicited FI in SWM drops to zero rapidly, as shown in Fig. 3a. As expected, BWM behaves robustly to DSE, and the FI grows consistently with increasing n. The primary limitation factors for this positive correlation between FI and n are the finite extinction ratio and pixel size in a practical experiment, which cause a small portion of photons to shine on the extinction point. Consequently, for BWM the extinction point eventually saturates for very large n, and the FI decays after reaching its maximum value, as shown in Fig. 3a. Nevertheless, the maximal FI of BWM is larger than those of CM and SWM by nearly three and two orders of magnitude, respectively. For a comprehensive comparison between BWM and SWM, further calculations are made for some small values of $\epsilon$ with m = 0 and ${{B}} = 1.3 \times 10^{ - 7}T$ as shown in Fig. 3b. Theoretically, the post-selection probability of SWM can be quadratically reduced by decreasing the value of $\epsilon$, and DSE can thus be effectively suppressed to acquire a better precision. However, BWM still exhibits a significant advantage in the achievable precision even for very small $\epsilon$.

    Fig. 3  FI against total incident photon number n to sense a static magnetic field.

    a By varying the incident photon number, the amount of elicited classical FI is calculated when the extinction ratio and post-selection strength $\epsilon$ are set to 90, 000 and 0.2, respectively. Compared to SWM and CW, BWM demonstrates excellent ability to subdue DSE and allows much more FI to be obtained. b For three values of $\epsilon$, the FI is calculated for SWM and BWM protocols with infinite extinction ratio. For each value of $\epsilon$, BWM detects much fewer photons than SWM, and the achievable FI significantly outperforms that of SWM
  • The advantage of BWM to mitigate DSE is demonstrated with the setup shown in Fig. 2. The static magnetic field produced by two electric coils is sensed through CM, SWM, and BWM schemes, and the change in the magnetic field can cause a spectral redistribution; thus, a more distinct redistribution results in greater measurement sensitivity. Figure 4a-c shows the normalized spectral distribution before and after applying the magnetic field for CM, SWM, and BWM, respectively. Note that the spectrum change in both CM and SWM is too subtle to be observed, and the spectrum change in BWM is transformed to a new pattern with perceptible distinguishability. These results indicate that BWM realizes a higher meter shift in measurement, as predicted by Eqs. (3), (5), and (7). The robustness to DSE can be revealed through the electron number distribution of CMOS for varying n, as shown in Fig. 4d-f for CM, SWM, and BWM, respectively. In CM, the pixels begin to saturate for 4.8 × 106 photons and completely saturate for 108 photons. In SWM, saturation begins when n = 107 and the profile completely flattens for n = 109 photons. Because of the ultra-sensitive extinction point in BWM, the electron number distribution of CMOS is not saturated up to 1010 photons and consistently provides a considerable amount of FI.

    Fig. 4  Spectrum distribution in CM, SWM, and BWM schemes.

    The spectrum re- distributions of a CM, b SWM, and c BWM are shown when the magnetic field is turned on. Here, blue and red curves represent the spectrum distribution before and after turning on the magnetic field. In d-f the electron number redistributions for CM, SWM, and BWM are shown for varying numbers of the incident photon. Typically, a larger photon number leads to more pixels becoming saturated in the CMOS. For CM and SWM, nearly all the pixels saturate when n approaches 108 and 109, and then the CMOS outputs a flattening profile, which limits the precision of the two protocols. By contrast, BWM is more robust to DSE because there is an ultra-sensitive extinction point, which is difficult to saturate even for 1010 photons

    By recording 6000 frames of electron number distribution for each value of n, maximum likelihood estimation (MLE) is utilized to estimate B. Briefly speaking, one estimation of B is given by MLE using v = 300 frames that are uniformly and randomly selected from 6000 frames recorded by CMOS. By repeating the MLE 100 times, we take the standard deviation of these 100 estimates as the precision ΔB (see "Materials and methods" for details). Figure 5 shows the precisions of CM, SWM, and BWM schemes with varying n. For a fixed value of n, the photon numbers participating in the system-meter coupling are identical for these three schemes. Eventually, almost half the incident photons are detected in CM for each outcome of the system, while in SWM and BWM only a small proportion of the photons are post-selected and thus detected. Although the post-selection parameter (the setting of polarizer) is the same for both BWM and SWM, much fewer photons are detected in BWM than that in SWM, since the system-meter joint states are different for these two schemes due to the pre-coupling in BWM. Consequently, for ${{n}} \ge 10^7$, BWM achieves a better precision than those of SWM and CM, and this advantage remains until the detector saturates in BWM.

    Fig. 5  Precision of magnetic sensing for CM, SWM, and BWM schemes.

    The precision ΔB obtained using CM, SWM, and BWM schemes are plotted versus the incident photon number n. CM and SWM achieve their best precisions of 2.2 × 104T and 3.63 × 105T with 5.9 × 106 and $9.7 \times 10^7$ photons, respectively. Afterward, the precisions degrade since DSE occurs. Because of the extinction effect, the precision of BWM keeps to be improved up to $5.08 \times 10^{10}$ photons and the best precision is $4.05 \times 10^{ - 6}T$. The theoretical precisions calculated from the FI for CM, SWM, and BWM schemes are shown as the solid lines in the same color with their respective experimental results

    Concerning the precision change with increasing n, when n is around $10^6$, and, as predicted, the precisions for all schemes improve with increasing n. When n approaches $5.9 \times 10^6$, the precision of CM reaches its minimum value $2.2 \times 10^{ - 4}T$ and then increases since DSE occurs. SWM scheme reaches its best precision $3.63 \times 10^{ - 5}T$ when ${{n}} = 9.7 \times 10^7$, and then degrades gradually; further increasing n will cause all pixels saturate for SWM; thus, the precision degrades rapidly, the MLE cannot converge and fails to give a reasonable estimate of B. By contrast, the precision for BWM continues to be enhanced with increasing n until $5.08 \times 10^{10}$, and the best precision $4.05 \times 10^{ - 6}T$ is obtained. Through SWM, the precision is improved by $\sim 6.1$ times compared to that of CM, and BWM further expands this superiority and achieves the best precision outperforming that of SWM by nearly one order of magnitude. The theoretical precision is calculated from the FI in Fig. 3a, and the resulting lines exhibit a trend that is similar to the experimental results with close values.

Discussion
  • Intuitively, by decreasing $\epsilon$, the post-selection probability in SWM can be infinitely minimized to produce an arbitrarily large weak value; therefore, it seems that DSE can be circumvented by SWM. Unfortunately, realistic optical elements can only achieve a limited $\epsilon$ on this normal post-selection and the surviving photons maintain a single peak structure, as shown in Fig. 1b. By introducing the pre-coupling, the number of photons surviving the post-selection is further reduced, and the spectrum is specially modified to generate an ultra-sensitive extinction point, which provides considerable FI even for very large incident photon numbers. Therefore, the advantage of BWM remains for very small values of $\epsilon$ (i.e., strong post-selection), which is clearly indicated by the theoretical calculation shown in Fig. 3b. Although stronger post-selection leads to more FI (better precision) in SWM, the FI of BWM protocol significantly exceeds that of SWM for each value of $\epsilon$.

    This advantage is also well confirmed by our experiment, in which BWM achieves precision that surpasses that of SWM by nearly one order of magnitude. Although our experiment demonstrates a magnetic-sensing scenario by measuring the spectrum, the proposed method is applicable to various optical measurement tasks because the extinction effect can occur in both frequency and time domains. In summary, our results pave the way to circumvent the limitation of DSE and realize higher precision in a low-cost manner, which explores the regime where current proposals fail in the presence of DSE.

Materials and methods
  • We primarily considered three effects in the response model ${{R}}_s(k_j|N_j)$. The first is the dark noise of CMOS. In our experiment, by taking 6000 frames without any incident light, we find that the distribution of dark noise satisfies the normal distribution ${{P}}_{\mathrm{d}}(k_{\mathrm{d}})\sim N(k_{{\mathrm{d}}0}, \sigma _{\mathrm{d}})$, where $k_{{\mathrm{d}}0} = 94.16$ and $\sigma _{\mathrm{d}} = 2.03$. The second effect is that the distribution of electrons excited by photons also satisfies a normal distribution ${{P}}_Q\left({k_j{\mathrm{|}}N_j} \right)\sim N({{\rm{\eta }}}N_j, \sigma _{N_j})$, where ${{\rm{\eta }}} = 31.3{\mathrm{\% }}$ is the quantum efficiency of the CMOS and ln $\sigma _{N_j} = 0.5908\, {\mathrm{ln}}\, N_j - 1.9986$. Thus, the electron distribution ${{R}}(k_j|N_j)$ is given by the convolution of the dark noise distribution and the electron distribution excited by photons:

    $$ {{R}}\left( {k_j{\mathrm{|}}N_j} \right) = \mathop {\sum }\limits_{k_{\mathrm{d}} = 58}^{k_{\mathrm{d}} = 140} {{P}}_Q\left( {k_j - d_{\mathrm{d}}{\mathrm{|}}N_j} \right){{P}}_{\mathrm{d}}(d_{\mathrm{d}}) $$ (9)

    here, we only take the sum over kd from 58 to 140 because the marginal probability of dark noise beyond this range is negligible. The third effect we must consider is the saturation effect. In our experiment, we set the saturation threshold to ${{k}}_s = 1200$. The overall response model is given as follows:

    $$ R_s\left( {k_j{\mathrm{|}}N_j} \right) = f\left( x \right) = \left\{ {\begin{array}{*{20}{c}} {R\left( {k_j{\mathrm{|}}N_j} \right),\;k_j < 1200} \\ {1 - \mathop {\sum }\nolimits_{k_{j\, < \,1200}} {{R}}\left( {k_j{\mathrm{|}}N_j} \right),\;k_j = 1200} \\ {0,\;k_j\, > \,1200} \end{array}} \right. $$ (10)

    In the saturation model, the saturation threshold is set to an artificial value of 1200 because the CMOS response becomes chaotic when the registered electron is above 1200 and cannot be described by a valid response model, which is required for the FI calculation and the use of MLE.

    With the probabilities of the readout electron numbers on each pixel, FI is calculated as follows:

    $$ {\mathrm{FI}} = \mathop {\sum }\limits_j \mathop {\sum }\limits_{k_j} \frac{{(\frac{{\partial P(k_j|B)}}{B})^2}}{{P(k_j|B)}} $$ (11)

    Several experimental parameters must be determined to calculate FI. Here, the wavelength of the laser is centered at ${{\rm{\lambda }}}_0 = 796\, {\mathrm{nm}}$ with ~12 nm full width at half maximum, and p0 is set to be $\frac{{2{{\rm{\pi }}}}}{{796}}\, {\mathrm{nm}}^{ - 1}$. Note that the Verdet constant V is approximately a constant in this 12 nm bandwidth, and it is measured as $70.35\; {\mathrm{rad}} \cdot T^{ - 1} \cdot m^{ - 1}$ for the utilized 1-cm-long FC. In addition, the dispersion relation must be calibrated because it determines ${\bar{n}}_j$. In this experiment, the dispersion relation is measured by determining the central wavelength of photons received by the jth pixel. To find this relation, we insert an etalon right before the grating and the wavelengths of the transmission peaks are measured by a fiber spectrometer. Knowing the wavelength of each peak which imposes on the jth pixel of the CMOS, we can find the relation between the wavelength λ and pixel number j is ${{\rm{\lambda }}} = 0.007331{\mathrm{j}} + 789.5$. The spectral profile $\left| { < {{\rm{\psi }}}\left({{p}} \right)} \right|{{\rm{\psi }}}\left({{p}} \right) > |^2$ is measured by the CMOS working in the dynamic range.

  • We employ a bootstrap method to obtain ΔB. We randomly select 300 frames from the 6000 taken frames. These 300 frames are used to obtain an estimate of B using the MLE method. To implement MLE, we must first define the loss function as follows:

    $$ L\left( {{B}} \right) = \mathop {\prod }\limits_{i = 1}^{300} \mathop {\prod }\limits_{j = 1}^{1920} R_s\left( {k_{ij}{\mathrm{|}}N_j} \right)P(N_j|\bar n_{ij}\left( B \right),\sigma _{ij}({{B}})) $$ (12)

    where j is the pixel number in a row on the CMOS and i is the frame number, and $\sigma _{ij}({{B}})$ is estimated as the standard deviation of electron counts of the 300 recorded frames. Here, the value of B is identified by maximizing this loss function. We repeat this process 100 times, and finally, we obtain 100 estimates of B and take the standard deviation as ΔB.

Acknowledgements
  • This work was supported by the National Key Research and Development Program of China (Nos. 2017YFA0304100, 2016YFA0302700), National Natural Science Foundation of China (Grant Nos. 11874344, 92065107, 61835004, 11774335, 91536219, 11821404), Key Research Program of Frontier Sciences, CAS (No. QYZDY-SSW-SLH003), Anhui Initiative in Quantum Information Technologies (AHY020100, AHY060300), the Fundamental Research Funds for the Central Universities (Grant Nos. WK2030020019, WK2470000026), Science Foundation of the CAS (No. ZDRW-XH-2019-1).

Conflict of interest
  • The authors declare no competing interests.

Reference (39)

Catalog

    /

    DownLoad:  Full-Size Img PowerPoint
    Return
    Return