# A Renormalisation Group Approach to the Universality of Wigner's Semicircle Law for Random Matrices with Dependent Entries.

1. Introduction

Random matrix theory (see the classical text [1]) first appeared in physics in Wigner's work on the level spacing in large nuclei. Since then, it has proven to have multiple applications to physics and other branches of science (see, e.g., [2]). Most of these applications rely on the universal behaviour of some of the observables for matrices of large size. A simple example is Wigner's semicircle law for the eigenvalue density that holds in the large N limit for matrices whose entries are independent and identically distributed.

Understanding the universal behaviour of eigenvalue distributions and correlations ranks among the major problems in random matrix theory. In this respect, the renormalisation group turns out to be a powerful technique. Introduced in the context of critical phenomena in statistical mechanics by K. Wilson to account for the universality of critical exponents, the latter has also been proven to be useful in understanding probability theory. For instance, it leads to an insightful proof of the central limit theorem (see the review by Jona-Lasinio [3] and references therein).

The renormalisation group has been used to derive the semicircle law for random matrices in the pioneering work of Brezin and Zee [4]. In the latter approach, the renormalisation group transformation consists in integrating over the last line and column of a matrix of size N+1 to reduce it to a matrix of size N. This leads to a differential equation for the resolvent G(z) = 1 /N<Tr[(z - M).sup.-1]> in the large N limit whose solution yields the semicircle law.

In this paper, we follow a different route: we first express the resolvent as an integral over replicas and introduce a differential equation for the replica effective action. This differential equation is a very simple analogue of Polchinski's exact renormalisation group equation [5]. It is used to derive inductive bounds on the various terms, ensuring that the semicircle law is obeyed provided the cumulants of the original matrix model fulfil some simple scaling bounds in the large N limit.

This paper is based on some work in collaboration with Krajewski et al. in which we extend Wigner's law to random matrices whose entries fail to be independent [6] to which we refer for further details. There have been other works on such an extension (see [7-9]).

2. What Are Random Matrices?

A random matrix is a probability law on a space of matrices, usually given by the joint probability density on its entries:

[rho](M) = [rho]([M.sub.11], [M.sub.12], ...). (1)

Thus, a random matrix of size N is defined as a collection of [N.sup.2] random variables. However, there is a much richer structure than this, relying notably on the spectral properties of the matrices.

Here, we restrict our attention to a single random matrix. Note that it is also possible to consider several random matrices, in which case the noncommutative nature of matrix multiplication plays a fundamental role, leading to the theory of noncommutative probabilities.

There are two important classes of probability laws on matrices:

(i) Wigner ensemble: the entries are all independent variables:

[mathematical expression not reproducible], (2)

up to the Hermitian condition [[bar.M].sub.ij] = [M.sub.ji].

(ii) Unitary ensemble: the probability law is invariant under unitary transformations:

[rho]([UMU.sup.[dagger]]) = [rho](M), (3)

for any unitary matrix U [member of] U(N).

The only probability laws that belong to both classes are the Gaussian ones:

[rho](M) [varies] exp - [1/2[[sigma].sup.2]] Tr[(M).sup.2], (4)

up to a shift of M by a fixed scalar matrix.

The main objects of interest are the expectation values of observables, defined as

<O> = [integral] dM[rho](M) O (M). (5)

Among the observables, the spectral observables defined as symmetric functions of the eigenvalues of M play a crucial role in many applications. This is essentially due to their universal behaviour: in the large N limit, for some matrix ensembles and in particular regimes, the expectation values of specific spectral observables do not depend on the details of the probability law [rho](M).

Universality is at the root of the numerous applications to physics and other sciences, since the results we obtain are largely model-independent. Among the applications to physics, let us quote the statistics of energy levels in heavy nuclei, disordered mesoscopic systems, quantum chaos, chiral Dirac operators, and so forth.

3. Wigner's Semicircle Law

In this paper, we focus on the eigenvalue density, defined as

[rho]([lambda]) = 1/N <[summation over (1[less than or equal to]i[less than or equal to]N)] [delta][[lambda] - [[lambda].sub.i](M/[square root of (N))]]>. (6)

In particular, a universal behaviour is expected in the large N limit for some ensembles.

For a Gaussian random Hermitian matrix [rho](M) [varies] exp - (1/2[[sigma].sup.2])Tr([M.sup.2]), the eigenvalue density obeys Wigner's semicircle law:

[mathematical expression not reproducible]. (7)

Empirically, [rho]([lambda]) may be determined by plotting the histogram of eigenvalue of a matrix taken at random with a given probability law (see Figure 1).

The derivation of Wigner's semicircle in the large N limit is based on the resolvent (also known as the Green function):

[mathematical expression not reproducible]. (8)

Then, the density of eigenvalues is recovered as

[rho]([lambda]) = G([lambda] - [i0.sup.+]) - G([lambda] + [i0.sup.+])/2i[pi], (9)

where we have used the relation

1/x [+ or -] [i0.sup.+] = p.v. 1/x [- or +] 2i[pi][delta](x). (10)

In the large N limit, for the Gaussian model, the resolvent obeys the self-consistency equation (also known as the Schwinger-Dyson equation) (see, e.g., [10], Section VII.4):

G(z) = [[infinity].summation over (k=0)] [[sigma].sup.k][G.sup.k](z)/[z.sup.k+1] = 1/z - [sigma]G(z). (11)

Its solution that behaves as 1/z for large z is

G(z) = z/2[[sigma].sup.2](1 - [square root of (1 - 4[[sigma].sup.2]/[z.sup.2]))]. (12)

Taking the cut of the square root on the negative real axis, we obtain Wigner's semicircle law (7) in the large N limit.

The semicircle law is not limited to the Gaussian case; it also holds for Wigner matrices in the large N limit. A random Hermitian N x N matrix is a Wigner matrix if

(i) real and imaginary parts of upper diagonal elements are independent and identically distributed (i.i.d.) with mean 0 and variance [sigma];

(ii) diagonal elements are i.i.d. with finite mean and variance and independent of the off-diagonal ones.

Then, in the limit N [right arrow] +[infinity], the eigenvalue distribution of M/[square root of N] is the semicircle law (7).

The original proof is of combinatorial nature and involves the expectation of the moments:

[mathematical expression not reproducible]. (13)

To derive this result, the idea is to first factorise [rho] for a Wigner ensemble as

[mathematical expression not reproducible], (14)

where [rho]' is the common probability density of the real diagonal terms and [rho]" is the common probability density of the real and imaginary parts of the off-diagonal terms.

Then, we expand the trace and integrate over the independent real variables [M.sub.ii], Re [M.sub.ij], and Im [M.sub.ij]. The power of N in the expectation of a given moment arises from the denominator 1/[N.sup.k/2+1] and from the number of independent indices in the summations. In the large N limit, the only configurations that survive are counted by Catalan numbers: [C.sub.l] = (2l)!/[(l!).sup.2](l + 1). Since the latter also appear in the Taylor expansion,

z/2[[sigma].sup.2](1 - [square root of (1 - 4[[sigma].sup.2]/[z.sup.2])] = [summation over (l[greater than or equal to]0)] [(2l)!/[(l!).sup.2](l + 1)] [[[sigma].sup.2l]/[z.sup.2l+1]], (15)

we conclude that

[mathematical expression not reproducible]. (16)

This is the form of the resolvent which leads to Wigner's semicircle law. Here, we see universality at work: in the large N limit, the eigenvalue density is given by the semicircle law, whatever the probability densities [rho]' and [rho]" are. However, this result relies on the independence of the matrix elements. In the next section, we will extend it to matrices whose entries are not necessary independent.

4. Wigner's Law beyond Wigner Ensembles

Let us introduce the cumulants, defined through their generating function:

[mathematical expression not reproducible]. (17)

In the physics terminology, these are the connected correlation functions. In particular, the Gaussian cumulants vanish beyond the quadratic term:

[mathematical expression not reproducible]. (18)

Therefore, cumulants of degree higher than 2 are a measure of the deviation from the Gaussian case.

Turning back to the general case, for each cumulant, we construct an oriented graph as follows (see Figure 2 for some examples):

(i) Vertices are distinct matrix indices in the cumulant.

(ii) There is an edge from i to j for every [M.sub.ij].

Since nonquadratic cumulants measure deviations from the Gaussian case, if the perturbation is small, it is reasonable to expect that the semicircle law is still obeyed.

To state this result, recall that an oriented graph is Eulerian if every vertex has an equal number of incoming and outgoing edges. Equivalently, it means that every connected component admits an Eulerian cycle, that is, an oriented cycle that passes through all edges, respecting the orientation. Furthermore, let us denote by v(G), e(G), and c(G) the numbers of vertices, edges, and connected components of G, respectively.

Theorem 1 (Wigner's law for matrices with dependent entries). Let [[rho].sub.N] be a probability law on the space of Hermitian N x N matrices M such that its cumulants can be decomposed as [C.sub.G] = [C'.sub.G] + [C".sub.G], with [C'.sub.G] being a Gaussian cumulant and [C".sub.G] being a perturbation such that, uniformly in the vertex indices [i.sub.1], ..., [i.sub.v(G)] (i.e., all constants involved should not depend on these indices),

(i) [lim.sub.N[right arrow][infinity]] [N.sup.v(G)-c(G)-e(G)/2][C.sub."G]([i.sub.1],...,[i.sub.v(G)]) = 0 if G is Eulerian,

(ii) [N.sup.v{G)-c{G)-e{G)/2] [C".sub.G]([i.sub.1],...,[i.sub.v(G)]) is bounded if G is not Eulerian.

Then, the moments of the eigenvalue distribution of the matrix M/[square root of (N)] converge towards the moments of the semicircle law, with a given by the Gaussian cumulant [<[M.sub.ij][M.sub.kl]>.sub.c] = [[sigma].sup.2][[delta].sub.il][[delta].sub.jk]:

[mathematical expression not reproducible]. (19)

For instance, for the graph

which is not Eulerian, with v = 3, e = 4, and c = 2, the cumulant should obey

1/N |[<[([M.sub.ij]).sup.2] [M.sub.jk][M.sub.ll]>.sub.c]] [less than or equal to] K, (20)

with K being a constant that does not depend on the indices i, j, k, and l. On the other hand, for the graph

which is Eulerian, with v = 2, e = 4, and c = 1, we impose

[mathematical expression not reproducible] (21)

uniformly in i and j.

As an illustration, we recover the case of Wigner matrices (with finite moments). Indeed,

(i) there is no graph with v [greater than or equal to] 3 (independence of off diagonal matrix elements);

(ii) for v = 1 and v = 2, e [greater than or equal to] 3, bounds are satisfied because of 1/[N.sup.e/2] and all moments are assumed to be finite;

(iii) [mathematical expression not reproducible] (off diagonal elements have mean value 0);

(iv) [mathematical expression not reproducible] (independence of diagonal elements);

(v) [mathematical expression not reproducible] (independence of diagonal and off diagonal elements);

(vi) [mathematical expression not reproducible] (independence of real and imaginary parts and equality of their distributions with mean value 0);

(vii) [mathematical expression not reproducible] is the Gaussian cumulant leading to the semicircle law.

The case of unitarily invariant matrices is critical since the bounds are saturated (see [6]). This is consistent since we know that the semicircle law is not obeyed by unitary non-Gaussian ensembles [11].

It is possible to give a combinatorial proof of this result based on the relation between moments and cumulants:

[mathematical expression not reproducible]. (22)

In the moment method, we have to estimate

[mathematical expression not reproducible]. (23)

Then, we express the moments in (23) in terms of cumulants using (22) and represent each cumulant as a graph. Because of the trace, one has to draw Eulerian cycles on the graphs after some vertex identifications. Then, the scaling bounds on the cumulants can be used to show that only Gaussian terms survive.

5. Proof Based on the Replica Effective Action

Let us give a renormalisation group proof of this result based on the replica effective action. The use of replicas in random matrix theory is a classical subject (see, e.g., [12] or [13]). To begin with, let us note that

Tr[(z - M/[square root of (N)]).sup.-1] = [partial derivative]/[partial derivative]z log det (z - M/[square root of (N)]). (24)

It is convenient to express the logarithm using the replica method. First, observe that

[mathematical expression not reproducible]. (25)

Then, we express the nth power of the determinant as a Gaussian integral over n replicas of a complex vector of size N (with a factor of [[pi].sup.nN] included in the measure):

1/[det.sup.n](z - M) = [integral] dX exp - Tr([X.sup.[dagger]](z - M)X), (26)

which fit into an N x n complex matrix [mathematical expression not reproducible].

The limit n [right arrow] 0 may be worrisome; its meaning is as follows. Because of U(n) invariance, any perturbative result in powers of 1/z is a polynomial in n, from which we retain only the linear term. Of course, this may not hold beyond perturbation theory, where replica symmetry breaking can occur.

After averaging over M with the random matrix density [rho](M), we obtain the following expression for the resolvent:

[mathematical expression not reproducible], (27)

where the replica potential is

[V.sub.0](X) = log <exp Tr ([X.sup.[dagger]] [M/[square root of (N)]]X)>. (28)

Because of the logarithm, the potential involves the cumulants and can be expanded over graphs as

[mathematical expression not reproducible], (29)

where s(e) is the source of edge e and t(e) is its target.

Let us introduce a replica effective action, obtained by a partial integration:

V(t, X) = log [integral] dY exp {- Tr([T.sup.[dagger]]Y)/t + [V.sub.0](X + Y)} - N n log t. (30)

The parameter t ranges between 0 (where we have no integration, V(t = 0, X) = [V.sub.0](X)) and t = 1/z.

The effective potential obeys a semigroup property that follows from Gaussian convolution (see, e.g., [14], Section A10.1):

V(t + s, X) = log [integral] dY exp {Tr([Y.sup.[dagger]]Y(/s + V(t, X + Y)} - N n log s. (31)

For small s = dt, it translates into the following renormalisation group equation, which is a simple version of Polchinski's exact renormalisation group equation [5]:

[mathematical expression not reproducible]. (32)

The first term on the RHS is referred to as the loop term, since it creates a new loop in the Feynman graph expansion of the effective action, while the second inserts a one-particle reducible line and is referred to as the tree term (see Figure 3).

Taking into account the boundary condition V(t = 0, X) = [V.sub.0](X), it is convenient to write (32) in integral form:

[mathematical expression not reproducible]. (33)

This allows us to derive inductive bounds in powers of t = 1/z.

From a physical point of view, we evaluate the effective potential by a large succession of small partial integrations, with a total weight given by t. Let us stress that, in our context, this differential equation is merely a tool to control the t dependence of the effective action after integrating with a t-dependent propagator.

The effective potential also admits an expansion over graphs:

[mathematical expression not reproducible]. (34)

This leads to a graphical interpretation of the action of the two differential operators in the renormalisation group equation (see Figure 4). Indeed, in the expansion (see (34)), an edge joining a vertex carrying label i to a vertex carrying j is equipped with a factor [[summation].sub.a] [X.sub.i,a] [[bar.X].sub.j,a], with a being a replica index. Then, the differential operator [partial derivative]/[partial derivative][X.sub.i,a] (resp., [partial derivative]/[partial derivative][[bar.X].sub.j,a]) removes the outgoing (resp., incoming) half edge. Finally, the remaining half edges are reattached and the vertices identified to yield a new graph on the RHS of (33), with one less edge. These operations are performed on the same graph for the loop term and on distinct ones for the tree term.

Let us decompose the effective cumulants appearing in (34) into Gaussian ones and perturbations and expand both in a power series in t = 1/z:

[mathematical expression not reproducible]. (35)

The Gaussian terms are those that are constructed using only the Gaussian term in the initial potential [V.sub.0](X). Even if [V.sub.0](X) is quartic in X, this does not hold for the Gaussian part of [V.sub.t](X), which contains terms of all orders. The perturbation collects all the remaining terms; they contain at least one non-Gaussian perturbation from [V.sub.0](X).

The renormalisation group equation (33) allows us to prove inductively on k that the perturbations [C".sup.(k).sub.G] obey the same scaling bound imposed on [C".sup.(0).sub.G] = [C".sup.(0).sub.G] and that the purely Gaussian terms do not grow too fast.

(i) [mathematical expression not reproducible] if G is Eulerian.

(ii) [mathematical expression not reproducible] is bounded if G is not Eulerian.

(iii) [mathematical expression not reproducible] is bounded for any G.

This involves a combinatorial discussion based on the graphical interpretation of Figure 4 which can be found in [6]. Let us simply mention that the terms that may violate the bounds are of higher order in n. Thus, they are harmless when taking the limit n [right arrow] 0 before the limit N [right arrow] +[infinity].

Finally, using (27) and the renormalisation group equation (33), the resolvent can be expressed as

[mathematical expression not reproducible]. (36)

The scaling bounds for the non-Gaussian cumulants impose, perturbatively in 1/z,

[mathematical expression not reproducible]. (37)

Therefore, only the Gaussian cumulants contribute and we recover Wigner's semicircle law.

6. Conclusion and Outlook

In this paper, we have argued that Wigner's semicircle law remains valid for matrices with dependent entries. The deviation from the independent case is measured by the joint cumulants of the entries, which are assumed to fulfil some scaling bound for large N. To establish this result, we have introduced an effective action for the replicas. This effective action obeys a renormalisation group equation that allowed us to prove perturbative bounds on the effective cumulants. As a consequence of these bounds, only the Gaussian terms contribute in the large N limit, thus establishing the validity of Wigner's semicircle law.

It may also be of interest to investigate the case of the sum of a random matrix M and a deterministic one A (see, e.g., [13], where such a model is discussed). In this case, the resolvent is expressed as

[mathematical expression not reproducible]. (38)

In our context, the deterministic matrix A induces a nontrivial kinetic for the replicas. In particular, if A is a discrete Laplacian, it yields a nontrivial renormalisation group flow that bears some similarities with the QFT renormalisation group. In this case, we expect to exploit the true power of the renormalisation group equation, with a discussion of fixed points and scaling dimensions.

https://doi.org/10.1155/2017/4098720

Disclosure

This paper was presented at the 5th Winter Workshop on Non-Perturbative Quantum Field Theory, Sophia Antipolis, March 2017.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author wishes to thank the organisers of the workshop, Mario Gattobigio, Thierry Grandou, and Ralf Hoffmann, for their kind invitation, as well as Stan Brodsky for insightful remarks. The author is also grateful to his collaborators, Dinh Long Vu and Adrian Tanasa.

References

[1] M. L. Mehta, Random Matrices, Academic Press, 2004.

[2] G. Akemann, J. Baik, and Ph. Di Francesco, Eds., Oxford Handbook of Random Matrix Theory, Oxford Univ. Press, 2011.

[3] G. Jona-Lasinio, "Renormalisation group and probability theory Phys," Physics Reports, vol. 352, article 0009219, 439 pages, 2001, https://arxiv.org/abs/cond-mat/0009219.

[4] E. Brezin and A. Zee, "Universality of the correlations between eigenvalues of large random matrices," Nuclear Physics B, vol. 402, no. 3, pp. 613-627, 1993.

[5] J. Polchinski, "Renormalization and effective lagrangians," Nuclear Physics B, vol. 231, no. 2, pp. 269-295, 1984.

[6] T. Krajewski, A. Tanasa, and D. L. Vu, "Wigner law for matrices with dependent entries: A perturbative approach," Journal of Physics A: Mathematical and Theoretical, vol. 50, no. 16, 2017, https://arxiv.org/abs/1609.01873.

[7] J. H. Schenker and H. Schulz-Baldes, "Semicircle law and freeness for random matrices with symmetries or correlations," Mathematical Research Letters, vol. 12, no. 4, pp. 531-542, 2005.

[8] K. Hofmann-Credner and M. Stolz, "Wigner theorems for random matrices with dependent entries: ensembles associated to symmetric spaces and sample covariance matrices," Electronic Communications in Probability, vol. 13, pp. 401-414, 2008.

[9] O. Friesen and M. Lowe, "The semicircle law for matrices with dependent entries".

[10] A. Zee, Quantum Field Theory in a Nutshell, Princeton University Press, Princeton, NJ, USA, 2nd edition, 2010.

[11] E. Brezin, C. Itzykson, G. Parisi, and J. B. Zuber, "Planar diagrams," Communications in Mathematical Physics, vol. 59, no. 1, pp. 35-51, 1978.

[12] E. Kanzieper, "Replica Approach in Random Matrix Theory," in The Oxford Handbook of Random Matrix Theory, Oxford Univ. Press, Oxford, UK, 2011.

[13] A. Zee, "Law of addition in random matrix theory," Nuclear Physics B, vol. 474, no. 3, pp. 726-744, 1996.

[14] J. Zinn-Justin, "Quantum field theory and critical phenomena," in International Series of Monographs on Physics, vol. 113, Oxford University Press, 2002.

Thomas Krajewski

Centre de Physique Theorique, Aix-Marseille Universite, Marseille, France

Correspondence should be addressed to Thomas Krajewski; krajewsk@gmail.com

Received 15 June 2017; Revised 1 September 2017; Accepted 1 October 2017; Published 31 December 2017

Caption: Figure 1: Spectrum of a large (N = 1000) Hermitian matrix with a Gaussian distribution ([sigma] = 1).

Caption: Figure 2: Examples of graph associated with cumulants.

Caption: Figure 3: Graphical interpretation of the renormalisation group equation.

Caption: Figure 4: Action of the differential operators on the vertices of the effective action.