arxiv:math-ph/ v1 14 Apr 2006

Transcription

1 Volume of the quantum mechanical state space arxiv:math-ph/643v 4 Apr 6 Attila Andai Department for Mathematical Analysis, Budapest University of Technology and Economics, H-5 Budapest XI. Sztoczek u., Hungary April 6, 6 Abstract The volume of the quantum mechanical state space over n-dimensional real, complex and quaternionic Hilbert-spaces with respect to the canonical Euclidean measure is computed, and explicit formulas are presented for the expected value of the determinant in the general setting too. The case when the state space is endowed with a monotone metric or a pullback metric is considered too, we give formulas to compute the volume of the state space with respect to the given Riemannian metric. We present the volume of the space of qubits with respect to various monotone metrics. It turns out that the volume of the space of qubits can be infinite too. We characterize those monotone metrics which generates infinite volume. Introduction The classical Jeffreys prior is the square root of the determinant of the classical Fisher information matrix, up to a normalization constant. Analogously, the quantum mechanical counterpart of the Jeffreys prior is the square root of the determinant of the quantum Fisher-information matrix. In the quantum mechanical case one can endow the state space with different Riemannian metrics. Some of them are the monotone metrics, which can be labeled by special operator monotone functions. In this general case the Jeffreys prior is still unknown explicitly. This prior was widely examined by Slater if the metric is the Bures metric [6, 7, 8, 9,, ]. In this paper we compute the volume of the state space with respect to the Lebesgue-measure, we give general formulas for the volume for monotone and pull-back metrics too and we characterize those monotone metrics which generate infinite volume of the state of qubits. In the first section we fix the notations for further computations and we mention some elementary lemmas which will be used in the sequel. In the second section we compute the volume of the quantum mechanical state space over n- dimensional real, complex and quaternionic Hilbert-spaces with respect to the canonical Euclidean measure. Before these general computations we compute the volume of the state space over the 3 and 4 dimensional real Hilbert-space keywords: state space, volume, monotone metrics; MSC: 53C, 8Q99 andaia@math.bme.hu

2 to give insight in to the general computational method. Moreover, we present explicitly the expected value of the determinant in the general setting. In the third section we consider the case when the state space is endowed with a monotone metric or a pull-back metric and we give formulas to compute the volume of the state space with respect to the given Riemannian metric. Finally in the fourth section we deal with the qubit case. We present the volume of this space with respect to various monotone metrics. It turns out that the volume of the space of qubits can be infinite too. We characterize those monotone metrics which generates infinite volume. Basic lemmas and notations The quantum mechanical state space consists of real, complex or quaternionic self-adjoint positive matrices with trace. We consider only the set of faithful states with real, complex and quaternionic entries. M R n = {X Mn, R X = X, X >, TrX = } M C n = {X Mn, C X = X, X >, TrX = } M H n = {X Mn, H X = X, X >, TrX = } The following lemmas will be our main tools, we will use them without mentioning it, and we also introduce some notations which will be used in the sequel. The first lemma is about some elementary properties of the gamma function Γ. Lemma. Consider the function Γ, which can be defined for z R + as Γz = t z e t dt. This function has the following properties for every natural number n and real parameter z R +. Γn = n! Γ + z = zγz Γ/ = π n!! n!! Γn + / = π Γn/ = π n n For an n n matrix A we set A i to be the left upper i i submatrix of A, where i =,...,n. The next two Lemmas can be found for example in Bhatiah s book [5]. Lemma. The n n self-adjoint matrix A is positive definite if and only if the inequality deta i > holds for every i =,...,n. Lemma 3. Assume that A is an n n matrix with entries x x R and B is a diagonal matrix with the elements B jj on the main diagonal, then deta + B = detb + x n n B jj. j= j i

3 Lemma 4. Assume that A is an n n self-adjoint, positive definite matrix with entries a ij i,j=,...,n and the vector x consists of the first n elements of the last column, that is x = a,n,..., a n,n. Then for the matrix T = deta n A n we have deta = a nn deta n x, Tx. Proof. Elementary matrix computation, one should expand deta by minors, with respect to the last row. Lemma 5. For parameters a, b R + and t R + the integral equalities hold. G a,b = t x a t x b +a+b Γa + Γb + dx = t Γa + b + x a x b dx = Γb + Γ a+ Γ a + b + 3 Proof. These are consequences of the beta integral x p x q dx = Γp + Γq +. Γp + q + Lemma 6. The surface of a unit sphere in an n dimensional space is F n = nπ n Γ n +. Proof. It follows from the well-known formula for the volume of the sphere in n dimension with radius r since F n = dvnr dr. r= V n r = rn π n Γ n +, When we integrate on a subset of the Euclidean space we always integrate with respect to the usual Lebesgue measure. The Lebesgue measure on R n will be denoted by λ n. Lemma 7. Consider the simplex { } n n = x,...,x n ], [ n x k =, k= then n n k x i dλ n x = Γk + n Γnk +. 3

4 Proof. The integral can be computed as n = n k x i dλ n x a n j= aj n a k i... [ ] n k a i a n da n... da da. Integrating with respect to a n, the integral is Γk + Γk + a k Γk + a k... a k k+ n a a n 3 a n and in general, the i-th integral is Γk + Γik + i Γi + k + ak a k...a ik+i ik+i. k a a n i a n i Thus the result is Γk + Γk + Γk + Γk + Γk + Γk + Γ3k Γ3k + 3 Γ4k + 4 Γk + Γn k + Γnk + = Γk + n Γnk +. Lemma 8. Assume that T is an n n self-adjoint, positive definite matrix and k, ρ R +. Set E R n T, ρ = {x Rn x, Tx < ρ}, T ij R; E C nt, ρ = {x C n x, Tx < ρ}, T ij C; E H n T, ρ = {x Hn x, Tx < ρ}, T ij H; then ρ x, Tx k ρ n dλ n x = +k F n G n,k, dett En RT,ρ ρ x, Tx k dλ n x = ρn+k dett F n G n,k, En C T,ρ ρ x, Tx k dλ 4n x = ρn+k dett F 4n G 4n,k. E4n H T,ρ Proof. We prove for the real case only, the other cases can be proved in the same way. The set EnT, ρ is a n dimensional ellipsoid, so to compute the integral first we transform our canonical basis to a new one, which is parallel to the axes of the ellipsoid, this is an orthogonal transformation so it s Jacobian is. After 4

5 we transform this ellipsis to a unit sphere, the Jacobian of this transformation is n ρ, µ k k= where µ k k=,...,n are the eigenvalues of T. Then we compute the integral in spherical coordinates. The integral with respect to the angles give the surface of the sphere F n r n and the radial part is ρ F n r n n dett ρ ρr k dr. Volume of the state space with respect to the Lebesgue measure Before investigating the general setting we compute the volume of the spaces M R 3 and M R 4. For a matrix with real entries A = a f h f b g h g c a f we set A = a, A =, A f b 3 = A and D denotes the matrix, which contains only the diagonal elements of A, that is D ij = δ ij A ii. The matrix A is in M R 3 if and only if the following set of inequalities hold deta = a > a + b + c = deta = ab f > deta 3 = abc + fgh h b g a f c >. These inequalities can be rewritten as a, b, c, f, T f < b deta, h, g, T h, g < c deta, where T i = deta i A i for i =,. It means that for fixed D and A the parameters h, g are in E R T, c deta. We set V A to be equal of the volume of the parameter space of h, g if D and A is given, that is V A = dg dh = c deta π dett Γ = πc deta E R T,c deta If D and A are fixed then we set V A to be equal of the volume of the parameter space f, g, h: V A = E R T,b deta V A df = 5 ab ab πc ab f df = π abc.

6 Finally the volume of the M R 3 space is V M R π 3 = V A dλ = a ab a bdb da = π 4. A 4 4 real, symmetric matrix with diagonal elements a, a, a 3 and a 4 b is an element of the space M R 4 if and only if deta = a > a > a 3 > a 4 > deta = a deta x, T x > deta 3 = a 3 deta x, T x > deta 4 = a 4 deta 3 x 3, T 3 x 3 >. 4 a k = k= We set T i = deta i A i V A 3 = for i =,, 3. For fixed parameters A 3 and D dλ 3 = a 3/ 4 F G, deta3, E R 3 T3,a4 deta3 where we used the notation of the previous example, in this case V A 3 is the volume the space of those parameters which do not belong to A 3 and D. Now assume that A and D is given, then V A = E R T,a3 deta a = F G, a 3/ 4 3/ 4 F G, deta3 dλ E R T,a3 deta = F F G, G,/ a 3/ 4 a 3/ 3 deta. a3 deta x, T x dλ x If D is given 3/ V A = F F G,G,/a4 a 3/ 3 deta dλ E R T,a deta Since = F F G, G,/ a 3/ 4 a 3/ 3 E R T,a deta = F F F G, G,/ G, a a a 3 a 4 3/. a deta x, T x dλ x a a a 3 a 4 3/ Γ3/ + 4 dλ 3 a = Γ 3 the volume of the 4 dimensional real state space is V M R 4 = F F F G, G,/ G, Γ3/ + 4 Γ = 3π4 8 9!. 6

7 Theorem. For every k N the volume of the state spaces M R k and MR k+ are k V M R k = k! πk k +k k!k i! + k! V M R k+ = π k +k k k! k!k i!. + 3k! Proof. A self-adjoint n n matrix with real entries A is in M R n if and only if i {,...,n} : deta i >, a i >, n a k =, where a i,...,n are the diagonal elements of A. First we assume that diagonal elements, D are given. If A n is fixed, then V A n = dλ n = a n / n F n G n, detan. E R n Tn,an detan If A n is fixed, then V A n = V A n dλ n k= E R n Tn,an detan = a n n F n G n, an deta n x, T n x dλ n x = a n E R n Tn,an detan n a n In general if A n k is fixed, then V A n k = n F n F n 3 G n, G n 3,/ deta n. k a n / n+ i F n i G n i,i / deta n k k, because this equation is correct for k = and by induction V A n k dλ n k E R n k T n k,a n k deta n k = k a n / n+ i F n i G n i,i / E R n k T n k,a n k deta n k = k a n / an k deta n k x, T n k x k dλ n k x n+ i F n i G n i,i / a n / n k F n k G n k,k/ deta n k k+ = V A n k. 7

8 It means that V A = n i= n F i So the volume of the real state space is n n / G n i,i / a n. V M R n = n F i G n i,i / n n n / a n dλ n a. The integral in this equation is and the product is ϕ = n If n = k +, then F i G n i,i / Γ n+ n Γ n +n = π n n 4 ϕ = πk + k k k! k! k+ n! n Γ n+ k Γ n n i Γ i+ Γ n i Γ i +. πi! which can be simplified using the equality ΓiΓi + / = i to i π k +k k kk! ϕ = k! k+ i!. If n = k, then using the same identity for the function Γ we have k k ϕ = π k k 3k k k! i!. k! Theorem. For every n N the volume of the state space M C n is V M C n = π nn n n i!.! Proof. The proof is similar to the real case, except that we have to take into account that the dimension of the parameter space of a matrix element A ij for i j indices is. If A n is fixed, then V A n = dλ n = an n F n 3 G n 3, deta n E C n Tn,an detan 8

9 and in general if A n k is fixed, then V A n k = k a n The volume of the complex state space is n+ i F n ig n i,i deta n k k. n V M C n = F i G n i,i n n n a n dλ n a, where the product is and the integral is π n n n n! n i! n! n n!. Theorem 3. For every n N the volume of the state space M H n is Proof. If A n is fixed, then V M H n = n!πn n n n! n i!. V A n = a n n F 4n 5 G 4n 5, deta n and in general if A n k is fixed, then V A n k = k a n The volume of the quaternionic state space is n+ i F 4n 4iG 4n 4i,i deta n k k. V M C n = n F 4i G 4n 4i,i n n n a n dλ n a, where the product is and the integral is π n n n n! n i! n! n n n!. A slight modification of the previous proofs gives the following Theorem. 9

10 Theorem 4. For every parameter α R + and n N the expected value of the function det α on the state spaces M R n, M C n and M H n with respect to the normalized Lebesgue measures µ R, µ C and µ H are M R n M C n M H n Γ deta α dµ R A = n +n Γ n+ deta α dµ C A = n! n! deta α dµ H A = Γn n Γn Γ n+ + α Γ n +n + nα Γn + α Γn + nα n n Γn + α Γn n + αn Γ i+ + α Γ i+ Γi + α Γi n Γi + α. Γi Proof. The proofs are similar, so we just prove the theorem for the real case only. First we compute the integral with respect to the Lebesgue measure, and we divide the result with the volume of the state space. The method is the same as in the previous theorems, so if A n is given, then V A n = deta n α dλ n E R n Tn,an detan = and the general formula is V A n k = E R n Tn,an detan an deta n x, T n x α dλn x = a n /+α n F n G n,α deta n +α k a n /+α n+ i F n i G n i,i /+α deta n k k +α. The integral of the function det α with respect to the Lebesgue measure is deta α dλ dimm R n A = M R n n n n +α F i G n i,i /+α a n dλ n a. n Dividing it with V M R n and after some simplification we get the formula which is in the Theorem. 3 Volume of the state space endowed with Riemannian metrics To simplify the notations the set of real or complex self-adjoint matrices will be denoted by M n, and the set of real or complex states by M n. Čencov and Morozova [6, ] studied first the monotone metrics on classical statistical manifolds, they proved that such metric is unique, up to normalization. The noncommutative extension of the Čencov Theorem was given by

11 Petz [3]. The stochastic maps are the counterpart of the Markovian maps in this setting. A linear map between matrix spaces T : M n M m is called a stochastic map if it is trace preserving and completely positive. Theorem 5. Consider the family of Riemannian-manifolds M n, g n n N. If for every stochastic map T : M n M m the following monotonicity property holds K TD TX, TX K D X, X D, X M n then there exists an operator monotone function f : R + R with the property fx = xfx, such that g D X, Y = Tr X R n,d fl n,d R n,d R Y n,d, for all n N where L n,d X = DX, R n,d X = XD for all D, X M n. These metrics are considered as the noncommutative generalizations of the Fisher-information. These metrics are called monotone metrics. It means that there exists a bijective mapping between the monotone family of metrics and some operator monotone functions. We use the normalization condition for the function f in the previous theorem f =. n Let D M n and choose a basis of R n such that D = µ j E jj is diagonal, where E jk j,k n is the usual system of matrix units. Let us define the following self-adjoint matrices. j= F jk = E jk + E kj H jk = i E jk i E kj j k n j < k n The set of matrices F ij i j n H ij i<j n form a basis of the tangent space at D for complex matrices and F ij i j n form a basis for real ones. We have for the metric from [] that if i < j n, k < l n : gdh ij, H kl = δ ik δ jl mµ i, µ j gdf ij, F kl = δ ik δ jl mµ i, µ j gdh ij, F kl =, if i < j n, k n : gdh ij, F kk = gdf ij, F kk =, if i n, k n : gdf ii, F kk = δ ik 4mµ i, µ i, where mµ i, µ j = µ j f µi The volume of a Riemannian manifold M, g is defined as detgx dλdimm x. M We use the canonical parametrization for the off-diagonal elements of M n and x,..., x n, x + +x n for the diagonal ones. The corresponding µ j.

12 tangent vectors for the diagonal coordinates are A i = E ii E nn for i =,...,n. Since gda i, A j = gde ii, E jj + gde nn, E nn = δ ij µ i + µ n, the determinant of the metric is detgd = m ij ϕ, if D M R n detgd = i<j n i<j n 4m ij ϕ, if D M C n, where ϕ is a determinant of an n n matrix ϕ = det µ + µ n µ n... µ n µ n µ + µ n... µ n µ n µ n... µ n + µ n = detd. Theorem 6. The volume of the real and complex state space endowed with a Riemannian metric which is generated by an operator monotone function f is V M R n, g f = nn 4 M R n V M C n, g f = nn M C n detd detd i<j n i<j n mµ i D, µ j D dλdimm R n D mµ i D, µ j D dλdimm C n D. We can endow the Riemannian space with a pull-back metric too. Consider the functions h : ], [ R with analytic continuation on a neighborhood of the ], [ interval and h x for every x ], [. We call such functions admissible functions. The space M n will geometrically be considered a Riemannian space R d, g E, where d R = n n+ for real matrices and d C = n for complex ones and g E is the canonical Riemannian metric on M n. That is, at every point D M n for every vectors X, Y M n in the tangent space at D the metric is g E DX, Y = Tr XY. For an admissible function h : ], [ R the pull back-geometry of the space M R n and MC n is the Riemannian geometry g h induced by the map φ h,n : M n M n D fd. This Riemannian space will be denoted by M n, g h. For example if the functions h are p p x if p or log x then we get the α-geometries [7, 9]. n If D M n is diagonal D = µ i E ii, then the metric can be computed see

13 [4] as if i < j n, k < l n : gdh ij, H kl = δ ik δ jl Mµ i, µ j gdf ij, F kl = δ ik δ jl Mµ i, µ j gdh ij, F kl =, if i < j n, k n : gdh ij, F kk = gdf ij, F kk =, if i n, k n : gdf ii, F kk = δ ik 4Mµ i, µ i, where hµ i hµ j if µ Mµ i, µ j = i µ j µ i µ j h µ i if µ i = µ j. Using the previous considerations about the volume of the state space we have the following Theorem. Theorem 7. The volume of the real and complex state space endowed with a pull-back metric g h are n n V M R n = nn 4 h µ j Mµ i D, µ j D dλdimm R n D M R n V M C n = nn M C n j= j i n n h µ j j= j i i<j n Mµ i D, µ j D dλ dimm C n D. i<j n 4 Volume of the state space of qubits In the space of qubits we choose the Stokes parametrization i.e. we write a state D M in the form D = + x y + i z. y + i z x Using these coordinates the spaces M C and M R can be identified with the unit ball in the Euclidean spaces R 3 and R. The metric g f which is generated by an operator monotone function in this coordinate system is 4λ λ mλ g f x, y, z =,λ mλ,λ g f x, y = 4λ λ mλ,λ The volume is an integral on the unit ball, which is in spherical and polar coordinates V M C r = 4π dr, r + rf respectively. V M R = π r r + r f r +r r +r dr,. 3

14 Corollary. The volume of the space M, g f where the metric g f is generated by an operator monotone function f is t V M C = π dt + t tft V M R = π t + t t + t ft dt. Here are some operator monotone functions which generate monotone metrics from [3, 4, 3] and the corresponding volumes. fx : V M C : V MR : + x π π x + x x π 8.98 log x x 4π 4 x + 4ππ 4π xx xlog x x π 4 + xlog x.5 x αx + α + αx + α x + βx + β βx + β β β π β? < β β x γ+ + x γ? < The parameters lie in the interval α ] ] [ [, ], β, and γ, ]. We have some open questions about this list. For every function in this list the complex state space has a greater volume; a natural question is if it is necessary? It seems that if for a function f the volume V M C is finite then for the transpose function f x = x fx the volume is infinity, except for the function fx = x, in this case f = f. Is it true in general? For some functions the difference between the volumes is infinity. What can be the statistical meaning of this phenomenon? The origin of the infinite volume of the space M C can be understood partially by the help of a representation theorem for operator monotone functions. This representation Theorem is due to Löwner [], but we use a modified version from [8]. Theorem 8. The map µ f, defined by fx = x tx + t dµt, for x >, 4

15 establishes a bijection between the class of positive Radon measures on [, ] and the class of operator monotone functions. The function f fulfills the condition fx = xfx for every positive x if and only if for every s [, ] the equality µ[, s] = µ[ s, ] holds. If f is an operator monotone function then it s transpose f is monotone too [5]. That is /fx can also be written in the form tx + t dµt, where µ is a probability measure on [, ] with the symmetric property µ[, s] = µ[ s, ]. Substituting this representation of f into the volume formula for M C we have that if µ is the corresponding symmetric measure for the function f, then the volume of the manifold M C,g f is V = z π arccosz + z z z z dµz. The integrand is continuous, monotonously decreasing and has a series expansion π 4 + π + 9π z 4 z 3 + π z +... near the origin. It s integral with respect to a symmetric probability measure is infinity if and only if z dµz = holds. So the volume of the complex state space of qubits is infinity if the metric is generated by a symmetric probability measure which is concentrated in some sense at the ends of the [, ] interval. If we consider the space of qubits with a pull-back metric, then using the above mentioned techniques we have the following corollary. Corollary. For an admissible function f let us consider the real and complex space M with the pull-back metric. The volume of this space is the following. V M R = π + r r f + f + r r V M C = π f + f Acknowledgements f f + r + r f f r dr r dr This work was supported by Hungarian Scientific Research Fund OTKA contract T46599, TS49835, and EU Network QP-Applications contract number HPRN CT 79. 5

16 References [] S. Amari, Differential-geometrical methods in statistics, Lecture Notes in Statistics 8, Springer-Verlag, Berlin-New York, 985. [] S. Amari, H. Nagaoka, Methods of information geometry, Translations of Mathematical Monographs 9, American Mathematical Society, Providence, R.I.. [3] A. Andai, Monotone Riemannian metrics on density matrices with nonmonotone scalar curvature, J. Math. Phys., 44 3, [4] A. Andai, On the curvature of the quantum state space with pull-back metrics, preprint, arxiv:math-ph/ [5] R. Bhatia, Matrix analysis, Graduate Texts in Mathematics 69, Springer- Verlag, New York, 997. [6] N. N. Čencov, Statistical decision rules and optimal inference, Translations of Mathematical Monographs 53, American Mathematical Society, Providence, R.I. 98. [7] P. Gibilisco, T. Isola, On the monotonicity of scalar curvature in classical and quantum information geometry, J. Math. Phys. 46 5, [8] P. Gibisco, T. Isola, Monotone metrics on statistica manifolds of density matrices by geometry of non-commutative l -space, Infinite Dimensional Anal. Quantum Probab. Related Topics, 9 4. [9] P. Gibilisco, G. Pistone, Analytical and geometrical properties of statistical connections in information geometry, in A. Beghi et al., editors, Mathematical Theory of Networks and Systems, Il Poigrafo, Padova 999, [] K. Löwner, Über monotone Matrixfunktionen, Math. Z , [] P. W. Michor, D. Petz, A. Andai, On the curvature of a certain Riemannian space of matrices, Quantum Probab. Infinite Dimensional Anal. 3, 99. [] E. A. Morozova, N. N. Čencov, Markov invariant geometry on state manifolds, in Current problems in mathematics. Newest results Russian, Akad Nauk SSSR. Vsesoyuz. Inst. Nauchn. i Tekhn. Inform. 36, Moscow 989, 69. [3] D. Petz, Monotone metrics on matrix spaces, Linear Algebra Appl , [4] D. Petz, Geometry of canonical correlation on the state space of a quantum system, J. Math. Phys , [5] C. R. Rao, Information and accuracy attainable in the estimation of statistical parameters, Rev. Bull. Calcutta Math. Soc , [6] P. B. Slater, Hall normalization constants for the Bures volumes of the n-state quantum systems, J. Phys. A 3 999,

17 [7] P. B. Slater, A priori probabilities of separable quantum states, J. Phys. A 3 999, [8] P. B. Slater, Comparative noninformativities of quantum priors based on monotone metrics, Phys. Lett. A , 8. [9] P. B. Slater, The quantum Jeffrey s prior/bures metric volume element for squeezed thermal states and a universal coding conjecture, J. Phys. A 9 996, L6 L65. [] P. B. Slater, Quantum Fisher-Bures information of two-level systems and a three-level extension, J. Phys. A 9 996, L7 L75. [] P. B. Slater, Applications of quantum and classical Fisher information to two-level complex and quaternionic and three-level complex systems, J. Math. Phys , [] R. F. Streater, Classical and quantum info-manifolds, Sūrikaisekikenkyūsho Kōkyūroku 96,