Quantum Fisher information

The quantum Fisher information is a central quantity in quantum metrology and is the quantum analogue of the classical Fisher information. The quantum Fisher information $F_{\rm {Q}}[\varrho ,A]$ of a state $\varrho$ with respect to the observable $A$ is defined as

$F_{\rm {Q}}[\varrho ,A]=2\sum _{k,l}{\frac {(\lambda _{k}-\lambda _{l})^{2}}{(\lambda _{k}+\lambda _{l})}}\vert \langle k\vert A\vert l\rangle \vert ^{2},$

where $\lambda _{k}$ and $\vert k\rangle$ are the eigenvalues and eigenvectors of the density matrix $\varrho ,$ respectively, and the summation goes over all $k$ and $l$ such that $\lambda _{k}+\lambda _{l}>0$.

When the observable generates a unitary transformation of the system with a parameter $\theta$ from initial state $\varrho _{0}$,

$\varrho (\theta )=\exp(-iA\theta )\varrho _{0}\exp(+iA\theta ),$

the quantum Fisher information constrains the achievable precision in statistical estimation of the parameter $\theta$ via the quantum Cramér–Rao bound as

$(\Delta \theta )^{2}\geq {\frac {1}{mF_{\rm {Q}}[\varrho ,A]}},$

where $m$ is the number of independent repetitions.

It is often desirable to estimate the magnitude of an unknown parameter $\alpha$ that controls the strength of a system's Hamiltonian $H=\alpha A$ with respect to a known observable $A$ during a known dynamical time $t$. In this case, defining $\theta =\alpha t$, so that $\theta A=tH$, means estimates of $\theta$ can be directly translated into estimates of $\alpha$.

Connection with Fisher information

Classical Fisher information of measuring observable $B$ on density matrix $\varrho (\theta )$ is defined as $F[B,\theta ]=\sum _{b}{\frac {1}{p(b|\theta )}}\left({\frac {\partial p(b|\theta )}{\partial \theta }}\right)^{2}$, where $p(b|\theta )=\langle b\vert \varrho (\theta )\vert b\rangle$ is the probability of obtaining outcome $b$ when measuring observable $B$ on the transformed density matrix $\varrho (\theta )$. $b$ is the eigenvalue corresponding to eigenvector $\vert b\rangle$ of observable $B$.

Quantum Fisher information is the supremum of the classical Fisher information over all such observables,

$F_{\rm {Q}}[\varrho ,A]=\sup _{B}F[B,\theta ].$

Relation to the symmetric logarithmic derivative

The quantum Fisher information equals the expectation value of $L_{\varrho }^{2}$, where $L_{\varrho }$ is the symmetric logarithmic derivative

Equivalent expressions

For a unitary encoding operation $\varrho (\theta )=\exp(-iA\theta )\varrho _{0}\exp(+iA\theta ),$, the quantum Fisher information can be computed as an integral,

$F_{\rm {Q}}[\varrho ,A]=-2\int _{0}^{\infty }{\text{tr}}\left(\exp(-\rho _{0}t)[\varrho _{0},A]\exp(-\rho _{0}t)[\varrho _{0},A]\right)\ dt,$

where $[\ ,\ ]$ on the right hand side denotes commutator. It can be also expressed in terms of Kronecker product and vectorization,

$F_{\rm {Q}}[\varrho ,A]=2\,{\text{vec}}([\varrho _{0},A])^{\dagger }{\big (}\rho _{0}^{*}\otimes {\rm {Identity}}+{\rm {Identity}}\otimes \rho _{0}{\big )}^{-1}{\text{vec}}([\varrho _{0},A]),$

where $^{*}$ denotes complex conjugate, and $^{\dagger }$ denotes conjugate transpose. This formula holds for invertible density matrices. For non-invertible density matrices, the inverse above is substituted by the Moore-Penrose pseudoinverse. Alternatively, one can compute the quantum Fisher information for invertible state $\rho _{\nu }=(1-\nu )\rho _{0}+\nu \pi$, where $\pi$ is any full-rank density matrix, and then perform the limit $\nu \rightarrow 0^{+}$ to obtain the quantum Fisher information for $\rho _{0}$. Density matrix $\pi$ can be, for example, ${\rm {Identity}}/\dim {\mathcal {H}}$ in a finite-dimensional system, or a thermal state in infinite dimensional systems.

Generalization and relations to Bures metric and quantum fidelity

For any differentiable parametrization of the density matrix $\varrho ({\boldsymbol {\theta }})$ by a vector of parameters ${\boldsymbol {\theta }}=(\theta _{1},\dots ,\theta _{n})$, the quantum Fisher information matrix is defined as

$F_{\rm {Q}}^{ij}[\varrho ({\boldsymbol {\theta }})]=2\sum _{k,l}{\frac {\operatorname {Re} (\langle k\vert \partial _{i}\varrho \vert l\rangle \langle l\vert \partial _{j}\varrho \vert k\rangle )}{\lambda _{k}+\lambda _{l}}},$

where $\partial _{i}$ denotes partial derivative with respect to parameter $\theta _{i}$. The formula also holds without taking the real part $\operatorname {Re}$, because the imaginary part leads to an antisymmetric contribution that disappears under the sum. Note that all eigenvalues $\lambda _{k}$ and eigenvectors $\vert k\rangle$ of the density matrix potentially depend on the vector of parameters ${\boldsymbol {\theta }}$.

This definition is identical to four times the Bures metric, up to singular points where the rank of the density matrix changes (those are the points at which $\lambda _{k}+\lambda _{l}$ suddenly becomes zero.) Through this relation, it also connects with quantum fidelity $F(\varrho ,\sigma )=\left(\mathrm {tr} \left[{\sqrt {{\sqrt {\varrho }}\sigma {\sqrt {\varrho }}}}\right]\right)^{2}$ of two infinitesimally close states,

$F(\varrho _{\boldsymbol {\theta }},\varrho _{{\boldsymbol {\theta }}+d{\boldsymbol {\theta }}})=1-{\frac {1}{4}}\sum _{i,j}{\Big (}F_{\rm {Q}}^{ij}[\varrho ({\boldsymbol {\theta }})]+2\!\!\sum _{\lambda _{k}({\boldsymbol {\theta }})=0}\!\!\partial _{i}\partial _{j}\lambda _{k}{\Big )}d\theta _{i}d\theta _{j}+{\mathcal {O}}(d\theta ^{4}),$

where the inner sum goes over all $k$ at which eigenvalues $\lambda _{k}({\boldsymbol {\theta }})=0$. The extra term (which is however zero in most applications) can be avoided by taking a symmetric expansion of fidelity,

$F\left(\varrho _{{\boldsymbol {\theta }}-d{\boldsymbol {\theta }}/2},\varrho _{{\boldsymbol {\theta }}+d{\boldsymbol {\theta }}/2}\right)=1-{\frac {1}{4}}\sum _{i,j}F_{\rm {Q}}^{ij}[\varrho ({\boldsymbol {\theta }})]d\theta _{i}d\theta _{j}+{\mathcal {O}}(d\theta ^{4}).$

For $n=1$ and unitary encoding, the quantum Fisher information matrix reduces to the original definition.

Quantum Fisher information matrix is a part of a wider family of quantum statistical distances.

Relation to fidelity susceptibility

Assuming that $\vert \psi _{0}(\theta )\rangle$ is a ground state of a parameter-dependent non-degenerate Hamiltonian $H(\theta )$, four times the quantum Fisher information of this state is called fidelity susceptibility, and denoted

$\chi _{F}=4F_{Q}(\vert \psi _{0}(\theta )\rangle ).$

Fidelity susceptibility measures the sensitivity of the ground state to the parameter, and its divergence indicates a quantum phase transition. This is because of the aforementioned connection with fidelity: a diverging quantum Fisher information means that $\vert \psi _{0}(\theta )\rangle$ and $\vert \psi _{0}(\theta +d\theta )\rangle$ are orthogonal to each other, for any infinitesimal change in parameter $d\theta$, and thus are said to undergo a phase-transition at point $\theta$.

Convexity properties

The quantum Fisher information equals four times the variance for pure states

$F_{\rm {Q}}[\vert \Psi \rangle ,H]=4(\Delta H)_{\Psi }^{2}$.

For mixed states, when the probabilities are parameter independent, i.e., when $p(\theta )=p$, the quantum Fisher information is convex:

$F_{\rm {Q}}[p\varrho _{1}(\theta )+(1-p)\varrho _{2}(\theta ),H]\leq pF_{\rm {Q}}[\varrho _{1}(\theta ),H]+(1-p)F_{\rm {Q}}[\varrho _{2}(\theta ),H].$

The quantum Fisher information is the largest function that is convex and that equals four times the variance for pure states. That is, it equals four times the convex roof of the variance

$F_{\rm {Q}}[\varrho ,H]=4\inf _{\{p_{k},\vert \Psi _{k}\rangle \}}\sum _{k}p_{k}(\Delta H)_{\Psi _{k}}^{2},$

where the infimum is over all decompositions of the density matrix

$\varrho =\sum _{k}p_{k}\vert \Psi _{k}\rangle \langle \Psi _{k}\vert .$

Note that $\vert \Psi _{k}\rangle$ are not necessarily orthogonal to each other. The above optimization can be rewritten as an optimization over the two-copy space as

$F_{Q}[\varrho ,H]=\min _{\varrho _{12}}2{\rm {Tr}}[(H\otimes {\rm {Identity}}-{\rm {Identity}}\otimes H)^{2}\varrho _{12}],$

such that $\varrho _{12}$ is a symmetric separable state and

${\rm {Tr}}_{1}(\varrho _{12})={\rm {Tr}}_{2}(\varrho _{12})=\varrho .$

Later the above statement has been proved even for the case of a minimization over general (not necessarily symmetric) separable states.

When the probabilities are $\theta$-dependent, an extended-convexity relation has been proved:

$F_{\rm {Q}}{\Big [}\sum _{i}p_{i}(\theta )\varrho _{i}(\theta ){\Big ]}\leq \sum _{i}p_{i}(\theta )F_{\rm {Q}}[\varrho _{i}(\theta )]+F_{\rm {C}}[\{p_{i}(\theta )\}],$

where $F_{\rm {C}}[\{p_{i}(\theta )\}]=\sum _{i}{\frac {\partial _{\theta }p_{i}(\theta )^{2}}{p_{i}(\theta )}}$ is the classical Fisher information associated to the probabilities contributing to the convex decomposition. The first term, in the right hand side of the above inequality, can be considered as the average quantum Fisher information of the density matrices in the convex decomposition.

Inequalities for composite systems

We need to understand the behavior of quantum Fisher information in composite system in order to study quantum metrology of many-particle systems. For product states,

$F_{\rm {Q}}[\varrho _{1}\otimes \varrho _{2},H_{1}\otimes {\rm {Identity}}+{\rm {Identity}}\otimes H_{2}]=F_{\rm {Q}}[\varrho _{1},H_{1}]+F_{\rm {Q}}[\varrho _{2},H_{2}]$

holds.

For the reduced state, we have

$F_{\rm {Q}}[\varrho _{12},H_{1}\otimes {\rm {Identity}}_{2}]\geq F_{\rm {Q}}[\varrho _{1},H_{1}],$

where $\varrho _{1}={\rm {Tr}}_{2}(\varrho _{12})$.

Relation to entanglement

There are strong links between quantum metrology and quantum information science. For a multiparticle system of $N$ spin-1/2 particles

$F_{\rm {Q}}[\varrho ,J_{z}]\leq N$

holds for separable states, where

$J_{z}=\sum _{n=1}^{N}j_{z}^{(n)},$

and $j_{z}^{(n)}$ is a single particle angular momentum component. The maximum for general quantum states is given by

$F_{\rm {Q}}[\varrho ,J_{z}]\leq N^{2}.$ Hence, quantum entanglement is needed to reach the maximum precision in quantum metrology.

Moreover, for quantum states with an entanglement depth $k$,

$F_{\rm {Q}}[\varrho ,J_{z}]\leq sk^{2}+r^{2}$

holds, where $s=\lfloor N/k\rfloor$ is the largest integer smaller than or equal to $N/k,$ and $r=N-sk$ is the remainder from dividing $N$ by $k$. Hence, a higher and higher levels of multipartite entanglement is needed to achieve a better and better accuracy in parameter estimation. It is possible to obtain a weaker but simpler bound

$F_{\rm {Q}}[\varrho ,J_{z}]\leq Nk.$

Hence, a lower bound on the entanglement depth is obtained as

${\frac {F_{\rm {Q}}[\varrho ,J_{z}]}{N}}\leq k.$

Relation to the Wigner–Yanase skew information

The Wigner–Yanase skew information is defined as

$I(\varrho ,H)={\rm {Tr}}(H^{2}\varrho )-{\rm {Tr}}(H{\sqrt {\varrho }}H{\sqrt {\varrho }}).$

It follows that $I(\varrho ,H)$ is convex in $\varrho .$

For the quantum Fisher information and the Wigner–Yanase skew information, the inequality

$F_{\rm {Q}}[\varrho ,H]\geq 4I(\varrho ,H)$

holds, where there is an equality for pure states.

Relation to the variance

For any decomposition of the density matrix given by $p_{k}$ and $\vert \Psi _{k}\rangle$ the relation

$(\Delta H)^{2}\geq \sum _{k}p_{k}(\Delta H)_{\Psi _{k}}^{2}\geq {\frac {1}{4}}F_{\rm {Q}}[\varrho ,H]$

holds, where both inequalities are tight. That is, there is a decomposition for which the second inequality is saturated, which is the same as stating that the quantum Fisher information is the convex roof of the variance over four, discussed above. There is also a decomposition for which the first inequality is saturated, which means that the variance is its own concave roof

$(\Delta H)^{2}=\sup _{\{p_{k},\vert \Psi _{k}\rangle \}}\sum _{k}p_{k}(\Delta H)_{\Psi _{k}}^{2}.$

Uncertainty relations with the quantum Fisher information and the variance

Knowing that the quantum Fisher information is the convex roof of the variance times four, we obtain the relation $$(\Delta A)^{2}F_{Q}[\varrho ,B]\geq \vert \langle i[A,B]\rangle \vert ^{2},$$ which is stronger than the Heisenberg uncertainty relation. For a particle of spin-$j,$ the following uncertainty relation holds $$(\Delta J_{x})^{2}+(\Delta J_{y})^{2}+(\Delta J_{z})^{2}\geq j,$$ where $J_{l}$ are angular momentum components. The relation can be strengthened as $$(\Delta J_{x})^{2}+(\Delta J_{y})^{2}+F_{Q}[\varrho ,J_{z}]/4\geq j.$$