Adaptive quantum metrology under general Markovian noise

Adaptive quantum metrology under general Markovian noise

Rafał Demkowicz-Dobrzański Faculty of Physics, University of Warsaw, ul. Pasteura 5, PL-02-093 Warszawa, Poland    Jan Czajkowski Faculty of Physics, University of Warsaw, ul. Pasteura 5, PL-02-093 Warszawa, Poland QuSoft, University of Amsterdam, Institute for Logic, Language and Computaion (ILLC),
P.O. Box 94242, 1090 GE Amsterdam, The Netherlands
Pavel Sekatski Institut fur Theoretische Physik, Universität Innsbruck, Technikerstr. 21a, A-6020 Innsbruck, Austria
Abstract

We consider a general model of unitary parameter estimation in presence of Markovian noise, where the parameter to be estimated is associated with the Hamiltonian part of the dynamics. In absence of noise, unitary parameter can be estimated with precision scaling as , where is the total probing time. We provide a simple algebraic condition involving solely the operators appearing in the quantum Master equation, implying at most scaling of precision under the most general adaptive quantum estimation strategies. We also discuss the requirements a quantum error-correction like protocol must satisfy in order to regain the precision scaling in case the above mentioned algebraic condition is not satisfied. Furthermore, we apply the developed methods to understand fundamental precision limits in atomic interferometry with many-body effects taken into account, shedding new light on the performance of non-linear metrological models.

pacs:
03.67a, 03.65Yz, 03.65.Ta, 06.20.-f, 42.50.Lc

I Introduction

With rapid advancements in quantum optical experimental techniques, the field of quantum metrology Giovannetti et al. (2011); Toth and Apellaniz (2014); Demkowicz-Dobrzanski et al. (2015); Pezzè et al. (2016); Schnabel (2016) is entering the stage where ubiquitous quantum features of light and matter are being harnessed to deliver ultra-sensitive measuring devices for real-life applications Schmidt et al. (2005); LIGO Collaboration (2011); Lücke et al. (2011); Taylor et al. (2013); Ockeloen et al. (2013); Kovachy et al. (2015); Barrett et al. (2016). Along with experimental advances, theoretical foundations for the field have been constantly developed. From the first proposals of utilizing squeezed states in optical interferometry Caves (1981), through identification of fundamental limits in decoherence-free metrology Braunstein (1992); Berry and Wiseman (2000); Giovannetti et al. (2006), general methods have been developed allowing to take into account the impact of realistic decoherence effects on the performance of metrological protocols Huelga et al. (1997); Fujiwara and Imai (2008); Genoni et al. (2011); Escher et al. (2011); Knysh et al. (2014); Demkowicz-Dobrzański et al. (2012, 2013); Jarzyna and Zwierz (2017) including the most general quantum adaptive strategies Demkowicz-Dobrzański and Maccone (2014); Sekatski et al. (2016a).

Most of the available general methods is based on the integrated form of the dynamics of a quantum system represented mathematically as a quantum channel Escher et al. (2011); Demkowicz-Dobrzański et al. (2012); Demkowicz-Dobrzański and Maccone (2014). This poses a serious difficulty when the dynamics is provided in terms of a Master differential equation. In this case obtaining the analytical form of the integrated dynamics is often impossible. This fact significantly limits the utility of the available methods making it often necessary to resort to numerical calculations instead of a more insightful analytical analysis. This deficiency has been successfully addressed in case of single qubit dynamics, where the full description of performance of the most general quantum metrological protocols has been given Sekatski et al. (2016a). In particular, it has been shown that provided the noise is represented by a single Pauli operator which is not proportional to the Hamiltonian itself, one can apply an error correction procedure allowing to reach the Heisenberg-like, , scaling of Quantum Fisher Information (QFI), where is the total evolution time of the probe system. On the other hand, for all other kind of Markovian noise processes the optimal QFI is limited by a classical-like scaling bound proportional to and hence results in a standard scaling of precision.

This paper provides a general solution to the problem of determining optimal performance of adaptive metrological schemes in a unitary parameter estimation problem for arbitrary Markovian dynamics. We present an explicit recipe that allows to obtain the formulas for the behavior of QFI in the optimal metrological protocol, based solely on the operators appearing explicitly in the Master equation in the standard Gorini-Kossakowski-Lindblad-Sudarshan Breuer and Petruccione (2002) form, with no need to integrate the dynamics whatsoever. The probe dynamics we consider is given by

 dρdt=−iω[H,ρ]+J∑j=1LjρL†j−12ρL†jLj−12L†jLjρ, (1)

where is the frequency-like parameter to be estimated associated with the unitary dynamics generated by the Hamiltonian , while are noise operators. In particular, we show that if

 (2)

where denote the hermitian and the anti-hermitian part of an operator, then the QFI scales at most linearly with , and the coefficient for the bound can be obtained from a solution of a simple semi-definite program. When restricted to a single qubit problem, this condition is equivalent to the one given in Sekatski et al. (2016a) requiring the noise not to be a single-rank Pauli linearly independent from the Hamiltonian. If the above linear dependence condition is not satisfied we discuss the possibility of implementing a “quantum error-correction”-like protocol that yields quadratic scaling of QFI in . In the qubit case this is always possible Sekatski et al. (2016a) using a scheme based on preparing a maximally entangled state of the probe system and an equally dimensional ancilla. Here we demonstrate that in higher dimensions this is in general no longer the case, and this approach is not always sufficient to overcome the effects of noise.

Further on we apply the newly developed quantitative methods to determine fundamental precision bounds in atomic metrological protocols involving many-body interactions. This allows us in particular to derive for the first time fundamental precision bounds on non-linear metrological protocols in presence of decoherence. In absence of decoherence it is known that in the case of the -body Hamiltonian non-linearity may help to improve the precision scaling of QFI to where is the number of atoms involved Boixo et al. (2007); Luis (2007); Boixo et al. (2008); Choi and Sundaram (2008); Napolitano and Mitchell (2010); Gross et al. (2010); Hall and Wiseman (2012); Joo et al. (2012); Sewell et al. (2014). We show that in the case of the -body Hamiltonian and -body noise the linear dependence condition implies QFI to scale no better than —a scaling formula identified in Beau and del Campo (2016) but only for GHZ states a and limited class of noise models. Apart from determining the scaling character of the bounds, we also provide explicit coefficients for the bounds in case of linear and non-linear atomic interferometry in presence of single and two-body losses. Note that we focus here on unitary parameter estimation in presence of noise and do not analyze the problem of estimating the noise parameter itself. This last problem, while interesting, does not enjoy equally spectacular quantum gains thanks to the use of entangled states as the unitary parameter estimation case. Often a completely uncorrelated state proves to be optimal, as for example in the problem of estimating losses or dephasing strength, while in other cases entanglement between a single probe and a passive ancilla is sufficient to reach optimality Hotta et al. (2005, 2006); Kołodyński and Demkowicz-Dobrzański (2013); Takeoka and Wilde (2016); Pirandola and Lupo (2017)

Ii Formulation of the problem

Considering the Master equation given in Eq. (1), let us denote by the integrated form of the dynamics so that

 ρωt=Eωt(ρ0)=∑jKωt,jρ0Kω†t,j(t), (3)

where are Kraus operators of the evolution.

The aim is to perform optimal estimation of parameter under the constraint of a fixed total evolution time under the most general adaptive quantum metrological scheme as depicted in Fig. 1.

Given the final state of the protocol , the fundamental limitation on the precision of estimating is given in terms of quantum Cramér-Rao bound:

 Δω≥1√FQ, FQ=2∑ab|⟨a|˙ϱωT||b⟩|2λa+λb (4)

where is QFI for unitary encoding, dot signifies , while , are eigevectors and eigenvalues of . In what follows we will use QFI as the figure of merit. Direct maximization of QFI of the final state over all tunable elements in the protocol, i.e. input state and controls is a virtually impossible task unless a decoherence-free case is considered where adaptiveness is useless Giovannetti et al. (2006).

Luckly, provided the integrated form of the dynamics in the form of (3) is given, one can apply the methods from Demkowicz-Dobrzański and Maccone (2014); Sekatski et al. (2016a) that allow to obtain a universally valid upper bound on QFI valid for arbitrary adaptive strategy, and hence a lower bound on uncertainty. The bound utilizes the observation that given a quantum evolution in the form of Kraus representation (3) one can consider equivalent Kraus representation consisting of operators related by a unitary matrix with the original ones, or written in a more concise notation , where .

The bound on QFI for the step adaptive protocol is then given in terms of minimization over different Kraus representations and reads:

 FQ≤4min~Kω,x{m∥α∥+m(m−1)∥β∥(x∥α∥+∥β∥+1/x)}, (5)

where

 α:=˙~K†˙~K,β:=−i˙~K†~K, (6)

is the operator norm and is a positive real parameter minimization over which helps to further tighten the bound. The crucial element here is that the unitary matrix can explicitly depend on the estimated parameter . Using the fact that for some hermitian matrix , one notes that this dependence enters the computation of and via the substitution:

 ˙~K=˙K−ihK,~K=K. (7)

In the most general adaptive approach it is always advantageous to take the limit as in this case the use of adaptive gates potentially provides the greatest benefits. Note that the limit is taken only in the duration of the sensing map and it does not influence the duration of the intermediate unitary control gates. In fact in our model we assume that the time needed to perform control gates is not included in the total resource count and hence the continuous limit does not affect the control gates time. At the same time this limit allows us to have arbitrary many control gates over the course of the sensing process and hence results in the most general adaptive strategy.

In order to keep this in mind, from now on we will therefore replace with . Taking this limit we may now utilize known relations between noise operators and the corresponding Kraus operators in order to write an explicit Kraus representation for the dynamics in the lowest order in :

 K0 =\openone−(12L†L+i% Hω)dt+O(dt2), (8) Kj =Lj√dt+O(dt32) (j=1,…,J) (9)

where . Note that Kraus operators’ index starts from while the noise opeartors’ index starts from . This approximation correctly recovers the dynamics in the linear order in and captures all the features of Markovian evolution.

Because of the different scaling appearing in and operators it will be convenient to introduce the following structure of the matrix :

 h=⎡⎢ ⎢ ⎢ ⎢ ⎢⎣h00h†hh⎤⎥ ⎥ ⎥ ⎥ ⎥⎦. (10)

Minimization over different Kraus representations in (5) amounts now to minimization over the matrix and taking into account that we consider the limit we get

 FQ≤4minh,x{T∥α∥dt−1+T2∥β∥dt−2(x∥α∥+∥β∥+1/x)}. (11)

The most interesting challenge posed be the above formula is to determine the situation where is limited by a bound scaling linearly in or where the bound scales as , in which case achieving the Heisenberg scaling may be possible.

Iii H∈S: T scaling of QFI

The bound will scale linearly in if we are able to find which makes as well as , since then by choosing we will get

 FQ≤4∥α(1)∥T (12)

in the limit . Let us write explicitly in terms of , , in leading orders in :

 β=Hdt+h00(\openone−L†Ldt)+(h†L+L†h)√dt+L†hLdt+O(dt32). (13)

Let us denote time expansion coefficients of as follows: . We now investigate the condition order by order in time. Making in orders and implies that as well as . The first non-trivial condition appears in order as this is the order where the Hamiltonian contributes and setting is not enough to get . With the above substitutions we may write the linear order coefficient of :

 β(1)=H+h(1)00\openone+h†(12)L+L†h(12)+L†h(0)L. (14)

Taking into account that is hermitian, this coefficient can be made zero if and only if the Hamiltonian , where subspace is defined in (2). Analysing the next order we get:

 β(32)=h(32)00\openone+h(1)†L+L†h(1)+L†h(1/2)L. (15)

We see that none of the coefficients appearing here appeared before when considering so we may put them all equal zero guaranteeing that , and proving the linear scaling of QFI.

In order to obtain a quantitative bound, as given in Eq. (12), let us now focus on the operator . Taking into account that as well as , the first nontrivial order is linear in and the corresponding coefficient reads:

 α(1)=(h(12)\openone+h(0)L)†(h(12)\openone+h(0)L). (16)

We now need to minimize the operator norm of the above coefficient over in order to get the tightest bound:

 FQ≤4Tmin{h(1)00,h(12),h(0)}∥α(1)∥,subject to: β(1)=0. (17)

Only in some particular cases this can be done analytically. Fortunately, the above problem can be implemented as a semi-definite program, as described explicitly in Appendix A. The implementation is similar to the one presented in Demkowicz-Dobrzański et al. (2012) for the Kraus operator formulation.

Since the bound (12) involves the operator norm it may not be immediate to apply it in infinite dimensional cases where the operators appearing in the Master equation are unbounded. This is for example the case when dealing with continuous variable systems. Following the original derivation of the bound (5), however, one can refine it by taking into account some properties of the state utilized in the protocol. In particular it might be that the states we deal with are restricted to consist of finite number of photons/atoms, or at least have the mean number of particles fixed. It might also, be the case that by various super-selection rules the whole Hilbert space is not available and the bound can be tightened by analysing the operator norm of separately in different super-selection sectors. This will amount to calculation of operator norms on subspaces. Moreover, provided sufficient information on the time evolution of the probe state is given, the bound may even be formulated as an integral over interrogation time of a time-dependent quantity. Namely

 FQ≤4∫T0min⟨α(1)⟩tdt,% subject to: β(1)=0, (18)

with setting a limit on the gain in QFI at a given instant of time, and is the state of the system at time — see Appendix B for the details how to tighten the bound in these cases.

In what follows, it will be convenient to assume that the Master equation (1) is given in the so called canonical form Pearle (2012); Hall et al. (2014), where all noise operators are traceless and orthogonal , see also C.1.

Let us now briefly comment on the single qubit case discussed in Sekatski et al. (2016a). Since operators in the canonical form are traceless there can be written as complex combinations of Pauli matrices. The condition required to get a better than linear scaling of QFI discussed there was that the noise is a single-rank Pauli linearly independent from the Hamiltonian. Mathematically this means that there is only one . Note that if we had two linearly independent , they would lead to being the full space of hermitian matrices, thanks to the fact that products appear in the definition of . Moreover, even with a single which is not proportional to a hermitian matrix, we would have two linear independent hermitian matrices from its hermitian and anti-hermitian part and hence again the generated would be equal to the whole space of hermitian matrices. Consequently, in the qubit case, the condition is equivalent to the one discussed in Sekatski et al. (2016a).

Iv H∉S: Possibility of T2 scaling of QFI

If and hence cannot be made zero then the second term in the bound (11) will not vanish (due to scaling as and canceling with term) and will result in an upper bound scaling as . This gives hopes for the Heisenberg scaling of precision. Below we discuss the possibility to construct an adaptive quantum error correction inspired strategy that allows to achieve a scaling of QFI and discuss some concrete examples.

Aside the probe system we allow for an additional ancillary system on which the evolution acts trivially. Let denote the input probe-ancilla state. The adaptive protocol we consider consists of intertwining of infinitesimal-time probe evolution

 Eωdt(ϱ)=ϱ+(−iω[H,ϱ]+L(ϱ))dt+O(dt2), (19)

where represents the noise part of the Master equation (1), and the error correction map applied after each infinitesimal time step . Hence, the final state of the probe and ancilla systems after the total evolution time , i.e. after applications of the evolution-correction step, reads:

 ϱωT=CωT(ϱ)=(C∘Eωdt)∘Tdt(ϱ). (20)

Formally, in the above formula we should write instead of , as the map acts also on the ancillary system in a trivial way. From now on, for simplicity of notation, we will assume that whenever an operator defined on the probe system alone acts on the extended probe-ancilla system it should be understood as extended in a trivial way.

To simplify the exposition we assume that estimation of is made around point. Otherwise by means of active control one can always compensate for the nonzero rotation term in the master equation—note that in this case the error correction operation may in general depend on . Since QFI depends on the state and its first derivative at , see Eq. (4), it is enough to consider the first order expansion of the final state in the estimated parameter : . The zeroth order is simply the action on the input state of the dynamics where the Hamiltonian part is dropped, while the first order term

 (21)

involves terms of where the Hamiltonian part of the dynamics enters linearly at different times.

Our goal is to design a protocol that protects the system from decoherence in a way that its QFI grows quadratically in and hence mimics the performance of noiseless frequency estimation. First of all, we demand that our protocol preserves the initial state in absence of the unitary evolution . Furthermore, let us define a decoherence-free qubit subspace , spanned by the input state as well as an orthogonal state such that , . The error-correction step projects the state after being acted on with back onto the subspace . We also assume that since otherwise the parameter would not be imprinted on the state at all. This allows us to simplify the the expression in Eq. (21)

 ˙ϱ0T =ic∫T0C0T−t(|ϕ⟩⟨ξ|−|ξ⟩⟨ϕ|)dt. (22)

In addition, we require that the control-assisted evolution preserves the coherence , which leads to . Calculating QFI with eigenvectors and yields quadratic as in the case of noiseless frequency estimation. Otherwise, had the evolution damped the coherence term, resulting in this would not be possible. Hence in summary, the requirement for the error correction map are the following:

 (i)C(|ϕ⟩⟨ϕ|+dtL(|ϕ⟩⟨ϕ|)) =|ϕ⟩⟨ϕ|+O(dt2) (23) (ii)C(|ξ⟩⟨ϕ|+dtL(|ξ⟩⟨ϕ|)) =|ξ⟩⟨ϕ|+O(dt2). (24)

In fact, the linearity and the trace preserving properties of , together with the above conditions, imply . As a result, the requirement for QFI to grow quadratically amounts to the requirement of existence of a two-dimensional subspace protected from decoherence up to the linear order in time. We may therefore utilize known results from approximate quantum error correction literature, which in this case reduce to the standard error-correction relation Knill and Laflamme (1997) for the set of error operators consisting of and the identity operator Bény (2011):

 (a)⟨ϕ|H|ξ⟩≠0, (25) (b)⟨ϕ|L†kLj|ξ⟩=⟨ϕ|Lj|ξ⟩=0, (26) (c)⟨ϕ|L†kLj|ϕ⟩=⟨ξ|L†kLj|ξ⟩, (27)

for all and , where (a) is an additional requirement that needs to be satisfied in order to keep non-trivial unitary evolution in the qubit subspace . Step by step derivation of the above conditions is provided in Appendix C.2.

Following the way the single-qubit error-correction protocols were applied in quantum metrology Arrad et al. (2014); Kessler et al. (2014); Dür et al. (2014); Sekatski et al. (2016a) a natural choice for is the maximally entangled state of probe+ancilla . Recall that in the canonical form of Master equation all are traceless and orthogonal , see C.1. Hilbert-Schmidt orthogonality of is automatically transferred to orthogonality of vectors (where should be in fact understood here as ). We then decompose the Hamiltonian such that while nonzero is orthogonal to all operators in . If we now take (note it is by construction orthogonal to as is in particular orthogonal to the identity operator), then one automatically satisfies the first two conditions. Condition (a) follows from , while (b) follows from

 ⟨ϕ|L†kLjH⊥|ϕ⟩∝TrH⊥L†kLj=0,⟨ϕ|LjH⊥|ϕ⟩=0, (28)

as is orthogonal to with respect to the Hilbert-Schmidt scalar product. For the qubit case Sekatski et al. (2016a) this construction also guarantees condition c) to be satisfied, as implies only one . Let , be the eigenbasis of : ( is orthogonal to and hence has eigenvalues). As a result , and consequently: . In this case one can state that Eq. (2) is an if and only if condition for impossibility of getting QFI scaling quadratically with .

In higher dimensions, however, an error-correction scheme based on the use of the maximally entangled state of system+ancilla will not work in general—see the Note Added at the end of the Conclusions section with a reference to the paper Zhou et al. (2017), where a universal construction of a quantum error-correction protocol satisfying all the required conditions has been provided whenever . This also, shows that the conditions which in our approach could be regarded as sufficient for quadratic scaling are actually also necessary.

V Atomic interferometry with one and two-body losses

We will now demonstrate an application of the developed methods to provide bounds in atomic interferometry models where two-body effects can be placed both in the noise part (two-body losses) or in the Hamiltonian part (the non-linear metrology model).

Let us consider a Bose-Einstein Condensate (BEC) system of two level atoms, where dynamics is described by the following master equation:

 dρdt=−iω[H(k),ρ]+L(1)(ρ)+L(2)(ρ), (29)

where

 L(1)(ρ) =2∑i=1γi(aiρa†i−12{a†iai,ρ}), (30) L(2)(ρ)=2∑i=1γii(a2iρa†2i−12{a†2ia2i,ρ})+γ12(a1a2ρa†1a†2−12{a†1a†2a1a2,ρ}), (31)

represent one-body and two-body loss processes with respective loss coefficients and representing anihilation operators removing an atom from the -th mode. The corresponding noise operators read: , . We will consider two different Hamiltonians that are associated with the sensing part of the dynamics

 H(k=1) =12(a†1a1−a†2a2), (32) H(k=2) =14:(a†1a1−a†2a2)2:, (33)

which correspond to linear and non-linear metrological scenarios. For the clarity of presentation, we have put the normal ordering operation in the definition of in order to make sure that we take into account only terms that appear due to interaction between two different particles.

Let us start with the linear Hamiltonian case , but keep both the single and two-body losses processes. This kind of model is well tailored to analyze matrological BEC experiments such as e.g. magnetometry experiment using spin-squeezed BEC Ockeloen et al. (2013). Let us calculate according to Eq .(14):

 β(1)=12(a†1a1−a†2a2)+h(1)00\openone+(∑i√γih(12)iai+h.c.)+∑ijh(0)ij√γiγja†iaj+…. (34)

We now ask about the possibility of choosing entries of in a way to set the above quantity to zero. Note that we have not written explicitly the noise operators related to two-body losses. The reason for that is that operators related to two-body losses appearing in the above equation would be of the form and would be linearly independent of the operators appearing in the Hamiltonian part. Hence trying to set we need to focus on one-body losses operators only. We are free to put all coefficients of in front of terms related with two-body losses equal to zero. If we succeed in setting using only one-body losses operators this will also imply that two-body losses are irrelevant in trying to assess the fundamental precision limit on frequency estimation in this case. By inspecting Eq.(34) it is clear that we can make by choosing , , , , .

We should remember, however, that when deriving the final bound using Eq. (16) we face the problem that operators appearing under the operator norm are unbounded and hence the bound formally will be infinite and hence useless. Physically this is due to the fact that we have not set any constraints on the number of atoms we use in the experiment. From now on we will assume we have atoms at our disposal and at every adaptive step we replace the lost atoms with fresh ones keeping the number of atoms constant. We discuss this approach in detail in Appendix B and argue that by doing so we do note lose generality of our bounds. Thanks to this, we are able to write when calculating the operator norm of , remembering that in the end we operate in a fixed particle number-subspace.

Let us now go back to Eq. (34). With fixed particle number constraint imposed, , and are no longer independent operators. This gives us an additional freedom in choosing coefficients of in order to keep , namely we can take , , , with being a free parameter. The QFI bound can be now obtained by minimizing over :

 FQ≤Tminξ∥γ−11(2ξ−1)2a†1a1+γ−12(2ξ+1)2a†2a2∥. (35)

Operators , commute and their common basis is (where is the number of atoms in mode ) hence in the above minimization we can replace with and with :

 FQ≤Tminξmax0≤n≤Nγ−11(2ξ−1)2n+γ−12(2ξ+1)2(N−n). (36)

The minimum is achieved for that satisfies and reads:

 FQ≤4TN(√γ1+√γ2)2. (37)

This bound indeed agrees with a known bound of particle interferometry with losses Kolodynski and Demkowicz-Dobrzanski (2010); Knysh et al. (2011); Escher et al. (2011); Demkowicz-Dobrzański et al. (2012)

 FQ≤4TtN(√1−η1η1+√1−η2η2)2, (38)

where after taking the limit . Note, however, that the derivation presented in this paper did not require any educated guess Escher et al. (2011) nor numerically indicated optimal form of Kraus representation Demkowicz-Dobrzański et al. (2012), but resulted in a purely algebraic analysis of the noise operators and the Hamiltonian appearing in the Master equation. Moreover, when deriving the bound we could clearly see that two-body losses do not have an impact on the bound.

We now move on to study the fundamental bounds in a non-linear metrological model with , in which case

 β(1)=14(a†21a21+a†22a22−2a†1a†2a1a2)+h(1)00\openone+h(0)11,11γ11a†21a21+h(0)22,22γ22a†22a22+h(0)12,12γ12a†1a2†a1a2+…, (39)

where we explicitly wrote only terms relevant for further discussion. In particular, we can ignore the one-body loss operators as they are linearly independent from the Hamiltonian and hence will not contribute to the bound. Similarly as in the linear case, we assume we deal with -atom states. Hence, we will utilize the fact that implying the following relation: , which allows us to introduce again a free parameter into coefficients of : , , , . Using Eq. (5) we arrive at the following bound:

 FQ≤minξT4∥∥∥2γ−112(ξ+1)2N(N−1)\openone+a†21a21(γ−111(ξ−1)2−2γ−112(ξ+1)2)+a†22a22(γ−122(ξ−1)2−2γ−112(ξ+1)2)∥∥∥. (40)

Since all operators under the operator norm commute and have common eigenbasis , we can write the above bound explicitly replacing with and with . Calculating the operator norm amounts now to maximization over . In case the above problem has an explicit solution which in the limit of large reads:

 FQ≤N2Tγ12⎧⎨⎩2(1+√λ)2,λ≥111+λ,λ<1. (41)

where . To the best of our knowledge this is the first example of a fundamental bound in the non-linear metrology model taking into account many-body decoherence effects. Details of the above derivation as well as the discussion of the case is presented in Appendix D.

Vi Quantum metrology with general many-body interactions

Let us now investigate what can be said in general concerning the fundamental bounds in non-linear metrology models with many-body interactions, without specifying the actual form of the dynamics but just its non-linear character and trying to identify the resulting characteristic scaling of QFI.

Consider a system of atoms, where the hamiltonian part is a result of interactions while the noise part is of type. To be more specific we consider the dynamics of the form:

 dρdt=−iω⎡⎣∑ν∈ΥkHν,ρ⎤⎦+γ∑μ∈Υl,jLμ,jρL†μ,j−12ρL†μ,jLμ,j−12L†μ,jLμ,jρ, (42)

where represents all element combinations of the element set, and the operator index denotes particles that a given operator acts on. We have also introduced a positive coefficient in order to be able to discuss effects of rescaling the noise strength. In what follows we will assume that . This is a general scenario considered in the field of non-linear metrology Boixo et al. (2007); Luis (2007); Boixo et al. (2008); Choi and Sundaram (2008); Napolitano and Mitchell (2010); Gross et al. (2010); Hall and Wiseman (2012); Joo et al. (2012); Sewell et al. (2014), but very often without the noise part. Including the noise part in this form in the analysis is extremely challenging and has only been analyzed for a very specific noise models and only in case of input GHZ states Beau and del Campo (2016).

Clearly, plugging all operators directly into formulas for and would make the problem intractable in case of large . We show that it is possible to apply our tools effectively to the dynamics of , , atoms and from this analysis infer the final scaling for the whole system of atoms. Recall that while investigating the fundamental limits to adaptive schemes we always consider the limit , and the only relevant order of the dynamics we need to take into account is the linear one. Hence, we may replace the original dynamics as represented by (42) with a scheme where each -tuple of particles experiences the dynamics sequentially. By a Trotter expansion argument this will introduce only a difference due to a potential lack of commutation of the operators acting on different subsets of particles, see Fig. 2.

This represents the dynamics in terms o elementary operations acting on particles only, which we refer to as subchannels. In the above scheme the -box denotes a free unitary evolution of particles for a time , while the -box represent the noisy part of the dynamics lasting also for a time . In order to keep the equivalence to the original problem we need to rescale the noise coefficient , where . The rescaling is necessary since the number of noisy gates applied is enhanced by a factor of compared with the original dynamics. Similarly, we need to modify the Hamiltonian evolution by rescaling —in the example of Fig. (2) this is not necessary since .

Let us assume that it is possible to find that makes , where should be understood as operator corresponding to the elementary dynamics . This again corresponds to the situation that the Hamiltonian belongs to which is constructed from noise operators entering according to Eq. (2). We can now apply the bound (11) treating as the fundamental building block for the adaptive strategy and since there are now such elementary blocks we arrive at: .

Let us inspect Eqs. (14,16) in order to understand the impact of the rescaling factors , on the value of the above bound. Rescaling of introduces an additional factor to all operators. Taking additionally into account that the Hamiltonian is rescaled by factor, then according to in order to satisfy constraint, needs to be rescaled by while by factors.

Together with (16) this implies that is rescaled by . Therefore we finally arrive at

 FQ≤4T∥α(1)εω,γ∥(Nn)χlχ2k∝TN2k−l. (43)

Note that the obtained scaling agrees with what we have obtained in the models discussed in Sec. V. In that case the situation where could be made zero corresponded to either or , and indeed we obtained respectively , and scalings. This shows that the general approach based on splitting the complex multiparticle dynamics into small subchannels involving only the number of particles required to model a given degree of nonlinearity is sufficient to obtain a proper scaling of the precision bounds. This general approach can also be utilized to obtain quantitative bounds which will be the subject of a separate publication Czajkowski and Demkowicz-Dobrzanski ().

The scaling of QFI or equivalently scaling of parameter estimation precision has been also observed in Beau and del Campo (2016), for protocols based on utilizing GHZ class of states and models where all and operators commute. Our approach proves the scaling for both the most general class of states and adaptive strategies. Once our bound can be derived (i.e. can be set to zero), we can claim that in models considered in Beau and del Campo (2016) indeed the GHZ as well as product states provide the optimal scaling. In the more general approach, however, where arbitrary states and adaptive strategies are allowed, this will not necessarily be the case.

To show both the power and simplicity of our approach, let us therefore consider the class of dynamics, which we will refer to as non-linear metrology with multi-particle dephasing, where and commute as considered in Beau and del Campo (2016), and try to apply our methods to derive general precision bounds in this case. Let us take and . Let us inspect the structure of the subspace , see Eq. (2), and invoke the representation of the dynamics in terms of subchannels as depicted in Fig. 2 and ask for what , we can satisfy the condition. The obvious case is when we simply consider subchannels in which case the operator is proportional to the operator . Similarly we can show that if , since if we take we can obtain from products of acting on two separate subsets of particles. More generally, provided is even and , we can consider particle subchannels and obtain by multiplying acting on two sets of particles where of them overlap. This way the product of two on overlapping particles produce the identity and we can obtain a fold tensor product of acting on the remaining particles.

In particular for non-linear metrology and linear dephasing we automatically get bound while for two-body Hamiltonian and nonlinear dephasing we get , proving fundamental character of the scaling obtained in Beau and del Campo (2016). However, if we take and , then according to our approach —as we cannot obtain a single from products of two or four that appear in the definition of . Note, that indeed in this case the subspace spanned by is actually immune to decoherence as all operators involve product of two acting on different sites and therefore yield a trivial factor when acting on states from this subspace. Still the single body Hamiltonian acts nontrivially, and we get showing the possibility of better than linear scaling both in and in and demonstrating that indeed in this case the bound is invalid.

Let us also comment here on the issue of optimization of quantum metrological protocols under the fixed definite particle number or just the mean number of particles fixed. In linear metrology QFI scales quadratically in the decoherence-free case and hence maximization of QFI under fixed mean particle number may lead to surprising conclusions of possibility of beating Heisenberg scaling or even reaching arbitrary high values, see e.g. Anisimov et al. (2010). This should be understood as deficiency of the QFI figure of merit which in general provides only a lower bound on achievable uncertainty via the Cramér-Rao bound whereas saturability of the bound requires more detailed arguments Tsang (2012); Berry et al. (2012); Giovannetti and Maccone (2012); Hall and Wiseman (2012). If decoherence makes QFI scale linearly with the number of particles, though, this issue becomes non-existent and one can replace fixed particle numbers with mean particle numbers in the formulas for the bounds Demkowicz-Dobrzański et al. (2013); Jarzyna and Demkowicz-Dobrzański (2015); Demkowicz-Dobrzanski et al. (2015). In the non-linear models discussed above, even in presence of decoherence, the bounds will in general scale super-linearly with particle numbers and hence again the task of maximizing QFI under the fixed mean particle number will be ill-posed and as a result meaningful discussion of such problems would require going beyond the QFI paradigm.

Vii Conclusions and discussion

Throughout this work, we have assumed Markovian semi-group dynamics. Validity of this approximation requires in particular coarse graining of the evolution of the system on time scales larger than the environmental characteristic relaxation time scales. Therefore the ,,continuous limit” , we have adopted while deriving our bounds, should be understood as the limit of very short times, but still in the regime where the Markovian approximation holds. Consequently, fundamental character of the bounds we derive, hinges upon the assumption that one cannot operate on time scales shorter than the ones characteristic for the Markovian approximation. If this assumption is dropped one may expect different time scalings of QFI, see e.g. Chin et al. (2012); Smirne et al. (2016); Sekatski et al. (2016b). Still, combining the framework of the most general adaptive protocols with non-Markovian or even Markovian non-semigroup dynamics, is a non-trivial task. In particular, focusing solely on the reduced system dynamics may not be in general sufficient to describe the effects of control operations acting on the system, and more information on system+environment interaction may be required. Moreover, in such scenarios, one cannot a priori justify the framework where the time used up by control operations is not relevant, as well as argue that a model with infinitely many control operations is the optimal one. We feel that providing a general framework to study the potential of adaptive quantum metrological protocols in presence of non-Markovian noise goes beyond direct generalization of the methods of this paper and leave this task for future research.

Another way of generalizing results of this work is to assume the environment can be partially monitored, which may open up completely new possibilities of fighting decoherence. It is also interesting to study multiparameter estimation scenarions were adaptiveness seems to play some role already at the decoherence-free level Yuan (2016).

To conclude. We have provided a simple algebraic criterion, Eq. (2), determining at most linear time scaling of QFI in estimation of a unitary parameter under general Markovian dynamics. Its high utility stems from the fact that it deals directly with the Hamiltonian and noise operators appearing in the Master equation and does not require integration of the dynamics. We have shown how the bounds can be derived in atomic interferometric models involving many-body interactions. In the qubit case it has been shown that the condition (2) is sufficient and necessary for fundamental linear time scaling of QFI. In case of arbitrary dynamics in arbitrary dimension, however, it has only been proven to be a sufficient condition. Whether it is also a necessary one remains an open question.

Note added. After completion of this work, a paper Zhou et al. (2017) appeared, where apart from independent derivation of the results presented in Sec. III of our paper, an explicit construction of quantum error-correction protocol has been provided yielding scaling of the QFI whenever condition (2) is not satisfied. As a result, this paper answered the open problem stated above on whether the condition (2) is indeed the if and only if condition for impossibility of preserving the scaling of QFI via the most general quantum protocols in the presence of noise for systems of arbitrary dimension. Note also that our conditions for an effective error correction protocol (25-27) are equivalent to the ones provided in Zhou et al. (2017). In particular by expressing projector from Zhou et al. (2017) using , from Sec.IV of our paper as , we see that their condition (15) is equivalent to ours (26-27) while their condition (17) is equivalent to ours (25)—on the one hand one can always find two states in this subspace that are coupled by the Hamiltonian, on the other if two states from the subspace are coupled the action of the Hamiltonian is nontrivial. Finally, let us briefly describe how the universal error correction code of Zhou et al. (2017) appears in our language. The key idea is to use the spectral decomposition of the part of the Hamiltonian perpendicular to the noise space , where and are orthogonal, positive semi-definite and since is traceless by construction. Next, on the probe plus ancilla system define the corresponding purifications states and satisfying

 P=h⊥TrA|p⟩⟨p|andQ=h⊥TrA|q⟩⟨q|, (44)

where since is non-zero. Then the virtual qubit states of Sec. IV can be defined as and .

Acknowledgements

We thank Wolfgang Dür, Janek Kołodyński, Animesh Datta, Philipp Treutlein and Krzysztof Pawłowski for fruitful discussions. This work was supported by the Polish Ministry of Science and Higher Education Iuventus Plus program for years 2015-2017 No. 0088/IP3/2015/73, National Science Center (Poland) grant No. 2016/22/E/ST2/00559 and Swiss National Science Foundation grant P300P2_167749.

Appendix A Calculating the bound via semi-definite programming

Given the condition is satisfied, we know that the scaling of QFI is bound to be linear in . In order to obtain the tightest bound of the form (12) we need to minimize the operator norm of keeping the constraint . First we construct a matrix

 A=⎡⎣√λ\openoneh(12)†\openone+L†h(0)h(12)\openone+h(0)L√λ\openone⊗J⎤⎦. (45)

Minimizing the operator norm is now equivalent to minimizing subject to with the additional constraint coming from the equation . The QFI bound can therefore be written as

 FQ≤4Tmin{h(1)00,h(12),h(0)}λ,subject to: A≥0, β(1)=0. (46)

The problem of determining the bounds is now fully specified as a semi-definite program using only the operators appearing in the Master equation.

Appendix B Tightening the bound using physical constraints

The bound in the form (11), which involves operator norms, may in some cases be tightened by taking into account additional physical constraints present in the considered problem.

b.1 Superselection rules

It is often the case in practice that the experimentalist does not have access to operations that create coherence between some sectors of a the total Hilbert space. For example, in no photonic experiment one is able to create, manipulate or detect coherence between sectors of different total number of particles in all modes (including local oscillators). This is a consequence of symmetries, like the time translation symmetry enforcing energy conservation in the example with particles. It is usually the case that the evolution, Eq. (1), does not create such coherences either, unless the evolution is “active” and describes the interaction of the system with a source. This can be formalized by introducing a map

 P(⋅)=∑kPk⋅Pk (47)

with , that projects the state of the systems onto sectors that are eigenspaces of the conserved quantities labeled by . Control operations and evolution are incoherent they commute with map , implying for instance that is block diagonal .

Under such circumstances QFI of a state might be overestimating the extractable classical Fisher information of the estimated parameter, as the information might be, for example, encoded in the coherence between different sectors. To get a tight expression one then has to explicitly account for the conservation law, by looking on QFI of the incoherent state . A similar treatment can be done in our case. Since the projection map commutes with both the controls and the evolution and has to be applied on the final state, we can explicitly account for the conservation law by considering the process where the free evolution induced by the master equation (1) is interjected with after each infinitesimal time step

 (48)

This composed channel