This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Error trade-off relations for two-parameter unitary model with commuting generators

Shin Funada    Jun Suzuki Graduate School of Informatics and Engineering, The University of Electro-Communications, 1-5-1 Chofugaoka, Chofu-shi, Tokyo, 182-8585 Japan
Abstract

We investigate whether a trade-off relation between the diagonal elements of the mean square error matrix exists for the two-parameter unitary models with mutually commuting generators. We show that the error trade-off relation which exists in our models of a finite dimension system is a generic phenomenon in the sense that it occurs with a finite volume in the spate space. We analyze a qutrit system to show that there can be an error trade-off relation given by the SLD and RLD Cramer-Rao bounds that intersect each other. First, we analyze an example of the reference state showing the non-trivial trade-off relation numerically, and find that its eigenvalues must be in a certain range to exhibit the trade-off relation. For another example, one-parameter family of reference states, we analytically show that the non-trivial relation always exists and that the range where the trade-off relation exists is up to about a half of the possible range.

I Introduction

An error trade-off relation upon estimating quantum parametric models is significantly different from the case of classical statistics. There were many examples exhibiting such genuine quantum effects helstrom ; holevo ; yl73 ; RDGill ; gibilisco ; watanabe ; suzuki:ijqi ; RJDD16 ; kull . The usual setting is to estimate the expectation values of two observables, thereby one can compare the derived error trade-off relation to the Heisenberg uncertainty relation. Those previous studies on the trade-off relation focused on the settings where the observables are non-commuting.

In the recent paper sf , we gave an example of physical system that shows a non-trivial trade-off relation between estimation error for the expectation values of two commuting observables. We investigated the uncertainty relation, or trade-off relation obtained by estimating the position of an electron in a uniform magnetic field as a parameter estimation problem of two-parameter unitary model. In this model, shifts in the position of the electron was generated by a unitary transformation with the canonical momenta, pxp_{x} and pyp_{y}. According to quantum mechanics, these generators pxp_{x} and pyp_{y} commute. As the main conclusion of our paper sf , we obtained a trade-off relation between xx and yy, even though the generators of the unitary transformation, pxp_{x} and pyp_{y} commute.

At first sight, this result came out counterintuitive, since two commuting generators should not give any correlation between two parameters of the quantum state. However, we clarified that two parameters are correlated, and hence, we cannot ignore a trade-off relation for our example. To be more precise, two symmetric logarithmic derivative (SLD) operators do not commute in this model, and hence it is genuine quantum mechanical. A natural question is then whether error trade-off relations of this kind exist or not in general, in particular, models of finite dimensional systems. In Ref. kull , for example, the trade-off relation of qubit systems and qutrit systems were investigated. However, neither the existence of the intersection of the SLD and RLD bounds nor its contribution to the trade-off was discussed.

A key observation in our study is that the SLD Cramér-Rao (CR) inequality does not give any trade-off relations, whereas the right logarithmic derivative (RLD) CR inequality does. Importantly, these two CR bounds need to intersect each other in order to show a meaningful trade-off relation unless the RLD CR bound dominates the SLD CR bound. In this way, we can characterize the shape of the error trade-off relation more accurately.

In this paper, we analyze finite dimensional systems to show that there can be such a trade-off relation given by the SLD and RLD CR bounds that intersect each other. As explicit examples, we study qutrit systems to demonstrate this counterintuitive result. We first disprove the existence of such error trade-off relations when the reference state is arbitrary pure states (Sec. III) or mixed qubit states (Sec. III). We then analyze error trade-off relations for a qutrit system numerically by randomly generating reference states. We observe that the occurrence of error trade-off relations is related to the eigenvalues of the reference state. In particular, they have to be distributed equally, otherwise a one particular large, or small, eigenvalue implies no error trade-off relation. To gain more insight into this kinds of trade-off relations, we show analytically that a non-trivial trade-off relation exists in a certain range of the reference state parameter which characterizes the reference state and that the region with the trade-off relation is up to about a half of the allowed region in one of those models.

II Model and error trade-off relation

II.1 Model

Let us consider arbitrary finite dimensional system. We consider the two-parameter unitary transformation with the generators XX and YY , i.e.,

U(θ1,θ2)=eiXθ1iYθ2.U(\theta^{1},\,\theta^{2})=\mathrm{e}^{-\mathrm{i}X\theta^{1}-\mathrm{i}Y\theta^{2}}. (1)

We denote the two-parameter family of states generated from the state ρ0\rho_{0} as ρθ\rho_{\theta}.

ρθ=U(θ1,θ2)ρ0U(θ1,θ2).\rho_{\theta}=U(\theta^{1},\,\theta^{2})\,\rho_{0}U^{\dagger}(\theta^{1},\,\theta^{2}). (2)

The state ρ0\rho_{0} is called as a reference state. In this paper, we mainly consider the case of the commuting generators, [X,Y]=0[X,\,Y]=0 unless stated explicitly.

II.2 Error trade-off relation

In the remaining of the paper, we consider unitary models only. We shall drop the parameter θ\theta to denote the quantum Fisher information matrices, since they are independent of the parameter. To derive a trade-off relation between the diagonal components of the mean square error (MSE) matrix V=[Vij]V=[V_{ij}], suppose we have a quantum CR inequality

VJQ1,V\geq J_{\mathrm{Q}}^{\>-1},

with JQJ_{\mathrm{Q}} a quantum Fisher information matrix. In Ref. sf , we derived a trade-off relation based on the inequality below,

[V11JQ 11][V22JQ 22]>|ImJQ 12| 2.\left[V_{11}-J_{\mathrm{Q}}^{\>11}\right]\left[V_{22}-J_{\mathrm{Q}}^{\>22}\right]>\left|\mathrm{Im}\,J_{\mathrm{Q}}^{\>12}\right|^{\,2}. (3)

where JQ1=[JQij]J_{\mathrm{Q}}^{-1}=[J_{\mathrm{Q}}^{\>ij}]. Im\mathrm{Im} denotes the imaginary part of a complex number. From this expression, we see that there exists a non-trivial error trade-off relation when ImJQ 120\mathrm{Im}\,J_{\mathrm{Q}}^{\>12}\neq 0. Note, however, that this inequality alone does not give a conclusive argument whether an error trade-off relation exists or not. This is because the quantum CR inequality is not tight unless certain special conditions are satisfied. The central idea of this paper is to consider two different CR inequalities set by the SLD and RLD Fisher information matrices. When combining two error trade-off relations, we can determine the shape of an error trade-off relation more accurately.

From the discussion above, we do not have the trade-off relation given by Eq. (3) when the SLD Fisher information matrix JSJ_{\mathrm{S}} is used. This is because it is a real symmetric matrix. The other candidate for giving rise to an error trade-off relation is the RLD Fisher information matrix JRJ_{\mathrm{R}}. In this case, the necessary condition to have an error trade-off relation is

|ImJR 12| 20.\left|\mathrm{Im}\,J_{\mathrm{R}}^{\>12}\right|^{\,2}\neq 0. (4)

as in Eq. (3). By defining δ:=JR, 12JR, 21\delta:=J_{\mathrm{R},\,12}-J_{\mathrm{R},\,21}, we have an equivalent condition,

Condition 1:δ=JR, 12JR, 210,\mathrm{Condition\,1}:\qquad\delta=J_{\mathrm{R},\,12}-J_{\mathrm{R},\,21}\neq 0, (5)

where JR=[JR,ij]J_{\mathrm{R}}=[J_{\mathrm{R},\,ij}].

(i,j)(i,\,j) component of the RLD Fisher information matrix, JR,ijJ_{\mathrm{R},\,ij} is defined as follows.

JR,ij=tr(ρ0LR,jLR,i),J_{\mathrm{R},\,ij}=\mathrm{tr}\left(\rho_{0}L_{\mathrm{R},\,j}\,L^{\dagger}_{\mathrm{R},\,i}\right), (6)

where iρθ|θ=0=ρ0LR,i\partial_{i}\rho_{\theta}\,|_{\theta=0}=\rho_{0}L_{\mathrm{R},\,i}. By using iρθ|θ=0=LR,iρ0\partial_{i}\rho_{\theta}\,|_{\theta=0}=L^{\dagger}_{\mathrm{R},\,i}\,\rho_{0}, Eqs (1), and  (2), we obtain

JR,ij=tr([Xjρ0][Xi,ρ0]ρ01),J_{\mathrm{R},\,ij}=-\mathrm{tr}\left([X_{j}\,\rho_{0}][X_{i},\,\rho_{0}]\,\rho_{0}^{-1}\right), (7)

where X1=XX_{1}=X and X2=YX_{2}=Y. With this, Condition 1 is

δ=tr([[X,ρ0],[Y,ρ0]]ρ01),\delta=\mathrm{tr}\left(\big{[}[X,\,\rho_{0}]\,,\,[Y,\,\rho_{0}]\big{]}\,\rho_{0}^{-1}\right), (8)

and thus, it is relatively easy to check this condition analytically. We stress that having commuting generators, [X,Y]=0[X,Y]=0 does not immediately imply δ=0\delta=0.

Now suppose Condition 1 is satisfied.

[V11JR 11][V22JR 22]>|ImJR 12| 2.\left[V_{11}-J_{\mathrm{R}}^{\>11}\right]\left[V_{22}-J_{\mathrm{R}}^{\>22}\right]>\left|\mathrm{Im}\,J_{\mathrm{R}}^{\>12}\right|^{\,2}. (9)

To give a conclusive argument for the existence of an error trade-off relation, we also consider consider the SLD CR inequality. Since the SLD Fisher information matrix is real, the diagonal components of the MSE matrix obey

V11JS 110,V22JS 220.V_{11}-J_{\mathrm{S}}^{\>11}\geq 0,\ V_{22}-J_{\mathrm{S}}^{\>22}\geq 0. (10)

Note that from the general relationship between the SLD and RLD Fisher information matrices, we have petz

JS1Re(JR1),J_{\mathrm{S}}^{\>-1}\geq\mathrm{Re}\,(J_{\mathrm{R}}^{\>-1}),

where Re\mathrm{Re} denotes the real part of a matrix. Since there exists a locally unbiased estimator such that its MSE for θ1\theta^{1} is arbitrary close to JS 11J_{\mathrm{S}}^{\>11}, (The same statement holds for θ2\theta^{2} as well.) intersections of two inequalities (9) and (10) imply the existence of the error trade-off relation. See Fig. 5 for the occurrence of intersections of the two bounds and the error trade-off relation as an example. Working out elementary algebra, we find that the following condition needs to be satisfied in order to have a non-trivial error trade-off relation sf .

Condition 2:Δ:=|ImJR 12|2[JR 11JS 11][JR 22JS 22]>0.\mathrm{Condition\,2}:\quad\Delta:=\left|\mathrm{Im}\,J_{\mathrm{R}}^{\>12}\right|^{2}-\left[J_{\mathrm{R}}^{\>11}-J_{\mathrm{S}}^{\>11}\right]\left[J_{\mathrm{R}}^{\>22}-J_{\mathrm{S}}^{\>22}\right]>0. (11)

And hence, we have a trade-off relation between V11V_{11} and V22V_{22} if these two conditions (5) and (11) are satisfied.

Next, let us make a remark about D-invariant models. It is known that, the RLD CR inequality is saturated when the model is D-invariant kahn ; hayashi ; kahn2 ; yamagata ; suzuki ; yang . This is true at least in the asymptotic setting. There is no intersection of the RLD and SLD CR bounds in the D-invariant models, because the RLD CR bound is dominant over the SLD CR bound. If the model is D-invariant and if the imaginary part of the off-diagonal elements of RLD Fisher information matrix are not zero, there is a trade-off relation that results from Condition 1 only. In the following, we mainly investigate the non-asymptotic setting unless stated explicitly. This is in contrast to the previous study kull , where the authors focused on the D-invariant model.

III Reference state: Pure state

We first consider the case that the reference state is a pure state, i.e., ρ0=|ψ0ψ0|\rho_{0}=\ket{\psi_{0}}\bra{\psi_{0}}. From Eqs. (1) and (2), ρθ\rho_{\theta} is expressed as

ρθ=eiXθ1iYθ2|ψ0ψ0|eiXθ1+iYθ2.\rho_{\theta}=\mathrm{e}^{-\mathrm{i}X\theta^{1}-\mathrm{i}Y\theta^{2}}\ket{\psi_{0}}\bra{\psi_{0}}\mathrm{e}^{\mathrm{i}X\theta^{1}+\mathrm{i}Y\theta^{2}}. (12)

Therefore, we have

1|ψθ=iX|ψθ,\displaystyle\partial_{1}\ket{\psi_{\theta}}=-\mathrm{i}X\ket{\psi_{\theta}}, (13)
2|ψθ=iY|ψθ,\displaystyle\partial_{2}\ket{\psi_{\theta}}=-\mathrm{i}Y\ket{\psi_{\theta}}, (14)

where |ψθ=eiXθ1iYθ2|ψ0\ket{\psi_{\theta}}=\mathrm{e}^{-\mathrm{i}X\theta^{1}-\mathrm{i}Y\theta^{2}}\ket{\psi_{0}}. In the pure state model, the RLD does not exist in general. Here, we use the generalized RLD instead fujiwara2 . The components of generalized RLD Fisher information matrix, J~R, 12\tilde{J}_{\mathrm{R},\,12} and J~R, 21\tilde{J}_{\mathrm{R},\,21} are given by fujiwara

J~R, 12\displaystyle\tilde{J}_{\mathrm{R},\,12} =4(ψ0|YX|ψ0ψ0|Y|ψ0ψ0|X|ψ0),\displaystyle=4(\braket{\psi_{0}}{Y\,X}{\psi_{0}}-\braket{\psi_{0}}{Y}{\psi_{0}}\braket{\psi_{0}}{X}{\psi_{0}}), (15)
J~R, 21\displaystyle\tilde{J}_{\mathrm{R},\,21} =4(ψ0|XY|ψ0ψ0|X|ψ0ψ0|Y|ψ0).\displaystyle=4(\braket{\psi_{0}}{X\,Y}{\psi_{0}}-\braket{\psi_{0}}{X}{\psi_{0}}\braket{\psi_{0}}{Y}{\psi_{0}}). (16)

We obtain δ\delta in Condition 1, or Eq. 5 as follows.

δ=J~R, 12J~R, 21=4ψ0|[Y,X]|ψ0.\delta=\tilde{J}_{\mathrm{R},\,12}-\tilde{J}_{\mathrm{R},\,21}=4\braket{\psi_{0}}{[Y,\,X]}{\psi_{0}}. (17)

If the reference state is a pure state and if the generators commute, Condition 1 does not hold. Therefore, there is no trade-off relation given by Eq. (3). See also  matsumoto .

IV Reference state: Qubit state

IV.1 General case

We consider the case of a single qubit in a mixed state. We first consider the general two-parameter unitary model to get insight into the problem. By using the Bloch vector, we can express the reference state ρ0\rho_{0} as

ρ0=12(I+s0σ).\rho_{0}=\frac{1}{2}(\mathrm{I}+\vec{s}_{0}\cdot\vec{\sigma}). (18)

where |s0|<1|\vec{s}_{0}|<1. ρθ\rho_{\theta} is given by Eq. (2). The generators X,YX,\,Y can also be expanded with using Pauli matrices.

X\displaystyle X =x0I+xσ,\displaystyle=x_{0}\,\mathrm{I}+\vec{x}\cdot\vec{\sigma}, (19)
Y\displaystyle Y =y0I+yσ.\displaystyle=y_{0}\,\mathrm{I}+\vec{y}\cdot\vec{\sigma}. (20)

The inverse of SLD and RLD Fisher information matrices, JS1J_{\mathrm{S}}^{\>-1} and JR1J_{\mathrm{R}}^{\>-1} are explicitly written as

JS1\displaystyle J_{\mathrm{S}}^{\>-1} =4detJS((y×s0)2(x×s0)(y×s0)(x×s0)(y×s0)(x×s0)2),\displaystyle=\frac{4}{\mathrm{det}\,J_{\mathrm{S}}}\begin{pmatrix}(\vec{y}\times\vec{s}_{0})^{2}&-(\vec{x}\times\vec{s}_{0})\cdot(\vec{y}\times\vec{s}_{0})\\ -(\vec{x}\times\vec{s}_{0})\cdot(\vec{y}\times\vec{s}_{0})&(\vec{x}\times\vec{s}_{0})^{2}\\ \end{pmatrix}, (21)
JR1\displaystyle J_{\mathrm{R}}^{\>-1} =JS1\displaystyle=J_{\mathrm{S}}^{-1}
+4detJS(0is0 2[s0(x×y)]is0 2[s0(x×y)]0),\displaystyle+\frac{4}{\mathrm{det}\,J_{\mathrm{S}}}\begin{pmatrix}0&-\mathrm{i}\vec{s}_{0}^{\>2}[\vec{s}_{0}\cdot(\vec{x}\times\vec{y})]\\ \mathrm{i}\vec{s}_{0}^{\>2}[\vec{s}_{0}\cdot(\vec{x}\times\vec{y})]&0\end{pmatrix}, (22)

where detJS\mathrm{det}\,J_{\mathrm{S}} is the determinant of JSJ_{\mathrm{S}}, and it is

detJS=16s0 2[s0(x×y)]2.\mathrm{det}\,J_{\mathrm{S}}={16\,\vec{s}_{0}^{\>2}[\vec{s}_{0}\cdot(\vec{x}\times\vec{y})]^{2}}. (23)

As shown in Eqs. (21) and (22), JS1=ReJR1J_{\mathrm{S}}^{\>-1}=\mathrm{Re}\,J_{\mathrm{R}}^{\>-1} holds. It follows that our qubit model is D-invariant. (See Lemma III-3 in Ref.suzuki .) Therefore, the RLD CR bound is asymptotically achievable and gives a trade-off relation. In this case, as explained earlier, the SLD and RLD CR bounds do not have intersections, but the trade-off relation exists in the asymptotic setting.

As for the Nagaoka bound, or the Gill-Massar (GM) bound for a two-parameter qubit model, which is known to be achievable nagaoka2 ; RDGill in the non-asymptotic setting, an inequality regarding the diagonal components of the MSE matrix can be derived. The inequality of the Nagaoka band is written as

[V11JS 11][V22JS 22]>1detJS.\left[V_{11}-J_{\mathrm{S}}^{\>11}\right]\left[V_{22}-J_{\mathrm{S}}^{\>22}\right]>\frac{1}{\mathrm{det}\,J_{\mathrm{S}}}. (24)

From Eqs. (10) and (22), we obtain the inequality of the RLD CR bound.

[V11JS 11][V22JS 22]>s0 2detJS.\left[V_{11}-J_{\mathrm{S}}^{\>11}\right]\left[V_{22}-J_{\mathrm{S}}^{\>22}\right]>\frac{\vec{s}_{0}^{\>2}}{\mathrm{det}\,J_{\mathrm{S}}}. (25)

We used JR 11=JS 11J_{\mathrm{R}}^{\>11}=J_{\mathrm{S}}^{\>11} and JR 22=JS 22J_{\mathrm{R}}^{\>22}=J_{\mathrm{S}}^{\>22}. Since s0 2<1\vec{s}_{0}^{\>2}<1, the Nagaoka bound is tighter than the RLD CR bound. This is because the Nagaoka bound is achieved by a separable measurement.

IV.2 Commuting generators’ case

Next, we derive a relationship between XX and YY, or x\vec{x} and y\vec{y} when XX and YY commute. From Eqs. (19) and (20), the commuting relation of XX and YY is given as

[X,Y]=[xσ,yσ]=2i(x×y)σ.[X,\,Y]=[\vec{x}\cdot\vec{\sigma},\,\vec{y}\cdot\vec{\sigma}]=2\mathrm{i}(\vec{x}\times\vec{y})\cdot\vec{\sigma}.

It immediately follows that the necessary and sufficient condition for XX and YY to commute is x×y=0\vec{x}\times\vec{y}=\vec{0}. There is no trade-off relation because the unitary transformation is no longer two-parameter model, because x\vec{x} and y\vec{y} are parallel when XX and YY commute.

V Reference state: Qutrit state

Let us consider a qutrit system, the three-dimensional system. To avoid non-regular models, we consider the full-rank model. Other regularity conditions are also imposed implicitly.

Since XX and YY commute, they are simultaneously diagonalizable. Without the loss of generality, for the calculation of δ\delta, we can use the representation so that both XX and YY can be diagonalized .

ρ0\displaystyle\rho_{0} =(ρ11ρ12ρ13ρ21ρ22ρ23ρ31ρ32ρ33),\displaystyle=\begin{pmatrix}\rho_{11}&\rho_{12}&\rho_{13}\\ \rho_{21}&\rho_{22}&\rho_{23}\\ \rho_{31}&\rho_{32}&\rho_{33}\\ \end{pmatrix}, (26)
X\displaystyle X =(x1000x2000x3),\displaystyle=\begin{pmatrix}x_{1}&0&0\\ 0&x_{2}&0\\ 0&0&x_{3}\\ \end{pmatrix},
Y\displaystyle Y =(y1000y2000y3).\displaystyle=\begin{pmatrix}y_{1}&0&0\\ 0&y_{2}&0\\ 0&0&y_{3}\\ \end{pmatrix}. (27)

By using Eq. (7), δ\delta is calculated as follows.

δ\displaystyle\delta =(detρ0)1(ρ12ρ23ρ31ρ21ρ32ρ13)[(y×x)1],\displaystyle=(\det\rho_{0})^{-1}\,\big{(}\rho_{12}\rho_{23}\rho_{31}-\rho_{21}\rho_{32}\rho_{13}\big{)}\,\left[\big{(}\vec{y}\times\vec{x}\big{)}\cdot\vec{1}\right], (28)

where x=(x1,x2,x3),y=(y1,y2,y3),\vec{x}=(x_{1},x_{2},x_{3}),\>\vec{y}=(y_{1},y_{2},y_{3}), and 1=(1,1,1)\vec{1}=(1,1,1). The condition of no trade-off relation, δ=0\delta=0 holds when

Im(ρ12ρ23ρ31)\displaystyle\mathrm{Im}\left(\rho_{12}\rho_{23}\rho_{31}\right) =0,\displaystyle=0, (29)
or\displaystyle\mathrm{or}
(y×x)1\displaystyle(\vec{y}\times\vec{x})\cdot\vec{1} =0.\displaystyle=0. (30)

Violation of these conditions together with Eq. (11) are the necessary and sufficient conditions to have a non-trivial error trade-off relation. In the case of qutrit, we cannot give an explicit expression of Δ\Delta in general. But, we can obtain Δ\Delta in a straightforward manner numerically.

In the following subsections, we give examples of reference states that give non-trivial error trade-off relations. One of them gives a relatively high possibility. Our main interest is to investigate the error trade-off relation for a given commuting XX and YY.

V.1 Example: reference state with multi-parameter

As one of the simplest examples, we pick an example with pure imaginary off-diagonal components as a reference state ρ0\rho_{0} with five reference state parameters v1v_{1}, v2v_{2}, v3v_{3}, u1u_{1}, u2u_{2}, and u3u_{3}. (v1+v2+v3=1v_{1}+v_{2}+v_{3}=1)

ρ0=13(v1iu1iu2iu1v2iu3iu2iu3v3).\rho_{0}=\frac{1}{3}\begin{pmatrix}v_{1}&-\mathrm{i}\sqrt{u_{1}}&\mathrm{i}\sqrt{u_{2}}\\ \mathrm{i}\sqrt{u_{1}}&v_{2}&-\mathrm{i}\sqrt{u_{3}}\\ -\mathrm{i}\sqrt{u_{2}}&\mathrm{i}\sqrt{u_{3}}&v_{3}\\ \end{pmatrix}. (31)

We choose the reference state ρ0\rho_{0} as above, because imaginary parts of the off-diagonal components of the reference state ρ0\rho_{0} are important to satisfy Condition 1 as seen in Eq. (29).

Refer to caption
Figure 1: Δ\Delta as a function of umaxu_{max}, the maximum of u1,u2u_{1},\,u_{2}, and u3u_{3} in Eq. (31). x=(1,2,3)\vec{x}=(1,~{}2,~{}3), y=(1.5,5,1)\vec{y}=(1.5,~{}5,~{}1).

We calculate Δ\Delta in Condition 2 with using the reference state ρ0\rho_{0} defined by Eq. (31) of which reference state parameters are generated by random numbers. We pick those which satisfy trρ0=1\mathrm{tr}\,\rho_{0}=1 and ρ0>0\rho_{0}>0 and calculate the RLD and SLD Fisher information matrices JSJ_{\mathrm{S}} and JRJ_{\mathrm{R}}. The RLD Fisher information matrix is obtained by using Eq. (7). The SLD Fisher information calculation is done in the standard method. (See for example, Refs. paris ; Liu .) The number of samples generated is on the order of 10610^{6}. Figure 1 shows Δ\Delta as a function of umaxu_{max}, the maximum of u1u_{1}, u2u_{2}, and u3u_{3}. There exists a region Δ>0\Delta>0. The ratio of obtaining Δ>0\Delta>0 out of all of the samples generated is 3.0%. Figures 2 and 3 show Δ\Delta as a function of λmin\lambda_{min} and λmax\lambda_{max}, respectively. λmin\lambda_{min} and λmax\lambda_{max} are the minimum and maximum of eigenvalues of ρ0\rho_{0}, respectively. For Δ\Delta to be positive, λmin\lambda_{min} and λmax\lambda_{max} must be in a certain range. λmin\lambda_{min} is more than about 0.13 and λmax\lambda_{max} is less than about 0.58.

Refer to caption
Figure 2: Δ\Delta as a function of λmin\lambda_{min}, the minimum of the eigenvalues of ρ0\rho_{0}
Refer to caption
Figure 3: Δ\Delta as a function of λmax\lambda_{max}, the maximum of the eigenvalues of ρ0\rho_{0}

V.2 Example: one-parameter family of reference states

Next, we set the reference state parameters in Eq. (31) as v1=v2=v3=1v_{1}=v_{2}=v_{3}=1 and u1=u2=u3=uu_{1}=u_{2}=u_{3}=u in order to investigate the model more in detail analytically. We pick the reference state parameters as above, because the result of Section V.1 indicates that the eigenvalues of the reference state Eq. (31) be roughly in the range 1/3±0.21/3\pm 0.2 to exhibit the non-trivial trade-off relation. The reference state ρ0\rho_{0} is, then explicitly written as

ρ0=13I+13u(0iii0iii0),\rho_{0}=\frac{1}{3}I+\frac{1}{3}\sqrt{u}\begin{pmatrix}0&-\mathrm{i}&\mathrm{i}\\ \mathrm{i}&0&-\mathrm{i}\\ -\mathrm{i}&\mathrm{i}&0\\ \end{pmatrix}, (32)

where II denotes 3×\times3 identity matrix. The reference state ρ0\rho_{0} is a sum of the completely mixed state of the qutrit system and a perturbation with one parameter uu. The parameter uu must be in the range, 0<u<1/30<u<1/3 for the reference state ρ0\rho_{0} to be positive. We exclude u=0u=0, because ρ0=I/3\rho_{0}=\mathrm{I}/3 at u=0u=0.

In the following, we show that the reference state ρ0\rho_{0} Eq. (31) always gives a non-trivial trade-off relation with a certain choice of the reference state parameter uu and that the possibility of seeing the non-trivial trade-off relation is not small.

V.2.1 Intersections of RLD and SLD CR bounds

From Condition 2 expressed by Eq. (11), Δ>0\Delta>0 needs to be satisfied in order to have a non-trivial error trade-off relation. We define a geometrical parameter, ζ\zeta as follows.

ζ=[1(x×y)]2(1×x)2(1×y)2.\zeta=\frac{[\vec{1}\cdot(\vec{x}\times\vec{y})]^{2}}{(\vec{1}\times\vec{x})^{2}(\vec{1}\times\vec{y})^{2}}. (33)

Let ξ=1×x\vec{\xi}=\vec{1}\times\vec{x} and η=1×y\vec{\eta}=\vec{1}\times\vec{y}. A vector analysis formula gives an expression,

ζ=13sin2θ13,\zeta=\frac{1}{3}\sin^{2}\theta\leq\frac{1}{3}, (34)

where sinθ=|ξ×η|/(|ξ||η|)\sin\theta=|\vec{\xi}\times\vec{\eta}|/(|\vec{\xi}||\vec{\eta}|). ζ=1/3\zeta=1/3 when θ=±π/2\theta=\pm{\pi/2}. ζ=0\zeta=0 is excluded, because ξ×η=0\vec{\xi}\times\vec{\eta}=\vec{0} gives δ=0\delta=0 from Eq. (30). Therefore, the possible range for the parameter ζ\zeta is 0<ζ1/30<\zeta\leq 1/3.

We introduce a function of uu at a given ζ\zeta, Fζ(u)F_{\zeta}(u) as

Fζ(u)=16ζ(3u27u+2)2u(3u29u+8)2.F_{\zeta}(u)=16\zeta(3u^{2}-7u+2)^{2}-u(3u^{2}-9u+8)^{2}. (35)

By using Fζ(u)F_{\zeta}(u), Δ\Delta is expressed as

Δ=916ζ2|ξ|2|η|2(2u)u(u27u+4)2Fζ(u).\Delta=\frac{9}{16\zeta^{2}|\vec{\xi}|^{2}|\vec{\eta}|^{2}(2-u)u(u^{2}-7u+4)^{2}}F_{\zeta}(u). (36)

The coefficient of Fζ(u)F_{\zeta}(u) in Eq. (36) is positive finite when 0<u<1/30<u<1/3. In order to investigate the range of uu that gives Δ>0\Delta>0, we can check the condition for Fζ(u)>0F_{\zeta}(u)>0 instead.

We can analytically show that Fζ(u)F_{\zeta}(u) is a monotonically decreasing function of uu and that there is always a unique solution u0u_{0} that satisfies Fζ(u0)=0F_{\zeta}(u_{0})=0 when 0<ζ1/30<\zeta\leq 1/3 and when ρ0>0\rho_{0}>0, i.e., 0<u<1/30<u<1/3. A detailed explanation is given in Appendix A. Figure 4 shows the solution u0u_{0} that satisfies Fζ(u)=0F_{\zeta}(u)=0. In the region where u<u0u<u_{0} at a given ζ\zeta, the non-trivial trade-off relation exists. We can regard u0u_{0} as the upper limit of uu that gives a non-trade off relation. It is worth noting that the upper limit of uu is almost a half of the maximum of uu, 1/31/3 at ζ=1/3\zeta=1/3. This means that the possibility of realizing non-trivial trade-off relation is not small.

Refer to caption
Figure 4: Solutoin u0u_{0} that satisfies Fζ(u0)=0F_{\zeta}(u_{0})=0 in the range 0<ζ1/30<\zeta\leq 1/3

Figure 5 shows an example in which the SLD and RLD CR bounds have two intersections. The parameters used are u=1/12u=1/12, x=(1,2,3)\vec{x}=(1,~{}2,~{}3), and y=(1.5,5,1)\vec{y}=(1.5,~{}5,~{}1).

We can see that Δ1\Delta_{1} and Δ2\Delta_{2} in Fig. 5 indicates the “strength” of trade-off relation by their definitions. They are calculated as

Δ1\displaystyle\Delta_{1} =ΔgS 22gR 22=34ζ|η|2u(u27u+4)(3u29u+8)Fζ(u),\displaystyle=\frac{\Delta}{g_{\mathrm{S}}^{\>22}-g_{\mathrm{R}}^{\>22}}=\frac{3}{4\,\zeta|\vec{\eta}|^{2}u(u^{2}-7u+4)(3u^{2}-9u+8)}F_{\zeta}(u),
Δ2\displaystyle\Delta_{2} =ΔgS 11gR 11=34ζ|ξ|2u(u27u+4)(3u29u+8)Fζ(u).\displaystyle=\frac{\Delta}{g_{\mathrm{S}}^{\>11}-g_{\mathrm{R}}^{\>11}}=\frac{3}{4\,\zeta|\vec{\xi}|^{2}u(u^{2}-7u+4)(3u^{2}-9u+8)}F_{\zeta}(u).

The strengths of trade-off relation is proportional to Δ\Delta. Figure 6 shows Δ1\Delta_{1} and Δ2\Delta_{2} as a function of the parameter uu. In the range where Δ1>0\Delta_{1}>0 or Δ2>0\Delta_{2}>0, the non-trivial trade-off relation exists. The strength of trade-off relation becomes stronger as uu approaches 0.

Refer to caption
Figure 5: Example of RLD and SLD CR bounds with the intersections: the reference state ρ0\rho_{0} defined by Eq. (32) with u=1/12u=1/12, x=(1,2,3)\vec{x}=(1,~{}2,~{}3), y=(1.5,5,1)\vec{y}=(1.5,~{}5,~{}1). The gray region is an allowed region.
Refer to caption
Figure 6: Δ1\Delta_{1} (solid line) and Δ2\Delta_{2} (dotted line) as a function of the parameter uu. x\vec{x} and y\vec{y} are the same as those used for Fig. 5. In the range where Δ1>0\Delta_{1}>0, therefore Δ2>0\Delta_{2}>0, the non-trivial trade-off relation exists.

V.3 Discussion

Unlike a qubit reference state or a pure state reference state, there exists a non-trivial trade-off relation for some qutrit reference states even when the generators commute. We show analytically that a non-trivial trade-off relation always exists in a certain range of the reference state parameter uu when the reference state ρ0\rho_{0} is defined by Eq. (32) that is a sum of the completely mixed state and a perturbation.

Furthermore, the strengths of trade-off relation Δ1\Delta_{1} and Δ2\Delta_{2} increase as uu approaches 0. This looks counterintuitive, because we can regard uu as a small perturbation from 3x3 identity matrix when u1u\ll 1 by the definition of ρ0\rho_{0}, Eq. (32). This reflects the fact that iρθ\partial_{i}\rho_{\theta} is not necessarily small when the perturbation itself is small. Since the (i,j)(i,j) component of the RLD Fisher information matrix is JR,ij=tr[jρθLR,i]J_{\mathrm{R},\,ij}=\mathrm{tr}[\partial_{j}\rho_{\theta}L_{R,\,i}^{\dagger}], the component gR,ijg_{\mathrm{R},\,ij} may not be small if iρθ\partial_{i}\rho_{\theta} is not small.

In a more general case when ρ0\rho_{0} is expressed by Eq. (31), we conducted numerical analysis. In this case also, there exists a non-trivial trade-off relation. Furthermore, in the case of four dimensional system with pure imaginary off-diagonal components, we also see a non-trivial trade-off relation by the same numerical analysis as well. With these, we conclude that the error trade-off relation is a generic phenomenon in the sense that it occurs with a finite volume in the spate space.

VI Conclusion

We have investigated whether the error trade-off relation exists in the generic two-parameter unitary models for finite dimensional systems with the commuting generators. By analyzing the necessary and sufficient conditions for the SLD and RLD CR bounds to intersect each other, we obtain the necessary and sufficient conditions for the existence of a non-trivial trade-off relation based on the SLD and RLD CR bounds for arbitrary finite dimensional system.

By using the conditions, we show two examples of the qutrit system with the non-trivial trade-off relation. The result of the reference state with multi-parameter indicates that the eigenvalues of the reference state be in a certain range. In the other model reference state with one-parameter, we show analytically that a non-trivial trade-off relation always exists in a certain range of the reference state parameter and that the region with the trade-off relation is up to about a half of the allowed region.

In our previous study about the trade-off relation of an infinite dimensional system sf , the bound is also given by both of the SLD and RLD CR bounds when the generators of the unitary transformation with the commuting generators. As shown in Figs. 5 and 6, we confirmed that what we saw in our previous study is not special, but generic. When the reference state is a pure state or a general qubit state, we disprove the existence of a non-trivial trade-off relation.

Acknowledgment

The work is partly supported by the FY2020 UEC Research Support Program, the University of Electro-Communications.

Appendix A Solution u0u_{0} of Fζ(u0)=0F_{\zeta}(u_{0})=0

In this section, we investigate the solution s0s_{0} of Fζ(u0)=0F_{\zeta}(u_{0})=0. We check up to the fourth partial derivative of Fζ(u)F_{\zeta}(u) with respect to ss to see Fζ(u)F_{\zeta}(u) in the allowed range for tt and ss.

Let Fζ(n)(u)=dnFζ(u)unF_{\zeta}^{(n)}(u)=\dfrac{d^{n}F_{\zeta}(u)}{\partial u^{n}}. Up to the fourth partial derivative of Fζ(u)F_{\zeta}(u) with respect to uu are as follows

Fζ(1)(u)\displaystyle F_{\zeta}^{(1)}(u) =45u464(1+7ζ)+72u3(3+8ζ)\displaystyle=-45u^{4}-64(1+7\zeta)+72u^{3}(3+8\zeta)
+32u(9+61ζ)9u2(43+224ζ),\displaystyle+32u(9+61\zeta)-9u^{2}(43+224\zeta),
Fζ(2)(u)\displaystyle F_{\zeta}^{(2)}(u) =2[90u3+108u2(3+8ζ)+16(9+61ζ)\displaystyle=2[-90u^{3}+108u^{2}(3+8\zeta)+16(9+61\zeta)
9u(43+224ζ)],\displaystyle-9u(43+224\zeta)],
Fζ(3)(u)\displaystyle F_{\zeta}^{(3)}(u) =18[43+30u2+224ζ24u(3+8ζ)],\displaystyle=-18[43+30u^{2}+224\zeta-24u(3+8\zeta)],
Fζ(4)(u)\displaystyle F_{\zeta}^{(4)}(u) =216(6+5u16ζ).\displaystyle=-216(-6+5u-16\zeta).

Fζ(3)(u)F_{\zeta}^{(3)}(u) is convex upward, because the coefficient of u2u^{2} in Fζ(3)(u)F_{\zeta}^{(3)}(u) is negative. Therefore, the extremum, in this case, the maximum of Fζ(3)(u)F_{\zeta}^{(3)}(u) is given by u0(4)u^{(4)}_{0} which is the solution of Fζ(4)(u0(4))=0F_{\zeta}^{(4)}(u^{(4)}_{0})=0. The solution u0(4)u^{(4)}_{0} is given by

u0(4)=25(8ζ+3).u^{(4)}_{0}=\frac{2}{5}(8\zeta+3).

u0(4)u^{(4)}_{0} which gives the maximum of Fζ(3)(u)F_{\zeta}^{(3)}(u) becomes minimum at ζ=0\zeta=0. At ζ=0\zeta=0, u0(4)=6/5=1.2>1/3u^{(4)}_{0}=6/5=1.2>1/3. Because of u0(4)>1/3u^{(4)}_{0}>1/3, Fζ(3)(u)F_{\zeta}^{(3)}(u) increases monotonically in the range 0<u<1/30<u<1/3.

Fζ(3)(u)F_{\zeta}^{(3)}(u) at u=1/3u=1/3 is

Fζ(3)(13)\displaystyle F_{\zeta}^{(3)}(\frac{1}{3}) =6(67+480ζ)<0when(0<ζ<13).\displaystyle=-6(67+480\zeta)<0\quad\mathrm{when}\quad(0<\zeta<\frac{1}{3}).

Then, we see Fζ(3)(u)<0F_{\zeta}^{(3)}(u)<0 when 0<u<1/30<u<1/3. Therefore, Fζ(2)(u)F_{\zeta}^{(2)}(u) decreases monotonically when 0<u<1/30<u<1/3.

Fζ(2)(u)F_{\zeta}^{(2)}(u) at u=1/3u=1/3 is

Fζ(2)(13)\displaystyle F_{\zeta}^{(2)}(\frac{1}{3}) =2863+800ζ>0when(0<ζ<13).\displaystyle=\frac{286}{3}+800\zeta>0\quad\mathrm{when}\quad(0<\zeta<\frac{1}{3}).

Therefore, Fζ(1)(u)F_{\zeta}^{(1)}(u) increases monotonically when 0<u<1/30<u<1/3.

Fζ(1)(13)\displaystyle F_{\zeta}^{(1)}(\frac{1}{3}) =329<0when(0<ζ<13).\displaystyle=-\frac{32}{9}<0\quad\mathrm{when}\quad(0<\zeta<\frac{1}{3}).

Therefore, Fζ(u)F_{\zeta}(u) decreases monotonically when 0<u<1/30<u<1/3. The values of Fζ(u)F_{\zeta}(u) at the both ends, u=0u=0 and u=1/3u=1/3 are

Fζ(0)\displaystyle F_{\zeta}(0) =64ζ,\displaystyle=64\zeta,
Fζ(13)\displaystyle F_{\zeta}(\frac{1}{3}) =25627.\displaystyle=-\frac{256}{27}.

With a given ζ\zeta in the range 0<ζ1/30<\zeta\leq 1/3, there always exists only one solution u0u_{0} that satisfies Fζ(u0)=0F_{\zeta}(u_{0})=0 in the range 0<u01/30<u_{0}\leq{1}/{3}.

References

  • (1) C. W. Helstrom Quantum Detection and Estimation Theory, Academic, New York (1976).
  • (2) A. S. Holevo, Probabilistic and Statistical Aspects of Quantum Theory, Edizioni della Normale, Pisa, 2nd ed (2011).
  • (3) H. Yuen and M. Lax, IEEE Trans. on Information Theory, IT19, 740 (1973).
  • (4) R. D. Gill, S. Massar Phys. Rev. Vol A 61, 042312 (2000). https://doi.org/10.1103/PhysRevA.61.042312
  • (5) P. Gibilisco, H. Hiai, D. Petz, IEEE Trans. Information Theory. A, Vol. 55, 439 (2009). https://doi.org/10.1109/TIT.2008.2008142
  • (6) Y. Watanabe, T. Sagawa, M. Ueda, Phys. Rev. A, Vol. 81, 042121 (2011). https://doi.org/10.1103/PhysRevA.84.042121
  • (7) J. Suzuki, Int. J. Quant. Inf. Vol. 9 1167 (2011). https://doi.org/10.1142/S0219749914500440
  • (8) S. Ragy, M. Jarzyna, and R. Demkowicz-Dobrzański, Phys. Rev. A, 94, 052108 (2016). https://doi.org/10.1103/PhysRevA.94.052108
  • (9) I. Kull, P. Allard Guérin, F. Verstraete, J. Phys. A: Math. Theor. 53 (2020) 244001. https://doi.org/10.1088/1751-8121/ab7f67
  • (10) S. Funada and J. Suzuki, Physica A Vol. 558, 124918 (2020). https://doi.org/10.1016/j.physa.2020.124918
  • (11) D. Petz, Linear Algebra Appl. Vol. 244, 81 (1996).
  • (12) M. Guţă and J. Kahn, Phys. Rev. A 73, 052108 (2006).
  • (13) M. Hayashi, K. Matsumoto, J. Math. Phys. Vol. 49, 102101 (2008). https://doi.org/10.1063/1.2988130
  • (14) J.  Kahn, M. Guţă 2009 Commun. Math. Phys. Vol. 289, 597 (2009). https://doi.org/10.1007/s00220-009-0787-3
  • (15) K. Yamagata, A. Fujiwara, R. D. Gill, Ann. Stat. Vol. 41, 2197 (2013). https://doi:10.1214/13-AOS1147
  • (16) J. Suzuki, J. Math. Phys., Vol. 57, 042201 (2016). https://doi.org/10.1063/1.4945086
  • (17) Y. Yang, G. Chiribella, and M. Hayashi, Comm. Math. Phys. Vol. 368, 223 (2019).
  • (18) A. Fujiwara, H. Nagaoka J. Math. Phys., Vol. 40, 4227 (1999). https://doi.org/10.1063/1.532962
  • (19) A. Fujiwara, H. Nagaoka Phys. Lett. A, Vol. 201, 119 (1995). https://doi.org/10.1016/0375-9601(95)00269-9
  • (20) K. Matsumoto J. Phys. A: Math. Gen., Vol. 35, 3111 (2002). https://doi.org/10.1088%2F0305-4470%2F35%2F13%2F307
  • (21) H. Nagaoka, in Asymptotic Theory of Quantum Statistical Inference, ed M Hayashi, Singapore: World Scientific p. 100–112 (2005)
  • (22) M.G.A Paris Int. J. Quantum Inf. 7 , 125 (2009). https://doi.org/10.1142/S0219749909004839
  • (23) J. Liu, H. Yuan, X.-M. Lu, and X. Wang, Q J. Phys. A: Math. Theor. 53 023001 (2020). https://doi.org/10.1088/1751-8121/ab5d4d