Monotonicity on homogeneous spaces

This paper presents a formulation of the notion of monotonicity on homogeneous spaces. We review the general theory of invariant cone fields on homogeneous spaces and provide a list of examples involving spaces that arise in applications in information engineering and applied mathematics. Invariant cone fields associate a cone with the tangent space at each point in a way that is invariant with respect to the group actions that define the homogeneous space. We argue that invariance of conal structures induces orders that are tractable for use in analysis and propose invariant differential positivity as a natural generalization of monotonicity on such spaces.

remain ordered at all subsequent times as the system evolves with the flow. Monotone flows and their discrete-time analogues, order-preserving maps, play an important role in the theory of dynamical systems and find applications to many biological, physical, chemical, and economic models [33,14]. These systems are closely related to linear dynamical systems with input and output channels, where the monotonicity of a nonnegative input is preserved by the output [43,3,21,2,23]. Recently, this type of input-output preserving system has been further extended to the notion of unimodality [22].
One area of application of monotonicity is the theory of consensus algorithms [35,44,45,27,49,48], where one is interested in designing and analyzing consensus protocols that define the interactions between a collection of agents exchanging information about their relative states via a communication network with the aim of achieving collective behavior. Monotone systems also arise naturally in many areas of biology [13,5,6]. A classical example of a monotone system arising from biology is described by a Kolmogorov model of interacting species where an increase in any population causes an increase in the growth rate of all other populations. Such systems are said to be cooperative [51]. Monotone systems also arise in areas of biology other than population dynamics. For instance, see [52] for an example concerning the dynamics of viral infections. Furthermore, monotone subsystems are often found as components of larger networks due to their robust dynamical stability and predictability of responses to perturbations. The decomposition of networks into monotone subsystems and the study of their interconnections using tools from control theory have also proven to be insightful [4,7,11].
This paper addresses the question of how to define monotonicity on a homogeneous manifold. The notion of order plays a defining role in monotonicity theory. In linear spaces, it is well-known that orders are intimately connected with the theory of pointed solid convex cones. In this paper, a solid convex cone K is said to be pointed if K ∩ −K = {0}. Every such cone K induces a partial order ≤ in a vector space, whereby a ≤ b if b − a lies in K. The simplest example is provided by the positive orthant R n + in R n consisting of vectors with nonnegative entries, which induces the standard vector order based on pairwise comparisons of vector entries. It is natural to generalize this approach by defining a field of cones on a manifold, whereby a cone is associated with the tangent space at each point on the manifold. A conal curve of a cone field is defined as a piecewise smooth curve whose tangent vector lies in the cone at every point along the curve wherever it exists. Cone fields induce the notion of conal orders, whereby a pair of points are said to be ordered if the first point can be joined to the second point by a conal curve. Conal orders locally define partial orders on a manifold. Whether the local partial order can be extended globally depends on the structure of the cone field and the underlying space.
In most applications in applied mathematics and engineering, we are interested in problems that are formulated on spaces with special geometries such as homogeneous spaces. These are manifolds that admit a transitive Lie group action and thus provide a way of systematically generating mathematical structures over the tangent bundle using constructs defined at a single point. This provides a methodology to incorporate the symmetries of the space in any additional structures that are endowed to the space for analysis and design purposes. The classical example of such a construction is that of a homogeneous Riemannian metric, which is entirely determined by the metric at a point. In a similar spirit, we define and characterize invariant cone fields on homogeneous spaces. In doing so we closely review elements of the general theory of homogeneous cone fields as outlined in the important work by Hilgert et al. in [25,24] and Neeb in [41]. We then present a number of examples of homogeneous spaces that arise in a variety of applications in information geometry, computational science and engineering, including Grassmann manifolds and spaces of symmetric positive definite matrices, and consider the existence of invariant cone fields on these spaces. A key theme of the paper is that geometric invariance yields 'tractability' in the analysis of orders and related concepts, which otherwise may appear daunting. In particular, we show that conality of geodesics on globally orderable Riemannian homogeneous spaces can be used to determine order relations between points on such spaces.
Cone fields and conal curves provide a local or differential way of thinking about order relations, which can be viewed as corresponding global concepts. Monotonicity itself is a global concept in the sense that it is classically defined in relation to some partial order. In extending any concept defined on vector spaces to manifolds, it is natural to seek the differential characterization of the property, which in turn will often provide a route for generalization to the nonlinear manifold setting. The local property that is equivalent to monotonicity in R n with respect to a partial order defined by a constant cone field is differential positivity [17]. We propose invariant differential positivity (i.e., differential positivity with respect to an invariant cone field) as a generalization of monotonicity to homogeneous spaces. We will show that invariant differential positivity is indeed equivalent to monotonicity when the cone field induces a global partial order. Furthermore, we discuss how the property remains useful in cases where the order is not a global partial order. Invariant differential positivity can be a powerful analytic tool for the study of monotonicity in a variety of contexts, including the theory of consensus of oscillators [36,40], nonlinear dynamical systems [16,15,17], and matrix monotone functions [32,10,37]. In this paper, we specify what we mean by invariant differential positivity on homogeneous spaces, including with respect to cones of rank k [46,47,19,38], which are generalizations of cones to structures that are closed and invariant under scaling by all real numbers. We also consider the strong implications that invariant differential positivity can have for the asymptotic behavior of dynamical systems.

Homogeneous spaces
A left action of a Lie group G on a manifold M is a smooth map Φ : G × M → M satisfying Φ(e, x) = x and Φ(g 1 g 2 , x) = Φ(g 1 , Φ(g 2 , x)) for all g 1 , g 2 ∈ G, x ∈ M, where e is the identity element in G. Note that for a given g ∈ G, the map x → Φ(g, x) is a diffeomorphism of M. The group G is referred to as a transformation group of the manifold M. We will use Φ(g, x) and g · x interchangeably in this paper. A homogeneous space is defined as a manifold M on which a Lie group G acts transitively.
Definition 1 A smooth manifold M is said to be a homogeneous space if there exists a Lie group G acting on M such that for all x 1 , x 2 ∈ M, there exists g ∈ G such that g · x 1 = x 2 .
Homogeneous spaces are closely connected to coset manifolds. For a given Lie group G and a closed subgroup H, consider the set G/H := {gH : g ∈ G} of left cosets of H in G. The set G/H is the set of equivalence classes for the equivalence relation ∼ on G defined by (1) The set G/H can be made into a manifold in a unique way if we require that the projection map π : G → G/H, π(g) := gH be a submersion; i.e., if we require that the differential map dπ| g is surjective for each g ∈ G. For each a ∈ G, define the left translation τ a : G/H → G/H by τ a (gH) := agH. Note that the left translations τ g are related to the left translations L g on the Lie group G by π•L g = τ g •π, for each g ∈ G.

Reductive homogeneous spaces
Let M = G/H be a homogeneous space and consider the natural projection π : G → G/H, π(g) = gH. The differential dπ| e : g → T o (G/H), where o = π(e) = eH, is given by for X ∈ g. As the map dπ| e : T e G → T o M is a vector space homomorphism, we have T e G/(ker dπ| e ) ∼ = im dπ| e . It follows from (2) that ker dπ| e = h, where h is the Lie algebra of H. Thus, we have the canonical isomorphism where g/h is the set of cosets X + h = {X + Y |Y ∈ h} for X ∈ g.
where m := h ⊥ with respect to an Ad-invariant inner product on g. Moreover, we note that the Killing form B : g × g → R, B(X, Y ) = tr(adX • adY ) of a Lie group G is always an Adinvariant symmetric bilinear form on g. Thus, if G is compact and semisimple so that −B is positive definite, then the Killing form defines a bi-invariant metric on G given by ·, · = −B(·, ·) [8].

Symmetric spaces
Symmetric spaces constitute an important class of homogeneous spaces that includes many of the spaces that are of interest in applications and discussed in this paper. A connected Riemannian manifold M is said to be a symmetric space if for each p ∈ M, there exists an isometry j p : M → M, such that where Id p is the identity map on T p M. The map j p has the property that it "reverses" the geodesics that pass through p ∈ M, in the sense that if . The Euclidean space R n is clearly symmetric. A less trivial example is the n-sphere S n embedded in R n+1 , where the symmetry at the north pole p = (1, 0, · · ·, 0) is given by j p (x 1 , x 2 , . . . , Denote the symmetry of the symmetric space M = G/H at o = eH by j. Now for each g ∈ G, define the map σ(g) : M → M by σ(g) = j • g • j. Since σ(g) is an isometry of M, it lies in G and thus we can define an automorphism σ : G → G by g → σ(g) = j • g • j −1 as j 2 = Id. Setting G σ = {g ∈ G : σ(g) = g} to be the fixed points of σ and G o σ its connected component, one can show that σ 2 = Id G and G σ is a closed subgroup of G that satisfies G o σ ⊆ H ⊆ G σ . The map σ is sometimes referred to as the involution map associated with the symmetric space G/H. Every symmetric space G/H with involution σ is a reductive homogeneous space with reductive decomposition g = h ⊕ m, where h = {X ∈ g : dσ| e X = X} and m = {X ∈ g : dσ| e X = −X}, and Ad(H)m ⊆ m [8].
3 Invariant cone fields on homogeneous spaces

Homogeneous cone fields
A wedge is a closed and convex subset of a vector space that is closed under multiplication by nonnegative scalars [25]. A wedge field W M on a manifold M smoothly assigns to each point x ∈ M a wedge W M (x) in the tangent space T x M.
for all g ∈ G and x ∈ M.
We now specialize to the case where the group action is transitive so that M is a homogeneous space. A homogeneous cone field on a homogeneous space M = G/H of a connected Lie group G assigns to each point x ∈ M a cone K M (x) in the tangent space T x M, such that the cone field is invariant under the action of G on M. Recall that elements of M can be identified with cosets gH := {gh : h ∈ H} in G/H and let o := eH denote the base-point in M. The left translations τ g : M → M are defined by τ g (x) = g · x for all g ∈ G and x ∈ M. Let g = T e G and h = T e H denote the Lie algebras of G and H, respectively. The canonical projection π : G → M defined by π(g) = gH induces a linear surjection dπ| e : g → T o M with ker dπ| e = h, so that we obtain the isomorphism g/h ∼ = T o M given by Thus, the surjection dπ| e is identified with the quotient map p : g → g/h, where p(X) = X + h.
for all h ∈ H and X ∈ g.
If we seek to describe invariant cone fields on M = G/H in terms of wedges defined in the Lie algebra g of the total space G, then there are some consistency requirements that must be satisfied. In particular, any cone K in then the associated pointed cone is well-defined. That is, for all g, g ∈ G corresponding to the same point x ∈ M (i.e. for all g, g ∈ G satisfying π(g) = π(g )), we have To see this, note that π(g) = π(g ) precisely if there exists h ∈ H such that g = gh. Thus, we have , whence the result follows from the η(H)-invariance of K. The following theorem from [25] describes the geometry of homogeneous cone fields on M = G/H. See figure 1.
Theorem 1 Let H be a closed subgroup of a Lie group G and W a wedge in where M := G/H, e is the identity element in G, o = eH is the base-point in M, and K is the pointed cone in T o M obtained as the projection of W onto g/h. Then, K G is an invariant wedge field on G and K M is a well-defined homogeneous or G-invariant cone field on M. Moreover, for each g ∈ G, where π : G → M is the canonical projection π(g) = gH.

Lie groups
Any Lie group G is itself a homogeneous space in at least two ways. First, it can be expressed as G ∼ = G/{e}. Alternatively, one can write G ∼ = G × G/G, where G×G acts on G by left and right translations and the isotropy subgroup is G diagonally embedded in G × G. Invariant cone fields can be defined on a Lie group as a homogeneous space using left-translation. Given a cone K in g, the corresponding left-invariant cone field K G is given by for all g ∈ G.

The quotient of the Heisenberg group by its center, G/Z(G)
The Heisenberg group G is a Lie group that arises in various fields including representation theory, sub-Riemannian geometry and quantum mechanics. It can be defined as the group of 3 × 3 upper triangular matrices with diagonal elements equal to 1 and group operation given by matrix multiplication. The Lie algebra g can be represented as the set of stricly upper triangular 3 × 3 matrices. That is, The center Z(G) of G is defined as the set Z(G) = {z ∈ G : zg = gz, ∀g ∈ G} and forms a subgroup H of G with Lie algebra h given by The quotient manifold G/H defines a homogeneous space of dimension 2. To construct an invariant cone field on G/H, we first look for a wedge W in g that satisfies condition (i) of Theorem 1; i.e., W ∩ −W = h. This is achieved precisely if W is of the form where K is any pointed convex solid cone in R 2 . For condition (ii) of Theorem 1, we consider Ad(H)W . Now since  for all c ∈ R and (α, β, γ) ∈ R 3 , Ad(H)(W ) = W trivially holds for any wedge W in g. Therefore, any wedge W of the form (17) uniquely defines an invariant cone field K M on M = G/H.

The n-spheres and Grassmannians
The n-sphere S n can be viewed as a homogeneous space with a transitive SO(n + 1) action, since any two points on S n embedded in R n+1 are related by a rotation. We fix the point o = (1, 0, . . . , 0) ∈ S n and note that the isotropy subgroup of o can be identified with SO(n) since it consists of matrices in SO(n + 1) of the form 1 0 0 R were R ∈ SO(n). Therefore, we can write S n = SO(n + 1)/SO(n). It is not possible to define homogeneous cone fields on every n-sphere. A direct way of proving this is to note that for all even-dimensional spheres S 2m , m ∈ N no global cone fields exist. This is a clear consequence of the Poincare-Brouwer theorem of algebraic topology, also known as the so-called hairy ball theorem, that the even-dimensional spheres do not admit any globally defined continuous and non-vanishing vector fields. Since a globally defined cone field on a manifold can be used to construct a continuous non-vanishing vector field by a continuous deformation of the cone field to a field of rays, the Poincare-Brouwer theorem implies the non-existence of global cone fields on S 2m . The set of all p-dimensional subspaces of R n is called the Grassmannian of dimension p in R n and is denoted by Gr(p, n) [31,1]. Grassmannians naturally arise in many applications including as parameter spaces in model estimation problems [53] and in computer vision applications including affine-invariant shape analysis, image matching, and learning theory [20,54]. The set Gr(p, n) can be endowed with a natural differentiable structure that turns it into a compact manifold of dimension p(n − p). The Grassmann manifold Gr(p, n) is a homogeneous space with a natural transitive O(n) action [12,9]: The Killing form B : g × g → R is non-degenerate for g = so(n). Thus, and m = h ⊥ with respect to the Killing form B of so(n). That is, As in the case of the n-spheres, the Poincare-Hopf theorem can be used to rule out the existence of homogeneous cone fields for most Grassmannians. Indeed, it can be shown using Schubert calculus [28] that unless p is odd and n is even, the real Grassmannian Gr(p, n) has a nonzero Euler characteristic and hence does not admit a continuous globally defined cone field as a corollary of the Poincare-Hopf theorem. In particular, Gr(p, n) does not in general admit a homogeneous cone field.

The space of positive definite matrices S + n
The space of positive definite matrices S + n of dimension n arises in many applications in information geometry and computational science. It is well known that S + n is a homogeneous space with a transitive GL(n)-action given by congruence transformations of the form The isotropy group of this action at Σ = I is precisely O(n), since τ Q : I → QIQ T = I if and only if Q ∈ O(n). Thus, we can identify any Σ ∈ S + n with an element of the quotient space GL(n)/O(n). That is The Lie algebra gl(n) of GL(n) consists of the set R n×n of all real n × n matrices equipped with the Lie bracket [ , as a sum of an antisymmetric part and a symmetric part, we have gl A cone field K on S + n is affine-invariant or homogeneous with respect to the quotient geometry S + for all Σ 1 , Σ 2 ∈ S + n . To generate such a cone field, we require a cone K(I) ⊂ m at identity that is Ad O(n) -invariant: Using such a cone, we uniquely generate a homogeneous cone field via The Ad O(n) -invariance condition (26) is satisfied if K(I) has a spectral characterization. That is, the characterization of X ∈ K(I) must only depend on the spectrum of X. For instance, tr(X) and tr(X 2 ) are spectral quantities because tr(X) is the sum of the eigenvalues of X and tr(X 2 ) is the sum of the squares of the eigenvalues of X. Therefore, these quantities are Ad O(n) -invariant. The following result gives a family of quadratic Ad O(n) -invariant cones in m, each of which generates a distinct homogeneous cone field on S + n [37].
The parameter µ controls the opening angle of the cone. If µ = 0, then (28) defines the half-space tr(X) ≥ 0. As µ increases, the opening angle of the cone becomes smaller and for µ = n (28) collapses to a ray. Now for any fixed µ ∈ (0, n), we obtain a unique well-defined affine-invariant cone field given by

Geodesics as conal curves
A cone field K M on a manifold M gives rise to a conal order ≺ on M. A continuous piecewise smooth curve γ : In the language of geometric control theory, the forward set of x is called the reachable set from x and the backward set of x is the set controllable to x. The closure ≤ K of this order is again an order and satisfies x ≤ K y if and only if y ∈ {z : x ≺ K z}. We say that M is globally orderable if ≤ K is a partial order. The conal order induced by a generic cone field on a path connected manifold is generally highly nontrivial. For instance, given a pair of points a, b ∈ M, the question of whether a and b are ordered or not is not at all straightforward to answer, since to rule out the existence of an order relation one has to demonstrate that none of an infinite collection of continuous piecewise smooth curves connecting the pair a, b is a conal curve. In this section, we discuss the significant role that geodesics play as conal curves on globally orderable Riemannian homogeneous spaces with respect to homogeneous cone fields, thereby reducing the search for an order relation between any two points a, b to a single statement on the pair of points. Thus, the use of invariant metric and conal structures on a globally orderable homogeneous space induces an order that is 'tractable' in the sense that we can check to see whether two points are ordered by checking a single condition.
First note that if G is a Lie group with a bi-invariant Riemannian metric, then the geodesics in G through the identity element e are precisely the oneparameter subgroups of G; i.e., curves γ of the form γ(t) = exp tX, where X ∈ g. That is, the Lie group exponential map coincides with the Riemannian exponential map in such cases. Geodesics through a point a ∈ G have the form γ(t) = a · exp(t a −1 · X), where X ∈ T a G. Now if G is equipped with a leftinvariant cone field K(g) = dL g | e K(e), then the geodesic γ(t) = a·exp(t a −1 ·X) through a in the direction of X is a conal curve if and only if X ∈ K(a), since Furthermore, all symmetric spaces are naturally reductive and thus g.o spaces [8].
Proposition 2 Let M = G/H be a naturally reductive homogeneous space with reductive decomposition g = h ⊕ m that is endowed with a homogeneous Riemannian metric. If K is a homogeneous cone field on M, then a geodesic γ = γ(t) through a point p ∈ M is a conal curve if and only if γ (0) ∈ K(p).
Proof By homogeneity, it is sufficient to consider geodesics through the basepoint o ∈ M, which are of the form γ(t) = exp(tX)·o, X ∈ m. Let τ g : M → M be the map τ g (x) = g · x for g ∈ G. We have as K is homogeneous. That is, γ is a conal curve if and only if its initial tangent vector lies in the cone at o. See figure 2.
We now consider the question of whether two given points a, b on a homogeneous space M = G/H equipped with a homogeneous cone field K are ordered. We first treat the example of R n endowed with the Euclidean metric and a constant cone field, i.e. invariant with respect to translations. Given a, b ∈ R n , we write a ≺ b if there exists a curve γ : [0, 1] → R n such that  γ(0) = a, γ(1) = b, and γ (t) ∈ K, for all t ∈ [0, 1]. Since K is closed and convex, we have as the integral can be thought of as the limit of a Riemann sum. But, of course, the integral is simply γ(1) − γ(0) = b − a. Thus, to check whether a and b are ordered with respect to a constant cone field, it is sufficient to check b − a ∈ K; i.e. a ≺ b if and only if the straight line from a to b is a conal curve. See figure  3. Now suppose that M = G/H is a reductive homogeneous space with a global order induced by a homogeneous cone field K. The following theorem is derived from [41].
It follows from the above result that any element of the set π −1 (↑ o) in G can be reached as the image of a vector Z in the wedge W under the exponential map. This observation can in turn be used to prove the following result. For simplicity, the theorem is formulated with respect to symmetric spaces, which include all of the examples considered in this paper.
Theorem 3 Let M = G/H be a globally orderable symmetric space with a homogeneous Riemannian metric and homogeneous cone field K. We have x 1 ≤ K x 2 if and only if the geodesic from x 1 to x 2 is a conal curve.
Proof By homogeneity, it is sufficient to consider the case where x 1 = o and x 2 = x. We define a wedge W in g by W := {Y + X : Note that W satisfies the relevant properties in Theorem 1 by construction. If o ≤ K x, it follows from Theorem 2 that there exists Z ∈ W such that By Lawson's polar decomposition theorem [30], any element g = exp Z of the semigroup S admits a unique decomposition as g = (exp X)h with X ∈ W ∩ m = K(o) and h ∈ H. Thus, we have since h · o = o for any h ∈ H. Thus, it follows that γ(t) = exp(tX) · o is a conal curve. Since all symmetric spaces are g.o. spaces, γ is precisely the Riemannian geodesic from o to x.

Corollary 1 Let M be a globally orderable symmetric space with metric and conal structures as in Theorem 3. A map F : M → M is monotone if and only if the geodesic from
is conal whenever the geodesic from x 1 to x 2 is conal.

Example: affine-invariant orders on S +
n Recall from Section 3.2 that S + n is a homogeneous space with quotient manifold structure GL(n)/O(n). The exponential map exp : m → S + n given by the usual matrix power series is a surjective map from the space of n × n symmetric matrices m onto S + n . Here m is identified with the tangent space of S + n at the identity I. Thus, the logarithm map log : S + n → m is well-defined on all of S + n . It is well-known that S + n can be equipped with a standard affineinvariant Riemannian metric ·, · Σ given by X, Y Σ = tr Σ −1 XΣ −1 Y for Σ ∈ S + n , X, Y ∈ T Σ S + n , which turns S + n = GL(n)/O(n) into a non-compact Riemannian homogeneous space with negative curvature [29]. The Riemannian distance between any two points Σ 1 , Σ 2 ∈ S + n is given by where λ i (Σ −1 1 Σ 2 ), (i = 1, . . . , n) denote the n real and positive eigenvalues of Σ −1 1 Σ 2 . It follows from (38) that d(Σ −1 1 , Σ −1 2 ) = d(Σ 1 , Σ 2 ). Thus, the inversion Σ → Σ −1 provides an involutive isometry on S + n , which shows that S + n is a Riemannian symmetric space and hence a g.o. space. Therefore, given Σ ∈ S + n , there exists a unique X = log Σ ∈ m, and the geodesic from I to Σ is given by γ(t) = exp(tX). Note that the domain of the injective curve γ can be extended to all of R.
In addition to the affine-invariant geometry of S + n , there is a natural 'flat' or translational geometry of S + n viewed as a cone embedded in the flat space of n × n symmetric matrices Sym n . The Löwner order ≥ L can be defined on Sym n by where A−B ≥ L O means that A−B is positive semidefinite. The restriction of (39) to S + n defines a partial order on S + n that coincides with the order induced on S + n ∼ = GL(n)/O(n) by the homogeneous cone field generated by the cone of positive semidefinite matrices in m ∼ = T I S + n . That is, in the special case where the cone K(I) at identity is itself the cone of positive semidefinite matrices, the translation-invariant and affine-invariant cone fields generated by K(I) agree on S + n . This is generally not the case for other choices of K(I), as shown in [37]. Now let K be an affine-invariant cone field on S + n . Such a cone field induces a global partial order on S + n [39]. Given Σ ∈ S + n , Theorem 3 implies that I ≤ K Σ if and only if the geodesic γ(t) = exp(t log Σ) from I to Σ is a conal curve; i.e., precisely if log Σ ∈ K(I). By homogeneity, it follows that Σ 1 ≤ K Σ 2 if and only if log(Σ −1/2 1 If K is the cone field corresponding to the Löwner order, then Σ 1 ≤ K Σ 2 if and only if log(Σ −1/2 1 If K is any of the quadratic affine-invariant cone fields in (29) which is equivalent to where .., n) denote the n real and positive eigenvalues of Σ 2 Σ −1 1 . Thus, the use of invariant causal structures has enabled us to answer the question of whether a pair of positive definite matrices Σ 1 and Σ 2 are ordered by checking a pair of inequalities involving the spectrum of Σ 2 Σ −1 1 .
(b) (a) Fig. 4 (a) A polyhedral cone of rank 2 in three dimensions. (b) The double cone is a quadratic cone of rank 1. The closure of its complement is a cone of rank 2. Note how a plane can be fitted within this set.

Cone fields of rank k
A closed set C in a vector space V is said to be a cone of rank k if (i) for any λ ∈ R, λ · C = C, and (ii) the maximum dimension of any subspace of V contained in C is k. Note that if K is a pointed convex cone, C = K ∪ −K is a cone of rank 1 according to the stated definition. A polyhedral cone of rank k is of the form C =K ∪ −K, whereK is given as the intersection of a collection of half-spaces. A second class of cones of rank k can be defined using quadratic forms. If P is a symmetric n × n matrix with k positive eigenvalues and n − k negative eigenvalues, then the set C = {x ∈ V : x, P x ≥ 0}, defines a cone of rank k. Note that the closure of the complement of a quadratic cone of rank k in R n is a quadratic cone of rank n − k. See figure 4 for an illustration of polyhedral and quadratic cones of rank 2 in three dimensions.
We can define homogeneous cone fields of rank k on a homogeneous space M = G/H in an analogous way to pointed and convex homogeneous cone fields. That is, a cone field C of rank k is homogeneous if it satisfies for all x, y ∈ M and g ∈ G such that y = g · x.

Monotonicity
A linear map T : V → V on a vector space V is positive with respect to a cone C of any rank k if T (C) ⊆ C. Strict positivity is characterized by T (C) ⊂ int C. A smooth map F : M → M is said to be differentially positive with respect to a cone field C = C(x) on M if dF | x C(x) ⊆ C(F (x)) for every x ∈ M [17].
Invariant differential positivity refers to differential positivity with respect to a homogeneous or invariant cone field on a homogeneous space M = G/H. Since a homogeneous cone field C satisfies C(g · x) = dτ g | x C(x), where τ g (x) = g · x is the left action of g ∈ G on x ∈ M, invariant differential positivity of F reduces to for any g 1 ∈ π −1 (x) and g 2 ∈ π −1 (F (x)), where π : G → M denotes the natural projection map. Note that (43) is a condition that is formulated in reference to a single cone C(o) ⊂ T o M. A continuous-time dynamical system with semiflow ψ = ψ(t, x) is differentially positive if the flow map ψ t : M → M, ψ t (x) = ψ(t, x) is differentially positive for any choice of t > 0. Notions of strict differential positivity and uniform strict differential positivity are defined in a natural way. In particular, uniform strict differential positivity is characterized by a cone contraction measure that is bounded below by some nonzero factor over a uniform time horizon [40]. Invariant differential positivity with respect to a pointed convex cone field can be thought of as a generalization of monotonicity. Recall that a system is monotone if it preserves a partial order. Monotonicity in a vector space with respect to a constant cone field is equivalent to differential positivity with respect to the same cone field. Similarly, if a homogeneous cone field on a homogeneous space defines a global partial order, then a map on this space is monotone if and only if it is differentially positive with respect to the same conal structure. To see this, note that a smooth map F : M → M is monotone with respect to a partial order ≥ on M if F (x 1 ) ≥ F (x 2 ) whenever x 1 ≥ x 2 . Let ≥ K denote the partial order induced by a homogeneous cone field K on M. If x 1 ≥ K x 2 , then there exists a conal curve γ : [0, 1] → M such that γ(0) = x 1 , γ(1) = x 2 and γ (t) ∈ K(γ(t)) for all t ∈ (0, 1). Now F •γ : [0, 1] → M is a curve in M with (F •γ)(0) = F (x 1 ), (F •γ)(1) = F (x 2 ), and (F • γ) (t) = dF | γ(t) γ (t). Thus, F • γ is a conal curve joining F (x 1 ) to F (x 2 ) if and only if dF | γ(t) K(γ(t)) ⊆ K(F (γ(t)); i.e., F is differentially positive with respect to K as expected. If an invariant cone field on a homogeneous space does not induce a partial order due to a failure of the antisymmetry condition arising from topological constraints, then monotonicity is not defined as it relies on the existence of a partial order. Nonetheless, invariant differential positivity provides the natural extension of the concept of monotonicity in this setting.

Example:
Order-preserving maps on S + n Invariant differential positivity can be a powerful tool for establishing monotonicity when the cone field characterization of a partial order is available. It can also be a particularly effective tool in proving monotonicity with respect to a family of partial orders corresponding to a collection of causal structures. A fundamental result in operator theory is the Löwner-Heinz theorem [10,32], which states that the map Σ → Σ r on S + n (n ≥ 2) is monotone with respect to the Löwner order if and only if r ∈ [0, 1]. The following theorem provides an extension of this important result to an infinite collection of affine-invariant causal structures on S + n , thereby highlighting the intimate connection of the Löwner-Heinz theorem to the affine-invariant geometry of S + n . The proof is based on differential positivity and can be found in [39].
Theorem 4 Let ≥ K denote the partial order induced by any quadratic affineinvariant cone field K on S + n . If Σ 1 ≥ K Σ 2 in S + n and r ∈ [0, 1], then Furthermore, if n ≥ 2 and r > 1, then the map Σ → Σ r is not monotone with respect to ≥ K .

Strict positivity
In linear positivity theory, strict positivity of a system with respect to a cone of rank k implies the existence of a dominant eigenspace of dimension k, which is an attractor for the system [19]. In the differential theory, the notion of a dominant eigenspace is replaced with that of a forward invariant distribution D of rank k corresponding to the dominant modes of the linearized system [38,18]. It is this distribution that shapes the asymptotic behavior of the dynamics. In particular, if D is involutive in the sense that for any pair of smooth vector fields X, Y defined near for all w ∈ D x , then there exists a unique integral manifold of D that is an attractor for all the trajectories from S.
Proof The strict differential positivity of Σ with respect to C on a bounded forward invariant region determines a splitting T x M = D x ⊕ D x , where D x and D x are distributions of rank k and n − k, respectively, and D x is forwardinvariant: and corresponds to the k dominant modes of the linearized system by arguments that can be found in the proofs of Theorem 1.2 of [42] and Theorem 1 of [18]. For each x ∈ M, define the map Φ x : where Π Dx , Π D x denote the linear projections onto the subspaces D x and D x . The map Φ(x) is clearly well-defined since Π Dx δx = 0 x for any δx ∈ C(x)\{0 x }. Strict differential positivity ensures that lim t→∞ Φ ψt(x) dψ t x δx = 0 for all , then for some α > 0 and w ∈ D x , we have δx + αw ∈ C(x) \ {0 x } and lim t→∞ Φ ψt(x) (δx + αw) = 0, which implies that lim t→∞ Π D ψ t (x) (dψ t | x δx) = 0 once again. Thus, in the limit of t → ∞, δx(t) = dψ t | x δx becomes parallel to D ψt(x) . It follows that any curve Γ = Γ (s) in S evolves so that ψ t (Γ (s)) asymptotically lies on an integral manifold of D. To prove uniqueness of the attractor, assume for contradiction that N 1 and N 2 are two distinct attractive integral manifolds of D and let x 1 ∈ N 1 , x 2 ∈ N 2 . By connectedness of S, there exists a smooth curve Γ in S connecting x 1 and x 2 . Since the curve ψ t (Γ (s)) converges to an integral manifold of D, N 1 and N 2 must be subsets of the same integral manifold of D, which provides the contradiction that completes the proof.
Note that condition (45) is necessary to ensure that vectors along the distribution D do not grow unbounded as they evolve by the variational flow, thereby ensuring that strict differential positivity results in contraction toward an integral manifold of D.
As noted earlier, invariant differential positivity is a generalization of monotonicity to homogeneous spaces. This generalization is made possible by the nonlinearity of the homogeneous space and allows for more complex asymptotic behavior to arise. For instance, as shown in the work of M. Hirsch [26,50], almost all bounded trajectories of a strongly monotone (strictly differentially positive) system converge to the set of equilibria. Moreover, under mild smoothness and boundedness assumptions, almost every trajectory converges to one equilibrium. On the other hand, for systems that are strictly differentially positive with respect to pointed convex invariant cone fields on homogeneous spaces, almost all bounded trajectories may converge to a limit cycle under similarly mild technical assumptions [17]. For example, invariant differential positivity on the cylinder has been used to establish convergence to a unique limit cycle in a nonlinear pendulum model [16]. Such a one-dimensional asymptotic behavior is possible for an invariantly differentially positive system due to the topology of the cylinder, which allows for the existence of closed conal curves.

Coset stabilization on Lie groups
In many problems in dynamical systems and control theory, we are interested in asymptotic convergence of trajectories to a submanifold of the state space. Such problems arise in numerous applications including consensus, synchronization, pattern generation, and path following. Here we consider a special class of such systems which are defined on a Lie group G and converge to submanifolds which arise as integrals of left-invariant distributions on G. A leftinvariant distribution D on G is a distribution that satisfies D g1g2 = dL g1 | g2 D g2 for all g 1 , g 2 ∈ G. Any such distribution uniquely determines a subspace D e of T e G ∼ = g and conversely every subspace of T e G defines a unique left-invariant distribution. Furthermore, the left-invariant distribution defined by a subspace U of the Lie algebra g is integrable if and only if U is a subalgebra h of g. Given such a left-invariant distribution D, its integral through the identity element e ∈ G is a subgroup H of G with Lie algebra h. The integral of D through any other point g ∈ G corresponds to a translation of the subgroup H on G and can be identified with the left coset gH.
In many applications, we are interested in stabilizing a submanifold corresponding to a coset gH in G. For instance, in satellite surveillance the attitude of the satellite is an element of the special orthogonal group SO(3) and we often seek to control the orientation of the satellite by ensuring that the telescope axis points to a fixed point on the Earth's surface. Since the set of attitudes which solve this problem can be identified with an SO(2) subgroup corresponding to rotations about the telescope axis, the problem is essentially a simple coset stabilization problem [34]. Another large class of coset stabilization problems involves consensus models involving N agents g k whose states evolve on a Lie group G and exchange information about their relative positions via a communication graph. The consensus manifolds in such problems take the form of a single copy of G repeated N times and diagonally embedded in the Cartesian product of G N , which is the (N dim G)-dimensional state space of the system. Such consensus manifolds correspond to fixed formations of the N agents evolving uniformly on G and correspond to (dim G)-dimensional cosets in G N . The synchronization manifold on which all N agents have the same state is a special case of such a coset in G N . A simple example of such a model is a network of N oscillators evolving on the N -torus, where the cosets are onedimensional and correspond to frequency synchronization and phase-locking behaviors. Now consider a homogeneous space M = G/H, where G is a semisimple and compact Lie group so that the negative of the Killing form −B of g is positive definite. Then M is reductive with reductive decomposition g = h⊕m, where m = h ⊥ . The negative of the Killing form −B induces a bi-invariant metric on G and a homogeneous metric on M. The tangent space at each point g ∈ G admits the decomposition T g G = H g ⊕H ⊥ g , where H = dL g | e h and H ⊥ = dL g | e m. We choose a basis for g of the form {e 1 , ..., e dim h , e dim h+1 , ..., e dim g }, where {e i : i = 1, ..., dim h} is a basis of h. Define a cone C(e) of rank k = dim h in g by where µ > 0 is a sufficiently small parameter to ensure that (48) is non-empty. C(e) can be uniquely extended to a left-invariant cone field C(g) = dL g | e C(e) of rank dim h on G.
Note that a dynamical systemġ =f (g) on G induces a well-defined projected flow on M generated byẋ = f (x) if for all g 1 ∼ g 2 , where ∼ denotes the equivalence relation induced by the coset manifold structure G/H; i.e., g 1 ∼ g 2 if and only if there exists h ∈ H such that g 1 = g 2 h.f H ⊥ (g) in (49) refers to the component off in H ⊥ g ⊂ T g G. The following theorem shows how strict differential positivity with respect to the left-invariant cone field C in G can induce a contractive dynamical system on M under suitable technical conditions. See figures 5 and 6 for illustrations of the relevant concepts.
Theorem 6 Let Σ :ġ =f (g) be a system on G that induces a well-defined projective dynamicsẋ = f (x) on M = G/H. Suppose that Σ is uniformly strictly differentially positive with respect to the left-invariant cone field C determined by (48) on a bounded, connected and forward invariant region S ⊂ G, with a distribution of dominant eigenspaces of the form D g = dL g | e h. If lim sup t→∞ dψ t | g w ψt(g) < ∞, for all w ∈ D g , then there exists a unique coset gH in S that is an attractor for all trajectories from S. Furthermore, the induced system on M is contractive with respect to any homogeneous metric and all trajectories from π(S) ⊂ M converge to the fixed point x = π(gH) ∈ π(S). set of vertices V and edges E evolving on G, we can define a class of consensus protocols on G as follows. For each g k ∈ G denote the Riemannian exponential and logarithm maps by exp g k : T g k G → G and log g k : U g k → T g k G, respectively, where U g k ⊂ G is the maximal set containing g k for which exp g k is a diffeomorphism. The system given bẏ defines a consensus protocol on G for constant vectors Ω k ∈ g and any collection of real-valued reshaping functions µ ki of the distance d that satisfy µ ki (0) = 0. Equation (50) defines a dynamical system on G × · · · × G = G N that yields a well-defined projected system on G N /G in the sense of Eq. (49) and Theorem 6, since log g·g k (g · g i ) = g · (log g k g i ) and d(g · g k , g · g i ) = d(g k , g i ), for all g, g k , g i ∈ G leave Eq. (50) invariant. For the sake of simplicity, we will consider the example of a network of N agents evolving on the circle S 1 . Equation (50) reduces to a system of the forṁ where θ k ∈ S 1 represents the phase of agent k, ω k ∈ R are prescribed 'intrinsic' frequencies, and µ ki now denotes an odd coupling function on the domain (−π, π) extended to R in such a way so as to make it 2π-periodic. Note that µ ki and µ ik need not be the same function. Let θ = (θ 1 , . . . , θ N ) denote an element of the N -torus T N and consider the N -tuple of vector fields ∂ ∂θ 1 , . . . , ∂ ∂θ N , which defines a basis of left-invariant vector fields on T N . Assuming that the coupling functions µ ki are differentiable and strictly monotonically increasing on (−π, π), then it can be shown that the linearizationδθ = A(θ)δθ of the system given by (52) is uniformly strictly differentially positive on the set T N π = {ϑ ∈ T N : |ϑ k − ϑ i | < π, (i, k) ∈ E} with respect to the invariant cone field for any strongly connected communication graph. Furthermore, the Perron-Frobenius vector field of the system on T N π is the left-invariant vector field 1(θ) = (1, . . . , 1) ∈ T θ T N , where the vector representation is given with respect to the invariant basis defined by ∂ ∂θ 1 , . . . , ∂ ∂θ N . Moreover, if we denote the flow of (52) by ψ t , then the condition A(θ)1(θ) = 0 implies that dψ t | θ 1 θ = 1 ψt(θ) , which ensures that lim sup t→∞ dψ t | θ 1(θ) ψt(θ) < ∞ for any flow confined to T N π , where · θ denotes the norm corresponding to the standard Riemannian metric on T N . If we add the requirement that the coupling functions µ ki be barrier functions on (−π, π) so that µ ki (α) → ∞ as α → π, then the flow ψ t will be forward-invariant on T N π , resulting in the following theorem.
Note that convergence to an integral curve of 1 on T N corresponds to a phase-locking behavior, whereby the collective motion asymptotically converges to movement in a fixed formation with frequency synchronization among the agents. Further details may be found in [40].

Conclusion
We have reviewed the notion of invariant cone fields on homogeneous spaces and presented examples of cone fields on a list of homogeneous spaces that are of special interest in applications in information science. Invariant differential positivity naturally arises as a generalization of monotonicity on homogeneous spaces within this context. Finally, we have illustrated the potential applications of monotone flows on homogeneous spaces in systems and control theory by reviewing the consensus problem on S 1 and discussing possible extensions of the approach to higher-dimensional spaces.