
The characteristic polynomial χ of a sequence of complex numbers z1,…,zn is defined by
χ(z):=n∏k=1(z−zk).
If z1,…,zn are now random, then one may ask about the roots z′1,…,z′n of the average characteristic polynomial
¯χ(z):=E(χ(z))=E(n∏k=1(z−zk))=n∏k=1(z−z′k).
Note that z′1,…,z′n are not necessarily real or distinct even if z1,…,zn are real and distinct. By expanding we get
¯χ(z)=zn+n∑k=1(−1)kk!zn−k∑i1≠⋯≠ikE(zi1⋯zik).
(Vieta's formulas: the coefficient of a polynomial are elementary symmetric functions of the roots, in honor of François Viète (1540-1603)). Now, if z1,…,zn are independent with common mean m then
¯χ(z)=zn+n∑k=1(−1)kk!zn−kn!(n−k)!mk=(z−m)n.
In this case z′1=⋯=z′n=m, and
1nn∑k=1δz′k=δm.
In contrast, note that by the law of large numbers,
1nn∑k=1δzka.s.⟶n→∞μ
where μ is the common law of the z′k. How about the case where z1,…,zn are dependent? Let us consider for instance the case where z1,…,zn are from the Gaussian Unitary Ensemble i.e. the eigenvalues of a n×n Gaussian Hermitian random matrix with density proportional to H↦exp(−nTr(H2)). In this case, it is well known that 1n∑nk=1δzk tend to the Wigner semi-circle law as n→∞. On the other hand, it is well known that ¯χ(z) is the n-th monic Hermite polynomial, and that 1n∑nk=1δz′k also tends to the semi-circle law as n→∞. Beyond the GUE, there are very nice answers when z1,…,zn follow a determinental point process, explored by Adrien Hardy in arXiv:1211.6564. One may ask the same for permanental point processes.
