Next Article in Journal
Revisiting Entanglement within the Bohmian Approach to Quantum Mechanics
Next Article in Special Issue
Relating Vertex and Global Graph Entropy in Randomly Generated Graphs
Previous Article in Journal
Coupled Node Similarity Learning for Community Detection in Attributed Networks
Previous Article in Special Issue
Divergence from, and Convergence to, Uniformity of Probability Density Quantiles
 
 
Correction published on 17 October 2018, see Entropy 2018, 20(10), 796.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantum Statistical Manifolds

Departement Fysica, Universiteit Antwerpen, Universiteitsplein 1, 2610 Wilrijk Antwerpen, Belgium
Entropy 2018, 20(6), 472; https://doi.org/10.3390/e20060472
Submission received: 26 May 2018 / Revised: 15 June 2018 / Accepted: 15 June 2018 / Published: 17 June 2018
(This article belongs to the Special Issue Entropy: From Physics to Information Sciences and Geometry)

Abstract

:
Quantum information geometry studies families of quantum states by means of differential geometry. A new approach is followed with the intention to facilitate the introduction of a more general theory in subsequent work. To this purpose, the emphasis is shifted from a manifold of strictly positive density matrices to a manifold of faithful quantum states on the C*-algebra of bounded linear operators. In addition, ideas from the parameter-free approach to information geometry are adopted. The underlying Hilbert space is assumed to be finite-dimensional. In this way, technicalities are avoided so that strong results are obtained, which one can hope to prove later on in a more general context. Two different atlases are introduced, one in which it is straightforward to show that the quantum states form a Banach manifold, the other which is compatible with the inner product of Bogoliubov and which yields affine coordinates for the exponential connection.

1. Introduction

The basic example of a quantum statistical system starts with a self-adjoint operator H on a finite-dimensional or separable Hilbert space H , with the property that the operator exp ( β H ) is trace-class for all β in an open interval D of the real line R . Then, for β D , the quantum expectation value A β of any bounded operator A in H is given by
A β = 1 Z ( β ) Tr e β H A , w i t h   Z ( β ) = Tr e β H .
Note that the quantum expectation is well-defined because the product of a trace-class operator with a bounded operator is again a trace-class operator. The operator H is the Hamiltonian. It defines a one parameter family of quantum states via Equation (1).
The quantum state (1) is a simple example of a model belonging to the quantum exponential family. In this case, the quantum states form a one-dimensional manifold. The goal of the present work is to search for a quantum exponential family formulated in a parameter-free way, similar to the formulation of Pistone and coworkers [1,2,3,4] in the non-quantum case, and to investigate a further generalization involving deformed exponential functions along the lines set out by the author in [5,6]. An alternative approach to parameter-free quantum information geometry is described in [7]. Which approach eventually will lead to a fully developed theory is hard to predict. Such a theory is expected to affect several domains of research, including Quantum Information Theory, Statistical Physics, in particular the study of phase transitions, and Complexity Theory. For a recent review of Information Geometry applied to complexity, see [8].
Early efforts to use geometric methods in the study of non-commutative information theory include the work by Ingarden and coworkers. See, for instance, [9,10]. The relation with Amari’s Information Geometry [11,12] was studied by Hasegawa [13,14,15]. He introduced an alpha-family of divergences D α ( ρ , σ ) where ρ and σ are any pair of density operators on a finite-dimensional Hilbert space H . The approach relies strongly on the properties of the trace.
The metric on the manifold of density matrices is the scalar product introduced by Bogoliubov and used by Kubo and Mori in the context of linear response theory. The generalization to an inner product for vector states on a von Neumann algebra is given in [16]. See also [17].
A more recent account on quantum information geometry is found in Chapter 7 of [12]. See also [18,19] and Example 3.8 of [20].
The parameter-free approach of Pistone and coworkers was generalized to the quantum context by Grasselli and Streater [21,22,23,24], See also [25]. Both the classical case and the quantum case need a regularizing condition on the allowed density functions, respectively density operators. Under this condition, they form a Banach manifold. Recently, Newton [26,27] proposed an alternative regularization based on a specific choice of a deformed logarithmic function. Part of the arguments in [27] can be transposed to the quantum setting [28].
The structure of the paper is as follows. In the next section, quantum states are labeled with operators belonging to the commutant of the GNS-representation rather than with density matrices. Section 3 describes the plane tangent at the reference state. Next, an atlas is introduced which contains a multitude of charts, one for each element of the manifold. Theorem 4 proves that the manifold is a Banach manifold and that the cross-over maps are linear operators. Section 5 introduces the inner product of Bogoliubov. The metric tensor is calculated. Next, alternative charts are introduced and their relation with the metric tensor is investigated. Section 8 and Section 9 discuss the mixture and the exponential connections. Proposition 4 proves that the alternative charts provide affine coordinates for the exponential connection. Section 10 contains a short presentation of the additional structure provided in quantum information geometry by the existence of modular automorphism groups. The final section discusses the results obtained so far. An Appendix about the GNS-representation and the modular operator is added for convenience of the reader.

2. Representation Theorems

In the present paper, the Hilbert space H is assumed to be finite dimensional. This solves the question of choosing an appropriate topology on the manifold of quantum states. In addition, all operators under consideration are bounded continuous. In fact, they are finite-dimensional matrices. In this way, the technical difficulties of working with unbounded operators are avoided.
A density matrix ρ is a self-adjoint operator with discrete spectrum consisting of non-negative eigenvalues which add up to one. This implies that the trace satisfies Tr ρ = 1 . The operator e β H / Z ( β ) , mentioned in the Introduction, is a density matrix of the kind we have in mind.
Introduce the notation A = B ( H ) for the C -algebra of bounded linear operators on the Hilbert space H . The notion of a quantum state coincides with the notion of a (mathematical) state ω on A . The latter is defined as a linear functional ω : A C , which satisfies the conditions of positivity and of normalization
ω ( A A ) 0 for all   A A , ω ( I ) = 1 ,
where I is the identity operator and A is the adjoint of A. In particular, any state ω belongs to the dual space of A as a Banach space.
The state ω is said to be faithful if ω ( A A ) = 0 implies A = 0 .
The Gelfand–Naimark–Segal (GNS) construction shows that given a state ω on a C -algebra A there exists a *-representation π of A as bounded linear operators on a Hilbert space H ω , together with an element Ω of H ω such that
ω ( A ) = ( π ( A ) Ω , Ω ) for all A A
and π ( A ) is dense in H ω . This representation is unique up to unitary equivalence. This representation is used here to make the transition from a situation where quantum states are described by a density matrix to the more general context of an arbitrary von Neumann algebra A of bounded operators on a separable Hilbert space H , together with a cyclic and separating vector Ω H of norm one.
In the case of the algebra of all N-by-N matrices a simple and explicit realization of the GNS-representation is possible. See the Appendix.
The relation between a density matrix ρ and the corresponding quantum state ω ρ , defined by
ω ρ ( A ) = Tr ρ A for all   A A
is a one-to-one relation. Indeed, if two density matrices ρ 1 and ρ 2 produce the same quantum expectations, then they coincide. Conversely, because the Hilbert space H is finite-dimensional, any quantum state ω determines a density matrix ρ such that ω = ω ρ . The state ω ρ is faithful if and only if the density matrix ρ is strictly positive.
For the sake of completeness, the proof of the following result is reproduced.
Theorem 1.
Let ρ and σ be two strictly positive density matrices operating in a finite-dimensional Hilbert space H . Let A denote the von Neumann algebra of linear operators on H . Let π ρ , H ρ , Ω ρ be the GNS-representation induced by ρ. Then, there exists a unique strictly positive operator X in the commutant π ( A ) such that
Tr σ A = ( π ρ ( A ) X 1 / 2 Ω ρ , X 1 / 2 Ω ρ ) f o r   a l l   A A .
Proof. 
Because H ρ is finite-dimensional and Ω ρ is cyclic and separating one has π ρ ( A ) Ω ρ = π ρ ( A ) Ω ρ = H ρ . Hence, there exists X in π ρ ( A ) such that
X Ω ρ = π ρ ( σ ρ 1 ) Ω ρ .
Then, one has for all A A
( π ρ ( A ) Ω ρ , X Ω ρ ) = ( π ρ ( A ) Ω ρ , π ρ ( σ ρ 1 ) Ω ρ ) = ( π ρ ( ρ 1 σ A ) Ω ρ , Ω ρ ) = Tr ρ ρ 1 σ A = Tr σ A .
In particular, take A = B B to obtain
( X π ρ ( B ) Ω ρ , π ρ ( B ) Ω ρ ) = Tr σ B B 0 ,
with equality if and only if B = 0 . This implies that X is a strictly positive operator.
X Ω ρ is the unique element of H ρ for which Equation (6) holds. Because Ω ρ is cyclic for π ρ , it is separating for the commutant. Hence, X is unique as well. ☐
Introduce the notation B ρ for the real Banach space formed by the self-adjoint elements K of π ρ ( A ) satisfying ( K Ω ρ , Ω ρ ) = 0 .
Theorem 2.
Let H , A and ρ be as in the previous Theorem. Let π ρ , H ρ , Ω ρ the GNS-representation induced by ρ. There is a one-to-one correspondence ξ ρ between faithful states ω on A and elements of B ρ . It satisfies
ω ( A ) = e α ρ ( K ) ( π ρ ( A ) e 1 2 K Ω ρ , e 1 2 K Ω ρ ) f o r   a l l   A A ,
with K = ξ ρ ( ω ) and the function α ρ given by
α ρ ( K ) = log ( e K Ω ρ , Ω ρ ) .
Proof. 
Let ω = ω σ . The previous theorem guarantees the existence of a unique strictly positive operator X in the commutant π ρ ( A ) . This operator X can be exponentiated. Let
K = log X ( log X Ω ρ , Ω ρ ) .
Then, ( K Ω ρ , Ω ρ ) = 0 holds by construction and Equation (8) is satisfied with α ρ ( K ) = ( log X Ω ρ , Ω ρ ) (remember that ( X Ω ρ , Ω ρ ) = 1 ).
Conversely, given K, the r.h.s. of Equation (8) defines a faithful state ω of A . ☐
The map ξ ρ is a chart which makes the manifold M of all faithful quantum states into a Banach manifold. The chart ξ ρ is said to be centered at ρ . It satisfies ξ ρ ( ρ ) = 0 .
All representations π ρ , H ρ , with ρ strictly positive, are unitary equivalent and can be identified. Therefore, in what follows, the index ρ of π ρ is dropped and the Hilbert space in which the representation π works is denoted H π .

3. The Tangent Plane at the Center

Let K = ξ ρ ( σ ) B ρ . Introduce the notation
Ψ K = e 1 2 ( K α ( K ) ) Ω ρ .
One has
d d t | t = 0 α ρ ( t K ) = ( K e t K Ω ρ , Ω ρ ) ( e t K Ω ρ , Ω ρ ) | t = 0 = ( K Ω ρ , Ω ρ ) = 0
and
d d t | t = 0 Ψ t K = d d t | t = 0 e 1 2 ( t K α ( t K ) ) Ω ρ = 1 2 K Ω ρ .
The density matrix σ t , defined by
Tr σ t A = ( π ( A ) Ψ t K , Ψ t K ) , A A
satisfies
d d t | t = 0 Tr σ t A = ( π ( A ) K Ω ρ , Ω ρ ) .
Hence, the linear functional f ρ , K defined by
A A f ρ , K ( A ) = ( π ( A ) Ω ρ , K Ω ρ )
is the derivative of the quantum state ω ρ in the direction ω σ , where the density matrix σ = σ 1 has the property that ξ ρ ( σ ) = K .
One concludes that the tangent plane T ρ M at the point ω ρ M consists of all linear hermitian functionals f ρ , K : A C of the form (10), with K B ρ . The functional f ρ , K belongs to the dual of A . In addition,
| | f ρ , K | | = sup A A f ρ , K ( A ) : | | A | | 1 = sup A A ( π ( A ) K Ω ρ , Ω ρ ) : | | A | | 1 = | | | K | 1 / 2 Ω ρ | | 2 | | | K | 1 / 2 | | 2 = | | K | | .
Hence, K f ρ , K is a bounded linear operator. This is a prerequisite for proving in the next Theorem that this map is the Fréchet derivative of the inverse of the chart ξ ρ . This bounded operator is denoted F ρ in what follows. One has F ρ K = f ρ , K . The inverse operator F ρ 1 satisfies F ρ 1 f ρ , K = K . It is well-defined. Indeed, f ρ , K = f ρ , L implies for all B A
0 = ( π ( B B ) ( K L ) Ω ρ , Ω ρ ) = ( ( K L ) π ( B ) Ω ρ , π ( B ) Ω ρ ) .
Because Ω ρ is a cyclic vector it follows that K = L .
Theorem 3.
The inverse of the map ξ ρ : M B ρ , defined in Theorem 2, is Fréchet-differentiable at ω = ω ρ . The Fréchet derivative is denoted F ρ . It maps K to f ρ , K , where the latter is defined by Expression (10).
Proof. 
Let K = ξ ρ ( ω σ ) . One calculates
| | ω σ ω ρ F ρ K | | = sup A A | ω σ ( A ) ω ρ ( A ) F ρ K ( A ) | : | | A | | 1 = sup A A | ( π ( A ) Ω ρ , [ e K α ( K ) I K ] Ω ρ ) | : | | A | | 1 | | e K α ( K ) I K | | | α ( K ) | + o ( | | K α ( K ) | | ) .
Note that
| α ( K ) | log | | e K | | | | K | |
and
| | K α ( K ) | | 2 | | K | | .
In addition, if | | K | | < 1 then one has
α ρ ( K ) log ( 1 + | | K Ω ρ | | 2 ) | K Ω ρ | | 2 .
This holds because λ ≤ 1 implies exp(λ) ≤ 1 + λ + λ2. One concludes that (11) converges to 0 faster than linearly as | | K | | tends to 0. This proves that F ρ K is the Fréchet derivative of ξ ρ ( ω σ ) ω σ at σ = ρ . ☐

4. The Atlas

Following the approach of Pistone and collaborators [1,3,4,27], we build an atlas of charts ξ ρ , one for each strictly positive density matrix ρ . The compatibility of the different charts requires the study of the cross-over map ξ ρ 1 ( σ ) ξ ρ 2 ( σ ) , where ρ 1 , ρ 2 , σ are arbitrary strictly positive density matrices.
Simplify notations by writing ξ 1 and ξ 2 instead of ξ ρ 1 , respectively ξ ρ 2 . Similarly, write Ω 1 and Ω 2 instead of Ω ρ 1 , respectively Ω ρ 2 , and F 1 , F 2 instead of F ρ 1 , respectively F ρ 2 .
Continuity of the cross-over map follows from the continuity of the exponential and logarithmic functions and from the following result.
Proposition 1.
Fix strictly positive density matrices ρ 1 and ρ 2 . There exists a linear operator Y such that for any strictly positive density matrix σ and corresponding positive operators X 1 , X 2 in the commutant A one has X 2 = Y X 1 Y * .
Proof. 
Using the notations of the Appendix one has
X i = J i ( ρ i - 1 / 2 σ ρ i - 1 / 2 I ) J i * , i = 1 , 2 .
Note that the isometry J depends on the reference state with density matrix ρ . Therefore it carries an index i. The above expression for X i implies that
X 2 = Y X 1 Y * w i t h Y = J 2 ( ρ 2 - 1 / 2 ρ 1 1 / 2 I ) J 1 * .
Theorem 4.
The set M of faithful states on the algebra A of square matrices, together with the atlas of charts ξ ρ , where ξ ρ is defined by Theorem 1, is a Banach manifold. For any pair of strictly positive density matrices ρ 1 and ρ 2 , the cross-over map ξ 2 ξ 1 - 1 is continuous.
Proof. 
The continuity of the map X 1 X 2 follows from the previous Proposition. The continuity of the maps K 1 X 1 and X 2 K 2 follows from the continuity of the exponential and logarithmic functions and the continuity of the function α . ☐

5. The Bogoliubov Inner Product

Umegaki’s divergence/relative entropy D ( σ , τ ) of a pair of strictly positive density matrices σ and τ is defined by [29,30,31]
D ( σ | | τ ) = Tr σ ( log σ log τ ) .
It can be used to define a metric tensor g σ , τ ( ρ ) , as explained below.
In the commutative context, Chentsov proved the uniqueness of the Fisher information matrix as a metric which is invariant under Markov morphisms. See, for instance, Theorem 2.1 of [32]. In the quantum case, the additional requirement of the existence of a dually-flat geometry is needed [21]. The notion of quantum relative entropy comes from Quantum Statistical Physics. In Quantum Information Theory, other quantities are being used as well. Alternatives include the trace distance, the Bures distance and the related fidelity function. See, for instance, Chapter 6 of [20].
Introduce σ s and τ t given by
σ s = 1 Z ( s ) exp log ρ + s ( log σ log ρ ) , τ t = 1 W ( t ) exp log ρ + t ( log τ log ρ ) ,
with
Z ( s ) = Tr exp log ρ + s ( log σ log ρ ) , W ( t ) = Tr exp log ρ + t ( log τ log ρ ) .
Both σ s and τ t are well-defined density matrices. The maps s ω σ s and t ω τ t describe two orbits in M , intersecting at ω ρ : σ 0 = τ 0 = ρ . For further use, note that Z ( 0 ) = W ( 0 ) = Z ( 1 ) = W ( 1 ) = 1 .
The metric tensor g σ , τ ( ρ ) is defined by
g σ , τ ( ρ ) = s t | s = t = 0 D ( σ s | | τ t ) .
With the help of the identity
d d t | t = 0 e H + t A = 0 1 d u e u H A e ( 1 u ) H ,
one obtains
d d t | t = 0 log τ t = log τ log ρ d d t | t = 0 log W ( t ) = log τ log ρ 0 1 d u Tr ρ u ( log τ log ρ ) ρ 1 u = log τ log ρ + D ( ρ | | τ ) ,
so that
g σ , τ ( ρ ) = s t | s = t = 0 Tr σ s log τ t = d d s | s = 0 Tr σ s log τ log ρ + D ( ρ | | τ ) = d d s | s = 0 Tr σ s log τ log ρ = 0 1 d u Tr ρ u ( log σ log ρ ) ρ 1 u log τ log ρ ) d d s | s = 0 Z ( s ) Tr ρ log τ log ρ = 0 1 d u Tr ρ u ( log σ log ρ ) ρ 1 u ( log τ log ρ ) D ( ρ | | σ ) D ( ρ | | τ ) .
This is the inner product of Bogoliubov. Its positivity is shown in the next section. It is straightforward to check that g σ , τ = g τ , σ .

6. Alternative Charts

The inner product (16) is expressed in terms of density matrices rather than tangent vectors. Let us therefore calculate the tangent vector of the orbit ω σ s defined by Definition (13).
Lemma 1.
For each A, self-adjoint element of A such that ω ρ ( A ) = 0 , there exists a unique element K of B ρ such that
0 1 d u π ρ u A ρ u Ω ρ = K Ω ρ .
Proof. 
An operator K in the commutant π ( A ) satisfying Equation (17) exists because π ( A ) Ω ρ = H π . It is unique because Ω ρ is separating for π ( A ) . It satisfies
( K Ω ρ , Ω ρ ) = 0 1 d u ( π ρ u A ρ u Ω ρ , Ω ρ ) = 0 1 d u Tr ρ 1 + u A ρ u = ω ρ ( A ) = 0 .
Finally, for any B in A , one has
( K Ω ρ , π ( B ) Ω ρ ) = ( π ( B ) Ω ρ , K Ω ρ ) = ( π ( B ) Ω ρ , 0 1 d u π ρ u A ρ u Ω ρ ) = 0 1 d u ( π ( ρ u A ρ u B ) Ω ρ , Ω ρ ) = 0 1 d u Tr ρ 1 u A ρ u B = 0 1 d u Tr ρ u A ρ 1 u B = 0 1 d u ( π ( B ρ u A ρ u ) Ω ρ , Ω ρ ) = ( K Ω ρ , π ( B ) Ω ρ ) .
This shows that K = K . One concludes that K belongs to B ρ . ☐
Lemma 2.
There exists a strictly positive operator G ρ on H π , which satisfies
G ρ 0 1 d u π ρ u A ρ u Ω ρ = π ( A ) Ω ρ f o r   a l l A A .
Proof. 
First, consider the operator X defined by
X π ( A ) Ω ρ = 0 1 d u π ρ u A ρ u Ω ρ .
It is well-defined because π ( A ) Ω ρ = 0 implies A = 0 . It is a positive operator. This follows from
( X π ( A ) Ω ρ , π ( A ) Ω ρ ) = 0 1 d u π ρ u A ρ u Ω ρ , π ( A ) Ω ρ = 0 1 d u Tr ρ A ρ u A ρ u = 0 1 d u Tr ρ ( 1 u ) / 2 A ρ u A ρ ( 1 u ) / 2 0 .
The latter expression vanishes if and only if ρ ( 1 u ) / 2 A ρ u A ρ ( 1 u ) / 2 = 0 for almost all u in [ 0 , 1 ] . Because ρ is strictly positive, this can happen only if A = 0 . This shows that the operator X is invertible. Take G ρ equal to the inverse of X to obtain the desired result. ☐
Theorem 5.
Let H , A and ρ be as in the previous theorems. Let π , H π , Ω ρ be the GNS-representation of A induced by ρ. Let G ρ be the positive operator defined by the previous lemma.
i) 
There exists a map χ ρ from the faithful states ω on A into the real Banach space B ρ , formed by the self-adjoint elements K of π ρ ( A ) satisfying ( K Ω ρ , Ω ρ ) = 0 , such that for any strictly positive density matrix σ one has
G ρ χ ρ ( ω σ ) Ω ρ = π ( A ρ , σ ) Ω ρ ,
with A ρ , σ in A given by
A ρ , σ = log σ log ρ + D ( ρ | | σ ) .
ii) 
The map χ ρ is injective.
iii) 
For each strictly positive density matrix σ is
d d s | s = 0 ω σ s = f ρ , K ,
with K = χ ρ ( ω σ ) , where σ s is defined by Definition (13) and f ρ , K is defined by (10).
Proof. 
i)
By the previous lemma, one has π ( A ρ , σ ) Ω ρ = G ρ X Ω ρ with X defined by
X = 0 1 d u π ρ u A ρ , σ ρ u .
Note that A ρ , σ is self-adjoint and satisfies ω ρ ( A ρ , σ ) = 0 . Hence, by Lemma 1, there exists a unique K in B ρ such that X Ω ρ = K Ω ρ . This shows that the map χ ρ which maps ω σ onto this element K of B ρ is well-defined.
ii)
Assume that χ ρ ( ω σ ) = χ ρ ( ω τ ) . This implies A ρ , σ = A ρ , τ and hence
log σ + D ( ρ | | σ ) = log τ + D ( ρ | | τ ) .
The latter implies
σ = τ e D ( ρ | | τ ) D ( ρ | | σ ) .
Because Tr σ = Tr τ , it follows that D ( ρ | | τ ) = D ( ρ | | σ ) and hence σ = τ . This shows that the map χ ρ is injective.
iii)
One has for all B A
d d s | s = 0 Tr σ s B = D ( ρ | | σ ) ω ρ ( B ) + 0 1 d u Tr ρ 1 u ( log σ log ρ ) ρ u B = 0 1 d u Tr ρ 1 u A ρ , σ ρ u B = 0 1 d u Tr ρ 1 u B ρ u A ρ , σ = 0 1 d u π B ρ u A ρ , σ ρ u Ω ρ , Ω ρ = π ( B ) G ρ 1 π ( A ρ , σ ) Ω ρ , Ω ρ = π ( B ) χ ρ ( ω σ ) Ω ρ , Ω ρ = f ρ , K ( B ) ,
with K = χ ρ ( ω σ ) .

7. The Riemannian Metric

Introduce an inner product · , · ρ on T ρ M defined by
f ρ , P , f ρ , Q ρ = G ρ P Ω ρ , Q Ω ρ ) ,
for any P , Q in B ρ . The matrix G ρ , which is introduced in Lemma 2, is strictly positive. Hence, the inner product is positive and non-degenerate. It defines a Riemannian geometry on the manifold M .
Theorem 6.
Let ρ , σ , τ be strictly positive density matrices. Let χ ρ be the chart defined in Theorem 5. Let · , · ρ be the inner product defined on the tangent plane T ρ M by (18). The inner product of Bogoliubov, defined by (15), satisfies
g σ , τ ( ρ ) = f ρ , P , f ρ , Q ρ ,
with P = χ ρ ( σ ) and Q = χ ρ ( τ ) .
Proof. 
From Expression (16) follows
g σ , τ ( ρ ) = 0 1 d u π ρ u ( log τ log ρ ) ρ u ( log σ log ρ ) Ω ρ , Ω ρ D ( ρ | | σ ) D ( ρ | | τ ) = 0 1 d u π ρ u ( A ρ , τ D ( ρ | | τ ) ) ρ u ( A ρ , σ D ( ρ | | σ ) ) Ω ρ , Ω ρ D ( ρ | | σ ) D ( ρ | | τ ) .
Use now that
( π ( A ρ , σ ) Ω ρ , Ω ρ ) = 0 1 d u π ( ρ u A ρ , τ ρ u ) Ω ρ , Ω ρ = 0
to obtain
g σ , τ ( ρ ) = 0 1 d u π ρ u ( A ρ , τ ) ρ u ( A ρ , σ ) Ω ρ , Ω ρ .
This can be written as
g σ , τ ( ρ ) = π A ρ , σ Ω ρ , G ρ 1 π ( A ρ , τ ) Ω ρ = G ρ χ ρ ( σ ) Ω ρ , χ ρ ( τ ) Ω ρ = f ρ , P , f ρ , Q ρ .

8. The Mixture Connection

Consider the situation in which the affine combinations of the form
t ρ t = ( 1 t ) ρ 0 + t ρ 1
are the geodesics of the geometry. Introduce the abbreviation ω t = ω ρ t , with the latter defined by (4). The derivative
d d t ω t = ω 1 ω 0
is a tangent vector, which is constant. This implies a vanishing connection.

9. The Exponential Connection

On the other hand, in the case of the exponential connection, the geodesics t ρ t are such that
log ρ t = ( 1 t ) log ρ 0 + t log ρ 1 ζ ( t ) = log ρ 0 + t H ζ ( t ) ,
where ζ ( t ) is a normalizing function and H is defined by H = log ρ 1 log ρ 0 . Note that
ω t ( H ) = Tr ρ t H = Tr ρ t log ρ 1 log ρ 0 = D ( ρ t | | ρ 0 ) D ( ρ t | | ρ 1 ) .
One has
d d t ρ t = d d ϵ | ϵ = 0 exp ( log ρ t + ϵ H ) d ζ d t ρ t = 0 1 d u ( ρ t ) u H ( ρ t ) 1 u d ζ d t ρ t .
Therefore, the derivative of the quantum state becomes
d d t ω t ( A ) = 0 1 d u Tr ( ρ t ) u H ( ρ t ) 1 u A d ζ d t ω t ( A ) .
Take A = I to find that
d ζ d t = ω t ( H ) .
Proposition 2.
The function ζ ( t ) is convex.
Proof. 
Let A = ρ t u / 2 H ρ t u / 2 . Then, one has
Tr ρ t 1 u H ρ t u H = Tr ρ t A A | Tr ρ t A | 2 = | ω t ( H ) | 2 .
Use this in
d 2 d t 2 ζ ( t ) = d d t Tr ρ t H = 0 1 d u Tr ρ t 1 u H ρ t u H | ω t ( H ) | 2 0 .
Because ζ ( 0 ) = ζ ( 1 ) = 0 , it follows that ζ ( t ) 0 on 0 < t < 1 . From Expression (20) and
d d t ω t ( H ) = d 2 d t 2 ζ ( t ) 0 ,
it follows that the expectation ω t ( H ) increases from ω 0 ( H ) = D ( ρ 0 | | ρ 1 ) 0 to ω 1 ( H ) = D ( ρ 1 | | ρ 0 ) 0 .
Proposition 3.
d ω t / d t is the derivative of ω t in the direction ω, with ω such that ξ t ( ω ) = χ t ( ρ 1 ) χ t ( ρ 0 ) .
Proof. 
From the definitions of χ ρ , G ρ and A σ , τ , it follows
( π ( A ) Ω t , χ t ( ω σ ) Ω t ) = ( π ( A ) Ω t , G t 1 π ( A σ t , σ ) Ω t ) = 0 1 d u π ( A ) Ω t , π ( ρ t u A σ t , σ ρ t u ) Ω t = 0 1 d u π ( ρ t u A σ t , σ ρ t u A ) Ω t , Ω t = D ( ρ t | | σ ) ρ t ( A ) + 0 1 d u π ( A ) Ω t , π ρ t u ( log σ log ρ t ) ρ t u Ω t .
Use this for σ = ρ 1 and for σ = ρ 0 and subtract. This gives
( π ( A ) Ω t , [ χ t ( ρ 1 ) χ t ( ρ 0 ) ] Ω t ) = D ( ρ t | | ρ 1 ) D ( ρ t | | ρ 0 ) ρ t ( A ) + 0 1 d u π ( A ) Ω t , π ρ t u ( log ρ 1 log ρ 0 ) ρ t u Ω t = ω t ( H ) ρ t ( A ) + 0 1 d u π ( A ) Ω t , π ρ t u H ρ t u Ω t = ω t ( H ) ρ t ( A ) + 0 1 d u Tr ρ t 1 u H ρ t u A = d d t ω t ( A ) = f t , K ,
with K = χ t ( ρ 1 ) χ t ( ρ 0 ) . ☐
The following result shows that χ ρ ( ω ) is an affine coordinate in the case of the exponential connection.
Proposition 4.
Let us be given a strictly positive density matrix ρ and a geodesic t ρ t of the form (19). Then,
χ ρ ( ω t ) = ( 1 t ) χ ρ ( ω 0 ) + t χ ρ ( ω 1 ) .
Proof. 
From (iii) of Theorem 5, it follows that for all A A
( π ( A ) Ω ρ , χ ρ ( ω t ) Ω ρ ) = d d s | s = 0 Tr σ s A ,
with σ = ρ t . The latter implies
σ s = 1 Z t ( s ) exp log ρ + s ( log ρ t log ρ ) ,
where
Z t ( s ) = Tr exp log ρ + s ( log ρ t log ρ ) .
Use Expression (19) to obtain
σ s = 1 Z t ( s ) exp ( 1 s ) log ρ + s ( ( 1 t ) log ρ 0 + t log ρ 1 ζ ( t ) ) .
There follows
( π ( A ) Ω ρ , χ ρ ( ω t ) Ω ρ ) = 0 1 d u Tr ρ u log ρ + ( 1 t ) log ρ 0 + t log ρ 1 ζ ( t ) ρ 1 u A d d s | s = 0 Z t ( s ) ω ρ ( A ) = ( 1 t ) 0 1 d u Tr ρ u log ρ + log ρ 0 ζ ( t ) ρ 1 u A + t 0 1 d u Tr ρ u log ρ + log ρ 1 ζ ( t ) ρ 1 u A d d s | s = 0 Z t ( s ) ω ρ ( A ) = ( 1 t ) ( π ( A ) Ω ρ , [ χ ρ ( ω 0 ) Ω ρ ) + t ( π ( A ) Ω ρ , χ ρ ( ω 1 ) ] Ω ρ ) d d s | s = 0 Z 0 ( s ) ω ρ ( A ) + d d s | s = 0 Z 1 ( s ) ω ρ ( A ) d d s | s = 0 Z t ( s ) ω ρ ( A ) = π ( A ) Ω ρ , ( 1 t ) χ ρ ( ω 0 ) + t χ ρ ( ω 1 ) .
Because A A is arbitrary and Ω ρ is cyclic and separating one concludes that Equation (24) holds. ☐

10. Modular Automorphisms

The quantum manifold M carries an additional structure, which is induced by the modular automorphism groups, one for each ρ M . In the commutative case, the automorphisms become trivial.
The Tomita–Takesaki theory [33] associates with each state ω ρ in M a self-adjoint operator Δ ρ on H π , which is called the modular operator. The one-parameter group of unitary operators Δ ρ i t defines a group of inner automorphisms of the algebra A . Indeed, for any A in A , the operator Δ ρ i t π ( A ) Δ ρ i t belongs again to π ( A ) . In particular, it induces a group of transformations of the manifold M by mapping any state ω σ onto the state ω σ , t defined by
ω σ , t ( A ) = Δ ρ i t π ( A ) Δ ρ i t Ω σ , Ω σ .
This group of transformations has ω ρ as a fixpoint because Δ Ω ρ = Ω ρ .
A useful property of the group of modular automorphims is the so-called KMS condition, named after Kubo, Martin and Schwinger. Given two elements A and B of A , the function F ( t ) , defined by
F ( t ) = π ( A ) Δ ρ i t π ( B ) Ω ρ , Ω ρ ,
has an analytic continuation in the complex plane such that
F ( t + i ) = π ( B ) Δ ρ i t π ( A ) Ω ρ , Ω ρ .
This property captures the essence of cyclic permutation under the trace and is helpful in the more general context when manipulating non-commuting pairs of operators.

11. Discussion

This paper reviews known and less known results of quantum information geometry. The Hilbert space is assumed to be finite-dimensional to avoid the technicalities coming with unbounded operators. They give rise to domain problems and require a specific choice of operator norm—see [23,28].
The present point of view differs from the usual one, which starts from the Hilbert space generated by the density matrices. Instead, the GNS-representation is used because it is more suited for later generalizations. The main goal of the present work is precisely to present those results for which one would like to find generalizations in the infinitely-dimensional case.
The manifold M of faithful quantum states can be parameterized in many ways. It is tradition to label each quantum state ω by a corresponding density matrix ρ . Here, the parameter-free approach of Pistone and coworkers [1,3,4,27] is followed. In particular, with each element ω ρ of M is associated a chart centered at ω ρ . Two atlases are introduced. The atlas with the charts ξ ρ , introduced in Section 4, is technically less complicated. It turns M into a Banach manifold. However, it is not linked in a straightforward manner with the Riemannian metric induced by Bogoliubov’s inner product. Therefore, another set of charts, denoted χ ρ , is introduced in Section 6. A link between the charts ξ ρ and χ ρ is found in Proposition 3.
The dually affine connections are shortly mentioned in Section 8 and Section 9. In the case of the exponential connection, the charts χ ρ provide affine coordinates.

Conflicts of Interest

The author declares no conflict of interest.

Appendix. The GNS-Representation of a Matrix Algebra

This appendix is added for convenience of the reader. Its content is well-known.
Choose an orthonormal basis of eigenvectors ψ n , n = 1 , 2 , , N , of the strictly positive density matrix ρ . It can be written as
ρ = n = 1 N p n E n where E n is the orthogonal projection onto C ψ n .
Introduce now the vector Ω ρ in H H defined by
Ω ρ = n = 1 N p n ψ n ψ n .
It satisfies | | Ω ρ | | 2 = n = 1 N p n = 1 . A short calculation shows that
Tr ρ A = ( ( A I ) Ω ρ , Ω ρ ) for all A A .
By assumption, all eigenvalues p n are strictly positive. Therefore, | | A Ω ρ | | = 0 implies A = 0 . This shows that Ω ρ is separating. Let E n , m be the orthogonal matrix which maps ψ m onto ψ n . It belongs to A and satisfies
E n , m Ω ρ = p m ψ n ψ m .
This shows that A Ω ρ equals all of H H . Hence, Ω ρ is a cyclic vector for A I . Because the GNS-representation is unique up to unitary equivalence, one concludes that A A I , together with the Hilbert space H H and the vector Ω ρ is equivalent.
Note that the commutant of A I equals I A .
Introduce now an anti-linear operator J defined by
J ( A I ) Ω ρ = ( I A ) Ω ρ
with A given by
A ψ n = m ( A ψ m , ψ n ) ψ m .
From | | A ψ n | | = | | A ψ n | | , it then follows that J is an isometry. A short calculation shows that for all A , B
( J ( A I ) Ω ρ , J ( B I ) Ω ρ ) = ( ( B I ) Ω ρ , ( A I ) Ω ρ ) .
Next, define an anti-linear operator S by S ( A I ) Ω ρ = ( A I ) Ω ρ . This is the modular conjugation operator. One verifies immediately that the conjugate F of S satisfies F ( I A ) Ω ρ = ( I A ) Ω ρ for all A.
The modular operator Δ by definition equals S S = F F . Let us verify that Δ = ρ ρ 1 . It suffices to show that ρ 1 / 2 ρ 1 / 2 J = S . One calculates
( ρ 1 / 2 ρ 1 / 2 ) J ( A I ) Ω ρ = ( ρ 1 / 2 ρ 1 / 2 ) ( I A ) Ω ρ = n p n ( ρ 1 / 2 ρ 1 / 2 ) ψ n A ψ n = n ψ n ρ 1 / 2 A ψ n = m , n ( A ψ m , ψ n ) ψ n ρ 1 / 2 ψ m = m , n p m ( A ψ m , ψ n ) ψ n ψ m = m p m A ψ m ψ m = ( A I ) Ω ρ = S ( A I ) Ω ρ .
From ( ρ A ) = ρ A , it follows that J ( ρ I ) = ( I ρ ) J and hence
S = ( ρ 1 / 2 ρ 1 / 2 ) J = J ( ρ 1 / 2 ρ 1 / 2 ) .
This implies that Δ = ρ ρ 1 .

References

  1. Pistone, G.; Sempi, C. An infinite-dimensional structure on the space of all the probability measures equivalent to a given one. Ann. Stat. 1995, 23, 1543–1561. [Google Scholar] [CrossRef]
  2. Gibilisco, P.; Pistone, G. Connections on non-parametric statistical manifolds by Orlicz space geometry. Inf. Dim. Anal. Quantum Prob. Rel. Topics 1998, 1, 325–347. [Google Scholar] [CrossRef]
  3. Pistone, G.; Rogantin, M. The exponential statistical manifold: Mean parameters, orthogonality and space transformations. Bernoulli 1999, 5, 721–760. [Google Scholar] [CrossRef]
  4. Pistone, G. Nonparametric information geometry. In Geometric Science of Information; Nielsen, F., Barbaresco, F., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 5–36. [Google Scholar]
  5. Naudts, J. Estimators, escort probabilities, and phi-exponential families in statistical physics. J. Ineq. Pure Appl. Math. 2004, 5, 102. [Google Scholar]
  6. Naudts, J. Generalised exponential families and associated entropy functions. Entropy 2008, 10, 131–149. [Google Scholar] [CrossRef]
  7. Ciaglia, M.F.; di Cosmo, F.; Laudato, M.; Marmo, G.; Mele, F.M.; Ventriglia, F.; Vitale, P. A pedagogical intrinsic approach to relative entropies as potential functions of quantum metrics: The q-z family. Ann. Phys. 2018, 395, 238–274. [Google Scholar] [CrossRef]
  8. Felice, D.; Cafaro, C.; Mancini, S. Information geometric methods for complexity. Chaos 2018, 28, 032101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ingarden, R.S. Information geometry in functional spaces of classical and quantum finite statistical systems. Int. J. Eng. Sc. 1981, 19, 1609–1633. [Google Scholar] [CrossRef]
  10. Ingarden, R.S.; Janyszek, H.; Kossakowski, A.; Kawaguchi, T. Information geometry of quantum statistical systems. Tensor 1982, 37, 105–111. [Google Scholar]
  11. Amari, S. Differential Geometric Methods in Statistics. In Lecture Notes in Statistics 28; Springer: New York, NY, USA, 1985. [Google Scholar]
  12. Amari, S.; Nagaoka, H. Methods of Information Geometry; Oxford University Press: Oxford, UK, 2000. [Google Scholar]
  13. Hasegawa, H. α-divergence of the non-commutative information geometry. Rep. Math. Phys. 1993, 33, 87–93. [Google Scholar] [CrossRef]
  14. Hasegawa, H. Exponential and mixture families in quantum statistics: Dual structure and unbiased parameter estimation. Rep. Math. Phys. 1997, 39, 49–68. [Google Scholar] [CrossRef]
  15. Hasegawa, H.; Petz, D. Non-commutative extension of information geometry II. In Quantum Communication, Computing and Measurement; Plenum Press: New York, NY, USA, 1997; pp. 109–118. [Google Scholar]
  16. Naudts, J.; Verbeure, A.; Weder, R. Linear response theory and the KMS condition. Commun. Math. Phys. 1975, 44, 87–99. [Google Scholar] [CrossRef]
  17. Petz, D.; Toth, G. The Bogoliubov inner product in quantum statistics. Lett. Math. Phys. 1993, 27, 205–216. [Google Scholar] [CrossRef]
  18. Petz, D. Quantum Information Theory and Quantum Statistics; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  19. Jenčová, A. Geometry of quantum states: Dual connections and divergence functions. Rep. Math. Phys. 2001, 47, 121–138. [Google Scholar] [CrossRef]
  20. Petz, D. Quantum Information Theory and Quantum Statistics; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  21. Grasselli, M.R.; Streater, R.F. On the uniqueness of the Chentsov metric in quantum information geometry. Infin. Dim. Anal. Quantum Prob. Rel. Top. 2001, 4, 173–182. [Google Scholar] [CrossRef]
  22. Streater, R.F. Duality in quantum information geometry. Open Syst. Inf. Dyn. 2004, 11, 71–77. [Google Scholar] [CrossRef]
  23. Streater, R.F. Quantum Orlicz spaces in information geometry. Open Syst. Inf. Dyn. 2004, 11, 359–375. [Google Scholar] [CrossRef]
  24. Grasselli, M.R. Dual connections in nonparametric classical information geometry. Ann. Inst. Stat. Math. 2010, 62, 873–896. [Google Scholar] [CrossRef]
  25. Jenčová, A. A construction of a nonparametric quantum information manifold. J. Funct. Anal. 2006, 239, 1–20. [Google Scholar] [CrossRef]
  26. Newton, N.J. An infinite-dimensional statistical manifold modeled on Hilbert space. J. Funct. Anal. 2012, 263, 1661–1681. [Google Scholar] [CrossRef]
  27. Montrucchio, L.; Pistone, G. Deformed exponential bundle: The linear growth case. In Geometric Science of Information, GSI 2017 LNCS Proceedings; Nielsen, F., Barbaresco, F., Eds.; Springer: Berlin/Heidelberg, Germany, 2017; pp. 239–246. [Google Scholar]
  28. Naudts, J. Quantum statistical manifolds: The linear growth case. arXiv 2018, arXiv:1801.07642. [Google Scholar]
  29. Umegaki, H. Conditional expectation in an operator algebra. IV. Entropy and information. Kodai Math. Sem. Rep. 1962, 14, 59–85. [Google Scholar] [CrossRef] [Green Version]
  30. Lindblad, G. Completely positive maps and entropy inequalities. Commun. Math. Phys. 1975, 40, 147–151. [Google Scholar] [CrossRef]
  31. Araki, H. Relative entropies for states of von Neumann algebras. Publ. RIMS Kyoto Univ. 1976, 11, 809–833. [Google Scholar] [CrossRef]
  32. Ay, N.; Jost, J.; Lê, H.V.; Schwachhöfer, L. Information Geometry; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  33. Takesaki, M. Tomita’s Theory of Modular Hilbert Algebras and Its Applications; Springer: Berlin/Heidelberg, Germany, 1970. [Google Scholar]

Share and Cite

MDPI and ACS Style

Naudts, J. Quantum Statistical Manifolds. Entropy 2018, 20, 472. https://doi.org/10.3390/e20060472

AMA Style

Naudts J. Quantum Statistical Manifolds. Entropy. 2018; 20(6):472. https://doi.org/10.3390/e20060472

Chicago/Turabian Style

Naudts, Jan. 2018. "Quantum Statistical Manifolds" Entropy 20, no. 6: 472. https://doi.org/10.3390/e20060472

APA Style

Naudts, J. (2018). Quantum Statistical Manifolds. Entropy, 20(6), 472. https://doi.org/10.3390/e20060472

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