Stiefel manifold

From Wikipedia the free encyclopedia

In mathematics, the Stiefel manifold is the set of all orthonormal k-frames in That is, it is the set of ordered orthonormal k-tuples of vectors in It is named after Swiss mathematician Eduard Stiefel. Likewise one can define the complex Stiefel manifold of orthonormal k-frames in and the quaternionic Stiefel manifold of orthonormal k-frames in . More generally, the construction applies to any real, complex, or quaternionic inner product space.

In some contexts, a non-compact Stiefel manifold is defined as the set of all linearly independent k-frames in or this is homotopy equivalent to the more restrictive definition, as the compact Stiefel manifold is a deformation retract of the non-compact one, by employing the Gram–Schmidt process. Statements about the non-compact form correspond to those for the compact form, replacing the orthogonal group (or unitary or symplectic group) with the general linear group.

Topology

[edit]

Let stand for or The Stiefel manifold can be thought of as a set of n × k matrices by writing a k-frame as a matrix of k column vectors in The orthonormality condition is expressed by A*A = where A* denotes the conjugate transpose of A and denotes the k × k identity matrix. We then have

The topology on is the subspace topology inherited from With this topology is a compact manifold whose dimension is given by

As a homogeneous space

[edit]

Each of the Stiefel manifolds can be viewed as a homogeneous space for the action of a classical group in a natural manner.

Every orthogonal transformation of a k-frame in results in another k-frame, and any two k-frames are related by some orthogonal transformation. In other words, the orthogonal group O(n) acts transitively on The stabilizer subgroup of a given frame is the subgroup isomorphic to O(nk) which acts nontrivially on the orthogonal complement of the space spanned by that frame.

Likewise the unitary group U(n) acts transitively on with stabilizer subgroup U(nk) and the symplectic group Sp(n) acts transitively on with stabilizer subgroup Sp(nk).

In each case can be viewed as a homogeneous space:

When k = n, the corresponding action is free so that the Stiefel manifold is a principal homogeneous space for the corresponding classical group.

When k is strictly less than n then the special orthogonal group SO(n) also acts transitively on with stabilizer subgroup isomorphic to SO(nk) so that

The same holds for the action of the special unitary group on

Thus for k = n − 1, the Stiefel manifold is a principal homogeneous space for the corresponding special classical group.

Uniform measure

[edit]

The Stiefel manifold can be equipped with a uniform measure, i.e. a Borel measure that is invariant under the action of the groups noted above. For example, which is isomorphic to the unit circle in the Euclidean plane, has as its uniform measure the natural uniform measure (arc length) on the circle. It is straightforward to sample this measure on using Gaussian random matrices: if is a random matrix with independent entries identically distributed according to the standard normal distribution on and A = QR is the QR factorization of A, then the matrices, are independent random variables and Q is distributed according to the uniform measure on This result is a consequence of the Bartlett decomposition theorem.[1]

Special cases

[edit]

A 1-frame in is nothing but a unit vector, so the Stiefel manifold is just the unit sphere in Therefore:

Given a 2-frame in let the first vector define a point in Sn−1 and the second a unit tangent vector to the sphere at that point. In this way, the Stiefel manifold may be identified with the unit tangent bundle to Sn−1.

When k = n or n−1 we saw in the previous section that is a principal homogeneous space, and therefore diffeomorphic to the corresponding classical group:

Functoriality

[edit]

Given an orthogonal inclusion between vector spaces the image of a set of k orthonormal vectors is orthonormal, so there is an induced closed inclusion of Stiefel manifolds, and this is functorial. More subtly, given an n-dimensional vector space X, the dual basis construction gives a bijection between bases for X and bases for the dual space which is continuous, and thus yields a homeomorphism of top Stiefel manifolds This is also functorial for isomorphisms of vector spaces.

As a principal bundle

[edit]

There is a natural projection

from the Stiefel manifold to the Grassmannian of k-planes in which sends a k-frame to the subspace spanned by that frame. The fiber over a given point P in is the set of all orthonormal k-frames contained in the space P.

This projection has the structure of a principal G-bundle where G is the associated classical group of degree k. Take the real case for concreteness. There is a natural right action of O(k) on which rotates a k-frame in the space it spans. This action is free but not transitive. The orbits of this action are precisely the orthonormal k-frames spanning a given k-dimensional subspace; that is, they are the fibers of the map p. Similar arguments hold in the complex and quaternionic cases.

We then have a sequence of principal bundles:

The vector bundles associated to these principal bundles via the natural action of G on are just the tautological bundles over the Grassmannians. In other words, the Stiefel manifold is the orthogonal, unitary, or symplectic frame bundle associated to the tautological bundle on a Grassmannian.

When one passes to the limit, these bundles become the universal bundles for the classical groups.

Homotopy

[edit]

The Stiefel manifolds fit into a family of fibrations:

thus the first non-trivial homotopy group of the space is in dimension n − k. Moreover,

This result is used in the obstruction-theoretic definition of Stiefel–Whitney classes.

See also

[edit]

References

[edit]
  1. ^ Muirhead, Robb J. (1982). Aspects of Multivariate Statistical Theory. John Wiley & Sons, Inc., New York. pp. xix+673. ISBN 0-471-09442-0.
  2. ^ Chikuse, Yasuko (1 May 2003). "Concentrated matrix Langevin distributions". Journal of Multivariate Analysis. 85 (2): 375–394. doi:10.1016/S0047-259X(02)00065-9. ISSN 0047-259X.
  3. ^ Pal, Subhadip; Sengupta, Subhajit; Mitra, Riten; Banerjee, Arunava (September 2020). "Conjugate Priors and Posterior Inference for the Matrix Langevin Distribution on the Stiefel Manifold". Bayesian Analysis. 15 (3): 871–908. doi:10.1214/19-BA1176. ISSN 1936-0975.