Skip to main content
9 of 13
deleted 9 characters in body
Olivier Bégassat
  • 2.8k
  • 4
  • 25
  • 30

Polynomial roots and convexity

A couple of years back, I came up with the following question, to which I have no answer to this day. I asked a few people about this, most of my teachers and some friends, but noone had ever heard of the question before, and noone knew the answer.

I hope this is an original question, but seeing how natural it is, I doubt this is the first time someone has asked it.

First, some motivation. Take $P$ any non zero complex polynomial. It is an easy and classical exercise to show that the roots of its derivative $P'$ lie in the convex hull of its own roots (I know this as the Gauss-Lucas property). To show this, you simply write $P=a\cdot\prod_{i=1}^{r}(X-\alpha_i)^{m_i}$ where the $\alpha_i~(i=1,\dots,r)$ are the different roots of $P$, and $m_i$ the corresponding multiplicities, and evaluate $\frac{P'}{P}=\sum_i \frac{m_i}{X-\alpha_i}$ on a root $\beta$ of $P'$ which is not also a root of $P$. You'll end up with an expression of $\beta$ as a convex combination of $\alpha_1,\dots,\alpha_r$. It is worth mentioning that all the convex coefficients are $>0$, so the new root cannot lie on the edge of the convex hull of $P$'s roots.

Now fix $P$ a certain non zero complex polynomial, and consider $\Pi$ it's primitive (antiderivative) that vanishes at $0:~\Pi(0)=0$ and $\Pi'=P$. For each complex $\omega$, write $\Pi_{\omega}=\Pi-\omega$ so that you get all the primitives of $P$. Also, define for any polynomial $Q$, $\mathrm{Conv}(Q)$, to be the convex hull of $Q$'s roots.

$\mathrm{MAIN~QUESTION}$: describe $\mathrm{Hull}(P)=\bigcap_{\omega\in\mathbb{C}}\mathrm{Conv}(\Pi_{\omega})$.

By the above quoted property, $\mathrm{Hull}(P)$ is a convex compact subset of the complex plane that contains $\mathrm{Conv}(P)$, but I strongly suspect that it is in general larger.

Here are some easy observations:

$1)$ replacing $P$ (resp. $\Pi$) by $\lambda P$ (resp. $\lambda \Pi$) will not change the result, and considering $P(aX+b)$ will change $\mathrm{Hull}(P)$ accordingly. Hence we can suppose both $P$ and $\Pi$ be unitary. The fact that $\Pi$ is no longer a primitive of $P$ is of no consequence.

$2)$ the intersection defining $\mathrm{Hull}(P)$ can be taken for $\omega$ ranging in a compact subset of $\mathbb{C}$: as $|\omega|\rightarrow\infty$, the roots of $\Pi_{\omega}$ will tend to become close to the $(\deg (P)+1)$-th roots of $\omega$, so for large enough $\omega$, their convex hull will always contain , say, $\mathrm{Conv}(\Pi)$.

$3)$ $\mathrm{Hull}(P)$ can be explicitely calculated in the following cases: $P=X^n$, $P$ of degree $1$ or $2$. There are only 2 kinds of degree $2$ polynomials: 2 simple roots or a double root. Using $z\rightarrow az+b$, one only has to consider $P=X^2$ and $P=X(X-1)$. The first one yields {$0$}, which equals $\mathrm{Conv}(X^2)$, the second one gives $[0,1]=\mathrm{Conv}(X(X-1))$.

Also if $\Pi$ is a real polynomial of odd degree $n+1$, that has all its roots real and simple, say $\lambda_1<\mu_1<\lambda_2<\dots<\mu_n<\lambda_{n+1}$, where I have also placed $P$'s roots $\mu_1,\dots,\mu_n$, and if you further assume that $\Pi(\mu_{2j})\leq\Pi(\mu_n)\leq\Pi(\mu_1)\leq\Pi(\mu_{2j+1})$ for all suitable $j$ (a condition that is best understood with a picture), then $\mathrm{Hull}(P)=\mathrm{Conv}(P)=[\mu_1,\mu_n]$: just vary $\omega$ between $[\Pi(\mu_n),\Pi(\mu_1)]$, the resulting polynomial $\Pi_{\omega}$ is always split over the real numbers and you get

$[\mu_1,\mu_n]=\mathrm{Conv}(P)\subset\mathrm{Hull}(P)\subset \mathrm{Conv}(\Pi_{\Pi(\mu_1)})\cap \mathrm{Conv}(\Pi_{\Pi(\mu_n)})=$ $[\mu_1,\dots]\cap [\dots,\mu_n]=[\mu_1,\mu_n]$

$4)$ The equation $\Pi_{\omega}(z)=\Pi(z)-\omega=0$ defines a Riemann surface, but I don't see how that could be of any use.

Computing $\mathrm{Hull}(P)$ for the next most simple polynomial $P=X^3-1$ has proven a challenge, and I can only conjecture what it might be.

Computing $\mathrm{Hull}(X^3-1)$ requires factorizing degree 4 polynomials, so one naturally tries to look for good values of $\omega$, the $\omega$ that allow for easy factorisation of $\Pi_{\omega}=X^4-4X-\omega$. For instance the $\omega$ that produce a double root. All that remains to be done afterwards is to factor a polynomial of degree $2$. The problem is symmetric, and you can focus on the case where 1 is the double root (i.e. $\omega=-3$). Plugging the result in the intersection, and rotating twice, you obtain the following superset of $\mathrm{Hull}(X^3-1)$: a hexagon that is the intersection of 3 similar isocele triangles with their main vertex located on the three third roots of unity $1,j,j^2$

$\mathrm{QUESTION}$: is this hexagon equal to $\mathrm{Hull}(X^3-1)$?

Here's why I think this might be.

Consider the question of how the convex hulls of the roots of $\Pi_{\omega}$ vary as $\omega$ varies. When $\omega_0$ is such that all roots of $\Pi_{\omega_0}$ are simple, then the inverse function theorem shows that the roots of $\Pi_{\omega}$ with $\omega$ in a small neighborhood of $\omega_0$ vary holomorphically $\sim$ linearly in $\omega-\omega_0$: $z(\omega)-z(\omega_0)\sim \omega-\omega_0$. If however $\omega_0$ is such that $\Pi_{\omega_0}$ has a multiple root $z_0$ of multiplicity $m>1$, then a small variation of $\omega$ about $\omega_0$ will split the multiple root $z_0$ into $m$ distinct roots of $\Pi_{\omega}$ that will spread out roughly as $z_0+c(\omega-\omega_0)^{\frac{1}{m}}$, where $c$ is some non zero coefficient. This means that for small variations, these roots will move at much higher velocities than the simple roots, and they will do the major contribution to the variation of $\mathrm{Conv}(\Pi_{\omega})$, also, they spread evenly out, and (at least if the multiplicity is greater or equal to $3$) they will tend to increase the convex hull around $z_0$. Thus it seems not too unreasonable to conjecture that the convex hull $\mathrm{Conv}(\Pi_{\omega})$ has what one can only describe as critical points at the $\omega_0$ that produce roots with multiplicities. I'm fairly certain there is a sort of calculus on convex sets that would allow one to make the above statement precise, but I don't know it.

Back to $X^3-1$: explicit calculations suggest that up to second order the double root $1$ of $X^4-4X+3-h$ for $|h|<<1$ splits in half nicely (here $\omega=-3+h$) and the convex hull will continue to contain the aforementioned hexagon.

$\mathrm{QUESTION\Conjecture}$: is it true that $\mathrm{Hull}(P)=\bigcap_{\omega\in\mathrm{MR}}\mathrm{Conv}(\Pi_{\omega})$, where $\mathrm{MR}$ is the set of all $\omega_0$ such that $\Pi_{\omega_0}$ has a multiple root, i.e. the set of all $\Pi(\alpha_i)$ where the $\alpha_i$ are the roots of $P$?

All previous examples of calculations agree with this, and I have tried as best I can to justify this guess heuristically.

Are you aware of a solution? Is this a classical problem? Is anybody brave enough to make a computer program that would compute some intersections of convex hulls obtained from the roots to see if my conjecture is any good?

Olivier Bégassat
  • 2.8k
  • 4
  • 25
  • 30