Next Article in Journal
SU(2) Decomposition for the Quantum Information Dynamics in 2d-Partite Two-Level Quantum Systems
Next Article in Special Issue
The Maximum Entropy Method in Ultrasonic Non-Destructive Testing—Increasing the Resolution, Image Noise Reduction and Echo Acquisition Rate
Previous Article in Journal
New Estimations for Shannon and Zipf–Mandelbrot Entropies
Previous Article in Special Issue
Thermodynamic Analysis of Irreversible Desiccant Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Information Geometry of Randomized Quantum State Tomography

1
Department of Mathematics, Osaka University, Toyonaka, Osaka 560-0043, Japan
2
Graduate School of Informatics and Engineering, The University of Electro-Communications, Chofu, Tokyo 182-8585, Japan
*
Author to whom correspondence should be addressed.
Submission received: 29 June 2018 / Revised: 5 August 2018 / Accepted: 13 August 2018 / Published: 16 August 2018
(This article belongs to the Special Issue Entropy: From Physics to Information Sciences and Geometry)

Abstract

:
Suppose that a d-dimensional Hilbert space H C d admits a full set of mutually unbiased bases | 1 ( a ) , , | d ( a ) , where a = 1 , , d + 1 . A randomized quantum state tomography is a scheme for estimating an unknown quantum state on H through iterative applications of measurements M ( a ) = | 1 ( a ) 1 ( a ) | , , | d ( a ) d ( a ) | for a = 1 , , d + 1 , where the numbers of applications of these measurements are random variables. We show that the space of the resulting probability distributions enjoys a mutually orthogonal dualistic foliation structure, which provides us with a simple geometrical insight into the maximum likelihood method for the quantum state tomography.

1. Introduction

Quantum state tomography is a method of estimating an unknown quantum state represented on some Hilbert space H , consisting of a fixed set of measurements that provides sufficient information about the unknown quantum state, as well as a data processing that maps each measurement outcome into the quantum state space S ( H ) on H [1]. A set of measurements that fulfils this requirement is sometimes called a measurement basis. For mathematical simplicity, we restrict ourselves to Hilbert spaces of finite dimensions.
To elucidate our motivation, let us treat the simplest case when H C 2 . It is well known that there is a one-to-one affine correspondence between the qubit state space
S ( C 2 ) : = { ρ C 2 × 2 | ρ 0 , Tr ρ = 1 }
and the unit ball (called the Bloch ball)
B : = x = ( x 1 , x 2 , x 3 ) R 3 x 2 : = ( x 1 ) 2 + ( x 2 ) 2 + ( x 3 ) 2 1 .
In fact, the correspondence is explicitly given by the Stokes parametrization
x ρ x = 1 2 ( I + x 1 σ 1 + x 2 σ 2 + x 3 σ 3 ) ,
where σ 1 , σ 2 , and σ 3 are the standard Pauli matrices. Since E ρ x [ σ i ] : = Tr ρ x σ i = x i for i { 1 , 2 , 3 } , the set σ = ( σ 1 , σ 2 , σ 3 ) of observables is regarded as an unbiased estimator [2,3,4] for the Stokes parameter x = ( x 1 , x 2 , x 3 ) . This is the basic idea behind the standard qubit state tomography, which runs as follows: suppose that, among N independent experiments, the ith Pauli matrix σ i was measured N / 3 times, and outcomes + 1 (spin-up) and 1 (spin-down) were obtained n i + and n i times, respectively. Then a naive estimate for the true value of the parameter x = ( x 1 , x 2 , x 3 ) is
x ^ = ( x ^ 1 , x ^ 2 , x ^ 3 ) : = n 1 + n 1 N / 3 , n 2 + n 2 N / 3 , n 3 + n 3 N / 3 .
When the estimate x ^ [ 1 , 1 ] 3 falls outside the Bloch ball B, it needs to be corrected so that the new estimate lies in the Bloch ball B. The maximum likelihood method is a canonical one to obtain a corrected estimate [2,5,6,7,8,9,10]. From the point of view of information geometry [11,12,13], the maximum likelihood estimate (MLE) is the orthogonal projection from the temporary estimate x ^ onto the Bloch ball B with respect to the standard Fisher metric along the ( m ) -geodesic [14], (cf., Appendix A).
Now let us deal with a slightly generalized situation: suppose that the ith Pauli matrix σ i was measured N i times and outcomes + 1 and 1 were obtained n i + and n i times, respectively, where { N i } i = 1 , 2 , 3 were random variables. Such a situation arises in an actual experiment due to unexpected particle loss [15]. We shall call such a generalized estimation scheme a randomized state tomography. A naive estimate in this case is the following:
x ^ = ( x ^ 1 , x ^ 2 , x ^ 3 ) : = n 1 + n 1 N 1 , n 2 + n 2 N 2 , n 3 + n 3 N 3 .
One may invoke the maximum likelihood method when x ^ falls outside the Bloch ball. It is then interesting to ask if there is also a useful geometrical picture for the MLE even when the numbers N i of measurements are random variables.
The above mentioned problem is naturally extended to quantum state tomography on an arbitrary Hilbert space that admits a full set of mutually unbiased bases [16,17]. In a d-dimensional Hilbert space H C d , k orthonormal bases
| α ( 1 ) α { 1 , , d } , | β ( 2 ) β { 1 , , d } , , | γ ( k ) γ { 1 , , d }
are called mutually unbiased if they satisfy
α ( a ) | β ( b ) 2 = 1 d
for all a , b { 1 , , k } with a b , and α , β { 1 , , d } . It is known that the number k of mutually unbiased bases (MUBs) is at most d + 1 [18]. If there are d + 1 MUBs, the Hilbert space H is said to admit a full set of MUBs. For example, when the dimension d of H is a power of a prime, H admits a full set of MUBs [19]. Whether or not any Hilbert space admits a full set of MUBs is an open question [16].
In what follows, unless otherwise stated, we assume that the Hilbert space H C d under consideration admits a full set of MUBs. As demonstrated in Appendix B (cf., [17,20]), each density operator ρ S ( H ) can be uniquely represented as
ρ = ρ ( ξ ) : = a = 1 d + 1 α = 1 d 1 ξ α ( a ) M α ( a ) + 1 β = 1 d 1 ξ β ( a ) M d ( a ) I ,
where
M ( a ) : = M 1 ( a ) , , M d ( a ) = | 1 ( a ) 1 ( a ) | , , | d ( a ) d ( a ) |
is the projection-valued measure (PVM) associated with the ath orthogonal basis in the MUBs, and
ξ : = ξ α ( a ) ( a , α ) { 1 , , d + 1 } × { 1 , , d 1 }
is a ( d 2 1 ) -dimensional real parameter that is chosen so that ρ ( ξ ) 0 . A simple calculation shows that, if the ath measurement M ( a ) is applied to the state ρ ( ξ ) , one obtains each outcome α { 1 , , d } with probability
p α ( a ) = Tr ρ ( ξ ) M α ( a ) = ξ α ( a ) , for α = 1 , , d 1 , 1 β = 1 d 1 ξ β ( a ) , for α = d .
This implies that the parametrization ξ ρ ( ξ ) establishes an affine isomorphism between the quantum state space
S ( C d ) : = ρ C d × d | ρ 0 , Tr ρ = 1
and the convex set
B : = ξ R d 2 1 | ρ ( ξ ) 0 .
Incidentally, the Stokes parametrization x ρ x for the qubit state space S ( C 2 ) is regarded as a special case of the above parametrization ξ ρ ( ξ ) for S ( C d ) . In fact, the eigenvectors of the Pauli matrices σ 1 , σ 2 , σ 3 form a full set of MUBs on C 2 , and the Stokes parametrization x = ( x 1 , x 2 , x 3 ) is related to the above parametrization ξ = ( ξ 1 ( 1 ) , ξ 1 ( 2 ) , ξ 1 ( 3 ) ) as
ξ 1 ( a ) = x a + 1 2 , ( a = 1 , 2 , 3 ) .
Now that a standard affine parametrization ξ ρ ( ξ ) has been established on an arbitrary Hilbert space H C d that admits a full set of MUBs, the scheme of randomized state tomography is naturally extended to H as follows. Suppose that the ath measurement M ( a ) was applied N ( a ) times and the outcome α { 1 , , d } was obtained n α ( a ) times, where N ( a ) a = 1 , , d + 1 were random variables. Then, due to (2), a naive estimate for the parameter ξ α ( a ) is
ξ ^ α ( a ) = n α ( a ) N ( a ) .
When the estimate ξ ^ : = ξ ^ α ( a ) [ 0 , 1 ] d 2 1 falls outside the parameter space B, one may invoke the maximum likelihood method to obtain a corrected estimate.
The objective of the present paper is to clarify that the ( m ) -projection interpretation for the MLE is still valid for the randomized state tomography by changing the standard Fisher metric into a deformed one depending on the realization of the random variables N ( a ) , which might as well be called a randomized Fisher metric. Such a novel geometrical picture will provide important insights into the quantum metrology.
The paper is organized as follows. In Section 2, we first introduce a statistical model on an extended sample space Ω that represents the randomized state tomography. We then clarify that the probability simplex P ( Ω ) is decomposed into mutually orthogonal dualistic foliation by means of certain ( m ) - and ( e ) -autoparallel submanifolds. In Section 3, we give a statistical interpretation of the above-mentioned dualistic foliation structure. In particular, we point out that the MLE is the ( m ) -projection with respect to a deformed Fisher metric that depends on the realization of the random variables N ( a ) . These results are demonstrated by several illustrative examples in Section 4. Finally, some concluding remarks are presented in Section 5. For the reader’s convenience, some background information is provided in Appendix A and Appendix B, including information geometry of the MLE and affine parametrization of a quantum state space S ( H ) .

2. Geometry of Randomized State Tomography

We identify the randomized state tomography on H C d with the following scheme [21]: at each step of the measurement, one chooses a PVM M ( a ) at random with probability s ( a ) , ( a = 1 , , d + 1 ), and applies the chosen PVM to yield an outcome α 1 , , d . The sample space Ω for this statistical picture is
Ω = ( a , α ) | a { 1 , , d + 1 } , α { 1 , , d } .
Suppose that the unknown state ρ is specified by the coordinate ξ B as (1). Then the corresponding probability distribution on Ω is represented by the d ( d + 1 ) -dimensional probability vector
p ( s , ξ ) : = s ( 1 ) ξ 1 ( 1 ) , , ξ d 1 ( 1 ) , 1 α = 1 d 1 ξ α ( 1 ) , , s ( d ) ξ 1 ( d ) , , ξ d 1 ( d ) , 1 α = 1 d 1 ξ α ( d ) , 1 a = 1 d s ( a ) ξ 1 ( d + 1 ) , , ξ d 1 ( d + 1 ) , 1 α = 1 d 1 ξ α ( d + 1 )
where the parameter s : = ( s ( 1 ) , , s ( d ) ) belongs to the domain
D : = s R d | s ( a ) > 0 for a { 1 , , d } , and a = 1 d s ( a ) < 1 .
Note that the family
p ( s , ξ ) | s D , ξ Ξ
with
Ξ : = ξ R d 2 1 | ξ α ( a ) > 0 for ( a , α ) { 1 , , d + 1 } × { 1 , , d 1 } , α = 1 d 1 and α = 1 d 1 ξ α ( a ) < 1 for a { 1 , , d + 1 }
forms a ( d 2 + d 1 ) -dimensional open probability simplex P ( Ω ) , and the parameters ( s , ξ ) form a coordinate system of P ( Ω ) . Since we are only interested in estimating the parameter ξ Ξ , the remaining parameter s D is understood as a set of nuisance parameters [2,12]. In what follows, we regard P ( Ω ) as a statistical manifold endowed with the standard dualistic structure g , ( e ) , ( m ) , where g is the Fisher metric, and ( e ) and ( m ) are the exponential and mixture connections [12].
Let us consider the following submanifolds of P ( Ω ) :
M ( s ) : = p ( s , ξ ) ξ Ξ
for each s D , and
E ( ξ ) : = p ( s , ξ ) s D
for each ξ Ξ . Since M ( s ) and E ( ξ ) are convex subsets of P ( Ω ) , they are both ( m ) -autoparallel. In addition, we have the following.
Proposition 1.
For each ξ Ξ , the submanifold E ( ξ ) is ( e ) -autoparallel. Furthermore, for each s D and ξ Ξ , the submanifolds M ( s ) and E ( ξ ) are mutually orthogonal with respect to the Fisher metric g.
Proof. 
Let us change the coordinate system ( s , ξ ) into ( η a , η b , α ) , where
η a : = s ( a )
for a { 1 , , d } , and
η b , α : = s ( b ) ξ α ( b )
for ( b , α ) { 1 , , d + 1 } × { 1 , , d 1 } . With this coordinate transformation, the probability vector p ( s , ξ ) is rewritten as
p η = a = 1 d + 1 η a , 1 , , η a , d 1 , η a α = 1 d 1 η a , α .
Here, η d + 1 is a function of { η a } a { 1 , , d } defined by
η d + 1 : = 1 a = 1 d η a ,
and is not a component of the coordinate system η : = ( η a , η b , α ) . We see from the representation (3) that the coordinate system η is ( m ) -affine. The potential function for η is given by the negative entropy
φ ( η ) : = ω Ω p η ( ω ) log p η ( ω ) = a = 1 d + 1 α = 1 d 1 η a , α log η a , α + η a β = 1 d 1 η a , β log η a β = 1 d 1 η a , β
and the dual ( e ) -affine coordinate system θ is given by
θ a = φ η a = log s ( a ) 1 b = 1 d s ( b ) + log 1 β = 1 d 1 ξ β ( a ) 1 β = 1 d 1 ξ β ( d + 1 )
for a { 1 , , d } , and
θ b , α = φ η b , α = log ξ α ( b ) 1 β = 1 d 1 ξ β ( b )
for ( b , α ) { 1 , , d + 1 } × { 1 , , d 1 } . Thus, fixing ξ is equivalent to fixing the coordinates θ b , α ( b , α ) { 1 , , d + 1 } × { 1 , , d 1 } , and the submanifold E ( ξ ) is generated by changing the remaining parameters θ a a { 1 , , d } . This implies that E ( ξ ) is ( e ) -autoparallel, proving the first part of the claim.
To prove the second part, let us introduce a mixed coordinate system [11]
η a ; θ b , α a { 1 , , d } , ( b , α ) { 1 , , d + 1 } × { 1 , , d 1 }
of P ( Ω ) . Since η a = s ( a ) , the submanifold M ( s ) is rewritten as
M ( s ) = p ( s , ξ ) | η a a { 1 , , d } are fixed and θ b , α ( b , α ) { 1 , , d + 1 } × { 1 , , d 1 } are arbitrary .
On the other hand, the submanifold E ( ξ ) is rewritten as
E ( ξ ) = p ( s , ξ ) | θ b , α ( b , α ) { 1 , , d + 1 } × { 1 , , d 1 } are fixed and η a a { 1 , , d } are arbitrary .
Thus, the orthogonality of M ( s ) and E ( ξ ) is an immediate consequence of the orthogonality of the dual affine coordinate systems θ and η with respect to the Fisher metric g. ☐
Proposition 1 implies that the manifold P ( Ω ) is decomposed into mutually orthogonal dualistic foliation based on the submanifolds M ( s ) and E ( ξ ) , as illustrated in Figure 1. We shall exploit this geometrical structure in the next section.

3. Estimation of the Parameter ξ

Let us proceed to the problem of estimating the unknown parameter ξ using the randomized tomography. Suppose that, among N independent repetitions of experiments, the ath measurement M ( a ) was applied N ( a ) times and outcomes α { 1 , , d } were obtained n α ( a ) times. Then temporary estimates ( s ^ , ξ ^ ) for the parameters ( s , ξ ) are given by
s ^ ( a ) : = N ( a ) N
for a { 1 , , d } , and
ξ ^ β ( b ) : = n β ( b ) N ( b )
for ( b , β ) { 1 , , d + 1 } × { 1 , , d 1 } . If ξ ^ has fallen outside the physical domain B, one may seek a corrected estimate by the maximum likelihood method. Observe that, due to (2), the empirical distribution q ^ N P ( Ω ) is represented as
q ^ N = p ( s ^ , ξ ^ ) .
On the other hand, the physical domain B in the parameter space Ξ corresponds to the subset
B : = { p ( s , ξ ) | s D , ξ B }
of P ( Ω ) , (see Figure 1). The MLE p * in P ( Ω ) is then given by
p * = arg min p B D ( q ^ N p ) ,
where D ( · · ) is the Kullback-Leibler divergence (cf., Appendix A). A crucial observation is the following.
Proposition 2.
The minimum in (5) is achieved on M ( s ^ ) B .
Proof. 
Let us take a point p ( s , ξ ) B arbitrarily. It then follows from the mutually orthogonal dualistic foliation of P ( Ω ) established in Proposition 1 that
D ( q ^ N p ( s , ξ ) ) = D ( p ( s ^ , ξ ^ ) p ( s , ξ ) ) = D ( p ( s ^ , ξ ^ ) p ( s ^ , ξ ) ) + D ( p ( s ^ , ξ ) p ( s , ξ ) ) D ( p ( s ^ , ξ ^ ) p ( s ^ , ξ ) ) .
In the second equality, the generalized Pythagorean theorem was used. Consequently,
min ξ B D ( p ( s ^ , ξ ^ ) p ( s , ξ ) ) min ξ B D ( p ( s ^ , ξ ^ ) p ( s ^ , ξ ) )
for all s D , and the right-hand side is achieved if and only if s = s ^ . ☐
The geometrical implication of Proposition 2 is illustrated in Figure 2. The MLE p * = p ( s ^ , ξ * ) is the ( m ) -projection from the empirical distribution p ( s ^ , ξ ^ ) to B , and is on the section M ( s ^ ) specified by the temporary estimate s ^ .
Now we arrive at a geometrical picture behind the parameter estimation based on randomized state tomography. Suppose we are given a temporary estimate ( s ^ , ξ ^ ) with ξ ^ B . Due to Proposition 2, we can restrict ourselves to section M ( s ^ ) as the search space for the MLE p * . Since each section M ( s ^ ) is affinely isomorphic to the parameter space Ξ , we can introduce a dualistic structure ( g ˜ , ˜ ( e ) , ˜ ( m ) ) on Ξ in the following way. Firstly, we identify the metric g ˜ with the Fisher metric g restricted on M ( s ^ ) , that is,
g ˜ ( s ^ , ξ ) ξ α ( a ) , ξ β ( b ) = η a , α ξ α ( a ) η b , β ξ β ( b ) g ( s , ξ ) η a , α , η b , β s = s ^ = s ( a ) s ( b ) 2 φ ( η ) η a , α η b , β s = s ^ = δ a b s ^ ( a ) 1 ξ d ( a ) + δ α β ξ α ( a ) ,
for a , b { 1 , , d + 1 } and α , β { 1 , , d 1 } , where s ^ ( d + 1 ) and ξ d ( a ) are formally defined as
s ^ ( d + 1 ) : = 1 a = 1 d s ^ ( a ) , ξ d ( a ) : = 1 α = 1 d 1 ξ α ( a ) .
Secondly, the mixture connection ˜ ( m ) on Ξ is defined through the natural affine isomorphism between M ( s ^ ) and Ξ . Finally, the dual connection ˜ ( e ) is defined by the duality
g ˜ ( ˜ X ( e ) Y , Z ) : = X g ˜ ( Y , Z ) g ˜ ( Y , ˜ X ( m ) Z ) .
Thus, the MLE ξ * in the parameter space Ξ is interpreted as the ˜ ( m ) -projection from ξ ^ to the physical domain B with respect to the metric g ˜ .

4. Examples

In this section, we present some examples that demonstrate the implication of Proposition 2 as well as the general diagram given in Figure 2.

4.1. When dim H = 2

Let us first study the simplest case when H = C 2 . A full set of MUBs is given by
| 1 ( 1 ) , | 2 ( 1 ) = 1 2 1 1 , 1 2 1 1 , | 1 ( 2 ) , | 2 ( 2 ) = 1 2 1 i , 1 2 1 i , | 1 ( 3 ) , | 2 ( 3 ) = 1 0 , 0 1 .
With these bases, the parameter representation (1) becomes
ρ = 1 2 1 + x 3 x 1 i x 2 x 1 + i x 2 1 x 3 ,
where x = ( x 1 , x 2 , x 3 ) is the standard Stokes parameter, which is related to ξ = ( ξ 1 ( 1 ) , ξ 1 ( 2 ) , ξ 1 ( 3 ) ) as x a = 2 ξ 1 ( a ) 1 for a = 1 , 2 , 3 .
Figure 3 demonstrates how the ˜ ( m ) -projection is realized. Here, the trajectories of ˜ ( m ) -projections that gives the MLE p * are plotted only on the x 1 x 2 -plane. The left and right panels correspond to the cases when N ( 1 ) : N ( 2 ) = 1 : 1 and N ( 1 ) : N ( 2 ) = 5 : 1 , respectively. The change of ξ 1 -coordinate relative to the change of x 2 -coordinate along each trajectory is less noticeable in the right panel than in the left panel. This is because a tomography with N ( 1 ) / N ( 2 ) = 5 provides us with more information about x 1 -coordinate, relative to x 2 -coordinate, as compared with the case when N ( 1 ) / N ( 2 ) = 1 .

4.2. When dim H = 3

The space H = C 3 admits a full set of MUBs; for example,
| 1 ( 1 ) , | 2 ( 1 ) , | 3 ( 1 ) = 1 0 0 , 0 1 0 , 0 0 1 , | 1 ( 2 ) , | 2 ( 2 ) , | 3 ( 2 ) = 1 3 1 1 1 , 1 3 1 ω ω 2 , 1 3 1 ω 2 ω , | 1 ( 3 ) , | 2 ( 3 ) , | 3 ( 3 ) = 1 3 ω 1 1 , 1 3 1 ω 1 , 1 3 1 1 ω , | 1 ( 4 ) , | 2 ( 4 ) , | 3 ( 4 ) = 1 3 ω 2 1 1 , 1 3 1 ω 2 1 , 1 3 1 1 ω 2 ,
where ω = ( 1 + i 3 ) / 2 is a primitive third root of unity. With these bases, the parameter representation (1) becomes
ρ = ξ 1 ( 1 ) a 12 i b 12 a 13 i b 13 a 12 + i b 12 ξ 2 ( 1 ) a 23 i b 23 a 13 + i b 13 a 23 + i b 23 1 ξ 1 ( 1 ) ξ 2 ( 1 ) ,
where
a 12 = 1 2 1 + ξ 1 ( 2 ) ξ 1 ( 3 ) ξ 2 ( 3 ) ξ 1 ( 4 ) ξ 2 ( 4 ) , a 13 = 1 2 1 + ξ 1 ( 2 ) + ξ 2 ( 3 ) + ξ 2 ( 4 ) , a 23 = 1 2 1 + ξ 1 ( 2 ) + ξ 1 ( 3 ) + ξ 1 ( 4 ) , b 12 = 3 6 1 ξ 1 ( 2 ) 2 ξ 2 ( 2 ) + ξ 1 ( 3 ) ξ 2 ( 3 ) ξ 1 ( 4 ) + ξ 2 ( 4 ) , b 13 = 3 6 1 + ξ 1 ( 2 ) + 2 ξ 2 ( 2 ) + 2 ξ 1 ( 3 ) + ξ 2 ( 3 ) 2 ξ 1 ( 4 ) ξ 2 ( 4 ) , b 23 = 3 6 1 ξ 1 ( 2 ) 2 ξ 2 ( 2 ) + ξ 1 ( 3 ) + 2 ξ 2 ( 3 ) ξ 1 ( 4 ) 2 ξ 2 ( 4 ) .
The physical domain B that corresponds to the state space S ( C 3 ) is a compact convex subset of the parameter space Ξ ( R 8 ) , and the extreme points of B form an algebraic variety with respect to the parameters
ξ = ( ξ 1 ( 1 ) , ξ 2 ( 1 ) , ξ 1 ( 2 ) , ξ 2 ( 2 ) , ξ 1 ( 3 ) , ξ 2 ( 3 ) , ξ 1 ( 4 ) , ξ 2 ( 4 ) ) .
A numerical example of a ˜ ( m ) -projection that gives the MLE is illustrated in Figure 4, where no probe particle is lost, that is, when
s ^ = 1 4 , 1 4 , 1 4 , 1 4 .
In Figure 4, the dot laid outside the greyish region indicates the empirical distribution, i.e., the temporary estimate
ξ ^ = ( 0.100 , 0.100 , 0.066 , 0.333 , 0.333 , 0.333 , 0.333 , 0.333 ) ,
and the corresponding MLE is
ξ * = ( 0.122 , 0.122 , 0.108 , 0.329 , 0.299 , 0.327 , 0.327 , 0.299 ) .
Furthermore, the greyish region represents the physical domain B cut by a two-dimensional affine subspace of Ξ specified by the equation
ξ = ( 1 s ) ξ ^ + s ξ * + t v .
The vector v was chosen randomly under the condition that
v ξ ^ ξ * and v = ξ ^ ξ * ,
where the orthogonality ⊥ and the norm · are understood relative to the standard Euclidean structure of R 8 . In Figure 4, the vector v was taken to be
v 1 = ( 0.036 , 0.038 , 0.012 , 0.026 , 0.038 , 0.011 , 0.002 , 0.005 )
in the left panel, and
v 2 = ( 0.028 , 0.000 , 0.006 , 0.034 , 0.024 , 0.022 , 0.034 , 0.030 )
in the right panel.
Figure 4 also demonstrates that the sections of the physical domain B show a variety of shapes. Unfortunately, due to this asymmetry of B, we were unable to find a (nontrivial) two-dimensional affine subspace on which every ˜ ( m ) -projection runs. Such a difficulty is in good contrast to the simplest case H C 2 , where the set B is rotationally symmetric and the ˜ ( m ) -projections can be displayed on any two-dimensional section of B that passes through the origin of B as Figure 3.

4.3. When dim H 4

The space H = C 4 is also known to admit a full set of MUBs since dim H = 4 is the second power of the prime number 2; for example [22],
| 1 ( 1 ) , | 2 ( 1 ) , | 3 ( 1 ) , | 4 ( 1 ) = 1 0 0 0 , 0 1 0 0 , 0 0 1 0 , 0 0 0 1 , | 1 ( 2 ) , | 2 ( 2 ) , | 3 ( 2 ) , | 4 ( 2 ) = 1 2 1 1 1 1 , 1 2 1 1 1 1 , 1 2 1 1 1 1 , 1 2 1 1 1 1 , | 1 ( 3 ) , | 2 ( 3 ) , | 3 ( 3 ) , | 4 ( 3 ) = 1 2 1 1 i i , 1 2 1 1 i i , 1 2 1 1 i i , 1 2 1 1 i i , | 1 ( 4 ) , | 2 ( 4 ) , | 3 ( 4 ) , | 4 ( 4 ) = 1 2 1 i i 1 , 1 2 1 i i 1 , 1 2 1 i i 1 , 1 2 1 i i 1 , | 1 ( 5 ) , | 2 ( 5 ) , | 3 ( 5 ) , | 4 ( 5 ) = 1 2 1 i 1 i , 1 2 1 i 1 i , 1 2 1 i 1 i , 1 2 1 i 1 i .
It is straightforward to calculate the parameter representation (1) of a state ρ S ( C 4 ) ; however, the corresponding density matrix is rather complicated, and we omit to display it here.
When H = C 6 , or more generally, when dim H is not a power of a prime, we do not know whether H admits a full set of MUBs. Let us touch upon a situation where a Hilbert space H , if it exists, does not admit a full set of MUBs. In this case, there is no measurement basis M ( a ) that allows a parametrization ξ of the state space S ( H ) having a direct connection to the probability distribution of the outcomes as (2). Such a situation could be comparable to the case when the Gell-Mann matrices [23] are used as the measurement basis for estimating an unknown state on H = C 3 . A state ρ S ( C 3 ) is represented as
ρ = ρ x : = 1 3 I + 3 i = 1 8 x i λ i ,
where λ 1 , , λ 8 are the Gell-Mann matrices, and x = ( x 1 , , x 8 ) is a set of real parameters. The physical domain
B = x R 8 ρ x 0
forms a compact convex subset of the unit ball in R 8 . With the state ρ x , the probability distribution of obtaining the eigenvalues ( 1 , 0 , 1 ) of the observable λ 1 is
1 3 x 1 + x 8 3 , 1 2 x 8 3 , 1 + 3 x 1 + x 8 3 ,
while the probability distribution of obtaining the eigenvalues ( 1 , 0 , 1 ) of the observable λ 2 is
1 3 x 2 + x 8 3 , 1 2 x 8 3 , 1 + 3 x 2 + x 8 3 .
Note that the probability of obtaining the eigenvalue 0 of λ 1 is identical to that of λ 2 . However, in a randomized estimation scheme in which λ i is measured N i times, the frequency of obtaining the eigenvalue 0 of λ 1 would be different from that of λ 2 . Consequently, one cannot assign a consistent temporary estimate x ^ 8 for the parameter x 8 in that case. Put differently, the empirical distribution q ^ N on the extended outcome space Ω does not in general have a coordinate representation (4). Thus, the existence of a full set of MUBs is crucial in our analysis.

5. Concluding Remarks

In the present paper, we explored an information geometrical structure of the randomized quantum state tomography, assuming that the Hilbert space under consideration admits a full set of MUBs. We first introduced a classical statistical model { p ( s , ξ ) } s , ξ on an extended sample space Ω , and found that the probability simplex P ( Ω ) was decomposed into mutually orthogonal dualistic foliation (Proposition 1). We then clarified that this geometrical structure had a statistical importance in estimating the coordinate ξ of an unknown quantum state ρ ( ξ ) under the existence of the nuisance parameter s (Proposition 2). This result gave a generalized insight into the ( m ) -projection interpretation for the MLE in that a similar interpretation was still valid for the randomized quantum state tomography by changing the standard Fisher metric into a deformed one. It also provided us with a new, convenient way of data processing in the actual quantum state tomography that may involve unexpected probe particle loss.
It should be noted that the existence of a full set of MUBs ensures the parametrization (1) of the quantum state space S ( H ) . Such a parametrization is distinctive in that it enables a direct correspondence between the parameter space and the probability simplex, realizing the coordinate representation (4) of the empirical distribution q ^ N . Thus, the use of a full set of MUBs is crucial in our analysis. Nevertheless, it is often the case that the Hilbert space under consideration takes the form H ( C p ) n for p = 2 or 3 because qubits or qutrits are often regarded as building blocks of various quantum protocols. Therefore, the existence of a full set of MUBs would not be too strong a requirement in applications.

Author Contributions

The authors contributed equally to this work.

Funding

The present study was supported by JSPS KAKENHI Grant Numbers JP22340019 and JP17H02861.

Acknowledgments

The authors are grateful to Ryo Okamoto and Shigeki Takeuchi for helpful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MUBsmutually unbiased bases
PVMprojection-valued measure
MLEmaximum likelihood estimate

Appendix A. Information Geometry of MLE

Let P ( Ω ) denote the set of probability distributions on a finite sample space Ω , i.e.,
P ( Ω ) : = p : Ω R p ( ω ) > 0 for all ω Ω , and ω Ω p ( ω ) = 1 .
This set may be identified with the ( | Ω | 1 ) -dimensional (open) simplex, where | Ω | denotes the number of elements in Ω , and thus it is sometimes referred to as the probability simplex on Ω . The set P ( Ω ) is also regarded as a statistical manifold endowed with the dualistic structure ( g , ( e ) , ( m ) ) , where g is the Fisher metric, and ( e ) and ( m ) are the exponential and mixture connections [11,12,13].
Suppose that the state of the physical system at hand belongs to a (closed) subset M of P ( Ω ) , but we do not know which is the true state. We further assume that the probability distributions of M are faithfully parametrized by a finite dimensional parameter θ as
M = { p θ ( ω ) | θ Θ } .
In this case, M is called a parametric model, and our task is to estimate the true value of the parameter θ that specifies the true state. Suppose that, by n independent experiments, we have obtained data ( x 1 , x 2 , , x n ) Ω n . This information is compressed into the empirical distribution, an element of P ( Ω ) defined by
q ^ n ( ω ) : = Number   of   occurrences   of   ω   in   data   ( x 1 , x 2 , , x n ) n = 1 n i = 1 n δ x i ( ω )
for each ω Ω , where δ x i ( ω ) is the Kronecker delta. If q ^ n belongs to the model M , then we have an estimate θ ^ n that satisfies p θ ^ n = q ^ n . However, the empirical distribution q ^ n does not always belong to the model M . When q ^ n M , we need to find an alternative estimate from the data. A canonical method of finding an alternative estimate p θ ^ n M is the maximum likelihood method, in which one seeks the maximizer of the likelihood function
θ p θ ( x 1 ) p θ ( x 2 ) p θ ( x n ) ,
within the domain Θ of the parameter θ , so that
θ ^ n : = arg max θ Θ p θ ( x 1 ) p θ ( x 2 ) p θ ( x n ) .
We can rewrite this relation as follows.
θ ^ n = arg max θ Θ 1 n i = 1 n log p θ ( x i ) = arg max θ Θ ω Ω q ^ n ( ω ) log p θ ( ω ) = arg min θ Θ ω Ω q ^ n ( ω ) log q ^ n ( ω ) log p θ ( ω ) = arg min θ Θ D ( q ^ n p θ ) ,
where
D ( q p ) : = ω Ω q ( ω ) log q ( ω ) p ( ω )
is the Kullback-Leibler divergence from q to p. In other words, the maximum likelihood estimate (MLE) p θ ^ n is the point on M that is “closest” from the empirical distribution q ^ n as measured by the Kullback-Leibler divergence:
p θ ^ n = arg min p M D ( q ^ n p ) .
Due to the generalized Pythagorean theorem, the MLE is geometrically understood as the ( m ) -projection from q ^ n to M or its boundary, as illustrated in Figure A1.
Figure A1. The maximum likelihood estimate p θ ^ n is the minimizer of the function p D ( q ^ n p ) with respect to p M , and is also understood as the ( m ) -projection from the empirical distribution q ^ n to M or its boundary.
Figure A1. The maximum likelihood estimate p θ ^ n is the minimizer of the function p D ( q ^ n p ) with respect to p M , and is also understood as the ( m ) -projection from the empirical distribution q ^ n to M or its boundary.
Entropy 20 00609 g0a1

Appendix B. Parametrization of S ( H )

Suppose that the Hilbert space H C d under consideration admits a full set of MUBs
| α ( a ) α { 1 , , d } , ( a = 1 , , d + 1 ) .
For each a { 1 , , d + 1 } , let
M ( a ) : = M 1 ( a ) , , M d ( a ) = | 1 ( a ) 1 ( a ) | , , | d ( a ) d ( a ) | .
Then, the operators
M α ( a ) I d ( a , α ) { 1 , , d + 1 } × { 1 , , d 1 }
are linearly independent, spanning the space of selfadjoint operators with zero trace. This is easily seen from the orthogonality relation:
Tr M α ( a ) I d M β ( b ) I d = δ a b δ α β 1 d .
Thus, given ρ S ( H ) , the operator ρ ( I / d ) is uniquely expanded as
ρ I d = a = 1 d + 1 α = 1 d 1 x α ( a ) M α ( a ) I d ,
where x α ( a ) are real numbers. We can regard x α ( a ) as a coordinate system of the state space S ( H ) . When d = 2 , this is identical to the Stokes parametrization, up to a factor of 2.
Now, let us change the coordinate system x α ( a ) into ξ α ( a ) as
x α ( a ) = ξ α ( a ) + β = 1 d 1 ξ β ( a ) 1 .
We then arrive at the parametrization (1), i.e.,
ρ = a = 1 d + 1 α = 1 d 1 ξ α ( a ) M α ( a ) + 1 β = 1 d 1 ξ β ( a ) M d ( a ) I .
This parametrization is useful in our analysis because it gives a direct connection to the probability distribution of outcomes of the measurement M ( a ) as
p α ( a ) : = Tr ρ M α ( a ) = ξ α ( a ) , for α = 1 , , d 1 , 1 β = 1 d 1 ξ β ( a ) , for α = d .

References

  1. Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  2. Lehmann, E.L.; Casella, G. Theory of Point Estimation, 2nd ed.; Springer: New York, NY, USA, 1998. [Google Scholar]
  3. Helstrom, C.W. Quantum Detection and Estimation Theory; Academic Press: New York, NY, USA, 1976. [Google Scholar]
  4. Holevo, A.S. Probabilistic and Statistical Aspects of Quantum Theory; North-Holland: Amsterdam, The Netherlands, 1982. [Google Scholar]
  5. Hradil, Z. Quantum-State Estimation. Phys. Rev. A 1997, 55, R1561–R1564. [Google Scholar] [CrossRef]
  6. Banaszek, K.; D’Ariano, G.M.; Paris, M.G.A.; Sacchi, M.F. Maximum-likelihood estimation of the density matrix. Phys. Rev. A 1999, 61, 010304. [Google Scholar] [CrossRef] [Green Version]
  7. Hradil, Z.; Summhammer, J.; Badurek, G.; Rauch, H. Reconstruction of the spin state. Phys. Rev. A 2000, 62, 014101. [Google Scholar] [CrossRef]
  8. James, D.F.V.; Kwiat, P.G.; Munro, W.J.; White, A.G. Measurement of qubits. Phys. Rev. A 2001, 64, 052312. [Google Scholar] [CrossRef] [Green Version]
  9. De Burgh, M.D.; Langford, N.K.; Doherty, A.C.; Gilchrist, A. Choice of measurement sets in qubit tomography. Phys. Rev. A 2008, 78, 052122. [Google Scholar] [CrossRef]
  10. Blune-Kohout, R. Optimal, reliable estimation of quantum states. New J. Phys. 2010, 12, 043034. [Google Scholar] [CrossRef] [Green Version]
  11. Amari, S.-I.; Nagaoka, H. Methods of Information Geometry; Translations of Mathematical Monographs 191; AMS and Oxford: Charles Street, RI, USA, 2000. [Google Scholar]
  12. Amari, S.-I. Differential-Geometrical Methods in Statistics; Lecture Notes in Statistics 28; Springer: Berlin, Germany, 1985. [Google Scholar]
  13. Murray, M.K.; Rice, J.W. Differential Geometry and Statistics; Chapman & Hall: London, UK, 1993. [Google Scholar]
  14. Fujiwara, A.; Yamagata, K. Data processing for qubit state tomography: An information geometric approach. arXiv, 2016; arXiv:1608.07983. [Google Scholar]
  15. Fraïsse, J.M.E.; Braun, D. Quantum channel-estimation with particle loss: GHZ versus W states. Quantum Meas. Quantum Metrol. 2016, 3, 53. [Google Scholar] [CrossRef]
  16. Durt, T.; Englert, B.-G.; Bengtsson, I.; Życzkowski, K. On mutually unbiased bases. Int. J. Quantum Inf. 2010, 8, 535–640. [Google Scholar] [CrossRef]
  17. Yuan, H.; Zhou, Z.; Guo, G. Quantum state tomography via mutually unbiased measurements in driven cavity QED systems. New J. Phys. 2016, 18, 043013. [Google Scholar] [CrossRef] [Green Version]
  18. Wootters, W.K.; Fields, B.D. Optimal state-determination by mutually unbiased measurements. Ann. Phys. 1989, 191, 363–381. [Google Scholar] [CrossRef]
  19. Bengtsson, I. Three ways to look at mutually unbiased bases. arXiv, 2006; arXiv:quant-ph/0610216. [Google Scholar]
  20. Ivonovic, I.D. Geometrical description of quantal state determination. J. Phys. A 1981, 14, 3241–3245. [Google Scholar] [CrossRef]
  21. Yamagata, K. Efficiency of quantum state tomography for qubits. Int. J. Quantum Inform. 2011, 9, 1167–1183. [Google Scholar] [CrossRef]
  22. Klappenecker, A.; Rötteler, M. Constructions of mutually unbiased bases. arXiv, 2003; arXiv:quant-ph/0309120. [Google Scholar]
  23. Gell-Mann, M. Symmetries of baryons and mesons. Phys. Rev. 1962, 125, 1067. [Google Scholar] [CrossRef]
Figure 1. Mutually orthogonal dualistic foliation of P ( Ω ) based on M ( s ) and E ( ξ ) . Each section M ( s ) is affinely isomorphic to the parameter space Ξ . The greyish cylindrical area indicates the subset B = { p ( s , ξ ) | s D , ξ B } of P ( Ω ) . In particular, for each s D , the intersection M ( s ) B is affinely isomorphic to the physical domain B that corresponds to the state space S ( H ) .
Figure 1. Mutually orthogonal dualistic foliation of P ( Ω ) based on M ( s ) and E ( ξ ) . Each section M ( s ) is affinely isomorphic to the parameter space Ξ . The greyish cylindrical area indicates the subset B = { p ( s , ξ ) | s D , ξ B } of P ( Ω ) . In particular, for each s D , the intersection M ( s ) B is affinely isomorphic to the physical domain B that corresponds to the state space S ( H ) .
Entropy 20 00609 g001
Figure 2. The maximum likelihood method in the framework of randomized tomography. Given a temporary estimate ( s ^ , ξ ^ ) with ξ ^ B , we can restrict ourselves to the section M ( s ^ ) as the search space for the MLE p * , and p * = p ( s ^ , ξ * ) is the ( m ) -projection from the empirical distribution p ( s ^ , ξ ^ ) to B on the section M ( s ^ ) .
Figure 2. The maximum likelihood method in the framework of randomized tomography. Given a temporary estimate ( s ^ , ξ ^ ) with ξ ^ B , we can restrict ourselves to the section M ( s ^ ) as the search space for the MLE p * , and p * = p ( s ^ , ξ * ) is the ( m ) -projection from the empirical distribution p ( s ^ , ξ ^ ) to B on the section M ( s ^ ) .
Entropy 20 00609 g002
Figure 3. The trajectories of ˜ ( m ) -projections on the Stokes parameter space when N ( 1 ) : N ( 2 ) = 1 : 1 (left) and N ( 1 ) : N ( 2 ) = 5 : 1 (right). The greyish disk represents the Bloch ball B.
Figure 3. The trajectories of ˜ ( m ) -projections on the Stokes parameter space when N ( 1 ) : N ( 2 ) = 1 : 1 (left) and N ( 1 ) : N ( 2 ) = 5 : 1 (right). The greyish disk represents the Bloch ball B.
Entropy 20 00609 g003
Figure 4. A trajectory of ˜ ( m ) -projection displayed on randomly chosen two-dimensional affine subspaces of Ξ to which both the empirical distribution (marked as a dot) and the MLE belong. The greyish region represents the physical domain B.
Figure 4. A trajectory of ˜ ( m ) -projection displayed on randomly chosen two-dimensional affine subspaces of Ξ to which both the empirical distribution (marked as a dot) and the MLE belong. The greyish region represents the physical domain B.
Entropy 20 00609 g004

Share and Cite

MDPI and ACS Style

Fujiwara, A.; Yamagata, K. Information Geometry of Randomized Quantum State Tomography. Entropy 2018, 20, 609. https://0-doi-org.brum.beds.ac.uk/10.3390/e20080609

AMA Style

Fujiwara A, Yamagata K. Information Geometry of Randomized Quantum State Tomography. Entropy. 2018; 20(8):609. https://0-doi-org.brum.beds.ac.uk/10.3390/e20080609

Chicago/Turabian Style

Fujiwara, Akio, and Koichi Yamagata. 2018. "Information Geometry of Randomized Quantum State Tomography" Entropy 20, no. 8: 609. https://0-doi-org.brum.beds.ac.uk/10.3390/e20080609

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop