On the zeros of random harmonic polynomials: the naive model
A complex harmonic polynomial is the sum of a complex polynomial and a conjugated complex polynomial, of degrees and respectively. Li and Wei present a formula for the expected number of zeros of a random harmonic polynomial in . In this paper we prove that if is a fixed number, this expectation is asymptotically as , and if we find a lower and upper bound of order for this expectation for sufficiently large .
A harmonic polynomial is a function of the form where and are analytic polynomials so that is harmonic over . We will refer to as a harmonic polynomial of degree and where has degree and has degree , and will use the convention . Since is not analytic, the Fundamental Theorem of Algebra does not tell us the number of zeros of . As with the number of real roots of real polynomials, there is a wide range of possibilities for the number of zeros of , which we will denote as or just . Even if the number cannot be completely determined by , some restrictions can yet be made on . Wilmhurst proved is bounded above by  if certain degenerate cases of are excluded (occuring with probability ), while the lower bound is found to be by the generalized argument principle. Both bounds are sharp, though for , the upper bound has been improved (see ) to , and it has been conjectured (see , , and ) that there is a general improvement that is linear in for each fixed .
Because of the range of values could take, one might ask for the average value of given a random harmonic polynomial, . A useful formula for determining the exact number of zeros of random real polynomials was provided by Kac  while Dunnage provided a similar formula for trigonometric polynomials . Shepp and Vanderbei go on to extend Kac’s result to determining the number of zeros of complex polynomials for a given complex domain .
The Kac formula and its generalizations have been used by Li and Wei to give an explicit formula for for a fixed and when the coefficients are independent Gaussian random variables. In the case when the variances of the coefficients are scaled by a binomial factor, they also determined the asymptotics of their formula as to be
An altered definition for the random coefficients, referred to as the “truncated model”, was studied by Lerario and Lundberg who found different asymptotics for their model, for , , only dependent on .
Although Li and Wei focused mainly on harmonic polynomials of the form
where and are independent Gaussians, with and , they also considered what I will refer to as the “naive” model of coefficients of random polynomials, namely so that the coefficients are i.i.d. They proved the following .
The expected number of zeros of
on an open domain , denoted by , is given by
where denotes the Lebesgue measure on the complex plane and
Also in the publication by the same authors,
Numerical analysis suggests that the asymptotic of above expectation for is for fixed , but a rigorous analytic asymptotic hasn’t been found.
We will prove this analytically, that is,
for fixed .
We will also explore the case when . We conjecture for some constant , because we were able to find both an upper bound and lower bound of this order.
When , there exists such that for sufficiently large
Acknowledgment. I wish to thank Erik Lundberg and Antonio Lerario for leading me in the direction of this problem and discussing it with me at length.
2. Preliminary Results
The main tool in proving Theorem 2 will be the generalization of a common theorem in Real Analysis, the General Lebegue Dominated Convergence Theorem; see  for proof.
Theorem 4 (General Lebegue Dominated Convergence Theorem).
Let be a sequence of measurable functions on that converges pointwise almost everywhere on to . Suppose there is a sequence of nonnegative measurable functions on that converges pointwise almost everywhere on to and dominates on in the sense that on for all . If
Thus in section 3 we search for the appropriate sequence in order to reach the result.
To maintain the clarity of the proof of the theorems, let us discuss a few of the inequalities used.
If is a positive real number and , , and , then , for .
That is a direct application of Cauchy-Schwarz, and the lemma follows.
Let and be non-negative functions over and let . Then and is integrable over if and are.
We first notice that the arithmetic mean of and is while the geometric mean is . The arithmetic and geometric mean inequality says that , which we relax for the purposes of this lemma. Since is a positive measurable function, integrability follows. ∎
3. Proof of Theorem 2
Beginning with the integral on the left side of (1), using polar coordinates and then a substitution of gives us
We use to be the function yielding 1 for values greater than and 0 otherwise and, for a fixed , define
We define the following functions of to simplify notation:
It follows from Lemma 5 and that
Using this, we can bound
and , by Lemma 6,
We can see when and when . When , we have , so we see that
The above bound for we will call .
Now the integrability of is important if we wish to use Theorem 4. We first note that and . From this we see
Then has an antiderivative when or ,
It follows that and clearly the limit as of this integral is .
Next we compute the point-wise limit of as , which we will call . Because when , then for , . Then
Then, to compute , we can use the following limits,
It is clear that for a fixed and , the limits as increases to infinity of , and are all as well as the limits of and . Then we can calculate
Meanwhile, if is the point-wise limit of , then
Then with a little computation we see and . Since
which proves Theorem 2.
4. Proof of Theorem 3
where , , and , again using the convention that
In the case where , then , , and . This allows us to simplify the integrand to
First we bound . By increasing all the coefficients of in from to , we can clearly see is bounded above by . Furthermore,
and for ,
A better lower bound of can be found using algebraically complicated calculus but we simply require a lower bound greater than .
Next we inspect another factor of the square of the integrand. Similar to the case with fixed , we can write as a product of and the derivative of , and similarly . As before, we use the identity that . Then
Here, can be bounded below by for all and above for for greater than some fixed . The lower bound of zero is clear. For the upper bound, we first note that it is less than for all . This can be seen by observing that and then by using the Cauchy Schwarz to see that
Then we see that
and . Thus is minimized at .
The limit of this lower bound as is and although this function seems to be decreasing as increases, we need not have a sharp lower bound and it will suffice that for large enough , we have , say, half of this limit as a lower bound. In summary, when
moreover, the upper bounds apply for all .
We will proceed by bounding the integrals over the following three intervals,
Let . Then in a similar way, for ,
Finally, for ,
We have determined the asymptotic order of growth of in two cases, namely when is fixed and when grows with . The outcome in the former case, where there are on average asymptotically zeros, is not completely unexpected. Indeed, one might expect that as increases, has a decreasing influence over the nature of when is fixed, so that may be heuristically treated as an analytic polynomial. Our result gives evidence to this conclusion since on average has zeros. On the other hand, when grows with we see a different behavior, implying that is quite distinguishable from an analytic polynomial of degree . Besides the additional factor in the total number of zeros, we also notice a disparity in the shape of the so-called “first intensity”, that is, in the integrand appearing in Theorem 1. In the case that is fixed the first intensity divided by converges to everywhere except along the unit circle. Moreover, after the change of variable, we find that the area to the left of under the curve is equal to the area to the right of , where corresponds to through the change of variables. This tells us that as increases, not only do the zeros concentrate on the unit circle, but they seem to do so in a symmetric manner (resembling the known symmetry of the first intensity of an analytic polynomial). In contrast, the first intensity when accumulates to the unit circle in an asymmetric manner when divided by , which can be seen when dividing the integrand by instead of and taking a limit. After the change of variables it is clear that this limit is for while it is for .
There are natural questions extending from this discussion. Though we know the order of growth of for to be , it is yet an open problem to prove the limit of the quotient exists and, if it does, to obtain the explicit value of the constant. One may also notice that we have only explored two cases for the behavior of . An obvious extension would be to investigate the case where , as was done in  and .
-  J. E. A. Dunnage, The number of real zeros of a random trigonometric polynomial, Proceedings of the London Mathematical Society s3-16 (1966), 53–84.
-  M. Kac, On the average number of real roots of a random algebraic equation, Bulletin of the American Mathematical Society 49 (1943), no. 4, 314–320.
-  D. Khavinson, S-Y. Lee, and A. Saez, Zeros of harmonic polynomials, critical lemniscates and caustics, preprint at arXiv:1508.04439.
-  D. Khavinson and G. Swiatek, On the maximal number of zeros of certain harmonic polynomials, Proceedings of the American Mathematical Society 131 (2003), 409–414.
-  S-Y. Lee, A. Lerario, and E. Lundberg, Remarks on Wilmshurst’s theorem, Indiana University Mathematics Journal 64 (2015), no. 4, 1153–1167.
-  A. Lerario and E. Lundberg, On the zeros of random harmonic polynomials: the truncated model, J. Math. Analysis Appl. 438 (2016), 1041–1054.
-  W. V. Li and A. Wei, On the expected number of zeros of random harmonic polynomials, Proceedings of the American Mathematical Society 137 (2009), 195–204.
-  H. L. Royden and P. M. Fitzpatrick, Real analysis, fourth ed., Pearson Education, Inc., Boston, MA, 2010.
-  Terry Sheil-Small, Complex polynomials, Cambridge University Press, 2002.
-  L. Shepp and R. Vanderbei, The complex zeros of random polynomials, Transactions of the American Mathematical Society 347 (1995), no. 11, 4365–4384.
-  A. S. Wilmshurst, The valence of harmonic polynomials, Proceedings of the American Mathematical Society 126 (1998), no. 7, 2077–2081.