Admissibility and Frame Homotopy for Quaternionic Frames

Abstract: We consider the following questions: when do there exist quaternionic frames with given frame spectrum and given frame vector norms? When such frames exist, is it always possible to interpolate between any two while fixing their spectra and norms? In other words, the first question is the admissibility question for quaternionic frames and the second is a generalization of the frame homotopy conjecture. We give complete answers to both questions. For the first question, the existence criterion is exactly the same as in the real and complex cases. For the second, the non-empty spaces of quaternionic frames with specified frame spectrum and frame vector norms are always path-connected, just as in the complex case. Our strategy for proving these results is based on interpreting equivalence classes of frames with given frame spectrum as adjoint orbits, which is an approach that is also well-suited to the study of real and complex frames.

\section{Introduction}\label{sec:intro}

In a finite-dimensional real or complex Hilbert space $\mathfrak{H}$, an (ordered) \emph{frame} is simply a collection $f_1, \dots , f_N \in \mathfrak{H}$ that spans $\mathfrak{H}$. When $N=d:=\dim \mathfrak{H}$, this is just a basis for $\mathfrak{H}$, but when $N > d$ a frame gives a redundant representation of a signal $v \in \mathfrak{H}$ by
\begin{equation}\label{eq:analysis}
	(\langle v, f_1 \rangle, \dots , \langle v, f_N \rangle)
\end{equation}
which can be more robust to erasures and other corruption of the data than a basis representation~\cite{Casazza:2003vp,Goyal:2001cd,Holmes:2004iv}.

Identifying the frame with the matrix $F = [f_1 | \dots | f_N]$ whose columns are the frame vectors, the redundant representation~\eqref{eq:analysis} corresponds to evaluation of the \emph{analysis operator} $v \mapsto F^\ast v$. Composition of the analysis operator with its adjoint \emph{synthesis operator} $w \mapsto F w$ gives the \emph{frame operator}
\[
	v \mapsto FF^\ast v.
\]
For orthonormal bases---or more generally, and by definition, \emph{Parseval frames}---the frame operator is the identity and the synthesis operator provides a simple method for reconstructing the signal $v$ from the data $F^\ast v$.

As in the case of Parseval frames, it is often desirable to choose frames with a fixed spectrum of the frame operator, for example to provide optimal reconstruction in a given noise model~\cite{Goyal:2001cd,Casazza:2011ev,Viswanath:2002bv}. Likewise, the squared norms $\|f_1\|^2,\dots,\|f_N\|^2$ of the frame vectors are often fixed for both practical and theoretical reasons~\cite{Kovacevic:2007fja,Viswanath:1999hf,Rupf:1994fl,Casazza:2006cx}.

In other words, given vectors $\boldsymbol{\lambda}=(\lambda_1, \dots , \lambda_d)$ and $\boldsymbol{r}=(r_1, \dots , r_N)$ of positive numbers, we are often interested in selecting frames from the space $\F^{\mathfrak{H},N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ of frames $f_1, \dots, f_N \in \mathfrak{H}$ so that $\boldsymbol{\lambda}$ is the spectrum of the frame operator $FF^\ast$ and $\|f_i\|^2 = r_i$ for $i=1, \dots , N$. In particular, two natural questions immediately present themselves:
\begin{enumerate}
	\item Does there exist a frame for $\mathfrak{H}$ with prescribed data $\boldsymbol{\lambda}$ and $\boldsymbol{r}$? In other words, is $\F^{\mathfrak{H},N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ non-empty?
	\item Can we interpolate between arbitrary elements of $\F^{\mathfrak{H},N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$? In other words, is $\F^{\mathfrak{H},N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ path-connected?
\end{enumerate}

Question 1 has been completely answered by Casazza and Leon~\cite{Casazza:2010ti}, who give a simple compatibility criterion for $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ which determines whether or not $\F^{\mathfrak{H},N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is empty. Moreover, the same criterion applies in both the real and the complex cases.

Question 2 is a generalization of the well-known \emph{frame homotopy conjecture}, which was posed by Larson in a 2002 REU and first appeared in the literature in Dykema and Strawn's 2006 paper~\cite{Dykema:2006ux}. This conjecture says that when $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ are constant vectors, the space $\F^{\mathfrak{H},N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is path-connected (when considered with the natural subspace topology). Up to scale one can assume that the constant vector $\boldsymbol{r} = (1, \dots , 1)$, and hence the frames under consideration are \emph{unit-norm tight frames}, so the frame homotopy conjecture says that the space of unit-norm tight frames in $\mathfrak{H}$ is path-connected.

The frame homotopy conjecture was proved for both real and complex frames by Cahill, Mixon, and Strawn in a 2017 paper~\cite{Cahill:2017gv}. In previous work~\cite{NeedhamSGC}, we gave a complete answer to Question 2 for complex frames, showing that the space $\mathcal{F}^{\C^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is always path-connected. For real frames with nonconstant $\boldsymbol{\lambda}$ or $\boldsymbol{r}$, Question 2 remains open.

There has been a recent flourishing of interest in \emph{quaternionic frames}~\cite{iverson_note_2021,Waldron:2020tp,Cohn:2016bz,EtTaoui:2020kf,Waldron:2020ti,Khokulan:2017ic,Sharma:2019js,Virender:2020ua}; that is, frames in $\mathbb{H}^d$ for $d \geq 1$ and $\mathbb{H}$ being the 4-dimensional skew field of quaternions. Since the group $\operatorname{Sp}(1)$ of unit quaternions is isomorphic to $\operatorname{Spin}(3)$ (the universal cover of the rotation group $\operatorname{SO}(3)$), the quaternions can be interpreted as a cone over $\operatorname{Spin}(3)$, so they are well-suited to parameterizing rotations in $\mathbb{R}^3$~\cite{Hanson:2006tr}. Just as complex numbers consist of a magnitude and a phase, quaternions consist of a magnitude and a versor, which determines a particular 3D rotation. Consequently, vectors $v \in \mathbb{H}^d$ can be used to record both magnitude and orientation information, for example in framed curves like polymers or inflatable elastic rods~\cite{Cantarella:2013bla,Hanson:2012vb,Howard:2011fj,Needham:2017vd}. Moreover, quaternions have become an increasingly popular tool for representing signals in data-driven applications. For example, quaternion-valued signals have recently been used to encode RGB images \cite{ell2006hypercomplex,fletcher2017development}, measurements in industrial machinery \cite{yi2017quaternion}, and multicomponent seismic measurements \cite{zhao2020quaternion}; see also the special issue of \emph{Signal Processing} on Hypercomplex Signal Processing~\cite{SPspecialissue}. This motivates the extension of signal processing techniques for classical (i.e., real- or complex-valued) signals to quaternionic signals.

Despite the difficulties introduced by the non-commutativity of the quaternions, it is still possible to define frames in $\mathbb{H}^d$, the norms of the frame vectors, and the spectrum of the frame operator (see \Cref{sub:quaternionic_frames}). In other words, the space $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ of length-$N$ frames in $\mathbb{H}^d$ with fixed frame spectrum $\boldsymbol{\lambda}$ and fixed frame vector norms $\|f_i\|^2 = r_i$ exists, and the goal of this paper is to give complete answers to Questions 1 and 2 in this case. Specifically, we show:

\begin{thm}\label{thm:existence}
	Let $N$, $d$, $\boldsymbol{\lambda}$, and $\boldsymbol{r}$ be as above, and additionally assume that $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ are sorted in non-increasing order: $\lambda_1 \geq \dots \geq \lambda_d > 0$ and $r_1 \geq \dots \geq r_N > 0$. Then $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is non-empty if and only if
	\[
		\sum_{i=1}^k r_i \leq \sum_{i=1}^k \lambda_i \text{ for all } k=1, \dots , d \qquad \text{and} \qquad \sum_{i=1}^N r_i  = \sum_{i=1}^d \lambda_i.
	\]
\end{thm}

\begin{thm}\label{thm:connected}
	For any $\boldsymbol{r}$ and $\boldsymbol{\lambda}$, the space $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is path-connected.
\end{thm}

Since the empty set is trivially path-connected, the substance of \Cref{thm:connected} is that all of the non-empty $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ spaces are path-connected.

Note that the condition in \Cref{thm:existence} is exactly the same as that given by Casazza and Leon in the real and complex cases~\cite{Casazza:2010ti}. In other words, for given $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ either all three of $\mathcal{F}^{\mathbb{R}^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$, $\mathcal{F}^{\C^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$, and $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ are empty, or all three are non-empty. Likewise, the result in \Cref{thm:connected} is the same as in the complex case~\cite{NeedhamSGC}, so that $\mathcal{F}^{\C^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ and $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ are both path-connected. In particular, \Cref{thm:existence,thm:connected} imply that all spaces of quaternionic unit-norm tight frames are non-empty and path-connected, so the standard frame homotopy conjecture is true in the quaternionic setting.

Our strategy for proving the complex version of \Cref{thm:connected} in~\cite{NeedhamSGC} was based on symplectic geometry, which is extremely well-suited to the complex setting, but not to either the real or quaternionic settings. In this paper, we instead identify (equivalence classes of) frames with fixed frame spectrum $\boldsymbol{\lambda}$ with an adjoint orbit of the action of the symplectic group\footnote{The (compact) symplectic group is the quaternionic analog of the unitary group; see \Cref{sub:quaternions} for the definition.} $\operatorname{Sp}(N)$ on the space of $N \times N$ Hermitian matrices, which is the natural home of the Gram matrix $F^\ast F$ associated to a frame $F$. This is not special to quaternionic frames: one can make an analogous identification for real and complex frames, so we view this as a unifying perspective on all three classes of frames.

The advantage is that such adjoint orbits are extremely nice geometrically. In particular, they are motivating special cases of both Kostant's convexity theorem~\cite{Kostant:1973ff} and of the theory of isoparametric submanifolds~\cite[Chapter~6]{Palais:1988ks}, and the basic strategy for proving \Cref{thm:existence,thm:connected} is to apply those very general tools to these specific problems.

Specifically, since $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ are essentially the spectrum and diagonal entries, respectively, of the Gram matrix $F^\ast F$, the real and complex versions of \Cref{thm:existence} can be proved using the Schur--Horn theorem~\cite{Antezana:2007ci,Tropp:2005do}. Likewise, we will prove \Cref{thm:existence} using a quaternionic Schur--Horn theorem. This quaternionic version of Schur--Horn is an easy consequence of Kostant's convexity theorem, but its statement does not seem to be readily accessible in the literature (though see~\cite[Example~8.6]{kobert_spectrahedral_2021}). Since it may be of some independent interest, we give a statement and proof in \Cref{thm:quaternionic Schur-Horn}. On the other hand, adjoint orbits (and, more generally, isoparametric submanifolds) share many of the nice features of Hamiltonian manifolds~\cite{Terng:1986hg,Mare:2005eo}. In particular, in this setting we have access to a version of Atiyah's connectedness theorem~\cite{Atiyah:1982ih} (see \Cref{thm:Mare}), and with that tool in hand the proof of \Cref{thm:connected} goes much as it did in the complex case.

\section{Quaternions and Quaternionic Frames}
\label{sec:background}
\subsection{Quaternions}
\label{sub:quaternions}
In thinking about $\mathbb{H}^m$ as a vector space over $\mathbb{H}$, some care is required. First, recall that $\mathbb{H}$ is a skew field whose elements can be written as
\[
	a + b \mathbf{i} + c \mathbf{j} + d \mathbf{k}
\]
with multiplication given by the identities
\[
	\mathbf{i}^2 = \mathbf{j}^2 = \mathbf{k}^2 = \mathbf{i} \mathbf{j} \mathbf{k} = -1,
\]
which in particular implies $\mathbf{i} \mathbf{j} = \mathbf{k}$, $\mathbf{j} \mathbf{i} = -\mathbf{k}$, etc. The products of distinct elements of $\{\mathbf{i},\mathbf{j},\mathbf{k}\}$ are encoded in the diagram
\[
	\begin{tikzpicture}[decoration = {markings,
	    mark = at position 0.999 with {\arrow[>=stealth]{>}}
	  }
	  ]
	  \path (90 :1cm) node (i) {$\mathbf{i}$}
	        (330:1cm) node (j) {$\mathbf{j}$}
	        (210:1cm) node (k) {$\mathbf{k}$};
	  \draw[postaction = decorate] (i) to[bend left=45] (j);
	  \draw[postaction = decorate] (j) to[bend left=45] (k);
	  \draw[postaction = decorate] (k) to[bend left=45] (i);
	\end{tikzpicture}
\]
which for example says that $\mathbf{k} \mathbf{i} = \mathbf{j}$ since the product follows arrows with their correct orientation and $\mathbf{k} \mathbf{j} = -\mathbf{i}$ since this product reverses directions of arrows.

Instead of thinking in terms of quadruples of real numbers, it is also often convenient to interpret quaternions as pairs of complex numbers:
\[
	q = a+b\mathbf{i}+c\mathbf{j}+d\mathbf{k} = (a+b\mathbf{i}) + (c+d\mathbf{i})\mathbf{j}.
\]
This makes it easy to identify an element of $\mathbb{H}$ with a $2 \times 2$ complex matrix:
\begin{equation}\label{eq:H to C map}
	z + w\mathbf{j} \mapsto \begin{bmatrix} z & w \\ -\overline{w} & \overline{z} \end{bmatrix}.
\end{equation}
This map is an algebra isomorphism between $\mathbb{H}$ and the collection of $2 \times 2$ complex matrices of the given form.

The \emph{conjugate} of a quaternion $q =a + b\mathbf{i}+c\mathbf{j}+d \mathbf{k} \in \mathbb{H}$ is defined by
\[
	\overline{q} = a - b\mathbf{i}-c\mathbf{j}-d\mathbf{k},
\]
and the \emph{modulus} by
\[
	|q| = \sqrt{q \overline{q}} = \sqrt{a^2+b^2+c^2+d^2}.
\]

Since quaternionic multiplication is non-commutative, we have to distinguish between left and right vector spaces over $\mathbb{H}$. In this paper we will follow Waldron's conventions~\cite{Waldron:2020ti} and consider $\mathbb{H}^m$ as a right vector space over $\mathbb{H}$, meaning that scalar multiplication happens on the right. This is done so that left matrix multiplication is linear: for $v_1,\dots , v_N \in \mathbb{H}^d$ thought of as column vectors, $\alpha_1,\dots , \alpha_N \in \mathbb{H}$, and a matrix $A \in \mathbb{H}^{m \times d}$, we have
\[
	A(v_1 \alpha_1 + \dots + v_N \alpha_N) = (A v_1) \alpha_1 + \dots + (A v_N) \alpha_N.
\]

Much of linear algebra over $\mathbb{H}$ goes through as over $\mathbb{R}$ or $\C$; see~\cite{Zhang:1997cd,Rodman:2014ue} for more, but in particular recall that $\mathbb{H}^m$ has a standard $\mathbb{H}$-valued Hermitian inner product given by
\begin{equation}\label{eq:inner product}
	\langle v, w \rangle = \sum_i \overline{w}_j v_j \in \mathbb{H}
\end{equation}
for all $v,w \in \mathbb{H}^m$. The \emph{Frobenius inner product} on matrices $A,B \in \mathbb{H}^{m \times k}$ is given by
\[
	\langle A, B \rangle_F := \operatorname{tr}(B^\ast A),
\]
where, as usual, $B^\ast$ is the conjugate transpose of the matrix $B$. This agrees with the Hermitian inner product~\eqref{eq:inner product} on vectorized versions of the matrices.

The quaternionic analog of the unitary group is usually called the \emph{symplectic group} (sometimes the \emph{compact symplectic group}) and denoted $\operatorname{Sp}(m)$. That is, a $m \times m$ quaternionic matrix $U \in \mathbb{H}^{m \times m}$ is symplectic if
\[
	\langle U v, U w \rangle = \langle v , w \rangle
\]
for all $v,w \in \mathbb{H}^m$; equivalently, $U^\ast U = \mathbb{I}_m$, the $m \times m$ identity matrix. $\operatorname{Sp}(m)$ is a compact, semisimple Lie group with type-$C$ Lie algebra $\mathfrak{sp}(m)$ consisting of the skew-Hermitian $m \times m$ quaternionic matrices; that is, those matrices $A$ so that $A^\ast = -A$. A simple parameter count shows that
\[
	\dim_\mathbb{R} \operatorname{Sp}(m) = \dim_\mathbb{R} \mathfrak{sp}(m) = 3m+4\frac{m(m-1)}{2} = 2m^2+m.
\]

Finally, we introduce the notation
\[
	\mathcal{H}(m) = \{A \in \mathbb{H}^{m \times m} : A^\ast = A\}
\]
for the $(2m^2-m)$-dimensional real vector space of $m \times m$ Hermitian quaternionic matrices, and let $\mathcal{H}_0(m)$ be the subspace of traceless Hermitian matrices.

\subsection{Quaternionic Frames}
\label{sub:quaternionic_frames}
Frames in $\mathbb{H}^d$ are defined in the usual way~\cite{Waldron:2020ti,Khokulan:2017ic,Sharma:2019js,Virender:2020ua}: a sequence $(f_i)$ of vectors in $\mathbb{H}^d$ is a \emph{frame} if there exist $0<A\leq B < \infty$ so that
\begin{equation}\label{eq:frame inequality}
	A\|v\|^2 \leq \sum_i |\langle v, f_i\rangle|^2 \leq B \|v\|^2.
\end{equation}
For finite collections of vectors the upper bound is automatically satisfied with $B = \sum_i \|f_i\|^2$ and the lower bound is satisfied when the $f_i$ span $\mathbb{H}^d$. In other words, a finite collection of vectors $f_1, \dots , f_N \in \mathbb{H}^d$ is a frame for $\mathbb{H}^d$ if and only if $\{f_1, \dots , f_N\}$ is a spanning set for $\mathbb{H}^d$. We will use $\F^{\quat^d,N}$ to denote the collection of all frames consisting of $N$ vectors in $\mathbb{H}^d$.

A frame is called \emph{tight} if we can choose $A=B$ in \eqref{eq:frame inequality}, and \emph{Parseval} if $A=B=1$. We get a nice alternative characterization of these frames by introducing some operators. To do so, we will usually identify a frame $f_1, \dots , f_N$ with the $d \times N$ matrix $F = [f_1 | \dots | f_N ] \in \mathbb{H}^{d \times N}$ whose columns are the frame vectors, and we will often just write $F \in \F^{\quat^d,N}$. The \emph{analysis operator} associated to the frame is the map $\mathbb{H}^d \to \mathbb{H}^N$ given by
\[
	v \mapsto (\langle v, f_1 \rangle , \dots , \langle v, f_N\rangle ) = F^\ast v,
\]
and the \emph{synthesis operator} is the map $\mathbb{H}^N \to \mathbb{H}^d$ given by
\[
	w \mapsto \sum_{i=1}^N f_i w_i = F w.
\]

Composing these two operators one way gives the \emph{frame operator} $S:\mathbb{H}^d \to \mathbb{H}^d$ given by
\[
	S(v) = \sum_{i=1}^N f_i \langle v, f_i \rangle = FF^\ast v.
\]
In this interpretation, a tight frame has frame operator $S = FF^\ast = A \mathbb{I}_d$, so in particular a Parseval frame satisfies $FF^\ast = \mathbb{I}_d$. Composing in the other way gives the \emph{Gram matrix} $F^\ast F$ whose entries are the pairwise inner products of the frame vectors. For real and complex frames, cyclic invariance of the trace immediately implies that $\operatorname{tr}(FF^\ast) = \operatorname{tr}(F^\ast F)$, which is an important and frequently-used identity. Over the quaternions, the trace of a product of matrices is not generally invariant under cyclic permutations of terms, but the real part is: $\Re(\operatorname{tr}(AB)) = \Re(\operatorname{tr}(BA))$ for quaternionic matrices $A$ and $B$ of dimensions for which these products make sense. Since the frame operator and Gram matrix are both Hermitian, and hence have real eigenvalues and real trace, it follows that
\begin{equation}\label{eq:cyclic trace}
	\operatorname{tr}(FF^\ast) = \Re(\operatorname{tr}(FF^\ast)) = \Re(\operatorname{tr}(F^\ast F)) = \operatorname{tr}(F^\ast F) = \sum_{i=1}^N \|f_i\|^2.
\end{equation}

In fact, since there is a singular value decomposition for quaternionic matrices~\cite{Zhang:1997cd,Rodman:2014ue}, the usual argument from the real and complex cases shows that the frame operator $FF^\ast \in \mathcal{H}(d)$ and the Gram matrix $F^\ast F \in \mathcal{H}(N)$ have the same nonzero (right) eigenvalues. In particular:

\begin{prop}\label{prop:spectrum}
	If $F \in \F^{\quat^d,N}$ has frame operator $FF^\ast$ with spectrum $\boldsymbol{\lambda} = (\lambda_1, \dots , \lambda_d) \in \mathbb{R}_+^d$,\footnote{Since the columns of $F$ are a spanning set for $\mathbb{H}^d$, the frame operator $FF^\ast$ must be nonsingular, so we know $\lambda_i > 0$ for all $i=1,\dots, d$.} then the Gram matrix $F^\ast F$ has spectrum $\widetilde{\boldsymbol{\lambda}} := (\lambda_1, \dots , \lambda_d, 0 , \dots , 0)$.
\end{prop}

Since we are interested in frames with specified spectrum of the frame operator, we introduce some terminology and notation:

\begin{definition}\label{def:frame spectrum}
	If $F \in \F^{\quat^d,N}$, call the spectrum of $FF^\ast$ the \emph{frame spectrum} of $F$. For given $\boldsymbol{\lambda} = (\lambda_1, \dots , \lambda_d) \in \mathbb{R}_+^d$, let $\F^{\quat^d,N}_{\boldsymbol{\lambda}}$ be the set of all frames in $\F^{\quat^d,N}$ with frame spectrum equal to $\boldsymbol{\lambda}$.
\end{definition}

Since Gram matrices are positive semidefinite and the spectral decomposition of Hermitian quaternionic matrices works as expected~\cite{Zhang:1997cd,Rodman:2014ue}, the same proof as in the real or complex case shows that every positive semidefinite Hermitian quaternionic matrix is the Gram matrix of a collection of vectors which is a frame for its span, and that this collection is uniquely determined up to the action of the symplectic group. In other words, if $M \in \mathcal{H}(N)$ is positive semidefinite of rank $d$, then there exists $F = [f_1 | \dots | f_N] \in \F^{\quat^d,N} \subset \mathbb{H}^{d \times N}$ so that $M = F^\ast F$. Notice that, for $U \in \operatorname{Sp}(d)$, the frame $UF$ has the same Gram matrix: $(UF)^\ast (UF) = F^\ast U^\ast U F = F^\ast F = M$. Conversely, if $G \in \mathbb{H}^{d \times N}$ so that $G^\ast G = M$, then $G = UF$ for some $U \in \operatorname{Sp}(d)$.

We summarize the above discussion in the following proposition:

\begin{prop}\label{prop:equivalence classes}
	$\operatorname{Sp}(d)$-equivalence classes of frames in $\F^{\quat^d,N}$ are uniquely determined by their Gram matrices, which consist of all rank-$d$ positive semidefinite elements of $\mathcal{H}(N)$.
\end{prop}

\section{Adjoint Orbits, Convexity, and Isoparametric Submanifolds}
\label{sec:isoparametric}

Our next goal is to connect the story of Gram matrices to adjoint orbits and isoparametric submanifolds. We will largely follow Mare's discussion~\cite[Example~5.4]{Mare:2005eo}; see \cite[\S III.7]{Helgason:1978vb} for a more general exposition of Cartan decompositions.

Thinking of a square quaternionic matrix $A \in \mathbb{H}^{m \times m}$ as
\[
	A = Z + W \mathbf{j}
\]
for $Z,W \in \C^{m \times m}$, we can define a map $\Psi_m:\mathbb{H}^{m \times m} \to \C^{2m \times 2m}$ analogously to~\eqref{eq:H to C map}:
\[
	\Psi_m(Z+W\mathbf{j}) :=  \begin{bmatrix} Z & W \\ -\overline{W} & \overline{Z} \end{bmatrix}.
\]

If $A \in \mathbb{H}^{m \times m}$ is invertible, then $\Psi_m(A)$ is as well. Moreover, the image $\Psi_m(\operatorname{Sp}(m))$ of $\operatorname{Sp}(m)$ under this map is a subgroup of $\operatorname{SU}(2m)$, the group of $2m \times 2m$ unitary matrices with determinant 1, and in fact $\Psi_m(\operatorname{Sp}(m))$ is precisely the fixed point set inside $\operatorname{SU}(2m)$ of the involution $\sigma:\C^{2m \times 2m} \to \C^{2m \times 2m}$ given by
\[
	\sigma(M) = \Omega^\ast \overline{M} \Omega,
\]
where $\Omega = \begin{bmatrix} 0 & \mathbb{I}_m \\ -\mathbb{I}_m & 0 \end{bmatrix}$.

Passing to the Lie algebra gives the Cartan decomposition
\[
	\mathfrak{su}(2m) = \mathfrak{k} \oplus \mathfrak{p},
\]
where $\mathfrak{k}$ is the $(+1)$-eigenspace of the linearization $d_{\mathbb{I}_m}\sigma: \mathfrak{su}(2m) \to \mathfrak{su}(2m)$ and $\mathfrak{p}$ is the $(-1)$-eigenspace.\footnote{Since $d_{\mathbb{I}_m}\sigma(A) = -A^\ast$ is an involution, its only eigenvalues are $\pm 1$.} The subalgebra $\mathfrak{k}$ is just the Lie algebra of $\Psi_m(\operatorname{Sp}(m)) \simeq \operatorname{Sp}(m)$, so $\mathfrak{k} \simeq \mathfrak{sp}(m)$, with an explicit isomorphism given by $\Psi_m|_{\mathfrak{sp}(m)}: \mathfrak{sp}(m) \to \mathfrak{k}$.

Notice that $[\mathfrak{k},\mathfrak{p}] \subseteq \mathfrak{p}$, and so the adjoint action of $\Psi_m(\operatorname{Sp}(m))$ on $\mathfrak{su}(2m)$ restricts to an action on $\mathfrak{p}$. At first glance this appears slightly esoteric, but we can use $\Psi_m$ to relate $\mathfrak{p}$ to a more familiar collection of matrices. Specifically, $\Psi_m|_{\mathcal{H}_0(m)}$ gives an $\operatorname{Sp}(m)$-equivariant linear isomorphism between the traceless Hermitian matrices $\mathcal{H}_0(m)$ and $\mathfrak{p}$, and the adjoint action of $\Psi_m(\operatorname{Sp}(m))$ on $\mathfrak{p}$ corresponds to the conjugation action of $\operatorname{Sp}(m)$ on $\mathcal{H}_0(m)$.

Under this isomorphism, we can identify the standard maximal abelian subspace of $\mathfrak{p}$ with the set $\mathfrak{a} \subset \mathcal{H}_0(m)$ of real diagonal $m \times m$ matrices with trace $0$. The orbit $\mathcal{O}_\Lambda := \operatorname{Sp}(m)\cdot \Lambda$ of a point $\Lambda = \operatorname{diag}(\lambda_1, \dots , \lambda_m) \in \mathfrak{a}$ is simply
\[
	\mathcal{O}_\Lambda = \{U \Lambda U^\ast : U \in \operatorname{Sp}(m)\},
\]
which is precisely the collection of Hermitian $m \times m$ quaternionic matrices with spectrum $(\lambda_1, \dots , \lambda_m)$. At this point, it is hopefully clear that we intend to use \Cref{prop:spectrum,prop:equivalence classes} to relate $\operatorname{Sp}(d)$-equivalence classes of frames to orbits of this form. Before doing so, let's see what desirable features these orbits have.

Needless to say, we could have introduced these orbits without recourse to Cartan decompositions or adjoint actions. The point of approaching things in this slightly roundabout way, though, is that we can now easily see that $\mathcal{O}_\Lambda$ fits into both \emph{Kostant's convexity theorem} and the story of so-called \emph{isoparametric submanifolds}.

\subsection{Kostant's Convexity Theorem}
\label{sub:Kostant}

Consider the Cartan decomposition
\[
	\mathfrak{G} = \mathfrak{K} \oplus \mathfrak{P}
\]
of the Lie algebra $\mathfrak{G}$ of a semisimple Lie group $G$.\footnote{In defiance of the usual convention, we are using capital fraktur letters here for the Lie algebra and its subspaces, so as not to confuse the general $\mathfrak{K}$ and $\mathfrak{P}$ discussed here with the specific $\mathfrak{k}$ and $\mathfrak{p}$ defined above.} Then, as above, $\mathfrak{K}$ is the Lie algebra of a compact Lie subgroup $K \subset G$, and the adjoint action of $K$ on $\mathfrak{G}$ restricts to an action on $\mathfrak{P}$. Let $\mathfrak{A} \subset \mathfrak{P}$ be a maximal abelian subspace and let $P:\mathfrak{P} \to \mathfrak{A}$ be the orthogonal (with respect to the Killing form) projection. The \emph{Weyl group} $W$ associated with the pair $(\mathfrak{A},\mathfrak{G})$ is the finite group $N_K(\mathfrak{A})/Z_K(\mathfrak{A})$, where $N_K(\mathfrak{A})$ is the normalizer of $\mathfrak{A}$ in $K$, and $Z_K(\mathfrak{A})$ is its centralizer.

Now, let $a \in \mathfrak{A} \subset \mathfrak{P}$ and let $\mathcal{O}_a$ be the orbit of $a$ under the adjoint action of $K$ on $\mathfrak{P}$. Then Kostant's convexity theorem characterizes the image of $\mathcal{O}_a$ under $P$:

\begin{thm}[Kostant~\cite{Kostant:1973ff}]\label{thm:Kostant}
	$P(\mathcal{O}_a) = \operatorname{conv}(W\cdot a)$, the convex hull of the Weyl orbit of $a \in \mathfrak{A}$.
\end{thm}

In the case where $G=\operatorname{SL}_n(\C)$ and $K = \operatorname{SU}(n)$, this is essentially the Schur--Horn theorem~\cite{schur1923uber,horn1954doubly}, which says that the diagonal entries of Hermitian matrices with fixed spectrum fill out the convex hull of the collection of vectors given by all possible re-orderings of the spectrum.

With $G=\operatorname{SU}(2m)$, $K=\operatorname{Sp}(m)$, and the Cartan decomposition $\mathfrak{su}(2m) = \mathfrak{k} \oplus \mathfrak{p}$ as above, \Cref{thm:Kostant} implies the following quaternionic analog of the Schur--Horn theorem (compare to~\cite[Example~8.6]{kobert_spectrahedral_2021}):

\begin{thm}\label{thm:quaternionic Schur-Horn}
	Let $\boldsymbol{\lambda}= (\lambda_1 , \dots , \lambda_m) \in \mathbb{R}^m$. Let $\mathcal{H}_{\boldsymbol{\lambda}}(m)$ be the collection of quaternionic Hermitian $m \times m$ matrices with spectrum $\boldsymbol{\lambda}$. Let $\Delta: \mathcal{H}(m) \to \mathbb{R}^m$ record the diagonal entries of a matrix. Then
	\[
		\Delta(\mathcal{H}_{\boldsymbol{\lambda}}(m)) = \operatorname{conv}(S_m \cdot \boldsymbol{\lambda}),
	\]
	the convex hull of the permutation orbit of $\boldsymbol{\lambda}$.
\end{thm}

\begin{proof}
	First of all, if $\lambda_1 + \dots + \lambda_m = 0$, then $\mathcal{H}_{\boldsymbol{\lambda}}(m) \subset \mathcal{H}_0(m)$ corresponds to the orbit of $\Lambda = \operatorname{diag}(\lambda_1, \dots , \lambda_m)$ under the adjoint action of $\operatorname{Sp}(m)$ on $\mathfrak{p} \simeq \mathcal{H}_0(m)$. Moreover, the Weyl group associated to $(\mathfrak{a},\mathfrak{su}(2m))$ is the symmetric group $S_m$ and the orthogonal projection $\mathfrak{p} \to \mathfrak{a}$ corresponds to the diagonal entry map $\Delta$, so \Cref{thm:Kostant} says exactly that
	\[
		\Delta(\mathcal{H}_{\boldsymbol{\lambda}}(m)) = \operatorname{conv}(S_m \cdot \boldsymbol{\lambda}).
	\]

	If $\lambda_1 + \dots + \lambda_m \neq 0$, then we just need to re-center, which we do using the map $\tau: \mathcal{H}_{\boldsymbol{\lambda}}(m) \to \mathcal{O}_\Lambda$ defined in \Cref{lem:re-center} below. In the notation of that lemma, \Cref{thm:Kostant} implies that
	\[
		(\Delta \circ \tau)(\mathcal{H}_{\boldsymbol{\lambda}}(m)) = \operatorname{conv}(t(S_m \cdot \boldsymbol{\lambda})).
	\]

	Since taking the convex hull is $t$-equivariant and since $\Delta \circ \tau = t \circ \Delta$, it follows that $\Delta(\mathcal{H}_{\boldsymbol{\lambda}}(m)) = \operatorname{conv}(S_m \cdot \boldsymbol{\lambda})$, as desired.
\end{proof}

\begin{lem}\label{lem:re-center}
	Let $\boldsymbol{\lambda} = (\lambda_1, \dots , \lambda_m) \in \mathbb{R}^m$, let $\sigma = \lambda_1 + \dots + \lambda_m$, and define $\Lambda := \operatorname{diag}(\lambda_1, \dots , \lambda_m)-\frac{\sigma}{m} \mathbb{I}_m$. Then the map $B \mapsto B - \frac{\sigma}{m} \mathbb{I}_m$ defines a diffeomorphism $\tau: \mathcal{H}_{\boldsymbol{\lambda}}(m) \to \mathcal{O}_\Lambda$ which is equivariant with respect to the $\operatorname{Sp}(m)$ actions on domain and range and makes the following diagram commute:
	\[\begin{tikzcd}
		{\mathcal{H}_{\boldsymbol{\lambda}}(m)} && {\mathcal{O}_\Lambda} \\
		\\
		{\mathbb{R}^m} && {\mathbb{R}^m}
		\arrow["\tau", from=1-1, to=1-3]
		\arrow["\Delta"', from=1-1, to=3-1]
		\arrow["t", from=3-1, to=3-3]
		\arrow["\Delta", from=1-3, to=3-3]
	\end{tikzcd}\]
	Here $t: \mathbb{R}^m \to \mathbb{R}^m$ is the translation map $t(\boldsymbol{x}) := \boldsymbol{x}-\left(\frac{\sigma}{m}, \dots , \frac{\sigma}{m}\right)$.
\end{lem}

\begin{proof}
	First, $\tau$ is certainly well-defined and smooth with smooth inverse, so it is a diffeomorphism.

	If $B \in \mathcal{H}_{\boldsymbol{\lambda}}(m)$ and $U \in \operatorname{Sp}(d)$, then
	\[
		U \tau(B) U^\ast = U\left(B - \frac{\sigma}{m}\mathbb{I}_m\right)U^\ast = UBU^\ast - \frac{\sigma}{m}\mathbb{I}_m = \tau(UBU^\ast),
	\]
	so $\tau$ is $\operatorname{Sp}(m)$-equivariant.

	Finally, for $B \in \mathcal{H}_{\boldsymbol{\lambda}}(m)$,
	\[
		(\Delta \circ \tau)(B) = \Delta\left(B - \frac{\sigma}{m}\mathbb{I}_m\right) = \Delta(B) - \left(\frac{\sigma}{m},\dots , \frac{\sigma}{m}\right) = (t \circ \Delta)(B),
	\]
	so the diagram commutes.
\end{proof}

\subsection{Isoparametric Submanifolds}
\label{sub:isoparametric submanifolds}

Isoparametric hypersurfaces were first studied by Cartan~\cite{Cartan:1938kh,Cartan:1939gv,Cartan:1939wj,Cartan:tw,Nomizu:1975bq}; for a more modern and general introduction see \cite[Chapter~6]{Palais:1988ks} or~\cite{Terng:1985eu}. First, we give the definition, the details of which will not concern us greatly:

\begin{definition}\label{def:isoparametric}
	A submanifold $M$ of a Riemannian manifold $N$ is \emph{isoparametric} if its normal bundle $\nu(M)$ is flat and the principal curvatures along any parallel normal field of $M$ are constant.
\end{definition}

Principal orbits of isotropy representations are classic examples of isoparametric submanifolds. In the context of our story, we have:

\begin{prop}[{see, e.g.,~\cite[Example~6.5.6]{Palais:1988ks}}]\label{prop:orbits are isoparametric}
	If $\Lambda \in \mathfrak{a}$ is generic (i.e., $\lambda_i \neq \lambda_j$ for all $i \neq j$), then $\mathcal{O}_\Lambda$ is an isoparametric submanifold of $\mathcal{H}_0(N)$, thought of as a Riemannian manifold by taking the Frobenius inner product on each tangent space. Moreover, the normal space to $\mathcal{O}_\Lambda$ at $\Lambda$ is just $\mathfrak{a}$.

	If $\Lambda$ is not generic, then $\mathcal{O}_\Lambda$ is parallel to an isoparametric submanifold of $\mathcal{H}_0(N)$.
\end{prop}

To explain the terminology in the previous sentence, if $M \subset \mathbb{R}^n$ is isoparametric and $\xi$ is a parallel section of the normal bundle $\nu(M)$, then $M_\xi := \{p + \xi(p) : p \in M\}$ is a \emph{parallel submanifold} to $M$, and these parallel submanifolds give a singular foliation of the ambient space. In the special case of isotropy orbits, the orbit foliation and the parallel foliation coincide.

Isoparametric submanifolds and their parallels have some of the nice features of symplectic manifolds admitting Hamiltonian torus actions without necessarily being symplectic. For example, the following fundamental result in symplectic geometry has an analog in the isoparametric setting.

\begin{thm}[{Atiyah~\cite{Atiyah:1982ih} and Guillemin--Sternberg~\cite{Guillemin:1982gx}}]\label{thm:AGS}
	Let $(M,\omega)$ be a compact symplectic manifold admitting a Hamiltonian action of a torus $T\simeq \operatorname{U}(1)^n$. Let $\mu:M \to \mathbb{R}^n$ be the associated momentum map.
	\begin{itemize}
		\item For any $v \in \mathbb{R}^n$, $\mu^{-1}(v)$ is either empty or connected.
		\item $\mu(M)$ is the convex hull of the images of the fixed points of the $T$-action.
	\end{itemize}
\end{thm}

\begin{example}
	Let $K$ be a compact Lie group of rank $n$ with Lie algebra $\mathfrak{K}$. Let $\mathcal{O} \subset \mathfrak{K}^\ast$ be an orbit of the coadjoint action of $K$ on $\mathfrak{K}^\ast$. This action of $K$ on $\mathcal{O}$ is Hamiltonian with momentum map being the inclusion $\mathcal{O} \hookrightarrow \mathfrak{K}^\ast$~\cite[Example~5.3.11]{mcduff2017introduction}.

	If $T \subset K$ is a maximal torus, then $T \simeq U(1)^n$ and the coadjoint action of $T$ on $\mathcal{O}$ is also Hamiltonian, with moment map given by the restriction of the projection $\mathfrak{K}^\ast \to \mathfrak{T}^\ast \simeq \mathbb{R}^n$ induced by the inclusion $T \hookrightarrow K$~\cite[Proposition~II.1.10]{Audin:2004bh}. Then \Cref{thm:AGS} implies that the image of this map is convex and its non-empty level sets are connected.

	In this case, convexity actually follows from \Cref{thm:Kostant}: let $G = K_\C$ be the complexification of $K$, with Lie algebra $\mathfrak{G} = \mathfrak{K} \oplus \mathbf{i} \mathfrak{K}$. This is a Cartan decomposition and the dual Lie algebra $\mathfrak{K}^\ast$ can be identified with the complementary subspace $\mathbf{i} \mathfrak{K}$ so that the coadjoint action of $K$ on $\mathfrak{K}^\ast$ corresponds to the adjoint action of $K$ on $\mathbf{i} \mathfrak{K}$ and $\mathfrak{T}^\ast \subset \mathfrak{K}^\ast$ corresponds to a maximal abelian subalgebra $\mathfrak{A} \subset \mathbf{i} \mathfrak{K}$. Hence, the projection map that \Cref{thm:Kostant} says has convex image is just the momentum map of the torus action. Of course, Kostant's convexity theorem was a major inspiration for Atiyah and Guillemin--Sternberg, whose result can be interpreted as generalizing this case of Kostant's theorem to arbitrary Hamiltonian torus actions.

	On the other hand, if $\mathcal{O}$ is a principal orbit, then $\mathcal{O}$ is an isoparametric submanifold of $\mathfrak{K}^\ast$ and the normal space at a point can be identified with $\mathfrak{A} \simeq \mathfrak{T}^\ast$. This gives yet another interpretation of the momentum map as the orthogonal projection onto the normal space at a point. Terng showed that, under this interpretation, the convexity part of \Cref{thm:AGS} generalizes to arbitrary isoparametric submanifolds:
\end{example}

\begin{thm}[{Terng~\cite{Terng:1986hg}}]\label{thm:Terng convexity}
	Let $M_\xi \subset \mathbb{R}^{n}$ be parallel to an isoparametric submanifold $M$. Let $p \in M$ so that $p + \xi(p) \in M_\xi$, let $\nu_p(M)$ be the normal space to $M$ at $p$, and let $P:M_\xi \to \nu_p(M)$ be orthogonal projection. Then $P(M_\xi)$ is a convex polytope.
\end{thm}

Given this, the projection map $P$ is like a momentum map for isoparametric submanifolds and their parallels, even though these submanifolds need not be symplectic. It is then reasonable to ask whether, as in the first part of \Cref{thm:AGS}, the level sets of $P$ are connected. They needn't be in general, but Mare gave a sufficient condition for all non-empty level sets of $P$ to be connected:

\begin{thm}[{Mare~\cite[Theorem~1.2 and Remark~1.3(a)]{Mare:2005eo}}]\label{thm:Mare}
	Let $M \subset \mathbb{R}^n$ be an isoparametric submanifold with all multiplicities $\geq 2$ and let $M_\xi$ be parallel to $M$. If $p \in M$ and $b \in \nu_p(M)$ is in the image of the projection $P: M_\xi \to \nu_p(M)$, then $P^{-1}(b)$ is connected.
\end{thm}

See Mare's paper for a general definition of multiplicities; in the setting of \Cref{thm:Kostant}, where we are considering orbits of the adjoint action of $K \subset G$ on the complementary subspace $\mathfrak{P} \subset \mathfrak{G}$, the multiplicities are the differences in dimension between a principal orbit and subprincipal orbits~\cite{Palais:1988ks,Hsiang:1988ka}. Recall from \Cref{prop:orbits are isoparametric} that the normal space can be identified with the maximal abelian subspace $\mathfrak{A}$.

\section{Existence}
\label{sec:admissibility}

Let $N \geq d \geq 1$ be integers and let $\boldsymbol{\lambda}=(\lambda_1, \dots, \lambda_d)$ and $\boldsymbol{r} = (r_1, \dots , r_N)$ be lists of positive numbers. Our main object of interest is the collection $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ of frames $F =[f_1 | \dots | f_N] \in \F^{\quat^d,N} \subset \mathbb{H}^{d \times N}$ so that $FF^\ast$ has spectrum $\boldsymbol{\lambda}$ and $\|f_i\|^2 = r_i$ for $i=1, \dots , N$. Our first goal is to answer the question: are there any such frames?

The analogous question for real and complex frames was answered by Casazza and Leon~\cite{Casazza:2010ti}, and the answer in the quaternionic case is essentially the same:

\begin{thm}\label{thm:admissibility}
	Let $N$, $d$, $\boldsymbol{\lambda}$, and $\boldsymbol{r}$ be as above. The space $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is non-empty if and only if $\boldsymbol{r} \in \operatorname{conv}(S_N \cdot \widetilde{\boldsymbol{\lambda}})$, the convex hull of the collection of vectors in $\mathbb{R}^N$ given by permuting the entries of $\widetilde{\boldsymbol{\lambda}} := (\lambda_1, \dots , \lambda_d,0,\dots,0)$.
\end{thm}

\begin{proof}
	\Cref{thm:quaternionic Schur-Horn} will be the key to the proof, so the goal is to relate $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ to a space of the form $\mathcal{H}_{\boldsymbol{\eta}}(m)$ for some choice of $\boldsymbol{\eta}$ and $m$.

	Recall from \Cref{prop:equivalence classes} that $\operatorname{Sp}(d)$-equivalence classes of frames in $\F^{\quat^d,N}$ are determined by their Gram matrices. Specifically, the space $\F^{\quat^d,N}_{\boldsymbol{\lambda}}/\operatorname{Sp}(d)$ of $\operatorname{Sp}(d)$-equivalence classes of frames with frame spectrum $\boldsymbol{\lambda}$ is diffeomorphic to the space of all possible Gram matrices of frames in $\F^{\quat^d,N}_{\boldsymbol{\lambda}}$.

	Now, if $F \in \F^{\quat^d,N}_{\boldsymbol{\lambda}}$, then definitionally $FF^\ast$ has spectrum $\boldsymbol{\lambda}$, and hence the Gram matrix $F^\ast F$ has spectrum $\widetilde{\boldsymbol{\lambda}}=(\lambda_1, \dots , \lambda_d, 0, \dots , 0) \in \mathbb{R}^N$. Conversely, any matrix in $\mathcal{H}_{\widetilde{\boldsymbol{\lambda}}}(N)$ can be realized as the Gram matrix of a frame in $\F^{\quat^d,N}_{\boldsymbol{\lambda}}$, so we see that
	\[
		\F^{\quat^d,N}_{\boldsymbol{\lambda}}/\operatorname{Sp}(d) \simeq \mathcal{H}_{\widetilde{\boldsymbol{\lambda}}}(N).
	\]

	In turn, if $\Delta:\mathcal{H}_{\widetilde{\boldsymbol{\lambda}}}(N) \to \mathbb{R}^N$ is the map which records diagonal entries, then the level set $\Delta^{-1}(\boldsymbol{r})$ is the collection of $\operatorname{Sp}(d)$-equivalence classes of frames with frame spectrum $\boldsymbol{\lambda}$ and squared frame norms $\boldsymbol{r}$; that is,
	\[
		\Delta^{-1}(\boldsymbol{r}) \simeq \F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})/\operatorname{Sp}(d).
	\]

	This means that $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is non-empty if and only if $\boldsymbol{r}$ is in the image of $\Delta$. But now \Cref{thm:quaternionic Schur-Horn} tells us that the image of $\Delta$ is exactly the convex hull of $S_N \cdot \widetilde{\boldsymbol{\lambda}}$, as desired.
\end{proof}

Since permuting entries doesn't change their sum, $\operatorname{conv}(S_N \cdot \widetilde{\boldsymbol{\lambda}})$ lies in the affine hyperplane $\{(x_1, \dots , x_N) : \sum_{i=1}^N x_i = \sum_{i=1}^d \lambda_i\}$. Hence, $\boldsymbol{r} \in \operatorname{conv}(S_N \cdot \widetilde{\boldsymbol{\lambda}})$ only if
\begin{equation}\label{eq:admissibility equality}
	\sum_{i=1}^N r_i = \sum_{i=1}^d \lambda_i.
\end{equation}
This is just a restatement of~\eqref{eq:cyclic trace}, which said the frame operator and the Gram matrix have the same trace.

In practice, the order of the frame vectors is no more than a bookkeeping convenience, so it is no problem to permute the frame vectors. Likewise, we can freely permute the numbers comprising the frame spectrum. In particular, we can get a more straightforward criterion for non-emptiness of $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ by sorting $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ in non-increasing order. Specifically, if we make the additional assumption that $\lambda_1 \geq \dots \geq \lambda_d > 0$ and $r_1 \geq \dots \geq r_N > 0$, then we see that $\boldsymbol{r} \in \operatorname{conv}(S_N \cdot \widetilde{\boldsymbol{\lambda}})$ if and only if it satisfies~\eqref{eq:admissibility equality} and
\begin{equation}\label{eq:admissibility inequality}
	\sum_{i=1}^k r_i \leq \sum_{i=1}^k \lambda_i \quad \text{for all }k=1, \dots , d.
\end{equation}

We call the (sorted) $\boldsymbol{r}$ satisfying~\eqref{eq:admissibility equality} and~\eqref{eq:admissibility inequality} \emph{$\boldsymbol{\lambda}$-admissible}. Notice that the admissibility criterion is exactly the same as that given by Casazza and Leon in the real and complex cases~\cite{Casazza:2010ti}. Thus, \Cref{thm:admissibility} is equivalent to \Cref{thm:existence}, which we restate in a more compact form:

\begin{existence}
	Let $N$, $d$, $\boldsymbol{\lambda}$, and $\boldsymbol{r}$ be as above, so that $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ are sorted in non-increasing order. Then $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is non-empty if and only if $\boldsymbol{r}$ is $\boldsymbol{\lambda}$-admissible.
\end{existence}

\section{Connectedness}
\label{sec:connectedness}

Next, we turn to the question of when $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is connected. As in the previous section, we will focus on $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})/\operatorname{Sp}(d)$, since the following lemma combined with the fact that $\operatorname{Sp}(d)$ is connected implies the quotient is connected if and only if $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is connected.

\begin{lem}\label{lem:connected quotient}
	Let $X$ be a topological space and let $G$ be a connected topological group acting continuously on $X$. If $X/G$ is connected, then $X$ is connected.
\end{lem}

This lemma follows from a standard point-set topology argument (see, e.g., \cite[Exercise~5.5]{Manetti:2015ep}) since connectedness of $G$ implies the fibers of the quotient map $X \to X/G$ are connected.

The strategy is to prove connectedness using \Cref{thm:Mare}, which applies to adjoint orbits of a compact group acting on the complementary subspace of the Lie algebra of some larger group $G$. Following the setup in \Cref{sec:isoparametric}, let $G=\operatorname{SU}(2N)$, so that $K=\operatorname{Sp}(N)$ and the complementary subspace $\mathfrak{p} \subset \mathfrak{su}(2N)$ can be identified with $\mathcal{H}_0(N)$, the space of traceless Hermitian quaternionic $N \times N$ matrices, and the standard maximal abelian subspace of $\mathfrak{p}$ corresponds to the subset $\mathfrak{a} \subset \mathcal{H}_0(N)$ of real diagonal $N \times N$ matrices with trace 0.

We saw in the proof of \Cref{thm:admissibility} that $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})/\operatorname{Sp}(d) \simeq \mathcal{H}_{\widetilde{\boldsymbol{\lambda}}}(N)$, the space of quaternionic Hermitian $N \times N$ matrices with spectrum $\widetilde{\boldsymbol{\lambda}} = (\lambda_1, \dots , \lambda_d, 0, \dots 0)$. Since $\lambda_1 + \dots + \lambda_d \neq 0$, this is not quite an orbit of the adjoint action of $\operatorname{Sp}(N)$ on $\mathcal{H}_0(N)$, but it is a simple translation of such an orbit.

Letting  $\sigma:=\lambda_1 + \dots + \lambda_d$ and
\[
	\widetilde{\Lambda} := \operatorname{diag}\left(\lambda_1 ,\dots , \lambda_d, 0, \dots, 0\right) - \frac{\sigma}{N} \mathbb{I}_N \in \mathcal{H}_0(N),
\]
define $\tau: \mathcal{H}_{\widetilde{\boldsymbol{\lambda}}}(N) \to \mathcal{O}_{\widetilde{\Lambda}}$ by $\tau(B) = B - \frac{\sigma}{N} \mathbb{I}_N$.

If $\boldsymbol{r}$ is $\boldsymbol{\lambda}$-admissible, then $\Delta^{-1}(\boldsymbol{r}) \subset \mathcal{H}_{\boldsymbol{\lambda}}(N)$ is non-empty, and hence so is
\[
	\tau(\Delta^{-1}(\boldsymbol{r})) = \Delta^{-1}(t(\boldsymbol{r})) \subset \mathcal{O}_{\widetilde{\Lambda}},
\]
where $t: \mathbb{R}^N \to \mathbb{R}^N$ is defined by $t(\boldsymbol{x}) := \boldsymbol{x} - \left(\frac{\sigma}{N},\dots , \frac{\sigma}{N}\right)$ and the above equality follows from the fact that the diagram in \Cref{lem:re-center} commutes. But now
\[
	\Delta^{-1}(t(\boldsymbol{r})) = P^{-1}(\operatorname{diag}(t(\boldsymbol{r}))),
\]
where $P: \mathcal{O}_{\widetilde{\Lambda}} \to \mathfrak{a}$ is the orthogonal projection. Since $\mathcal{O}_{\widetilde{\Lambda}}$ is an isotropy orbit\footnote{In fact, it is an example of a quaternionic flag manifold, and much is known about its cohomology~\cite{Mare:2008iz,Hsiang:1988ka,Mare:2006bf}.} and hence also a parallel submanifold to some principal orbit, and since all multiplicities in this setting are equal to 4~\cite[\S 3.3]{Hsiang:1988ka}, connectedness of $\Delta^{-1}(t(\boldsymbol{r}))$ follows from \Cref{thm:Mare}. Since $\tau$ is a diffeomorphism, we conclude that
\[
	\Delta^{-1}(\boldsymbol{r}) \simeq \F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})/\operatorname{Sp}(d)
\]
is connected.

Finally, applying \Cref{lem:connected quotient} shows that $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is connected whenever it is non-empty; i.e., whenever $\boldsymbol{r}$ is $\boldsymbol{\lambda}$-admissible. In turn, since $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is a real algebraic set in $\mathbb{H}^{d \times N} \simeq \mathbb{R}^{4dN}$, it is locally path-connected (in fact, triangulable by \L{}ojasiewicz's triangulation theorem \cite{lojasiewicz1964triangulation}) so that connectivity implies path-connectivity.

Thus, we have proved \Cref{thm:connected}, which we now restate:

\begin{connectedness}
	For any $\boldsymbol{r}$ and $\boldsymbol{\lambda}$, the space $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is path-connected.
\end{connectedness}

If we prefer to focus on the space
\[
	\F^{\quat^d,N}_S(\boldsymbol{r}) := \{F \in \F^{\quat^d,N} : FF^\ast = S \text{ and } \|f_i\|^2 = r_i \text{ for all } i = 1 , \dots , N\}
\]
of frames with fixed frame operator (rather than fixed frame spectrum) and fixed frame vector norms, connectivity still holds:

\begin{cor}\label{cor:fixed frame operator}
	If $S \in \mathcal{H}(N)$ is positive definite with spectrum $\boldsymbol{\lambda}$ and $\boldsymbol{r} \in \mathbb{R}_+^N$, the space $\F^{\quat^d,N}_S(\boldsymbol{r})$ is (i) non-empty if and only if $\boldsymbol{r}$ is $\boldsymbol{\lambda}$-admissible and (ii) path-connected.
\end{cor}

\begin{proof}
	If $\F^{\quat^d,N}_S(\boldsymbol{r}) \subset \F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is non-empty, then certainly $\boldsymbol{r}$ is $\boldsymbol{\lambda}$-admissible by \Cref{thm:admissibility}. Conversely, if $\boldsymbol{r}$ is $\boldsymbol{\lambda}$-admissible, then $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is non-empty, so that there exists $F \in \F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$, and $FF^\ast = V S V^\ast$ for some $V \in \operatorname{Sp}(d)$. But then $V^\ast F \in \F^{\quat^d,N}_S(\boldsymbol{r})$, so we see that $\F^{\quat^d,N}_S(\boldsymbol{r})$ is non-empty as well.

	Since the empty set is trivially path-connected, the only thing to prove is that $\F^{\quat^d,N}_S(\boldsymbol{r})$ is path-connected whenever it is non-empty. In this case, we know from \Cref{thm:connected} that $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ is path-connected. So if $F_0, F_1 \in \F^{\quat^d,N}_S(\boldsymbol{r}) \subset \F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$, then there exists a path $\widetilde{\gamma}: [0,1] \to \F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ so that $\widetilde{\gamma}(0) = F_0$ and $\widetilde{\gamma}(1) = F_1$.

	Of course, the path $\widetilde{\gamma}$ need not stay in $\F^{\quat^d,N}_S(\boldsymbol{r})$, but we can easily fix it up to do so. Indeed, the frame operator
	\[
		\widetilde{\gamma}(t)\widetilde{\gamma}(t)^\ast = U_t S U_t^\ast
	\]
	determines a continuous path $\{U_t: t \in [0,1]\} \subset \operatorname{Sp}(d)$, and hence
	\[
		\gamma(t) := U_t^\ast \widetilde{\gamma}(t)
	\]
	is a path in $\F^{\quat^d,N}_S(\boldsymbol{r})$ connecting $F_0$ and $F_1$.
\end{proof}

\section{Discussion}
\Cref{cor:fixed frame operator} is the exact analog of the main theorem in our previous paper~\cite{NeedhamSGC} on complex frames. While our argument in the complex case was based on symplectic geometry, the strategy employed here for proving \Cref{thm:connected} and \Cref{cor:fixed frame operator} can be adapted to the complex setting, so we view this approach as somewhat more general.

Indeed, much of the setup goes through in the real case as well: orthogonal equivalence classes of real frames with fixed frame spectrum correspond to adjoint orbits by way of their Gram matrices, and hence are parallel to isoparametric submanifolds of the space of symmetric matrices. Unfortunately, all multiplicities in this case are equal to 1~\cite[\S 3.2]{Hsiang:1988ka}, so \Cref{thm:Mare} does not apply. Indeed, while some of the real frame spaces $\mathcal{F}^{\mathbb{R}^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$ are connected~\cite{Cahill:2017gv}, others are not~\cite{Kapovich:1995wg} (compare with~\cite[Theorem~2.7]{Goyal:2001cd}), and it seems challenging (but interesting!) to characterize which $\boldsymbol{\lambda}$ and $\boldsymbol{r}$ lead to which outcome. Since the language of isoparametric submanifolds provides a common framework for understanding real, complex, and quaternionic frames, this perspective seems ripe for further exploration.

\Cref{thm:existence,thm:connected} have the same statements as the corresponding results in the complex case, which gives some reason to hope that direct translations of other results about complex frames might give true statements about quaternionic frames. For example, it seems likely that there is a quaternionic extension of the eigenstep method~\cite{Cahill:2013jv} which would construct all elements of $\F^{\quat^d,N}_{\boldsymbol{\lambda}}(\boldsymbol{r})$, and more generally give a constructive proof of \Cref{thm:quaternionic Schur-Horn} (cf.~\cite{Fickus:2013ki}).

\subsection*{Acknowledgments}

We are very grateful to Augustin-Liviu Mare, Emily King, and Colin Roberts for providing inspiration and sharing their knowledge and insight. This work was partially supported by grants from the National Science Foundation (DMS--2107808, Tom Needham; DMS--2107700, Clayton Shonkwiler).

Functional Analysis (math.FA); Differential Geometry (math.DG)


On generalization of D'Aurizio-Sándor trigonometric inequalities with a parameter

Abstract: In this work, we generalize the D'Aurizio-Sándor inequalities (\cite{D'Aurizio,Sandor}) using an elementary approach. In particular, our approach provides an alternative proof of the D'Aurizio-Sándor inequalities. Moreover, as an immediate consequence of the generalized D'Aurizio-Sándor inequalities, we establish the D'Aurizio-Sándor-type inequalities for hyperbolic functions.

\section{Introduction}

Based on infinite product expansions and inequalities on series and the Riemann's zeta function, D'Aurizio (\cite{D'Aurizio}) proved the following inequality:
\begin{equation}\label{eqn: D'Aurizio ineq}
\frac{1-\displaystyle\frac{\cos x}{\cos \frac{x}{2}}}{x^2}<\frac{4}{\pi^2},
\end{equation}
where $x\in(0,\pi/2)$. Using an elementary approach, S\'andor (\cite{Sandor}) offered an alternative proof of \eqref{eqn: D'Aurizio ineq} by employing trigonometric inequalities and an auxiliary function. In the same paper, S\'andor also provided the converse to \eqref{eqn: D'Aurizio ineq}:
\begin{equation}
\frac{1-\displaystyle\frac{\cos x}{\cos \frac{x}{2}}}{x^2}>\frac{3}{8},
\end{equation}
where $x\in(0,\pi/2)$. In addition, S\'andor found the following analogous inequality \eqref{eqn: Sandor ineq}  holds true for the case of sine functions:

\begin{theorem}[D'Aurizio-S\'andor inequalities (\cite{D'Aurizio,Sandor})]\label{thm: Daurizio-Sandor ineq original}
The two double inequalities
\begin{equation}
\frac{3}{8}<\frac{1-\displaystyle\frac{\cos x}{\cos \frac{x}{2}}}{x^2}<\frac{4}{\pi^2}
\end{equation}
and
\begin{equation}\label{eqn: Sandor ineq}
\frac{4}{\pi^2}(2-\sqrt{2})<\frac{\displaystyle 2-\frac{\sin x}{\sin \frac{x}{2}}}{x^2}<\frac{1}{4}
\end{equation}
hold for any $x\in(0,\pi/2)$.
\end{theorem}

Throughout this paper, we denote $\frac{1-\frac{\cos x}{\cos \frac{x}{p}}}{x^2}$ and $\frac{p-\frac{\sin x}{\sin \frac{x}{p}}}{x^2}$ by $f_c(x)$ and $f_s(x)$, respectively:

\begin{align}
f_{p}^{c}(x)= & \frac{1-\displaystyle\frac{\cos x}{\cos \frac{x}{p}}}{x^2}, \\
f_{p}^{s}(x)= & \frac{p-\displaystyle\frac{\sin x}{\sin \frac{x}{p}}}{x^2}.
\end{align}

Our aim is to generalize the D'Aurizio-S\'andor inequalities for the case of $f_{p}^{c}(x)$ and $f_{p}^{s}(x)$ as follows:

\begin{theorem}[Generalized D'Aurizio-S\'andor inequalities]\label{thm: Daurizio-Sandor ineq generalized}
Let $0<x<\pi/2$. Then the two double inequalities
\begin{equation}\label{eqn: cos ineq p>2}
\frac{4}{\pi^2}<\frac{1-\displaystyle\frac{\cos x}{\cos \frac{x}{p}}}{x^2}<\frac{p^2-1}{2\,p^2}
\end{equation}
and
\begin{equation}\label{eqn: sin ineq p>2}
\frac{4}{\pi^2}\left(p-\csc \left(\frac{\pi}{2p}\right)\right)<\frac{\displaystyle p-\frac{\sin x}{\sin \frac{x}{p}}}{x^2}<\frac{p^2-1}{6\,p}
\end{equation}
hold for $p=3, 4 , 5,\cdots$.  In particular, the double inequality \eqref{eqn: sin ineq p>2} remains true when $p=2$ while the double inequality \eqref{eqn: cos ineq p>2} is reversed when $p=2$.
\end{theorem}

The remainder of this paper is organized as follows. Section~\ref{sec: proof of the main results} is devoted to the proof of Theorem~\ref{thm: Daurizio-Sandor ineq generalized} and an alternative proof of Theorem~\ref{thm: Daurizio-Sandor ineq original}. In Section~\ref{sec: D'Aurizio-Sandor inequalities for hyperbolic functions}, we establish analogue of Theorem~\ref{thm: Daurizio-Sandor ineq generalized} for hyperbolic functions. As an application of Theorem~\ref{thm: Daurizio-Sandor ineq generalized}, we apply in Section~\ref{sec: application} inequality \eqref{eqn: sin ineq p>2} to the Chebyshev polynomials of the second kind and establish a trigonometric inequality.

\section{Proof of the main results}\label{sec: proof of the main results}

At first we will prove the following lemma. The lemma provides expressions of the higher-order derivative  $\frac{d^2}{dx^2}(x^3\frac{d}{dx} f_{p}^{\triangle}(x))$ involving $f_{p}^{\triangle}(x)$ $(\triangle=c, s)$, which are helpful in proving Theorem~\ref{thm: Daurizio-Sandor ineq generalized}. We note that the sign of $\frac{d^2}{dx^2}(x^3\frac{d}{dx} f_{p}^{\triangle}(x))$ plays a crucial role in proving Theorem~\ref{thm: Daurizio-Sandor ineq generalized}.

\begin{lemma}\label{lem: lots of f'''}
Let $0<x<\pi/2$ and $k=1,2,3,\cdots$. Then when $p\in\mathbb{R}$ and $p\neq0$, we hav
\begin{description}
  \item[$(i)$]
\begin{align}\label{eqn: f''' cos} \nonumber
\frac{d^2}{dx^2}& \left(x^3\frac{d}{dx} f_{p}^{c}(x)\right)
=-\frac{x\,\csc ^4\left(\frac{x}{p}\right)}{8\,p^3}
\bigg(
\left(p+1\right)^3
   \sin \Big(x-\frac{3 x}{p}\Big)+\left(p-1\right)^3 \sin
   \Big(x+\frac{3 x}{p}\Big)\\
+& \left(3 p^3+3 p^2-15 p-23\right)
   \sin \Big(x-\frac{x}{p}\Big)+\left(3 p^3-3 p^2-15 p+23\right)
   \sin\Big(x+\frac{x}{p}\Big)\bigg);
\end{align}
  \item[$(ii)$]
\begin{align}\label{eqn: f''' sin} \nonumber
\frac{d^2}{dx^2}& \left(x^3\frac{d}{dx} f_{p}^{s}(x)\right)
=\frac{x\,\csc ^4\left(\frac{x}{p}\right)}{8\,p^3}
\bigg(
\left(p+1\right)^3
   \sin \Big(x-\frac{3 x}{p}\Big)-\left(p-1\right)^3 \sin
   \Big(x+\frac{3 x}{p}\Big)\\
+& \left(-3 p^3-3 p^2+15 p+23\right)
   \sin \Big(x-\frac{x}{p}\Big)+\left(3 p^3-3 p^2-15 p+23\right)
   \sin\Big(x+\frac{x}{p}\Big)\bigg).
\end{align}
In particular,
  \item[$(iii)$] when $p=2\,k$,
\begin{equation}\label{eqn: 2k sin}
\frac{d^2}{dx^2} \left(x^3\frac{d}{dx} f_{p}^{s}(x)\right)=-\frac{x}{4\,k^3}\sum_{j=0}^{k-1} (2\,j+1)^3\,\sin\left(\frac{2\,j+1}{2\,k}\,x\right);
\end{equation}
  \item[$(iv)$] when $p=2\,k+1$,
\begin{align}
\label{eqn: 2k+1 cos}
\frac{d^2}{dx^2} \left(x^3\frac{d}{dx} f_{p}^{c}(x)\right)=&-\frac{16\,x}{(2\,k+1)^3}\sum_{j=1}^{k} j^3\,\sin\left(\frac{2\,j}{2\,k+1}\,x\right)\,(-1)^{j-1}, \\
\label{eqn: 2k+1 sin}
\frac{d^2}{dx^2} \left(x^3\frac{d}{dx} f_{p}^{s}(x)\right)=&-\frac{16\,x}{(2\,k+1)^3}\sum_{j=1}^{k} j^3\,\sin\left(\frac{2\,j}{2\,k+1}\,x\right).
\end{align}
 \item[$(v)$] For $\triangle=c, s$ and $p\in\mathbb{R}\setminus\{0\}$,
\begin{equation}
 \displaystyle \lim_{x\to0} \frac{d}{dx}\left(x^3\frac{d}{dx} f_{p}^{\triangle}(x)\right)=\lim_{x\to0} x^3\frac{d}{dx} f_{p}^{\triangle}(x)=0.
\end{equation}
\end{description}
\end{lemma}

\begin{proof}
$(i)$, $(ii)$ and $(v)$ follows directly from calculations using elementary Calculus. In particular, trigonometric addition formulas are used in proving $(i)$ and $(ii)$. To prove \eqref{eqn: 2k sin}, we claim
\begin{equation
-\frac{x}{4\,k^3}\sum_{j=0}^{k-1} (2\,j+1)^3\,\sin\left(\frac{2\,j+1}{2\,k}\,x\right)=-\,x\,\frac{d^3}{dx^3} \left( \frac{\sin x}{\sin\left(\frac{x}{2\,k}\right)}\right).
\end{equation}
Indeed, we rewrite
\begin{equation}\label{eqn: trans f''' p=2k}
\frac{1}{4\,k^3}\sum_{j=0}^{k-1} (2\,j+1)^3\,\sin\left(\frac{2\,j+1}{2\,k}\,x\right)
= 2\frac{d^3}{dx^3} \left(\sum_{j=0}^{k-1}\cos\left(\frac{2\,j+1}{2\,k}\,x\right)\right).
\end{equation}
On the other hand, making use of Euler's formula $e^{i\,z}=\cos z+i\,\sin z$ leads to an alternative expression of the left-hand side of \eqref{eqn: trans f''' p=2k}:
\begin{align}
\label{}
\sum_{j=0}^{k-1}\cos\left(\frac{2\,j+1}{2\,k}\,x\right)
&= \sum_{j=0}^{k-1} \displaystyle\Re\left\{ e^{i(\frac{x}{2\,k}+\frac{x}{k}\,j)}\right\}
  = \Re\left\{e^{i\frac{x}{2\,k}}\sum_{j=0}^{k-1} \displaystyle \left(e^{i\frac{x}{k}}\right)^j\right\}   \\
&= \Re\left\{e^{i\frac{x}{2\,k}}\displaystyle\frac{1-e^{i\,x}}{1-e^{i\frac{x}{k}}} \right\}
  = \Re\left\{e^{i\frac{x}{2\,k}}\displaystyle\frac{e^{\frac{i\,x}{2}}(e^{-\frac{i\,x}{2}}-e^{\frac{i\,x}{2}})}{e^{\frac{i\,x}{2\,k}}(e^{-\frac{i\,x}{2\,k}}-e^{\frac{i\,x}{2\,k}})} \right\}\\
&= \Re\left\{ e^{\frac{i\,x}{2}}\,\frac{\sin \left(\frac{x}{2}\right)}{\sin \left(\frac{x}{2\,k}\right)}\right\}
  =\cos \left(\frac{x}{2}\right)\frac{\sin \left(\frac{x}{2}\right)}{\sin \left(\frac{x}{2\,k}\right)}=\frac{\sin x}{2\,\sin \left(\frac{x}{2\,k}\right)},
\end{align}
where $\Re\left\{z\right\}$ is the real part of $z$ and $i=\sqrt{-1}$. Now it suffices to show
\begin{equation
\frac{d^2}{dx^2} \left(x^3\frac{d}{dx} f_{p}^{s}(x)\right)=-x\frac{d^3}{dx^3} \left( \frac{\sin x}{\sin\left(\frac{x}{2\,k}\right)}\right).
\end{equation}
Using \eqref{eqn: f''' sin} in $(ii)$, this can be achieved by straightforward calculations. Thus $(iii)$ is true. The proof of $(iv)$ is similar, and we omit the details. We complete the proof of Lemma~\ref{lem: lots of f'''}.\qed

\end{proof}

We provide here an alternative proof of the two double inequalities in Theorem~\ref{thm: Daurizio-Sandor ineq original}.

\begin{proof}[Proof of Theorem~\ref{thm: Daurizio-Sandor ineq original}]
 To this end, we show that for $x\in(0,\pi/2)$, $f_{2}^{c}(x)=\frac{1-\frac{\cos x}{\cos \frac{x}{2}}}{x^2}$ is strictly increasing while $f_{2}^{s}(x)=\frac{2-\frac{\sin x}{\sin \frac{x}{2}}}{x^2}$ is strictly decreasing. These lead to the desired inequalities since it is easy to see that
\begin{alignat}{4
\lim_{x\to0}f_{2}^{c}(x) & =\frac{3}{8},
& \quad
\lim_{x\to\pi/2}f_{2}^{c}(x) & =\frac{4}{\pi^2},
& \qua
\\
\lim_{x\to0}f_{2}^{s}(x) & =\frac{1}{4},
& \quad
\lim_{x\to\pi/2}f_{2}^{s}(x) & =\frac{4}{\pi^2}(2-\sqrt{2}).
& \qua
\end{alignat}
To see $f_{2}^{c}(x)$ is strictly increasing, we employ \eqref{eqn: f''' cos} in Lemma~\ref{lem: lots of f'''} to obtain
\begin{align}\label{eqn: f''' lemma1}\nonumber
\frac{d^2}{dx^2}\left(x^3\frac{d}{dx} f_{2}^{c}(x)\right)
&=-\frac{x}{64} \sec^4\left(\frac{x}{2}\right) \left(-44 \sin\left(\frac{x}{2}\right)+5 \sin \left(\frac{3 x}{2}\right)+\sin \left(\frac{5 x}{2}\right)\right)     \\
&=-\frac{x}{16} \sec^4\left(\frac{x}{2}\right)
\sin \left(\frac{x}{2}\right) (\cos x-2)(\cos x+5)>0.
\end{align}
As $\lim_{x\to0} \frac{d}{dx}\left(x^3\frac{d}{dx} f_{2}^{c}(x)\right)=0$, it follows that $\frac{d}{dx}\left(x^3\frac{d}{dx} f_{2}^{c}(x)\right)>0$. We are led to $x^3\frac{d}{dx} f_{2}^{c}(x)>0$ or $\frac{d}{dx} f_{2}^{c}(x)>0$ since $\lim_{x\to0} \left(x^3\frac{d}{dx} f_{2}^{c}(x)\right)=0$. This that shows $f_{2}^{c}(x)$ is strictly increasing.

By using \eqref{eqn: 2k sin} in Lemma~\ref{lem: lots of f'''}, we have
\begin{equation}
\frac{d^2}{dx^2}\left(x^3\frac{d}{dx} f_{2}^{s}(x)\right)=-\frac{x}{4}\,\sin \left(\frac{x}{2}\right)<0,
\end{equation}
from which we infer that $\frac{d}{dx}\left(x^3\frac{d}{dx} f_{2}^{s}(x)\right)<0$ since $\lim_{x\to0} \frac{d}{dx}\left(x^3\frac{d}{dx} f_{2}^{s}(x)\right)=0$ by $(v)$ of Lemma~\ref{lem: lots of f'''}. Then
\begin{equation}
\displaystyle \frac{d}{dx}\left(x^3\frac{d}{dx} f_{2}^{s}(x)\right)<0
\end{equation}
together with the fact $\lim_{x\to0} \left(x^3\frac{d}{dx} f_{2}^{s}(x)\right)=0$ from $(v)$ of Lemma~\ref{lem: lots of f'''} yields $x^3\frac{d}{dx} f_{2}^{s}(x)<0$ or $\frac{d}{dx} f_{2}^{s}(x)<0$. Thus we have shown that $f_{2}^{s}(x)$ is strictly decreasing. This completes the proof of the theorem.\qed
\end{proof}

We are now in the position to give the proof of Theorem~\ref{thm: Daurizio-Sandor ineq generalized}.

\begin{proof}[Proof of Theorem~\ref{thm: Daurizio-Sandor ineq generalized}]
The proof of the case when $p=2$ has been given in Theorem~\ref{thm: Daurizio-Sandor ineq generalized}. For $p\ge3$, we prove the desired inequalities by showing that $\frac{d}{dx} f_{p}^{\triangle}(x)<0$ for $\triangle=c, s$. Due to $(i)$ of Lemma~\ref{lem: lots of f'''}, we see that $\frac{d^2}{dx^2}\left(x^3\frac{d}{dx} f_{p}^{c}(x)\right)<0$ for $p\ge3$. Instead of employing $(ii)$ of Lemma~\ref{lem: lots of f'''}, we use \eqref{eqn: 2k sin} and \eqref{eqn: 2k+1 sin} in Lemma~\ref{lem: lots of f'''} to conclude that $\frac{d^2}{dx^2}\left(x^3\frac{d}{dx} f_{p}^{s}(x)\right)<0$. Thus we have for $\triangle=c, s$,
\begin{equation}
\frac{d^2}{dx^2}\left(x^3\frac{d}{dx} f_{p}^{\triangle}(x)\right)<0.
\end{equation}
Because of the first vanishing limit in $(v)$ of Lemma~\ref{lem: lots of f'''}, it follows that
\begin{equation}
\frac{d}{dx}\left(x^3\frac{d}{dx} f_{p}^{\triangle}(x)\right)<0,
\end{equation}
which, together with the fact that the second limit in $(v)$ of Lemma~\ref{lem: lots of f'''} vanishes, implies that $x^3\frac{d}{dx} f_{p}^{\triangle}(x)<0$ or $\frac{d}{dx} f_{p}^{\triangle}(x)<0$ for $\triangle=c, s$. It remains to find the following limits:
\begin{alignat}{4
\lim_{x\to0}f_{p}^{c}(x) & =\frac{p^2-1}{2\,p^2},
& \quad
\lim_{x\to\pi/2}f_{p}^{c}(x) & =\frac{4}{\pi^2},
& \qua
\\
\lim_{x\to0}f_{p}^{s}(x) & =\frac{p^2-1}{6\,p},
& \quad
\lim_{x\to\pi/2}f_{p}^{s}(x) & =\frac{4}{\pi^2}\left(p-\csc \left(\frac{\pi}{2p}\right)\right).
& \qua
\end{alignat}
We immediately have
\begin{equation
\frac{4}{\pi^2}=\lim_{x\to\pi/2}f_{p}^{c}(x)<\frac{1-\displaystyle\frac{\cos x}{\cos \frac{x}{p}}}{x^2}<\lim_{x\to0}f_{p}^{c}(x)=\frac{p^2-1}{2\,p^2}
\end{equation}
and
\begin{equation
\frac{4}{\pi^2}\left(p-\csc \left(\frac{\pi}{2p}\right)\right)=\lim_{x\to\pi/2}f_{p}^{s}(x)<\frac{\displaystyle p-\frac{\sin x}{\sin \frac{x}{p}}}{x^2}<\lim_{x\to0}f_{p}^{s}(x)=\frac{p^2-1}{6\,p}.
\end{equation}
The proof is completed.\qed

\end{proof}

\section{Generalized D'Aurizio-S\'andor inequalities for hyperbolic functions}\label{sec: D'Aurizio-Sandor inequalities for hyperbolic functions}

In this section, we show an analogue of Theorem~\ref{thm: Daurizio-Sandor ineq generalized} for the case of hyperbolic functions holds true. Let
\begin{align}
h_{p}^{c}(x)& = \frac{1-\displaystyle\frac{\cosh x}{\cosh \frac{x}{p}}}{x^2}, \\
h_{p}^{s}(x)& = \frac{p-\displaystyle\frac{\sinh x}{\sinh \frac{x}{p}}}{x^2}.
\end{align}
Following the same arguments for proving Lemma~\ref{lem: lots of f'''}, it can be shown that Lemma~\ref{lem: lots of f'''} with $\cos x$, $\sin x$ and $f_{p}^{\triangle}(x)$ $(\triangle=c, s)$ replaced by $\cosh x$, $\sinh x$ and $h_{p}^{\triangle}(x)$ $(\triangle=c, s)$ respectively, remains true. It follows that we can prove $\frac{d}{dx} h_{p}^{\triangle}(x)<0$ for $\triangle=c, s$ as in the proof of Theorem~\ref{thm: Daurizio-Sandor ineq generalized}. It remains to calculate the following limits:
\begin{alignat}{4
\lim_{x\to0}f_{p}^{c}(x) & =\frac{1-p^2}{2\,p^2},
& \quad
\lim_{x\to\pi/2}f_{p}^{c}(x) & =\frac{4}{\pi^2}\left(1-\cosh \left(\frac{\pi }{2}\right)\sech\left(\frac{\pi }{2 p}\right)\right),
& \qua
\\
\lim_{x\to0}f_{p}^{s}(x) & =\frac{1-p^2}{6\,p},
& \quad
\lim_{x\to\pi/2}f_{p}^{s}(x) & =\frac{4}{\pi ^2}\left(p-\sinh \left(\frac{\pi }{2}\right)\csch\left(\frac{\pi }{2p}\right)\right).
& \qua
\end{alignat}
Thus, we have the following analogue of Theorem~\ref{thm: Daurizio-Sandor ineq generalized} for $\cosh x$ and $\sinh x$.
\begin{theorem}\label{thm: Daurizio-Sandor ineq generalized hyperbolic fcn}
Let $0<x<\pi/2$. Then the two double inequalities
\begin{equation}\label{eqn: cosh ineq p>2}
\frac{4}{\pi^2}\left(1-\cosh \left(\frac{\pi }{2}\right)\sech\left(\frac{\pi }{2 p}\right)\right)<\frac{1-\displaystyle\frac{\cosh x}{\cosh \frac{x}{p}}}{x^2}<\frac{1-p^2}{2\,p^2}
\end{equation}
and
\begin{equation}\label{eqn: sinh ineq p>2}
\frac{4}{\pi ^2}\left(p-\sinh \left(\frac{\pi }{2}\right)\csch\left(\frac{\pi }{2p}\right)\right)<\frac{\displaystyle p-\frac{\sinh x}{\sinh \frac{x}{p}}}{x^2}<\frac{1-p^2}{6\,p}
\end{equation}
hold for $p=3, 4 , 5,\cdots$. In particular, the double inequality \eqref{eqn: cosh ineq p>2} is reversed when $p=2$ while the double inequality \eqref{eqn: sinh ineq p>2} remains true when $p=2$.
\end{theorem}

\section{Application of the generalized D'Aurizio-S\'andor inequalities to the Chebyshev polynomials of the second kinds}\label{sec: application}

The first few Chebyshev polynomials of the second kind $U_n(x)$ $(n=0,1,2,\cdots)$ are (\cite{AS2,Chebyshev})
\begin{align}
\label{eqn: Chebyshev poly n=0}
U_0(x)&=1, \\
\label{eqn: Chebyshev poly n=1}
U_1(x)&=2\,x, \\
\label{eqn: Chebyshev poly n=1}
U_2(x)&=4\,x^2-1, \\
\label{eqn: Chebyshev poly n=1}
U_3(x)&=8\,x^3-4\,x, \\
\label{eqn: Chebyshev poly n=1}
U_4(x)&=16\,x^4-12\,x^2+1, \\
\label{eqn: Chebyshev poly n=1}
U_5(x)&=32\,x^5-32\,x^3+6\,x, \\
\label{eqn: Chebyshev poly n=1}
U_6(x)&=64\,x^6-80\,x^4+24\,x^2-1.
\end{align}
In this section, we apply Theorem~\ref{thm: Daurizio-Sandor ineq generalized} to $U_n(x)$ with $x=\cos \theta$. By means of the formula $U_{n}(\cos \theta)=\frac{\sin ((n+1)\,\theta)}{\sin \theta}$, we obtain the following corollary.

\begin{corollary}\label{cor: Daurizio-Sandor ineq applied to Chebyshev poly}
Let $y\in(0,\frac{\pi}{2\,p})$. The double inequality
\begin{equation
\frac{p}{6}\left((1-p^2 )y^2+6\right)<U_{p-1}(\cos y)<p-\frac{4}{\pi^2}\left(p-\csc \left(\frac{\pi}{2p}\right)\right)\,p\,y^2
\end{equation}
holds for $p=2, 3, 4 , 5,\cdots$.
\end{corollary}

\begin{proof}
The double inequality \eqref{eqn: sin ineq p>2} in Theorem~\ref{thm: Daurizio-Sandor ineq generalized} can be written as
\begin{equation
p-\frac{p^2-1}{6\,p}\,x^2<\frac{\sin x}{\sin \frac{x}{p}}<p-\frac{4}{\pi^2}\left(p-\csc \left(\frac{\pi}{2p}\right)\right)\,x^2, \quad x\in(0,\pi/2).
\end{equation}
Letting $x/p=y$, we have
\begin{equation
\frac{p}{6}\left((1-p^2 )y^2+6\right)<\frac{\sin (p\,y)}{\sin y}<p-\frac{4}{\pi^2}\left(p-\csc \left(\frac{\pi}{2p}\right)\right)\,p\,y^2, \quad y\in(0,\frac{\pi}{2\,p}).
\end{equation}
Since $\frac{\sin (p\,y)}{\sin y}=U_{p-1}(\cos y)$, the proof is completed.\qed
\end{proof}

\begin{example}
Letting $p=7$ in Corollary~\ref{cor: Daurizio-Sandor ineq applied to Chebyshev poly} results in the following inequality
\begin{equation
7-56\,y^2<64\,\cos ^6 y-80\,\cos ^4 y+24\,\cos ^2 y-1<7-\frac{196 \left(7-\csc \left(\frac{\pi}{14}\right)\right)}{\pi ^2}\,y^2,
\end{equation}
where $y\in(0,\frac{\pi}{14})\approx (0,0.2244)$ and $\frac{196 \left(7-\csc \left(\frac{\pi}{14}\right)\right)}{\pi ^2}\approx 49.7673$.
\end{example}

\textit{Acknowledgements}. The authors wish to express sincere gratitude to Tom Mollee for his careful reading of the manuscript and valuable suggestions to improve the readability of the paper. Thanks are also due to Chiun-Chuan Chen and Mach Nguyet Minh for the fruitful discussions. The authors are grateful to the anonymous referee for many helpful comments and valuable suggestions on this paper.

\Refs

\bibitem{AS2}
        \by M. Abramowitz and I. A. Stegun (Eds)
        \book Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables
        \publ National Bureau of Standards, Applied Mathematics Series {\bf 55}, 9th printing
        \publaddr Washington
        \year 1970
\endref

\bibitem{D'Aurizio}
        \by J. D'Aurizio
        \paper Refinements of the {S}hafer-{F}ink inequality of arbitrary
              uniform precision
        \jour Math. Inequal. Appl.
        \vol 17
        \issue 4
        \yr 2014
        \pages 1487--1498
\endref

\bibitem{Chebyshev}
        \by T. J. Rivlin
        \book Chebyshev Polynomials
        \publ Wiley
        \publaddr New York
        \year 1990
\endref

\bibitem{Sandor}
        \by J. S\'andor
        \paper On D'Aurizio's trigonometric inequality
        \jour J. Math. Inequal.
        \vol 10
        \issue 3
        \yr 2016
        \pages 885--888
\endref

\endRefs

\end{document}