in this short introduction we would like to introduce the two players in the title , random matrix theory ( rmt ) and quantum chromodynamics ( qcd ) , on a superficial level .
this gives a motivation why and where it will be beneficial to relate these two seemingly unrelated theories , and what the reader may expect to learn from these notes .
let us begin with qcd , the theory of the strong interactions .
it is part of the standard model of elementary particles that also describes the weak and electromagnetic interactions among all elementary particles .
the latter two interactions will however not play any role here .
the particle content of qcd are the @xmath0 quarks with masses @xmath1 , @xmath2 , and the gluons . in nature @xmath3 flavours have been observed named up , down , strange , charm , top and bottom quark . in the following
we will keep @xmath0 as a free parameter , and often consider only the lightest up and down quark ( @xmath4 ) as they constitute the most common particles as proton ( p ) , neutron ( n ) and pions ( @xmath5 ) .
the quarks interact through the gluons , the carriers of force , with a coupling constant @xmath6 .
both come in 3 colours and the interaction is described through a field strength and covariant derivative carrying an @xmath7 lie group structure .
qcd is a strongly interacting , relativistic quantum field theory ( qft ) which is very difficult .
fortunately we will only need to know a few features and some of its global symmetries that will be described in section [ symm ] .
for a standard textbook on perturbative qcd we refer for example to @xcite .    roughly speaking qcd has two different phases that are schematically depicted in figure [ fig : phasediag ] .
they are characterised by an order parameter , the chiral condensate @xmath8 . at high energies corresponding to high temperatures t ,
@xmath9 and quarks and gluons form a plasma , that has been observed in collision experiments e.g. at the relativistic heavy ion collider rhic in brookhaven .
this phase also existed in the early universe , with the cooling down happening close to the temperature axis at low baryon density parametrised by @xmath10 . in these situations as well as in single particle collisions produced at collider experiments a perturbative expansion in powers of the coupling @xmath6 typically applies , as described e.g. in @xcite .
the second phase with @xmath11 at low temperature and density is the one in which we live . here
quarks condense into colourless objects , that is into baryons made of 3 quarks of different colour adding up to white ( like p or n ) , or mesons made of a quark and an anti - quark with its anti - colour ( like the @xmath12 s ) .
gluons are also confined to these objects , and a proof of confinement of quarks and gluons from qcd is still considered to be an open millennium problem . in qcd
it happens that confinement goes along with the spontaneous breaking of chiral symmetry - hence the name of @xmath8 .
this global symmetry will be described in detail later , and eventually lead us to a rmt description .
( 12.2,5.5 ) ( 1.0,0.3)(1,0)8 ( 9.2,0.2)@xmath10 ( 1.0,0.3)(0,1)5 ( 0.1,5.1)@xmath13 ( 0.1,4.3)@xmath14 ( 7.1,4.0)@xmath9 ( 1.0,-0.1)@xmath15 ( 2.9,1.8)@xmath11 ( 3.9,4.25)@xmath16 ( 8.0,-0.1)@xmath17 ( 1.0,4.4)(4.7,4.3)(4.0,4.3 ) ( 4.0,4.3)(7.0,4.2)(8.0,0.3 )    in the phase with @xmath11 perturbation theory breaks down and one has to apply other methods .
if one wants to continue to work with first principles and the qcd action , one possibility is to study qcd on a finite space - time lattice of volume @xmath18 and lattice spacing @xmath19 , equipped with a euclidean metric .
the numerical solution of this theory nowadays reproduces the masses of particles that are composed of elementary ones to a very high precision , as testified in the particle data booklet @xcite . to that aim
two limits have to be made , the continuum limit sending @xmath20 and the thermodynamical limit sending @xmath21 . for a standard textbook on lattice
qcd we refer to @xcite .
a second possibility is to approximate qcd in the confined phase by effective theories , that describe for example only the low energy excitation .
one of these is chiral perturbation theory ( chpt ) , see @xcite for a review , that describes the low momentum modes that appear after the spontaneous breaking of chiral symmetry , the so - called goldstone bosons .
another such theory is rmt .
it should be clear by now that rmt will not solve qcd ( nor chpt ) , in the continuum or on the lattice .
however , it will describe certain aspects in an analytic fashion , namely the spectral properties of the small eigenvalues of the qcd dirac operator @xmath22 . in a finite volume
the density of eigenvalues of @xmath22 satisfies the banks - casher relation @xcite @xmath23 which relates it to the described setup
. in particular rmt will predict the detailed dependence on the parameters @xmath1 , @xmath10 , @xmath19 , and @xmath18 after being appropriately rescaled , as well as on a the zero - eigenvalues of @xmath22 which relate to a topological index @xmath24 .
and most remarkably these predictions have been verified in comparison to lattice qcd by many groups .
we will not repeat these findings here and refer to the literature at the end of each section .
the idea to apply rmt to qcd goes back to the seminal works @xcite that started from the simplest case with @xmath25 which is called quenched approximation , at @xmath26 .
the field has since developed enormously , leading to the detailed analytical knowledge that we will describe .
there exist a number of excellent reviews already , @xcite , notably the lecture notes from the les houches session in 2004 by j. verbaarschot @xcite .
this brings me to the main goals of these lecture notes , to be addressed in the order of the subsequent sections .
what can we predict from rmt that can be compared to lattice qcd ?
in section [ opdirac ] we will start with the solution of the corresponding rmt using the theory of orthogonal polynomials .
this section is rather mathematical and contains no further physics input .
the second question , what the limit is , in which qcd reduces to rmt , is addressed in section [ symm ] . here
the global symmetries of qcd are explained , leading to chpt and eventually to a rmt description .
the more physics inclined reader may jump to this section first .
one of my personal motivations why to add another review on this topic is answered in section [ recent ] where the development in the last 10 years is reviewed , including some mathematical aspects .
this covers the dependence on finite lattice spacing @xmath19 and on chemical potential @xmath10 , for an earlier review on the latter see @xcite .
let us move to rmt , where i assume that the reader has already some working knowledge .
rmt is a much older topic , beginning in the late 1920 s with wishart in mathematical statistics and in the late 1950 s with wigner and dyson in nuclear physics .
the number of its applications is huge and still increasing , and we refer to @xcite for applications to quantum physics and to @xcite for a recent compilation containing physics ( including qcd in chapter 32 ) , mathematics and more . in these lecture notes
we will focus only on the theory of orthogonal polynomials .
for the standard gaussian ensembles they are covered in the last edition of mehta s classical book @xcite , for more recent monographs see @xcite .
the symmetry class we will be interested in for the application to qcd is the chiral gaussian unitary ensemble ( chgue ) , which is also known as wishart or laguerre unitary ensemble , and extensions thereof . in the simplest case of @xmath25 and
when all the abovementioned parameters from qcd are absent it is defined by the probability space of complex @xmath27 matrices @xmath28 .
all its matrix elements are independent and share the same complex normal distribution .
the mean density @xmath29 of the positive eigenvalues of @xmath30 , where @xmath31 denotes the hermitian conjugate , is well known in the limit @xmath32 .
after an appropriate rescaling it is given by the marchenko - pastur law depicted in figure [ fig : mpdensity ] , cf . the lectures by bouchaud where this was derived .
( 12.2,4.5 )   which is the limiting global spectral density of the chgue .
we also indicate the locations where different local statistics applies : they are given by i ) the bessel - kernel close to the hard edge , ii ) the sine - kernel in the bulk and iii ) the airy - kernel at the soft edge .
[ fig : mpdensity],title="fig:",width=264 ] ( 0.2,-0.2)@xmath33 ( -7.8,4.0)@xmath29 ( -1.9,-0.2)iii ) ( -6.6,-0.2)i ) ( -4.5,-0.2)ii )    the same law is obtained when taking real or quaternion valued matrix elements , constituting the chiral gaussian orthogonal and symplectic ensemble ( chgoe and chgse ) , respectively .
however , the local statistics of the 3 ensembles differs , and depends on the location in the spectrum we consider .
for finite @xmath34 the eigenvalue statistics follows a determinantal ( chgue ) or pfaffian point process ( chgoe , chgse ) , and the limiting kernels in the various locations of the spectrum shown in figure [ fig : mpdensity ] are given in terms of bessel , sine or airy functions .
in view of the relation ( [ bcrel ] ) we will be interested in the local statistics for the chgue at the origin also called hard edge .
the other two symmetry classes at the hard edge will also find applications in qft as pointed out in @xcite , cf .
@xcite , but will not be discussed here . in the limit where rmt applies the eigenvalues @xmath35 of the qcd dirac operator @xmath22
will be given by @xmath36 and thus come in pairs , where @xmath37 are the eigenvalues of @xmath30 ( the positive part @xmath38 are the singular values of @xmath28 ) .
after this change of variables the marchenko - pastur law becomes a semi - circle , which is why the global density of @xmath22 at the origin will become constant ( and not divergent as in figure [ fig : mpdensity ] ) . for finite @xmath34
the chgue can be solved in terms of the classical laguerre polynomials
. it will be the content of sections [ opdirac ] and [ recent ] to include more structure from qcd into this rmt , while maintaining its exact solvability .
in this section we go directly to the rmt of qcd in an eigenvalue representation that generalises the chgue , by including the mass terms of @xmath0 quark flavours . in subsection [ defs ]
we briefly recall the orthogonal polynomial ( op ) formalism and define the quantities of our interest .
we then attack the mathematical problem of computing these by deriving properties of op with a general weight , in particular including averages of characteristic polynomials in subsection [ genweight ] .
this allows us to determine all quantities explicitly in terms of standard laguerre polynomials in subsection [ alllaguerre ] , which facilitates the large-@xmath34 limit to be taken in subsection [ lim ] .
we begin by stating the rmt that we want to solve in its eigenvalue representation , defined by the following partition function @xmath39 it gives the normalisation constant of the unnormalised joint probability density function ( jpdf ) @xmath40 of all eigenvalues . here
we have included the dependence on the following parameters to be related to qcd later : @xmath41 called the topological index taking fixed values , the masses @xmath42 , @xmath43 , with @xmath0 counting the number of inserted characteristic polynomials .
in the quenched theory with @xmath25 their product is absent .
these parameters are all collected in the weight function @xmath44 and finally we also have introduced the vandermonde determinant @xmath45=\prod_{1\leq i < j\leq n}(x_j - x_i)\ . \label{delta}\ ] ] it depends on the set of all eigenvalues @xmath46 . the partition function eq .
( [ partitionfunct ] ) is relevant for qcd for @xmath47 , representing the chgue in the presence of @xmath0 characteristic polynomials ( mass terms ) . for completeness and in order to state some open problems later
we have also introduced the chgoe and chgse with @xmath48 and 4 , respectively .
all three rmt are also called wishart - laguerre ensembles and can be written in matrix representation , @xmath49\prod_{f=1}^{n_f}\det \left [ \begin{array}{cc } m_f 1_n&i w\\ iw^\dag & m_f 1_{n+\nu}\\ \end{array } \right ] \exp[-\mbox{tr}(ww^\dag)]\ . \label{matrixrep}\ ] ] here @xmath28 is an @xmath50 matrix taking values @xmath51 for @xmath52 respectively , and @xmath53 is its hermitian conjugate .
the integration @xmath54 $ ] is over the flat lebesgue measure of all independent matrix elements .
the identity matrices in front of the scalar mass variables @xmath55 of dimensions
@xmath34 and @xmath56 are denoted by @xmath57 and @xmath58 , respectively .
here we see that @xmath24 counts the number of zero eigenvalues of the wishart matrix @xmath59 , and the @xmath37 , @xmath60 denote its non - zero positive eigenvalues we always have an even number @xmath61 of doubly degenerate masses in eq .
( [ partitionfunct ] ) .
we will come back to this at the end of subsection [ alllaguerre ] . ] .
the change of variables from matrix elements to eigenvalues leads to the vandermonde determinant to the power @xmath52 times the @xmath24-dependent part as the jacobian , and we refer to standard text books on rmt for the derivation @xcite . in the lectures of j .-
bouchaud these ensembles were discussed for @xmath25 , in particular when @xmath24 is of the order of @xmath34 .
this matrix representation eq .
( [ matrixrep ] ) will play an important role later in section [ symm ] when relating the rmt to qcd based on symmetries . in the following
we will mainly focus on the chgue with @xmath47 which is the rmt relevant for qcd .
the values @xmath48 and 4 apply to other qfts and are very briefly discussed at the end of subsection [ alllaguerre ] .    what are the quantities we would like to calculate as functions of the masses @xmath55 and @xmath24 , apart from the normalising partition function ? the jpdf represents a determinantal ( or pfaffian ) point process introduced in the lectures by a. borodin .
therefore we can determine all @xmath62-point eigenvalue density correlations defined for all @xmath63 as @xmath64\ , \ \ \mbox{for } \ \ \beta=2\ . \label{rkb2}\end{aligned}\ ] ] in the second line valid for @xmath47 only they can be expressed in terms of the determinant of the kernel @xmath65 that does not contain the weights in our convention .
it contains op with respect to the weight eq .
( [ weight ] ) , @xmath66 for @xmath67 similar expressions exist in terms of pfaffian determinants of a matrix kernel of op with respect to a skew symmetric inner product replacing eq .
( [ op ] ) .
these are called skew op and we refer to @xcite for details .
we choose the op in eq .
( [ op ] ) to be monic , @xmath68 , with squared norms @xmath69 .
in particular eq .
( [ rk ] ) leads to the spectral density @xmath70 ( normalised to @xmath34 ) and to the determinantal expression for the jpdf , @xmath71 .
two further quantities can be defined that are of our interest .
the @xmath62-th gap probability defined for all @xmath63 reads @xmath72 it is the probability that @xmath62 eigenvalues are in @xmath73 $ ] and @xmath74 in @xmath75 .
for example @xmath76 gives the gap probability that the interval @xmath73 $ ] is empty .
if we order the eigenvalues @xmath77 the @xmath62-th eigenvalue distribution defined as @xmath78 gives the probability to find the @xmath62-th eigenvalue at @xmath79 , @xmath80 eigenvalues in @xmath73 $ ] and @xmath74 in @xmath75 .
these individual eigenvalue distributions are normalised @xmath81 , as can be easily checked ( see e.g. the appendix of @xcite )
. we will be mostly interested in the example of the distribution of the smallest non - zero eigenvalue , @xmath82 .
equations ( [ gapk ] ) and ( [ pk ] ) are not independent , in fact the latter follow by differentiation : @xmath83 which can be solved for @xmath84 as follows : @xmath85 the gap probabilities can also be expressed in terms of the kernel @xmath86 as a fredholm determinant , cf . the lectures by a. borodin and eqs .
( [ e0andreiev ] ) and ( [ expand ] ) , but we will use a different route to compute them explicitly .
our goal is now to determine all the quantities that we have introduced in this subsection , including the partition function . for @xmath25 this would be an easy task as the polynomials orthogonal with respect to @xmath87 $ ] are the well known laguerre polynomials .
it turns out that also for @xmath88 we can express all quantities in terms of laguerre polynomials , including the dependence on the parameter @xmath55 and @xmath24 that we are after .
this will be the goal of the following 2 subsections .
the first part of this subsection is standard material about op on the real line and can be found in any standard reference , e.g. in @xcite including also op in the complex plane needed later . for the second part with weights containing characteristic polynomials we refer to @xcite for more details .
let us consider a general measurable weight function on the positive real half line with the condition that all moments exist , @xmath89 .
the monic polynomials @xmath90 orthogonal with respect to @xmath91 can then be recursively constructed using gram - schmidt , expressing them as the ratio of two determinants : @xmath92\right)^{-1},\ ] ] see @xcite for details .
it obviously satisfies @xmath93 .
these polynomials obey a three - step recurrence relation , @xmath94 which can be see as follows . with the @xmath90 forming a complete set we can expand @xmath95 using that in the last integral @xmath96 is a polynomial of degree @xmath97 and the orthogonality eq .
( [ op ] ) it follows that @xmath98 for @xmath99 and @xmath100 . for the orthonormal polynomials @xmath101 the recursion eq .
( [ recurrence ] ) becomes more symmetric , @xmath102 and we obtain the christoffel - darboux formula for the kernel @xmath103 this simply follows by multiplying the sum by @xmath104 with @xmath105 , using the recursion and the fact that this is a telescopic sum . from lhpital s rule for the kernel at equal arguments we see that in the asymptotic limit @xmath32 we only need to evaluate the asymptotic of @xmath106 , @xmath107 , their derivatives , and the ratio of norms @xmath108 to determine the spectral density : @xmath109 let us give two examples .
for the gue with weight @xmath110 $ ] on @xmath111 instead of @xmath112 we have @xmath113 the probabilist s hermite polynomials .
second , for the laguerre weight @xmath87 $ ] , eq .
( [ weight ] ) with @xmath25 on @xmath112 , the monic op are given by the appropriately rescaled generalised laguerre polynomials @xmath114 : @xmath115 where @xmath24 can also be real here .
the asymptotic analysis of the kernel ( and density ) can now be made , and we will come back to this in subsection [ lim ] .    in the following
we will express the partition function , corresponding ops , kernel and gap probability @xmath76 in terms of expectation values of characteristic polynomials in a first step , and then reduce them to determinants of laguerre polynomials in a second step . in order to prepare this
the following exercise is very useful .
consider the 2 vandermode determinants eq .
( [ delta ] ) inside the partition function ( [ partitionfunct ] ) for @xmath47 . by adding columns inside the vandermonde determinant
we do not change its value . starting with the last column of highest degree we can thus express @xmath116 in terms of arbitrary monic polynomials
@xmath117 : @xmath118 \ .
\label{delta - id1}\ ] ] choosing those monic polynomials @xmath119 that precisely satisfy the orthogonality relation ( [ op ] ) we obtain for the partition function : @xmath120\right)^2\nonumber\\ & = & \sum_{\sigma,\sigma^\prime\in s_n}(-1)^{\sigma+\sigma^\prime } \prod_{l=1}^n\int_0^\infty dx_l w(x_l){p}_{\sigma(l)-1}(x_l){p}_{\sigma^\prime(l)-1}(x_l ) \ = \
\label{znorm}\end{aligned}\ ] ] in the second step we have laplace expanded both determinants into sums over permutations @xmath121 and used the orthogonality ( [ op ] ) .
the computation of the partition function thus amounts to determine the squared norms of the ops with respect to the weight eq .
( [ weight ] ) .
let us define the expectation value of an operator @xmath122 that only depends on the set of eigenvalues @xmath123 : @xmath124 then the following identity holds for the monic ops : @xmath125 it is given by the expectation value of a characteristic polynomial with respect to @xmath126 eigenvalues for an arbitrary degree , where @xmath127 .
it can also be expressed in terms of matrices of size @xmath126 , choosing @xmath128 $ ] .
this relation is known as the heine formula , dating back to the 19th century .
the proof uses the simple identity for the vandermonde determinant , @xmath129 \ .
\label{delta - id2}\ ] ] using this identity in the expectation value we arrive at the same situation as in eq .
( [ znorm ] ) , with one larger determinant : @xmath130 here the orthogonality of the polynomials enforces @xmath131 , and the remaining norms cancel due to eq .
( [ znorm ] ) . as a check also the leading coefficient of eq .
( [ opid ] ) can be compared , taking the asymptotic limit @xmath132 which leads to @xmath133 on both sides .
note that eq .
( [ opid ] ) provides an @xmath126-fold integral representation of the @xmath134 for an arbitrary weight function .
compared to the known single integral representation , e.g. for the hermite polynomials in the example in eq .
( [ hermite ] ) , we thus have the following duality relation @xmath135\ .
% = \langle\det[x - ih]^n\rangle_1 .\ ] ] here we take a single average over @xmath136 copies of characteristic polynomials .
the existence of such relations between 1- and @xmath136-fold averages of characteristic polynomials are known only for gaussian weights .
they can also be derived using the supersymmetric method as reviewed in chapter 7 of @xcite by thomas guhr , see also @xcite for further dualities using op techniques .    in the next step
we will express the kernel itself as an expectation value of two characteristic polynomials , following the work of p. zinn - justin @xcite .
namely it holds that @xmath137 the proof uses the same idea as before , now including each product into a different vandermonde determinant of size @xmath138 .
we have @xmath139 compared to eq .
( [ proofpn ] ) the orthogonality only implies @xmath140 , and its value can still run over all possible values from 1 to @xmath138 . in the last line
we give the expression for the kernel @xmath141 in terms of the christoffel - darboux formula ( [ cd ] ) , using monic polynomials instead .
once again the leading coefficients of left and right hand side can be see to agree @xmath142 in the limit @xmath143 .
an immediate question arises : what is the expectation value for more than two products of characteristic polynomials ? for that purpose let us introduce some more notation following baik , deift and strahov @xcite , to where we refer for more details .
denote by @xmath144}(x)$ ] the op with respect to weight @xmath145}(x)=(\prod_{j=1}^l(y_j - x))w(x)$ ] for @xmath146 , with @xmath147}(x)=p_n(x)$ ] .
the following relation holds which is called christoffel formula : @xmath148}(x)=\frac{1}{(x - y_1)\cdots(x - y_l ) } \left| \begin{array}{lll } p_n(y_1)&\ldots&p_{n+l}(y_1)\\ \vdots & & \vdots\\ p_n(y_l)&\ldots&p_{n+l}(y_l)\\ p_n(x)&\ldots&p_{n+l}(x)\\ \end{array } \right|\left(\det_{1\leq i , j\leq l}[p_{n+j-1}(y_i)]\right)^{-1}. \label{christoffel}\ ] ] as explained in @xcite this can be seen as follows .
it is clear that the determinant in the middle , which we call @xmath149}(x)$ ] , is a polynomial of degree @xmath150 in @xmath33 .
it has zeros at @xmath151 , and thus @xmath149}(x)/(x - y_1)\cdots(x - y_l)$ ] is a polynomial of degree @xmath136 in @xmath33 .
thus we have that @xmath152}(x)}{(x - y_1)\cdots(x - y_l)}w^{[l]}(x)=0\ , \mbox{for}\ \
j=0,1,\ldots , n-1\ , \ ] ] due to @xmath153 being op with respect to weight @xmath91 , after cancelling the factors from @xmath145}(x)$ ] .
the last factor in eq .
( [ christoffel ] ) ensures that @xmath144}(x)$ ] is monic , as can bee seen from taking the limit @xmath132 .
the repeated product of the christoffel formula ( [ christoffel ] ) leads to the following theorem ( see theorem 2.3 @xcite ) for the average of products of @xmath97 characteristic polynomials with @xmath146 , which is due to brzin and hikami @xcite : @xmath154 = \prod_{j=0}^lp_n^{[j]}(y_{j+1})= \left\langle   \prod_{i=1}^{n}\left(\prod_{j=1}^{l+1}(y_j - x_i ) \right)\right\rangle_n . \label{prodchar}\ ] ] for the last equality we have used eq .
( [ opid ] ) that leads to @xmath155}(y_{j+1})=\frac{\left\langle   \prod_{i=1}^{n}\left((y_{j+1}-x_i)\prod_{p=1}^{j}(y_p - x_i ) \right)\right\rangle_n}{\left\langle
\prod_{i=1}^{n}\left(\prod_{p=1}^{j}(y_p - x_i ) \right)\right\rangle_n}. \label{pnvev}\ ] ] in the product almost all factors cancel out . the case for a single characteristic polynomial with @xmath156 in eq .
( [ prodchar ] ) was already stated eq .
( [ opid ] ) .    in @xcite in theorem 3.2
a further identity was derived for the expectation value of the product of an even number of characteristic polynomials , expressing it through a determinant of kernels divided by two vandermonde determinants , @xmath157= \left\langle   \prod_{i=1}^{n}\left(\prod_{j=1}^{k}(\lambda_j - x_i ) ( \mu_j - x_i)\right)\right\rangle_n . \label{prod2char}\ ] ] the simplest example for this relation was derived in eq .
( [ kernelid ] ) . in @xcite
this set of identities was further generalised , expressing the expectation value of arbitrary products of characteristic polynomials by a determinant containing both kernels and polynomials divided by two vandermonde determinants in many different and equivalent ways , cf .
( [ cprodchar ] ) .
the simplest example with 3 products reads @xcite @xmath158    let us mention that in @xcite also expectation values of ratios of characteristic polynomials were determined .
they can be expressed in terms of the cauchy transforms @xmath159 of the monic polynomials @xmath90 , @xmath160 as well as through mixed chistoffel - darboux kernels containing both cauchy transforms and polynomials .
the simplest example is given by @xmath161 the fact that both sides agree @xmath162 for @xmath132 can be easily seen from the definition ( [ cauchy ] ) , expanding the geometric series and using the orthogonality of the polynomial @xmath163 to monic powers less than @xmath62 .    a second particularly important example is the expectation value of the ratio of two characteristic polynomials , due to the following relation to the resolvent or stieltjes transform @xmath164 : @xmath165 the resolvent can be used to obtain the spectral density through the following relation @xmath166 \ . \label{inversion}\ ] ] there are many examples where the op technique is not available , but where the expectation value of ratios of characteristic polynomials can be found by other means , e.g. by using supersymmetry , replicas or loop equations .
this then leads to an alternative way to determine the spectral density , or higher @xmath62-point correlation functions , by taking the average over @xmath62 resolvents and then taking the imaginary parts with respect to each of the @xmath62 arguments .
we refer to chapters 7 , 8 and 16 in @xcite for further details and references .      using the results from the previous subsection we are now ready to express the partition function , ops , kernel , and consequently all @xmath62-point correlation functions with @xmath0 characteristic polynomials , as well as the corresponding gap probabilities through expectation values of characteristic polynomials with respect to the quenched weight
@xmath87 $ ] with @xmath25 .
because the ops of this quenched weight are laguerre polynomials , everything will be finally expressed through these , which yields explicit and exact expressions for any finite @xmath34 and @xmath0 . in the next subsection [ lim ]
will use them to take the large-@xmath34 limit based on the known asymptotic of the laguerre polynomials .
we shall adopt the notation from the previous section , labelling all the above listed quantities by superscript @xmath167 $ ] compared to the quenched weight without superscript .
we have for the partition function eq .
( [ partitionfunct ] ) with masses @xmath168}(\{m\})&=&\left(\prod_{j=1}^n\int_0^\infty dx_j\ x_j^\nu e^{-x_j}\prod_{f=1}^{n_f}(x_j+m_f^2)\right)\delta_n(\{x\})^2 = z_n\left\langle \prod_{j=1}^n\prod_{f=1}^{n_f}(x_j+m_f^2)\right\rangle_n \nonumber\\ & = & n!(-1)^{n_f(n+(n_f-1)/2)}\prod_{j=1}^{n}\gamma(j+\nu)\prod_{j=1}^{n+n_f}\gamma(j)\   \frac{\det_{1\leq i , j \leq n_f } \left[l_{n+j-1}^\nu(-m_i^2)\right]}{\delta_{n_f}(\{-m^2\})}. \label{zmassive}\end{aligned}\ ] ] this result simply follows from choosing @xmath169 in eq .
( [ prodchar ] ) , together with eqs .
( [ znorm ] ) and ( [ laguerre ] ) for the quenched polynomials .
it is clear that if some of the masses become zero , say @xmath126 out of @xmath0 , from looking at eq .
( [ partitionfunct ] ) for @xmath47 without further calculation this leads to the shift @xmath170 in the remaining determinant of size @xmath171 .
this property is called flavour - topology duality . in the case
when all masses are degenerate , @xmath172 , lhopital s rule eventually leads to @xcite @xmath173}(m)=\frac{(-1)^{nn_f}}{\prod_{l=0}^{n_h-1}l!}\det_{0\leq i , j \leq n_f-1 } \left[(-1)^{n+i}(n+i)!(l_{n+i}^\nu(-m^2))^{(j)}\right]\ , \label{zmassivedeg}\ ] ] where the superscript @xmath174 denotes the @xmath175-th derivative with respect to the argument @xmath176 .    the unquenched op @xmath177}(x)$ ] follow directly from eq .
( [ pnvev ] ) at @xmath178 , with @xmath179 . for example using eq .
( [ laguerre ] ) the op for the weight @xmath180}(x)=(x+m^2)x^\nu\exp[-x]$ ] the corresponding op read @xcite @xmath181}(x)=\frac{h_nk_{n+1}(x ,- m^2)}{p_n(-m^2)}=\frac{(-1)^{n+1}(n+1)!}{(x+m^2)l_n^\nu(-m^2)}(l_{n+1}^\nu(x)l_{n}^\nu(-m^2)-l_{n+1}^\nu(-m^2)l_{n}^\nu(x))\ . \label{pmassive}\ ] ] however , it is much simpler to directly use the kernel from eq .
( [ kernelid ] ) , rather than expressing it through the polynomials @xmath182}(x)$ ] using the christoffel - darboux identity ( [ cd ] ) .
indeed the kernel is given by eq .
( [ kernelid ] ) as @xmath183}(x , y)&= & \frac{1}{h_{n-1}^{[n_f]}}\left\langle \prod_{i=1}^{n-1}(x - x_i)(y - x_i ) \right\rangle_{n-1}^{[n_f ] } = \frac{1}{h_{n-1}}\frac { \left\langle\prod_{i=1}^{n-1}\left((x - x_i)(y - x_i)\prod_{f=1}^{n_f}(m_f^2+x_i)\right ) \right\rangle_{n-1 } } { \left\langle\prod_{i=1}^{n}\prod_{f=1}^{n_f}(m_f^2+x_i ) \right\rangle_{n } } \nonumber\\ & = & \frac{(-1)^{n_f-1}(n+n_f)!}{\gamma(n+\nu)(y - x)\prod_{f=1}^{n_f}(y+m_f^2)(x+m_f^2 ) } \frac { \left| \begin{array}{lll } l_{n-1}^\nu(-m_1 ^ 2)&\ldots&l_{n+n_f}^\nu(-m_1 ^ 2)\\ \vdots & & \vdots\\ l_{n-1}^\nu(-m_{n_f}^2)&\ldots&l_{n+n_f}^\nu(-m_{n_f}^2)\\ l_{n-1}^\nu(x)&\ldots&l_{n+n_f}^\nu(x)\\ l_{n-1}^\nu(y)&\ldots&l_{n+n_f}^\nu(y)\\ \end{array } \right| } { \det_{1\leq f , g\leq n_f}[l_{n+g-1}^\nu(-m_f^2)]}. \label{kernelmassive}\end{aligned}\ ] ] in the first step after inserting eq .
( [ kernelid ] ) we have included the mass dependent norm @xmath184}$ ] into the massive partition function in the denominator .
it is given by the product of its norms times @xmath185 .
thus we have increasing the average from @xmath186 to @xmath34 . in the final result in the second line we have already taken out all signs and factorials from the two determinants and cancelled them partly .
also the ratio of the two vandermonde determinants , one of which containing the arguments @xmath33 and @xmath187 has been simplified using the definition ( [ delta ] ) .
the kernel eq .
( [ kernelmassive ] ) determines all @xmath62-point eigenvalue correlation functions from eq .
( [ rkb2 ] ) including their mass dependence .
in particular for @xmath188 we obtain the spectral density with @xmath0 masses : @xmath189}(x)=x^\nu\exp[-x]\prod_{f=1}^{n_f}(x+m_f^2)\ k_n^{[n_f]}(x , x)\ , \label{r1massive}\ ] ] after applying lhpital s rule once .
this leads to @xmath190 in the last row in the numerator .
note that the product from the weight in front of the kernel cancels part of its denominator .
let us now turn to the gap probabilities and distribution of smallest eigenvalues , starting with @xmath76 from eq .
( [ gapk ] ) . the following andriev integral identity ( cf .
borodin s lectures ) can be used to get a first expression : @xmath191\det_{1\leq i , j \leq n}[\psi_i(x_j ) ] = n!\det_{1\leq i , j \leq n}\left[\int dx\phi_i(x)\psi_j(x)\right]\ .
\label{andreiev}\ ] ] the only condition to hold is that all integrals of the functions @xmath192 and @xmath193 exist .
the proof merely uses the laplace expansion of the left hand side : @xmath194 and integration over common arguments of @xmath195 and @xmath196 .
this implies @xmath197 , which is yet another permutation .
following eq .
( [ gapk ] ) we can thus write @xmath198}(s ) = \frac{1}{n!%\prod_{l=0}^{n-1}h_l } \prod_{p=1}^n\int_s^\infty dx_p   \det_{1\leq i , j \leq n}\left[\sqrt{w^{[n_f]}(x_j)}\hat{p}_{i-1}^{[n_f]}(x_j)\right]^2 = \det_{1\leq i , j \leq n}\left[\delta_{ij}-\int_0^s dx w^{[n_f]}(x)\hat{p}_{i-1}^{[n_f]}(x)\hat{p}_{j-1}^{[n_f]}(x)\right]\ ! , \label{e0andreiev}\ ] ] which can be interpreted as a fredholm determinant .
in the first step we have replaced the vandermonde determinants by determinants of monic polynomials , and then included the square root of the weight functions and of the norms stemming from the normalising partition function into the determinants , in order to make the polynomials orthonormal , @xmath199}(x_j)$ ] . in the next step
we have used the andrief formula and the orthonormality on @xmath200 .
knowing both the monic polynomials from eq .
( [ pnvev ] ) and their squared norms from @xmath201}=z_j^{[n_f]}/(j z_{j-1}^{[n_f]})$ ] in terms of expectation values of laguerre polynomials , eq .
( [ e0andreiev ] ) is a valid expression for the gap @xmath76 .
however , we will use a more direct expression to be derived below . as a further remark we note that replacing @xmath202 in the first equality in ( [ e0andreiev ] ) we obtain an equivalent expansion of the fredholm determinant in the second equation : @xmath198}(s ) = 1+\sum_{k=1}^n\frac{(-1)^k}{k!}\int_0^s dx_1\ldots \int_0^s
dx_k r_k^{[n_f]}(x_1\ldots , x_k ) .
\label{expand}\ ] ] in order to express the gap probability directly as an expectation value let us go back to its definition ( [ gapk ] ) , which we state for general @xmath63 : @xmath203}(s)&=&\frac{1}{z_n^{(\beta , n_f,\nu)}}\left(\prod_{j=1}^n\int_s^\infty dx_j\ x_j^{\frac\beta2(\nu+1)-1}e^{-x_j}\prod_{f=1}^{n_f}(x_j+m_f^2)\right ) |\delta_n(\{x\})|^\beta\nonumber\\ & = & \frac{e^{-n s}}{z_n^{(\beta , n_f,\nu)}}\left(\prod_{j=1}^n\int_0^\infty dy_j\ y_j^0 ( y_j+s)^{\frac\beta2(\nu+1)-1}e^{-y_j}\prod_{f=1}^{n_f}(y_j+s+m_f^2)\right ) |\delta_n(\{y\})|^\beta .
\label{gap - shift}\end{aligned}\ ] ] here we have shifted all integration domains , substituting @xmath204 , @xmath60 .
the shift can be simply worked out , in particular it leaves the vandermonde determinant invariant , as a consequence the numerator is again given by a partition of @xmath205 masses with values @xmath206 for the first @xmath0 , and @xmath207 for the remaining ones , with an effective topological charge @xmath208 parameter satisfying @xmath209 .
this idea has been used in several papers @xcite to give closed form expressions as an alternative to the fredholm determinant above .
let us first consider @xmath47 , with @xmath210 .
consequently eq .
( [ gap - shift ] ) can be written as @xmath211}(s)= e^{-ns } \frac{z_n^{(2,n_f+\nu,0)}}{z_n^{(2,n_f,\nu ) } } \left\langle \prod_{i=1}^n\left ( ( s+y_j)^\nu\prod_{f=1}^{n_f}(m_f^2+s+y_j)\right)\right\rangle_{n,\,\nu_{ef\!f}=0}. \label{gap - vevb2}\ ] ] the simplest examples with @xmath25 and @xmath212 thus read @xmath213 in the second line we have replaced the expectation value of a single characteristic polynomial eq .
( [ opid ] ) by the laguerre polynomial eq .
( [ laguerre ] ) and fixed the normalisation by the requirement @xmath214 .
the quenched gap probability ( @xmath25 ) with arbitrary @xmath24 can be computed from the massive partition functions eq .
( [ zmassivedeg ] ) at complete degeneracy .
the smallest eigenvalue distribution @xmath82 easily follows for the examples we have just given by differentiation of the gap probability , see eq .
( [ p1 ] ) . for higher gap probabilities and for computing @xmath84
directly the same trick from eq .
( [ gap - shift ] ) can be used .
first consider only the integrals @xmath215 in the definitions ( [ gapk ] ) and ( [ pk ] ) .
then do the shift @xmath204 , and perform the remaining integrals @xmath216 over the obtained expectation value of characteristic polynomials in the end .
we refer to @xcite for more details .
let us now briefly comment on @xmath67 , mainly because of recent developments and open questions in these symmetry classes .
we begin with @xmath48 . here
the effective topological charge is @xmath217 to satisfy @xmath218 in eq .
( [ gap - shift ] ) .
the additional mass terms originating from the shift by @xmath219 in eq .
( [ gap - shift ] ) appear with multiplicity @xmath220 : @xmath221}(s)= e^{-ns } \frac{z_n^{(1,n_f+(\nu-1)/2,1)}}{z_n^{(1,n_f,\nu ) } } \left\langle \prod_{i=1}^n\left ( ( s+y_j)^{(\nu-1)/2}\prod_{f=1}^{n_f}(m_f^2+s+y_j)\right)\right\rangle_{n,\,\nu_{ef\!f}=1}. \label{gap - vevb1}\ ] ] that is for @xmath222 odd , @xmath223 we have @xmath224 extra masses @xmath207 compared to the @xmath0 shifted masses . here
the simplest case is @xmath225 with @xmath25 , for which we obtain @xmath226 for arbitrary @xmath34 .
it agrees with the result for @xmath47 at @xmath227 , eq .
( [ e0nu0 ] ) . for general odd
@xmath24 the gap probabilities have been computed and we refer to @xcite for results . however , for @xmath228 even we obtain a half integer number @xmath229 of additional mass terms .
this leads to the question of calculating expectation values including products of square roots of characteristic polynomials in eq .
( [ gap - vevb1 ] ) . in @xcite
the expectation values needed to determine @xmath76 for arbitrary even @xmath228 and general @xmath0 have been computed for the chgoe , answering at least the question for the gap probability .
the cases @xmath227 @xcite and @xmath230 @xcite were previously known from different considerations .
namely in @xcite a recursive construction was made for the smallest eigenvalue for @xmath48 valid for all @xmath24 .
however , in this approach the pfaffian structure appearing for @xmath231 is not at all apparent .    independently in @xcite the question about expectation values of square roots of determinants has been asked and answered for special cases for the goe in the context of scattering in chaotic quantum systems .
it is an open problem if the most general expectation value of such products ( or ratios ) of square roots has a determinantal or pfaffian structure as in eq .
( [ prodchar ] ) .
let us turn to @xmath232 .
the effective topological charge we obtain is @xmath233 to satisfy @xmath234 in eq .
( [ gap - shift ] ) .
the fact that it is not integer does not pose a problem as the index of the generalised laguerre polynomials eq .
( [ laguerre ] ) in terms of which we have expressed expectation values of characteristic polynomials can be chosen accordingly .
the more sever problem is the following . due to
kramer s degeneracy the eigenvalues @xmath37 of a self - dual quaternion valued @xmath27 matrix @xmath235 always come in pairs , implying @xmath236=\prod_{j=1}^n(x - x_j)^2 $ ] .
consequently the mass terms generated from eq .
( [ matrixrep ] ) always occur with an _ even _ power @xmath0 of two - fold degenerate masses in eq .
( [ partitionfunct ] ) .
in contrast in eq .
( [ gap - shift ] ) we always need to evaluate an odd number of powers @xmath237 .
thus as for @xmath48 with even @xmath24 we need to evaluate square roots of determinants in order to compute the gap probability for @xmath232 .
this is an open problem so far and only the gap probability ( and smallest eigenvalue distribution ) with @xmath227 is known explicitly @xcite .
a taylor series expansion exists though for @xmath238 , following from group integrals of kaneko - type @xcite .      after having exhaustively presented the solution of the eigenvalues model eq .
( [ partitionfunct ] ) and its correlation functions for a finite number of eigenvalues @xmath34 at arbitrary @xmath0 for @xmath47 , we will now turn to the large-@xmath34 limit . as it is true in general in rmt one
has to distinguish between different large-@xmath34 scaling regimes .
the global spectral statistics , which includes the global macroscopic density given by the semi - circle law for the gue , or the marchenko - pastur law for the chgue , see figure [ fig : mpdensity ] , is concerned with correlations between eigenvalues that have many other eigenvalues ( in fact a finite fraction of all ) in between them . for a discussion of this large-@xmath34 limit alternative techniques
are available such as loop - equations , and we refer to @xcite for a standard reference where all correlation functions including subleading contributions are computed recursively for @xmath47 . in this limit the fluctuations of the eigenvalues on a local scale ( of a few eigenvalues ) are averaged out , and typically expectation values factorise .
in contrast when magnifying the fluctuations among eigenvalues at a distance of @xmath239 one speaks of microscopic limits .
the value of @xmath240 and the form of the limiting kernel depend on the location in the spectrum , see figure [ fig : mpdensity ] .
at the so - called soft edge of the marchenko - pastur law one finds the airy - kernel , whereas in the bulk of the spectrum the sine - kernel is found , see e.g. @xcite . here
we will be interested in the eigenvalues in the vicinity of the origin presenting a hard edge with @xmath241 and the limiting kernel to be computed below is the bessel - kernel .
why this limit is relevant for the application of rmt to qcd has been indicated already in the introduction after eq .
( [ bcrel ] ) .    as a further consequence of this comparison to qcd
we will change variables from wishart eigenvalues @xmath37 of @xmath30 to dirac operator eigenvalues @xmath36 ( more details follow in section [ symm ] ) .
looking at the partition function eq .
( [ partitionfunct ] ) the dirac eigenvalues and masses have to be rescaled with the same power in @xmath34 .
otherwise we would immediately loose the dependence on the masses . on the rmt side
the hard - edge scaling limit zooming into the vicinity of the origin is conventionally defined including factors by taking @xmath32 and @xmath242 , such that the following product remains finite - independent here in contrast to other conventions . ]
@xmath243 in the last subsection we have eventually expressed all quantities of interest in terms of the laguerre polynomials of the quenched weight . therefore the only piece of information we need to take the asymptotic limit of all these quantities is the well known limit @xcite : @xmath244 consequently all laguerre polynomials of positive argument turn into @xmath245-bessel functions .
those polynomials containing negative arguments will turn into @xmath246-bessel functions , due to the relation @xmath247 for integer @xmath24 . note that both types of bessel functions also appear together ,
see e.g. in eqs .
( [ pmassive ] ) and ( [ kernelmassive ] ) .
the asymptotics of the laguerre op eq .
( [ laguerrelim ] ) in fact holds for op with respect to much more general weight functions .
this phenomenon is called universality . after first results in @xcite more sophisticated rigorous mathematical methods including the riemann - hilbert approach
were developed , see chapter 6 by a. kuijlaars in @xcite for a detailed discussion and for references .    in the following
we will give a few examples for the hard edge limit , starting with the partition function .
looking at eq .
( [ zmassive ] ) it is clear that we have to normalise it differently , such that the limit exists and gives a function of the limiting masses from eq .
( [ rmtmass ] ) . for @xmath248
this is very easy to do by dividing out the quenched partition function , and we define : @xmath249}(m)=m^{\nu}n^{-\frac{\nu}{2}}\frac{z_n^{[n_f=1]}(m)}{z_n n!}= m^{\nu}n^{-\frac{\nu}{2 } } l_n^\nu(-m^2)\ , \label{zrmtredef}\ ] ] leading to the limit @xmath250}(m)={\cal z}^{[n_f=1]}(\tilde{m})= i_\nu(\tilde{m})\ . \label{zrmtnf1}\ ] ] note that our conventions are chosen such that for small mass the limiting partition function behaves as @xmath251}(\tilde{m}\approx0)\sim \tilde{m}^\nu$ ] , rather than unity .
the reason for this choice will also become transparent in the next section . for general @xmath0 we do not spell out the normalisation constant explicitly and directly give the limiting result for non - degenerate masses @xmath252}(\{m\})={\cal z}^{[n_f]}(\{\tilde{m}\})=   2^{\frac{n_f(n_f-1)}{2}}\prod_{j=0}^{n_f-1}\gamma(j+1)\   \frac{\det_{1\leq f ,
g\leq n_f}\left[\tilde{m}_f^{g-1}i_{\nu+g-1}(\tilde{m}_f)\right]}{\delta_{n_f}(\{\tilde{m}^2\})}\ , \label{zrmtnf}\ ] ] where we adopted the normalisation from @xcite . in the completely degenerate limit this leads to : @xmath253}\right|_{m_f = m } { \cal z}^{[n_f]}(\tilde{m } ) = { \det_{1\leq f , g\leq n_f}\left[i_{\nu+g - f}(\tilde{m})\right]}\ .
\label{zrmtnfdeg}\ ] ] the direct limit @xmath32 of eq .
( [ zmassivedeg ] ) ( or lhopital s ) rule from eq .
( [ zrmtnf ] ) ) leading to a determinant of derivatives of @xmath246-bessel functions can be brought to this form of a toeplitz determinant using identities for bessel functions .    as the next step we will directly jump to the limiting kernel , the reason being two - fold .
first , it is shorter to directly construct all correlation function form the kernel , eq .
( [ kernelmassive ] ) , rather than constructing the kernel through the massive op s first .
second and more importantly , the limit of the individual op s does not necessarily exist . while here in the hard edge limit this is not the case , the well - known sine- and cosine - asymptotic of the hermite polynomials in the bulk of the spectrum is only achieved after multiplying them with the square root of the weight function .
we obtain for the limiting kernel from eq .
( [ kernelmassive ] ) , after appropriately normalising and changing to dirac eigenvalues , @xmath254}(\tilde{y}_1,\tilde{y}_2)&\sim & \lim_{\stackrel{\mbox{\small $ n\to\infty$}}{\mbox{\small $ x_1,x_2,m_f\to0 $ } } } ( w(x_1)w(x_2))^{\frac12}k_n^{[n_f]}(x_1,x_2)\nonumber\\ & = & \frac{1}{(\tilde{y}_1 ^ 2-\tilde{y}_2 ^ 2)\prod_{f=1}^{n_f}\sqrt{(\tilde{y}_1 ^ 2+\tilde{m}_f^2)(\tilde{y}_2 ^ 2+\tilde{m}_f^2 ) } } \frac { \left| \begin{array}{lll } i_\nu(\tilde{m}_1)&\ldots&\tilde{m}_1^{n_f+1}i_{n_f+\nu+1}(\tilde{m}_1)\\ \vdots & & \vdots\\ i_\nu(\tilde{m}_{n_f})&\ldots&\tilde{m}_{n_f}^{n_f+1}i_{n_f+\nu+1}(\tilde{m}_{n_f})\\ j_\nu(\tilde{y}_1)&\ldots&(-\tilde{y}_1)^{n_f+1}j_{\nu+n_f+1}(\tilde{y}_1)\\ j_\nu(\tilde{y}_2)&\ldots&(-\tilde{y}_2)^{n_f+1}j_{\nu+n_f+1}(\tilde{y}_2)\\ \end{array } \right| } { \det_{1\leq f , g\leq n_f}[m_f^{g-1}i_{\nu+g-1}(\tilde{m}_f)]}.\ \ \ \   \label{kernelmassivelim}\end{aligned}\ ] ] at @xmath25 it reduces to the standard bessel - kernel : @xmath255 here we have given two equivalent forms of the same kernel using an identity @xcite for @xmath245-bessel functions , @xmath256 .
following eq .
( [ rkb2 ] ) we have already included the weight function into the limiting kernel above . while the exponential part cancels out due to @xmath257\to1 $ ] the remaining @xmath24-dependent part
is contributing .
we can thus immediately write down an example , the quenched microscopic density @xmath258 we have omitted here the @xmath24 delta - functions of the exact zero - eigenvalues .
they will enter the discussion later in subsection [ wrmt ] where they can be seen .
the density is plotted in figure [ fig : besseldensity ] for @xmath227 and 1 .
the next step is to obtain the limiting distribution of the smallest eigenvalue . because we change from wishart to dirac eigenvalues @xmath259 , we obtain the following scaling limit for the examples of the gap probability ( [ e0nu0 ] ) and ( [ e0nu1 ] ) at @xmath25 , using the same scaling for the dirac eigenvalues as in eq .
( [ rmtmass ] ) : @xmath260 note that the extra power of @xmath34 from the shift ( [ gap - shift ] ) leads to the fact that here the exponential is not vanishing .
the liming quenched gap probability for general @xmath24 follows from the completely degenerate massive partition function ( [ zrmtnfdeg ] ) : @xmath261\ .
\label{equenched}\ ] ] the limiting distribution of the corresponding smallest dirac eigenvalue follows according to eq .
( [ p1 ] ) , after changing to dirac eigenvalues : @xmath262 where for the last step we have used an identity for bessel functions .
the last two formulas are compared to the quenched bessel density eq .
( [ besseldensity ] ) in figure [ fig : besseldensity ] .
it turns out that for general @xmath24 the derivative of the determinant in eq .
( [ equenched ] ) can be written again as a determinant , following the ideas sketched after eq .
( [ e0nu1 ] ) that the distribution of the smallest eigenvalue can itself be expressed in term of an expectation value of characteristic polynomials , @xmath263\ , \label{p1quenched}\ ] ] see @xcite for a detailed derivation , including masses .
( 12.2,5.5 ) ) and the distribution of the smallest eigenvalue for @xmath227 ( left ) and @xmath264 ( right ) using eqs .
( [ p1nu0lim ] ) and ( [ p1nu1lim ] ) .
it can be nicely seen how the distribution of the smallest eigenvalues follows the microscopic density for small @xmath265 . the local maxima of the density further to the right mark the locations of the second , third etc .
eigenvalues , cf .
the limiting value @xmath266 that equals the macroscopic density at the origin is marked as a horizontal line .
[ fig : besseldensity],title="fig:",width=264 ] ) and the distribution of the smallest eigenvalue for @xmath227 ( left ) and @xmath264 ( right ) using eqs .
( [ p1nu0lim ] ) and ( [ p1nu1lim ] ) .
it can be nicely seen how the distribution of the smallest eigenvalues follows the microscopic density for small @xmath265 . the local maxima of the density further to the right mark the locations of the second , third etc .
eigenvalues , cf .
the limiting value @xmath266 that equals the macroscopic density at the origin is marked as a horizontal line .
[ fig : besseldensity],title="fig:",width=264 ] ( 0.1,0.2)@xmath265 ( -7.1,0.2)@xmath265 ( -2.5,3.5)@xmath267 ( -10.8,3.5)@xmath267 ( -3.5,1.4)@xmath268 ( -11.5,1.4)@xmath268    it should be mentioned that alternative representations exist for the gap probability eq .
( [ p1quenched ] ) that are based on fredholm determinant analysis @xcite .
they are given in terms of an exponential of an integral of the solution of a painlev v equation , which is valid for any real @xmath269 .
for integer values of @xmath24 this expression reduces to the simpler forms that we have given here .
furthermore , in @xcite representations in terms of hypergeometric functions with matrix argument were given that are valid for all @xmath270 .
in @xcite the equivalence to the above expressions for @xmath47 was proved .
the distributions that we have computed so far , as for example the density and smallest eigenvalue in figure [ fig : besseldensity ] have been extensively compared to numerical solutions from lattice qcd .
for example the quenched densitiy in figure [ fig : besseldensity ] was matched with qcd in @xcite , whereas the quenched distribution of the smallest eigenvalue also given in figure [ fig : besseldensity ] was compared in @xcite to qcd for different values of @xmath24 ( and also to other qfts corresponding to @xmath67 ) . a matching to unquenched results for the smallest eigenvalues
can be found in @xcite , and we refer to @xcite for a more exhaustive discussion of the lattice results .
apart from giving non - specialists some idea about the theory of quantum chromodynamics ( qcd ) we will have to answer three major questions in this section :    * how is the dirac operator defined the spectrum of which can be described by rmt ? * what are the global symmetries of qcd that determine the rmt to be used ? * what is the limit that leads from qcd to this rmt ?
let us give first short answers that also indicate how this section will be organised .
the eigenvalues of the qcd dirac operator will be partly described by rmt , actually only its smallest eigenvalues , and we will introduce this operator first in subsection [ qcd ] .
the anti - hermiticity properties of the dirac operator and the concept of chiral symmetry will allow us to introduce some formal aspects and the definition of the qcd partition function .
the breaking of chiral symmetry in subsection [ xsb ] then leads us to the first approximation of qcd , to chiral perturbation theory that only contains the lightest excitations , the goldstone bosons . in a last step in subsection [ echpt ] this interacting non - linear , effective qft
is simplified in the epsilon - regime , that agrees with rmt to leading order .    in what follows it
is important that we will talk about space - time coordinates @xmath271 with components @xmath272 , @xmath273 equipped with the euclidean metric @xmath274 given by the kronecker symbol .
this so - called wick rotation from minkowski to euclidean metric is necessary for two reasons : only in this setting does the dirac operator have defined anti - hermiticity properties , and second only for qcd on a euclidean space - time lattice numerical solutions of qcd can be easily performed to which we can compare our rmt .
we use einstein s summation conventions meaning that greek indices appearing twice are automatically summed over from 1 to 4 .
the qcd dirac operator is a linear , matrix valued differential operator defined as follows : @xmath275 where @xmath276 .
the @xmath277 are the dirac - gamma matrices in euclidean space - time , satisfying the clifford algebra @xmath278 the curly brackets denote the anti - commutator . a standard representation as @xmath279 matrices is given as follows in terms of the pauli matrices @xmath280 , @xmath281 : @xmath282 all matrices are hermitian , @xmath283 , including @xmath284 .
this is only true in the euclidean setting
. obviously one can construct projection operators as @xmath285 the real parameter @xmath6 in the dirac operator is the coupling constant of the strong interactions , thus @xmath286 corresponds to the free dirac operator without interactions .
it acts on 4-vectors @xmath287 called spinors in minkowski space , due to their transformation properties there which we do nt need to specify here .
the real vector potential @xmath288 is an element of this lie group for the three colours ( the subscript @xmath289 stands for colour ) .
it is a scalar with respect to the @xmath277 s , just as the partial derivatives in eq .
( [ dirac ] ) .
the @xmath287 come in 3 copies with colours blue @xmath290 , green @xmath291 and red @xmath292 , and so the interacting dirac operator acts on the full vector @xmath293 . for each quark
flavour up , down etc .
we have such a vector @xmath294 with @xmath295 . in the standard model of elementary particle physics
we have @xmath3 , but often we will keep @xmath0 fixed as a free parameter
. we could have written the dirac operator with tensor product notation to underline these structures , but the precise @xmath296 structure will not be crucial in the following .
we are now ready to state some important global symmetries of the dirac operator .
first , it follows from the algebra ( [ clifford ] ) that @xmath297 the dirac operator is block off - diagonal also called chiral , where the symbol @xmath298 still denotes a differential operator that depends on @xmath33 and @xmath299 .
suppose we can find the eigenfunctions @xmath300 of the dirac operator labelled by @xmath62 : @xmath301 in a suitably regularised setting ( e.g. a finite box ) . because of the anti - hermiticity the eigenvalues @xmath302 are purely imaginary .
furthermore , because of eq .
( [ dchiral ] ) multiplying equation ( [ eigenfunct ] ) by @xmath284 leads to different eigenfunctions @xmath303 with eigenvalues @xmath304 , provided that @xmath305 .
thus the non - zero eigenvalues of the dirac operator come in pairs @xmath306 .
the euclidean qcd partition function or path integral as it is called in qft can be formally written down as follows : @xmath307\int[d\psi]\exp\left [   -\int d^4x \sum_{q=1}^{n_f}\overline{\psi}^q(x)(d\!\!\!\!/+m_q)\psi^q(x)-\int d^4x \frac12\mbox{tr}(f_{\mu\nu}f_{\mu\nu } ) \right]\nonumber\\ & = & \int[da_\mu]\prod_{q=1}^{n_f}\det[d\!\!\!\!/+m_q]\exp\left [ -\int d^4x \frac12\mbox{tr}(f_{\mu\nu}f_{\mu\nu})\right]\ . \label{zqcd}\end{aligned}\ ] ] here @xmath308 is the dirac conjugate and the @xmath296 field strength tensor is defined by the commutator @xmath309=-ig_sf_{\mu\nu}$ ] .
the integration over @xmath310 and @xmath311 are formal ( path ) integrals that we will not specify in the following , see standard textbooks on qcd for a discussion as @xcite .
in the second equation however we have formally integrated out the quarks due to the following observation . for a complex @xmath34-vector @xmath312 and @xmath27
matrix @xmath313 we have @xmath314\sim \frac{1}{\det[b]}\ , \ \
\int d^{2n}\psi\exp[-\psi^\dag b
\psi]\sim { \det[b]}\ , \label{detintegrals}\ ] ] whereas in the second identity we have integrated over a complex vector @xmath315 of anticommuting variables .
integration and differentiation over such fermionic variables can be defined in a precise way , and we refer to chapter 7 in @xcite for details .
the reason for fermions to be represented by anti - commuting variables is that in the canonical quantisation of free dirac fields one has to use anti - commutators for these fields .
comparing eq .
( [ dchiral ] ) and the second line of eq .
( [ zqcd ] ) with eq .
( [ matrixrep ] ) we already get a first idea about the random matrix approximation of qcd , as it was suggested initially in @xcite : the dirac operator from eq .
( [ dchiral ] ) is replaced by a constant matrix with the same block structure , @xmath316 , and the average over the fields strength tensor is replaced by a gaussian average over @xmath28 .
we further learn in this approximation that the eigenvalues @xmath317 of the gaussian random matrix @xmath30 correspond to the squared eigenvalues of the dirac operator rotated to the real line , and that the number of eigenvalues @xmath34 is proportional to the dimension of the qcd dirac operator truncated in this way .
however , as we will explain in the next subsection the precise form of the limit from qcd to rmt was better understood later , representing a controlled approximation .    furthermore , in the case when all masses vanish @xmath318 the action in the exponent in the first line of eq .
( [ zqcd ] ) has an additional symmetry . defining the vector @xmath319 as well as the projections @xmath320 onto left ( l ) and right ( r )
chirality we can write the fermionic part of the action as @xmath321 it is invariant under the global rotations @xmath322 with @xmath323 and likewise for @xmath324 .
in contrast , even when all masses are equal , @xmath325 , the mass term @xmath326 is only invariant under the diagonal transformation with @xmath327 . in qcd the resulting global symmetry
is @xmath328 , with the @xmath329-factors called @xmath330 for vector and @xmath331 for axial - vector split off .
the latter is broken through quantum effects also called an anomaly , whereas the former remains unbroken through the mass term .
this leads to the explicit chiral symmetry breaking pattern in qcd : @xmath332      in this subsection we will describe the first step of simplifying qcd to what is called its low - energy effective theory .
we begin with the following observation .
it is found that the qcd vacuum @xmath333 , the ground state with the lowest energy , breaks chiral symmetry _
spontaneously _ , in precisely the same way as the theory with equal masses in eq .
( [ xsbpattern ] ) breaks it explicitly : @xmath334 the parameter @xmath8 is called the chiral condensate and it is the order parameter for the spontaneous breaking of the chiral symmetry ( [ xsbpattern ] ) .
it turns out that for qcd the low energy phase with non - vanishing @xmath11 coincides with the phase or different representations this is not necessarily the case . ] where the quarks and the gluons , the carriers of the strong interaction represented by @xmath299 , are confined to colourless objects , see the phase diagram in figure [ fig : phasediag ] . in the phase with @xmath11 we need lattice qcd or other effective field theory description ( such as rmt ) because a perturbative expansion in the coupling constant @xmath6 will not work here , because the coupling is large .
an important consequence was drawn by banks and casher @xcite , relating the global continuum density @xmath335 of the dirac operator eigenvalues in eq .
( [ eigenfunct ] ) at the origin , suitably regularised in a finite volume @xmath18 , to @xmath8 : @xmath336 here we added the second order term to eq .
( [ bcrel ] ) derived by stern and smilga @xcite , containing the pion decay constant @xmath337 as an additional parameter .
the first term leads to the following interpretation : at the origin the global density is constant , and thus the average distance between eigenvalues there is @xmath338 .
this is very different from free particles in a box where @xmath339 with @xmath340 in 4 euclidean dimensions .
one can thus say that the smallest dirac operator eigenvalues make an important contribution to build up the chiral condensate @xmath11 , by piling up very closely spaced at the origin .
we are thus lead to define rescaled , dimensionless eigenvalues and masses , as they appear on the same footing inside the determinant in eq .
( [ zqcd ] ) , together with a rescaled microscopic density of the dirac operator @xmath341 : @xmath342 this is called the microscopic limit of qcd .
it has to be compared e.g. to eq .
( [ besseldensity ] ) for @xmath25 , after identifying rescaled rmt eigenvalues @xmath343 and masses @xmath344 from eq .
( [ rmtmass ] ) .
the breaking of a global continuous symmetry in any theory leads to goldstone bosons .
this can be visualised by analogy in figure [ fig : goldstone ] , showing the classical potential of a particle before ( left ) and after symmetry breaking ( right ) . in the left picture with a convex potential any excitation to make a particle move in the potential costs energy .
after symmetry breaking the curvature at the origin has changed . consequently there
now exists the possibility to excite a mode in angular direction along the valley of the potential , without costing potential energy .
modes that do need energy to be excited still exist in radial direction .
( 12.2,5.5 )   before symmetry breaking ( @xmath345 , left ) and after breaking ( @xmath346 , right ) , obtained simply by switching the sign of the quadratic term .
[ fig : goldstone],title="fig:",width=264 ]   before symmetry breaking ( @xmath345 , left ) and after breaking ( @xmath346 , right ) , obtained simply by switching the sign of the quadratic term .
[ fig : goldstone],title="fig:",width=264 ] ( -5.0,0.7)@xmath33 ( -12.0,0.7)@xmath33 ( -7.0,4.8)@xmath347 ( -0.9,1.4)@xmath187 ( -7.9,1.4)@xmath187 ( -14.1,4.8)@xmath348    in our analogy particles that require energy to be excited correspond to massive particles .
consequently after symmetry breaking there appear massless modes - the goldstone bosons - and modes that remain massive . if in qcd chiral symmetry was an exact symmetry we would expect to observe exactly massless goldstone bosons in the phase with broken symmetry .
however , it is only an approximate symmetry due to @xmath349 , leading to approximately massless bosons after symmetry breaking .
these are the 3 light pions observed in nature , @xmath350 and @xmath351 , with masses of about 135 mev compared to the heavier non - goldstone particles with masses of 800 -1000 mev .
recall that in the confined phase we only observe colourless bound states made from two quarks , mesons , which are bosons , from three quarks , baryons , which are fermions , or more .
the pions are unstable and decay mostly into a muon and its anti - neutrino for @xmath352 , the corresponding anti - particles for @xmath353 , or 2 photons for @xmath351 .
the number of goldstone bosons equals the dimension of the coset @xmath354 of the groups before and after breaking , which is @xmath355 , the dimension of @xmath356 in our case . to consider the pions as goldstone bosons thus means that we only keep the lightest quarks up and down , so that @xmath4 in our theory . for many purposes at low energy
this is a good approximation .
sometimes also the strange quark is included , with @xmath357 leading to consider the 8 lightest particles as goldstone bosons .
if we neglect all non - goldstone bosons in our theory - as they will not be relevant at low enough energies - and formally integrate them out we arrive at the low - energy effective theory of the goldstone bosons , @xmath358\exp\left[-\int d^4x \frac{f_\pi^2}{4}\mbox{tr}\left({\partial_\mu u(x)}{\partial_\mu u(x)^\dag}\right ) + \int d^4x\frac{\sigma}{2}\mbox{tr}\left(m(u(x)+u^\dag(x))\right )   \right ] \label{zchpt}\ ] ] where @xmath359 , \   m=\mbox{diag}(m_1,\ldots , m_{n_f } ) .
\label{u(x)}\ ] ] here the implicit sum over @xmath290 runs from 1 to @xmath355 . for @xmath4 the 3 fields @xmath360 actually denote the pion fields , with @xmath361 and @xmath362 . in eq .
( [ u(x ) ] ) we have already split off the constant , @xmath33-independent modes @xmath363 of the fields @xmath360 for later convenience . usually at this stage @xmath363 is set to the identity .
( [ zchpt ] ) represents the partition function or path integral of chiral perturbation theory ( chpt ) , which replaces that of qcd eq .
( [ zqcd ] ) in the low - energy regime . in eq .
( [ zchpt ] ) we only give the leading order ( lo ) lagrangian , all higher order terms in powers of @xmath364 and its derivatives that are invariant exist and we will come back to the question when they contribute in the next section . for more details about chpt we refer to @xcite .
chpt is an interacting , non - linear qft . when expanding the exponential in eq .
( [ u(x ) ] ) in powers of @xmath360 and considering only quadratic terms we obtain the standard lagrangian for scalar fields that leads to the klein - gordon equation . from that expansion
we can make the following identification for the masses @xmath365 of the pion fields @xmath360 , the gell - mann  oakes  renner ( gor ) relation : @xmath366 the qft in eq .
( [ zchpt ] ) is not yet related to rmt , and we turn to this limit in the next subsection .      the so - called epsilon - regime of chpt ( echpt ) was introduced by gasser and leutwyler @xcite before a rmt of qcd was constructed .
their motivation was to introduce a regime where analytic computations in chpt were feasible , in particular to understand the mass dependence of the partition function and the role of topology .
in particular leutwyler and smilga later computed the partition function analytically in this regime and deduced so - called sum rules for the dirac operator eigenvalues @xcite . working in a finite volume @xmath340 of linear size @xmath126
, the authors of @xcite introduced the following scaling : @xmath367 this scaling implies in particular @xmath368 where the former follows from the microscopic scaling limit of the banks - casher relation ( [ microlim ] ) , and the latter ensures that the kinetic term of the lo action is dimensionless .
the physical interpretation of the epsilon - regime follows from the scaling of the pion mass , which is implied from multiplying the gor relation eq .
( [ gor ] ) with @xmath18 .
it has to hold that @xmath369 , such that both sides are of the order @xmath370 . in other words ,
the compton wavelength of the pion fields is much larger than the size of the box : @xmath371 which is clearly an unphysical regime .
normally one would consider this to be undesirable , because in order to observe physical , propagating pions one would choose a scaling such that the pions do fit into the box @xmath18 , with @xmath372 .
in the so - called @xmath373-regime both momenta and pion mass are scaled in the same way , @xmath374 , where such an inequality can be satisfied .
we will come back to this issue very briefly below . in the epsilon - regime introduced above
we obtain the following limit @xcite @xmath375 \exp\left[-\int d^4x \frac12\partial_\mu\pi^b(x)\partial_\mu\pi^b(x ) \right]\int_{su(n_f)}du_0\exp\left[\frac12
\sigma v\mbox{tr}(m(u_0+u_0^\dag))+{\cal o}(\epsilon^2)\right ] .
\label{zechpt}\ ] ] here the propagating modes @xmath360 completely factorise from the zero - modes @xmath363 and contribute to the partition function as a free gaussian field theory - which is just a constant factor .
the reason that we do not expand in the constant modes of the pions as well , but keep them as a full matrix @xmath363 as in eq . ( [ u(x ) ] ) is that their quantum fluctuations are of @xmath376 , see @xcite for a more detailed discussion .
all higher order terms including those not given in the chiral lagrangian ( [ zchpt ] ) are subleading in this limit .
it is in this limit that the mass - dependence of the limiting partition function eq .
( [ zechpt ] ) and of rmt agree , and in fact also all dirac operator eigenvalue correlation functions . before we make this statement more precise
let us add further comments .
first , it is not only possible to make analytical computations in the epsilon - regime , but also to compare them to lattice qcd . on the lattice
the simulation parameters can be freely chosen , such that the epsilon - regime is reached and a quantitative comparison is possible .
second , it is possible to compute so - called one or higher loop corrections to the partition function ( and correlation functions ) eq . ( [ zechpt ] ) .
for example writing down the @xmath377 term explicitly , expanding the exponential and using wick s theorem with the finite volume propagator corresponds to computing the one - loop corrections .
third , one may perform also perturbative computations in the physical @xmath373-regime , or even chose an extrapolation between the two regimes , cf .
@xcite for a discussion and references .
let us now make the map between rmt and echpt more precise .
for that purpose it is convenient to introduce an extra term in the qcd lagrangian , the so - called theta - vacuum term containing the dual field strength tensor @xmath378 .
it is defined using the totally antisymmetric epsilon - tensor @xmath379 in 4 dimensions .
the extra term reads in appropriate normalisation : @xmath380 it is a topological term , and here @xmath24 is the number of left handed @xmath381 minus right handed @xmath382 zero modes of the dirac operator eq .
( [ dirac ] ) - which so far has not made an appearance in this section about qcd . due to the atiah - singer index theorem
@xmath24 is equal to the winding number of the gauge fields @xmath299 which is also a topological quantity . in the partition function eq .
( [ zqcd ] ) we have thus implicitly summed over all the topologies , which can now be written as follows : @xmath383 here @xmath384 denotes the qcd partition function at fixed topology .
inverting this equation it can be written as @xmath385 what is the advantage ?
in rmt in the previous section we always had a fixed @xmath386 - to which we should now compare for non - negative values . ] . furthermore , the unitary group integral over @xmath356 that we found in the epsilon - regime eq .
( [ zechpt ] ) can now be promoted to an integral over the full unitary group , according to @xmath387^\nu \label{groupint}\ ] ] defining the rescaled quark masses in the epsilon - regime as follows , @xmath388 we obtain the limiting relation between the mass - dependent part of the qcd partition function in the epsilon - regime at fixed topology and rmt : @xmath389}(\{\hat{m}\})&=&\int_{u(n_f ) } du_0\det[u_0]^\nu\exp\left[\frac12 \mbox{tr}(\hat{m}(u_0+u_0^\dag))\right]\nonumber\\ & = & 2^{\frac{n_f(n_f-1)}{2}}\prod_{j=0}^{n_f-1}\gamma(j+1 ) \frac{\det_{1\leq i , j\leq n_f}[\hat{m}_i^{j-1}i_\nu(\hat{m}_i)]}{\delta_{n_f}(\{\hat{m}^2\})}= % \lim_{\stackrel{\mbox{\small $ n\to\infty$}}{\mbox{\small $ m_f\to0 $ } } }   { \cal z}^{[n_f]}(\{\hat{m}\})\ . \label{lsint}\end{aligned}\ ] ] in the last equality with the limiting rmt partition function from eq .
( [ zrmtnf ] ) we have identified rescaled random matrix and qcd quark masses @xmath344 for all flavours . the computation of the integral over the appropriately normalised haar measure of the unitary group @xmath390 in eq .
( [ lsint ] ) goes back to @xcite , see also @xcite and @xcite for our normalisation .
a concise derivation in terms of group characters was presented in @xcite . in the simplest case with @xmath248
it is easy to see that : @xmath391}(\hat{m})=\int_{-\pi}^\pi\frac{d\theta}{2\pi } e^{i\nu\theta}\exp\left [ \frac12 \hat{m}(e^{i\theta}+e^{-i\theta})\right]=i_\nu(\hat{m})\ , \label{zechptnf1}\ ] ] which is one of the standard definitions of the modified bessel - function .
furthermore , we can also verify the completely degenerate case @xmath392 , by applying the andrief formula ( [ andreiev ] ) .
diagonalising @xmath363 by a unitary transformation @xmath393 , @xmath394 , we obtain again the squared vandermonde determinant as a jacobian : @xmath395^\nu\exp\left[\frac{\hat{m}}{2}\mbox{tr}(u_0+u_0^\dag)\right]=const .
\prod_{j=1}^{n_f}\left(\int_{-\pi}^\pi\frac{d\theta_j}{2\pi}\exp[i\nu\theta_j+\hat{m}\cos(\theta_j)]\right ) |\delta_{n_f}(\{e^{i\theta}\})|^2\nonumber\\ & = & n_f!\ const.\det_{1\leq j , k\leq n_f}\left [ \int_{-\pi}^\pi\frac{d\theta}{2\pi}\exp[i(\nu+k - j)\theta+\hat{m}\cos(\theta)]\right]=\det_{1\leq j , k\leq n_f}[i_{\nu+k - j}(\hat{m})]\ . \label{zdeg}\end{aligned}\ ] ] it agrees with eq .
( [ zrmtnfdeg ] ) from rmt , after appropriately normalising .
we observe that both for degenerate and non - degenerate masses the partition functions with @xmath0 flavours eqs .
( [ lsint ] ) and ( [ zdeg ] ) are determinants of single flavour partition functions with @xmath248 .
the matching of the limiting rmt partition function and the group integral from the epsilon - regime of chpt at fixed topology is only part of their equivalence . in particular
this matching does not imply that all dirac operator eigenvalue correlation functions agree .
for example at @xmath248 the echpt partition function agrees for all three symmetry classes @xmath52 @xcite , while their densities differ . at first sight it is not clear how to access the dirac operator spectrum , once we have moved from the quarks to the goldstone bosons as fundamental low energy degrees of freedom .
however , looking at the alternative construction of the spectral density via the resolvent eq .
( [ resolvent ] ) helps .
consider adding an auxiliary pair of quarks with masses @xmath33 and @xmath187 respectively to the qcd lagrangian eq .
( [ zqcd ] ) .
if we arrange for the second quark to be bosonic , we obtain the following average with respect to the unperturbed qcd partition function , at fixed topology : @xmath396}{\det[y - d\!\!\!\!/\ ] } \right\rangle_\nu^{\rm qcd}= \frac{\int_{u(n_f+1|1 ) } du_0\ \mbox{sdet}[u_0]^\nu\exp\left[\frac12 \mbox{str}(\hat{{\cal m}}(u_0+u_0^\dag))\right ] } { { \cal z}_{\rm echpt}^{[n_f]}(\{\hat{m}\})}\ , \label{qcdresolvent}\ ] ] with @xmath397 , in the limit of the epsilon - regime .
because of the bosonic nature of the second auxiliary quark in the denominator the integral is now over the supergroup @xmath398 with the corresponding supersymmetric trace and determinant , str and sdet respectively , see e.g. chapter 7 of @xcite .
its evaluation and differentiation to obtain the resolvent has to be done carefully . in @xcite to where we refer for details
it was shown in this fashion that the microscopic density eq .
( [ besseldensity ] ) and its extension to arbitrary @xmath0 follows from eq .
( [ qcdresolvent ] ) , thus establishing the equivalence to rmt on the level of the microscopic density .
in order to know the distribution of the smallest eigenvalues the knowledge of all density correlation functions is necessary , see eq .
( [ expand ] ) , as was pointed out in @xcite .
the equivalence proof for all @xmath62-point density correlation functions was achieved in @xcite by directly matching the @xmath62-point resolvent generating functions with @xmath399 additional ratios of quark determinants of both theories , using superbosonisation techniques .    a few comments regarding universality and corrections to the rmt regime
are in place here .
already the fact that two apparently different theories as echpt and rmt yield the same dirac operator correlation functions is a striking result which is based on symmetries and universality .
it is thus surprising that even when taking into account one - loop corrections to echpt the same universal correlation functions were found in @xcite .
the only corrections to be made are subsumed in a modification @xmath400 , and @xmath401 , where @xmath402 and @xmath403 are constants that depend on @xmath0 , @xmath337 and the geometry of the discretisation , see e.g. @xcite for details and references .
so far @xmath337 did not appear in eigenvalue correlation functions , but later when including a chemical potential it will .
this modification means that the chiral condensate ( and pion decay constant ) and thus the rescaled masses in eq .
( [ hatm ] ) and eigenvalues simply get renormalised accordingly . only two - loop corrections within echpt contain non - universal terms that lead out of the rmt universality class @xcite .
however , even at the lo the agreement between the rmt and echpt rapidly breaks down when going to higher energies , meaning to higher rescaled dirac operator eigenvalues . at some point
the fluctuating modes will start to contribute , and eventually also the chpt approach breaks down as contributions from full qcd will appear .
the scale where this happens has been called thouless energy @xcite in analogy to applications in condensed matter physics , to where we refer for details .
the point where this happens scales with @xmath404 and @xmath405 . beyond this scale one sees a rise of the global spectral density as in eq .
( [ bcss ] ) , away from the rmt prediction .
the scaling of the thouless energy has been confirmed from lattice data in a different theory , in qft with two colours @xcite .
for the same theory deep in the bulk regime the local statistics was seen to agree again with rmt bulk statistics of the gse @xcite , which is the non - chiral @xmath232 symmetry class .
here we have already mentioned one of the two different symmetry classes corresponding to @xmath48 and 4 in eq .
( [ partitionfunct ] ) .
they are also relevant but for different qfts other than qcd , e.g. with only two colours @xmath406 .
these have chiral symmetry patters different from eq .
( [ xsbpattern ] ) and we refer to @xcite for a detailed discussion . the matching with the 3 chiral rmts in eq .
( [ partitionfunct ] ) was initially proposed in @xcite . while the agreement with lattice simulations of the respective theories ( cf .
@xcite ) leave little doubt that these rmt are equivalent to the epsilon - regime of the corresponding chpt , such an equivalence has not been established even on the level of the partition function for general @xmath0 .
this is due to the more complicated structure of the corresponding group integral , see however @xcite and @xcite for a discussion .
let us briefly comment also on the influence of the dimension . in 3 space - time dimensions
there is no chiral symmetry , and the rmt for qcd is given by the gue instead @xcite , including mass terms .
in 2 dimensions which again has chiral symmetry the classification is much richer , and we refer to @xcite for details .
this makes contact with the symmetry classification of topological insulators @xcite and the so - called bott - periodicity in any dimension .
in this section we will cover some recent developments in the application of rmt to qcd , where the detailed dependence of the partition function and dirac operator eigenvalue correlation functions on finite lattice spacing @xmath19 or chemical potential @xmath10 are computed .
in the first subsection [ wrmt ] we will study the influence of a finite space - time lattice on the spectrum as it is expected to be seen in numerical simulation of qcd on a finite lattice far enough from the continuum . the corresponding rmt and effective field theory that we will consider go under the name of wilson ( w)rmt as introduced in @xcite , and wilson ( w)chpt , see @xcite for a standard reference .
we will only sketch the ideas beginning with the symmetries and wchpt , and then write down the joint density of eigenvalues of the hermitian wilson dirac operator .
its solution and its non - hermitian part will be referred to the literature .
the addition of a chemical potential term for the quarks to the dirac operator renders its eigenvalues complex .
the solution of the corresponding rmt using op in the complex plane is sketched in the second subsection [ rmtmu ] where we will give some details .
the first application of such a rmt was already introduced by stephanov in 1996 @xcite to explain the difference between the quenched and unquenched theory .
the development of the corresponding theory of op in the complex plane started later in @xcite and @xcite .
it is already partly covered in the les houches lecture notes from 2005 @xcite , for a slightly more recent review on rmt of qcd with chemical potential see @xcite . here
we will comment on some recent developments relating to products of random matrices , see @xcite for a detailed review ( and the lectures of y. tourigny in this volume ) .
the discretisation of derivatives on a space - time lattice is not a unique procedure , which also applies to the dirac operator in eq .
( [ dirac ] ) . we would thus first like to motivate why in this subsection we will study the modification proposed by wilson , the wilson dirac operator @xmath407 which is no longer anti - hermitian . here
@xmath19 is the lattice spacing and @xmath408 denotes the laplace operator which is hermitian .
of course both operators still have to be discretised .
we can not possibly give justice to the vast literature on this subject and will only focus on aspects relevant for the application of rmt .
the reason why wilson proposed to add the laplacian is the so - called doubler problem . in the continuum with @xmath409 the relativistic energy - momentum relation for a particle with mass @xmath410 , energy @xmath411 , and 4-momentum @xmath412 reads @xmath413 here we have spelled out the sum explicitly .
particles satisfying this relation are called on - shell ( on the energy shell ) . after standard discretisation
this relation becomes @xmath414 which in the limit @xmath20 leads back to eq .
( [ emom ] ) .
in contrast to the continuum relation , here with any @xmath412 satisfying this equation also @xmath415 , @xmath416 etc .
fulfil eq .
( [ emoma ] ) . in total
we have @xmath417 possibilities to add @xmath418 to the components of @xmath412 , which all correspond to the same on - shell particle .
this is the doubler problem .
wilson s modification to add the laplacian to the dirac operator amounts to adding the term @xmath419 to the right hand side of eq .
( [ emoma ] ) . in the limit
@xmath20 this makes the extra 15 particles heavy and thus removes the doublers .
the modification comes with a prices , as the laplacian explicitly breaks chiral symmetry .
we will not discuss other possibilities here , for example in @xcite for a rmt of the so - called staggered dirac operator , but rather stick to wilson s choice . despite the non - hermiticity of @xmath420 eq .
( [ wdirac ] ) it satisfies the so - called @xmath284-hermiticity : @xmath421 here we have simply used that @xmath422 and the anti - commutator from eq .
( [ dchiral ] ) .
consequently one can define the following hermitian operator called @xmath423 @xmath424 where traditionally the quark mass @xmath425 is added to the definition .
looking at eq .
( [ dchiral ] ) for the dirac operator and eq .
( [ gammas ] ) for @xmath284 we can immediately give the chiral block structure of @xmath420 and @xmath423 : @xmath426 here @xmath427 and @xmath428 are hermitian operators .
note the change in sign in the mass term in the lower right block of @xmath423 . as previously for @xmath409 this immediately leads to a good ansatz to make for a wrmt , as it was made in @xcite .
before we turn to analyse this in some detail let us directly go to the epsilon - regime of qcd in the wilson formulation at fixed topology , wechpt .
it turns out that at lo 3 extra terms appear in the chiral lagrangian as was discussed in @xcite ( note the different sign convention ) .
they contain 3 new low energy constants ( lec ) @xmath429 , @xmath430 and @xmath431 , and all terms are proportionally to @xmath432 : @xmath433}(\{\hat{m}\})&= & \int_{u(n_f ) } du_0\det[u_0]^\nu\exp\big[+\frac12 \sigma v\mbox{tr}({m}(u_0+u_0^\dag ) ) -a^2v w_8\mbox{tr}(u_0 ^ 2+u_0^{\dag\,2})\nonumber\\ & & \quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad -a^2 vw_6(\mbox{tr}(u_0+u_0^\dag))^2-a^2 vw_7(\mbox{tr}(u_0-u_0^\dag))^2 \big]\ .
\label{zwechpt}\end{aligned}\ ] ] it is clear that this implies the scaling @xmath434 or @xmath435 , in addition to the scaling of the quark masses with the volume in eq .
( [ exponenteps ] ) .
let us discuss some simplifications first . without fixing topology in the special case of @xmath436 the term tr@xmath437 vanishes , and the terms proportional to @xmath431 and @xmath429 are equivalent . in general , in the second line of eq .
( [ zwechpt ] ) the 2 squares of the traces can be linearised by so - called hubbart - stratonovich transformations , at the expense of 2 extra gaussian integrals .
the linearised terms can be included into a shift of the mass for @xmath429 , and into so - called axial mass terms @xmath438 for @xmath430 , see @xcite for more details . also for the latter term
the corresponding group integral generalising eq .
( [ lsint ] ) is known , see @xcite .
therefore in the following we will set @xmath439 , having in mind that we need to perform 2 extra integrals for the full partition function ( [ zwechpt ] ) for non - zero values of these lecs .
let us define @xmath440 and compute the partition function only including this term .
for @xmath248 we obtain @xmath441}(\hat{m})=\int_{-\pi}^\pi\frac{d\theta}{2\pi } e^{i\nu\theta+\hat{m}\cos(\theta)-4\hat{a}_8 ^ 2\cos(\theta)^2 + 2\hat{a}_8 ^ 2 }   = e^{2\hat{a}_8 ^ 2}\int_{-\infty}^\infty \frac{dx}{\sqrt{\pi}}e^{-x^2}\left(\frac{\hat{m}-4ix\hat{a}_8}{\hat{m}+4ix\hat{a}_8}\right)^{-\frac{\nu}{2 } } i_\nu\left(\sqrt{\hat{m}^2 + 16x^2\hat{a}_8 ^ 2}\right ) .
\label{zwechptnf1}\ ] ] it is no longer elementary , but after linearising the @xmath442 term it can be written as a gaussian integral over the one - flavour partition function @xmath443 in the continuum at shifted mass . for
degenerate masses and general @xmath0 we can again diagonalise @xmath363 and use andrief s integral formula , as in the derivation of eq .
( [ zdeg ] ) , leading to : @xmath444}(\hat{m})=\det_{1\leq j , k\leq n_f}\left[e^{2\hat{a}_8 ^ 2 } \int_{-\pi}^\pi\frac{d\theta}{2\pi}\exp\left[i(\nu+k - j)\theta+\hat{m}\cos(\theta)-4\hat{a}_8 ^ 2\cos(\theta)^2\right]\right ] .
\label{zwechptnf}\ ] ] once again it is given by the determinant of @xmath248 flavour partition functions . in order to determine the spectral density of @xmath423 or @xmath420 one then has to introduce additional auxiliary quark pairs as described in eq .
( [ qcdresolvent ] ) ( or use replicas ) .
prior to the rmt calculation this was done for the spectral density from wechpt in @xcite , for both @xmath423 or @xmath420 .
let us now turn to the rmt side as introduced in @xcite . for simplicity we will only consider the hermitian operator @xmath423 , for the discussion of the complex eigenvalue spectrum of @xmath420 we refer to @xcite and to @xcite for its complete solution .
actually two slightly different rmt have been proposed , which have turned out to be equivalent on the level of the jpdf . in @xcite
the partition function ( [ matrixrep ] ) was generalised using eq .
( [ dwchiral ] ) in an obvious way , by adding two hermitian gaussian matrices @xmath445 and @xmath313 of dimensions
@xmath34 and @xmath56 , respectively : @xmath446}\sim\int[dw][da][db]\prod_{f=1}^{n_f}\det \left [ \begin{array}{cc } m_f 1_n+aa&i w\\ -iw^\dag & -m_f 1_{n+\nu}-ab\\ \end{array } \right ] \exp\left[-\frac14\mbox{tr}(a^2+b^2 + 2ww^\dag)\right ] .
\label{wimatrixrep}\ ] ] in @xcite the following modification was made : rather than adding a block - diagonal hermitian matrix to @xmath423 as the last term in eq .
( [ dwchiral ] ) , a full hermitian matrix @xmath447 was added to it , @xmath448 , filling the off - diagonal blocks with a rectangular complex matrix @xmath449 ( and trivially changing the sign of @xmath450 ) : @xmath451}&=&\int[dw][dh]\prod_{f=1}^{n_f}\det \left [ \begin{array}{cc } m_f 1_n+a&i w+\omega\\ -iw^\dag + \omega^\dag&-m_f 1_{n+\nu}+b\\ \end{array } \right ] \exp\left[-\frac{1}{2(1-a^2)}\mbox{tr}ww^\dag-\frac{1}{4a^2}\mbox{tr}h^2\right],\nonumber\\ & & h=\left(\begin{array}{cc } a&\omega\\ \omega^\dag&b\\ \end{array}\right),\ \ a\in[0,1]\ . \label{wiimatrixrep}\end{aligned}\ ] ] the different rescaling of the variances with @xmath19 was made in @xcite to underline that on the level of rmt this corresponds to a two - matrix model that interpolates between the chgue and the gue , in the limits @xmath20 and @xmath452 , respectively .
it was shown in @xcite and @xcite that in the limit @xmath32 , identifying @xmath453 in addition to the rescaled mass terms the two rmt partition functions and eq . ( [ zwechpt ] ) with @xmath454 agree : @xmath444}=\lim_{\stackrel{\mbox{\small $ n\to\infty$}}{\mbox{\small $ a , m_f\to0 $ } } } { \cal z}_{wi,\,n}^{[n_f]}= \lim_{\stackrel{\mbox{\small $ n\to\infty$}}{\mbox{\small $ a , m_f\to0 $ } } } { \cal z}_{wi\!i,\,n}^{[n_f]}\ .
\label{wequiv}\ ] ] as we have said before the presence of the @xmath429- and @xmath430-terms can be achieved by 2 extra gaussian integrals , see @xcite .
it is quite remarkable that rmt allows to include such a detailed structure of wchpt .
the fact that we have `` only '' considered rmts for @xmath423 does not affect the equivalence argument of the partition functions . because of eq .
( [ d5 ] ) and @xmath455=1 $ ] the partition functions of @xmath420 and @xmath423 agree .
let us just state the jpdf of the eigenvalues @xmath456 of @xmath423 for degenerate masses @xmath425 without derivation , see @xcite and @xcite for details , @xmath457}(m ) & \sim & \exp\big[\frac{nm^{2}}{2a^2(1-a^2)}\big ] \int_{-\infty}^\infty\prod_{j=1}^{2n+\nu } dd_jd_j^{n_f } \exp\big[{-\frac{d_j^2}{4a^2}}\big ] \delta_{2n+\nu}(\{d\})\ \ \ \
\nonumber\\ & & \times\mbox{pf}_{1\leq i , j\leq 2n+\nu;\,1\leq q\leq \nu } \left [ \begin{array}{cc } f(d_j - d_i ) & \!d_i^{\ , q-1}\exp\left[-\frac{d_im}{2a^2}\right]\\   -d_j^{\ , q-1}\exp\left[-\frac{d_jm}{2a^2}\right ] & { 0}\\ \end{array } \right]\ , \label{zdevpf}\end{aligned}\ ] ] where we have defined the antisymmetric weight @xmath458   \sinh\big[\frac{xu}{2a^2}\big].\ ] ] this jpdf has two new features .
first , there appears a vandermonde times a pfaffian ( pf ) defined as @xmath459=\sqrt{\det[a]}$ ] for antisymmetric matrices @xmath445 of even dimension , see e.g. @xcite for a definition in terms of permutations . this is quite common for such interpolating two - matrix models , see @xcite for classical works regarding the gue - goe and the gue - gse transitions .
second , and most importantly the solution for the correlation functions involves both skew op and op at the same time , where the op are hermite polynomials of a gue of size @xmath460 . we refer to @xcite for details , in particular that this sub - structure of a finite size gue remains valid in the large-@xmath34 limit . in @xcite where only the cases @xmath227 and 1 were solved this structure
was not observed , because only hermite polynomials of degree 0 and 1 occur .    instead of giving any further details of the solution of wrmt we show plots in figure [ fig : wbessel ] below for the limiting quenched microscopic density @xmath461 of @xmath423 for @xmath212 , illustrating the effect of the rescaled finite lattice spacing @xmath462 .
in order to compare to the results from the chgue in figure [ fig : besseldensity ] let us first discuss the difference between @xmath423 and @xmath22 at @xmath409 , where @xmath463 , see eqs .
( [ wdirac ] ) and ( [ d5 ] ) .
the rotation with @xmath284 merely makes the dirac operator hermitian rather than anti - hermitian , which corresponds to rotating the eigenvalues from the imaginary to the real axis ( as we had already done in writing the jpdf in eq .
( [ partitionfunct ] ) of real eigenvalues ) .
the shift by @xmath464 however introduces a gap such that there are no eigenvalues of @xmath423 in @xmath465 $ ] .
furthermore , the @xmath24 exact zero eigenvalues which are not shown in figure [ fig : besseldensity ] get moved away from the origin to one of the edges of this gap , which in our convention is at @xmath466 .
therefore we now have to give the density on @xmath111 instead of @xmath112 , and it holds for @xmath25 : @xmath467 the effect of @xmath468 is both to broaden the @xmath24 delta - functions and to eventually fill the gap @xmath469 $ ] . because eigenvalues repel each other the @xmath24 zero - modes will now spread , and form approximately a gue of finite size @xmath24 , see @xcite for a further discussion . in order to make the comparison more transparent
let us also choose @xmath470 , so that we can plot the density on @xmath112 only , as is shown in figure [ fig : wbessel ] .
the effect of @xmath468 is particularly striking on the zero - modes .
( 12.2,5.5 )   for @xmath227 ( left ) and @xmath264 ( right ) as a function of @xmath471 , taken from @xcite . because of choosing @xmath470 it is still symmetric around the origin . for @xmath472
we are back to the microscopic density of bessel functions as in figure [ fig : besseldensity ] , whereas for @xmath473 ( blue ) and @xmath474 ( red ) the density starts to spread into the region around the origin . for @xmath264
the delta - function of the one zero - eigenvalue is broadened by increasing @xmath471 - which is why this is plotted on a different scale compared to @xmath227 . at these small values of @xmath471 the effect on the density is remarkably localised to the origin .
only in the limit @xmath475 we will arrive at the microscopic density of the gue @xcite which is completely flat ( and normalised to @xmath476 here ) .
[ fig : wbessel],title="fig:",width=264 ]   for @xmath227 ( left ) and @xmath264 ( right ) as a function of @xmath471 , taken from @xcite .
because of choosing @xmath470 it is still symmetric around the origin .
for @xmath472 we are back to the microscopic density of bessel functions as in figure [ fig : besseldensity ] , whereas for @xmath473 ( blue ) and @xmath474 ( red ) the density starts to spread into the region around the origin . for @xmath264 the delta - function of the one zero - eigenvalue is broadened by increasing @xmath471 - which is why this is plotted on a different scale compared to @xmath227 .
at these small values of @xmath471 the effect on the density is remarkably localised to the origin . only in the limit @xmath475
we will arrive at the microscopic density of the gue @xcite which is completely flat ( and normalised to @xmath476 here ) .
[ fig : wbessel],title="fig:",width=264 ] ( 0.1,-0.2)@xmath265 ( -7.1,-0.2)@xmath265 ( -10.8,3.5)@xmath461 ( -3.5,1.4)@xmath461    the density of the hermitian wilson dirac operator @xmath423 has been compared to quenched lattice simulations in @xcite . because this involves multi - parameter fits several quantities
have been proposed that are simple to measure @xcite .
a comparison to unquenched data remains a difficult question , see however @xcite .
the rmt approach presented here has also lead to insights about the sign and constraints amongst the lecs @xmath477 that were previously controversial @xcite .
when the lattice artefacts become very strong additional unphysical phases appear , and the possibility to access these as a function of the lec and quark content @xmath0 was clarified in @xcite .
the derivation of all @xmath62-point correlation functions including @xmath478 based on wechpt has been pushed forward in @xcite , leading to alternative representations compared to the op approach @xcite .
for the density of @xmath423 in the @xmath373-regime see @xcite .
the computation of individual eigenvalue distribution is not known beyond the expansion of the fredholm expansion proposed in @xcite .      in this subsection
we will study the influence of a chemical potential @xmath10 on the dirac operator spectrum .
we will first state its global symmetry in qcd , then construct the rmt and at the end link to the epsilon - regime of the corresponding chiral lagrangian .
the addition of a chemical potential @xmath479 for each quark flavour to the qcd action in eq .
( [ zqcd ] ) amounts to add the terms @xmath480 . in classical statistical mechanics
the addition of a chemical potential allows the number of the corresponding particles to fluctuate . here
we will add the same chemical potential @xmath481 to all flavours ( which makes it the baryon chemical potential ) .
in analogy to eq .
( [ dchiral ] ) the global symmetry of the dirac operator plus @xmath482 changes to @xmath483 and thus becomes complex non - hermitian .
it holds that @xmath484 , and thus for purely imaginary chemical potential @xmath485 it is still anti - hermitian . because the operator remains chiral @xmath486 the non - zero complex eigenvalues continue to come in pairs @xmath487 , using the same argument as after eq .
( [ eigenfunct ] ) .
the fact that the spectrum of the dirac operator with chemical potential is complex ( and thus the qcd action too ) has severe repercussions for lattice qcd . because the action can no longer be interpreted as a probability weight , the standard method on the lattice , weighting a given configuration by its action breaks down .
several methods have been invented to circumvent this so - called sign problem , and we refer to @xcite for a review .
in contrast the rmt for qcd with chemical potential remains analytically solvable with op in the complex plane , as we will see .    in @xcite
a rmt was constructed by stephanov where @xmath298 was replaced by a complex gaussian random matrix @xmath488 , cf .
( [ matrixrep ] ) for @xmath47 compared to ( [ dchiral ] ) . while this rmt was very useful to illuminate the difference between @xmath25 and @xmath88 , which is more dramatic than for @xmath489 , this model was not analytically tractable with op at finite-@xmath34 .
therefore osborn @xcite introduced a two - matrix model where @xmath10 is multiplied by a second random matrix @xmath490 , see eq .
( [ 2matrixrep ] ) below , implying that the chemical potential term is not diagonal in that basis . while at first sight more complicated this model turns out to be exactly solvable using op in the complex plane , as we will demonstrate .
a further equivalent representation was chosen in @xcite , multiplying @xmath491 by @xmath492 $ ] instead .
this parametrisation allowed to understand , why the sign problem in the numerical solution of this rmt for qcd with @xmath493 can be circumvented @xcite . prior to introducing an rmt with @xmath10 as in eq .
( [ dchiralmu ] ) a similar model was introduced @xcite describing the effect of temperature @xmath13 .
such a model can only be solved using the supersymmetric method , see @xcite , with the effect of renormalising @xmath494 to be temperature dependent , keeping the correlation functions otherwise unchanged .
a further modification of stephanov s model by including both parameters @xmath495 in eq .
( [ dchiralmu ] ) was very successful as a schematic model for the phase diagram of qcd @xcite , predicting the qualitative features depicted in figure [ fig : phasediag ] .
let us turn to the rmt @xcite defined as @xmath496[dw_2]\prod_{f=1}^{n_f}\det \left [ \begin{array}{cc } m_f 1_n&i w_1+\mu w_2\\ iw_1^\dag + \mu w_2^\dag & m_f 1_{n+\nu}\\ \end{array } \right ] \exp[-\mbox{tr}(w_1w_1^\dag + w_2w_2^\dag)]\nonumber\\ & \sim & \int[dx_1][dx_2]\prod_{f=1}^{n_f}\det \left [ \begin{array}{cc } m_f 1_n&x_1\\ x_2&m_f 1_{n+\nu}\\ \end{array } \right ] \exp\left[-a(\mu)\mbox{tr}(x_1x_1^\dag + x_2x_2^\dag)\right]\nonumber\\ & & \quad\quad\quad \times\exp\left[b(\mu)\mbox{tr}(x_1x_2+x_2^\dag x_1^\dag)\right]\ , \ \
a(\mu)=\frac{1+\mu^2}{4\mu^2}\ , \ \ b(\mu)=\frac{1-\mu^2}{4\mu^2}\ . \label{2matrixrep}\end{aligned}\ ] ] in the second step
we have simply changed variables to the matrices @xmath497 and @xmath498 which are now coupled .
while at @xmath499 the matrices @xmath500 and @xmath501 are again independent gaussian random matrices - this is called maximal non - hermiticity - at @xmath489 we have @xmath502 with total correlation ( and we are back to the chgue ) .
a look at the characteristic equation for the massless rmt dirac operator , @xmath503=\det[z^2-x_1x_2]\ , \label{muchar}\ ] ] reveals that we are after the complex eigenvalues in this setting see @xcite .
] of the product of the two correlated matrices @xmath504 . as in eq .
( [ matrixrep ] ) we could have allowed for real or quaternion valued matrix elements with @xmath67 .
however , the joint density of complex eigenvalues does not enjoy a closed form for all 3 @xmath63 s as in eq .
( [ partitionfunct ] ) .
in particular real matrices are special as they can have real eigenvalues and complex conjugated eigenvalues pairs , as the characteristic equation ( [ muchar ] ) is real .
for the solution of the respective rmts see @xcite and @xcite , respectively .
the parametrisation of a generic complex non - hermitian matrix @xmath505 in terms of its complex eigenvalues is usually done using the schur decomposition @xmath506 , where @xmath507 is unitary , @xmath508 contains the complex eigenvalues , and @xmath13 is a strictly upper triangular complex matrix ( it can be chosen to be strictly lower triangular instead ) . for two matrices however , in order to compute the jacobian we have to use the following generalised schur decomposition @xcite @xmath509 with @xmath510 and @xmath511 . here
are unitary , @xmath513 are diagonal matrices with complex entries ( which are not the complex eigenvalues of @xmath500 and @xmath501 ) , and @xmath514 are strictly upper triangular .
the complex eigenvalues of the product matrix @xmath505 are thus given by the elements of the diagonal matrix @xmath515 .
let us directly give the result for the joint density from @xcite ( see @xcite for an alternative derivation ) : @xmath516}&= & \left(\prod_{j=1}^n\int d^2z_j|z_j|^{\nu}k_\nu(a(\mu)|z_j|)\exp\left[-\frac{b(\mu)}{2}(z_j+z_j^*)\right]\prod_{f=1}^{n_f}(z_j+m_f^2)\right ) |\delta_n(\{z\})|^2 \label{zrmtmu } \\
& = & \left(\prod_{j=1}^n\int d^2z_j\right ) { \cal p}_{jpdf}(z_1,\ldots , z_n)\ .
\label{jpdfmu}\end{aligned}\ ] ] note that in view of eq .
( [ muchar ] ) we should later take the square root of the eigenvalues @xmath517 of @xmath504 to change to dirac operator eigenvalues . in the complex plane
this is not unique and without loss of generality one can then restrict to the half - plane , also in view of the @xmath518 pairs of dirac eigenvalues . before we discuss the properties of eq .
( [ zrmtmu ] ) and its correlation functions let us try to understand why the jacobian of the transformation ( [ schur ] ) does not only include the modulus squared vandermonde determinant ( as in the complex ginibre ensemble , see @xcite ) and the zero - modes @xmath519 , but also a modified @xmath520-bessel function , even though we started from ( coupled ) gaussian matrices .
the same phenomenon happens when looking at products of complex gaussian random variables @xmath521 , and asking for the distribution @xmath522 of the product random variable @xmath523 : @xmath524\delta^{(2)}(z - z_1z_2)=(2\pi)^2 \int_0^\infty \frac{drr}{r^2}\exp[-r^2-|z|^2/r^2]= 2\pi k_0(2|z|)\ . \label{kweight}\ ] ] obviously we are dealing with the analogue of quadratic matrices here .
multiplying @xmath525 such random variables we obtain an @xmath526-fold integral representation for the distribution of @xmath527 now given by a meijer g - function : @xmath528 see @xcite for its standard definition .
these weights appear when studying the complex eigenvalues of the product of @xmath136 independent gaussian matrices @xmath529 and we refer to @xcite for the solution of this model .
the complex eigenvalue correlation functions of the model ( [ zrmtmu ] ) are defined as in eq .
( [ rk ] ) and it holds @xmath530\ . \label{rkc}\ ] ] they are again given in terms of the kernel @xmath531 of polynomials @xmath532 and @xmath533 which are bi - orthogonal in the complex plane with respect to the weight @xmath534\prod_{f=1}^{n_f}(z_j+m_f^2)\ ,   \label{muweight}\\ & & \ \
\int d^2z\ w(z ) p_k(z)q_l(z)=h_k\delta_{kl}\ .
\label{opc}\end{aligned}\ ] ] we will not discuss gap probabilities or individual eigenvalue distributions here and refer to @xcite .
unlike in the real eigenvalue case they will now depend on the choice of the geometry of the region containing @xmath535 eigenvalues , and only in the rotationally invariant case @xmath499 a radial ordering of eigenvalues according to their modulus provides a natural choice .
before we solve the rmt eq .
( [ zrmtmu ] ) a word of caution is in order . in standard op theory in the complex plane only real positive weights are considered , which is only true for @xmath25 in eq .
( [ muweight ] ) . in this case
@xmath536 and standard techniques as gram - schmidt apply to construct these polynomials , see e.g. @xcite .
when @xmath88 in our case the weight is complex , and thus a priori the partition function is as well as the @xmath537 are .
therefore the density no longer has a probabilistic interpretation , it has a real and imaginary part , and we refer to @xcite for a more detailed discussion
.    we will first provide the solution of the simplest quenched case with @xmath25 . then the partition function and general eigenvalue correlation functions including @xmath0 mass terms can be obtained along the same lines as in section [ opdirac ] , expressing them in terms of the op at @xmath25 .
surprisingly , despite the complicated non - gaussian weight function in the complex plane the op are again given by laguerre polynomials of complex arguments @xcite . )
was made in @xcite . ]
the following holds , where we drop the argument @xmath10 of the parameters @xmath538 in eq .
( [ 2matrixrep ] ) : @xmath539 l_j^\nu(cz)l_k^\nu(cz)=\delta_{jk}h_k\ , \
\mbox{with } \ \
c=\frac{a^2-b^2}{2b}\ , \ \
\label{opmu}\ ] ] here @xmath540 is an integer , and the laguerre polynomials can be made monic as in eq .
( [ laguerre ] ) .
the proof we sketch follows @xcite and uses induction of depth 2 , for an earlier longer proof see @xcite .
insert the following integral representations into the orthogonality relation ( [ opmu ] ) for @xmath212 , @xmath541\ , \nonumber\\ l_j^{\nu}(z)&=&\oint_\gamma \frac{du}{2\pi } \frac{\exp\left[-\frac{zu}{1-u}\right]}{(1-u)^{\nu+1}u^{j+1}}\ , \label{cintrep}\end{aligned}\ ] ] where @xmath542 is a contour enclosing the origin in positive direction , but not the point @xmath543 .
then the gaussian integrals over the real and imaginary parts of @xmath544 can be solved , and subsequently the integral over @xmath545 from the @xmath520-bessel function .
the orthogonality then easily follows applying the residue theorem twice . for the induction we use the 3-step recurrence relation for the laguerre polynomials and @xmath520-bessel function @xcite ,
see @xcite for details .
we thus have solved the quenched case as we know all correlation functions from the kernel @xmath546 in terms of the quantities in eq .
( [ opmu ] ) , by using eq .
( [ rkc ] ) .
a similar solution exists @xcite for the so - called elliptic complex ginibre ensemble of a single complex matrix @xmath547 , where the hermitian part @xmath548 and anti - hermitian part @xmath549 of the matrix have different variances that depend on a parameter .
the solution is given in term of hermite polynomials in the complex plane and allows to interpolate between the gue and the complex ginibre ensemble at maximal non - hermiticity , where @xmath447 and @xmath445 have the same variance .
we refer to @xcite for a detailed discussion .
let us now turn to the unquenched solution with @xmath88 described in great detail in @xcite , expressing everything in terms of the laguerre polynomials of the quenched solution from eq .
( [ opmu ] ) .
because we follow the same strategy as in subsections [ genweight ] and [ alllaguerre ] using expectation values of characteristic polynomials we can be very brief . the corresponding proofs for manipulating expectation values of characteristic polynomials and their complex conjugates defined as in eq . ( [ ave ] )
can be found in @xcite . in particular for the monic op in the complex plane satisfying eq .
( [ opc ] ) for a weight @xmath522 that can be complex we have @xcite @xmath550 the proof goes exactly the same way as in eq .
( [ proofpn ] ) , where for @xmath551 we multiply the product into @xmath552 , increasing it to @xmath553 .
for @xmath554 we simply choose @xmath555 , increasing it to @xmath556 .
likewise for the kernel we obtain @xmath557 modifying the proof in eq .
( [ proofkernel ] ) accordingly by multiplying the first product into @xmath558 , and the second into @xmath559 .
the only difference here is that in general no christoffel - darboux formula is available .
this is despite the fact that the laguerre polynomials in our example eq . ( [ opmu ] ) still satisfy the standard 3-step recurrence ( which is not true for op in the complex plane in general ) .
the step in the proof of eq .
( [ cd ] ) that breaks down is the cancellation of the summands after multiplying eq .
( [ clagkernel ] ) by @xmath560 .
the most general result corresponding to eq .
( [ prodchar ] ) now contains both products of @xmath520 characteristic polynomials and of @xmath126 products of complex conjugated ones . without loss of generality
let us choose @xmath561 ( @xmath562 can be obtained by complex conjugation ) , to state the result from @xcite : @xmath563 . \label{cprodchar}\ ] ] here the left block of kernels inside the determinant is of size @xmath564 .
the proof is straightforward and uses induction . using this result we are now ready to compute the partition function eq .
( [ zrmtmu ] ) for general @xmath565 . in the quenched case
it is given again by @xmath185 times the products of the norms from eq .
( [ opmu ] ) , as derived in eq .
( [ znorm ] ) . for @xmath88 we can express it in terms of the laguerre polynomials from eq .
( [ opmu ] ) , and obtain an almost identical answer to eq .
( [ zmassive ] ) : @xmath173}(\{m\})=%z_n^{[n_f=0]}\left\langle \prod_{l=1}^n\prod_{f=1}^{n_f}(z_l+m_f^2)\right\rangle_n&= & \frac{(-1)^{nn_f}n!\prod_{j=0}^{n-1}h_j}{\delta_{n_f}(\{m\ } ) } \det_{1\leq j , g\leq n_f}\left[(-1/c)^{n+g-1}(n+g-1)!l^\nu_{n+g-1}(-cm_f^2)\right ] .
\label{zmassivec}\ ] ] as in eq .
( [ zmassive ] ) it is given by a determinant of laguerre polynomials of real arguments which are real .
the only difference to @xmath489 is the scale factor @xmath566 and the @xmath10-dependence inside the norms @xmath567 from eq .
( [ opmu ] ) .
we will not give further details regarding the kernel @xmath568}(z , u^*)$ ] constituting the @xmath62-point correlation functions in eq .
( [ rkc ] ) . following the first line of eq .
( [ kernelmassive ] ) it can be expressed through the quenched average of @xmath569 characteristic polynomials times one conjugated characteristic polynomial , divided by the partition function from eq .
( [ zmassivec ] ) .
the spectral density with @xmath0 masses then follows as in eq .
( [ r1massive ] ) , by multiplying @xmath568}(z , z^*)$ ] with the weight eq .
( [ muweight ] ) , see @xcite for more details . clearly it is not real in general .
we will now very briefly discuss the large-@xmath34 limit and matching with the corresponding echpt .
it turns out that for non - hermitian matrices that are coupled as in eq .
( [ 2matrixrep ] ) the possibility of taking such limits is much richer than in the hermitian ( or maximally non - hermitian ) case .
in addition to the standard macroscopic and microscopic statistics a further class of limits exist , where the matrix @xmath504 becomes almost hermitian , taking the parameter @xmath570 at such a rate that that the following products @xmath571 are constant . here
the scaling of the complex eigenvalues ( and masses ) is unchanged compared to eq .
( [ rmtmass ] ) .
this limit was first observed in the elliptic complex ginibre ensemble already mentioned @xcite and called weakly non - hermitian @xcite .
this is in contrast to the limit at strong non - hermiticity where @xmath572 $ ] is kept fixed and a different scaling of @xmath527 is chosen . here
the statistics of the product of two independent complex matrices with @xmath499 is found , see @xcite and @xcite for details .
looking at our result eq .
( [ zmassivec ] ) it is clear that the limiting partition function is the same as for @xmath489 in eq .
( [ zrmtnf ] ) and is thus independent of the rescaled @xmath573 .
in contrast already the limiting quenched density does depend on @xmath573
. it can be easily obtained from eq .
( [ clagkernel ] ) , replacing the sum by an integral and using eq .
( [ laguerrelim ] ) for the laguerre polynomials , @xmath574\int_0 ^ 1 dtt \exp[-2 t^2\hat{\mu}^2]|j_\nu(\tilde{z}t)|^2\ . \label{cdensity}\ ] ] it is shown in figure [ fig : besselc ] . here
we have again changed to dirac eigenvalues following @xcite .
this result was obtained earlier using replicas and the toda equation in @xcite .
for @xmath575 the integral becomes elementary and matches eq .
( [ besseldensity ] ) , while the prefactors turn into a delta - function .
( 12.2,5.5 ) ) as a function of @xmath576 for @xmath227 ( left ) and @xmath264 ( right ) , at @xmath577 .
the density of real eigenvalues from figure [ fig : besseldensity ] can still be recognised , it now spreads into the complex plane .
because in the limit @xmath578 it will become proportionally to a delta - function @xmath579 times the bessel density the vertical scale is different .
[ fig : besselc],title="fig:",width=264 ] ) as a function of @xmath576 for @xmath227 ( left ) and @xmath264 ( right ) , at @xmath577 .
the density of real eigenvalues from figure [ fig : besseldensity ] can still be recognised , it now spreads into the complex plane .
because in the limit @xmath578 it will become proportionally to a delta - function @xmath579 times the bessel density the vertical scale is different .
[ fig : besselc],title="fig:",width=264 ] ( -5.1,0.8)@xmath265 ( -1.0,1.3)@xmath580 ( -12.1,0.8)@xmath265 ( -8.0,1.3)@xmath580 ( -7.0,4.4)@xmath581 ( -14.0,4.4)@xmath581    the fact that the partition function is @xmath573-independent and the density is not has been called the `` silver blaze problem '' @xcite , alluding to a novel by a.c . doyle .
how can we resolve this puzzle ?
let us look again at chpt , the theory of goldstone bosons . because these are mesons they do nt carry baryon charge and
thus do not feel the effect of adding @xmath10 to the action ( mesons consist of a quark and anti - quark whereas baryons , as the proton , of 3 quarks ) .
we can still make the chpt partition function @xmath10-dependent by adding an integer number @xmath582 of quarks with the opposite chemical potential @xmath583 . as discussed after eq .
( [ dchiralmu ] ) this amounts to add complex conjugated determinants of the dirac operator to the action , leading to @xmath584}(\{m\},\{n\})= z_n \left\langle \prod_{l=1}^n\left(\prod_{i=1}^{n_f}(z_l+m_f)\right)\left(\prod_{j=1}^{n_f^*}(z_l^*+n_f^2)\right)\right\rangle_n\ , \label{zrmtmix}\ ] ] which we can also compute using eq .
( [ cprodchar ] ) .
the corresponding limiting partion function of echpt now does depend on @xmath10 as follows : @xmath585}(\{m\},\{n\})=\int_{u(n_f+n_f^ * ) } du_0\det[u_0]^\nu\exp\left[\frac12 \sigma v\mbox{tr}({m}(u_0+u_0^\dag))-\frac14\mu^2 vf_\pi^2[u_0,b][u_0,b]\right ] .
\label{zechptmu}\ ] ] here @xmath586 contains both kinds of masses and @xmath587 is the metric with the signature of the chemical potentials .
the commutator comes from the fact that in the dirac operator @xmath10 acts as an additional vector potential .
it obviously vanishes at @xmath588 . once again in the limit the partition functions from rmt and echpt agree , when identifying the rescaled masses and chemical potentials as @xmath589 : @xmath590}(\{m\},\{n\})={\cal z}_{\rm echpt}^{[n_f+n_f^*]}(\{\hat{m}\},\{\hat{n}\}).\ ] ] they are given by the determinant from eq .
( [ cprodchar ] ) in terms of the limiting kernel and the limiting polynomials which are again @xmath443 .
the corresponding group integral eq .
( [ zechptmu ] ) was computed independently in @xcite using an explicit parametrisation of @xmath591 , see @xcite for the general @xmath0 case .    what about the equivalence of rmt and echpt for the density , in other words how can we generate the spectral density eq .
( [ cdensity ] ) from a chiral lagrangian ?
although one may still define and generate the resolvent as in eq .
( [ resolvent ] ) , the relation to the spectral density changes compared to eq .
( [ inversion ] ) : @xmath592 due to @xmath593 on @xmath594 .
while it is subtle to extract the non - holomorphic part from the expectation value of the ratio of characteristic polynomials two alternative ways exist to formulate the generation of a two - dimensional delta - function . one way is to use so - called replicas which amounts to introduce @xmath62 extra pairs of complex conjugated quarks in the chiral lagrangian . using the relation to the toda lattice equation , the replica limit can be made well defined as was developed in @xcite .
the generating function for the density with @xmath595 , @xmath596}$ ] on the echpt side thus depends on @xmath573 , which is why the resulting density also depends on it .
the second possibility is to use supersymmetry , which will also involve pairs of complex conjugated bosons @xcite , @xmath597}{\det[(u - d\!\!\!\!/(\mu))(u - d\!\!\!\!/(\mu))^\dag+\kappa^2]}\right\rangle_\nu^{\rm qcd}\right|_{u = z}\right)\ , \label{susyresmu}\ ] ] and both methods have been shown to be equivalent @xcite . in either case
the silver blaze puzzle is resolved .
we would also like to mention the papers @xcite where averages of ratios of characteristic polynomials and their complex conjugates were evaluated in a purely complex op framework .
the most prominent example for such an average is the distribution of the phase of @xmath598 @xmath599}\ , \label{phase}\ ] ] and we refer to @xcite and @xcite for a detailed discussion including its physical interpretation .    when trying to compare the rmt predictions to qcd lattice simulations the presence of the chemical potential also has a virtue , despite the sign problem : it couples to the second lec in the chiral lagrangian , @xmath337 , in eq .
( [ zechptmu ] ) . just as for @xmath26
a fit of the density to lattice data helps to determine @xmath8 , a two - parameter fit at @xmath600 allows to measure @xmath337 in this way . because for imaginary chemical potential @xmath601 the dirac operator spectrum remains real as discussed after eq .
( [ dchiralmu ] ) , this was proposed to determine @xmath337 using standard lattice simulations @xcite .
there is also a corresponding two - matrix model @xcite which in the large-@xmath34 with rescaled @xmath601 as in eq .
( [ weaklim ] ) becomes equal to eq .
( [ zechptmu ] ) , with @xmath602 ( which is of course still convergent ) . turning back to real chemical potential
most comparisons were done within the quenched approximation , starting with @xcite , cf .
@xcite for the influence of topology .
despite the difficulty of defining individual eigenvalue distributions these were compared with in @xcite .
there are many issues that we have not discussed about the application to qcd with chemical potential , apart from methods different from the op approach mentioned above .
this includes the different origin of chiral symmetry breaking at @xmath603 that is strongly related to the oscillatory behaviour of the unquenched density and we refer to @xcite and to @xcite for a review .
thanks to its analytical solution the rmt eq .
( [ zrmtmu ] ) has been used as a testing ground for numerical ways to solve the sign problem .
this has been successful for the rmt itself and is reported in @xcite .
furthermore , the study of the singular values of the qcd dirac operator with chemical potential has been proposed in @xcite , that could also be used to analyse different phases ( not depicted in figure [ fig : phasediag ] ) for larger values of @xmath10 . here
both the lattice realisation and the rmt have many open questions , and a first step in rmt has been taken in @xcite . here
a different kernel is found at the origin , the so - called meijer g - kernel @xcite that generalises the bessel kernel .
the mathematical question of universality in the limit of weak non - hermiticity for non - gaussian rmts is still open at the moment , despite first heuristic attempts @xcite . for a further discussion including all known limiting kernels in this weak non - hermiticity limit resulting from gaussian models we refer to @xcite and references therein .
let us briefly summarise the lessons we have learned in these lectures .
the theory of orthogonal polynomials on the real line and in the complex plane and their properties has been extensively discussed .
a particular focus was put on weights including characteristic polynomials , as this allowed us to determine the dirac operator spectrum of qcd including the contribution of @xmath0 quarks flavours in a random matrix approach .
based on the global symmetries of qcd we have taken a limit to the low energy regime that eventually leads to a random matrix description in a controlled approximation , to which corrections can be computed . even though this is an unphysical limit it can be compared to numerical solutions of qcd on the lattice .
this leads to a determination of the low energy constants of the theory and tests of lattice algorithms to be used then in the physical regime .
we have shown how the introduction of more details in the theory of qcd such as the effect of finite lattice spacing or of chemical potential , that make the eigenvalues complex , modify the spectrum .
and we have mentioned several open problems relevant for rmt and for its application like the computation of new correlation functions or the proof of universality for existing results .
they hopefully contribute to keep this an interesting and lively area of research for both mathematicians and theoretical physicists .
+ acknowledgements : it has been a great pleasure to be back in les houches , now as a lecturer , and i would like to express my sincere gratitude to the organisers for inviting me .
my understanding of this subject would not have been possible without the many collaborations on this topic , and i would like to thank all my coworkers for sharing their insights and compassion .
i am particularly thankful to jac verbaarschot for many valuable comments on a first draft of this manuscript .
furthermore i would like to thank the simons center for geometry and physics at stony brook university ( programme on statistical mechanics and combinatorics 2016 ) for the kind hospitality where part of these lecture notes were written up .
this work was partly supported by dfg grant ak 35/2 - 1 .
f.  basile and g.  akemann , jhep * 0712 * ( 2007 ) 043 [ arxiv:0710.0376 [ hep - th ] ] . c.  lehner and t.  wettig , jhep
* 0911 * ( 2009 ) 005 [ arxiv:0909.1489 [ hep - lat ] ] . g.  akemann , f.  basile and l.  lellouch , jhep * 0812 * ( 2008 ) 069 [ arxiv:0804.3809 [ hep - lat ] ] .