Constructing Subspace Packings from Other Packings
Emily J. King

TL;DR
This paper explores methods to construct subspace packings from existing configurations, analyzing their algebraic and geometric properties, and extends the theory to include various optimal packings like equiisoclinic and orthoplex-bound arrangements.
Contribution
It generalizes known construction techniques for subspace packings, characterizes preserved properties, and analyzes new configurations including those outside traditional parameter regimes.
Findings
Provides a unified framework for constructing subspace packings.
Characterizes algebraic and geometric properties preserved in constructions.
Includes analysis of optimal packings like orthoplex-bound arrangements.
Abstract
The desirable properties when constructing collections of subspaces often include the algebraic constraint that the projections onto the subspaces yield a resolution of the identity like the projections onto lines spanned by vectors of an orthonormal basis (the so-called tightness condition) and the geometric constraint that the subspaces form an optimal packing of the Grassmannian, again like the one-dimensional subspaces spanned by vectors in an orthonormal basis. In this article a generalization of related constructions which use known packings to build new configurations and which appear in numerous forms in the literature is given, as well as the characterization of a long list of desirable algebraic and geometric properties which the construction preserves. Another construction based on subspace complementation is similarly analyzed. While many papers on subspace packings focus…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Constructing Subspace Packings from Other Packings
Emily J. King
Colorado State University
Abstract
The desirable properties when constructing collections of subspaces often include the algebraic constraint that the projections onto the subspaces yield a resolution of the identity like the projections onto lines spanned by vectors of an orthonormal basis (the so-called tightness condition) and the geometric constraint that the subspaces form an optimal packing of the Grassmannian, again like the one-dimensional subspaces spanned by vectors in an orthonormal basis. In this article a generalization of related constructions which use known packings to build new configurations and which appear in numerous forms in the literature is given, as well as the characterization of a long list of desirable algebraic and geometric properties which the construction preserves. Another construction based on subspace complementation is similarly analyzed. While many papers on subspace packings focus only on so-called equiisoclinic or equichordal arrangements, attention is also given to other configurations like those which saturate the orthoplex bound and thus are optimal but lie outside of the parameter regime where equiisoclinic and equichordal packings can occur. Keywords: fusion frame, Grassmannian packing, simplex bound, orthoplex bound, equichordal, strongly simplicial, equiisoclinic MSC 2010: 42C15, 14M15
1 Introduction
1.1 Motivation
The goal is to find optimal configurations of subspaces which are of interest mathematically but also in applications such as coding theory (see, e.g., [Cre08, PWTH18, XZG05, KP03, KPCL09]), quantum information theory (see, e.g., [FHS17, AFZ15, SS98, GR09]), and more. The usefulness of such configurations often comes from whether the projections onto the subspaces (approximately) yield a resolution of the identity and whether the angles between the subspaces are as large as possible. Configurations which satisfy the former condition are called (tight) fusion frames and the latter condition are called Grassmannian packings.
In general, Grassmannian fusion frames, which are fusion frames which correspond to optimal packings of Grassmannian spaces with respect to the chordal distance, are optimally robust against noise and erasures. Under certain models of noise and erasures, a type of Grassmannian fusion frame called equichordal is shown to be optimal [KPCL09]. Under other models, the subclass of Grassmannian fusion frames which are called equiisoclinic have been proven to be the best [Bod07], while for certain coding theory regimes such packings are not optimal [PWTH18]. Equiisoclinic tight fusion frames are also proven to be deterministically optimal in block sparse recovery [SAH14, EKB10].
We complete this section by introducing Grassmannian fusion frames and basic notation (Section 1.2). In Section 2, we give generalizations, classifications, and examples of constructions from the literature of various Grassmannian packings that can be constructed from other such packings. In particular, a construction of Grassmannian packings and fusion frames from ones known to exist which generalizes different constructions based on Kronecker products in [LS73b, BCP*+*13, Cre08, SAH14, CFM*+*11, CTX15, EKB10, MKGH20] is presented in Theorem 2.2 and Corollary 2.3, coupled with proofs of which desirable algebraic and geometric properties are inherited from the original packings. Then subspace complementation, a construction method found in, e.g. [BCP*+*13], is characterized including a focus on optimal Grassmannian packings which are not equichordal (Proposition 2.5).
Throughout the paper, will always either denote or . Further for , we define , to be the set of matrices with entries in , and to be the identity matrix. Finally, for , we write for the column span of .
1.2 Fusion Frames and Grassmannian Packings
Our objects of interest are collections of subspaces, which may be viewed as points in a Grassmannian.
Definition 1.1**.**
For , set to be the collection of dimensional subspaces of . is called a Grassmannian. is endowed with a metric space structure induced by the chordal distance (see, e.g., [CHS96])
[TABLE]
for , where is the orthogonal projection onto .
We do not consider the trivial cases and as they are not interesting. We also note that there are many other metrics that one can endow the Grassmannian with, like the spectral distance, the geodesic distance, and the Fubini-Study distance (see, e.g., [CHS96, DHST08]). Let us consider a set of vectors for and define for each . Then is an orthonormal basis precisely when is the orthogonal projection onto for each , , and
[TABLE]
We now generalize these traits to systems which may be overcomplete and consist of subspaces of dimension greater than . Collections of possibly redundant systems of subspaces which yield algebraic reconstruction properties akin to orthonormal bases have appeared in the literature for many years with many names, including stable space splittings of Hilbert spaces [Osw94, Osw97], systems of bounded quasi-projectors [For03, For04], (weighted projective) resolutions of the identity [For03, Bod07], -frames [Sun06], and frames of subspaces [CK04]. We will use the name fusion frames and the related terminology, which were introduced in [CKL08]. See [CK13, Chapter 13] for a general overview of fusion frames.
Definition 1.2**.**
A finite collection of subspaces is a tight fusion frame of -dimensional subspaces with unit weights for if there exists an (called the fusion frame bound) satisfying
[TABLE]
where is the orthogonal projection onto . The map is called the fusion frame operator.
One may loosen this definition by allowing non-equidimensional subspaces, non-equal weights, and the fusion frame operator to only be an approximation of the identity; however, we will not be concerned with these cases in this paper. To avoid being verbose, we shall refer to tight fusion frames of -dimensional subspaces with unit weights as tight fusion frames. Given , we fix for each an orthonormal basis for the subspace and denote by the matrix . We further define
[TABLE]
Then it is clear that (3) holds, i.e., that is a tight fusion frame of -dimensional subspaces with unit weights for with fusion frame bound precisely with the rows of are orthogonal with norm .
We will consider three definitions of “equal” geometric spread between subspaces. See [LS73b, ET07, BCP*+*13, BG73, Cre08], in particular [Theorem 2.3] of [LS73b], [Theorem 1] of [BG73], and [p. 4] of [Cre08].
Definition 1.3**.**
Let (not necessarily a fusion frame) with corresponding orthonormal bases as the columns of . Then we say
- •
is equichordal when for all with , is constant;
- •
is strongly simplicial when for all with , has the same set of eigenvalues; and
- •
is equiisoclinic when there exists an such that for all with , .
If are the eigenvalues of , then for all the which satisfy are called the principal angles between and .
We note that given a set of parameters, the only variable in the definition of the chordal distance (1) is ; thus, being equichordal means that is constant for , as one would hope given the name.
If we fix , , and , then the Grassmannian packing problem concerns finding elements in so that the minimal distance between any two subspaces is as large as possible, just like in (2). An algorithm to approximate solutions is in [DHST08], while [Slo] has a (somewhat dated) list of best known packings when . For fixed parameters , , and , the maximizers of are called Grassmannian fusion frames (which may or may not be tight fusion frames).
Definition 1.4 and Theorem 1.5 summarize results concerning the Grassmannian packing problem found in [Ran55, CHS96, KPCL09, LS73a, BH16, Hen05, FJMW17] and many more sources.
Definition 1.4**.**
Define
[TABLE]
This is known as Gerzon’s bound and comes from the dimension of the smallest vector subspace of which contains the symmetric / self-adjoint matrices.
Theorem 1.5**.**
Let , then
[TABLE]
The bound in (4) is saturated if and only if is an equichordal tight fusion frame. Further, the bound (4) can only be saturated if . When ,
[TABLE]
Thus if a tight fusion frame is equichordal, strongly simplicial, or equiisoclinic it is a Grassmannian fusion frame since all of those configurations are equichordal. The bound in (4) is known as the simplex bound, and the bound in (5) is the orthoplex bound. If the orthoplex bound is saturated and , then we call an orthoplectic Grassmannian packing. Such a packing need not be a tight fusion frame. We further note that when , it is possible for to saturate the orthoplex bound but not be a Grassmannian fusion frame.
There are some special terms for the above-defined concepts when (see, e.g., [JKM19, CK13, Wal18]). Namely, tight fusion frames of -dimensional subspaces with unit weights are called finite unit norm tight frames. Also, the definitions of equichordal, strongly simplicial, and equiisoclinic coincide and are jointly known as equiangular.
Some constructions and characterizations of equiisoclinic packings are in [ET06, ET07, ET18, Hog77, LS73b, CHR*+*99, SS98], of strongly simplicial packings are in [Cre08], of equichordal packings are in [Cre08, CHR*+*99, KPCL09, BP15, Kin13, Kin19], and of orthoplectic Grassmannian packings are in [SS98, BH16, GR09].
2 Constructing New Grassmannian Packings from Old
A common method for constructing new optimal Grassmannian packings out of already known ones is to use a Kronecker product or Kronecker-like product. We make note of the following standard definition and properties.
Definition 2.1**.**
For and , with the th row and th column of denoted by , we define the Kronecker product as
[TABLE]
A certainly incomplete summary of such constructions making use of the Kronecker product in the literature follows. We will often use the more general term tensor instead of Kronecker product to circumvent awkward phrasing. In [LS73b], the authors tensor real equiangular lines (not necessarily a frame) with orthogonal matrices to obtain equiisoclinic packings (not necessarily a fusion frame). They achieve this indirectly by tensoring the Gram matrix of the equiangular lines with an identity matrix ([proof of Theorem 3.7], [Remark 4.2]). The same construction was repeated in [SAH14, Theorem 4] but specifically for the case of equiangular tight frames (although the authors falsely assume that all Grassmannian frames are equiangular tight frames). Similarly, in [MKGH20], orthonormal basis vectors for equichordal tight fusion frames are tensored with an identity to generate equichordal tight fusion frames. The constructions in [CTX15] involve tensoring vectors in an equiangular tight frame [Theorem 3] or mutually unbiased bases [Theorem 6] (see also [EKB10]) with a fixed unitary and then considering the spectral distance – which encourages equiisoclinic packings – between the subspaces. In [Hog77], the concept of isoclinic covariant functors as a method of creating new sets of equiisoclinic subspaces from pre-existing ones is introduced. Tensoring vectors with an identity matrix would be one example of such a functor. Instead of starting with an equiangular tight frame, the author of [Cre08] tensors the projections of an equichordal Grassmannian fusion frame with an identity matrix to obtain another equichordal Grassmannian fusion frame in [Proposition 12]. Finally, in [BCP*+*13, Corollary 5], each column vector in a truncated orthonormal basis is tensored with a unitary matrix representing the different subspaces of a particular type of tight fusion frame to obtain a new tight fusion frame, called the Naimark complement, and in [Theorem 7], a specific application of this starting with a collection of orthonormal bases for the entire space is used to obtain a equiisoclinic Grassmannian fusion frame. Furthermore, [CFM*+*11] essentially has the construction from Corollary 2.3 but in the specific case of building new tight fusion frames from old ones, without a focus on the geometric properties. We now present a construction which subsumes the ones listed above. The theorem also shows that the construction preserves various desired properties separately.
Theorem 2.2**.**
Let be a collection of -dimensional subspaces in . For each , fix an orthonormal basis of . Further let be a collection of unitaries in . For each and define
[TABLE]
Then the following statements hold.
For each , the columns of are a set of orthonormal vectors in ; 2. 2.
* is equichordal (not necessarily a fusion frame) if and only if is;* 3. 3.
* is strongly simplicial (not necessarily a fusion frame) if and only if is;* 4. 4.
* is equiisoclinic (not necessarily a fusion frame) if and only if is;* 5. 5.
* saturates the orthoplex bound if and only if does, but they cannot both be orthoplectic Grassmannian packings if ; and* 6. 6.
* is a tight fusion frame if and only if is. In this case, they have the same fusion frame bound.*
Proof.
We begin by computing blocks of for , making use of basic properties of the Kronecker product. We note that for ,
[TABLE]
When , Equation 6 simplifies to
[TABLE]
where is the zero matrix. In this case, the off-diagonal blocks of the Gram matrix are the zero matrix and the diagonal blocks are the identity. Thus the columns of are orthonormal, and Statement 1 is proven.
To prove Statements 2 – 5, we are concerned with the properties of , or more specifically , when . As usual, we define to be the matrix with columns (from the original packing). Then the entry of is precisely . Thus, it follows from Equation 6 that
[TABLE]
It follows from Definition 1.3, that is equichordal if and only if is constant for all , which happens if and only if is constant for all . Thus, Statement 2 is proven. Similarly, has spectrum independent of (resp. is for all equal to a constant multiple times the identity) if and only if has spectrum independent of (resp. is for all equal to a constant multiple times the identity). Hence, Statements 3 and 4 are proven.
If (resp., ) saturates the orthoplex bound, then the maximum value of the chordal distance satisfies for some
[TABLE]
We can see from Equation 8 that . Thus has at least one pair of subspaces at the orthoplex bound if and only if does as well. However, such a configuration is only optimal when (resp., ). Unless , cannot fall in both ranges.
To prove Statement 6, we define
[TABLE]
We will make use of the fact that (resp., ) is a tight fusion frame if and only if (resp., ) is a constant multiple of the identity. We calculate the respective matrix products as the sum of the rank one tensors formed from the columns.
[TABLE]
It follows that if and only if . ∎
We may generalize Theorem 2.2.
Corollary 2.3**.**
Let be a collection of -dimensional subspaces in and be a collection of -dimensional subspaces in . For each , fix an orthonormal basis of and of . Further define for each . Finally for each and define
[TABLE]
Then the following statements hold.
For each , the columns of are a set of orthonormal vectors in ; 2. 2.
* is equichordal (not necessarily a fusion frame) if and only if and are;* 3. 3.
* is strongly simplicial (not necessarily a fusion frame) if and only if and are;* 4. 4.
* is equiisoclinic (not necessarily a fusion frame) if and only if and are; and* 5. 6.
* is a tight fusion frame if and only if and are. In this case, the fusion frame bound of is the product of the fusion frame bounds of and .*
Proof.
The proof of the corollary follows the proof of Theorem 2.2 quite closely. We note that the s are partial isometries, so when , Equation 6 still simplifies, and the off-diagonal blocks of the Gram matrix are the zero matrix and the diagonal blocks are the identity. Thus the columns of are orthonormal, and Statement 1 is proven.
For the proof of Statements 2 – 4, we can only simplify down to Equation 7. That is,
[TABLE]
We now make use of the fact that for arbitrary square matrices and , and for arbitrary matrices the singular values of are simply the products of all of the singular values of with each of the singular values of .
Statement 5 from Theorem 2.2 does not carry over well.
Finally, we note that Statement 6 is simply [Theorem 4] from [CFM*+*11]. ∎
Another way to generate a fusion frame from another fusion frame is via subspace complementation [BCP*+*13].
Lemma 2.4** ([MBI92, QZL05]).**
Let and be subspaces of . The nonzero principal angles between and are equal to the nonzero principal angles between their orthogonal complements and .
Proposition 2.5**.**
Let . Then the following statements about hold.
* is a tight fusion frame with bound if and only if is a tight fusion frame with bound ;* 2. 2.
* is equichordal (not necessarily a fusion frame) if and only if is;* 3. 3.
* is strongly simplicial (not necessarily a fusion frame) if and only if is;* 4. 4.
If is equiisoclinic (not necessarily a fusion frame), then is equiisoclinic (not necessarily a fusion frame) if and only if or the subspaces are trivially all the same; and 5. 5.
* is an orthoplectic Grassmannian packing if and only if is.*
Proof.
Statement 1 is [BCP*+*13, Theorem 5]. We note that for tight fusion frames, the fusion frame bound is ; hence, . Statements 2 and 3 follow immediately from Lemma 2.4. If is equiisoclinic, then each pair of subspaces has equal principal angles. If the principal angles are all [math], then the are all equal and hence so are the . Otherwise the principal angles are in . By Lemma 2.4, the non-zero principal angles of the are the same, meaning that one must have for equiisoclinicity to be preserved, showing Statement 4. Finally, for Statement 5, is an orthoplectic Grassmannian packing by Theorem 1.5 if and there exist with such that
[TABLE]
Since the number of subspaces and the dimension of the base space does not change when taking orthogonal complements, being an orthoplectic Grassmannian packing implies that is as well and vice versa. ∎
There are two huge differences when comparing the results of subspace complementation with the tensor construction (Theorem 2.2). Initially, in contrast to the tensor construction which always destroys the optimality of an orthoplectic Grassmannian packing, subspace complementation preserves it. We also note that subspace complementation does not in general preserve equiisoclinicity. We can see that in the following example.
Example 2.6**.**
Let be the subspaces spanned by the columns of the Sylvester-Hadamard matrix with the first row removed:
[TABLE]
This fusion frame of -dimensional subspaces is actually an equiangular tight frame and thus trivially equiisoclinic. Then the columns of the following matrices serve as orthonormal bases for with and
[TABLE]
Although is a strongly simplicial tight fusion frame, it is not equiisoclinic since, for example
[TABLE]
On the other hand, Theorem 2.5.4 is not vacuously true, as seen in the following example.
Example 2.7**.**
We may use Theorem 2.2 to generate equiisoclinic systems in any even dimensional space consisting of up to subspaces of dimension . To do this, we start with any equiangular system of vectors in and tensor them with unitaries for any . It follows from Theorem 2.2.4 that the result is an equiisoclinic collection of -dimensional subspaces of . Gerzon’s bound (Theorem 1.5) implies that if and if ; it ends up that equiangular lines exist for each possible . Any pair of non-orthogonal unit vectors will generate an equiangular system in that is not a tight frame, while a pair of orthogonal unit vectors yields an equiangular system that is a tight frame. The Mercedes-Benz frame is an example of an equiangular tight frame of vectors in :
[TABLE]
and the symmetric informationally complete operator-valued measure in is an example of an equiangular tight frame of vectors in :
[TABLE]
A particularly nice presentation of an equiisoclinic tight fusion frame of -dimensional subspaces of resulting from (the complement of) a Singer difference set may be found in [Example 3.7] of [FS20].
Acknowledgements
The author was supported in part by the Explorationsprojekt “Hilbert Space Frames and Algebraic Geometry” funded by the Zentrum für Forschungsförderung der Univeristät Bremen. The author would like to thank the anonymous referee for their helpful suggestions.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[AFZ 15] D. M. Appleby, Christopher A. Fuchs, and Huangjun Zhu, Group theoretic, Lie algebraic and Jordan algebraic formulations of the SIC existence problem , Quantum Inf. Comput. 15 (2015), no. 1-2, 61–94. MR 3362195
- 2[BCP + 13] Bernhard G. Bodmann, Peter G. Casazza, Jesse D. Peterson, Ihar Smalyanau, and Janet C. Tremain, Excursions in harmonic analysis , vol. 1, ch. Fusion Frames and Unbiased Basic Sequences, pp. 19–34, Birkhäuser Boston, 2013.
- 3[BG 73] Ȧke Björck and Gene H. Golub, Numerical methods for computing angles between linear subspaces , Math. Comp. 27 (1973), 579–594. MR 0348991
- 4[BH 16] Bernhard G. Bodmann and John Haas, Achieving the orthoplex bound and constructing weighted complex projective 2-designs with Singer sets , Linear Algebra Appl. 511 (2016), 54–71.
- 5[Bod 07] Bernhard G. Bodmann, Optimal linear transmission by loss-insensitive packet encoding , Appl. Comput. Harmon. Anal. 22 (2007), no. 3, 274–285. MR MR 2311854 (2008 c:42031)
- 6[BP 15] Irena Bojarovska and Victoria Paternostro, Gabor fusion frames generated by difference sets , Wavelets and Sparsity XVI (Manos Papadakis and Vivek K. Goyal, eds.), Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 9597, 2015.
- 7[CFM + 11] Peter G. Casazza, Matthew Fickus, Dustin G. Mixon, Yang Wang, and Zhengfang Zhou, Constructing tight fusion frames , Appl. Comput. Harmon. Anal. 30 (2011), 175–187.
- 8[CHR + 99] A. R. Calderbank, R. H. Hardin, E. M. Rains, P. W. Shor, and N. J. A. Sloane, A group-theoretic framework for the construction of packings in Grassmannian spaces , J. Algebraic Combin. 9 (1999), no. 2, 129–140.
