From the Jordan Product to Riemannian Geometries on Classical and Quantum States

The Jordan product on the self-adjoint part of a finite-dimensional C*-algebra A is shown to give rise to Riemannian metric tensors on suitable manifolds of states on A, and the covariant derivative, the geodesics, the Riemann tensor, and the sectional curvature of all these metric tensors are explicitly computed. In particular, it is proved that the Fisher–Rao metric tensor is recovered in the Abelian case, that the Fubini–Study metric tensor is recovered when we consider pure states on the algebra B(H) of linear operators on a finite-dimensional Hilbert space H, and that the Bures–Helstrom metric tensors is recovered when we consider faithful states on B(H). Moreover, an alternative derivation of these Riemannian metric tensors in terms of the GNS construction associated to a state is presented. In the case of pure and faithful states on B(H), this alternative geometrical description clarifies the analogy between the Fubini–Study and the Bures–Helstrom metric tensor.


Introduction
The study of geometrical structures on the space of classical and quantum states is a well-developed and constantly growing subject. In the context of classical probability theory, it is almost impossible to overestimate the impact of the Fisher-Rao metric tensor G FR introduced by Rao in [1] based on the Fisher information matrix introduced in [2] (see also [3][4][5][6][7]). From the theoretical point of view, this metric tensor on the space of probability measures is characterized by a universality property, that is, it is the only Riemannian metric tensor (up to a constant number) that is equivariant with respect to the family of Markov morphisms between probability spaces (see [6,[8][9][10][11] for more detailed discussions and for the proofs of this statement in various contexts). In the quantum context, that is, when probability distributions are "replaced" by density operators on the Hilbert space of the system, the situation becomes more involved because, as Petz showed in [12] finishing the work started by Cencov and Morozowa in [13], there is an infinite number of metric tensors on the manifold of (invertible) density operators satisfying the quantum version of the equivariance property of the Fisher-Rao metric tensor, namely, the equivariance with respect to the class of completely-positive, trace-preserving maps. Furthermore, much has been discovered (and is still being discovered), on the link between these metric tensors and families of quantum relative entropies (see [14][15][16][17][18]). Quite interestingly, but not completely surprising, all these metric tensors reduce to the Fisher-Rao metric tensor if a suitable classical limit is performed.
From a purely theoretical point of view, there is no need to restrict our attention to geometrical structures on quantum states associated with covariant tensor fields as in the case of metric tensors discussed above. Indeed, in [19][20][21][22][23][24][25], the associative product of the algebra B(H) of linear operators on the finite-dimensional Hilbert space H associated with a quantum system has been suitably exploited to define two contravariant tensor fields on the space of self-adjoint operators on H, and these tensor fields have been used to give a geometrical description of the Gorini-Kossakowski-Lindblad-Sudarshan (GKLS) equation describing the dynamical evolution of open quantum systems (see [19,21,22,24,[26][27][28]). These two tensor fields, named Λ and R, are associated with the antisymmetric part (the Lie product) and the symmetric part (the Jordan product) of the associative product in B(H), respectively. In particular, the tensor field Λ turns out to be the Poisson tensor associated with the coadjoint action of the unitary group U (H) on the space of self-adjoint operators on H, when the latter is thought of as the dual space of the Lie algebra u(H). Consequently, it makes sense to study the symplectic foliation associated with Λ, and it turns out that the symplectic leaves of Λ passing through quantum states are the manifolds of isospectral density operators, in particular, the leaf passing through a pure state is diffeomorphic to the complex projective space of the Hilbert space H of the system. From the point of view of open quantum dynamics, the manifolds of isospectral quantum states are not "big enough" in the sense that the dynamical evolution associated with the GKLS equation will cross these manifolds transversally unless it coincides with the unitary evolution of a closed system. Furthermore, from the point of view of quantum information theory, the manifold on which the metric tensors appearing in Petz's classification are defined is the manifold of invertible quantum states, which is the union of all the manifolds of isospectral quantum states with maximal rank. Consequently, the geometry of Λ is not enough to capture all the relevant features of quantum states.
In [29][30][31], it has been shown that the relevant manifolds of quantum states may be described as submanifolds of homogeneous spaces of the group GL(H) of invertible operators on H. Specifically, if ξ denotes a self-adjoint operator on H, the group GL(H) acts on the space of self-adjoint operators according to the map ξ → g ξ g † , where † denotes the adjoint operator on B(H). This action does not preserve the spectrum, and in particular the trace, of ξ unless g is unitary, however, it preserves the positivity of ξ in the sense that g ξ g † is positive semi-definite if ξ is positive semi-definite, and it preserves the rank of ξ. Furthermore, the orbit of GL(H) through a positive semi-definite operator turns out to be a submanifold of the space of self-adjoint operators on H, which is a homogeneous space of GL(H) consisting of all the positive semi-definite operators with the same rank. Then, on each of these orbits, it is possible to select all those elements with unit trace and prove that this set is actually a submanifold of the orbit which, by construction, consists of all the quantum states with fixed rank. In particular, we have two extreme cases: the manifold of pure quantum states consisiting of quantum states with rank equal to 1; the manifold of invertible quantum states consisting of quantum states with maximal rank. The manifold of quantum states with maximal rank thus obtained coincides with the manifold of invertible quantum states on which Petz's metric tensors are defined. Furthermore, it is possible to link the tensor fields Λ and R with the action of GL(H) by showing that the vector fields associated with linear functions by means of Λ and R provide the representation of the Lie algebra of GL(H) integrating to the action of GL(H) described above (see [20,22]).
In this work, we show that the orbits of GL(H) through positive semi-definite operators behave, with respect to R, as the symplectic leaves (e.g., the manifolds of isospectral states) behave for the Poisson tensor Λ. Specifically, we show that R is invertible on each of these orbits and its inverse gives rise to a Riemannian metric tensor the geometry of which we characterize by computing the covariant derivative, the geodesics, the Riemann tensor, and the sectional curvature. Then, we study the Riemannian geometries on the manifolds of quantum states with the same rank arising from the fact that each of these manifolds is a submanifold of a given orbit of GL(H) through positive semi-definite operators.
To be more precise, in Section 2, we review the geometrical aspects of the space of states of a finite-dimensional C * -algebra that we will need in the rest of the work. In Sections 3 and 4, we actually prove the aforementioned results not only in the case of the quantum states associated with the algebra B(H), but in the more general case of states of a finite-dimensional C * -algebra A . This framework allows us to deal with classical and quantum states with the same formalism because the states of a finite-dimensional, Abelian C * -algebra A are in one-to-one correspondence with the probability distributions on a finite-outcome space. Furthermore, it nicely fits into the recently developed groupoidal approach to quantum theories developed in [32][33][34][35][36][37].
Then, in Section 6, we will take inspiration from Uhlmann's geometric construction of the Bures-Helstrom metric tensor (see [38][39][40][41][42][43][44][45][46][47]) to show that the Riemannian geometries on the states associated with the Jordan product may be realized as "projected shadows" of the Riemannian geometries of suitable spheres in suitable Hilbert spaces by means of the GNS construction.
The geometrical picture that will eventually emerge from the present work is that the Jordan product derived from the associative product in A generates Riemann metric tensors on the manifolds of states on A that are associated with the action of the group G of invertible elements in A . In particular, in Sections 5, 7, and 8 respectively, these Riemannian metric tensors are shown to coincide with the Fisher-Rao metric tensor when A is Abelian, with the Fubini-Study metric tensor when A = B(H) and we consider the manifold of pure quantum states (rank-one projectors in B(H)) , and with the Bures-Helstrom metric tensor when A = B(H) and we consider the manifold of invertible quantum states. Consequently, with the help of some imagination, we may interpret all these seemingly different metric tensors as being different faces of the same object, namely, the contravariant tensor R determined by the Jordan product. Section 9 contains some concluding remarks.

Geometrical Aspects of Positive Linear Functionals and States
Let A be a finite-dimensional, unital C * -algebra in which the involution is denoted by † and unit by I. We refer to [48][49][50] for the basic definitions concerning C * -algebras. Let A sa be the self-adjoint part of A . The associative product of A gives rise to a commutative product {, } and to a non-commutative product [ where ı is the imaginary unit. Note that both {, } and [[, ]] are non-associative, and [[, ]] satisfies the Jacobi identity. These two products make A sa into a Banach-Lie-Jordan algebra (see [51][52][53][54][55][56] [57]). Let V be the self-adjoint part of the Banach dual A * of A , that is, the set of all the linear functionals where · denotes complex conjugation.
A positive linear functional ω is called faithful if The space of positive linear functionals on A (excluding the null functional) is denoted by P, while P + denotes the space of faithful, positive linear functionals, which is an open submanifold of V .
For future reference, we need to briefly recall the so-called Gelfand-Naimark-Segal (GNS) construction associated with a positive linear functional ω (see [48,50,58] for more details). Given ω, we define the set This is a left ideal in A called the Gel'fand ideal of ω. Then, we consider the bilinear form on A given by (a, b) ω := ω(a † b).
It is easily seen that (, ) ω is a pre-inner product that descends to the quotient By completing H ω with respect to (, ) ω , we obtain a complex Hilbert space H ω the elements of which are written as ψ a to emphasize that they are associated with (equivalence classes of) elements of A . The Hilbert space product on H ω is written as , , and there is a representation The vector ψ I is cyclic with respect to r ω and separating for the commutant of r ω (A ) in B(H ω ) (see prop. 2.5.3 of [49]). Moreover, every vector in H ω gives rise to a positive linear functional on A by means of If ω is faithful, then N ω = {0} and H ω = A . Therefore, in the finite-dimensional case, H ω coincides with A and r ω becomes the left regular representation a → L a of A on itself.

Proposition 1.
Let ω be a positive linear functional on A , and let a ∈ A be an element of the Gel'fand ideal N ω (see Equation (5)). Then, we have Proof. We just need to apply the Cauchy-Schwarz inequality (see prop. 2.3.10.b of [49]) to the positive-semidefinite sesquilinear form defined by ω in Equation (6). Specifically, we have where, in the last equality, we used the fact that a is in N ω .
The space P is not a smooth manifold in the usual sense of differential geometry (see [59,60] for the appropriate definition of smooth manifold). However, there is a linear left action α of the Lie group G of invertible elements in A on the space V given by (see [29][30][31]) This action preserves P, and every orbit of α is a smooth submanifold of V . That is, even though P itself does not have a smooth structure, it is stratified by the orbits of this action all of which are homogeneous spaces and hence smooth manifolds; in fact, the manifold structure of these homogeneous spaces coincides with that induced by the inclusion into V as the action of G is defined on all of V . The top stratum, i.e., that of maximal dimension, is easily seen to be the space P + of faithful, positive linear functionals.
Recall that every a ∈ A sa may be identified with a real-valued, linear function on V , that we denote by l a , by means of the expression l a (ξ) := ξ(a). (13) Since V is a finite-dimensional Banach space, the map a → l a is an isomorphism between A sa and V * = A * * sa , and thus the differential of the linear functions on V associated with elements in A sa generate the cotangent space T * ξ V at each ξ. Now, given a ∈ A sa , we introduce the vector fields X a and Y a given by For reasons that will be clear later, we call Y a a gradient vector field, and X a a Hamiltonian vector field. Next, we define the vector field A direct computation based on Equation (14) shows that the Lie bracket [, ] between V ab and V cd reads In particular, it follows that the Hamiltonian vector fields {X a } a∈A sa provide an anti-representation of the Lie algebra u of the unitary group U of A . According to what will be proved below, the left action of U generated by the Hamiltonian vector fields is just α restricted to U .
We will now prove that the vector fields V ab are the fundamental vector fields of the action α. Specifically, we consider an element g ∈ G , and write it as g = e 1 2 (a+ıb) (17) with a, b ∈ A sa . This is always possible because A = A sa ⊕ ıA sa , where ı is the imaginary unit, is the Lie algebra of the Lie group G of invertible elements in A . Then, we may consider the smooth curve g(t) = e t 2 (a+ıb) (18) starting at g(0) = I, and compute the fundamental vector field F of the action α associated with g(t). This vector field is defined as the infinitesimal generator of the one-parameter group of diffeomorphisms of V generated by g(t) by means of α (see [59] (p. 331)). Specifically, it is for all ξ ∈ V . Comparing Equation (19) with Equations (15) and (14), we conclude that F = V ab as claimed.
Let us fix an orbit O ⊂ P of G . The tangent space T ω O is thus identified with the subspace of V ∼ = T ω V written as In this work, the symbol ∼ = denotes the identification of two different sets. Then, the cotangent space From the practical point of view, we define the functions for every b ∈ A sa , where i is the canonical immersion of the orbit O ⊂ P in V , and we obtain that the cotangent vector dl + b (ω) at every ω ∈ O is identified with b. Clearly, since the set {dl b (ξ)} b∈A sa is an overcomplete basis for T * ξ V for every ξ ∈ V , we have that the set {dl + b (ω)} b∈A sa is an overcomplete basis for T * ω O for every ω ∈ O. Now, we will pass from positive linear functional, to states. A positive linear functional ρ is called a state if where I is the identity element in A . The space of states S is the convex body in P which is the intersection of P with the affine hyperplane determined as the inverse image of 1 through the linear function l I , with I ∈ A being the identity element. Consequently, if O is an orbit of G in P through ω, we may consider the inverse image of 1 through the function l + I and obtain a smooth manifold, denoted by O 1 , of states as a closed submanifold of O. Clearly, we may do that for every orbit O in P, and thus we obtain a stratification of S into the disjoint union of smooth manifolds, where the top stratum, denoted by S + , is the space of faithful states. Note that some of these manifolds can be degenerate, i.e., points. In particular, this happens when A is Abelian and O 1 contains a pure state (recall that pure states are the extremal points of S ). In the following, whenever we consider a manifold O 1 of states, we will always implicitly assume that O 1 is not a single point.
Concerning pure states, it is worth mentioning that, according to [61], the functional representation of a commutative C * -algebra in terms of complex-valued functions on the space of pure states may be extended to any noncommutative C * -algebra by looking at the space of pure states as a bundle of Kähler manifolds, and using the Kähler metric to define a noncommutative product between complex-valued functions on the pure states.
We will now see how the group G acts on every O 1 making it a homogeneous space. To this scope, we first note that, if ρ is a state sitting inside O ⊂ P, then α(g, ρ) is in general not a state. From the infinitesimal point of view, this is related to the fact that the gradient vector fields on O do not preserve O 1 because Y a (l + I ) in general does not vanish. However, if we set the canonical immersion map given by identification. Furthermore, every Hamiltonian vector field X a is tangent to O 1 , and we denote by X a the vector field on O 1 which is i 1+ -related with X a . For reasons that will be clear later, we call Y a a gradient vector field, and X a a Hamiltonian vector field. Now, we define the vector fields {Υ ab } a,b∈A sa by means of Quite interestingly, a direct computation shows that which means that we also have Comparing Equation (26) with Equation (16), we conclude that the vector fields {Υ ab } a,b∈A sa provide a representation of the Lie algebra g of G which is tangent to O 1 . Furthermore, if we define the map Φ : O 1 −→ O 1 given by (see also [29][30][31]) it is not hard to show that it is a left action of G on O 1 which is transitive (essentially because α is transitive on O). In particular, the space S + of faithful states is an orbit of G . The flow of the vector field Υ ab is just Φ(g(t), ρ), where g(t) is defined as in Equation (18), and thus the fundamental vector fields of Φ are precisely the Υ ab 's. Note that the map Φ is well-defined because the denominator term is always strictly positive since g is an invertible element. However, note that Φ does not preserve the convex structure of S , that is, we have in general.
Let us now fix the orbit O 1 ⊂ S . Defining the function for every a ∈ A sa , it is immediate to check that The tangent space at ρ ∈ O 1 is identified with the subspace of V ∼ = T ρ V written as while the cotangent space From the practical point of view, given ρ ∈ O 1 , just as it happens for the orbit O ⊂ P, we obtain that the cotangent vector de b (ρ) is identified with b. Clearly, since the set {dl a (ξ)} a∈A sa is an overcomplete basis for T * ξ V for every ξ ∈ V , we have that the set {de b (ρ)} b∈A sa is an overcomplete basis for T * ρ O 1 for every ρ ∈ O 1 .

Remark 1.
An identification similar to that given in Equation (31) (with b = 0) may be found also in [5], under the name of e-representation, and in [62] for faithful states. However, in these works, the authors consider only the case A = B(H) (with dim(H) < ∞) so that they identify the space of faithful states S + with the space of faithful density operators in B(H) by means of the isomorphism between B(H) and its dual induced by the trace on H, and no mention is made of the gradient and Hamiltonian vector fields nor of the associated action of G on S + . On the other hand, here we want to stress that the identification of T ρ O 1 with a linear subspace of V given by Equation (31) works for every orbit O 1 ⊂ S and it is part of the "internal geometry" of the space of states of A .

From the Jordan Product to Riemannian Geometries
We will now exploit the Jordan-Lie-algebra structure of A sa introduced above to obtain geometric tensor fields on V , specifically, we obtain a symmetric, contravariant bivector field R associated with the Jordan product {, }, and a Poisson bivector field Λ associated with the Lie product [[, ]] on A sa . This is the generalization to a generic (finite-dimensional) C * -algebra A of what is done in [19][20][21][22][23][24][25] for the specific case A = B(H) for a finite-dimensional Hilbert space H. Then, we will show how the manifolds of positive linear functionals introduced in the previous section may be interpreted as a sort of analogs of symplectic leaves for the symmetric tensor R in a sense that will be specified later. This will allow us to define Riemannian geometries on the orbits of G in P that will be studied in some detail.
In order to define R and Λ, we recall that the differentials of the linear functions l a with a ∈ A sa generate the cotangent space T * ξ V at every ξ ∈ V , so that we may set and extend these objects by linearity obtaining two contravariant tensor fields The antisymmetry of [[, ]] implies that Λ is antisymmetric, while the symmetricity of {, } implies that R is symmetric. Moreover, note that both tensors have non-constant rank, and Λ = 0 if A is Abelian. If A = B(H) for some finite-dimensional Hilbert space H, it is a matter of direct inspection to show that tensor fields R and Λ defined above coincide with those introduced in [19][20][21][22][23][24][25].
The Lie algebra of the unitary group U may be identified with the space A sa of self-adjoint elements in A (see Equation (17)), and thus Λ may be interpreted as the Kostant-Kirillov-Souriau Poisson tensor associated with the coadjoint action of the unitary group U . Since Λ is a Poisson tensor, we may introduce the Hamiltonian vector field X f associated with a smooth function f on V by means of Λ by setting In particular, it is immediate to check that the Hamiltonian vector field associated with the linear function l a (with a ∈ A sa ) coincides with the fundamental vector field X a = V 0a for α introduced in the previous section. Analogously, we may (improperly) define the gradient vector field Y f associated with a smooth function f on V by means of R by setting Again, it is immediate to check that the gradient vector field associated with the linear function l a (with a ∈ A sa ) coincides with the fundamental vector field Y a = V a0 for α introduced in the previous section. This gives an intimate connection between the tensor field Λ and R and the action α of G on V . In particular, the Hamiltonian vector fields X a generate the action of the unitary group U ⊂ G and the orbits of this action, which are embedded, compact submanifolds of V because U is a compact group (only in finite dimensions), are the leaves of the symplectic foliation associated with the Poisson tensor Λ. When A coincides with the C * -algebra B(H) of linear operators on the finite-dimensional, complex Hilbert space H, it is not hard to see that the orbit through ξ ∈ V is in one-to-one correspondence with the set of self-adjoint operators on H that are isospectral with the self-adjoint operatorξ, which is uniquely associated to ξ by means of the isomorphism between A sa and V induced by the trace on H. When A is Abelian, then Λ = 0 and the action of U on V is trivial. The orbits of U are such that the Poisson tensor Λ is invertible on them and thus gives rise to a symplectic form on every orbit. We may try to obtain a similar construction for the tensor field R, that is, we may try to find suitable submanifolds of V on which R is invertible. In a certain sense, we are looking for analogs, for R, of what would be the symplectic leaves of Λ.
Quite interestingly, we will see that every orbit O of positive linear functionals provides an example of such an analog of a symplectic leaf. In particular, we will see that the inverse G of R on O determines a Riemannian metric tensor, and compute its associated covariant derivative, sectional curvature, and Riemann tensor. Then, we will study the Riemannian geometry on the orbit O 1 ⊂ O of states arising from the canonical immersion i 1+ : O 1 −→ O given by the identification map.
Let us fix the orbit O ⊂ P. It is an immersed submanifold of V , consequently, the set {dl + a (ω)} a∈A sa , where l + a = i * l a , is an overcomplete basis for the cotangent space T ω O. Therefore, we may define the symmetric, (0, 2) tensor field R on O by setting and then extend by linearity just as we did for the definition of R. By construction, we have that for all 1-forms θ 1 , θ 2 on V . The tensor R may be thought of as the restriction of R to O.
The contravariant tensor R is symmetric, invertible and positive.

Proof. By definition of R, a cotangent vector dl
because ω is a positive linear functional. Recalling the definition of the GNS ideal N ω of ω given in Equation (5), and recalling that c is self-adjoint, it is immediate to see that the equality in the previous equation holds if and only if dl where the last equality follows from the fact that c is in N ω and from Proposition 1. Then, since the tangent vector v ω was arbitrary, we conclude that the cotangent vector dl + c (ω) must be the zero cotangent vector, and thus it follows that R is positive and invertible on O ⊂ P.
Because of Proposition 2, the covariant tensor is a Riemannian metric tensor on the orbit O ⊂ P. We can immediately compute the gradient vector field W a associated with the function l + a by means of G. In order to characterize W a , it is sufficient to obtain its action on all the functions from which we obtain that W a = Y a (44) with Y a the fundamental vector field V a0 of the action α of G . Consequently, we have and the fact that Y a is the gradient vector field associated with l + a by means of G explains why we already called it a gradient vector field when we defined it in Section 2. Furthermore, since the By directly applying the definition of a gradient vector field, we have Furthermore, we have where L X a denotes the Lie derivative (see [59,60]), and where we used Equations (45) and (16). However, it also holds where we used Equation (45), the first equality in Equation (14), and the fact that the Lie product is a derivation of the Jordan product. Therefore, comparing the previous two equations, and recalling that the Jordan product is symmetric, we obtain for all a ∈ A sa . This means that the metric G is invariant under the action of the unitary group U defined by the Hamiltonian vector fields associated with the functions l + a with a ∈ A sa . Regarding the evaluation of G on Hamiltonian vector fields, we can say the following. First of all, we note that the set {Y b (ω)} b∈A sa is an overcomplete basis for T ω O for every ω ∈ O, therefore, given the Hamiltonian vector field X b , we can express the tangent vector X b (ω) in terms of gradient tangent vectors, that is, we have where B b ω is an element of A sa that depends on both b and ω. Specifically, B b ω is the element in A sa such that for all c ∈ A sa . In general, B b ω is not unique, but it is defined up to an element in the isotropy algebra g ω of ω. However, since the gradient tangent vector Y a (ω) vanishes for every a ∈ A sa ∩ g ω , we conclude that the non-uniqueness of B b ω does not affect Y B b ω (ω). Consequently, we may write Note that, unlike Equations (45) and (46), Equation (52) expresses a pointwise relation. Now, we will compute the covariant derivative associated with the metric G. By applying the Koszul formula (see [60] (thm. 4.3.1)), we obtain where we used Equations (46) and (48) from which it follows that where we used Equation (14). Eventually, we may write the covariant derivative of a gradient vector field with respect to another gradient vector field as At this point, we may compute the Riemann curvature tensor of G as A direct application of Equation (56) gives and thus In order to compute the scalar products in the last line of the previous equation, we start noting that, since ∇ is compatible with G, we have for every Y a , X b , Y c , which implies Therefore, after exploiting Equation (61) to rewrite the scalar products in the last line of Equation (59), and after a simple manipulation involving the Jacobi identity for the Lie product [[, ]] (see [53]), we obtain Now, a tedious but straightforward computation based on the expansion of all the Jordan and Lie products according to Equation (1) shows that Now, the sectional curvature K G of G (see [60] for the definition) reads where and thus we obtain Note that, when A is Abelian, both the sectional curvature K G and the Riemann curvature R G identically vanish for every orbit O.

Riemannian Metrics on Manifolds of States
Since we have the Riemannian metric G on the orbit O ⊂ P, we may consider the orbit O 1 ⊂ O and pull back G to O 1 by means of the immersion map i 1+ obtaining the Riemannian metric tensor In the following, we will see that this metric tensor determined by the Jordan product "becomes" the Fisher-Rao metric tensor, the Fubini-Study metric tensor, or the Bures-Helstrom metric tensor when we make suitable choices for the explicit form of A and O.
Recalling the definition of the vector fields Y a and X a given in Section 2 below Equation (23), we immediately obtain The set {Y b (ρ)} b∈A sa is an overcomplete basis for T ρ O 1 for every ρ ∈ O 1 , and from the first line of Equation (68) and the second relation in Equation (30), we have that which means that Y a is the gradient vector field associated with e a by means of the Riemannian metric G 1 , and this explains why we already called them gradient vector fields in Section 2. By adapting the proof used for G, we may prove that the Hamiltonian vector fields preserve G 1 , that is, we may prove that L X a G 1 = 0 for all a ∈ A sa . We therefore conclude that G 1 is invariant with respect to the action of the unitary group U on O 1 obtained by restricting Φ to U ⊂ G . Now, we will compute the covariant derivative ∇ 1 associated with G 1 . By applying the Koszul formula (see [60] (thm. 4.3.1)) and proceeding as we did in obtaining Equation (56), we obtain  (23) and (71), it immediately follows that while, from Equations (23) and (56), it immediately follows that This means that O 1 is not a totally geodesic submanifold of O. Then, concerning the Riemann curvature tensor R G 1 associated with G 1 , we have the standard formula which becomes where N abcd , the sectional curvature K G 1 of G 1 is easily computed to be and we see that the sectional curvature of every O 1 is constant and equal to 1 4 when A is Abelian.

The Fisher-Rao Metric Tensor
Here, we will study the case in which A is Abelian, and provide a link between the Riemannian metric tensor G 1 and the Fisher-Rao metric tensor G FR . First of all, we may consider the Abelian C * -algebra A = C n (with n < ∞) endowed with component-wise algebraic operations without loss of generality because, up to isomorphism, this is the unique finite-dimensional, Abelian C * -algebra. Then, we take the canonical basis {e j } j=1,...,n in C n , so that every a ∈ A can be written as a = a j e j with a j ∈ C for all j = 1, ..., n. In particular, a is a self-adjoint element if and only if a j is real for all j = 1, ..., n. By considering the dual basis {e j } j=1,...,n of {e j } j=1,...,n , we introduce its associated Cartesian coordinate system {p j } j=1,...,n on V , and it is immediate to check that V may be identified with R n , the cone P may be identified with the positive orthant P n in R n (without the zero), and the space of states S may be identified with the standard n-simplex ∆ n in R n . Concerning the orbits of G , we first fix a reference element ω * ∈ P such that p j (ω * ) = 0 only for a subset J * ⊆ 2 {1,...,n} , and then note that the orbit O of G through ω * ∈ P is given by all those positive linear functionals ω such that p j (ω) = 0 if and only if j ∈ J * . In particular, every orbit O may be identified with the open interior of the positive cone in R m with m = card(J * ) (cardinality of J * ).
Similarly, if we fix a reference state ρ * ∈ S such that p j (ρ * ) = 0 only for a subset J * ⊆ 2 {1,...,n} (i.e., a probability vector supported on J * ), we have that the orbit O 1 of G through ρ * ∈ S is given by all those states ρ such that p j (ρ) = 0 if and only if j ∈ J * . In particular, every orbit O 1 may be identified with the open interior of the m-simplex in R m with m = card(J * ).
The linear function l a associated with a ∈ A sa reads l a = a j p j , and we have that because A is Abelian, and because A is Abelian and its product operation is defined component-wise. This means that the tensor field Λ vanishes, while the tensor field R may be written as It is then clear that the "inverse" G of R on the orbit O may be written as and the pullback tensor on O 1 coincides with the canonical Fisher-Rao metric tensor G FR (see [5] (sec. 2.2)) on O 1 when the latter is considered as the open interior of the m-simplex in R m with m = card(J * ). Quite provocatively, we may say that the Fisher-Rao metric tensor is a built-in feature of the C * -algebraic approach to classical probability theory on finite sample spaces, and that, in this context, the natural object to start with is the contravariant tensor field R from which the Fisher-Rao metric tensor (a covariant tensor field) may be recovered as explained above.

From the Gelfand-Naimark-Segal Construction to Riemannian Geometries
In this section, we will take inspiration from the reduction procedure adopted to define the Fubini-Study metric tensor on the manifold of pure quantum states (see [64,65]), as well as from Uhlmann's purification procedure adopted to define the Bures-Helstrom metric tensor on the manifold of faithful quantum states (see [44][45][46][47]66,67]), to give a more appealing geometrical picture of the metric tensors on the manifolds of states on A introduced in the previous section. Essentially, we will build a geometrical picture that is somehow dual to the one presented before. Indeed, in Section 4, the Riemannian manifold (O 1 , G 1 ) was thought of as the source of a Riemannian immersion into the Riemannian manifold (O, G), while here, the Riemannian manifold (O 1 , G 1 ) will be shown to be the target of a Riemannian submersion from an open submanifold of a suitably big sphere.
In order to develop our ideas, we need to briefly recall some aspects of the geometry of a complex Hilbert space H as a real Kähler manifold (see [64,68,69]). First of all, every complex Hilbert space H may always be thought of as a real, smooth Hilbert manifold (much of what we will say applies also to infinite-dimensional Hilbert spaces, but we will confine the discussion to the finite-dimensional case). Indeed, we may always consider the realification H R as the model space. The realification H R is obtained by restricting linear combinations of elements in H to have only real coefficients, and by defining a real Hilbert product as the real part of the complex Hilbert product on H. In the following, we will write H to denote both the complex Hilbert space and the real, smooth Hilbert manifold modeled on H R . This should not induce confusion since the context will always clarify which is the mathematical aspect of H we are referring to. The tangent space T ψ H at ψ ∈ H may be identified with H itself because H is a vector manifold. Consequently, we may set and so that a direct computation shows that Ω is a symplectic form and that E is a Riemannian metric tensor. Now, let us fix a state ρ. We denote by H the GNS Hilbert space associated with ρ, by ψ I the GNS vector associated with ρ, and with r the GNS representation of A in B(H). Every vector ψ a ∈ H gives rise to a positive linear functional on A given by in particular, if the vector is normalized, the associated positive linear functional is a state.
For every ψ a ∈ H 0 = H − {0} we set so that ψ a is automatically on the unit sphere S 1 in H, and every element in S 1 may be written as ψ a for some ψ a with a / ∈ N ρ . Note that the condition ψ a |ψ a = ρ(a † a) = 0 means that a is in the ideal N ρ , which means that ψ a is the null vector in H. Then, the map π : S 1 −→ S given by is easily seen to be continuous with respect to the topology underlying the standard differential structure of S 1 and the weak* topology on S (which in the finite-dimensional case coincides with the norm topology). Furthermore, the image of S 1 through π consists of all those states ρ ψ a acting as In particular, from this last equation we immediately see that the smooth homogeneous space O 1 containing ρ is in the image of S 1 through π. Now, in the finite-dimensional case, the Hilbert space H is just the quotient A /N ρ , where N ρ is the Gel'fand ideal associated with ρ. There is a natural projection map pr : A −→ H, and this projection map is an open map because H = A /N ρ is the quotient by a group action. Consequently, the space is open in H because it is the image pr(G ) of the open set G ⊂ A , and if we set we immediately conclude that S 1 (G ) is an open submanifold of the unit sphere S 1 , and that, essentially by definition, every ψ ∈ S 1 (G ) is such that there exists an invertible element g ∈ G such that ψ = pr(g) = ψ g . This means that the image π(S 1 (G )) in S through the map π defined in Equation (90) coincides with the orbit O 1 ⊂ S and vice versa. This is analogous to the correspondence between the positive octant of the sphere and the unit simplex, as exploited in [9]. We note that there is a left action β of G on the sphere S 1 given by This action is smooth, and its fundamental vector fields Ψ ab are easily seen to be where we implicitly assumed that g = e 1 2 (a+ıb) with a, b ∈ A sa . Essentially by definition, this action preserves the open submanifold S 1 (G ) of the unit sphere, and it is actually transitive on it as can be checked by direct inspection. Consequently, the vector fields Ψ ab provide an overcomplete basis for the tangent space at each ψ g ∈ S 1 (G ). Proof. To prove the proposition, we will show that T ψ g π(Ψ ab ( ψ g )) = Υ ab (π( ψ g )), (96) where Υ ab is the fundamental vector field of the (transitive) action Φ of G on O 1 (see Equation (24)), which proves that T ψ g π is surjective for every ψ g ∈ S 1 (G ). The proof of Equation (96) is obtained by noting that and then directly computing (a−ıb) ) r(c) r(e t 2 (a+ıb) )|ψ g ψ g |r(e t 2 (a−ıb) ) r(e t 2 (a+ıb) )|ψ g c]] (π( ψ g )) − e a (π( ψ g )) e c (π( ψ g )) .
(98) Remark 2. From this proposition, we conclude that every Ψ ab is π-related with the fundamental vector field Υ ab = Y a + X b of the action Φ of G on O 1 given in Equation (27), and thus the action of G on O 1 may be seen as the projected shadow of the action of G on S 1 (G ). In particular, the same is true for the action of the unitary group U on O 1 . Furthermore, the validity of Equation (96) implies that the kernel of T ψ g π coincides with the isotropy algebra g Φ π( ψ g ) of the action Φ at π( ψ g ). This is consistent with the fact that the equality is equivalent to h = k g (100) with k in the isotropy subgroup G π( ψ g ) of π( ψ g ) ∈ O 1 ⊂ S with respect to Φ, as can easily be checked.
On the unit sphere S 1 there is the action of another relevant Lie group, namely, the Lie group U ⊂ B(H), which consists of unitary elements in the commutant A of r(A ) in B(H). The action of U on S 1 is just the restriction of its natural action on H, and this action is proper because both U and S 1 are compact (in the finite-dimensional case). The fundamental vector fields of this action will be denoted by Ξ b , where B ∈ A ⊂ B(H) is skew-adjoint (i.e., B † = −B), and it is possible to prove that every Ξ b commutes with every Ψ ab . Indeed, the flow Fl Ψ ab t of Ψ ab on ψ c is given by (see [59] for the definition of the flow of a vector field) with g t = e t 2 (a+ıb) (see Equation (94)), while the flow Fl with U t = e tB , and, since U t ∈ A , it immediately follows that the flows of these vector fields commute, i.e., the vector fields themselves commute.

Proposition 4.
The group U acts freely and properly on S 1 (G ), and there is a diffeomorphism between the quotient space M = S 1 (G )/U endowed with the canonical smooth structure coming from the free and proper action and the smooth manifold O 1 endowed with the smooth structure recalled in Section 2.
Proof. First of all, a direct computation shows that for every U ∈ U and every ψ g ∈ S 1 (G ), and thus the action of U preserves S 1 (G ) because S 1 (G ) is the preimage of O 1 through π. Then, we note that the action of U in S 1 (G ) is proper for the group being compact. To show that it is free, note that, by construction of the GNS representation, the vector ψ I is cyclic for r(A ), and it is separating for A because of a standard result (see [49] (prop. 2.5.3)). Consequently, if a ∈ A is such that a (ψ I ) = 0, then a is the zero element in B(H). This means that the isotropy group U ψ I of ψ I with respect to the action of U ⊂ A is the trivial group consisting only of the identity operator on H. Furthermore, the isotropy group U ψ g of every ψ g ∈ S 1 (G ) is the trivial subgroup. Indeed, every and thus ψ g is separating for A (see [49] (prop. 2.5.3)) and we may proceed as before. From this, we have that the quotient space M = S 1 (G )/U is a smooth manifold. The group G acts on M by means of the projection of the action β on S 1 (G ) introduced before because the fundamental vector fields Ψ ab generating β commute with the fundamental vector fields Ξ b generating the action of U giving rise to M. Furthermore, this action is transitive on M because β is transitive on S 1 (G ). We denote this action by Φ, and we have From this it follows that h is in the isotropy group of [ ψ g ] with respect to Φ if and only if (recall that the isotropy subgroup of [ ψ g ] with respect to Φ is the set of all group elements in U leaving that is, if and only if there exists a unitary element U h such that Clearly, this means that and thus h is in the isotropy subgroup G π( ψ g ) of π( ψ g ) ∈ O 1 ⊂ S with respect to Φ (see Remark 2).
Therefore, [ ψ g ] ∈ M and π( ψ g ) ∈ O 1 have the same isotropy group for every ψ g ∈ S 1 (G ), and this implies that there is a diffeomorphism between M endowed with the differential structure and the G -action coming from the quotient procedure on S 1 (G ) and the manifold O 1 ⊂ S endowed with the smooth structure and the G -action recalled in Section 2. Note, however, that this diffeomorphism does not extend to the boundaries where it is only a homeomorphism.
At this point, we consider the metric tensor g on S 1 (G ) which is the pullback of the Euclidean tensor E on H by means of the canonical immersion of S 1 (G ) into H given by the identification map, and we will prove that the projection map π : S 1 (G ) −→ O 1 is a Riemannian submersion between (S 1 (G ), g) and (O 1 , G 1 ). First of all, we note that the vector fields Ξ b generating the action of U on S 1 (G ) span the kernel of π because of Proposition 4. Then, a direct computation shows that for every a ∈ A sa and every B in the Lie algebra of U . Indeed, we have because B is skew-adjoint and commutes with r(a). This means that the linear span of the Ψ a0 's is in the orthogonal complement of the vertical distribution. However, π being a submersion, we have and since the Ψ a0 's are π-related with the gradient vector fields Y a (see Remark 2), and these vector fields provide an overcomplete basis of tangent vectors at each point in O 1 , we conclude that the linear span of the Ψ a0 's generates the whole orthogonal complement of the vertical distribution at each point in S 1 (G ). What is left to prove is that for all a, b ∈ A sa . For this purpose, recalling first Equation (95), then Equation (97), and then Equation (68), we have which proves the validity of Equation (113), which implies that π is a Riemannian submersion between (S 1 (G ), g) and (O 1 , G 1 ) as claimed. The fact that π is a Riemannian submersion implies that every geodesic of G 1 on O 1 may be written as the projection of a geodesic of g on S 1 (G ) having initial tangent vector in the horizontal distribution. Therefore, since S 1 (G ) is an open submanifold of the unit sphere in H, and since g is the restriction of four times the round metric on S 1 , we have that the explicit expression of the geodesic γ of g on S 1 (G ) starting at ψ g with initial (constant) tangent vector φ = 0 is given by where |φ| 2 = φ|φ . The tangent vector φ is horizontal if and only if there exists a ∈ A sa such that φ = Ψ a0 ( ψ g ), and it is different from the null vector if and only if a is not in the Gel'fand ideal generated by ρ ψ g = π( ψ g ). Therefore, assuming a not to be in the Gel'fand ideal generated by ρ ψ g = π( ψ g ), a direct computation shows that the geodesic σ starting at ρ g ∈ O 1 with initial tangent vector Y a (ρ g ) reads σ(t) = cos 2 (N g,a t) ρ g + sin 2 (N g,a t) where we have set ρ g := ρ ψ g , In general, this geodesic "leaves" O 1 remaining in S , and after some time it returns in O 1 essentially because the geodesic γ of which σ is the projection is a great circle on a sphere. A case in which O 1 is geodesically complete is when A is the algebra of linear operators on a finite-dimensional, complex Hilbert space and O 1 is the (compact) manifold of pure states. As we will see below, in this case G 1 corresponds to the Fubini-Study metric tensor. Note that, when A = C n , Equation (116) gives an explicit form for the geodesics of the Fisher-Rao metric tensor.

Remark 3.
Note that the procedure applied here to S 1 (G ) may be adapted in the obvious way to the open submanifold H(G ) of H introduced in Equation (92). The result is that we obtain a Riemannian submersion between (H(G ), E ) (where E is the Euclidean metric tensor on H given in Equation (87)) and (O, G) where O is the orbit of G in the space of positive linear functionals P containing the reference state ρ thought of as an element of P. Accordingly, it is possible to obtain an explicit expression also for the geodesics of G.

The Fubini-Study Metric Tensor
We will now explicitly perform the construction presented above in the case where A is the algebra B(H) of linear operators on the finite-dimensional, complex Hilbert space H, and we consider the orbit O 1 of pure states. In doing this, we will essentially recover the standard construction of the diffeomorphism of the complex projective space CP(H) associated with H with the manifold O 1 of pure states, and we will see that the Riemannian metric G 1 on O 1 ∼ = CP(H) is (a multiple of) the Fubini-Study metric tensor.
First of all, let us recall that the space of pure states on B(H) is the space of rank one projectors on H (this is true only in the finite-dimensional case; if H is infinite-dimensional, then the space of rank one projectors on H is the space of normal pure states), that is, a pure state ρ ψ on B(H) may always be written as for some non-zero vector ψ ∈ H. Let us fix a normalized vector ψ ∈ H, and let us introduce an orthonormal basis {e j } j=1,...,dim(H) in H such that e 1 = ψ. Then, the Gel'fand ideal N ρ ψ of the reference state ρ ψ is easily seen to be given by all those linear operators a on H that can be written as From this, we obtain that the GNS Hilbert space associated with ρ ψ can be identified with H itself, and thus the GNS representation of B(H) on the GNS Hilbert space may be identified with B(H) itself.
Moreover, the GNS representation is irreducible (as it must be because we are considering the GNS representation associated with a pure state) and its commutant is given by the multiples of the identity operator on H.
In the case we are considering, it is immediate to check that the open submanifold H(G ) defined in Equation (92) coincides with H 0 = H − {0}, and we conclude that The unitary group U of the commutant of A = B(H) is just the action of the phase group U(1) consisting of elements of the form e ıθ I, with θ ∈ R. Therefore, in the case at hand, the quotient space S 1 (G )/U ∼ = O 1 appearing in the general construction presented in the previous section is just S 1 /U(1), which is precisely the complex projective space associated with H, and we obtain the well-known diffeomorphism between the manifold O 1 of pure states on B(H) thought of as rank one projectors with the complex projective space CP(H). Under this isomorphism, the action of the unitary group U = U (H) on O 1 coincides with the canonical action of the unitary group on the complex projective space, and this is enough to conclude that the pullback to CP(H) of the Riemannian metric G 1 on O 1 is a multiple of the Fubini-Study metric tensors. Indeed, we know that G 1 is invariant with respect to the action of the unitary group on O 1 (see Section 4), and thus its pullback on CP(H) will be invariant with respect to the canonical action of the unitary group on the complex projective space, and this forces this metric to be a multiple of the Fubini-Study metric tensor, since the latter as a metric of a symmetric space is characterized by that property, see for instance [60].

The Bures-Helstrom Metric Tensor
In this subsection, we will explore the case where A is again the algebra B(H) of linear operators on the finite-dimensional, complex Hilbert space H, but the orbit O 1 is the orbit of faithful states. As will be clear, in this case we obtain Uhlmann's construction of the Bures-Helstrom metric tensor (see [44][45][46][47]). Sometimes, this metric tensor is also called the Bures metric, or the Quantum Fisher Information Matrix.
First of all, recall that the Hilbert space trace Tr(·) gives an isomorphism between A and its dual A * . Essentially, every ξ ∈ A * may be identified with an element in A , denoted again by ξ with an evident abuse of notation, such that ξ(a) = Tr(ξ a) ∀ a ∈ A . (121) In view of the literature on the quantum information of finite-level systems, in the remainder of this subsection, we will always maintain a bipolar attitude and think of ξ as either an element of A or of A * , hoping that no serious confusion arises. Accordingly, the vector space V is also thought of as the space A sa of self-adjoint elements in A , P is also thought of as the space of positive elements in A , S is also thought of as the space of density operators in A . Now, we fix the reference state τ to be the maximally mixed state associated with the element I n , where n = dim(H). Since the reference state is faithful, then its Gel'fand ideal contains only the zero element in the algebra A = B(H), and thus the vector space underlying the GNS Hilbert space is B(H) itself. Therefore, the Hilbert product in the GNS Hilbert space H τ reads and we conclude that the GNS Hilbert space is essentially the Hilbert space of Hilbert-Schmidt operators on H. Then, we easily obtain that the open submanifold H(G ) defined in Equation (92) coincides with G itself, and thus S 1 (G ) = H(G ) ∩ S 1 coincides with the set of all the invertible elements in the algebra satisfying the normalization condition 1 n Tr(g † g) = 1 .
The GNS representation of B(H) coincides with the left action of B(H) on itself, and its commutant coincides with B(H) acting by means of the right action on itself. Consequently, if ρ is a faithful state on B(H), that is, an invertible, positive operator on H with unit trace, we have that all the vectors ψ g ∈ S 1 (G ) such that their projections π( ψ g ) coincide with ρ may be written as g u, where g ∈ G satisfies the normalization condition given in Equation (123), and u is a unitary element in G , that is, a unitary operator on H. Comparing what we have just obtained with the results in [44][45][46][47], it follows that the general construction presented in Section 6 essentially reduces to Uhlmann's purification procedure used to define the Bures-Helstrom metric tensor. However, we will now give a more explicit proof of the equivalence of G 1 with the Bures-Helstrom metric tensor, which will also stress the fact that some of the computational difficulties usually associated with the expression of the Bures-Helstrom metric tensor (as explicitly stated for instance in [39] below equation 9.43) may be attributed to a particular realization of the tangent space at each T ρ O 1 , which is not well-suited for the Bures-Helstrom metric tensor.
To better appreciate this instance, we recall that, in the standard approach to the definition of the Bures-Helstrom metric tensor, the manifold S + of faithful states on H is identified with the manifold of invertible density operators, the tangent space T ρ S + is identified with the affine hyperplane A 0 sa of self-adjoint elements in A = B(H) with vanishing trace. This realization of T ρ S + is clearly different from the one used here in terms of the gradient vector fields, and we will now analyze their relation. The gradient vector fields on V allow to identify a tangent vector at ξ ∈ V with an element in A sa by means of Y a (ξ) = {ξ, a}, where ξ is thought of as an element of A sa . It is worth noting that, when we consider a state ρ, the tangent vector Y a (ρ) provides a geometrical version of the Symmetric Logarithmic Derivative at ρ widely used in quantum estimation theory [42,43,[70][71][72]. On the other hand, since V ∼ = A sa , we may also define a constant vector field associated with every a ∈ A sa by setting If we restrict our considerations to the open submanifold P + ⊂ V of faithful positive linear functionals, we may relate gradient vector fields with constant vector fields as follows. First of all, recall that every ω ∈ P + is identified with an invertible element in A sa , and we may write where L ω : A sa −→ A sa is the linear operator given by L ω (a) := ω a, and R ω : A sa −→ A sa is the linear operator given by R ω (a) := a ω. Clearly, both L ω and R ω are positive, invertible linear operators on A sa because ω is a positive, invertible operator on H, and thus A ω also is an invertible linear operator. It clearly holds and thus at every faithful positive linear functional ω ∈ P + . Furthermore, the gradient vector fields on the manifold S + ⊂ P + of faithful states allow to identify a tangent vector at ρ ∈ S + with an element in A sa by means of Y a (ρ) = {ρ, a} − Tr(ρa) ρ.
However, since S + is the submanifold of P + determined by the inverse image of 1 with respect to the linear function l + I , and since a tangent vector at a point in P + may be identified with a self-adjoint element in A sa , we obtain the identification of T ρ S + with the linear subspace A 0 sa ⊂ A sa consisting of traceless elements mentioned before. Specifically, every constant vector field Z a in Equation (125) is tangent to S + whenever a ∈ A sa is such that Tr(a) = 0. Consequently, there is a vector field Z a on S + to which Z a is i 1+ -related.
Then, it is possible to prove that (see [40,41,46]), given arbitrary tangent vectors a, b ∈ A 0 sa ∼ = T ρ S + , the Bures-Helstrom metric tensor acts as where A −1 ρ is the inverse of A ρ (see Equation (126)), which is the invertible (because ρ is invertible) linear operator on B(H) given by A ρ (b) := 1 2 (ρ b + b ρ). Note that the expression of G BU given in Equation (131) crucially depends on the identification of the tangent space T ρ S + at ρ ∈ S + with the space of traceless, self-adjoint elements in B(H). That is, G BH is expressed in terms of the basis of vector fields Z a on S + that are i 1+ -related with the constant vector fields Z a on P + defined in Equation (125), and, in this way, it becomes necessary to find the explicit form of the operator A −1 ρ at every ρ making the explicit computation of the action of the metric tensor not straightforward. On the other hand, the metric tensor G 1 in Equation (68) is "visualized" in terms of the gradient vector fields, and, with respect to these vector fields, its explicit expression is very easy to compute. We will now show that G BH and G 1 coincide because we have (G 1 (Z a , Z b )) (ρ) = (i * 1+ G) ρ (Z a (ρ), Z b (ρ)) = G ρ T ρ i 1+ (Z a (ρ)), T ρ i 1+ (Z b (ρ)) where we used the first equality in Equation (128) in the fourth equality. From Equation (132), we conclude that G 1 = i * 1+ G is precisely the Bures-Helstrom metric tensor as claimed. This shows that with respect to the gradient vector fields, the explicit expression of the Bures-Helstrom metric is relatively easy to compute. The fact that we no longer have to find the explicit expression of the operator A −1 ρ at every ρ is due to the fact that we work with gradient vector fields, and we believe that this instance should be interpreted as an argument in favor of using the gradient vector because these vector fields better express the geometrical properties of the manifold of quantum states.

Concluding Remarks
In the context of quantum information theory, it is well-known that there is an infinite number of metric tensors on the manifold of faithful quantum states satisfying a property, which is the quantum analog of the monotonicity under Markov maps characterizing the Fisher-Rao metric tensor (up to a constant factor) in the classical case. The act of choosing which one of these metric tensors to use is thus an additional freedom that the quantum framework provides.
In this work, we presented a geometrical description of one of these quantum metric tensors, the so-called Bures-Helstrom metric tensor, which is rooted in the C * -algebraic nature of the space of quantum states. Indeed, the theoretical framework of C * -algebras allows to deal with classical probability distributions and quantum states "at the same time" because both of them can be realized as concrete realizations of the space of states on suitable C * -algebras, and from this point of view, the algebraic structures on the algebra may be exploited to investigate the geometrical properties of the space of states. We focused on the finite-dimensional case, and we studied the geometrical structures associated with the symmetric part (Jordan product) of the associative product of the algebra, and the result is the definition of Riemannian metric tensors on submanifolds of states. In particular, we obtain that the Jordan product determines the Fisher-Rao metric tensor in the classical case, the Fubini-Study metric tensor in the case of pure quantum states, and the Bures-Helstrom metric tensor in the case of faithful quantum states, thus providing a theoretical framework in which all these seemingly different Riemannian metric tensors actually appear as different realizations of the "same" conceptual entity.
Finally, we want to mention that the geometrical picture outlined in this work heavily relies on the Jordan product of A sa . As it is known, any associative algebra A over R gives rise to a Jordan algebra A J with Jordan product given by a b := 1 2 (a · b + b · a), where · is the associative product in A (see [51,52,54]). Consequently, it would be interesting to understand if and how it is possible to build a similar picture for the space of states of a Jordan algebra associated with an associative algebra that is not a C * -algebra. In particular, "natural candidates" would be the many types of Geometric Algebras [73,74].