Zeros of fractional derivatives of polynomials
We investigate the behavior of fractional derivatives of polynomials. In particular, we consider the locations and the asymptotic behaviour of their zeros and give bounds for their Mahler measure.
\section{Introduction}

Questions concerning finding exact or approximate values of the zeros of polynomial functions $p(x) = c_n x^n + c_{n-1}x^{n-1} + \cdots + c_1 x + c_0$ are classical, and (for the case of real coefficients $c_0, c_1, \cdots c_n$) properties of the distribution of these zeros have been studied since at least 1637, when Descartes established his fundamental Rule of Signs (in {\em La G{\' e}om{\' e}trie} \cite{descartes1637}). This important result was refined to finite intervals by Budan in 1807 and by Fourier in 1820 (see \cite{bf1} or \cite{bf2}). By then, thanks to the work of Euler, Gauss and Argand (see \cite{fine1997} for more details), the Fundamental Theorem of Algebra had been established, guaranteeing that a polynomial of degree $n$ has exactly $n$ complex zeros (counted with multiplicity). Not much later,
in 1829, Cauchy \cite{cauchy1829} was able to prove that all zeros of a monic polynomial $p(x)$, with complex coefficients, must lie inside the disk $|z| < 1 + \max_{0 \leq k \leq n-1} |c_k|$, 
the bounds that were eventually generalized by Landau \cite{landau1907}, Fej{\' e}r \cite{fejer1908} and others.  

In another direction,
one could ask about the relation between locations of the zeros of a polynomial $p(x)$ and the zeros of its derivative $p^{\prime}(x)$, as Rolle has done in his {\em  Trait{\' e} d’alg{\`e}bre} of 1690 (see \cite{shain1937}); the well-known theorem bearing his name -- that states that between any two zeros of a real polynomial there lies at least one zero of the derivative -- was proved rigorously by Cauchy \cite{cauchy1823} in 1823. In the complex plane, the situation becomes even more interesting. As Gauss noted in 1836, all zeros of $p^{\prime}(x)$ lie in the convex hull of the zeros of $p(x)$.
The first proof of this proposition was published by Lucas \cite{lucas1874} in 1874; it is now known as the Gauss-Lucas Theorem
(also see \cite{marden1966}). At the beginning of the 20th century it was refined by B{\^o}cher \cite{bocher1904}, Jensen \cite{jensen1912} and Walsh \cite{walsh1920}, and in more recent times, several related extensions and generalizations of it have been considered by Dimitrov \cite{dimitrov1998}, Brown \& Xiang \cite{brown1999}, Sendov \cite{sendov2021}, Tao \cite{tao2022} and others.

 

The main aim of our work is to investigate connections between these two central themes. We will try to show that their key ideas can be combined in a very natural way, but to quite surprising effects, if one considers the fractional 
derivatives $p^{(\alpha)} (x)$, where $\alpha \in \mathbb{R}$ is a variable $0 \leq \alpha \leq n = \deg p(x)$. 
Our main goal in this paper will be to answer one of the most intriguing questions that arises as soon as one begins to study these topics: since obviously $\deg p^{\prime}(x) =  \deg p(x) - 1$, and the Fundamental Theorem asserts than the same reduction must occur for the total number of zeros, 
what happens to the zeros of fractional derivatives, as the real $\alpha$ increases continuously from $0$ to $n$? How 
do the zeros of polynomials vanish, and why? As it turns out, these questions have remarkably simple and elegant answers. Namely, for a polynomial $p(x)$ of degree $n$, each of its $n$ zeros will belong to a path of unique length that connects it to the origin, where the ``length'' of the path can be measured by the number of zeros of its derivatives it contains; in other words, for each $0 \leq k \leq n-1$ there will be a unique path (originating at one of the zeros of $p(x)$) that will contain exactly $k$ zeros of its higher derivatives. (Figure \ref{fig two} shows this general property for a generic cubic polynomial).  

 

\begin{figure}[ht]
\includegraphics[width=0.45\textwidth]{paper-deg-3-center-0.pdf}
\includegraphics[width=0.45\textwidth]{paper-deg-3-center+1-I.pdf}
\caption{Paths \(z(\alpha)\) of zeros of the Riemann-Liouville fractional derivatives 
\(\RL{0}{\alpha}p\) and
\(\RL{1-i}{\alpha}p\)
of the polynomial \(p(x)=x^3 + (1+i)x^2 + x - i\).
We end the paths when \(z(\alpha)\) reaches the ``origin'' (\(a=0\) in the first case, and \(a=1-i\) in the second).
}\label{fig two}
\end{figure}

Another goal of this paper will be to try to understand some of the particulars of the paths the fractional zeros take, their dynamical properties. In order to state our results concerning this general flow of polynomial zeros more precisely, first we will need to recall some basic definitions and properties 
of fractional derivatives. There exist a multitude of different definitions of fractional derivatives, each with its own particular advantages and disadvantages. Unlike in 
our study of the Riemann zeta function and the Stieltjes constants where the Gr\"unwald-Letnikov fractional derivative
\cite{g:1867,l:1869-1} worked 
really well (see \cite{fps} and \cite{fps2}), in the case of polynomials their divergence forces us instead to employ 
the  Riemann-Liouville fractional derivative \cite{riefrac}, which also can be thought of as a truncated version of the Gr\"unwald-Letnikov fractional derivative \cite{ortigueira2011fractional}. 


The rest of the paper is structured as follows. Our Section \ref{sec derv} contains a short introduction to the Riemann-Liouville differintegral  focusing on results most applicable to fractional derivatives of polynomials. 
In Section \ref{sec low} this theory is applied to the two simplest cases: polynomials of degree one and two. These are the two cases where, thanks to the manageable classical formulas for the zeros, all the main questions can be conclusively answered. With the cubic polynomials things become somewhat murky, but general convergence trends can still be established. In Section \ref{sec flow} we do just that: we examine how the zeros of integral derivatives are connected to the zeros of fractional derivatives in the most general setting, and we look at the paths of zeros and investigate their convergence and the overall flow. In Section \ref{sec bounds} we consider the behaviour of the zeros on a larger scale and we prove bounds for
the Mahler measure of the fractional derivatives, which are then, in Section \ref{caputo}, also established for the Caputo fractional derivative. Finally, in Section \ref{conclusion}, we discuss some intriguing open problems and unsolved questions. 

\section{Riemann-Liouville Fractional 
Derivatives}\label{sec derv}
In full generality, the \(\alpha\)-th Riemann-Liouville fractional derivative is defined as follows, see \cite{oldham-spanier}, for example: 

\begin{definition}[Riemannn-Liouville fractional derivative]\label{def rl}
Let \(f\) be analytic on a convex open set \(C\) let and \(a\in C\). 
For \(\alpha>0\) the Riemann-Liouville integral is
\[
\RLI{a}{\alpha} f(t)=\frac{1}{\Gamma(\alpha)}
\int_a^t f(t)(t-\tau)^{\alpha-1} \,d\tau.
\]
Set \(\RLI{a}{0}f(t)=f(t)\) and for \(\alpha>0\) define the Riemann-Liouville fractional derivative as
\[
\RL{a}{\alpha} f(t) = \frac{d^m}{dt^m} \RLI{a}{m-\alpha} f(t),
\]
where $m=\lceil \alpha \rceil$.  For \(\alpha<0\) set
\(\RL{a}{\alpha}f(t)=\RLI{a}{-\alpha} f(t)\).
\end{definition}
In what follows, we will only consider the special case of polynomials, composed of the simple power functions $p(x)=(x-a)^{\beta}$, where $\beta \in \mathbb{R}$, $a\in\mathbb{C}} \newcommand{\N}{\mathbb{N}$. For these, the $\alpha$-th Riemann-Liouville fractional derivative can be computed 
using the Power Rule: 
\begin{align}\label{eq poly derv}
\RL{a}{\alpha}(x-a)^{\beta}=\begin{cases}
0 & \text{if }\alpha-\beta-1 \in \N\\
\dfrac{\Gamma(\beta+1)}{\Gamma(\beta-\alpha+1)}(x-a)^{\beta-\alpha} & \text{otherwise}
\end{cases}
\end{align}
where, for $z \in \mathbb{C}$, the gamma function $\Gamma(z)$ is defined as $\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}dt$; it satisfies $\Gamma(1+n)=n\Gamma(n)$, which implies that $\Gamma(n) = (n-1)!$ for $n \in \mathbb{N}$. The function $\Gamma(z)$ has no zeros in the complex 
plane $\mathbb{C}$, and has poles at all the negative integers (see \cite{artin1964}).

\begin{remark}
The Riemann-Liouville fractional derivative of a monomial \(f(x)=x^n\) is multivalued.  When changing the branch of the complex logarithm in the computation of the fractional derivative all coefficients of the derivative are changed by the same factor.
So choosing a different branch of the complex logarithm does not change the zeros of the derivative, which means that we
can fix the branch in our consideration of zeros of derivatives of polynomials.
We use the principal branch of the complex logarithm.
\end{remark}

It should be noted that the constant $a$ that centers the expansion (\ref{eq poly derv}) plays a key role in all our computations below, as the ``origin,'' or the limit of convergence, of the flow of zeros of derivatives (Figure \ref{fig two} illustrates its role). Also noteworthy is the fact that the Riemann-Liouville fractional derivative satisfies all properties expected of a regular derivative, with the exception of the composition rule. The following example shows why it fails: 

\begin{example} 
From (\ref{eq poly derv}),  the \(1.5\)-th derivative of \(p(x)=1\) is
\(
\textstyle\RL{0}{1.5}p(x)=\frac{2}{ \csqrt{\pi}}x^{-1.5}
\)
and 
\(\textstyle\RL{0}{1}\left(\RL{0}{0.5}p(x)\right) =\frac{2}{ \csqrt{\pi}}x^{-0.5}.
\)
However, \(\RL{0}{0.5}\left(\RL{0}{1}p(x)\right)=\RL{0}{0.5}0=0.
\)
\end{example}

When \(\beta\in\mathbb{R}} \newcommand{\Z}{\mathbb{Z}\setminus\Z\) we still have
\(\RL{a}{\alpha}\left(\RL{a}{\beta}p(x)\right)=\RL{a}{\alpha+\beta}p(x)\).



\begin{remark} It is possible to go beyond the standard values of $0 \leq \alpha$, and consider what happens for $\alpha < 0$. Here, there will be \(n+1\) complex roots, because the first term in Equation \ref{eq lem derv} below has the root \(0\). Just like in the standard case, the extended curves \(z(\alpha)\) of zeros of the differintegral \(\RL{a}{\alpha}p\) will be continuous for \(\alpha<0\) unless
\(\RL{a}{\alpha}p\) has a double root; however, they will not be smooth at integral \(\alpha > 0\).  More on this will be said in Section \ref{sec flow} below.
\end{remark}




In what follows, we consider the zeros of the the fractional derivatives of polynomials \(p\in\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\) of degree \(n\), and we investigate the implicit functions \(z:(0,n)\to \mathbb{C}} \newcommand{\N}{\mathbb{N}\) given by
\[
\RL{a}{\alpha}f(z(\alpha))=0.
\]
If \(\left(\RL{a}{\alpha}p(x)\right)'\ne 0\), for \(\alpha\in(0,n)\), then 
\(z(\alpha)\) is differentiable on \((0,n)\).
We denote the roots of the polynomial \(p(x)\) by \(z_1, z_2, \dots,z_n\)
and for \(1\le k\le n\) we define the the implicit function \(z_k:[0,n)\to\mathbb{C}} \newcommand{\N}{\mathbb{N}\)  by \(z_k(0)=z_k\)
and \(\RL{a}{\alpha}p(z(\alpha))=0\).


 

The following representation of the fractional derivatives of a general monic polynomial will be most useful.

\begin{lemma}\label{lem derv}
Let \(p(x)=(x-a)^n+\sum_{j=0}^{n-1} c_j (x-a)^{j}\in\mathbb{C}[x]\) and \(\alpha\not\in\N\).  Then
\begin{equation}\label{eq lem derv}
\RL{a}{\alpha}p(x)=\frac{n!}{\Gamma(n+1-\alpha)} (x-a)^{-\alpha}
\left[(x-a)^{n}+
\sum_{j=0}^{n-1} 
\left(
\prod_{k=j+1}^n\!(k-\alpha)\right)\cdot
\frac{j!}{n!} c_j (x-a)^{j}
\right].
\end{equation}
\end{lemma}


\begin{proof}
Applying the Power Rule (\ref{eq poly derv}), we obtain \begin{align*}
\RL{a}{\alpha}p(x)  &=
\frac{\Gamma(n+1)}{\Gamma(n+1-\alpha)} (x-a)^{n-\alpha}+
\sum_{j=0}^{n-1} \frac{\Gamma(j+1)}{\Gamma(j+1-\alpha)} c_j (x-a)^{j-\alpha}\\
&=\frac{\Gamma(n+1)}{\Gamma(n+1-\alpha)} (x-a)^{-\alpha}
\left[(x-a)^{n}+
\sum_{j=0}^{n-1} 
\frac{\Gamma(n+1-\alpha)}{\Gamma(j+1-\alpha)}
\frac{\Gamma(j+1)}{\Gamma(n+1)} c_j (x-a)^{j}
\right]
\\
&=\frac{n!}{\Gamma(n+1-\alpha)} (x-a)^{-\alpha}
\left[(x-a)^{n}+
\sum_{j=0}^{n-1} 
\left(
\prod_{k=j+1}^n\!(k-\alpha)\right)\cdot
\frac{j!}{n!} c_j (x-a)^{j}
\right],
\end{align*}
as wanted. \qedhere
\end{proof}

\begin{remark} The representation of the fractional derivatives of polynomials given in the above  Lemma \ref{lem derv} has the property that their roots only depend on the factors in square brackets, which in turn implies the useful fact that the branch cut of the complex logarithm does not affect the paths \(z(\alpha)\) of
zeros of these fractional derivatives.

We also see that for \(p(x)=(x-a)^n\) we have:
\begin{equation}\label{eq at a}
\RL{a}{\alpha}p(a)=
\begin{cases}
0 &\mbox{if }\alpha<n\\
\Gamma(n+1) & \mbox{if }\alpha=n\\
\mbox{undefined} &\mbox{if }\alpha>n
\end{cases}
\end{equation}
\end{remark}


\begin{remark} A few words should be said about our plots of the implicit functions \(z:\mathbb{R}} \newcommand{\Z}{\mathbb{Z}\to\mathbb{C}} \newcommand{\N}{\mathbb{N}\)
with \(z(\alpha)\) given by \(\RL{a}{\alpha}p(z(\alpha))=0\).  
The dots labelled `\(\bullet\)\textsf{k}' represent zeros of the \(k\)th Riemann-Liouville differintegral (thus `\(\bullet \)\textsf{0}' 
represent the zeros of the polynomial \(p(x)\) itself), while circles
`\(\circ\)\textsf{k}' represent points that are limits of \(z(\alpha)\) as \({\alpha\to k}\) but are not zeros of \(\RL{a}{k}p\).  These occur for integral \(k\) with \(k\ge\deg p(x)\), where \(\RL{a}{k}p(x)\) is constant, for example at \(x=a\), see Equation (\ref{eq at a}).
Moreover, when a point is either a zero or a limit point of zeros of both the \(j\)th and the \(k\)th differintegrals, then it is represented by `\(\bullet\) \textsf{j \& k}' or `\(\circ\) \textsf{j \& k}' respectively.  
In Figures \ref{fig two}, \ref{fig double}, and \ref{fig longest} we let all paths of zeros of \(\RL{a}{\alpha}p(x)\) end when they reach the origin \(a\).  In Figures \ref{fig one},  \ref{fig double}, \ref{fig p1},
\ref{fig quintic paths}, and \ref{fig quintic abs}
we continue the paths past \(a\). 
In Figures \ref{fig one}, \ref{fig double}, and \ref{fig p1} we display the path of zeros \(z(\alpha)\) for \(\alpha<0\) and \(\alpha>n\) in lighter colors than for \(0\le \alpha < n \) where \(n\) is the degree of the polynomial.  
The point \(a\) is only labeled with values for \(\alpha\) for \(\alpha\ge 0\).
\end{remark}



\begin{figure}
\includegraphics[width=0.45\textwidth]{paper-deg-1.pdf}
\includegraphics[width=0.45\textwidth]{paper-deg-2.pdf}
\caption{Path $z(\alpha)$ of zeros of differintegrals of the polynomials \(p(x)=x-2-3i\) and $p(x)=(x-2-3i)(x+2+i)$.}
\label{fig one}
\end{figure}


\section{Low Degrees}\label{sec low}
Formulas for finding zeros of polynomials of low degrees have been known for centuries. Applying the Riemann-Liouville derivative to these low-degree cases
has proved to be a simple but informative exercise. In this section we summarize some of these results, stating the most useful ones as lemmas. They are 
examples of a 
dynamic that shares certain key characteristic with most high-degree cases, but some aspects of which are often unique. For example, in the linear case, 
the path the zeros take is also linear, while already in the quadratic case one observes a considerably more complex behavior.

Let us start with the linear polynomials. Here the situation is simple. The paths of zeros will always be linear, too, and they can be completely described. We have: 

\begin{lemma}[Linear Polynomials]\label{lem linear}
Let \(p(x)=(x-a)+c_0\). As $\alpha$ increases from $0$ to $1$, the path of zeros of \(\RL{a}{\alpha}p(x)\) is given by \(z(\alpha)=(\alpha-1)c_0+a\).
\begin{proof}
From (\ref{eq poly derv}), with $\beta = 1$, we get: \(\RL{a}{\alpha}p(x)=(x-a)^{-\alpha}\left((x-a)+(1-\alpha)c_0\right)\).
\end{proof}
\end{lemma}

\begin{remark}
    In addition to considering the $\alpha$-th derivatives in the usual range $0 \leq \alpha \leq n$, one could also look at what happens when $\alpha < 0$ and when $\alpha>n$.   
    In the linear case this is, again, simple.   From Lemma \ref{lem linear} we can deduce that, as with $\alpha<0$, the line of roots of the derivatives continues. Similarly, for $\alpha>1$, $\alpha\not\in\Z$ lie on the same line, see Figure \ref{fig one} below. 
\end{remark}











Let us now consider the quadratic case. This case is considerably more interesting, since there are
now two paths of zeros of the fractional derivatives, and they exhibit a much more complex and intricate behavior.
We first notice that the path of the zero closest to $a$ will directly connect with $a$, while the path of the farther zero will in the process of reaching $a$ pass through the zero of the first derivative. 

\begin{lemma}\label{lem quad root order}
Let $p(x)=x^2+c_1(x-a)+c_0$ with discriminant \(d=\) and roots \(s_1\) and \(s_2\).
\begin{enumerate}
\item If \(d\in\mathbb{C}} \newcommand{\N}{\mathbb{N}\setminus\mathbb{R}} \newcommand{\Z}{\mathbb{Z}^{<0}\), denote by \( \csqrt{d}\) the 
complex number $r$ with \(r^2=d\) and \(\Re(r)> 0\).  Then for the 
roots \(s_1=\frac{-b+ \csqrt{b^2-4c}}{2}\) and \(s_2=\frac{-b- \csqrt{b^2-4c}}{2}\)
of \(g(x)=x^2+bx+c\in\mathbb{C}[x]\) we have \(|s_1|\le |s_2|\).
\item If \(d\in\mathbb{R}} \newcommand{\Z}{\mathbb{Z}^{<0}\), then \(|s_1|=|s_2|\).
\end{enumerate}
\end{lemma}

\begin{proof}
\begin{enumerate}
\item Write \(d=1-\frac{4c}{b^2}\).
Then
\(|s_1|=\left|\frac{b}{2}\right|\cdot \left|-1+ \csqrt{d}\right|\)
and \(|s_2|=\left|\frac{b}{2}\right|\cdot \left|-1- \csqrt{d}\right|\).
Now
\begin{align*}
\left|-1+ \csqrt{d}\right|^2 & = \left(-1+ \csqrt{d}\right)\left(-1+\overline{ \csqrt{d}}\right)=1- \csqrt{d}-\overline{ \csqrt{d}}+|d| \text{ and }\\
\left|-1- \csqrt{d}\right|^2 &= \left(-1- \csqrt{d}\right)\left(-1-\overline{ \csqrt{d}}\right)=1+ \csqrt{d}+\overline{ \csqrt{d}}+|d|
\end{align*}
Because \(\Re\left( \csqrt{d}\right)> 0 \) we have
\(|s_1|\le |s_2|\). 
\item Now assume $\Re(r)=0$, meaning $r$ is completely imaginary, i.e. $r=it$. Then using the information above we have, 
\begin{align*}
    \left|-1+ \csqrt{d}\right|^2 & =1- \csqrt{d}-\overline{ \csqrt{d}}+|d| = 1-r-\overline{r}+|r^2| =1-it-\overline{it}+\left|(it)^2\right| \\
    & =1-t^2 =1+it-it-t^2 =1+it+\overline{it}-t^2 =\left|-1- \csqrt{d}\right|^2.
\end{align*}
Therefore, we have \(|s_1|= |s_2|\).\qedhere
\end{enumerate}
\end{proof}

\begin{remark}
    As of yet, there is no known reliable ordering of the zeros for any of the higher degree polynomials. In fact, there exist examples of cubic polynomials for which the standard Euclidean distance (which works so well for the linear and the quadratic cases) can be shown to fail: see Figure \ref{fig longest}.  
\end{remark}


 

 In addition to the natural ordering on the quadratic roots, another question that seems to be of interest is the one that concerns the trends of descent of their paths, especially since it had such a nice answer in the linear case. As it turns out, the asymptotes of the two quadratic paths exist, and the quadratic formula alone is enough to help us find them. 
 
  
 

\begin{theorem}\label{prop roots quad}
For the quadratic polynomial $p(x)=(x-a)^2+c_1(x-a)+c_0$, the paths of zeros of the fractional derivatives \(\RL{a}{\alpha} p(x)\) are given as
\begin{equation}\label{eq roots quad}
z_{1,2}(\alpha)=a+\frac{-(2-\alpha) c_1\pm \csqrt{(2-\alpha)^2\cdot (c_1)^2 - 8(2-\alpha)(1-\alpha) \cdot c_0}}{4},
\end{equation}
with \(|z_1(\alpha)|\ge|z_2(\alpha)|\), for {\(\alpha\in[0,1]\)}, and \(\lim_{\alpha\to 1}z_1(\alpha)=a\)
and 
\(\lim_{\alpha\to 2}z_{1,2}(\alpha)=a\).
\begin{proof}
With the help of (\ref{eq poly derv}), the fractional derivatives of \(p(x)\) can be written as 
\[
\RL{a}{\alpha}p(x)=\frac{2}{\Gamma(3-\alpha)}x^{-\alpha}
\left(
(x-a)^{2} +\frac{(2-\alpha)\cdot c_1}{2} (x-a) +\frac{(2-\alpha)(1-\alpha)\cdot c_0}{2}\right).
\]
Now, set $y=x-a$. Then the roots of $y^{2} +\frac{(2-\alpha)\cdot c_1}{2} y +\frac{(2-\alpha)(1-\alpha)\cdot c_0}{2}$ . 
are 
\[
z_{1,2}(\alpha)=a+\frac{-(2-\alpha) c_1\pm \csqrt{(2-\alpha)^2\cdot (c_1)^2 - 8(2-\alpha)(1-\alpha) \cdot c_0}}{4}
\]
The ordering of the roots \(|z_1(\alpha)|\ge|z_2(\alpha)|\), for {\(\alpha\in[0,1]\)},  follows with Lemma \ref{lem quad root order}. Furthermore, we have 
\[
\lim_{\alpha\to 1} z_{1,2}(\alpha)=a+\frac{- c_1\pm \csqrt{ (c_1)^2}}{4}=a+\frac{-c_1\pm c_1}{4}.
\]
Thus $\lim_{\alpha\to 1} z_1(\alpha)=a$ and $\lim_{\alpha\to 1} z_2(\alpha)= a+\frac{c_1}{2}$.
Similarly
\[    
\lim_{\alpha\to 2}\left[a+\frac{-(2-\alpha) c_1\pm \csqrt{(2-\alpha)^2\cdot (c_1)^2 - 8(2-\alpha)(1-\alpha) \cdot c_0}}{4}\right]=a.\qedhere
\]
\end{proof}
\end{theorem}

 

\begin{corollary}
For \(\alpha \to \pm\infty\), the asymptotes of the quadratic paths are
\[\textstyle z_{1,2}(\alpha)\approx\frac{(2-\alpha) c_1}{4}\left[-1\pm \csqrt{1-\frac{8c_0}{c_1}}\right]. \]
\end{corollary}

\begin{proof}
By (\ref{eq roots quad}) 
we have,   
\begin{align*}
        z_{1,2}(\alpha)&=\frac{-(2-\alpha) c_1\pm(2-\alpha) c_1 \csqrt{1 - 8\frac{(1-\alpha)}{(2-\alpha)} \cdot \frac{c_0}{c_1}}}{4}\\
        &=\frac{(2-\alpha) c_1}{4}\left[-1\pm \csqrt{1-\frac{8c_0(1-\alpha)}{c_1^2(2-\alpha)}}\right],\\
\end{align*}
and since $\frac{1-\alpha}{2-\alpha}\to 1$, as $\alpha\to \pm\infty$, this yields linear asymptotes for $z_1(\alpha)$ and $z_2(\alpha)$.
\end{proof}



An noteworthy special case occurs when the polynomial has a double root. Then the paths display an interesting symmetry, see Figure \ref{fig double}.
In fact, it is easy to see that specializing our Proposition \ref{prop roots quad} to the case of a double zero of the polynomial itself yields:

\begin{corollary}
If $p(x)=(x-z_0)^2$ then the zeros of \(\RL{0}{\alpha}p\) are
\[z_{1,2}(\alpha)=\frac{z_0\left(-(2-\alpha) \pm  \csqrt{\alpha(2-\alpha)}\right)}{2}. \]
\end{corollary}




\begin{figure}
\includegraphics[width=0.45\textwidth]{paper-deg-2-double-totale.pdf}
\includegraphics[width=0.45\textwidth]{paper-deg-2-double-root-abs.pdf}
\caption{Paths \(z(\alpha)\) of zeros of the fractional derivatives of 
\(p(x)=x^2+(2+6i)x-12+9i\) where \(\RL{0}{1/2}p\) has a double root.  In the plot on the right the wide, light green graph represents the real part of \(z(\alpha\)), while the thin, dark red graph represents its imaginary part. 
}\label{fig double}
\end{figure}


Another natural question to ask is whether, given that a quadratic polynomial has  distinct zeros, can its fractional derivative have a double zero. Setting \(z_1(\alpha)=z_2(\alpha)\) one gets:

\begin{corollary}\label{lem double zero}
Let $p(x)=(x-a)^2+c_1(x-a)^1+c_0(x-a)^0$. Then the fractional derivative \(\RL{a}{\alpha}(p(x))\) will have a double zero precisely for one \(\alpha\in\mathbb{R}} \newcommand{\Z}{\mathbb{Z}\setminus\N\), namely: 
\(\alpha = 1-\frac{c_1^2}{8c_0- c_1^2}\).
\end{corollary}




\section{Flow of Zeros}\label{sec flow}

As stated above, one of our main goals was to consider the paths of zeros of the fractional derivatives of polynomials \(p(x) \in\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\) of arbitrary degrees. Unfortunately, unlike in the linear and quadratic cases, already for the cubics we find that the situation becomes considerably more complicated. This can be seen from the fact that one of the nicest properties -- the natural ordering of zeros -- fails already for degree 3: in other words, 
it is not true in general that zeros furthest away from the origin will yield the longest paths of zeros of fractional derivatives on its way to the origin. Figure
\ref{fig longest} shows a notable counterexample.

\begin{figure}
\includegraphics[width=0.45\textwidth]{paper-furthest-longest-counterexample-paths.pdf}
\includegraphics[width=0.45\textwidth]{paper-furthest-longest-counterexample-abs.pdf}
\caption{The paths \(z(\alpha)\) of zeros of the fractional derivatives of 
the cubic \(p(x)=(x-3-2i)(x-2+5i)(x-4+4i)\)
and the absolute values \(|z(\alpha)|.\)
In the latter case, the zero of \(p(x)\) with the greatest absolute value is
not the starting point of the longest path.
}
\label{fig longest}
\end{figure}


However, certain convergence properties of the paths can be established in general. For example, the following theorem shows that all the paths terminate in the origin $a$. 

\begin{theorem}\label{theo roots lim}
Let \(p(x) \in \mathbb{C}[x]\) of degree \(n\) such that for all \(\alpha\in [0,n]\)
the fractional derivative \(p^{(\alpha)}(x)\) has no double zeros.
Then there is an ordering of the roots \(z_1,\,z_2,\,...,\,z_n\) or \(p\) such that for 
\(p^{(\alpha)}(z_j(\alpha)) = 0\) and \(z_j(0)  = z_j\) for \(j\in\{1,...,n\}\) we have  
\(\lim_{\alpha \to j} z_j(\alpha) = a\)
\end{theorem}





\begin{proof}
We  denote the coefficients in the expansion proved in Lemma \ref{lem derv} by
\begin{equation}
d^\alpha_j:=\prod_{k=j+1}^n\!(k-\alpha)\cdot \frac{j!}{n!}. \label{eq dalphaj}
\end{equation}
Let \(0\le j\le n\) and \(m>j\).  Here, clearly
\[
\lim_{\alpha\to m} 
d^\alpha_j=
\lim_{\alpha\to m} 
\prod_{k=j+1}^n\!(k-\alpha)\cdot \frac{j!}{n!}=0. 
\]
Write the coefficients of the derivatives  \(\RL{a}{\alpha}p(x)\) by
\(d^\alpha_j\,c_j\) 
as symmetric functions of the roots 
\(z_1(\alpha),z_2(\alpha), \dots, z_n(\alpha)\) of  \(\RL{a}{\alpha}p(x)\).
Because 
\[
0=\lim_{\alpha\to m} d^\alpha_0\,c_0= \lim_{\alpha\to m} \prod_{k=1}^n z_k(\alpha)
\]
we have
\(\lim_{\alpha\to m} z_{k_1}(\alpha)=0\), for at least one \(1\le k_1\le n\).
Inductively, continuing with the next coefficient we get:
\begin{equation}\label{eq sym 1}
0=\lim_{\alpha\to m} d^\alpha_1\,c_1= \lim_{\alpha\to m} 
\sum_{l=1}^n \prod_{k\ne l} z_k(\alpha).
\end{equation}
There is one summand in (\ref{eq sym 1}) that does not contain \(z_{k_1}(\alpha)\).
So we need to have \(\lim_{\alpha\to m} z_{k_2}(\alpha)=0\) for at least one \(k_2\ne k_1\).
This argument also holds for all \(d^\alpha_j\,c_j\) with \(j<m\).
Therefore, we get \(\lim_{\alpha\to m} z_{k}(\alpha)=0\) for at least \(m\) distinct 
\(k\in\{1,\dots, n\}\).
\end{proof}



\begin{figure}
\includegraphics[width=0.45\textwidth]{paper-deg-3-p1-totale.pdf}
\includegraphics[width=0.45\textwidth]{paper-deg-3-p1-totale-mahler.pdf}
\caption{Paths \(z(\alpha)\) of zeros of \(\RL{0}{\alpha}(x^3+x^2+x+1+i)\).
illustrating the growth of the absolute value of \(z(\alpha)\) as
\(\alpha\to\infty\) and \(\alpha\to-\infty\) and the Mahler measure
of \(\RL{0}{\alpha}p\) along with the bounds from Theorems \ref{theo upper bound} and \ref{theo lower bound}.
}
\label{fig p1}
\end{figure}

\section{Bounds}\label{sec bounds}

As stated in the introduction, for \(f\in\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\) the Gauss-Lucas theorem states that all zeros of \(f'\) lie in the convex hull of the set of zeros of \(f\), see \cite[Theorem 6.1]{marden1966}. By induction this generalizes to all integral derivatives. Unfortunately, although all roots of the fractional derivatives converge to the origin, by our Theorem \ref{theo roots lim}, the analogue of the Gauss-Lucas theorem does not hold for the fractional derivatives.
This is an immediate consequence of a result by Genchev and Sendov
{\cite{genchev-sendov}}, which is also stated as {\cite[Theorem 2]{nikolov-sendov}}:

\begin{theorem}
Let \(L:\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\to\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\) be a linear operator, such that \(L(p)\ne 0\) implies that
the convex hull of the set of roots of \(p\) contains the roots of \(L(p)\). 
Then \(L\) is a linear functional or there are \(c\in\mathbb{C}} \newcommand{\N}{\mathbb{N}\setminus\{0\}\) and \(k\in\N\) such that \(L(p)=cp^{(k)}\).
\end{theorem}

Figure \ref{fig longest} illustrates this result by giving a specific counterexample to the Gauss-Lucas property for the case of the Riemann-Liouville fractional derivatives. Nevertheless, it is possible to make some useful statements about how the absolute values of zeros
\(z_k(\alpha)\)
of the fractional derivatives \(\RL{0}{\alpha} f\)
decrease as \(\alpha\) increases in terms of the Mahler measure of \(f\).

Let \(p(x)=x^n+\sum_{j=1}^{n-1} c_j x^j=\prod_{j=1}^n (x-z_j)\in\mathbb{C}[x]\).
For the Mahler measure \(M(f)\) \cite{mahler} we have
\[
M(p) = \exp\left(\int_0^1 \log(|f(e^{2\pi i\theta})|)\, d\theta \right)= \prod_{j=1}^n \max\{1,|z_j|\}.
\]
Denote the height of \(f\) by
\(
\nheight{p} = \max\{c_0,\dots,c_n\}
\)
and the length of \(f\) by
\(
\nlength{p}=|c_0|+\dots+|c_n|.
\)
Recall that Mahler was able to prove the bounds 
\begin{equation}\label{eq mahler height}
\binom{n}{\lfloor n/2 \rfloor}^{-1} \nheight{p} \le M(p) \le \nheight{p}  \csqrt{n+1}
\end{equation}
and
\begin{equation}\label{eq mahler length}
2^{-n} \nlength{p} \le M(p) \le  \nlength{p}.
\end{equation}
For \(\ntwo{p} = \left(\sum_{j=1}^n |c_j|^2\right)^\frac{1}{2}\) we have Landau's inequality
\cite{landau1905}
\begin{equation}\label{eq landau}
M(p)\le ||p||_2.
\end{equation}
We generalize the definition of the Mahler measure to fractional derivatives.
Let \(Z(\alpha)\) be the set of zeros of \(\RL{0}{\alpha}p\).  We set
\[
M\left(\RL{0}{\alpha}p\right)=\prod_{z\in Z(\alpha)}\max\{1,|z|\}.
\]
and prove bounds similar to Equations \ref{eq mahler height},
\ref{eq mahler length}, and \ref{eq landau}
for the fractional cases. 
We first estimate the coefficients \(d^\alpha_j\) from Proposition \ref{theo roots lim} and then use them to derive bounds for \(M\left(\RL{0}{\alpha}p\right)\).

\begin{figure}
\includegraphics[width=0.45\textwidth]{paper-mahler-paths.pdf}
\includegraphics[width=0.45\textwidth]{paper-caputo-mahler-paths.pdf}
\caption{Comparing paths \(z(\alpha)\) of zeros of the Riemann-Liouville (left) and Caputo (right) fractional derivatives 
of the  quintic 
\(p(x)=(x+1-2i)(x-3-2i)(x-2+5i)(x-3+3i)(x+1+5i)\), for \(0\le\alpha\le 5\).
}
\label{fig quintic paths}
\end{figure}

\begin{lemma}\label{lem dj}
Let \(n\in\N\) and \(\alpha\in\mathbb{R}} \newcommand{\Z}{\mathbb{Z}\setminus\N\).
Let
\(d^\alpha_j=\left(\prod_{k=j+1}^n\!(k-\alpha)\right)\cdot \frac{j!}{n!}
\) where \(0 \le j\le n-1\).  Then
\begin{enumerate}
\item\label{dj gt 0}
\(|d^\alpha_j| \le \frac{n-\alpha}{n}\), for  \(0<\alpha< n\).
\item\label{dj lt}
\(|d^\alpha_j| \le \frac{\ff{|n-\alpha|}{n}}{n!}\), 
where \(\ff{|n-\alpha|}{n}=\prod_{k=0}^{n-1} |n-\alpha-k|\),
for \(\alpha<0\) and \(\alpha>n\).
\item\label{dj lt 0 gt 2n}
\( |d^\alpha_j|\ge\frac{n-\alpha}{n}\),  for \(\alpha<0\)
and \(\alpha>2n\).
\item\label{dj gt n}
\(|d^\alpha_j|\ge\frac{\alpha-n}{n} \binom{n-1}{(n-1)/2}^{-1}\), for \(n<\alpha\le 2n\).
\end{enumerate}
\end{lemma}
\begin{proof}
\begin{enumerate}
\item
For \(0\le\alpha < j+1\), we have
\[
|d^\alpha_j| = \frac{n-\alpha}{n} \prod_{k=j+1}^{n-1}\frac{k-\alpha}{k} \cdot \frac{j!}{j!}< \frac{n-\alpha}{n}.
\]
When \(j+1 < \alpha\) set \(h:=\lfloor\alpha\rfloor\).  We get
\begin{align*}
|d^\alpha_j| 
&=
\frac{n-\alpha}{n}\cdot\frac{(n-1-\alpha)\cdots(h+1-\alpha)\cdot(\alpha-j-1)\cdots(\alpha-h)\cdot j!}{(n-1)!}
\le \frac{n-\alpha}{n}.
\end{align*}
\item For \(\alpha<0\) and \(\alpha>n\), we have
\begin{align*}
|d^\alpha_j| 
&=
\frac{|n-\alpha|\cdots|j+1-\alpha|\cdot j!}{n!}\\
&=
\frac{|\alpha-n|\cdots|\alpha-j-1|\cdot j!}{n!}\\
&\le
\frac{|\alpha-n|\cdots|\alpha-1|}{n!}=\frac{|n-\alpha|_n}{n!}.
\end{align*}
For \(\alpha<0\) and \(\alpha>2n\), we have
\begin{align*}
|d^\alpha_j| 
&=\frac{n-\alpha}{n}\cdot\frac{(n-1-\alpha)\cdots(j+1-\alpha)\cdot j!}{(n-1)!} \\
&\geq\frac{n-\alpha}{n}\cdot\frac{(n-1)\cdots(j+1)\cdot j!}{(n-1)!} 
=\frac{n-\alpha}{n}.
\end{align*}
\item
For \(\alpha>2n\), we have
\begin{align*}
|d^\alpha_j| 
=\frac{\alpha-n}{n}\cdot\frac{(\alpha-n+1)\cdots(\alpha-j-1)\cdot j!}{(n-1)!} \ge \frac{\alpha-n}{n}.
\end{align*}
\item
For \(2n>\alpha>n\), we have
\begin{align*}
|d^\alpha_j| 
&=\frac{(\alpha-n)\cdots(\alpha-j-1)\cdot j!}{n!}\\
&\ge\frac{\alpha-n}{n}\cdot\frac{1\cdots(n-j-1)\cdot j!}{(n-1)!}\\
&\ge\frac{\alpha-n}{n} \binom{n-1}{(n-1)/2}^{-1}.
\end{align*}
\end{enumerate}
\end{proof}

\begin{figure}
\includegraphics[width=0.45\textwidth]{paper-mahler-abs.pdf}
\includegraphics[width=0.45\textwidth]{paper-caputo-mahler-abs.pdf}
\caption{Comparing the absolute value of \(z(\alpha)\) of zeros of the Riemann-Liouville (left) and Caputo (right) fractional derivatives 
of the  quintic 
\(p(x)=(x+1-2i)(x-3-2i)(x-2+5i)(x-3+3i)(x+1+5i)\), for \(0\le\alpha\le 5\).
}
\label{fig quintic abs}
\end{figure}


We now show that \(M\left(\RL{0}{\alpha} p\right)\) 
is bounded from above for \(0<\alpha<\deg p\) where the bound linearly 
decreases with \(\alpha\).  Furthermore, in Theorem \ref{theo lower bound},
we see that for \(\alpha<0\) and \(\alpha>\deg p\) the Mahler measure
\(M\left(\RL{0}{\alpha} p\right)\) increases at least linearlily with \(\alpha\).

In our proofs  we use the following notation. We set
\begin{align}\label{eq equivalent poly}
\RLs{0}{\alpha}f(x) = \frac{\Gamma(n+1-\alpha)}{n!}x^\alpha\cdot \RL{0}{\alpha}f(x).
\end{align}
Because \(\RLs{0}{\alpha}f(x)\)
has the same roots as \(\RL{0}{\alpha}\) except for the additional root \(0\).
we have 
\begin{align*}
M\left(\RLs{0}{\alpha}f\right)  
= M\left(\RL{a}{\alpha}f)\right).
\end{align*}

\begin{theorem}\label{theo upper bound}
Let \(f(x)\in\mathbb{C}[x]\) be monic and let \(0<\alpha<n\).  Then
\begin{enumerate}
\item $M(\RL{0}{\alpha} f)  \le \frac{n-\alpha}{n}\nlength{f}+1$
\item $M(\RL{0}{\alpha} f)  \le  \csqrt{n+1} \max\left\{\frac{n-\alpha}{n} \nheight{f},1\right\}$
\item $M(\RL{0}{\alpha} f)  \le \frac{n-\alpha}{n} \ntwo{f}+1$
\end{enumerate}
\end{theorem}
\begin{proof}
Write \(f(x)=\sum_{j=1}^n c_j x^j\).
With the notation from Lemma \ref{lem derv} we have
\begin{align*}
\RLs{0}{\alpha}(x)
= x^{n}+ \sum_{j=0}^{n-1}
\left(
\prod_{k=j+1}^n\!(k-\alpha)\right)\cdot
\frac{j!}{n!} c_j x^{j} = x^n+\sum_{j=0}^{n-1} d^{\alpha}_j\,c_j\, x^j.
\end{align*}
This yields the following three bounds for the Mahler measure of \(\RL{a}{\alpha} f \).
\begin{enumerate}
\item
With Lemma \ref{lem dj} (\ref{dj gt 0}), we get
\begin{align*}
\nlength{\RLs{0}{\alpha}}
&=|d_0^{(\alpha)}c_0|+|d_1^{(\alpha)}c_1|+\dots+|d_{n-1}^{(\alpha)}c_{n-1}|+1\\
&\le \frac{n-\alpha}{n}\left(|c_0|+|c_1|+\dots+|c_{n-1}|\right)+1\\
&\le \frac{n-\alpha}{n} \nlength{f}+1.
\end{align*}
With Mahler's (\ref{eq mahler length}), we get
\[
M\left(\RL{a}{\alpha} f\right) =
M\left(\RLs{0}{\alpha}\right) \le
\nlength{\RLs{0}{\alpha}} \le \frac{n-\alpha}{n}\nlength{f}+1.
\]
\item
By Lemma \ref{lem dj} (\ref{dj gt 0}), we have
\begin{align*}
\nheight{\RLs{0}{\alpha}}
&=\max\{|d_0^{(\alpha)}c_0|,|d_1^{(\alpha)}c_1|,\dots,|d_{n-1}^{(\alpha)}c_{n-1}|,1\}\\
&\le \max\left\{\frac{n-\alpha}{n}\max\{|c_0|,|c_1|,\dots,|c_{n-1}|\},1\right\}
\\
&\le \frac{n-\alpha}{n}\nheight{f}+1.
\end{align*}
With Mahler's (\ref{eq mahler height}),
we get
\[
M\left(\RL{a}{\alpha} f\right) =
M\left(\RLs{a}{\alpha} f\right) \le
 \csqrt{n+1} \nheight{\RLs{0}{\alpha}} \le  \csqrt{n+1}  \max\left\{\frac{n-\alpha}{n} \nheight{f},1\right\}.
\]
\item
By Lemma \ref{lem dj} (\ref{dj gt 0}), we have
\begin{align*}
\ntwo{\RLs{0}{\alpha}}
&= \csqrt{(d_0^{(\alpha)}c_0)^2+\dots+(d_{n-1}^{(\alpha)}c_{n-1})^2+1}\\
&\le\frac{n-\alpha}{n} \csqrt{c_0^2+\dots+c_{n-1}^2}+1\\
&\le\frac{n-\alpha}{n}\ntwo{f}+1.
\end{align*}
With Landau's (\ref{eq landau}), we get
\[
M(\RL{a}{\alpha} f )=
M\left(\RLs{0}{\alpha}\right) \le \ntwo{\RLs{0}{\alpha}} \le\frac{n-\alpha}{n} ||f||_2+1.
\]
\end{enumerate}
\end{proof}


Now we are able to prove:


\begin{theorem}\label{theo lower bound}
Let \(f(x)\in\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\) be monic of degree \(n\).  Then
\begin{enumerate}
\item\label{tlb1}\(M(\RL{0}{\alpha} f)\ge 2^{-n}\frac{n-\alpha}{n}(\nlength{f}-1)+1\), for  \(\alpha<0\) and \(\alpha>2n\)
\item\label{tlb2}\(M(\RL{0}{\alpha} f)\le \frac{\ff{|n-\alpha|}{n}}{n!}\nlength{f}+1\), for  \(\alpha<0\) and \(\alpha>n\)
\item\label{tlb3}\(M(\RL{0}{\alpha} f)\ge 2^{-n}\binom{n-1}{\lceil(n-1)/2\rceil}^{-1}\frac{\alpha-n}{n}\nlength{f}\), for \(n<\alpha\le 2n\)
\end{enumerate}
\begin{proof}
\begin{enumerate}
\item[\ref{tlb1}.]
By Lemma \ref{lem dj} (\ref{dj lt 0 gt 2n}), we have 
\begin{align*}
\nlength{\RLs{0}{\alpha}f}
&=|d_0^{(\alpha)}c_0|+|d_1^{(\alpha)}c_1|+\dots+|d_{n-1}^{(\alpha)}c_{n-1}|+1\\
&\ge \frac{n-\alpha}{n}\left(|c_0|+|c_1|+\dots+|c_{n-1}|+1-1\right)+1\\
&=\frac{n-\alpha}{n}(\nlength{f}-1)+1.
\end{align*}
With Mahler's (\ref{eq mahler length}), we get
\[
M(\RL{0}{\alpha} f) =M\left(\RLs{0}{\alpha}\right) \ge 2^{-n} \nlength{\RLs{0}{\alpha}} \ge 2^{-n}\frac{n-\alpha}{n}(\nlength{f}-1)+1.
\]
\item[\ref{tlb2}.]
By Lemma \ref{lem dj} (\ref{dj lt 0 gt 2n}), we have
\begin{align*}
\nlength{\RLs{0}{\alpha}f}
&= \csqrt{(d_0^{(\alpha)}c_0)^2+\dots+(d_{n-1}^{(\alpha)}c_{n-1})^2+1}\\
&\le\frac{\ff{|n-\alpha|}{n}}{n!}\csqrt{c_0^2+\dots+c_{n-1}^2}+1\\
&\le\frac{\ff{|n-\alpha|}{n}}{n!}\ntwo{f}+1.
\end{align*}




With Landau's (\ref{eq landau}),
we get
\[
M\left(\RL{a}{\alpha} f\right) =
M\left(\RLs{a}{\alpha} f\right) \le
\nlength{\RLs{a}{\alpha} f}
\le \frac{\ff{|n-\alpha|}{n}}{n!}\ntwo{f}+1.
\]
\item[\ref{tlb3}.]
By Lemma \ref{lem dj} (\ref{dj gt n})
\begin{align*}
\nlength{\RLs{0}{\alpha}}
&=|d_0^{(\alpha)}c_0|+|d_1^{(\alpha)}c_1|+\dots+|d_{n-1}^{(\alpha)}c_{n-1}|+1\\
&\ge \frac{\alpha-n}{n}\binom{n-1}{\lceil(n-1)/2\rceil}^{-1}\left(|c_0|+|c_1|+\dots+|c_{n-1}|+1\right)\\
&=\frac{\alpha-n}{n}\binom{n-1}{\lceil(n-1)/2\rceil}^{-1}\nlength{f}
\end{align*}
With Mahler's (\ref{eq mahler length}), we get
\[
M(\RL{0}{\alpha} f) =M\left(\RLs{0}{\alpha}\right) \ge 2^{-n} \nlength{\RLs{0}{\alpha}} \ge 2^{-n}\frac{\alpha-n}{n}\binom{n-1}{(n-1)/2}^{-1}\nlength{f}.
\]
\end{enumerate} 
\end{proof}
\end{theorem}


In Figure \ref{fig p1} we present the paths of zeros and the
Mahler measures of the fractional derivatives of a degree 3 polynomial along with the bounds from Theorems \ref{theo lower bound} and \ref{theo upper bound}.
Furthermore Figures \ref{fig one}, \ref{fig double}, and \ref{fig p1} show 
the growth of $M(\RL{0}{\alpha} p)$ for \(\alpha<0\) and \(\alpha>0\).

\begin{figure}
\begin{center}
\includegraphics[width=0.45\textwidth]{paper-caputo-rieliou-mahler-bound.pdf}
\end{center}
\caption{Mahler measure of the Riemann-Liouville and Caputo fractional derivatives of the quintic \(p(x)=(x+1-2i)(x-3-2i)(x-2+5i)(x-3+3i)(x+1+5i)\)  for \(0\le\alpha\le 5\) along with the bound from
Theorem \ref{theo upper bound} and Corollary \ref{cor upper bound caputo}.}
\label{fig bound mahler}
\end{figure}


\section{Caputo Fractional Derivatives} \label{caputo}

 
Switching the order of differentiation and integration in Definition \ref{def rl} yields the Caputo fractional derivative, see \cite{diethelm2010} for example.

\begin{definition}[Caputo fractional derivative]\label{def caputo}
Let \(f\) be analytic on a convex open set \(C\) let and \(a\in C\).  Let $\alpha>0$ and $m=\lceil \alpha \rceil$ then the \(\alpha\)-th Caputo fractional derivative is
\[
\CA{a}{\alpha} f(t)
=
\RLI{a}{m-\alpha} 
\frac{d^m}{dt^m} 
f(t).
\]
\end{definition}
The Caputo fractional derivative has the advantage that the composition rule holds, while the power rule is very similar to that of the Riemann-Liouville fractional derivative: 
\[
    \CA{a}{\alpha}x^\beta= \begin{cases}
        0 & 
        \begin{array}{l}\text{if } \beta \in \{0,1,\dots, m-1\}
        \end{array}\\
        \dfrac{\Gamma(\beta+1)}{\Gamma(\beta-\alpha+1)}(x-a)^{\beta-\alpha} & 
\begin{array}{l}        
        \text{if } \beta\in\N \text{ and }\beta \ge  m\text{ or }\\[-0.5ex]
        \text{if } \beta\not\in\N \text{ and }\beta > m-1,
\end{array}
    \end{cases}
\]
where $m=\lceil \alpha \rceil$.

These derivatives differ, but obey some common general trends. Figures \ref{fig quintic paths} and \ref{fig quintic abs} compare the paths of zeros of the Riemann-Liouville and Caputo derivatives of a degree 5 polynomial.

The upper bounds from Theorem \ref{theo upper bound} easily transfer 
to the Caputo fractional derivative.  Because the coefficients of the Caputo fractional derivatives of a polynomial \(p\in\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\) are either the 
same as those of the Riemann-Liouville fractional derivative or zero, 
we have
\[
\norm{\CAs{0}{\alpha}}{k} \le \norm{\RLs{0}{\alpha}}{k},
\]
for  \(k\in\{1,2,\infty\}\).  This yields the bounds:
\begin{corollary}
\label{cor upper bound caputo}
Let \(p\in\mathbb{C}} \newcommand{\N}{\mathbb{N}[x]\) be monic of degree $n$ and let \(0<\alpha<n\).  Then 
\begin{enumerate}
\item $M(\CA{0}{\alpha} p)  \le \frac{n-\alpha}{n}\nlength{p}+1$,
\item $M(\CA{0}{\alpha} p)  \le  \csqrt{n+1} \max\left\{\frac{n-\alpha}{n} \nheight{p},1\right\}$,
\item $M(\CA{0}{\alpha} p)  \le \frac{n-\alpha}{n} \ntwo{p}+1$.
\end{enumerate}
\end{corollary}





\section{Conclusion} \label{conclusion}

The bounds we have established in Section 5 and Section 6 were sufficient for our purposes, but they are 
far from best possible. 
Figure \ref{fig bound mahler} illustrates the decline of $M(\RL{0}{\alpha} p)$ and $M(\CA{0}{\alpha} f)$, when $\alpha$ approaches \(n\), as described by Theorem \ref{theo upper bound} and Corollary \ref{cor upper bound caputo} and Theorem \ref{theo roots lim}. In a future work, it would be interesting to consider the true growth of the paths of zeros, for \(\alpha<0\) and \(\alpha>n\). Also, the maximal extent of loops of paths, after traversing the origin \(a\), is something that could be worth looking at. Moreover, as Figure \ref{fig p1} clearly shows, the paths exhibit very distinct linear asymptotes in both directions \(\alpha\to\infty\) and \(\alpha\to-\infty\). For polynomials of higher degrees, their exact directions are not yet known. 

In addition to this, as we have noted earlier, the useful Gauss-Lucas property does not hold universally for the fractional derivatives of polynomials. However, some of the dynamical properties we have observed could be investigated with insights related to those that play a key role in the integral case. In particular, Gauss himself suggested a very intriguing physical interpretation of the nontrivial critical points of a polynomial (the critical points which are not zeros) as the equilibrium points in certain force fields, generated by particles placed at the zeros of the polynomial, with masses equal to the multiplicity of the zeros and repelling with a force inversely proportional to the distance. This amazing physical application of a purely theoretical polynomial concept is exceedingly intriguing and should be investigated further. It could go a long way in explaining the profound intricacies of the paths of zeros, and their seemingly chaotic local behavior.  








2023-01-20T02:15:14.000
2301
2301.08149
en
https://arxiv.org/abs/2301.08149
We investigate the behavior of fractional derivatives of polynomials. In particular, we consider the locations and the asymptotic behaviour of their zeros and give bounds for their Mahler measure.
General Mathematics (math.GM)
Zeros of fractional derivatives of polynomials
Qwen/Qwen-72B
1. YES
2. YES