
2
on the target frequency, f, in question and its interac-
tion with the sensor as expressed in the Hamiltonian,
H(f). This pulse sequence will use the same exponen-
tially growing sensing time as in the quantum PEA only
in sequential order, from the shortest to the longest, and
not simultaneously, similar to the Kitaev’s iterative PEA
[22]. After each pulse sequence with one sensing time,
the sensor is measured, and the outcome (u) can be one
of the two states of the sensor - zero or one. This out-
come is used in the second step to update the probability
function, Pm(u|f), to measure the sensor state, |0ior |1i,
given that there is interaction with the unknown param-
eter with frequency f. The nature of the sensor’s inter-
action with the target parameter in the pulse sequence is
encoded in the probability function.
The critical step of the algorithm is in step 3 , where
one applies a Bayesian update to estimate the unknown
parameter [23–25]
Pposterior(f|u)∝Pm(u|f)Pprior(f|u) (2)
where P(f|u) is the probability function of the mea-
surement outcome given the target parameter, subscript
posterior is the new probability after each Bayesian
update and prior is the old one from the last update.
P(u|f) is the probability function of the target parame-
ter given the outcome of the measurement is u, the sub-
script mdenotes a single outcome. Since the adaptive
1
2
4
3
uniform distribuon
Measurement
Calculate
Adapt variables
Bayesian update
Threshold
Single shot readout
number of posive results
Threshold
number of posive results
non single-shot readout
FIG. 1. (a) Graphical illustration of the adaptive phase esti-
mation algorithm comprising four steps: (1) A pulse sequence
suitable for the estimation of the unknown parameter, given
the nature of the interaction between the sensor and the pa-
rameter. This pulse sequence will be applied with exponen-
tially growing sensing times. The state of the sensor is mea-
sured after every sensing time. (2) Calculating the proba-
bility function for the state of the sensor given the unknown
parameter. (3) Using Bayes’ Theorem to update the proba-
bility function for the parameter. (4) Calculating the optimal
variables for extracting maximal information from the next
iteration. After Mkiterations for each sensing time, the final
distribution will be the estimate of the unknown parameter.
(b-c) Schematic illustration of the measurement outcome of a
single-shot (b) or averaged (c) sensor.
PEA applies the sensing scheme with different sensing
times sequentially, each measurement holds less infor-
mation about the phase than the quantum PEA. The
penalty in the full scheme is that each sensing time is
measured multiple times by changing one of the sensing
variables, as is illustrated in step 4 . The number of it-
erations Mk=G+ (K−k)Ffor each sensing time grows
as the sensing time gets shorter, where Gand Fare op-
timized parameters, and kis the index of the sensing
time [26]. The adaptive character of the scheme is es-
tablished in step 4 . In this step, the optimal variables
for gaining maximal information are calculated based on
the last probability function and then transferred to the
pulse sequence of the next iteration.
Adaptive PEA has been studied extensively. Theoret-
ical works suggested controlling the sensing phase or the
sensing time [27] to enhance sensitivity. Others used nu-
merical simulations [28, 29], and several did experimental
studies with different sensors [23, 30] to prove the feasi-
bility and benefits of this protocol. All of these studies
were performed with a single-shot readout (SSR) sen-
sor, where the state of the sensor can be measured after
one measurement with high fidelity (Fig. 1b). Neverthe-
less, in some cases, non-SSR sensors are the only possible
sensing approach, for instance, for imaging nanoscale bi-
ological samples with high special resolution and in am-
bient conditions. These sensors are characterized by the
high ratio of classical noise added in the measurement,
for example, low photon collection efficiency in optically
read-out systems, compared to the quantum projection
noise of the system [11]. This causes the histogram of
the measurement outcomes to mix ‘0’ and ‘1’. Therefore,
assigning the measurement outcome to one state of the
sensor with high fidelity, i.e., in one shot, is impossible
(Fig. 1c).
For a non-SSR sensor, the pulse-sequence and the mea-
surement should be applied for many repetitions to assess
the sensor state, still with a non-negligible error. This sit-
uation requires adjusting the probability function P(u|f)
used in the Bayesian update to the averaged measure-
ment result. The most common and simple solution is
to use a threshold that is calculated based on the prob-
ability to measure a positive outcome from the sensor at
each of the states, which can be a collection of photons
for an optically measurable sensor (See Appendix V E).
In this method the measurement is repeated for Rtimes
and the number of positive outcomes ris assigned to a
state of the sensor, u, based on the calculated threshold;
we call this method “majority voting”. This approach
results in a binary outcome from a large batch of size R
repetitions of the measurement. This method’s benefit
is the possibility of using the probability function and
Bayesian update as in the SSR sensor scheme [31]. How-
ever, it disregards most of the possible outcomes from
the Rrepetitions by using only a binary span of results.
Therefore, it is also more prone to noise, where a noisy
measurement can be assigned to the wrong binary option
[32].