Hong-Ou-Mandel interferometry on a biphoton beat note

TL;DR

This paper demonstrates that high-precision Hong-Ou-Mandel interferometry can be achieved using biphoton states with only two well-separated frequencies, simplifying quantum sensing requirements.

Contribution

It challenges the assumption that large bandwidth is necessary for precise measurements, showing that discrete color entanglement suffices.

Findings

Achieved high-precision delay measurements with biphotons of two frequencies.

Validated the approach experimentally by detecting delays in optical fibers.

Used Fisher Information to optimize measurement points.

Abstract

Hong-Ou-Mandel interference, the fact that identical photons that arrive simultaneously on different input ports of a beam splitter bunch into a common output port, can be used to measure optical delays between different paths. It is generally assumed that great precision in the measurement requires that photons contain many frequencies, i.e., a large bandwidth. Here we challenge this well-known assumption and show that the use of two well-separated frequencies embedded in a quantum entangled state (discrete color entanglement) suffices to achieve great precision. We determine optimum working points using a Fisher Information analysis and demonstrate the experimental feasibility of this approach by detecting thermally-induced delays in an optical fiber. These results may significantly facilitate the use of quantum interference for quantum sensing, by avoiding some stringent conditions such as the requirement for large bandwidth signals.

Equations

$$\begin{array}{r}\delta \tau \ge \frac{1}{2\sqrt{NQ}}=\delta {\tau }_{QCR},\end{array}$$

$$\begin{array}{r}Q=⟨\frac{\partial \Psi (\tau )}{\partial \tau }|\frac{\partial \Psi (\tau )}{\partial \tau }⟩-|⟨\Psi (\tau )|\frac{\partial \Psi (\tau )}{\partial \tau }⟩{|}^2\end{array}$$

$$\begin{array}{rl}|\Psi (\tau )⟩=& \frac{1}{\sqrt{2}}\int d\Omega f(\Omega )\times \\ & [e^{i(\Delta +2\Omega )\tau }{a}_1^{†}({\omega }_1^0+\Omega )a_2^{†}({\omega }_2^0-\Omega )-\\ & a_1^{†}({\omega }_2^0+\Omega )a_2^{†}({\omega }_1^0-\Omega )]|vac⟩,\end{array}$$

$$\begin{array}{r}\delta {\tau }_{QCR}=\frac{1}{N^{1/2}}\frac{1}{({\Delta }^2+4{\sigma }^2{)}^{1/2}},\end{array}$$

$$|\Psi (\tau )⟩\to |{\Psi }_A(\tau )⟩+|\Phi (\tau )⟩,$$

$$\begin{array}{r}{P}_c(\tau )=\frac{1}{2}[1+cos(\Delta \tau +\varphi )exp(-2{\sigma }^2{\tau }^2)],\end{array}$$

$$\begin{array}{rl}{P}_2(\tau )& =\frac{1}{2}(1-\gamma {)}^2[1+\alpha cos(\Delta \tau )e^{-2{\sigma }^2{\tau }^2}]\\ P_1(\tau )& =\frac{1}{2}(1-\gamma {)}^2[\frac{1+3\gamma }{1-\gamma }-\alpha cos(\Delta \tau )e^{-2{\sigma }^2{\tau }^2}]\\ P_0(\tau )& ={\gamma }^2,\end{array}$$

$$\begin{array}{r}\delta {\tau }_{CR}=\frac{1}{(N{F}_{\tau }{)}^{1/2}}\ge \delta {\tau }_{QCR},\end{array}$$

$$F_{\tau }=\frac{({\partial }_{\tau }{P}_2(\tau ){)}^2}{P_2(\tau )}+\frac{({\partial }_{\tau }{P}_1(\tau ){)}^2}{P_1(\tau )}+\frac{({\partial }_{\tau }{P}_0(\tau ){)}^2}{P_0(\tau )}.$$

$$\begin{array}{r}\underset{\tau \to 0}{\lim }{F}_{\tau }={\Delta }^2+4{\sigma }^2.\end{array}$$

$$\begin{array}{rl}0& =:({\partial }_{\tau }log\mathcal{L}{)}_{{\tilde{\tau }}_{MLE}}\\ & =\frac{N_0{P}_0(\tau {)}^{\prime }}{P_0(\tau )}+\frac{N_1{P}_1(\tau {)}^{\prime }}{P_1(\tau )}+\frac{N_2{P}_2(\tau {)}^{\prime }}{P_2(\tau )},\end{array}$$

$$\begin{array}{r}|\psi {⟩}_{\omega }^{-}\to (\frac{|{\omega }_s{⟩}_1|{\omega }_i{⟩}_2-|{\omega }_i{⟩}_1|{\omega }_s{⟩}_2}{\sqrt{2}})\otimes |H{⟩}_1|V{⟩}_2.\end{array}$$

$$\frac{d\beta }{dT}\approx (\frac{2\pi }{{\lambda }_s}-\frac{2\pi }{{\lambda }_i})\frac{dN}{dT}{L}_o+(\frac{2\pi }{{\lambda }_s}{N}_o^{{\lambda }_s}-\frac{2\pi }{{\lambda }_i}{N}_o^{{\lambda }_i})\frac{dL}{dT},$$

$$\begin{array}{r}{\hat{a}}_1^{†}({\omega }_1)=\frac{1}{\sqrt{2}}[{\hat{a}}_3^{†}({\omega }_1)+{\hat{a}}_4^{†}({\omega }_1)]\\ {\hat{a}}_2^{†}({\omega }_2)=\frac{1}{\sqrt{2}}[{\hat{a}}_3^{†}({\omega }_2)-{\hat{a}}_4^{†}({\omega }_2)],\end{array}$$

$$|\Psi (\tau )⟩\to |{\Psi }_A(\tau )⟩+|\Phi (\tau )⟩,$$

$$\begin{array}{rl}|{\Psi }_A(\tau )⟩=& \frac{1}{2}\int d\Omega f(\Omega )(1+e^{-i\Delta \tau })\\ & [{\hat{a}}_3^{†}({\omega }_1^0+\Omega ){\hat{a}}_4^{†}({\omega }_2^0-\Omega )\\ & -{\hat{a}}_3^{†}({\omega }_2^0-\Omega ){\hat{a}}_4^{†}({\omega }_1^0+\Omega )]|vac⟩\\ |\Phi (\tau )⟩=& \frac{1}{2}\int d\Omega f(\Omega )(1-e^{-i\Delta \tau })\\ & [{\hat{a}}_3^{†}({\omega }_1^0+\Omega ){\hat{a}}_3^{†}({\omega }_2^0-\Omega )\\ & -{\hat{a}}_4^{†}({\omega }_1^0+\Omega ){\hat{a}}_4^{†}({\omega }_2^0-\Omega )]|vac⟩.\end{array}$$

$$\begin{array}{r}{F}_{\tau }=\frac{(1-{\gamma }^2)[\alpha \Delta sin(\Delta \tau )+4\alpha {\sigma }^2\tau cos(\Delta \tau ){]}^2{e}^{-4{\sigma }^2{\tau }^2}}{4P_1(\tau )P_2(\tau )/(1-\gamma {)}^4}.\end{array}$$

$$\begin{array}{rl}0& =:({\partial }_{\tau }log\mathcal{L}{)}_{{\tilde{\tau }}_{MLE}}=\frac{N_1{P}_1(\tau {)}^{\prime }}{P_1(\tau )}+\frac{N_2{P}_2(\tau {)}^{\prime }}{P_2(\tau )}.\end{array}$$

$$N_1{P}_2(\tau ){|}_{{\tilde{\tau }}_{MLE}}=N_2{P}_1(\tau ){|}_{{\tilde{\tau }}_{MLE}}.$$

$$\begin{array}{r}{\tilde{\tau }}_{MLE}=\arccos (\frac{\frac{1+3\gamma }{1-\gamma }{N}_2-N_1}{\alpha (N_1+N_2)e^{-2{\sigma }^2{\tau }_s^2}})/\Delta ,\end{array}$$

$$\begin{array}{rl}\beta =& \frac{2\pi }{{\lambda }_s}(\frac{dN}{dT}{L}_o+\frac{dL}{dT}{N}_o^{{\lambda }_s}+\frac{dN}{dT}\frac{dL}{dT}T)T\\ & -\frac{2\pi }{{\lambda }_i}(\frac{dN}{dT}{L}_o+\frac{dL}{dT}{N}_o^{{\lambda }_i}+\frac{dN}{dT}\frac{dL}{dT}T)T.\end{array}$$

Full text

I Introduction

The exploitation of quantum interference promises to enhance sensing technologies beyond the possibilities of classical physics. Hong-Ou-Mandel (HOM) interference is a prototypical example of such a quantum phenomenon, that lacks any counterpart in classical optics. When two identical photons in a global pure state impinge on a beam splitter from separate input modes, they both leave the beam splitter through the same output port, as a consequence of their bosonic nature hong1987measurement. On the other hand, if the input photons are not identical, or they are independent but not in a pure state, the “bunching” probability is directly related to the photons’ level of indistinguishability, or its degree of purity mosley2008heralded. This effect enables a wide range of quantum information processing tasks, ranging from the characterization of ideally identical single-photon emitters aharonovich2016solid, the implementation of photonic Bell state measurements for entanglement swapping and quantum teleportation pan2012multiphoton, or in tailoring high-dimensional entangled states of light zhang2016engineering; ndagano2019entanglement.

HOM interferometry also holds great promise for sensing schemes that require precise knowledge of optical delays. When the relative arrival time of two photons is varied, the coincidence rate exhibits a characteristic dip with a width that is related to the photons’ coherence time. Notably, and unlike other interferometric approaches based on first-order interference, HOM interference is not affected by variations in the optical phase. As a consequence, a HOM interferometer maintains its ability to measure time delays, even when fluctuations of path length difference are on the order of the wavelength. This feature has resulted in proposals for HOM-based time delay sensors with an ultrahigh timing resolution lyons2018attosecond and novel protocols such as Quantum Optical Coherence Tomography (QOCT) that benefit from other quantum features, such as the cancellation of some deleterious dispersion effects nasr2003demonstration.

In the context of such applications, the broad consensus has been that the width of the dip, i.e. the coherence time, imposes the ultimate limit on the precision. As a consequence, ultra-broad-band photon sources have long been hailed as a vital prerequisite for ultra-precise HOM interferometry.

Here we embark on an alternative route towards ultra-precise HOM interferometry using superpositions of two well-separated and entangled discrete frequency modes and coincidence detection on the bi-photon beat note. The manifestation of this fourth-order spatial beating effect is an oscillation within a typically Gaussian envelope that is determined by the coherence time of the two photons, as a direct result of relative phase shift between distinct colors ou1988observation; rarity1990two. We explore the sensitivity limits as a function of the difference frequency of color-entangled states, as imposed by the Quantum Cram$\acute{e}$r-Rao (QCR) bound, and find that the precision with which the delays can be measured is mainly determined not by the coherence time of photons, but by the separation of the center frequencies of the state. Aside from promising improved precision, the approach allows to increase the dynamic range of a HOM-based sensor, provided the required frequency non-degenerate states can be generated in a tunable manner.

We show how suitable frequency entangled states are readily obtained with comparatively little technological effort by employing a variation of the source scheme recently developed in Ref. chen2018polarization. Building on the measurement and estimation strategy by analyzing the Fisher information to determine the optimum working points for frequency-degenerate HOM interference, recently proposed in Ref. lyons2018attosecond, we experimentally demonstrate an optimized HOM sensor that we use to detect delays introduced by temperature drifts in an optical fiber.

The results obtained in this proof of concept experiment show that quantum interference of unconventional frequency states on a beam splitter provides a simple way of enhancing the timing resolution in HOM-based sensors and may also indicate a new direction towards fully harnessing HOM interference in quantum sensing and quantum information processing.

II Results

II.1 HOM interfereometry with frequency entangled states

Let us search for the ultimate limits to the precision of a HOM-based sensor. We consider the generic task of estimating an unknown parameter $\tau$ of a physical system. We prepare a probe state $\lvert\Psi_{0}\rangle$ that is transformed as $\lvert\Psi_{0}\rangle\rightarrow\lvert\Psi(\tau)\rangle$ upon interaction with the physical system. The transformed state is then subjected to a particular measurement strategy to obtain an estimator of $\tau$. Irrespective of the specifics of the final measurement step, we may already state a fundamental limit for the precision of estimation $\delta\tau$ helstrom1969quantum; fujiwara1995quantum:

[TABLE]

where

[TABLE]

and $N$ is the number of independent trials of the experiment. The generality of this statement, known as the Quantum Cram$\acute{e}$r-Rao bound, is remarkable: no matter what ingenious measurement procedure the experimenter may contrive, she will never achieve a precision better than $\delta\tau_{QCR}$. Since the QCR bound is attached to a particular quantum state, it is clear that the appropriate choice of the probe state is of the utmost importance pirandola2018advances.

Let us now consider an experimental configuration where paired photons (signal and idler), with central frequencies $\omega_{1}^{0}$ and $\omega_{2}^{0}$, originate from a parametric down-conversion process (SPDC) pumped by a CW pump with frequency $\omega_{p}=\omega_{1}^{0}+\omega_{2}^{0}$. Each photon of the pair is injected into one of the two arms of a HOM interferometer. The time delay of interest is one that may occur due to an imbalance between the two arms of the interferometer. Even though the common case in HOM interferometry is to consider signal and idler photons with the same central frequency, in the following we allow for a more general configuration where the state of interest is a discrete or continuous frequency entangled state ramelow2009discrete:

[TABLE]

where $\Delta=\omega_{1}^{0}-\omega_{2}^{0}$ is the difference frequency of two well-separated center frequency bins, $\lvert vac\rangle$ is the vacuum state, and $f(\Omega)$ is an the Gaussian spectral amplitude function with $\int d\Omega|f(\Omega)|^{2}=1$. For this state, the Quantum Cram$\acute{e}$r-Rao limit on the estimation of time delays writes

[TABLE]

where $\sigma=\sqrt{\langle\Omega^{2}\rangle-\langle\Omega\rangle^{2}}$ is the RMS (root mean square) bandwidth of SPDC photons. The dependence of the QCR bound on frequency detuning $\Delta$ gives us a first indication to the potential use of non-degenerate frequency entanglement as an alternative to large bandwidth for enhanced resolution HOM interferometry.

Up until now, we have only considered limitations that are inherent to the particular choice of the quantum state. We must confirm that we can experimentally realize this potential benefit using an appropriate measurement strategy, i.e. one that allows us to saturate equation (4). As we shall see in the following, this can be accomplished via coincidence detection in the output ports of a balanced beam splitter. The beam splitter transforms the bi-photon state (see Methods for details) to

[TABLE]

where $\lvert\Psi_{A}(\tau)\rangle$ and $\lvert\Phi(\tau)\rangle$ correspond to the events that two photons emerge in opposite and identical outports, respectively. The normalized coincidence detection probability $P_{c}(\tau)=|\langle\Psi(\tau)|\Psi_{A}(\tau)\rangle|^{2}$ reads

[TABLE]

where $\phi$ is a relative phase factor.

In the case of a real HOM interferometer, that is subject to photon loss $\gamma$ and imperfect experimental visibility $\alpha$, there are three possible measurement outcomes; either both photons are detected, one photon is detected, or no photon detected. The corresponding probability distributions read

[TABLE]

where subscripts 0, 1 and 2 denote the number of detectors that click, corresponding to total loss, bunching and coincidence, respectively. For a more detailed discussion refer to lyons2018attosecond. The outcome probabilities in this measurement can now be used to construct an estimator for the value of $\tau$.

An estimator $\tilde{\tau}$ is a function of the experimental data that allows us to infer the value of the unknown time delay using a particular statistical model for the probability distribution of measurement outcomes. It is thus itself a random variable, that can be constructed from the probability distributions $P_{i}(\tau)$ as a function of time delay. The average of an unbiased estimator corresponds to the real time delay. For any such estimator, classical estimation theory states standard deviation is lower bounded by

[TABLE]

where the Fisher information $F_{\tau}$ reads

[TABLE]

This limit is known as the Cram$\acute{e}$r Rao bound. It is tied to a particular quantum state and a specific measurement strategy. Evaluating the Fisher information for this set of probabilities, we find that its upper bound is achieved in ideal case ($\gamma=0$, $\alpha=1$) at position of $\tau\rightarrow 0$ as

[TABLE]

In the case of zero loss and perfect visibility we recover the Quantum Cram$\acute{e}$r Rao Bound, thus confirming that the measurement strategy is indeed optimal.

While equation (8) provides an ultimate bound on the achievable precision of estimation that can be achieved, the approach does not yet tell us how to construct a suitable estimator for $\tau$. To this end a widely used analytical technique is maximum likelihood estimation (MLE). The likelihood function $\mathcal{L}(\tau)$ is defined from measurement outcomes, whose logarithm can be maximized by using optimization algorithm such as Gradient Descent to predict the parameter $\tau$ that we want to infer. In our framework, the likelihood function is a multinomial distribution as $\mathcal{L}(N_{0},N_{1},N_{2}|\tau)\propto P_{0}(\tau)^{N_{0}}P_{1}(\tau)^{N_{1}}P_{2}(\tau)^{N_{2}}$, where $N_{0}$, $N_{1}$ and $N_{2}$ denote the numbers of events that no, only one and two detector(s) click(s), respectively. Note that $P_{0}(\tau)$, being independent of $\tau$, results in a constant scale factor that is of no relevance to the final calculation of the Fisher information and parameter estimation. The likelihood is extremized as lyons2018attosecond:

[TABLE]

and solving this equation enables us to predict an optimal estimator as $\tilde{\tau}_{MLE}$.

II.2 Experiment

Experimental realization of bi-photon beat note.

We generate photon pairs via spontaneous parametric down-conversion pumped with a continuous-wave pump laser. The experimental configuration implemented to generate the desired frequency entangled state of distant frequency modes (i.e. signal and idler frequencies that are separated by more than the spectral bandwidth $\omega_{s}-\omega_{i}\gg\Delta\omega$) is a modified crossed-crystal configuration kwiat1999ultrabright; steinlechner2012high shown in the inset of Fig. 1. In this configuration, two nonlinear crystals for type-II SPDC are placed in sequence, whereby the optical axis of the second crystal is rotated by 90∘ with respect to the first. Balanced pumping of the two crystals ensures equal probability amplitudes for SPDC emission $\lvert V,\omega_{p}\rangle\rightarrow\lvert V,\omega_{s}\rangle\lvert H,\omega_{i}\rangle$ in the first-, or $\lvert H,\omega_{p}\rangle\rightarrow\lvert H,\omega_{s}\rangle\lvert V,\omega_{i}\rangle$ in the second crystal, where $H/V$ denote horizontal and vertical polarizations. The photons are guided to a PBS, which maps the orthogonally polarized photon pairs into two distinct spatial modes ($1,2$) in the desired frequency entangled state

[TABLE]

The frequency entangled photons are routed to the input ports of a beam splitter. After operation of HOM interference, we only focus on the situation that two detectors indiscriminately register coincidence, i.e., exiting via different ports, as a direct consequence of anti-bunching effect of photons entangled in the form of anti-symmetric state.

As the central wavelengths of down-converted photons are related to the phase-matching temperature of nonlinear crystals, our source has the ability to produce color tunable frequency entangled photon pairs. We analyze the HOM signal for various frequency detunings to demonstrate this flexibility (see Fig. 2). By fitting these interference fringes to normalized coincidence probability as equation (6), we are able to estimate single photon frequency bandwidth to be $0.253\,\mathrm{THz}$, which corresponds to a bandwidth in wavelength of $0.55\,\mathrm{nm}$ and a coherence time of $3.5\,\mathrm{ps}$. Frequency detunings are much larger than single photon bandwidth such that two frequency bins could be separated completely. The visibilities of these experimentally measured frequency entangled photon pairs can reach $0.85\pm 0.05$. The maximal frequency detuning we have measured is $17.08\,\mathrm{THz}$ at temperature of $100\,\mathrm{{}^{\circ}C}$, which is about 68 times the single photon frequency bandwidth.

Fisher information analysis.

Figure 3 demonstrates the explicit procedure of parameter estimation and their corresponding Fisher information in experiment, from which we see that frequency detuning can facilitate the achievement of higher resolution and precision. The oscillation of Fisher information within two-photon coherence time is a key signature of discrete frequency entanglement ramelow2009discrete. Here the maximal Fisher information we have obtained is $245\,\mathrm{ps^{-2}}$ for frequency detuning of $5.34\,\mathrm{THz}$, which inversely reveals the highest precision of $639\,\mathrm{as}$, i.e., relative path delay of $192\,\mathrm{nm}$, for experimental trials of $O(10^{4})$. It is noticed that the quadratic dependence of Fisher information as a function of frequency detuning could be used to further enhance the Fisher information with respect to the frequency degenerate case, where values of $\sim 8\,\mathrm{ps^{-2}}$ have already been reported lyons2018attosecond.

Experimental application as a temperature sensor.

In order to demonstrate the viability principle of employing our HOM sensor, we performed a proof of concept experiment in which we estimate the time delay due to linear expansion of a jacket optical fiber. In order to verify the conclusion that quantum metrology based on frequency entanglement with larger frequency detuning has higher precision, we experimentally measure two-fold coincidence probabilities and predict the thermal coefficient by heating the sensing fiber to vary relative phase shifts (see Fig. 4).

In principle, the relative phase shift varies almost linearly with fiber length and is described as $\beta=N_{g}kL$, where $L$ is sensing fiber length, $N_{g}$ is the material group index and $k$ is the light wave number lagakos1981temperature. Since the input frequency entangled state of HOM sensor is highly sensitive to transmission time, the relative phase shift, introduced by the length extension of fiber, can be expressed as a function of heating temperature, and resulting in the thermal coefficient as

[TABLE]

where $\lambda_{s/i}$ is center wavelength of signal or idler photons, $N_{o}^{\lambda_{s/i}}$ and $L_{o}$ are the corresponding parameters at room temperature, $T$ is heated temperature, $\frac{dN}{dT}$ and $\frac{dL}{dT}$ are thermal coefficients of material group index and fiber length, respectively.

We notice that the thermal coefficient of shifted phase is related to frequency detuning, which agrees well with the experimental measurement results (see Fig. 4), and results in the coincidence probability varies as cosine function (see Fig. 4). The measured thermal coefficients is $0.13\,\mathrm{rad/deg}$, $0.2\,\mathrm{rad/deg}$, $0.3\,\mathrm{rad/deg}$ and $0.48\,\mathrm{rad/deg}$ for frequency detunings of $3.7\,\mathrm{THz}$, $7.4\,\mathrm{THz}$, $11.2\,\mathrm{THz}$ and $17.1\,\mathrm{THz}$, respectively. The refractive index of pure silica is wavelength dependent, and its first derivative with respect to temperature is about $1\times 10^{-5}\,\mathrm{/deg}$ bruckner1970properties. Then we are able to estimate the thermal coefficient of linear expansion of jacket optical fiber to be $\frac{dL}{dT}\sim 4.8\times 10^{-7}\,\mathrm{m/deg}$, which agrees well with the results reported in Ref. tateda1980thermal; priest1997thermal. Accordingly the maximal frequency detuning that we observed in this proof-of-principle experiment enables us to achieve temperature resolution of $0.12\,\mathrm{deg}$.

III Discussion

We have demonstrated a new approach to HOM interferometry based on discrete frequency entanglement of well separated frequency modes and detection of a beat note coincidence signal.

Previous HOM-interferometric sensing schemes required perfect frequency degenerate and ultra-broad-band SPDC emission. Any wavelength distinguishability decreases visibility of the HOM dip and correspondingly diminishes the resolution. Providing suitable quantum sources for this case is a significant challenge, as it either requires the engineering of aperiodic poling structures or the use of very short nonlinear crystals, at the cost of efficiency. In contrast, the approach outlined here requires only a sufficiently large non-degeneracy, whereby the spectral bandwidth can be small. We have experimentally demonstrated how to generate suitable discrete frequency-entangled states, in a manner that can be readily extended to larger wavelength separations. For example, $\lambda_{s}$= $1,500\,\mathrm{nm}$ and $\lambda_{2}$ = $800\,\mathrm{nm}$ ($1,000\,\mathrm{THz}$ angular difference frequency) a timing sensitivity of 9 as could already be obtained for only $N=10^{4}$ detection events. Backed by the results of our proof-of-concept experiment, this shows that the approach can provide higher resolution and highly sensitive measurement, and makes it an ideal candidate for more quantum enhanced metrology applications.

Although this work only reports the advantages of our approach in estimating delays, similar great enhancement can also be achieved for a variety of applications like state discrimination or hypothesis testing.

In conclusion, we believe that fully harnessing HOM interference and frequency entanglement will provide additional tools, e.g. for frequency shaping of photons and interference phenomena in general, ultimately broadening the path towards practical quantum applications.

IV Methods

Entangled Photon source. In our experimental realization of flexible frequency entanglement source chen2018polarization, two mutually orthogonally oriented 10-mm-long ppKTP crystals are manufactured to provide collinear phase matching with pump (p), signal (s) and idler (i) photons at center wavelengths of $\lambda_{p}\approx 405\,\mathrm{nm}$ and $\lambda_{s,i}\approx 810\,\mathrm{nm}$. They are pumped with a $405\,\mathrm{nm}$ continuous wave grating-stabilized laser diode. To achieve the desired diagonally and anti-diagonally polarization states for simplifying alignment, we design an oven with V-groove such that two crossed crystals are oriented at $45\,\mathrm{{}^{\circ}}$. Since the pump beam is horizontally polarized, it is equally likely to generate a photon pair in the first or second crystal, resulting in a state of equation (12). The relative phase factor is compensated by tilting a half wave plate. Long pass filter is set to block pump beam. Then PBS routes a pair of photons into two distinct spatial modes according to orthogonal polarizations.

HOM interferometer. In spatial mode 1, a translation stage introduces a relative path delay to accomplish the task of scanning HOM interference fringes. Polarization controllers are required to compensate polarization difference of biphotons such that only frequency entanglement can make contributions to the interference effect. Finally the anti-bunched photons are detected by silicon avalanche photon diodes, and two-fold events are identified using a fast electronic AND gate when two photons arrive at the detectors within a coincidence window of $\sim 3\,\mathrm{ns}$.

Coincidence signal with frequency-entangled states. An optimal measurement procedure may allow us to saturate the limit set by equation (8). In the ideal, lossless with perfect visibility, case, such a measurement can be accomplished by interference on a balanced beam splitter. The beam splitter transformation on the input modes can be expressed by

[TABLE]

where $\omega_{1}$ and $\omega_{2}$ denote the signal or idler frequency mode that are incident from opposite ports, and subscripts 1/2 (3/4) represent two input (output) ports of that beam splitter. Accordingly the state is transformed as

[TABLE]

where these state contributions can be expressed as

[TABLE]

Due to HOM interference on the beam splitter coincidence detection in distinct spatial modes projects onto the state component $\lvert\Psi_{A}(\tau)\rangle$.

Fisher information. In a specific experiment (measurement strategy), with outcomes $x_{i}$, and corresponding probability distributions $P_{i}(\tau)$, any unbiased estimator will fulfill equation (8), where the Fisher information $F_{\tau}$ quantifies the information that a particular measurement can reveal about the unknown parameter of interest. Note that optimizing over all probability distributions results we recover the QCR bound. The outcomes of this measurement are sufficient to obtain an estimator for the value of $\tau$.

By substituting equation (9) with the corresponding probabilities from equation (7), we could calculate the Fisher information as

[TABLE]

We note that the Fisher information is undefined at position of $\tau=0$ in ideal case since the denominator will be zero.

Maximum-likelihood estimator. Since no prior knowledge is provided, we can apply maximum likelihood estimation approach to predict the target parameter. We extremized the likelihood function as

[TABLE]

Based on the calculation in equation (7), we know $P_{1}(\tau)^{\prime}=-P_{2}(\tau)^{\prime}$ such that

[TABLE]

For the sake of simplicity, $\tilde{\tau}$ in term of $e^{-2\sigma^{2}\tau^{2}}$ can be considered as a constant value, i.e., coarse sensing position $\tau_{s}$ where Fisher information is highest. Thus we get an optimal estimator to variable relative time delay as

[TABLE]

and the values of parameters $\tau_{s}$, $\gamma$, $\alpha$, $\sigma$ and $\Delta$ need to be separately estimated before the measurements begin.

Shifted phase of temperature sensor. The introduced phase shift in frequency entangled state can be expressed as a function of heating temperature as

[TABLE]

Since $\frac{dN}{dT}\frac{dL}{dT}$ is in the order of much smaller magnitude, we ignore the term of $\frac{dN}{dT}\frac{dL}{dT}T$ in equation (21).

See supplementary materials for more information.

Data availability

Data available on request from the authors.

Acknowledgements

We thank Thomas Scheidl, Sebastian Ecker, Soeren Wengerowsky, Johannes Handsteiner, Siddarth Joshi, and Lukas Bulla for experimental support and helpful conversations. YC thanks Lijun Chen for support. The research leading to these results has received funding from the H2020 European Programme under Grant Agreement 801060 Q-MIC, the Austrian Research Promotion Agency (FFG) Projects - Agentur für Luft- und Raumfahrt (FFG-ALR contract 6238191 and 866025), the European Space Agency (ESA contract 4000112591/14/NL/US) as well as the Austrian Academy of Sciences. YC acknowledges personal funding from Major Program of National Natural Science Foundation of China (No. 11690030, 11690032), National Key Research and Development Program of China (2017YFA0303700); the National Natural Science Foundation of China (No.61771236), and from a Scholarship from the China Scholarship Council (CSC) and the program B for Outstanding PhD candidate of Nanjing University. This work was supported by the Fraunhofer Internal Programs under Grant No. Attract 066-604178. JPT acknowledges financial support from Fundacio Cellex, from the Government of Spain through the Severo Ochoa Programme for Centres of Excellence in R&D (SEV-2015-0522), from Generalitat de Catalunya under the programs ICREA Academia and CERCA, and from the project 17FUN01 BeCOMe within the Programme EMPIR, and initiative co-founded by the European Union and the EMPIR Participating Countries.

Contributions

F.S. developed the initial idea for this work. Y.C. conducted the experiment under supervision from F.S. and R.U. Theoretical analysis was carried out by J.T. and Y.C.. Y.C. and F.S. wrote the first draft and all authors contributed to the final version of the manuscript.

Competing interests

The authors declare that there are no competing interests.