{"id":40,"date":"2010-05-12T07:52:16","date_gmt":"2010-05-12T05:52:16","guid":{"rendered":"http:\/\/djalil.chafai.net\/blog\/?p=40"},"modified":"2019-11-10T19:13:33","modified_gmt":"2019-11-10T18:13:33","slug":"circular-law-theorem-for-random-markov-matrices","status":"publish","type":"post","link":"https:\/\/djalil.chafai.net\/blog\/2010\/05\/12\/circular-law-theorem-for-random-markov-matrices\/","title":{"rendered":"Circular law theorem for random Markov matrices"},"content":{"rendered":"<p style=\"text-align: justify;\">I have uploaded on arXiv a <a title=\"Circular Law Theorem for Random Markov Matrices\" href=\"http:\/\/arxiv.org\/abs\/0808.1502\">paper<\/a> written with <a title=\"Charles Bordenave\" href=\"\/scripts\/search.php\/?q=Charles+Bordenave\">Charles Bordenave<\/a> and <a title=\"Pietro Caputo\" href=\"\/scripts\/search.php\/?q=Pietro+Caputo\">Pietro Caputo<\/a><em> entitled Circular Law Theorem for Random Markov Matrices<\/em>. We provide, probably for the first time, an almost sure\u00a0 circular law theorem for a model of random matrices with dependent entries and finite positive variance. This work leads us to formulate some open problems of variable difficulty, listed at the end of the introduction. Feel free to solve them!<\/p>\n<p style=\"text-align: justify;\">To be a bit more precise, let $(X_{jk})_{jk\\geq1}$ be i.i.d. nonnegative random variables with bounded density, mean $m$, and finite positive variance $\\sigma^2$. Let $M$ be the $n\\times n$ random Markov matrix with i.i.d. rows defined by $M_{jk}=X_{jk}\/(X_{j1}+\\cdots+X_{jn})$.\u00a0 Let $\\lambda_1,\\ldots,\\lambda_n$ be the eigenvalues of $\\sqrt{n}M$ i.e. the roots in $\\mathbb{C}$ of its characteristic polynomial. Our main result states that with probability one, the counting probability measure $\\frac{1}{n}\\delta_{\\lambda_1}+\\cdots+\\frac{1}{n}\\delta_{\\lambda_n}$ converges weakly as $n\\to\\infty$ to the uniform law on the disk $\\{z\\in\\mathbb{C}:|z|\\leq m^{-1}\\sigma\\}$ (circular law).<\/p>\n<p style=\"text-align: justify;\">In particular, when $X_{11}$ follows an exponential law, the  random matrix $M$ belongs to the Dirichlet Markov Ensemble of  random stochastic matrices, and our main result gives a positive answer to the question  raised in a <a title=\"The Dirichlet Markov Ensemble \" href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2575404\">previous  paper<\/a>.<\/p>\n<p style=\"text-align: justify;\">The matrix $M$ is Markov: its entries belong to $[0,1]$ and each row sums up to $1$. It has equally distributed dependent entries. However, the rows of $M$ are i.i.d. and follow an exchangeable law on $\\mathbb{R}^n$ supported by the simplex. Since $\\sigma&lt;\\infty$, the uniform law of large numbers of Bai and Yin states that almost surely<\/p>\n<p style=\"text-align: justify;\">$\\displaystyle\\max_{1\\leq i\\leq n}\\left|X_{i1}+\\cdots+X_{in}-nm\\right|=o(n)$.<\/p>\n<p style=\"text-align: justify;\">This suggests that $m^{-1}\\sqrt{n}M$ is approximately equal to $n^{-1\/2}X$ for $n$ large enough. One can then expect that almost surely the counting probability measure of the singular values of $\\sqrt{n}M$ will converge as $n\\to\\infty$ to\u00a0 a quartercircular law while the counting probability measure of the eigenvalues of $\\sqrt{n}M$ will converges to a circular law. We show that this heuristics is valid. There is however a complexity gap between these two statements due to the fact that for nonnormal operators such as $M$, the eigenvalues are less stable than the singular values under perturbations (here the least singular value enters the scene).<\/p>\n<p style=\"text-align: justify;\">The  proof of the main result is inspired from the<a title=\"Random matrices:  Universality of ESDs and the circular law\" href=\"http:\/\/arxiv.org\/abs\/0807.4898\"> Tao and Vu  version of the  Girko Hermitization via logarithmic potentials<\/a>. More  precisely, the starting point is the following Hermitization identity  valid for any $n\\times n$ complex matrix $A$ and any $z$ outside the spectrum of $A$:<\/p>\n<p style=\"text-align: justify;\">$$\\displaystyle  \\frac{1}{n}\\sum_{k=1}^n\\log|\\lambda_k(A)-z|=\\frac{1}{n}\\sum_{k=1}^n\\log(\\lambda_k(\\sqrt{(A-zI)(A-zI)^*})).$$<\/p>\n<p style=\"text-align: justify;\">The left hand side is the logarithmic potential at point $z$ of the counting probability measure of the eigenvalues of $A$. The right hand side is the integral of $\\log(\\cdot)$ for the counting probability measure of the singular values of $A-zI$. Thanks to the <a title=\"Logarithmic potential and Hermitization\" href=\"\/blog\/2010\/05\/14\/logarithmic-potential-and-hermitization\/\">Girko Hermitization theorem<\/a>, the problem boils down then to show that for  Lebesgue almost all $z\\in\\mathbb{C}$, almost surely,<\/p>\n<ol>\n<li>the counting probability measure $\\nu_{\\sqrt{n}M-zI}$ of the  singular values of $\\sqrt{n}M-zI$ tends weakly as $n\\to\\infty$ to a deterministic probability measure $\\nu_z$ related  to the circular law<\/li>\n<li>the function $\\log(\\cdot)$ is uniformly integrable for the  sequence $(\\nu_{\\sqrt{n}M-zI})_{n\\geq1}$<\/li>\n<\/ol>\n<p style=\"text-align: justify;\">For the  first statement, we use a result of <a title=\"On the  empirical distribution of eigenvalues of large dimensional  information-plus-noise-type matrices.\" href=\"http:\/\/www.ams.org\/mathscinet-getitem?mr=2322123\">Dozier and  Silverstein<\/a>, which is based on Hermitian techniques: truncation, centralization,\u00a0 recursion for the Cauchy-Stieltjes transform (trace of the resolvent) via Sherman-Morrison bloc inversion, producing finally a cubic equation for the Cauchy-Stieltjes transform of the limiting law.<\/p>\n<p style=\"text-align: justify;\">The second statement\u00a0 involves as essential ingredients a Talagrand  concentration inequality (for the control of the distance of rows to the span of other rows), and a polynomial control of the operator norm  of the  resolvent (i.e. the least singular value of $\\sqrt{n}M-zI$). The  bounded density assumption  is purely technical and comes from the way we  control the operator  norm of the resolvent. A possible route for the  removal of this  assumption is to control the <a title=\"Previous post on  least singular  values for certain matrices with independent rows\" href=\"\/blog\/2010\/04\/27\/least-singular-value-of-random-matrices-with-independent-rows\/\">least   singular value for a certain kind of matrices with independent rows<\/a>.<\/p>\n<p style=\"text-align: justify;\"><img title=\"Superposition of la_2,...,la_n for 10 simulations of n^(1\/2)M with n=250 and X11 symmetric Bernoulli.\" src=\"\/blog\/wp-content\/uploads\/cirmar\/bern.jpg\" alt=\"\" \/><\/p>\n<p style=\"text-align: justify;\"><img title=\"Voronoi tessellation of la_2,...,la_n for a single simulation of n^(1\/2)M with n=250 and X11 symmetric Bernoulli.\" src=\"\/blog\/wp-content\/uploads\/cirmar\/bernv.jpg\" alt=\"\" \/><\/p>\n","protected":false},"excerpt":{"rendered":"<p>I have uploaded on arXiv a paper written with Charles Bordenave and Pietro Caputo entitled Circular Law Theorem for Random Markov Matrices. We provide, probably&#8230;<\/p>\n<div class=\"more-link-wrapper\"><a class=\"more-link\" href=\"https:\/\/djalil.chafai.net\/blog\/2010\/05\/12\/circular-law-theorem-for-random-markov-matrices\/\">Continue reading<span class=\"screen-reader-text\">Circular law theorem for random Markov matrices<\/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":93},"categories":[1],"tags":[],"_links":{"self":[{"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/posts\/40"}],"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=40"}],"version-history":[{"count":4,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/posts\/40\/revisions"}],"predecessor-version":[{"id":11758,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/posts\/40\/revisions\/11758"}],"wp:attachment":[{"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/media?parent=40"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/categories?post=40"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/djalil.chafai.net\/blog\/wp-json\/wp\/v2\/tags?post=40"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}