Approximating monomials using Chebyshev polynomials

Abstract: This paper considers the approximation of a monomial $x^n$ over the interval $[-1,1]$ by a lower-degree polynomial. This polynomial approximation can be easily computed analytically and is obtained by truncating the analytical Chebyshev series expansion of $x^n$. The error in the polynomial approximation in the supremum norm has an exact expression with an interesting probabilistic interpretation. We use this interpretation along with concentration inequalities to develop a useful upper bound for the error.

\section{Motivation and Introduction}\label{sec:intro}
We are interested in approximating the monomial $x^n$ by a polynomial of degree $0 \leq k < n$ over the interval $[-1,1]$. The monomials $1,x,x^2,\dots$ form a basis for $C[-1,1]$, so it seems unlikely that we can represent a monomial in terms of lower degree polynomials. In Figure~\ref{fig:monomialbasis}, we plot a few functions from the monomial basis over $[0,1]$; the basis function look increasingly alike as we take higher and higher powers, i.e., they appear to ``lose independence.'' Numerical analysts often avoid the monomial basis in polynomial interpolation since they result in ill-conditioned Vandermonde matrices, leading to poor numerical performance in finite precision arithmetic. This loss of independence means that it is reasonable to approximate the monomial $x^n$ as a linear combination of lower order monomials, i.e., a lower order polynomial approximation. The natural question to ask, therefore, is: how small can $k$ be so that a well-chosen polynomial of degree $k$ can accurately approximate $x^n$?

\begin{figure}[!ht]
    \centering
    \includegraphics[scale=0.4]{monomial_basis}
    \caption{Visualization of a few monomials in the interval $[0,1]$.}
    \label{fig:monomialbasis}
\end{figure}

The surprising answer to this question is that we can approximate the monomial $x^n$ over $[-1,1]$ by a polynomial of small degree, which we will make precise. Let $\|f\|_\infty = \max_{x \in [-1,1]} | f(x)|$ denote the supremum norm on $C[-1,1]$  and let  $\pi_k^*(\cdot)$ be the best polynomial approximation to $x^n$ in this norm; that is
\[ E_{n,k} := \min_{\pi \in \mathcal{P}_k} \|x^n- \pi(x) \|_\infty= \|x^n - \pi_k^*(x)\|_\infty,   \]
where $\mathcal{P}_k$ is a vector space of polynomials with real coefficients of degree at most $k$.  The minimizer $\pi_k^*(\cdot)$ exists and is unique~\cite[Chapter 10]{trefethen2013approximation}, but does not have a closed form expression. Newman and Rivlin~\cite[Theorem 2]{newman1976approximation} showed that\footnote{We briefly mention that the notation our manuscript differs from~\cite{newman1976approximation} in that we reverse the roles of $n$ and $k$.}
\begin{equation}\label{eqn:newmanrivlin} \frac{p_{n,k}}{4e} \leq \|x^n -  \pi_k^*(x)\|_\infty \leq p_{n,k},  \end{equation}
where the term $p_{n,k}$ is given by the formula  \[ p_{n,k}  = \frac{1}{2^{n-1}} \sum_{j = \lfloor (n+k)/2 \rfloor +1}^n \binom{n}{j}. \]
Since $p_{n,k}$ involves the sum of binomial coefficients, it has a probabilistic interpretation which we explore in Section~\ref{sec:prob}.

To see why a small $k$ is sufficient, consider the upper bound $p_{n,k}$. In Section~\ref{sec:prob} we use the probabilistic interpretation to obtain the following bound  $p_{n,k} \leq 2\exp\left(-{k^2}/{2n}\right)$. Suppose we are given a user-defined tolerance $\epsilon > 0$. To ensure
\[ \|x^n -  \pi_k^*(x)\|_\infty \leq \epsilon,\]
we need to  choose $k \geq \sqrt{2n\log(2/\epsilon)}$. The accuracy of the polynomial approximation is visualized in Figure~\ref{fig:approx}, where in the left panel we plot the monomial $x^{n}$ for $n=75$ and the best polynomial approximation $\pi_{k}^*$ for $k=5,15,25$. The polynomial $\pi_k^*$ is computed using the Remez algorithm, implemented in chebfun~\cite{Driscoll2014}. We see that for $k=25$, the polynomial approximation looks very accurate. In the right panel, we display $p_{n,k}$, which is the upper bound of the best polynomial approximation, as well as the upper bound for $p_{n,k}$. We see that $p_{n,k}$ and its upper bound both have sharp decay with increasing $k$.  Numerical evidence in~\cite{nakatsukasa2018rational} further confirms this analysis; the authors show that the error $E_{n,k}$ behaves approximately like $\frac12 \text{erfc}(k/\sqrt{n})$, where erfc is the complementary error function. .

\begin{figure}[!ht]
    \centering
    \includegraphics[scale=0.5]{monomial_approx.eps}
    \caption{(left) Approximation of the monomial $x^{n}$ for $n=75$ by $\pi_k^*$, (right) $p_{n,k}$ and its upper bound $2\exp(-k^2/2n)$. The visualization is restricted to the interval $[0,1]$.}
    \label{fig:approx}
\end{figure}

 Polynomial and rational approximations to the monomial has received considerable attention in the approximation theory community, and surveys of various results can be found in~\cite{reddy1987approximations,nakatsukasa2018rational}. Polynomial approximations to high order monomials have many applications in numerical analysis.   This key insight was exploited by Cornelius Lanczos~\cite{lanczos1988applied} in his ``$\tau$-method'' for the numerical solution of differential equations. For a simulating discussion on this topic, please see~\cite{nakatsukasa2018rational}. In numerical linear algebra, this has been exploited to efficiently compute matrix powers and the Schatten p-norm of a matrix~\cite{avron2011randomized,dudley2020monte}.

 In this short note, we show to construct a polynomial approximation $x^n \approx \phi_k(x)$ using a truncated Chebyshev polynomial expansion. The error in the truncated representation equals the sum of the discarded coefficients and is precisely $p_{n,k}$. The polynomial $\phi_k$ and the resulting error can both be computed analytically and, therefore, is of great practical use.  We briefly review Chebyshev polynomials in Section~\ref{sec:cheby} and state and prove the main result in Section~\ref{sec:main}. In Section~\ref{sec:prob}, we explore probabilistic interpretations of $p_{n,k}$ and obtain bounds for partial sums of binomial coefficients.

\section{Chebyshev polynomials}\label{sec:cheby}
The Chebyshev polynomials of the first kind $T_n(x)$ for $n=0,1,2,\dots$ can be represented as
\[T_n(x) = \cos(n\arccos{x}) \qquad x \in [-1,1]. \]
Starting with $T_0(x) = 1$, $T_1(x) = x$, the Chebyshev polynomials satisfy a recurrence relationship of the form  $T_{n+1}(x) = 2xT_n(x) - T_{n-1}(x)$ for $n\geq 1$. The Chebyshev polynomials are orthogonal with respect to the weighted inner product \[\langle u,v \rangle = \int_{-1}^1 w(x) u(x) v(x) dx \]
where the weight function takes the form $w(x) = (1-x^2)^{-1/2}$. Any function $f \in C[-1,1]$ that is Lipschitz continuous can be represented in terms of a Chebyshev polynomial expansion of the form
\[ f(x) = \frac12 c_0 + \sum_{j=1}^\infty c_jT_j(x) \qquad x \in [-1,1], \]
where the coefficients $c_j$ are obtained as  $c_j = \frac{2}{{\pi}}\langle f(x), T_j(x)\rangle$ and the series is uniformly convergent. The monomial $x^n$ is rather special since it has the following exact representation in terms of the Chebyshev polynomials~\cite[Section 4]{cody1970survey} \begin{equation}\label{eqn:xncheby} x^n = \sum_{j=0}^{n}{}^{'} c_jT_j(x),\end{equation}
where ${}^{'}$ means the summand corresponding to $j=0$ is halved (if it appears) and  the coefficients $c_j$ for $j=0,\dots,n$ are
\begin{equation}\label{eqn:cj} c_j = \left\{ \begin{array}{ll} 2^{1-n}\binom{n}{(n-j)/2}  & n-j \text{ even} \\ 0 & \text{otherwise}.\end{array}  \right.\end{equation}
Equation~\eqref{eqn:xncheby} takes a more familiar form, when we consider the trigonometric perspective of Chebyshev polynomials. For example, the well-known trigonometric identity $\cos(3\theta) = 4\cos^3 \theta - 3\cos \theta$, can be arranged as \[ \cos^3 \theta = \frac{3}{4}\cos\theta + \frac{1}{4} \cos(3\theta) = \frac{1}{2^2} \left( \binom{3}{1} \cos\theta + \binom{3}{0} \cos(3\theta)\right). \]
With $x=\cos\theta$, we get $x^3 = 2^{-2} \binom{3}{1} T_1(x) + 2^{-2} \binom{3}{0} T_3(x)$. For completeness, we provide a derivation of~\eqref{eqn:xncheby} in Appendix~\ref{app:der}. It is important to note here that the series in~\eqref{eqn:xncheby} is finite, but can be truncated to obtain an accurate approximation; see Section~\ref{sec:main}.

Chebyshev polynomials have many applications in approximation theory and numerical analysis~\cite{trefethen2013approximation} but we limit ourselves to two such examples here. First, if the function is differentiable $r$ times or analytic, the Chebyshev coefficients exhibit decay (algebraic or geometric respectively). Therefore, the Chebyshev series can be truncated to obtain an  polynomial approximation of the function and the accuracy of the approximation depends on the rate of decay of the coefficients. Another application of Chebyshev polynomials is in the theory and practice of polynomial interpolation. The polynomial ${q}_{n-1}^*(x) = x^{n} - 2^{1-n}T_{n}(x)$ solves the minimax problem
\begin{equation}\label{eqn:minmax} \min_{q \in \mathcal{P}_{n-1}} \|x^{n} - q(x)\|_\infty = 2^{1-n}. \end{equation}
Based on the minimax characterization, to interpolate a function over $[-1,1]$ by a polynomial of degree $n-1$, the function to be interpolated  should be evaluated at the roots of the Chebyshev polynomial $T_{n}$ given by the points $x_j = \cos\left(\frac{2j+1}{2n }\pi\right)$ for $j=0,\dots,n-1$.

\section{Main result}\label{sec:main}
We construct the polynomial approximation $ x^n \approx \phi_k(x)$ by truncating the Chebyshev polynomial expansion in~\eqref{eqn:xncheby} beyond the term $j=k$. That is \[\phi_k(x) := \sum_{j=0}^k{}^{'}c_jT_j(x).\]
Our main result is the following theorem, which quantifies the error in the polynomial approximation. The proof of this theorem is based on the expression in~\eqref{eqn:xncheby}. We believe this result is new.
\begin{theorem}\label{thm:main}
The error in the polynomial approximation $\phi_k(x)$ satisfies
\[  \|x^n -  \phi_k(x)\|_\infty = p_{n,k}. \]
\end{theorem}

\begin{proof}
 From~\eqref{eqn:xncheby}, $x^{n} -  \phi_k(x) = \sum_{j=k+1}^nc_jT_j(x)$. Using triangle inequality, we find that $\|x^n-\phi_k(x)\|_\infty \leq \sum_{j=k+1}^n c_j$ since the coefficients are nonnegative and the Chebyshev polynomials are bounded as $|T_j(x)| \leq 1$. Substituting the coefficients $c_j$ from~\eqref{eqn:cj}, to get
\begin{equation}
    \label{eqn:inter}
    \|x^n-\phi_k(x)\|_\infty \leq \frac{1}{2^{n-1}}\sum_{\substack{j=k+1 \\ n-j \text{ even}}}^n \binom{n }{(n-j)/2}.
\end{equation}
Using the properties of the binomial coefficients, the summation simplifies as
\[
\sum_{\substack{j=k+1 \\ n-j \text{ even}}}^n \binom{n }{(n-j)/2} =  \sum_{\substack{j=k+1 \\ n+j \text{ even}}}^n \binom{n }{(n+j)/2} =\sum_{j=\lfloor (n+k)/2\rfloor+1}^n \binom{n }{j}. \]
Plug this identity into~\eqref{eqn:inter} to get $\|x^n -  \phi_k(x)\|_\infty \leq p_{n,k}$. The bound is clearly achieved at $x= 1$, where all the Chebyshev polynomials take the value $1$.
\end{proof}
This theorem shows that the polynomial approximation $\phi_k$ is nearly optimal, and the error due to this approximation is $p_{n,k}$. However, it is the optimal polynomial for the special case $k=n-1$. It is easy to see that $x^n - \phi_{n-1}(x) = 2^{1-n}T_n(x)$ and so $\phi_{n-1}$ is the same as the best polynomial approximation $q^*_{n-1}$ in~\eqref{eqn:minmax}. For $k < n-1$, from~\eqref{eqn:newmanrivlin} and Theorem~\ref{thm:main}
\[ \frac{ \|x^n-\phi_k(x)\|_\infty}{4e} \leq  \|x^n-\pi_k^*(x)\|_\infty \leq  \|x^n-\phi_k(x)\|_\infty, \]
so that the error in the Chebyshev polynomial approximation is suboptimal by at most the factor $4e \approx 10.87$. Therefore, by using $\phi_k$ we lose only one significant digit of accuracy compared to $\pi^*_k$.

\section{A probabilistic digression}\label{sec:prob}
In Section~\ref{sec:intro}, we saw that the error in the monomial approximation depends on  $p_{n,k}$. Since $p_{n,k}$ depends on the sum of binomial coefficients, it has a probabilistic interpretation. Newman and Rivlin~\cite{newman1976approximation} observed that if a fair coin is tossed $n$ times, $p_{n,k}$ is the probability that the magnitude of the difference between the number of heads and the number of tails exceeds $k$. They used this insight along with the de Moivre-Laplace theorem~\cite[Section 1.3]{vershynin2018high} (which is a special case of the Central Limit Theorem) to obtain the approximation $p_{n,k} \approx 2 \text{erfc}(k/\sqrt{n})$.

To convert this into a rigorous inequality for $p_{n,k}$ we use a different tool from probability, namely, concentration inequalities. The inequalities are useful in quantifying how much a random variable deviates from its mean.  We start with the following alternative interpretation for $p_{n,k}$: it is twice the probability that greater than   $\lfloor (n+k)/2 \rfloor$ coin tosses result in heads (or equivalently tails). We associate each coin toss with an independent Bernoulli random variable $X_i$ with parameter $p=1/2$ since the coin is fair. The random variable $X = \sum_{i=1}^nX_i$ has the Binomial distribution with parameters $n$ and $p$. Then,
\[ p_{n,k} =  2\mathbb{P}\left( \lfloor (n+k)/2 \rfloor +1 \leq X \leq n\right)  \leq 2\mathbb{P}\left(X  \geq (n+k)/2 \right). \]
Since $X$ has the Binomial distribution, we can once again use the de Moivre-Laplace theorem, to say that as $n\rightarrow \infty$,
\[ \frac{X - np}{\sqrt{np(1-p)}} \longrightarrow \mathcal{N}(0,1), \qquad \text{in distribution}.\]
  Roughly speaking, this theorem says that $X$ behaves as a normal random variable with mean $n/2$ and variance $n/4$. Since the tails of normal distributions decay exponentially, we expect that $X$ lies in the range $\frac{n}{2} \pm 1.96 \sqrt{\frac{n}{4}}$ with nearly $95\%$ probability; alternatively, the probability that it is outside this range is very small. To make this more precise, we apply Hoeffding's concentration inequality~\cite[Theorem 2.2.6]{vershynin2018high}, to obtain
\[  \mathbb{P}\left(X  \geq (n+k)/2 \right) = \mathbb{P}\left(X - \mathbb{E}[X] \geq k/2 \right) \leq \exp\left(-\frac{k^2}{2n}\right). \]
This gives our desired bound $p_{n,k} \leq 2\exp(-{k^2}/{2n})$.

We can use a similar technique to prove the following result which may be of independent interest. If $0 \leq k \leq n/2$, then
\[ \sum_{j=0}^k \binom{n}{j} \leq 2^n \exp\left( - \frac{(n-2k)^2}{2n}\right).\]
 Other concentration inequalities such as Chernoff and Bernstein (see~\cite[Chapter 2]{vershynin2018high}) also give equally interesting bounds. We invite the reader to explore such results.

\section{Acknowledgements}
The author would like to thank Alen Alexanderian, Ethan Dudley, Ivy Huang, Ilse Ipsen, and Nathan Reading for comments and feedback. The work was supported by the National Science Foundation through the grants DMS-1745654 and DMS-1845406.
\section{Declaration of Interest}
The author has no relevant financial or non-financial competing interests to report.


The quaternionic Gauss-Lucas Theorem

Abstract: The classic Gauss-Lucas Theorem for complex polynomials of degree $d\ge2$ has a natural reformulation over quaternions, obtained via rotation around the real axis. We prove that such a reformulation is true only for $d=2$. We present a new quaternionic version of the Gauss-Lucas Theorem valid for all $d\geq2$, together with some consequences.

\section{Introduction}

Let $p$ be a complex polynomial of degree $d\geq2$ and let $p'$ be its derivative. The Gauss-Lucas Theorem asserts that the zero set of $p'$ is contained in the convex hull $\mathcal K(p)$ of the zero set of $p$. The classic proof uses the logarithmic derivative of $p$ and it strongly depends on the commutativity of ${\mathbb{C}}$.

This note deals with a quaternionic version of such a classic result.
We refer the reader to \cite{GeStoSt2013} for the notions and properties concerning the algebra ${\mathbb{H}}$ of quaternions we need here.
The ring ${\mathbb{H}}[X]$ of quaternionic polynomials is defined by fixing the position of the coefficients w.r.t.\  the indeterminate $X$ (e.g.\ on the right) and by imposing commutativity of $X$ with the coefficients when two polynomials are multiplied together (see e.g.\ \cite[\S 16]{Lam}). Given two polynomials $P,Q\in{\mathbb{H}}[X]$, let $P\cdot Q$ denote the product obtained in this way. If $P$ has real coefficients, then $(P\cdot Q)(x)=P(x)Q(x)$.
In general, a direct computation (see \cite[\S 16.3]{Lam}) shows that if $P(x)\ne0$, then
\begin{equation}\label{product}
(P\cdot Q)(x)=P(x)Q(P(x)^{-1}xP(x)),
\end{equation}
while $(P\cdot Q)(x)=0$ if $P(x)=0$.
In this setting, a {(left) root} of a polynomial $P(X)=\sum_{h=0}^dX^h a_h$ is an element $x\in{\mathbb{H}}$ such that $P(x)=\textstyle\sum_{h=0}^dx^h a_h=0$.

Given  $P(X)=\sum_{k=0}^dX^ka_k\in {\mathbb{H}}[X]$,
consider the polynomial $P^c(X)=\sum_{k=0}^dX^k\bar a_k$ and the \emph{normal polynomial} $N(P)=P \cdot P^c=P^c\cdot P$. Since $N(P)$ has real coefficients, it can be identified with a polynomial in ${\mathbb{R}}[X]\subset{\mathbb{C}}[X]$. We recall that a subset $A$ of ${\mathbb{H}}$ is called \emph{circular} if, for each $x\in A$, $A$ contains the whole set (a 2-sphere if $x\not\in{\mathbb{R}}$, a point if $x\in{\mathbb{R}}$)
\begin{equation}\label{sx}
{\mathbb{S}}_x=\{pxp^{-1}\in{\mathbb{H}}\;|\;p\in{\mathbb{H}}^*\},
\end{equation}
where ${\mathbb{H}}^*:={\mathbb{H}}\setminus\{0\}$. In particular, for any imaginary unit $I\in{\mathbb{H}}$, ${\mathbb{S}}_I={\mathbb{S}}$ is the 2-sphere of all imaginary units in ${\mathbb{H}}$. For every subset $B$ of ${\mathbb{H}}$ we define its \emph{circularization} as the set $\bigcup_{x\in B}{\mathbb{S}}_x$. It is well-known (\cite[\S3.3]{GeStoSt2013}) that the zero set $V(N(P))\subset{\mathbb{H}}$ of the normal polynomial is the circularization of the zero set $V(P)$, which consists of isolated points or isolated 2-spheres of the form \eqref{sx} if $P\neq0$.

Let the degree $d$ of $P$ be at least 2 and let $P'(X)=\sum_{k=1}^dX^{k-1}ka_k$ be the derivative of $P$.
It is known (see e.g.\ \cite{GGSarxiv}) that the Gauss-Lucas Theorem does not hold directly for quaternionic polynomials. For example, the polynomial $P(X)=(X-i)\cdot (X-j)=X^2-X(i+j)+k$ has zero set $V(P)=\{i\}$, while $P'$ vanishes at $x=(i+j)/2$.

Since the zero set $V(P
$ of $P$ is contained in the set $V(N(P)
$, a natural reformulation in ${\mathbb{H}}[X]$ of the classic Gauss-Lucas Theorem is the following: $V(N(P'))\subset\mathcal K(N(P))$ or equivalently
\begin{equation}\label{eq:g-l}
V(P')\subset\mathcal K(N(P)),
\end{equation}
where $\mathcal K(N(P))$ denotes the convex hull of $V(N(P))$ in ${\mathbb{H}}$. This set is equal to the circularization of the convex hull of the zero set of $N(P)$ viewed as a polynomial in ${\mathbb{C}}[X]\subset{\mathbb{H}}[X]$. Recently two proofs of the above inclusion \eqref{eq:g-l} were presented in \cite{VlacciGL,GGSarxiv}.

Our next two propositions prove that inclusion \eqref{eq:g-l} is correct in its full generality only when $d=2$.

\section{Gauss-Lucas polynomials}

\begin{definition}
Given a polynomial $P\in{\mathbb{H}}[X]$ of degree $d\geq 2$ we say that $P$ is a \emph{Gauss-Lucas polynomial} if $P$ satisfies \eqref{eq:g-l}.
\end{definition}

\begin{proposition}
If $P$ is a polynomial in ${\mathbb{H}}[X]$ of degree $2$, then $V(N(P))={\mathbb{S}}_{x_1}\cup{\mathbb{S}}_{x_2}$ for some $x_1,x_2\in{\mathbb{H}}$ (possibly with ${\mathbb{S}}_{x_1}={\mathbb{S}}_{x_2}$) and
\[
V(P')\subset\bigcup_{y_1\in{\mathbb{S}}_{x_1},y_2\in{\mathbb{S}}_{x_2}}\left\{\frac{y_1+y_2}{2}\right\}.
\]
In particular every polynomial $P\in{\mathbb{H}}[X]$ of degree $2$ is a Gauss-Lucas polynomial.
\end{proposition}
\begin{proof}
Let $P(X)=X^2a_2+Xa_1+a_0\in{\mathbb{H}}[X]$ with $a_2\neq0$. Since $P\cdot a_2^{-1}=Pa_2^{-1}$, $(P\cdot a_2^{-1})'=P'\cdot a_2^{-1}=P'a_2^{-1}$, we can assume $a_2=1$. Consequently, $P(X)=(X-x_1)\cdot(X-x_2)=X^2-X(x_1+x_2)+x_1x_2$ for some $x_1,x_2\in{\mathbb{H}}$. Then $x_1\in V(P)$ and $\bar x_2\in V(P^c)$, since $P^c(X)=(X-\bar x_2)\cdot(X-\bar x_1)$. Therefore $x_1,x_2\in V(N(P))$. On the other hand $V(P')=\{(x_1+x_2)/2\}$ as desired.
\end{proof}

\begin{remark}
Let $P(X)=\sum_{k=0}^dX^ka_k\in{\mathbb{H}}[X]$ of degree $d\geq2$ and let $Q:=P\cdot a_d^{-1}$ be the corresponding monic polynomial. Since $V(P)=V(Q)$ and $V(P')=V(Q')$, $P$ is a Gauss-Lucas polynomial if and only if $Q$ is.
\end{remark}

\begin{proposition} \label{prop:-4}
Let $P\in{\mathbb{H}}[X]$ of degree $d\geq3$. Suppose that $N(P)(X)=X^{2e}\cdot(X^2+1)^{d-e}$ for some $e<d$ and that $N(P')(X)=\sum_{k=0}^{2d-2}X^kb_k$ contains a unique monomial of odd degree, that is, $b_k\neq0$ for a unique odd $k$. Then $P$ is not a Gauss-Lucas polynomial.
\end{proposition}
\begin{proof}
Since $V(N(P))\subset\{0\}\cup{\mathbb{S}}$, $\mathcal K(N(P))\subset\im({\mathbb{H}})=\{x\in{\mathbb{H}}\,|\,\re(x)=0\}$. Then it suffices to show that $N(P')$ has at least one root in ${\mathbb{H}}\setminus\im({\mathbb{H}})$. Let $G,L\in{\mathbb{R}}[X]$ be the unique real polynomials such that $G(t)+iL(t)=N(P')(it)$ for every $t\in{\mathbb{R}}$. Since $N(P')$ contains a unique monomial of odd degree, say $X^{2\ell+1}b_{2\ell+1}$, then $L(t)=(-1)^\ell b_{2\ell+1}t^{2\ell+1}$ and hence $V(N(P'))\cap i{\mathbb{R}}\subset\{0\}$. Being $V(N(P'))$ a circular set, it holds $V(N(P'))\cap\im({\mathbb{H}})\subset\{0\}$.
Since $N(P')$ contains at least two monomials, namely those of degrees $2\ell+1$ and $2d-2$, we infer that it must have a nonzero root. Therefore
$V(N(P'))\not\subset\{0\}$ and hence $V(N(P'))\not\subset\im({\mathbb{H}})$, as desired.

\end{proof}

\begin{corollary}\label{counterexample}
Let $d\geq3$ and let
\[
P(X)=X^{d-3}\cdot(X-i)\cdot(X-j)\cdot(X-k).
\]
Then $N(P)(X)=X^{2d-6}\cdot(X^2+1)^3$ and $N(P')$ contains a unique monomial of odd degree, namely $-4X^{2d-5}$. In particular $P$ is not a Gauss-Lucas polynomial.
\end{corollary}
\begin{proof}
By a direct computation we obtain:
\begin{align}
P(X)&=X^d-X^{d-1}(i+j+k)+X^{d-2}(i-j+k)+X^{d-3},\\
P'(X)&=dX^{d-1}-(d-1)X^{d-2}(i+j+k)+(d-2)X^{d-3}(i-j+k)+(d-3)X^{d-4},\label{P1}
\\
N(P')(X)&=d^2X^{2d-2}+3(d-1)^2X^{2d-4}-4X^{2d-5}+3(d-2)^2X^{2d-6}+(d-3)^2X^{2d-8}.
\end{align}
Proposition \ref{prop:-4} implies the thesis.
\end{proof}

Let $I\in{\mathbb{S}}$ and let ${\mathbb{C}}_I\subset{\mathbb{H}}$ be the complex plane generated by 1 and $I$. Given a polynomial $P\in{\mathbb{H}}[X]$, we will denote by $P_I:{\mathbb{C}}_I\to{\mathbb{H}}$ the restriction of $P$ to ${\mathbb{C}}_I$. {\it If $P_I$ is not constant, we will denote by $\mathcal K_{{\mathbb{C}}_I}(P)$ the convex hull in the complex plane ${\mathbb{C}}_I$ of the zero set $V(P_I)=V(P)\cap{\mathbb{C}}_I$. If $P_I$ is constant, we set $\mathcal K_{{\mathbb{C}}_I}(P)={\mathbb{C}}_I$.}

If ${\mathbb{C}}_I$ contains every coefficient of $P\in{\mathbb{H}}[X]$, then we say that $P$ is a \emph{${\mathbb{C}}_I$-polynomial}.

\begin{proposition}
The following holds:
\begin{enumerate}
 \item[(1)] Every ${\mathbb{C}}_I$-polynomial of degree $\geq2$ is a Gauss-Lucas polynomial.
 \item[(2)] Let $d\geq3$, let ${\mathbb{H}}_d[X]=\{P\in{\mathbb{H}}[X]\,|\,\deg(P)=d\}$ and let $E_d[X]$ be the set of all elements of ${\mathbb{H}}_d[X]$ that are not Gauss-Lucas polynomials.

 Identify each $P(X)=\sum_{k=0}^dX^ka_k$ in ${\mathbb{H}}_d[X]$ with $(a_0,\ldots,a_d)\in{\mathbb{H}}^d\times{\mathbb{H}}^*\subset{\mathbb{R}}^{4d+4}$ and endow ${\mathbb{H}}_d[X]$ with the relative Euclidean topology.
 Then $E_d[X]$ is a nonempty open subset of ${\mathbb{H}}_d[X]$. Moreover $E_d[X]$ is not dense in ${\mathbb{H}}_d[X]$, being $X^d-1$ an interior point of its complement.
\end{enumerate}
\end{proposition}
\begin{proof}
If $P$ is a ${\mathbb{C}}_I$-polynomial, then $P_I$ can be identified with an element of ${\mathbb{C}}_I[X]$. Consequently, the classic Gauss-Lucas Theorem gives $V(P')\cap{\mathbb{C}}_I=V(P_I')\subset\mathcal K_{{\mathbb{C}}_I}(P)$.
The zero set of the ${\mathbb{C}}_I$-polynomial $P'$ has a particular structure (see \cite[Lemma~3.2]{GeStoSt2013}): $V(P')$ is the union of $V(P')\cap{\mathbb{C}}_I$ with the set of spheres ${\mathbb{S}}_x$ such that $x,\bar x\in V(P'_I)$. It follows that
\[
V(P')\subset\mathcal K(N(P)).
\]

This proves (1).

Now we prove (2). By Corollary \ref{counterexample} we know that $E_d[X]\neq\emptyset$. If $P\in E_d[X]$, then $V(N(P'))\not\subset\mathcal K(N(P))$.
$N(P)$ and $N(P')$
are polynomials with real coefficients. Since the roots of $N(P)$ and of $N(P')$ depend continuously on the coefficients of $P$ and $\mathcal K(N(P))$ is closed in ${\mathbb{H}}$, for every $Q\in{\mathbb{H}}_d[X]$ sufficiently close to $P$, $V(N(Q'))$ is not contained in $\mathcal K(N(Q))$, that is $Q\in E_d[X]$.

To prove the last statement, observe that $P(X)=X^d-1$ is not in $E_d[X]$ from part (1). Since $V(P')=V(N(P'))=\{0\}$ is contained in the interior of the set $\mathcal K(N(P))$, for every $Q\in{\mathbb{H}}_d[X]$ sufficiently close to $P$, $V(Q')$ is still contained in $\mathcal K(N(Q))$.
\end{proof}

\section{A quaternionic Gauss-Lucas Theorem}

Let $P(X)=\sum_{k=0}^dX^ka_k\in{\mathbb{H}}[X]$ of degree $d\geq2$. For every $I\in{\mathbb{S}}$, let $\pi_I:{\mathbb{H}}\to{\mathbb{H}}$ be the orthogonal projection onto ${\mathbb{C}}_I$ and $\pi_I^\bot=id-\pi_I$. Let $P^I(X):=\sum_{k=1}^dX^ka_{k,I}$ be the ${\mathbb{C}}_I$-polynomial with coefficients $a_{k,I}:=\pi_I(a_k)$.

\begin{definition}
We define the \emph{Gauss-Lucas snail of $P$} as the following subset $\mk{sn}(P)$ of ${\mathbb{H}}$:
\[
\mk{sn}(P):=\bigcup_{I\in{\mathbb{S}}}\mathcal K_{{\mathbb{C}}_I}(P^I).
\]
\end{definition}

Our quaternionic version of the Gauss-Lucas Theorem reads as follows.

\begin{theorem}\label{thm}
For every polynomial $P\in{\mathbb{H}}[X]$ of degree $\geq2$,
\begin{equation}\label{snail}
V(P')\subset\mk{sn}(P).
\end{equation}
\end{theorem}
\begin{proof}
Let $P(X)=\sum_{k=0}^dX^ka_k$ in ${\mathbb{H}}_d[X]$ with $d\geq2$. We can decompose the restriction of $P$ to ${\mathbb{C}}_I$ as $P_I=\pi_I\circ P_I+\pi_I^\bot\circ P_I={P^I}_{|{\mathbb{C}}_I}+\pi_I^\bot\circ P_I$. If $x\in{\mathbb{C}}_I$, then $P^I(x)\in{\mathbb{C}}_I$ while $(\pi_I^\bot\circ P_I)(x)\in{\mathbb{C}}_I^\bot$. The same decomposition holds for $P'$. This implies that $V(P')\cap{\mathbb{C}}_I\subset V((P^I)')\cap{\mathbb{C}}_I$. The classic Gauss-Lucas Theorem applied to $P^I$ on ${\mathbb{C}}_I$ gives $V(P')\cap{\mathbb{C}}_I\subset \mathcal K_{{\mathbb{C}}_I}(P^I)$. Since $V(P')=\bigcup_{I\in{\mathbb{S}}}(V(P')\cap{\mathbb{C}}_I)$, the inclusion \eqref{snail} is proved.
\end{proof}

If $P$ is monic Theorem \ref{thm} has the following equivalent formulation: {\it For every monic polynomial $P\in{\mathbb{H}}[X]$ of degree $\geq2$, it holds}
\begin{equation}
\mk{sn}(P')\subset\mk{sn}(P).
\end{equation}

\begin{remark}\label{rem:monic}
If $P$ is a nonconstant monic polynomial in ${\mathbb{H}}[X]$, then two properties hold:
\begin{itemize}
\item[(a)]
$\mathcal K_{{\mathbb{C}}_I}(P^I)$ is a compact subset of ${\mathbb{C}}_I$ for every $I\in{\mathbb{S}}$.
\item[(b)]
$\mathcal K_{{\mathbb{C}}_I}(P^I)$ depends continuously on $I$.
\end{itemize}
Let $I\in{\mathbb{S}}$. Since $P$ is monic, also $P^I$ is a monic, nonconstant polynomial and then $\mathcal K_{{\mathbb{C}}_I}(P^I)$ is a compact subset of ${\mathbb{C}}_I$. This prove property (a). To see that (b) holds, one can apply the Continuity theorem for monic polynomials (see e.g.~\cite[Theorem~1.3.1]{RahmanSchmeisser}). The roots of $P^I$ depend continuously on the coefficients of $P^I$, which in turn depend continuously on $I$. Therefore the convex hull $\mathcal K_{{\mathbb{C}}_I}(P^I)$ depends continuously on $I$.
Observe that (a) and (b) can not hold for polynomials that are not monic. For example, let $P(X)=X^2i$. Then, given $I=\alpha_1i+\alpha_2j+\alpha_3k\in{\mathbb{S}}$,  $\mathcal K_{{\mathbb{C}}_I}(P^I)=\{0\}$ if $\alpha_1\ne0$ and $\mathcal K_{{\mathbb{C}}_I}(P^I)={\mathbb{C}}_I$ if $\alpha_1=0$, since in this case $P^I$ is constant.

\end{remark}

\begin{remark}\label{rem:closed}
If $P$ is a monic polynomial in ${\mathbb{H}}[X]$ of degree $\geq2$, then its Gauss-Lucas snail is a closed subset of ${\mathbb{H}}$.
To prove this fact, consider $q\in{\mathbb{H}}\setminus\mk{sn}(P)$ and choose $I\in{\mathbb{S}}$ such that $q\in{\mathbb{C}}_I$. Write $q=\alpha+I\beta\in{\mathbb{C}}_I$ for some $\alpha,\beta\in{\mathbb{R}}$ and define $z:=\alpha+i\beta\in{\mathbb{C}}$. Since $P$ is monic,
$\mk{sn}(P)\cap{\mathbb{C}}_I=\mathcal K_{{\mathbb{C}}_I}(P^I)$ is a compact subset of ${\mathbb{C}}_I$. Moreover,
$\mathcal K_{{\mathbb{C}}_I}(P^I)$ depends continuously on $I$, and then there exist an open neighborhood $U_I$ of $z$ in ${\mathbb{C}}$ and an open neighborhood $W_I$ of $I$ in ${\mathbb{S}}$ such that the set
\[
\textstyle
[U_I,W_I]:=\bigcup_{J\in W_I}\{a+Jb\in{\mathbb{C}}_J\;|\; a+ib\in U_I\}
\]
is an open neighborhood of $q$ in $\bigcup_{J\in W_I}{\mathbb{C}}_J$, and it is disjoint from $\mk{sn}(P)$. If $q\not\in{\mathbb{R}}$ then $q$ is an interior point of ${\mathbb{H}}\setminus\mk{sn}(P)$, because $[U_I,W_I]$ is a neighborhood of $q$ in ${\mathbb{H}}$ as well.
Now assume that $q\in{\mathbb{R}}$. Since ${\mathbb{S}}$ is compact there exist $I_1,\ldots,I_n\in{\mathbb{S}}$ such that $\bigcup_{\ell=1}^nW_{I_\ell}={\mathbb{S}}$. It follows that $[\bigcap_{\ell=1}^nU_\ell,{\mathbb{S}}]$ is a neighborhood of $q$ in ${\mathbb{H}}$, which is  disjoint from $\mk{sn}(P)$. Consequently $q$ is an interior point of ${\mathbb{H}}\setminus\mk{sn}(P)$ also in this case. This proves that $\mk{sn}(P)$ is closed in ${\mathbb{H}}$.

In Proposition \ref{pro:sn-compact} below we will show that the Gauss-Lucas snail of a monic polynomial in ${\mathbb{H}}[X]$ of degree $\geq2$ is also a compact subset of ${\mathbb{H}}$.
\end{remark}

If all the coefficients of $P$ are real, then $\mk{sn}(P)$ is a circular set.
In general, $\mk{sn}(P)$ is neither closed nor bounded nor circular, as shown in the next example.

\begin{example}
Let $P(X)=X^2i+X$. Given $I=\alpha_1i+\alpha_2j+\alpha_3k\in{\mathbb{S}}$, $P^I(x)=X^2I\alpha_1+X$ and then $\mathcal K_{{\mathbb{C}}_I}(P^I)=\{0\}$ if $\alpha_1=0$ while $\mathcal K_{{\mathbb{C}}_I}(P^I)$ is the segment from 0 to $I{\alpha_1}^{-1}$ if $\alpha_1\ne0$. It follows that
$\mk{sn}(P)=\{x_1i+x_2j+x_3k\in\im({\mathbb{H}})\,|\,0<x_1\le1\}\cup\{0\}$.
Finally, observe that the monic polynomial $Q(X)=-P(X)\cdot~i=X^2-Xi$ corresponding to $P$ has compact Gauss-Lucas snail $\mk{sn}(Q)=\{x_1i+x_2j+x_3k\in\im({\mathbb{H}})\,|\,(x_1-1/2)^2+x_2^2+x_3^2\le 1/4\}$.
\end{example}

\begin{remark}
Even for ${\mathbb{C}}_I$-polynomials, the Gauss-Lucas snail of $P$ can be strictly smaller than the circular convex hull $\mathcal K(N(P))$. For example, consider the ${\mathbb{C}}_i$-polynomial $P(X)=X^3+3X+2i$, with zero sets $V(P)=\{-i,2i\}$ and $V(P')={\mathbb{S}}$. The set $\mathcal K(N(P))$ is the closed three-dimensional disc in $\im({\mathbb{H}})$, with center at the origin and radius 2. The Gauss-Lucas snail
$\mk{sn}(P)$ is the subset of $\im({\mathbb{H}})$ obtained by rotating around the $i$-axis the following subset of the coordinate plane $L=\{x=x_1i+x_2j\in\im({\mathbb{H}})\,|\,x_1,x_2\in{\mathbb{R}}\}$:
\[\{x=\rho\cos(\theta)i+\rho\sin(\theta)j\in L\;|\;0\le\theta\le\pi,\,0\le\rho\le2\cos(\theta/3)\}.\]
Therefore $\mk{sn}(P)$ is a proper subset of $\mathcal K(N(P))$ (the boundaries of the two sets intersect only at the point $2i$). Its boundary is obtained by rotating a curve that is part of the \emph{lima\c con trisectrix} (see Figure \ref{fig:limacon}).
\end{remark}

\begin{figure}
\begin{center}
\includegraphics[width=6cm]{GaussLucas.png}\\
\caption{Cross-sections of $\mk{sn}(P)$ (gray), of $V(P')$ and $\mathcal K(N(P))$ (dashed).
\label{fig:limacon}
\end{center}
\end{figure}

\subsection{Estimates on the norm of the critical points}

Let $p(z)=\sum_{k=0}^da_kz^k$ be a complex polynomial of degree $d\ge1$. The norm of the roots of $p$ can be estimated making use of the norm of the coefficients $\{a_k\}_{k=0}^d$ of $p$. There are several classic results in this direction
(see e.g.\ \cite[\S8.1]{RahmanSchmeisser}). For instance the estimate \cite[(8.1.2)]{RahmanSchmeisser} (with $\lambda=1,p=2$) asserts that
\begin{equation}\label{eq:cauchy}
\textstyle
\max_{z\in V(p)}|z|\leq|a_d|^{-1}\sqrt{\sum_{k=0}^d|a_k|^2}\;.
\end{equation}

\begin{proposition}\label{pro:sn-compact}
For every monic polynomial $P\in{\mathbb{H}}[X]$ of degree $d\geq2$, the Gauss-Lucas snail $\mk{sn}(P)$ is a compact subset of ${\mathbb{H}}$.
\end{proposition}
\begin{proof}
Since $P=\sum_{k=0}^dX^ka_k$ is monic, every polynomial $P^I$ is monic. From \eqref{eq:cauchy} it follows that $\max_{x\in V(P^I)}|x|^2\le\sum_{k=0}^d|\pi_I(a_k)|^2\le\sum_{k=0}^d|a_k|^2$ and hence $\mk{sn}(P)\subset \{x\in{\mathbb{H}}\,|\, |x|^2\le \sum_{k=0}^d|a_k|^2\}$ is bounded. Since $\mk{sn}(P)$ is closed in ${\mathbb{H}}$, as seen in Remark \ref{rem:closed}, it is also a compact subset of ${\mathbb{H}}$.
\end{proof}

Define a function $C:{\mathbb{H}}[X]\to{\mathbb{R}}\cup\{+\infty\}$ as follows: $C(a):=+\infty$ if $a$ is a quaternionic constant and
\[
\textstyle
C(P):=|a_d|^{-1}\sqrt{\sum_{k=0}^d|a_k|^2} \qquad \text{if $P(X)=\sum_{k=0}^dX^ka_k$ with $d\ge1$ and $a_d\neq0$.}
\]

\begin{proposition}\label{pro:estimate}
For every polynomial $P\in{\mathbb{H}}[X]$ of degree $d\geq1$, it holds
\begin{equation}\label{eq:C}
\max_{x\in V(P)}|x|\leq C(P).
\end{equation}
\end{proposition}
\begin{proof}
We follow the lines of the proof of estimate \eqref{eq:cauchy} for complex polynomials given in \cite{RahmanSchmeisser}. Let $P(X)=\sum_{k=0}^dX^ka_k$ with $d\ge1$ and $a_d\neq0$.
We can assume that $P(X)$ is not the monomial $X^d a_d$, since in this case the thesis is immediate.
Let $b_k=|a_ka_d^{-1}|$ for every $k=0,\ldots,d-1$. The real polynomial
$h(z)=z^d-\sum_{k=0}^{d-1}b_kz^k$ has exactly one positive root $\rho$ and is positive for real $z>\rho$ (see \cite[Lemma~8.1.1]{RahmanSchmeisser}). Let $S:=\sum_{k=0}^{d-1}b_k^2=C(P)^2-1$. From the Cauchy-Schwartz inequality, it follows that
\[\left(\sum_{k=0}^{d-1}b_kC(P)^k\right)^2\le
S\sum_{k=0}^{d-1}C(P)^{2k}=
(C(P)^2-1)\frac{C(P)^{2d}-1}{C(P)^2-1}<C(P)^{2d}.
\]
Therefore $h(C(P))>0$ and then $C(P)>\rho$.  Let $x\in V(P)$.
It remains to prove that $|x|\le\rho$. Since $x^d=-\sum_{k=0}^{d-1}x^ka_ka_d^{-1}$, it holds
\[|x|^d\le\sum_{k=0}^{d-1}|x|^k\left|a_ka_d^{-1}\right|=\sum_{k=0}^{d-1}|x|^kb_k.\]
This means that $h(|x|)\le0$, which implies $|x|\le\rho$.
\end{proof}

From Proposition~\ref{pro:estimate} it follows that for every polynomial $P\in{\mathbb{H}}[X]$ of degree $d\geq2$, it holds
\begin{equation}\label{eq:C}
\max_{x\in V(P')}|x|\leq C(P').
\end{equation}
Theorem \ref{thm} allows to obtain a new estimate.

\begin{proposition}
Given any polynomial $P\in{\mathbb{H}}[X]$ of degree $d\geq2$, it holds:
\begin{equation}\label{eq:ijk'}
\max_{x\in V(P')}|x|\leq\sup_{I\in{\mathbb{S}}}\{C(P^I)\}.
\end{equation}
\end{proposition}
\begin{proof}
If $x\in V(P')\cap{\mathbb{C}}_I$, Theorem \ref{thm} implies that $x\in\mathcal K_{{\mathbb{C}}_I}(P^I)$. Therefore

\begin{equation*}\label{eq:ijk}
\max_{x\in V(P')\cap{\mathbb{C}}_I}|x|\leq C(P^I) \quad \text{ for every $I\in{\mathbb{S}}$ with $V(P')\cap{\mathbb{C}}_I\ne\emptyset$},
\end{equation*}
from which inequality \eqref{eq:ijk'} follows.
\end{proof}

Our estimate \eqref{eq:ijk'} can be strictly better than classic estimate \eqref{eq:C}, as explained below.

\begin{remark}
Let $d\geq3$ and let $P(X)=X^{d-3}\cdot(X-i)\cdot(X-j)\cdot(X-k)$. Using \eqref{P1}, by a direct computation we obtain:
\[
C(P')=d^{-1}\sqrt{8d^2-24d+24}.
\]
Moreover, given $I=\alpha_1i+\alpha_2j+\alpha_3k\in{\mathbb{S}}$ for some $\alpha_1,\alpha_2,\alpha_3\in{\mathbb{R}}$ with $\alpha_1^2+\alpha_2^2+\alpha_3^2=1$, we have
$ \pi_I(i+j+k)=\langle I,i+j+k\rangle I=(\alpha_1+\alpha_2+\alpha_3)I$ and
$\pi_I(i-j+k)=\langle I,i-j+k\rangle I=(\alpha_1-\alpha_2+\alpha_3)I$
and hence
\[
 C(P^I)=\sqrt{4+4\alpha_1\alpha_3}\leq\sqrt{4+2(\alpha_1^2+\alpha_3^2)}\leq\sqrt{6}.
\]
This implies that
\[
\sup_{I\in{\mathbb{S}}}\{C(P^I)\}\leq\sqrt{6}.
\]
For every $d\geq11$ it is easy to verify that $\sqrt{6}<C(P')$ so
\[
\sup_{I\in{\mathbb{S}}}\{C(P^I)\}<C(P'),
\]
as announced.
\end{remark}

\begin{remark}
Some of the results presented here can be generalized to real alternative *-algebras, a setting in which polynomials can be defined and share many of the properties valid on the quaternions (see \cite{GhPe_AIM}).
The polynomials given in Corollary \ref{counterexample} can be defined every time the algebra contains an Hamiltonian triple $i,j,k$. This property is equivalent to say that the algebra contains ${\mathbb{H}}$ as a subalgebra (see \cite[\S8.1]{Numbers}). For example, this is true for the algebra of octonions and for the Clifford algebras with signature $(0,n)$, with $n\ge2$. Therefore in all such algebras there exist polynomials for which the zero set $V(P')$ (as a subset of the \emph{quadratic cone}) is not included in the circularization of the convex hull of $V(N(P))$ viewed as a complex polynomial.
\end{remark}

The quaternionic Gauss-Lucas Theorem