{"id":969,"date":"2010-11-02T00:27:40","date_gmt":"2010-11-01T22:27:40","guid":{"rendered":"http:\/\/djalil.chafai.net\/blog\/?p=969"},"modified":"2017-04-21T15:25:19","modified_gmt":"2017-04-21T13:25:19","slug":"aspects-of-the-complex-ginibre-ensemble","status":"publish","type":"post","link":"https:\/\/djalil.chafai.net\/blog\/2010\/11\/02\/aspects-of-the-complex-ginibre-ensemble\/","title":{"rendered":"Aspects of the Complex Ginibre Ensemble"},"content":{"rendered":"<figure id=\"attachment_9489\" aria-describedby=\"caption-attachment-9489\" style=\"width: 289px\" class=\"wp-caption aligncenter\"><a href=\"http:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/ginibre.jpg\"><img loading=\"lazy\" class=\"wp-image-9489 size-medium\" src=\"http:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/ginibre-289x300.jpg\" alt=\"Ginibre Ensemble 1000 x 1000\" width=\"289\" height=\"300\" srcset=\"https:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/ginibre-289x300.jpg 289w, https:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/ginibre-144x150.jpg 144w, https:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/ginibre.jpg 553w\" sizes=\"(max-width: 289px) 100vw, 289px\" \/><\/a><figcaption id=\"caption-attachment-9489\" class=\"wp-caption-text\">Ginibre Ensemble of size 1000 x 1000. GNU-R simulation courtesy of A. Hardy.<\/figcaption><\/figure>\n<p style=\"text-align: justify;\">This post is devoted to the Complex Ginibre Ensemble, the subject of an expository talk that I gave few months ago. Let \\( {(G_{i,j})_{i,j\\geq1}} \\) be an infinite table of i.i.d. random variables on \\( {\\mathbb{C}\\equiv\\mathbb{R}^2} \\) with \\( {G_{11}\\sim\\mathcal{N}(0,\\frac{1}{2}I_2)} \\). The Lebesgue density of the \\( {n\\times n} \\) random matrix \\( {\\mathbf{G}=(G_{i,j})_{1\\leq i,j\\leq n}} \\) in \\( {\\mathcal{M}_n(\\mathbb{C})\\equiv\\mathbb{C}^{n\\times n}} \\) is <a id=\"eqginden\"><\/a><\/p>\n<p style=\"text-align: center;\">\\[ \\mathbf{A}\\in\\mathcal{M}_n(\\mathbb{C}) \\mapsto \\pi^{-n^2}e^{-\\sum_{i,j=1}^n|\\mathbf{A}_{ij}|^2} = \\pi^{-n^2}e^{-\\mathrm{Tr}(\\mathbf{A}\\mathbf{A}^*)} \\ \\ \\ \\ \\ (1) \\]<\/p>\n<p style=\"text-align: justify;\">where \\( {\\mathbf{A}^*} \\) the conjugate-transpose of \\( {\\mathbf{A}} \\). We say that \\( {\\mathbf{G}} \\) belong to the <em>Complex Ginibre Ensemble<\/em>. This law is unitary invariant, in the sense that if \\( {\\mathbf{U}} \\) and \\( {\\mathbf{V}} \\) are \\( {n\\times n} \\) unitary matrices then \\( {\\mathbf{U}\\mathbf{G}\\mathbf{V}} \\) and \\( {\\mathbf{G}} \\) are equally distributed. The entry-wise and matrix-wise real and imaginary parts of \\( {\\mathbf{G}} \\) are independent and belong to the Gaussian Unitary Ensemble (GUE) of Gaussian random Hermitian matrices (the converse is also true).<\/p>\n<p style=\"text-align: justify;\">The <em>eigenvalues<\/em> \\( {\\lambda_1(\\mathbf{A}),\\ldots,\\lambda_n(\\mathbf{A})} \\) of a matrix \\( {\\mathbf{A}\\in\\mathcal{M}_n(\\mathbb{C})} \\) are the roots in \\( {\\mathbb{C}} \\) of its characteristic polynomial \\( {P_\\mathbf{A}(z)=\\det(\\mathbf{A}-z\\mathbf{I})} \\), labeled so that \\( {|\\lambda_1(\\mathbf{A})|\\geq\\cdots\\geq|\\lambda_n(\\mathbf{A})|} \\) with growing phases.<\/p>\n<blockquote style=\"background: white; border: solid thick #e4e5e7; text-align: justify; padding-left: 1em;\"><p><b>Lemma 1 (Diagonalizability)<\/b> <em><a id=\"lediag\"><\/a> For every \\( {n\\geq1} \\), the set of elements of \\( {\\mathcal{M}_n(\\mathbb{C})} \\) with multiple eigenvalues has zero Lebesgue measure in \\( {\\mathbb{C}^{n\\times n}} \\). In particular, the set of nondiagonalizable elements of \\( {\\mathcal{M}_n(\\mathbb{C})} \\) has zero Lebesgue measure in \\( {\\mathbb{C}^{n\\times n}} \\).<\/em><\/p><\/blockquote>\n<p style=\"text-align: justify;\"><em>Proof:<\/em> If \\( {\\mathbf{A}\\in\\mathcal{M}_n(\\mathbb{C})} \\) has characteristic polynomial<\/p>\n<p style=\"text-align: center;\">\\[ P_\\mathbf{A}(z)=z^n+a_{n-1}z^{n-1}+\\cdots+a_0, \\]<\/p>\n<p style=\"text-align: justify;\">then \\( {a_0,\\ldots,a_{n-1}} \\) are polynomials functions of \\( {\\mathbf{A}} \\). The <em>resultant<\/em> \\( {R(P_\\mathbf{A},P'_\\mathbf{A})} \\) of \\( {P_\\mathbf{A},P'_\\mathbf{A}} \\), called the <em>discriminant<\/em> of \\( {P_\\mathbf{A}} \\), is the determinant of the \\( {(2n-1)\\times(2n-1)} \\) Sylvester matrix of \\( {P_\\mathbf{A},P'_\\mathbf{A}} \\). It is a polynomial in \\( {a_0,\\ldots,a_{n-1}} \\). We have also the Vandermonde formula<\/p>\n<p style=\"text-align: center;\">\\[ |R(P_\\mathbf{A},P'_\\mathbf{A})|=\\prod_{i{&lt;}j}|\\lambda_i(\\mathbf{A})-\\lambda_j(\\mathbf{A})|^2. \\]<\/p>\n<p style=\"text-align: justify;\">Consequently, \\( {\\mathbf{A}} \\) has all eigenvalues distinct if and only if \\( {\\mathbf{A}} \\) lies outside the polynomial hyper-surface \\( {\\{\\mathbf{A}\\in\\mathbb{C}^{n\\times n}:R(P_\\mathbf{A},P'_\\mathbf{A})=0\\}} \\). \u2610<\/p>\n<p style=\"text-align: justify;\">Since the law of \\( {\\mathbf{G}} \\) is absolutely continuous, from lemma <a href=\"#lediag\">1<\/a>, we get that a.s. \\( {\\mathbf{G}\\mathbf{G}^*\\neq \\mathbf{G}^*\\mathbf{G}} \\) (nonnormality) but \\( {\\mathbf{G}} \\) is diagonalizable with distinct eigenvalues. Note that if \\( {\\mathbf{G}=\\mathbf{U}\\mathbf{T}\\mathbf{U}^*} \\) is the Schur unitary decomposition where \\( {\\mathbf{U}} \\) is unitary and \\( {\\mathbf{T}} \\) is upper triangular, then \\( {\\mathbf{T}=\\mathbf{D}+\\mathbf{N}} \\) with \\( {\\mathbf{D}=\\mathrm{diag}(\\lambda_1(\\mathbf{G}),\\ldots,\\lambda_n(\\mathbf{G}))} \\) and<\/p>\n<p style=\"text-align: center;\">\\[ \\mathrm{Tr}(\\mathbf{G}\\mathbf{G}^*) =\\mathrm{Tr}(\\mathbf{D}\\mathbf{D}^*) +\\mathrm{Tr}(\\mathbf{N}\\mathbf{N}^*). \\]<\/p>\n<p style=\"text-align: justify;\">Following <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=0173726\">Ginibre<\/a>, one may compute the joint density of the eigenvalues by integrating <a href=\"#eqginden\">(1)<\/a> over the non spectral variables \\( {\\mathbf{U}} \\) and \\( {\\mathbf{N}} \\). The result is stated in Theorem <a href=\"#thcomplex-ginibre-density\">2<\/a> below. The law of \\( {\\mathbf{G}} \\) is invariant by the multiplication of the entries with a common phase, and thus the law spectrum of \\( {\\mathbf{G}} \\) is rotationally invariant in \\( {\\mathbb{C}^n} \\). In the sequel we set<\/p>\n<p style=\"text-align: center;\">\\[ \\Delta_n:=\\{(z_1,\\ldots,z_n)\\in\\mathbb{C}^n:|z_1|\\geq\\cdots\\geq|z_n|\\}. \\]<\/p>\n<blockquote style=\"background: white; border: solid thick #e4e5e7; text-align: justify; padding-left: 1em;\"><p><b>Theorem 2 (Spectrum law)<\/b> <em><a id=\"thcomplex-ginibre-density\"><\/a> \\( {(\\lambda_1(\\mathbf{G}),\\ldots,\\lambda_n(\\mathbf{G}))} \\) has density \\( {n!\\varphi_n\\mathbf{1}_{\\Delta_n}} \\) where<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\varphi_n(z_1,\\ldots,z_n)=\\frac{\\pi^{-n}}{\\prod_{k=1}^nk!} \\exp\\left(-\\sum_{k=1}^n|z_k|^2\\right)\\prod_{1\\leq i{&lt;}j\\leq n}|z_i-z_j|^2. \\]<\/em><\/p>\n<p><em>In particular, for every symmetric Borel function \\( {F:\\mathbb{C}^n\\rightarrow\\mathbb{R}} \\),<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\mathbb{E}[F(\\lambda_1(\\mathbf{G}),\\ldots,\\lambda_n(\\mathbf{G}))] =\\int_{\\mathbb{C}^n}\\!F(z_1,\\ldots,z_n)\\varphi_n(z_1,\\ldots,z_n)\\,dz_1\\cdots dz_n. \\]<\/em><\/p>\n<\/blockquote>\n<p style=\"text-align: justify;\">We will use Theorem <a href=\"#thcomplex-ginibre-density\">2<\/a> with symmetric functions of the form<\/p>\n<p style=\"text-align: center;\">\\[ F(z_1,\\ldots,z_n) =\\sum_{i_1,\\ldots,i_k \\text{ distinct}}f(z_{i_1})\\cdots f(z_{i_k}). \\]<\/p>\n<p style=\"text-align: justify;\">The Vandermonde determinant comes from the Jacobian of the diagonalization, and can be interpreted as an electrostatic repulsion. The spectrum is a Gaussian determinantal process. One may also take a look at the fourth chapter of the <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2552864\">book by Hough, Krishnapur, Peres, and Virag<\/a> for a generalization to zeros of Gaussian analytic functions. Theorem <a href=\"#thcomplex-ginibre-density\">2<\/a> is reported in chapter 15 of <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2129906\">Mehta's book<\/a>. <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=0173726\">Ginibre<\/a> considered additionally the case where \\( {\\mathbb{C}} \\) is replaced by \\( {\\mathbb{R}} \\) or by the quaternions. These two cases are less studied than the complex case, due to their peculiarities. For instance, for the real Gaussian case, and following <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=1437734\">Edelman<\/a> (Corollary 7.2), the probability that the \\( {n\\times n} \\) matrix has all its eigenvalues real is \\( {2^{-n(n-1)\/4}} \\), see also <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2363393\">Akemann and Kanzieper<\/a>. The whole spectrum does not have a density in \\( {\\mathbb{C}^n} \\) in this case!<\/p>\n<blockquote style=\"background: white; border: solid thick #e4e5e7; text-align: justify; padding-left: 1em;\"><p><b>Theorem 3 (\\( {k} \\)-points correlations)<\/b> <em><a id=\"thcomplex-ginibre-k-points\"><\/a> Let \\( {z\\in\\mathbb{C}\\mapsto\\gamma(z)=\\pi^{-1}e^{-|z|^2}} \\) be the density of the standard Gaussian law \\( {\\mathcal{N}(0,\\frac{1}{2}I_2)} \\) on \\( {\\mathbb{C}} \\). Then for every \\( {1\\leq k\\leq n} \\), the ``\\( {k} \\)-point correlation''<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\varphi_{n,k}(z_1,\\ldots,z_k) := \\int_{\\mathbb{C}^{n-k}}\\!\\varphi_n(z_1,\\ldots,z_k)\\,dz_{k+1}\\cdots dz_n \\]<\/em><\/p>\n<p><em>satisfies to<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\varphi_{n,k}(z_1,\\ldots,z_k) = \\frac{(n-k)!}{n!}\\gamma(z_1)\\cdots\\gamma(z_k) \\det\\left[K(z_i,z_j)\\right]_{1\\leq i,j\\leq k} \\]<\/em><\/p>\n<p><em>where<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ K(z_i,z_j) :=\\sum_{\\ell=0}^{n-1}\\frac{(z_iz_j^*)^\\ell}{\\ell!} =\\sum_{\\ell=0}^{n-1}H_\\ell(z_i)H_\\ell(z_j)^* \\quad\\text{with}\\quad H_\\ell(z):=\\frac{1}{\\sqrt{\\ell!}}z^\\ell. \\]<\/em><\/p>\n<p><em>In particular, by taking \\( {k=n} \\) we get<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\varphi_{n,n}(z_1,\\ldots,z_n) =\\varphi_n(z_1,\\ldots,z_n) =\\frac{1}{n!}\\gamma(z_1)\\cdots\\gamma(z_n)\\det\\left[K(z_i,z_j)\\right]_{1\\leq i,j\\leq n}. \\]<\/em><\/p>\n<\/blockquote>\n<p style=\"text-align: justify;\"><em>Proof:<\/em> Calculations made by <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2129906\">Mehta<\/a> (chapter 15 page 271 equation 15.1.29) using<\/p>\n<p style=\"text-align: center;\">\\[ \\prod_{1\\leq i{&lt;}j\\leq n}|z_i-z_j|^2 =\\prod_{1\\leq i{&lt;}j\\leq n}(z_i-z_j)\\prod_{1\\leq i{&lt;}j\\leq n}(z_i-z_j)^* \\]<\/p>\n<p style=\"text-align: justify;\">and<\/p>\n<p style=\"text-align: center;\">\\[ \\det\\left[z_j^{i-1}\\right]_{1\\leq i,j\\leq k}\\det\\left[(z_j^*)^{i-1}\\right]_{1\\leq i,j\\leq k} =\\frac{\\prod_{j=1}^kj!}{n!}\\det\\left[K(z_i,z_j)\\right]_{1\\leq i,j\\leq k}. \\]<\/p>\n<p style=\"text-align: justify;\">\u2610<\/p>\n<p style=\"text-align: justify;\">Recall that the <em>empirical spectral distribution<\/em> of an \\( {n\\times n} \\) matrix \\( {A} \\) is given by<\/p>\n<p style=\"text-align: center;\">\\[ \\mu_\\mathbf{A}:=\\frac{1}{n}\\sum_{k=1}^n\\delta_{\\lambda_k(\\mathbf{A})}. \\]<\/p>\n<blockquote style=\"background: white; border: solid thick #e4e5e7; text-align: justify; padding-left: 1em;\"><p><b>Theorem 4 (Mean Circular Law)<\/b> <em><a id=\"thcomplex-ginibre-circle-weak\"><\/a> For every continuous bounded function \\( {f:\\mathbb{C}\\rightarrow\\mathbb{R}} \\),<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\lim_{n\\rightarrow\\infty}\\mathbb{E}\\left[\\int\\!f\\,d\\mu_{\\frac{1}{\\sqrt{n}}\\mathbf{G}}\\right] =\\pi^{-1}\\int_{|z|\\leq 1}\\!f(z)\\,dxdy. \\]<\/em><\/p>\n<\/blockquote>\n<p style=\"text-align: justify;\"><em>Proof:<\/em> From Theorem <a href=\"#thcomplex-ginibre-k-points\">3<\/a>, with \\( {k=1} \\), we get that the density of \\( {\\mathbb{E}\\mu_{\\mathbf{G}}} \\) is<\/p>\n<p style=\"text-align: center;\">\\[ \\varphi_{n,1}: z\\mapsto \\gamma(z)\\left(\\frac{1}{n}\\sum_{\\ell=0}^{n-1}|H_\\ell|^2(z)\\right) =\\frac{1}{n\\pi}e^{-|z|^2}\\sum_{\\ell=0}^{n-1}\\frac{|z|^{2\\ell}}{\\ell!}. \\]<\/p>\n<p style=\"text-align: justify;\">Following <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2129906\">Mehta<\/a> (chapter 15 page 272), by elementary calculus, for every compact \\( {C\\subset\\mathbb{C}} \\),<\/p>\n<p style=\"text-align: center;\">\\[ \\lim_{n\\rightarrow\\infty}\\sup_{z\\in C} \\left|n\\varphi_{n,1}(\\sqrt{n}z)-\\pi^{-1}\\mathbf{1}_{[0,1]}(|z|)\\right|=0. \\]<\/p>\n<p style=\"text-align: justify;\">The \\( {n} \\) in front of \\( {\\varphi_{n,1}} \\) is due to the fact that we are on the complex plane \\( {\\mathbb{C}=\\mathbb{R}^2} \\) and thus \\( {d\\sqrt{n}xd\\sqrt{n}y=ndxdy} \\). Here is the start of the elementary calculus: for \\( {r^2&lt;n} \\),<\/p>\n<p style=\"text-align: center;\">\\[ e^{r^2}-\\sum_{\\ell=0}^{n-1}\\frac{r^{2\\ell}}{\\ell!} =\\sum_{\\ell=n}^\\infty\\frac{r^{2\\ell}}{\\ell!} \\leq \\frac{r^{2n}}{n!}\\sum_{\\ell=0}^\\infty\\frac{r^{2\\ell}}{(n+1)^\\ell} =\\frac{r^{2n}}{n!}\\frac{n+1}{n+1-r^2} \\]<\/p>\n<p style=\"text-align: justify;\">while for \\( {r^2&gt;n} \\),<\/p>\n<p style=\"text-align: center;\">\\[ \\sum_{\\ell=0}^{n-1}\\frac{r^{2\\ell}}{\\ell!} \\leq\\frac{r^{2(n-1)}}{(n-1)!}\\sum_{\\ell=0}^{n-1}\\left(\\frac{n-1}{r^2}\\right)^\\ell \\leq \\frac{r^{2(n-1)}}{(n-1)!}\\frac{r^2}{r^2-n+1}. \\]<\/p>\n<p style=\"text-align: justify;\">This leads to the result by taking \\( {r^2:=|\\sqrt{n}z|^2} \\). \u2610<\/p>\n<p style=\"text-align: justify;\">The real Gaussian version of Theorem <a href=\"#thcomplex-ginibre-circle-weak\">4<\/a> was established by <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=1231689\">Edelman<\/a> (theorem 6.3). The sequence \\( {(H_k)_{k\\in\\mathbb{N}}} \\) forms an orthonormal basis of square integrable analytic functions on \\( {\\mathbb{C}} \\) for the standard Gaussian on \\( {\\mathbb{C}} \\). The uniform law on the unit disc (known as the circular law) is the law of \\( {\\sqrt{V}e^{2i\\pi W}} \\) where \\( {V} \\) and \\( {W} \\) are i.i.d.\u00a0uniform random variables on the \\( {[0,1]} \\). <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2446909\">Ledoux<\/a> makes use of this point of view for his interpolation between complex Ginibre and GUE via the Girko elliptic laws (see also <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2288065\">Johansson<\/a>).<\/p>\n<blockquote style=\"background: white; border: solid thick #e4e5e7; text-align: justify; padding-left: 1em;\"><p><b>Theorem 5 (Strong Circular Law)<\/b> <em><a id=\"thcomplex-ginibre-circle-strong\"><\/a> Almost surely, for all continuous bounded \\( {f:\\mathbb{C}\\rightarrow\\mathbb{R}} \\),<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\lim_{n\\rightarrow\\infty}\\int\\!f\\,d\\mu_{\\frac{1}{\\sqrt{n}}\\mathbf{G}} = \\pi^{-1}\\int_{|z|\\leq 1}\\!f(z)\\,dxdy. \\]<\/em><\/p>\n<\/blockquote>\n<p style=\"text-align: justify;\"><em>Proof:<\/em> We reproduce Silverstein's argument, published by <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=841088\">Hwang<\/a>. The argument is similar to the proof of the strong law of large numbers for i.i.d. random variables with finite fourth moment. It suffices to establish the result for compactly supported continuous bounded functions. Let us pick such a function \\( {f} \\) and set<\/p>\n<p style=\"text-align: center;\">\\[ S_n:=\\int_{\\mathbb{C}}\\!f\\,d\\mu_{\\frac{1}{\\sqrt{n}}\\mathbf{G}} \\quad\\text{and}\\quad S_\\infty:=\\pi^{-1}\\int_{|z|\\leq 1}\\!f(z)\\,dxdy. \\]<\/p>\n<p style=\"text-align: justify;\">Suppose for now that we have <a id=\"eqmom4\"><\/a><\/p>\n<p style=\"text-align: center;\">\\[ \\mathbb{E}[\\left(S_n-\\mathbb{E} S_n\\right)^4]=O(n^{-2}). \\ \\ \\ \\ \\ (2) \\]<\/p>\n<p style=\"text-align: justify;\">By monotone convergence,<\/p>\n<p style=\"text-align: center;\">\\[ \\mathbb{E}\\sum_{n=1}^\\infty\\left(S_n-\\mathbb{E} S_n\\right)^4 =\\sum_{n=1}^\\infty\\mathbb{E}[\\left(S_n-\\mathbb{E} S_n\\right)^4]&lt;\\infty \\]<\/p>\n<p style=\"text-align: justify;\">and consequently \\( {\\sum_{n=1}^\\infty\\left(S_n-\\mathbb{E} S_n\\right)^4&lt;\\infty} \\) a.s. which implies \\( {\\lim_{n\\rightarrow\\infty}S_n-\\mathbb{E} S_n=0} \\) a.s. Since \\( {\\lim_{n\\rightarrow\\infty}\\mathbb{E} S_n=S_\\infty} \\) by theorem <a href=\"#thcomplex-ginibre-circle-weak\">4<\/a>, we get \\( {\\lim_{n\\rightarrow\\infty}S_n=S_\\infty} \\) a.s. Finally, one can swap the universal quantifiers on \\( {\\omega} \\) and \\( {f} \\) thanks to the separability of \\( {\\mathcal{C}_c(\\mathbb{C},\\mathbb{R})} \\). To establish <a href=\"#eqmom4\">(2)<\/a>, we set<\/p>\n<p style=\"text-align: center;\">\\[ S_n-\\mathbb{E} S_n=\\frac{1}{n}\\sum_{i=1}^nZ_i \\quad\\text{with}\\quad Z_i:=f\\left(\\lambda_i\\left(\\frac{1}{\\sqrt{n}}\\mathbf{G}\\right)\\right). \\]<\/p>\n<p style=\"text-align: justify;\">Next, we obtain, with \\( {\\sum_{i_1,\\ldots}} \\) running over distinct indices in \\( {1,\\ldots,n} \\),<\/p>\n<p style=\"text-align: center;\">\\[ \\mathbb{E}\\left[\\left(S_n-\\mathbb{E} S_n\\right)^4\\right] =\\frac{1}{n^4}\\sum_{i_1}\\mathbb{E}[Z_{i_1}^4] \\]<\/p>\n<p style=\"text-align: center;\">\\[ +\\frac{4}{n^4}\\sum_{i_1,i_2}\\mathbb{E}[Z_{i_1}Z_{i_2}^3] \\]<\/p>\n<p style=\"text-align: center;\">\\[ +\\frac{3}{n^4}\\sum_{i_1,i_2}\\mathbb{E}[Z_{i_1}^2Z_{i_2}^2] \\]<\/p>\n<p style=\"text-align: center;\">\\[ +\\frac{6}{n^4}\\sum_{i_1,i_2,i_3}\\mathbb{E}[Z_{i_1}Z_{i_2}Z_{i_3}^2] \\]<\/p>\n<p style=\"text-align: center;\">\\[ +\\frac{1}{n^4}\\sum_{i_1,i_2,i_3,i_3,i_4}\\!\\!\\!\\!\\!\\!\\mathbb{E}[Z_{i_1}Z_{i_3}Z_{i_3}Z_{i_4}]. \\]<\/p>\n<p style=\"text-align: justify;\">The first three terms of the right hand side are \\( {O(n^{-2})} \\) since \\( {\\max_{1\\leq i\\leq n}|Z_i|\\leq\\Vert f\\Vert_\\infty} \\). Finally, some calculus using the expressions of \\( {\\varphi_{n,3}} \\) and \\( {\\varphi_{n,4}} \\) provided by Theorem <a href=\"#thcomplex-ginibre-k-points\">3<\/a> allows to show that the remaining two terms are also \\( {O(n^{-2})} \\). See <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=841088\">Hwang<\/a> (page 151). \u2610<\/p>\n<p style=\"text-align: justify;\">Following <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=1148410\">Kostlan<\/a>, the integration of the phases in the joint density of the spectrum given by Theorem <a href=\"#thcomplex-ginibre-density\">2<\/a> leads to theorem <a href=\"#thcomplex-ginibre-module\">6<\/a> below. See also <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=1986426\">Rider<\/a> and the <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2552864\">book by Hough, Krishnapur, Peres, and Virag<\/a> for a generalization to determinental processes.<\/p>\n<blockquote style=\"background: white; border: solid thick #e4e5e7; text-align: justify; padding-left: 1em;\"><p><b>Theorem 6 (Layers)<\/b> <em><a id=\"thcomplex-ginibre-module\"><\/a> If \\( {Z_{1},\\ldots,Z_{n}} \\) are independent with \\( {Z_k^2\\sim\\Gamma(k,1)} \\) for every \\( {k} \\), then<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\left\\{|\\lambda_1(\\mathbf{G})|,\\ldots,|\\lambda_n(\\mathbf{G})|\\right\\} \\overset{d}{=} \\left\\{Z_{1},\\ldots,Z_{n}\\right\\}. \\]<\/em><\/p>\n<\/blockquote>\n<p style=\"text-align: justify;\">Note that \\( {(\\sqrt{2}Z_k)^2\\sim\\chi^2(2k)} \\) which is useful for \\( {\\sqrt{2}\\mathbf{G}} \\). Since \\( {Z_k^2\\overset{d}{=}E_1+\\cdots+E_k} \\) where \\( {E_1,\\ldots,E_k} \\) are i.i.d. exponential random variables of unit mean, we get, for every \\( {r&gt;0} \\),<\/p>\n<p style=\"text-align: center;\">\\[ \\mathbb{P}\\left(\\rho(\\mathbf{G})\\leq \\sqrt{n}r\\right) =\\prod_{1\\leq k\\leq n}\\mathbb{P}\\left(\\frac{E_1+\\cdots+E_k}{n}\\leq r^2\\right) \\]<\/p>\n<p style=\"text-align: justify;\">where \\( {\\rho(\\mathbf{G})=\\max_{1\\leq i\\leq n}|\\lambda_i(\\mathbf{G})|=|\\lambda_1(\\mathbf{G})|} \\) is the <em>spectral radius<\/em> of \\( {\\mathbf{G}} \\).<\/p>\n<p style=\"text-align: justify;\">The law of large numbers suggests that \\( {r=1} \\) is a critical value. The central limit theorem suggests that \\( {n^{-1\/2}\\rho(\\mathbf{G})} \\) behaves when \\( {n\\gg1} \\) as the maximum of a standard Gaussian i.i.d\u00a0sample, for which the fluctuations follow the Gumbel law. A quantitative central limit theorem and the Borel-Cantelli lemma provides the follow result. The full proof is in <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=1986426\">Rider<\/a>.<\/p>\n<blockquote style=\"background: white; border: solid thick #e4e5e7; text-align: justify; padding-left: 1em;\"><p><b>Theorem 7 (Convergence and fluctuation of the spectral radius)<\/b> <a id=\"thcomplex-ginibre-spectral-radius\"><\/a><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\mathbb{P}\\left(\\lim_{n\\rightarrow\\infty}\\frac{1}{\\sqrt{n}}\\rho(\\mathbf{G})=1\\right)=1. \\]<\/em><\/p>\n<p><em>Moreover, if \\( {\\gamma_n:=\\log(n\/2\\pi)-2\\log(\\log(n))} \\) then<\/em><\/p>\n<p style=\"text-align: center;\"><em>\\[ \\sqrt{4n\\gamma_n} \\left(\\frac{1}{\\sqrt{n}}\\rho(\\mathbf{G})-1-\\sqrt{\\frac{\\gamma_n}{4n}}\\right) \\overset{d}{\\underset{n\\rightarrow\\infty}{\\longrightarrow}} \\mathcal{G}. \\]<\/em><\/p>\n<p><em>where \\( {\\mathcal{G}} \\) is the Gumbel law with cumulative distribution function \\( {x\\mapsto e^{-e^{-x}}} \\) on \\( {\\mathbb{R}} \\).<\/em><\/p><\/blockquote>\n<p style=\"text-align: justify;\">The convergence of the spectral radius was obtained by <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2129906\">Mehta<\/a> (chapter 15 page 271 equation 15.1.27) by integrating the joint density of the spectrum of Theorem <a href=\"#thcomplex-ginibre-density\">2<\/a> over the set \\( {\\bigcap_{1\\leq i\\leq n}\\{|\\lambda_i|&gt;r\\}} \\). The same argument is reproduced by <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=841088\">Hwang<\/a> (pages 149-150). Let us give now an alternative derivation of Theorem <a href=\"#thcomplex-ginibre-circle-weak\">4<\/a>. From Theorem <a href=\"#thcomplex-ginibre-spectral-radius\">7<\/a>, the sequence \\( {(\\mathbb{E}\\mu_{n^{-1\/2}\\mathbf{G}})_{n\\geq1}} \\) is tight and every adherence value \\( {\\mu} \\) is supported in the unit disc. From Theorem <a href=\"#thcomplex-ginibre-density\">2<\/a>, such a \\( {\\mu} \\) is rotationally invariant, and from Theorem <a href=\"#thcomplex-ginibre-module\">6<\/a>, the image of \\( {\\mu} \\) by \\( {z\\in\\mathbb{C}\\mapsto|z|} \\) has density \\( {r\\mapsto 2r\\mathbf{1}_{[0,1]}(r)} \\) (use moments!). Theorem <a href=\"#thcomplex-ginibre-circle-weak\">4<\/a> follows immediately.<\/p>\n<p style=\"text-align: justify;\">Note that the <a href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2641363\">recent book of Forrester<\/a> contains many computations for the Ginibre Ensemble. See also the <a href=\"http:\/\/arxiv.org\/abs\/0911.5645\">chapter by Khoruzhenko and Sommers in the Oxford Handbook of Random Matrices<\/a>.<\/p>\n<figure id=\"attachment_9485\" aria-describedby=\"caption-attachment-9485\" style=\"width: 225px\" class=\"wp-caption alignright\"><a href=\"https:\/\/en.wikipedia.org\/wiki\/Jean_Ginibre\"><img loading=\"lazy\" class=\"wp-image-9485 size-medium\" src=\"http:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/jean.ginibre-225x300.jpg\" alt=\"Photo of Jean Ginibre\" width=\"225\" height=\"300\" srcset=\"https:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/jean.ginibre-225x300.jpg 225w, https:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/jean.ginibre-113x150.jpg 113w, https:\/\/djalil.chafai.net\/blog\/wp-content\/uploads\/2010\/11\/jean.ginibre.jpg 426w\" sizes=\"(max-width: 225px) 100vw, 225px\" \/><\/a><figcaption id=\"caption-attachment-9485\" class=\"wp-caption-text\">Jean Ginibre - Mathematical physicist<\/figcaption><\/figure>\n","protected":false},"excerpt":{"rendered":"<p>This post is devoted to the Complex Ginibre Ensemble, the subject of an expository talk that I gave few months ago. Let \\( {(G_{i,j})_{i,j\\geq1}} \\)&#8230;<\/p>\n<div class=\"more-link-wrapper\"><a class=\"more-link\" href=\"https:\/\/djalil.chafai.net\/blog\/2010\/11\/02\/aspects-of-the-complex-ginibre-ensemble\/\">Continue reading<span class=\"screen-reader-text\">Aspects of the Complex Ginibre Ensemble<\/span><\/a><\/div>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"iawp_total_views":2419},"categories":[1],"tags":[],"_links":{"self":[{"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/posts\/969"}],"collection":[{"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/comments?post=969"}],"version-history":[{"count":28,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/posts\/969\/revisions"}],"predecessor-version":[{"id":9491,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/posts\/969\/revisions\/9491"}],"wp:attachment":[{"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/media?parent=969"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/categories?post=969"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/tags?post=969"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}