Hilbert-Schmidt speed as an efficient figure of merit for quantum estimation of phase encoded into the initial state of open n-qubit systems.

Hilbert-Schmidt speed (HSS) is a special type of quantum statistical speed which is easily computable, since it does not require diagonalization of the system state. We find that, when both HSS and quantum Fisher information (QFI) are calculated with respect to the phase parameter encoded into the initial state of an n-qubit register, the zeros of the HSS dynamics are actually equal to those of the QFI dynamics. Moreover, the signs of the time-derivatives of both HSS and QFI exactly coincide. These findings, obtained via a thorough investigation of several paradigmatic open quantum systems, show that HSS and QFI exhibit the same qualitative time evolution. Therefore, HSS reveals itself as a powerful figure of merit for enhancing quantum phase estimation in an open quantum system made of n qubits. Our results also provide strong evidence for both contractivity of the HSS under memoryless dynamics and its sensitivity to system-environment information backflows to detect the non-Markovianity in high-dimensional systems, as suggested in previous studies.

Introduction

The process of quantum parameter estimation typically consists in preparing the probe system in an initial state, letting it interact with the target system to encode the information about the unknown parameter, then measuring the probe to extract the information and finally estimate the parameter. It should be noted that in this process, the system also may be affected by different noises. Provided that the physical mechanism governing the system dynamics is known, one may deduce a value of the parameter, with a given precision, by comparison between the input and the output states of the probe[1,2].

Phase estimation is at the heart of quantum metrology[2-15] due to the fact that, in many technological areas, the parameter estimation problem can be reduced to determining a phase shift induced in the quantum state of the probe. Measuring an unknown phase has many significant applications in different scenarios, some of them being gravitational wave observation[16-19], test-mass position measurements[20], frequency standards[21], gravimeters[22,23] and detection of weak signals or defects resulting in the design of very sensitive sensors[24]. In these scenarios an interferometric scheme is typically used to implement the quantum phase estimation, as happens for optical interferometry in gravitational wave detectors[25], atomic interferometry in Ramsey spectroscopy and physical law testing[13], optical imaging or laser gyroscopes[2] to name a few. All of these applications usually aim at optimal estimation of a relative phase gathered by the signal along one arm of the interferometer[26].

According to the quantum Cramér-Rao theorem, the precision of quantum phase estimation is bounded by the inverse of the quantum Fisher information (QFI)[4,27-29], which thus represents a central quantity in quantum metrology. In fact, evaluation of the QFI provides the ultimate quantum limits to measurement precision and consequently a general benchmark to assess quantum metrological protocols. The QFI is also a measure of quantum statistical speed that quantifies the sensitivity of an initial state with respect to changes of the parameter which should be estimated. In general, each measure of statistical distance naturally leads to a statistical speed for parametric evolutions of classical probability distributions or quantum states. This statistical speed can be obtained by the change in distance originated from a small variation of this parameter (i.e., the derivative of the distance). The quantum statistical speed is then obtained by maximizing over the classical statistical speed over all quantum measurements[30].

Inspired by the fact that the QFI is a measure of quantum statistical speed, obtained from the Hellinger distance[31], here we investigate the application of the Hilbert–Schmidt speed (HSS), another quantifier of quantum statistical speed, in the process of quantum phase estimation. Because calculating the QFI for high-dimensional quantum systems is typically complicated[32], it is useful to inquire the efficiency of the HSS, which is an easily computable quantity having the advantage of avoiding diagonalization of the evolved density matrix, in the quantum estimation theory. In this paper, we show that the HSS can be indeed exploited as a powerful and convenient figure of merit in quantum metrology for n-qubit systems. This result gains particular attention considering the fact that many quantum information protocols are designed by n-qubit registers.

Results

Preliminaries: quantum Fisher information (QFI)

We start our analysis by recalling the general formulation leading to defining a kind of quantum statistical speed by which the quantum Fisher information (QFI) can be characterized.

First, we consider the (classical) Hellinger distance[31]where and represent the probability distributions, i.e., probabilities of occurrence of different possible outcomes for two different experiments. Here the random variable x is assumed to take only discrete values. Formally, in order to achieve the statistical speed from a given statistical distance, one should quantify the distance between infinitesimally close distributions taken from a one-parameter family with parameter . In other words, the statistical speed versus a special parameter is defined as the change in the corresponding distance originated from a small variation of the parameter. Following this prescription and performing a Taylor expansion at for small values of , we find that the classical statistical speed associated with the (classical) Hellinger distance is given bywheredenotes the Fisher information[4,28].

Extending these classical notions to the quantum case and considering a given pair of quantum states and , one may write and representing the measurement probabilities associated with the positive-operator-valued measure (POVM) defined by the set which satisfies . The associated quantum distance can be obtained by maximizing the (classical) Hellinger distance over all possible choices of POVMs[33], that iscalled the Bures distance[28] in which denotes the state fidelity[34]. Now, the quantum statistical speed[28] can be defined asin which the QFI appears and is given by[1,5,6,28]where and denote, respectively, the eigenvectors and eigenvalues of the density matrix . In fact, the QFI is obtained by maximizing the Fisher information over all possible POVMs[4], that iswhere and . The fundamental relationship between the QFI and its corresponding quantum bound is expressed by the quantum Cramér-Rao bound[28]which sets the precision limit for quantum estimation of unknown parameter : .

Preliminaries: Hilbert–Schmidt speed (HSS)

We now report the definition of the Hilbert–Schmidt speed (HSS). Introducing the distance measure[30]where as well as are probability distributions, and subsequently considering the classical statistical speedone can define a special kind of quantum statistical speed which is called HSS. Following the procedure discussed in the previous subsection for obtaining the corresponding quantum relations, one can determine the Hilbert–Schmidt distance[35]not requiring the diagonalization of the argument operator. The corresponding quantum statistical speed is just the HSS given by[30]which can be computed without diagonalizing . This is a significant advantage from a computational perspective.

Quantum estimation through HSS

Because both QFI and HSS are quantum statistical speeds associated, respectively, with the Bures and Hilbert–Schmidt distances, it is reasonable to investigate how they can be related to each other. Based on analytical and numerical calculation, we find that an important relationship exists between them. We summarize our main result in the following.

Main Result. Suppose that we are given a pure initial state of an n-qubit quantum register such aswhere is the normalization factor and denotes the computational basis of the collective Hilbert space of dimension d. Then, let this state be affected by an arbitrary quantum channel giving the output state . Under these conditions, we find that and , computed with respect to phase parameter encoded into the input state of Eq. (13), exhibit the same qualitative dynamics, that is: Notice that the relation (i) above immediately implies that maxima and minima times of QFI and HSS coincide. Therefore, by investigating the HSS dynamics we can faithfully detect the instants when the optimal phase estimation is achieved.

if , one has , where indicates the sign of the argument x.

.

The sanity check of these relationships is performed by presenting various physical examples in the following subsections. It should be noted that using the general hierarchy between the HSS and QFI discussed in the context of quantum statistical speeds[30], one can show that , hence leads to . However, the reverse (i.e., detecting the QFI zeros through the HSS zeros, which is the study of this paper) cannot be extracted from the above general inequality. Our subsequent analysis just confirms that if HSS vanishes then QFI also vanishes, leading to our main result. We also point out that, in the following, the special choice of the initial state of the system is chosen in order to simplify the analytical expressions of the QFI and the HSS. It is emphasized that our main result is independent of which one of the phase parameters appearing in Eq. (13) is estimated. The results along the manuscript are general and hold for any other choices of the initial state where a relative phase parameter appears.

One-qubit system

First we focus on the paradigmatic one-qubit system where the qubit interacts with a dissipative reservoir (cavity), as depicted in Fig. 1, through the Hamiltonian ()[36]where is the transition frequency of the qubit, represent the system raising and lowering operators, is the frequency of the k-th field mode of the reservoir, () denotes the k-mode annihilation (creation) operator, and with being the coupling constant with the k-th mode. The cavity is initially in the vacuum state and uncorrelated from the qubit. Consider then the qubit prepared in the initial statewhere the parameter is the one to be estimated, while just fixes the probability amplitudes in the quantum superposition.

Figure 1

Sketch of a single qubit of transition frequency embedded in a dissipative zero-temperature cavity, with spectral width (rate of photon losses). The qubit is initially prepared in a given superposition of its two states, with a relative phase parameter . This figure was created using Keynote, version 10.3.9 (7029.9.8), URL: https://apps.apple.com/it/app/keynote/id409183694?mt=12.

In the limit of a continuum of reservoir modes , where denotes the reservoir spectral density[36,37]. Restricting ourselves to the case of a single excitation in the atom-cavity system, we find that the cavity mode can be eliminated in favour of an effective spectral density of the following form[36]where parameter , connected to the reservoir correlation time by , represents the spectral width for the qubit-reservoir coupling. Moreover, the decay rate is related to the system (qubit) relaxation time scale , over which the state of the system changes, by means of .

From the Hamiltonian of Eq. (14), at zero temperature and in the strong-coupling regime with the above Lorentzian spectral density for the cavity modes, one finds that the dynamics of the qubit in the basis is described by the evolved reduced density matrix[36,38]where the coherence characteristic function P(t) iswith .

Inserting the time-dependent density matrix of Eq. (17) into Eqs. (6) and (12), we find that the QFI and HSS associated with an initial phase are, respectively, given byleading to the relationsAccordingly, we see that when the HSS vanishes, the QFI also equals zero. Moreover, at all times when , the signs of and coincide (HSS is a nonnegative quantity) and hence they exhibit the same qualitative dynamics. In particular, the times when the optimal estimation is achieved, i.e., (maximum points of the QFI), can be easily detected by looking at the HSS dynamics.

Two-qubit systems

The validity of our main result is here discussed for two-qubit systems in three different scenarios: local coupling to independent environments, coupling to a common environment, teleportation of entanglement between the two qubits.

Coupling to independent environments

We consider a composite quantum system which consists of two separated identical qubits independently interacting with their own zero-temperature dissipative reservoir, as illustrated in Fig. 2. The two local qubit-reservoir interactions are assumed to be the same, for the sake of simplicity and without losing any generality. The two cavities (reservoirs), each exhibiting a Lorentzian spectral density, are initially in the vacuum state and uncorrelated from the two-qubit system. Knowing the evolved density matrix of the single qubit discussed in previous subsection, one can easily obtain the density matrix evolution of the two independent qubits[38-40]. In the computational basis , we take the two qubits initially prepared in the statewith the phase parameter , giving rise to the evolved reduced density matrixwhere is the coherence characteristic function of Eq. (18).

Figure 2

Sketch of two separated identical qubits each one interacting with its own zero-temperature dissipative cavity (). The initial state of the two qubits is such that a relative phase parameter appears between two-qubit basis states. This figure was created using Keynote, version 10.3.9 (7029.9.8), URL: https://apps.apple.com/it/app/keynote/id409183694?mt=12.

Computing the QFI of Eq. (6) and the HSS of Eq. (12) with respect to the phase parameter , one promptly gets, respectively, the expressionsleading toThe expressions above show once again our main result. In fact, assuming that , we have that the two time derivatives and have the same sign. In addition, at times when , the QFI also vanishes and hence no information can be extracted from the system. Similar expressions would be obtained for different choices of the two-qubit initial state.

Coupling to a common environment

We study two separated qubits interacting with a common zero-temperature bosonic reservoir, as depicted in Fig. 3. In such a case, the interaction between the two qubits is not direct but mediated by the common quantum reservoir. The total Hamiltonian of the system plus the reservoir is written as , with[41]where and denote, respectively, inversion operator and transition frequency of the j-th qubit (), () represents the k-mode creation (annihilation) operator of quanta (photons) of the environment, and in which is the coupling constant with the k-th mode. Moreover, the interaction of the j-th qubit with the reservoir is quantified by the dimensionless constant depending on the value of the cavity field at the qubit position, which can be effectively controlled by means of dc Stark shifts tuning the atomic transition in and out of resonance. We investigate the case when the two atomic qubits interact resonantly with the reservoir described by a Lorentzian spectral density and they have the same transition frequency, i.e., .

Figure 3

Illustration of two separated qubits both interacting with a zero-temperature common dissipative cavity. The initial state of the qubits is such that a relative phase parameter appears between the two-qubit basis states. This figure was created using Keynote, version 10.3.9 (7029.9.8), URL: https://apps.apple.com/it/app/keynote/id409183694?mt=12.

It is useful to introduce a collective coupling constant , the relative strengths such that , and mutually orthogonal quantum statesThe cavity is initially in its vacuum state , uncorrelated from the two-qubit system. With these definitions, for an initial two-qubit state of the formone finds that the evolved reduced density matrix of the two-qubit system in the computational basis is given by[41]where, defining , one hasMoreover, introducing the dimensionless quantities and , where is the reservoir correlation time and denotes the vacuum Rabi frequency, one obtainsThe above equations allow us to calculate the relevant quantities for our study, that is QFI and HSS. The QFI can be computed analytically but its expression is too cumbersome to be presented here. On the other hand, we find that the HSS is given byThe dynamical comparison between the two figures of merit is then performed by plots. In particular, we show in Fig. 4 that both the QFI and HSS dynamics simultaneously exhibit an oscillatory behavior such that their maximum and minimum points exactly coincide. This plot qualitatively verifies our result that the HSS can detect exactly the times at which the best phase estimation occurs (maximum of the QFI). Therefore, also for this system, the HSS, similarly to QFI, can be used as a distinguishability metric on the space of quantum states which quantifies the maximum amount of information on an unknown phase parameter attainable by a given probe state.

Figure 4

Dynamics of quantum Fisher information (red solid line) and Hilbert–Schmidt speed (amplified by 1.4 times for comparison, blue dashed line), as a function of the dimensionless time for the two-qubits system coupled to a common reservoir, with and . Notice that approximately oscillates between 0.3 and 1 due to the special choice of the initial parameters.

Two-qubit entanglement teleportation

We now study the relationship between QFI and HSS within the scenario of noisy two-qubit teleportation, where a two-qubit entangled state is to be teleported from a place to another[42,43]. Let us first describe the noisy two-qubit system to this aim, which shall constitute the resource channel for the teleportation.

One of the most important noise models used in the low-temperature regime is described by a spin environment[44,45]. In particular, in order to achieve suitable operations in experiments performed to study macroscopic quantum coherence and decoherence, one requires temperatures close to absolute zero. Here we consider a two-qubit system interacting with an external environment composed of N spins. The general Hamiltonian is written as where the system, environment and interaction Hamiltonians are given, respectively, by[44,45]where and denote, respectively, the characteristic frequency of the i-th qubit and the coupling strength between the two spin qubits. Moreover, denotes the tunneling matrix element for the i-th environmental spin, while () represents the coupling between qubit 1 (qubit 2) and the spins of the environment.

Initializing the two-qubit system in an extended Werner-like state[40]where quantifies the degree of purity of the state, is the identity matrix andone finds that the evolved reduced density matrix in the two-qubit computational basis is given by[45]where the decoherence factor Q(t) is

Assuming that the two qubits are shared between Alice (location A) and Bob (location B), we use two copies of this system state as a resource channel for the teleportation of an unknown entangled state from A to B, according to a basic scheme[42] illustrated in Fig. 5. It is useful to introduce the Bell states ’s associated with the Pauli matrices ’s bywhere ( being the identity matrix), , , and . Now, one can generalize the standard teleportation protocol[42] and find that the output state of the two-qubit teleportation is[43]where , being the resource state for the teleportation which is equal to the reduced density matrix of Eq. (35) in our model. Accordingly, for the input state withwhere , we find that the output state of the teleportation is given bywith and . Using this expression for the output state, we find that the QFI and HSS associated to the phase parameter encoded into the input state sent into the teleportation channel are, respectively,From these expressions we straightforwardly obtainleading to our main result, i.e., the possibility of extracting the QFI dynamics through the HSS dynamics due to the fact that the two figures of merit have the same qualitative time behavior.

Figure 5

Sketch of the noisy two-qubit teleportation protocol. Two copies of the evolved noisy two-qubit state are employed as resource channel for the protocol, shared between the two locations A and B. The input state to be teleported, containing information about the phase parameter , is manipulated in location A. After the protocol, the state is teleported at the location B encoded in the two qubits coming from the resource state[42]. This figure was created using Keynote, version 10.3.9 (7029.9.8), URL: https://apps.apple.com/it/app/keynote/id409183694?mt=12.

General case: n-qubit system ()

Now the validity of the main result is investigated for high-dimensional systems in two scenarios: local coupling to independent environments and coupling to a common environment.

Here we analyze the case of an open quantum system made of n noninteracting qubits locally subject to their own reservoir, a typical configuration for quantum networking with quantum registers. To this aim, one has to choose the type of qubit and the local qubit-environment interaction.

We specifically consider the dynamics of a topological qubit realized by two Majorana modes which are generated at the endpoints of some nanowire with strong spin-orbit interaction, placed on top of an s-wave superconductor and driven by an external magnetic field B along the wire axis direction[46,47]. We also assume that each Majorana mode is coupled to the metallic nanowire via a tunnel junction in the way that the tunneling strength is controllable by an external gate voltage. The total Hamiltonian is written aswhere denotes the Hamiltonian of the topological qubit, V represents the system-environment interaction Hamiltonian, while is the environment Hamiltonian whose elementary constituents can be considered as electrons or holes. The decoherence which affects the topological qubit is modeled as a fermionic Ohmic-like environment described by spectral density with . The Ohmic, super-Ohmic and sub-Ohmic environments are characterized by , and , respectively. Since these Majorana modes used as the topological qubit are zero-energy modes, we have . Moreover, the interaction Hamiltonian V, constructed by the electron creation (annihilation) operators with Majorana modes and , satisfies the propertiesBefore turning on the interaction V, the two Majorana modes construct a topological (non-local) qubit with states and related to each other bywhere , are chosen to be represented by , , , in which ’s are the usual Pauli matrices.

Let us now assume, as usual, that the state of the total qubit-environment system is uncorrelated initially: , where and are the initial density matrices of the topological qubit and of its environment, respectively. Supposing that the initial state of the Majorana qubit is written asone finds that the reduced density matrix of the topological qubit at time t can be obtained by a dynamical map such that (for details, see Ref.[46])wherewith indicating the high-frequency cutoff for the linear spectrum of the edge state, representing the Gamma function, and being the generalized hypergeometric function . From the eigenvalues and eigenvectors of the Choi matrix[48] of the map , we obtain the corresponding operator-sum representation with Kraus operators given byIn the following, we take a system formed by n noninteracting topological qubits (each one defined by two Majorana modes as above) such that each qubit locally interacts with the environment described above by the interaction (), as displayed in Fig. 6. Note that the effects of the environment on each of the qubits can be canceled by setting the corresponding external magnetic field to zero. With this in mind, we focus on the scenarios in which m of the qubits are affected by the noise, while the remaining ones () are noiseless. Since the environments are independent in our model, the Kraus operators are just tensor products of Kraus operators acting on each of the qubits, noting that the Kraus operators of the noiseless qubits are set to identity operator[38,49].

Figure 6

Illustration of a system of n noninteracting topological qubits (), each one built by two Majorana modes and (only , are explicitly indicated for simplicity), locally interacting by to the environment (constituted by electrons and holes). The n qubits are initially prepared in a state containing relative phase parameters between the computational basis states. This figure was created using Keynote, version 10.3.9 (7029.9.8), URL: https://apps.apple.com/it/app/keynote/id409183694?mt=12.

We start by considering a three-qubit system with , initially prepared in the W-like state[50]where two phase parameters , appear. We find that the QFI and HSS associated with phase parameter is obtained asIt is easily found that

As clear from Eq. (51), the HSS is always nonzero and we conclude that the QFI dynamics can be completely determined by analyzing the HSS dynamics: our main result is once again proved. Now we compute the measures associated with the phase parameter , leading to the expressionsfrom which we obtainso that our main result is manifestly confirmed also in this case.

In Fig. 7, we numerically illustrate the same qualitative time behavior of the HSS and QFI for two different initial states: (a) a four-qubit system and (b) a six-qubit system (with ), initially prepared in the following W-like statesandrespectively.

Figure 7

Dynamics of quantum Fisher information F(t) (red solid line) and Hilbert–Schmidt speed HSS(t) (amplified by 1.3 times for comparison, blue dashed line) which are computed with the same phase parameter as a function of time t for the (a) four-qubit system prepared in initial state (55) and (b) six-qubit system prepared in initial state (56) described in Fig. 6 when one of the qubits is affected by the noise.

As another example, we choose n and m arbitrary with the n-qubit register prepared in a Greenberger-Horne-Zeilinger (GHZ)-like state[50,51] written aswhere the relevant phase parameter appears. Calculating the evolved state of the system, we find that the corresponding QFI and HSS are given, respectively, byHence, we can writeexplicitly leading to our main result also in this general case of n qubits, as already confirmed by the other previous examples with a fixed number of qubits.

Coupling to a common environment: depolarizing channel

Now we investigate the phase estimation problem for a natural parametrization of arbitrary pure states affected by the (isotropic) depolarizing channel[49] given by the mapin which and d denote the reliability of the channel and dimension of the Hilbert space, respectively. This form of collective channel, mixing the initial n-qubit pure state with white noise and illustrated in Fig. 8, plays an important role in various quantum information tasks, such as quantum cloning machines[52], NMR quantum computing[53], entanglement optimization[54], quantum error correction[55], and quantum repeaters[56]. We focus on the following parametrization of input pure statesAfter lengthy calculation one can find that the general expression of QFI associated with any parameter is given by[57]Moreover, we obtain that the HSS corresponding to the parameter can be represented asleading to the following relationship between QFI and HSS:Because , our main result can be immediately confirmed. Combining this result with the fact that (see Eq. (63)), we see that when the reliability parameter of the depolarizing channel oscillates with time, an indication of non-Markovianity[58], the oscillations are in turn detected easily by the HSS.

Figure 8

Illustration of an n-qubit register collectively subject to white noise under the action of a common depolarizing channel. This channel turns an initial pure state into the mixture . This figure was created using Keynote, version 10.3.9 (7029.9.8), URL: https://apps.apple.com/it/app/keynote/id409183694?mt=12.

Coupling to a common environment: dephasing channel

Here, the validity of our main result is checked for the dynamics of an n-qubit system with total angular momentum in the presence of collisional dephasing, as depicted in Fig. 9. Using a unitary evolution, we can prepare the probes in the following GHZ-like initial statewhere for denotes the collective basis constructed from the eigenvectors of , representing the z component of the total angular momentum for all qubits, i.e., . The effect of the collisional dephasing channel on input state (65), can be modeled as the interaction between the qubits and the common thermal reservoir[59,60]. In the interaction picture, the master equation of the system can be expressed as[61]where represents the dephasing rate which may be time-dependent, and denotes the reduced density operator of the system. In the case of single-qubit systems, Eq. (66) reduces tocorresponding to a single-qubit dephasing channel. Because the QFI and the HSS remain invariant under the unitary evolution being independent of the phase parameter, working in the interaction picture does not affect the result of the calculation. Solving Eq. (82), one can find that the time evolution of the density matrix elements is expressed aswhere .

Figure 9

Illustration of an n-qubit register subject to the action of a common dephasing channel. This channel turns an initial pure GHZ state into a mixed state . This figure was created using Keynote, version 10.3.9 (7029.9.8), URL: https://apps.apple.com/it/app/keynote/id409183694?mt=12.

We are interested to compare the dynamics of the QFI and HSS when computed with respect to initial phase imprinted on the input state (65). The corresponding QFI is given by[60]In addition, we find that the HSS has the following simple exponential form:resulting in the relationby which our main result can be verified straightforwardly.

Discussion

Quantum information processing based on n-qubit registers provides a paradigmatic playground for fundamental research on composite quantum systems with potential technological advances[50,62-65]. For this reason, finding useful tools for the convenient characterization of quantum features within quantum registers is particularly desired.

In this paper, we have constructed a strong relationship for n-qubit systems between the Hilbert–Schmidt speed (HSS), which is a special case of quantum statistical speed, and the quantum Fisher information (QFI), a key concept in parameter estimation theory and in quantum metrology in general. The idea underlying this relationship stems from the fact that the QFI, quantifying the sensitivity of an initial state with respect to changes of the parameter of a dynamical evolution, is itself a quantum statistical speed extracted from the Bures distance. In fact, our detailed analysis, carried out for several paradigmatic open quantum systems ranging from one qubit to n qubits under different environmental conditions, shows that these two quantities exhibit the same qualitative dynamical behaviors. Therefore, in contrast to the typical computational complication of the QFI, especially for multipartite systems, our findings provide evidence that the HSS can be instead adopted as an efficient tool in the process of quantum phase estimation implemented by n-qubit systems, because of its straightforward determination.

We point out that the present results are an advance of previously reported analysis[66], where the HSS has been shown to be useful for metrological resolution of phase shift induced by a unitary transformation without noisy dynamics. Here, instead, we provide evidence that the HSS is a powerful figure of merit for the enhanced estimation of a phase encoded in the initial state of an open system made of n qubits subject to noisy channels. As an interesting outlook, these results suggest to extend the analysis to problems in noisy quantum phase estimation, which encompass cases when both parameter encoding and noisy dynamics happen simultaneously (as occurs, for instance, in frequency estimation processes under pure dephasing[67]).

On the one hand, it is known that the QFI monotonically decreases under Markovian dynamics, since it cannot increase under completely positive maps[68-70], so that it can be used as a witness of non-Markovianity. Originally, introducing a flow of QFI as , it has been suggested[71] that if for some t, then the time evolution is non-Markovian. Nevertheless, the efficiency of the QFI flow to detect the non-Markovianity in various scenarios has not been yet analytically compared with the other faithful witnesses of non-Markovianity (see Ref.[72] for a numerical study on the interplay between the dynamics of the QFI and the trace distance an important indicator of non-Markovianity). On the other hand, the HSS flow has been recently proposed as a faithful witness of non-Markovianity in low dimensional systems[73,74]. As a useful consequence, our results also provide a sanity check of the QFI flow as a witness of non-Markovianity. Moreover, because the QFI is always contractive under Markovian dynamics, our results supply strong evidence for contractivity of the HSS under memoryless evolution of high-dimensional systems and pave the way to further studies on its applications in measuring the non-Markovianity in open quantum systems made of qudits.

Finally, it should be noted that much attention has been recently devoted to exploring the use of continuous variable (CV) systems in quantum information processing[75-77]. It originates from the fact that the continuous-spectrum quantum variables can be easier to manipulate than quantum bits for performing various quantum information processes. Therefore, generalizing our results to the applicability of HSS in the CV context remains of interest and will be addressed elsewhere.