Characteristic polynomials of certain random symmetric matrices and the complexity of random Morse functions - MathOverflow most recent 30 from http://mathoverflow.net 2013-05-25T19:28:52Z http://mathoverflow.net/feeds/question/84687 http://www.creativecommons.org/licenses/by-nc/2.5/rdf http://mathoverflow.net/questions/84687/characteristic-polynomials-of-certain-random-symmetric-matrices-and-the-complexi Characteristic polynomials of certain random symmetric matrices and the complexity of random Morse functions Liviu Nicolaescu 2012-01-01T16:51:00Z 2012-01-09T10:16:41Z <p>Investigations concerning random Morse functions led me to the following problem. Consider the classical GOE of $m\times m$ real symmetric matrices $A$ with independent Gaussian entries with zero means and variances</p> <p>$$\boldsymbol{E}(a_{ii}^2)=2 \boldsymbol{E}(a_{ij}^2)= 2$$</p> <p>for all $i \neq j$. Consider the function</p> <p>$$F_m(x, y) = \boldsymbol{E}_{GOE}\bigl( |\det(y+A)| e^{ -x(tr A)^2 } \bigr),$$</p> <p>$x,y$ real, $x>0$. What can one say about the behavior of $F_m(x,y)$ as $m\rightarrow \infty$.</p> <p>Equivalently we can consider the Gaussian ensemble $\mathcal{S}(m,x)$ of symmetric $m\times m$ real matrices with probability density</p> <p>$$dP(A)=\frac{1}{Z_{m,x}} e^{-\frac{1}{2}tr(A^2)-x(tr A)^2} \prod_{i\leq j} da_{ij},$$</p> <p>$x>0$, and then ask for the bevavior as $m\rightarrow \infty$ of the expectation</p> <p>$$\boldsymbol{E}_{\mathcal{S}(m,x)}\left( |\det(A+y)|\right).$$</p> <p>Observe that </p> <p>$$GOE= \mathcal{S}(m,x)_{x=0}.$$</p> <p>The normalizing constant $Z_{m,x}$ can be explicitly computed for any $x$ and thus</p> <p>$$F_m(x,y)= \frac{Z_{m,x}}{Z_{m,0}} \boldsymbol{E}_{\mathcal{S}(m,x)}\left( |\det(A+y)|\right).$$</p> <p>In the geometric problem I am interested $x=\frac{1}{8}$. In this case the ensemble $\mathcal{S}_m:=\mathcal{S}(m, \frac{1}{4})$ can be described as the ensemble of real, symmetric $m\times m$ matrices whose entries are mean zero Gaussian variables satisfying the covariance equalities</p> <p>$$\boldsymbol{E}\left( a_{ij} a_{k\ell}\right)=-\frac{2}{2+m}\delta_{ij}\delta_{k\ell} +\left( \delta_{ik}\delta_{j\ell}+ \delta_{i\ell}\delta_{jk}\right).$$</p> <p>Note that as $m\rightarrow \infty$ this ensemble resembles more and more the classical GOE which satisfies the covariance equalities</p> <p>$$\boldsymbol{E}\left( a_{ij} a_{k\ell}\right)= \left(\delta_{ik}\delta_{j\ell}+ \delta_{i\ell}\delta_{jk}\right).$$</p> <p>Finally, I want to explain how is this related to Morse theory. To put things in perspective observe that if $A$ is a symmetric $m\times m$ matrix, then its spectrum can be identified with the set of critical values of the restriction to the unit sphere in $\mathbb{R}^m$ of the quadratic polynomial</p> <p>$$\mathbb{R}^m\ni x\mapsto q_A(X)=(Ax,x).$$</p> <p>To a Morse function $f$ on a compact smooth manifold $M$ of dimension $m$ we can associate two measures.</p> <p>(a) A measure $K_f$ on $M$ defined as the sum of Dirac delta's concentrated at the critical points of $f$</p> <p>$$K_f=\sum_{df(p)=0}\delta_p.$$</p> <p>(b) A measure $\Delta_f$ on $\mathbb{R}$ supported on the set of critical values of $f$ and defined as the pushforward of $K_f$ via $f$, </p> <p>$$\Delta_f:=f_*(K_f).$$</p> <p>In other words, $\Delta_f$ counts the critical values with multiplicity. Note that when $f$ is the restriction to the unit sphere of the quadratic form $q_A$ then $\Delta_f$ coincides with the spectral measure of $A$. </p> <p>Fix a Riemann metric $g$ on $M$ and an orthonormal $(\Psi_k)_{k\geq 0}$ basis of $L^2(M)$ consisting of eigenfunctions of the Laplacian</p> <p>$$\Delta \Psi_k=\lambda_k \Psi_k.$$</p> <p>Fix i.i.d. standard Gaussian random variables $(x_k)_{k\geq 0}$ and for every $L >0$ define the random function</p> <p>$$f_L=\sum_{\lambda_k\leq L^2}x_k\Psi_k.$$</p> <p>The function $f_L$ is roughly speaking a random polynomial of large degree. Equivalently one should think of $f_L$ as a random element in the space $U_L$ spanned by the eigenfunctions corresponding to eigenvalues $\leq L^2$ and equipped with the standard Gaussian measure. The large $L$ behavior of $\dim U_L$ is governed by Weyl's asymptotic formula</p> <p>$$\dim U_L \sim const. L^m.$$</p> <p>To $f_L$ we associate two random measures</p> <p>$$K_{f_L},\;\; \Delta_{f_L}$$</p> <p>that have normalized expectations</p> <p>$$K_L:=\frac{1}{\dim U_L} \boldsymbol{E}( K_{f_L} ),$$</p> <p>$$\Delta_L:=\frac{1}{\dim U_L} \boldsymbol{E}( \Delta_{f_L} ).$$</p> <p>Above, $K_L$ is a measure on $M$ and $\Delta_L$ is a measure on $\mathbb{R}$. I can show that as $L\to\infty$ the measure $K_L$ converges weakly to $C_m dV_g$, where $dV_g$ denotes the volume measure determined by the metric $g$, and $C_m$ is a certain explicit constant that depends only on $m$ but not on $(M,g)$. Thus, the critical points of a random $f_L$, $L\gg 0$, is uniformly distributed on average.</p> <p>As $L\to \infty$ the measure a suitable rescaled version of $\Delta_L$ converges to a measure $d\mu_m(y)$ on $\mathbb{R}$ that is absolutely continuous with respect to the Lebesgue measure. More precisely</p> <p>$$d\mu_m(y)=\rho_m(y) dy= Const_m \times \boldsymbol{E}_{\mathcal{S}(m,1/8)}\left( \;|\det(A-s_my )|\;\right) e^{-\frac{y^2}{2 }} dy,$$</p> <p>$$s_m=\sqrt{\frac{m+4}{m+2}}$$</p> <p><strong>Remark.</strong> The measure $d\mu_m(y)$ can also be given a description as a conditional expectation. To explain this I need to introduce another Gaussian ensemble of symmetric $m\times m$ matrices.</p> <p>To describe it observe that to any such matrix $A$ we can associate a quadratic form $q_A$ on $\mathbb{R}^m$,</p> <p>$$q_A(x)=(Ax, x).$$</p> <p>We have a unique, centered Gaussian probability measure on the space of symmetric $m\times m$ matrices with variance</p> <p>$$V(A)=\int_{\mathbb{R^m}} q_A(x)^2 \frac{e^{-\frac{|x|^2}{2}}}{(2\pi)^{\frac{m}{2}}} dx.$$</p> <p>Denote by $\mathcal{U}_m$ this Gaussian ensemble of symmetric matrices. (I use the symbol $\mathcal{U}_m$ because this ensemble has a <a href="http://www.nd.edu/~lnicolae/CritSetStat.pdf" rel="nofollow">remarkable universality property.</a>)</p> <p>Now fix a standard (scalar) Gaussian r.v. $Y$ such that the pair $(A,Y)$ is a Gaussian vector satisfying the correlation equalities</p> <p>$$\boldsymbol{E}(a_{ij} Y)=s_m\delta_{ij}$$.</p> <p>Then for any Borel subset of $\mathbb{R}$ we have</p> <p>$$\mu_m(B)=\boldsymbol{E}_{\mathcal{U}_m}\Bigl( |\det A|\;\Bigl|\; Y\in B\Bigr).$$</p>