BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Peter J. Forrester

Peter J. Forrester appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPmath.PRcond-mat.stat-mechmath.CAmath.CVnlin.SIComputationhep-thmath.COmath.HOphysics.data-anquant-ph

Catalog footprint

What is connected

67works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

67 published item(s)

preprint2026arXiv

Integrability enabled computations relating to the fixed trace Laguerre ensemble

Studies of density matrices for random quantum states lead naturally to the fixed trace Laguerre ensemble in random matrix theory. Previous studies have uncovered explicit rational function formulas for moments of purity statistic (trace of the squared density matrix), and also a third order linear differential equation satisfied by the eigenvalue density. We further probe the origin of these results from the viewpoint of integrability, which is taken here to mean wider classes of recursions and differential equations, and give extensions. Prominent in our study are first order linear matrix differential equations. One application given is to the derivation of the third order scalar equation for the density. Another is to obtain the explicit rational function formula for the variance of the purity statistic in the $β$ generalised fixed trace Laguerre ensemble. In the original case ($β= 2$), the purity cumulants are expressed in terms of the large argument expansion of a particular $σ$-Painlevé IV transcendent. In a different but related direction, the exact computation of the two-point correlation for the fixed determinant circular unitary ensemble SU$(N)$ is given the Appendix.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Progress on the study of the Ginibre ensembles I: GinUE

The Ginibre unitary ensemble (GinUE) consists of $N \times N$ random matrices with independent complex standard Gaussian entries. This was introduced in 1965 by Ginbre, who showed that the eigenvalues form a determinantal point process with an explicit correlation kernel, and after scaling they are supported on the unit disk with constant density. For some time now it has been appreciated that GinUE has a fundamental place within random matrix theory, both for its applications and for the richness of its theory. Here we review the progress on a number of themes relating to the study of GinUE. These are eigenvalue probability density functions and correlation functions, fluctuation formulas, sum rules and asymptotic behaviours of correlation functions, and normal matrix models. We discuss too applications in quantum many body physics and quantum chaos, and give an account of some statistical properties of the eigenvectors.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Progress on the study of the Ginibre ensembles II: GinOE and GinSE

This is part II of a review relating to the three classes of random non-Hermitian Gaussian matrices introduced by Ginibre in 1965. While part I restricted attention to the GinUE (Ginibre unitary ensemble) case of complex elements, in this part the cases of real elements (GinOE, denoting Ginibre orthogonal ensemble) and quaternion elements represented as $2 \times 2$ complex blocks (GinSE, denoting Ginibre symplectic ensemble) are considered. The eigenvalues of both GinOE and GinSE form Pfaffian point processes, which are more complicated than the determinantal point processes resulting from GinUE. Nevertheless, many of the obstacles that have slowed progress on the development of traditional aspects of the theory have now been overcome, while new theoretical aspects and new applications have been identified. This permits a comprehensive account of themes addressed too in the complex case: eigenvalue probability density functions and correlation functions, limit formulas for correlation functions, fluctuation formulas, sum rules, gap probabilities and eigenvector statistics, among others. Distinct from the complex case is the need to develop a theory of skew orthogonal polynomials corresponding to the skew inner product associated with the Pfaffian. Another distinct theme is the statistics of real eigenvalues, which is unique to GinOE. These appear in a number of applications of the theory, coming from areas as diverse as diffusion processes and persistence in statistical physics, topologically driven parametric energy level crossings for certain quantum dots, and equilibria counting for a system of random nonlinear differential equations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Finite size corrections relating to distributions of the length of longest increasing subsequences

Considered are the large $N$, or large intensity, forms of the distribution of the length of the longest increasing subsequences for various models. Earlier work has established that after centring and scaling, the limit laws for these distributions relate to certain distribution functions at the hard edge known from random matrix theory. By analysing the hard to soft edge transition, we supplement and extend results of Baik and Jenkins for the Hammersley model and symmetrisations, which give that the leading correction is proportional to $z^{-2/3}$, where $z^2$ is the intensity of the Poisson rate, and provides a functional form as derivates of the limit law. Our methods give the functional form both in terms of Fredholm operator theoretic quantities, and in terms of Painlevé transcendents. For random permutations and their symmetrisations, numerical analysis of exact enumerations and simulations gives compelling evidence that the leading corrections are proportional to $N^{-1/3}$, and moreover provides an approximation to their graphical forms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rank $1$ perturbations in random matrix theory -- a review of exact results

A number of random matrix ensembles permitting exact determination of their eigenvalue and eigenvector statistics maintain this property under a rank $1$ perturbation. Considered in this review are the additive rank $1$ perturbation of the Hermitian Gaussian ensembles, the multiplicative rank $1$ perturbation of the Wishart ensembles, and rank $1$ perturbations of Hermitian and unitary matrices giving rise to a two-dimensional support for the eigenvalues. The focus throughout is on exact formulas, which are typically the result of various integrable structures. The simplest is that of a determinantal point process, with others relating to partial differential equations implied by a formulation in terms of certain random tridiagonal matrices. Attention is also given to eigenvector overlaps in the setting of a rank $1$ perturbation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Joint moments of a characteristic polynomial and its derivative for the circular $β$-ensemble

The problem of calculating the scaled limit of the joint moments of the characteristic polynomial, and the derivative of the characteristic polynomial, for matrices from the unitary group with Haar measure first arose in studies relating to the Riemann zeta function in the thesis of Hughes. Subsequently, Winn showed that these joint moments can equivalently be written as the moments for the distribution of the trace in the Cauchy unitary ensemble, and furthermore relate to certain hypergeometric functions based on Schur polynomials, which enabled explicit computations. We give a $β$-generalisation of these results, where now the role of the Schur polynomials is played by the Jack polynomials. This leads to an explicit evaluation of the scaled moments for all $β> 0$, subject to the constraint that a particular parameter therein is equal to a non negative integer. Consideration is also given to the calculation of the moments of the singular statistic $\sum_{j=1}^N 1/x_j$ for the Jacobi $β$-ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotic correlations with corrections for the circular Jacobi $β$-ensemble

Previous works have considered the leading correction term to the scaled limit of various correlation functions and distributions for classical random matrix ensembles and their $β$ generalisations at the hard and soft edge. It has been found that the functional form of this correction is given by a derivative operation applied to the leading term. In the present work we compute the leading correction term of the correlation kernel at the spectrum singularity for the circular Jacobi ensemble with Dyson indices $β= 1,2$ and 4, and also to the spectral density in the corresponding $β$-ensemble with $β$ even. The former requires an analysis involving the Routh-Romanovski polynomials, while the latter is based on multidimensional integral formulas for generalised hypergeometric series based on Jack polynomials. In all cases this correction term is found to be related to the leading term by a derivative operation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computable structural formulas for the distribution of the $β$-Jacobi edge eigenvalues

The Jacobi ensemble is one of the classical ensembles of random matrix theory. Prominent in applications are properties of the eigenvalues at the spectrum edge, specifically the distribution of the largest (e.g. Roy's largest root test in multivariate statistics) and smallest (e.g.~condition numbers of linear systems) eigenvalues. We identify three ranges of parameter values for which the gap probability determining these distributions is a finite sum with respect to particular bases, and moreover make use of a certain differential-difference system fundamental in the theory of the Selberg integral to provide a recursive scheme to compute the corresponding coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linear Differential Equations for the Resolvents of the Classical Matrix Ensembles

The spectral density for random matrix $β$ ensembles can be written in terms of the average of the absolute value of the characteristic polynomial raised to the power of $β$, which for even $β$ is a polynomial of degree $β(N-1)$. In the cases of the classical Gaussian, Laguerre, and Jacobi weights, we show that this polynomial, and moreover the spectral density itself, can be characterised as the solution of a linear differential equation of degree $β+1$. This equation, and its companion for the resolvent, are given explicitly for $β=2$ and $4$ for all three classical cases, and also for $β=6$ in the Gaussian case. Known dualities for the spectral moments relating $β$ to $4/β$ then imply corresponding differential equations in the case $β=1$, and for the Gaussian ensemble, the case $β=2/3$. We apply the differential equations to give a systematic derivation of recurrences satisfied by the spectral moments and by the coefficients of their $1/N$ expansions, along with first-order differential equations for the coefficients of the $1/N$ expansions of the corresponding resolvents. We also present the form of the differential equations when scaled at the hard or soft edges.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rate of convergence at the hard edge for various Pólya ensembles of positive definite matrices

The theory of Pólya ensembles of positive definite random matrices provides structural formulas for the corresponding biorthogonal pair, and correlation kernel, which are well suited to computing the hard edge large $N$ asymptotics. Such an analysis is carried out for products of Laguerre ensembles, the Laguerre Muttalib-Borodin ensemble, and products of Laguerre ensembles and their inverses. The latter includes as a special case the Jacobi unitary ensemble. In each case the hard edge scaled kernel permits an expansion in powers of $1/N$, with the leading term given in a structured form involving the hard edge scaling of the biorthogonal pair. The Laguerre and Jacobi ensembles have the special feature that their hard edge scaled kernel -- the Bessel kernel -- is symmetric and this leads to there being a choice of hard edge scaling variables for which the rate of convergence of the correlation functions is $O(1/N^2)$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Co-rank 1 projections and the randomised Horn problem

Let $\hat{\boldsymbol x}$ be a normalised standard complex Gaussian vector, and project an Hermitian matrix $A$ onto the hyperplane orthogonal to $\hat{\boldsymbol x}$. In a recent paper Faraut [Tunisian J. Math. \textbf{1} (2019), 585--606] has observed that the corresponding eigenvalue PDF has an almost identical structure to the eigenvalue PDF for the rank 1 perturbation $A + b \hat{\boldsymbol x} \hat{\boldsymbol x}^\dagger$, and asks for an explanation. We provide this by way of a common derivation involving the secular equations and associated Jacobians. This applies too in related setting, for example when $\hat{\boldsymbol x}$ is a real Gaussian and $A$ Hermitian, and also in a multiplicative setting $A U B U^\dagger$ where $A, B$ are fixed unitary matrices with $B$ a multiplicative rank 1 deviation from unity, and $U$ is a Haar distributed unitary matrix. Specifically, in each case there is a dual eigenvalue problem giving rise to a PDF of almost identical structure.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Schwinger-Dyson and loop equations for a product of square Ginibre random matrices

In this paper, we study the product of two complex Ginibre matrices and the loop equations satisfied by their resolvents (i.e. the Stieltjes transform of the correlation functions). We obtain using Schwinger-Dyson equation (SDE) techniques the general loop equations satisfied by the resolvents. In order to deal with the product structure of the random matrix of interest, we consider SDEs involving the integral of higher derivatives. One of the advantage of this technique is that it bypasses the reformulation of the problem in terms of singular values. As a byproduct of this study we obtain the large $N$ limit of the Stieltjes transform of the $2$-point correlation function, as well as the first correction to the Stieltjes transform of the density, giving us access to corrections to the smoothed density. In order to pave the way for the establishment of a topological recursion formula we also study the geometry of the corresponding spectral curve. This paper also contains explicit results for different resolvents and their corrections.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Multiplicative Convolution of Real Asymmetric and Real Antisymmetric Matrices

The singular values of products of standard complex Gaussian random matrices, or sub-blocks of Haar distributed unitary matrices, have the property that their probability distribution has an explicit, structured form referred to as a polynomial ensemble. It is furthermore the case that the corresponding bi-orthogonal system can be determined in terms of Meijer G-functions, and the correlation kernel given as an explicit double contour integral. It has recently been shown that the Hermitised product $X_M \cdots X_2 X_1A X_1^T X_2^T \cdots X_M^T$, where each $X_i$ is a standard real complex Gaussian matrix, and $A$ is real anti-symmetric shares exhibits analogous properties. Here we use the theory of spherical functions and transforms to present a theory which, for even dimensions, includes these properties of the latter product as a special case. As an example we show that the theory also allows for a treatment of this class of Hermitised product when the $X_i$ are chosen as sub-blocks of Haar distributed real orthogonal matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A. Hurwitz and the origins of random matrix theory in mathematics

The purpose of this article is to put forward the claim that Hurwitz's paper "Uber die Erzeugung der Invarianten durch Integration." [Gott. Nachrichten (1897), 71-90] should be regarded as the origin of random matrix theory in mathematics. Here Hurwitz introduced and developed the notion of an invariant measure for the matrix groups $SO(N)$ and $U(N)$. He also specified a calculus from which the explicit form of these measures could be computed in terms of an appropriate parametrisation - Hurwitz chose to use Euler angles. This enabled him to define and compute invariant group integrals over $SO(N)$ and $U(N)$. His main result can be interpreted probabilistically: the Euler angles of a uniformly distributed matrix are independent with beta distributions (and conversely). We use this interpretation to give some new probability results. How Hurwitz's ideas and methods show themselves in the subsequent work of Weyl, Dyson and others on foundational studies in random matrix theory is detailed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The scaling of boson sampling experiments

Boson sampling is the problem of generating a quantum bit stream whose average is the permanent of a $n\times n$ matrix. The bitstream is created as the output of a prototype quantum computing device with $n$ input photons. It is a fundamental challenge to verify boson sampling, and the question of how output count rates scale with matrix size $n$ is crucial. Here we apply results from random matrix theory to establish scaling laws for average count rates in boson sampling experiments with arbitrary inputs and losses. The results show that, even with losses included, verification of nonclassical behaviour at large $n$ values is indeed possible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Volumes for ${\rm SL}_N(\mathbb R)$, the Selberg integral and random lattices

There is a natural left and right invariant Haar measure associated with the matrix groups GL${}_N(\mathbb R)$ and SL${}_N(\mathbb R)$ due to Siegel. For the associated volume to be finite it is necessary to truncate the groups by imposing a bound on the norm, or in the case of SL${}_N(\mathbb R)$, by restricting to a fundamental domain. We compute the asymptotic volumes associated with the Haar measure for GL${}_N(\mathbb R)$ and SL${}_N(\mathbb R)$ matrices in the case of that the operator norm lies between $R_1$ and $1/R_2$ in the former, and this norm, or alternatively the 2-norm, is bounded by $R$ in the latter. By a result of Duke, Rundnick and Sarnak, such asymptotic formulas in the case of SL${}_N(\mathbb R)$ imply an asymptotic counting formula for matrices in SL${}_N(\mathbb Z)$. We discuss too the sampling of SL${}_N(\mathbb R)$ matrices from the truncated sets. By then using lattice reduction to a fundamental domain, we obtain histograms approximating the probability density functions of the lengths and pairwise angles of shortest length bases vectors in the case $N=2$ and 3, or equivalently of shortest linearly independent vectors in the corresponding random lattice. In the case $N=2$ these distributions are evaluated explicitly.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analogies between random matrix ensembles and the one-component plasma in two-dimensions

The eigenvalue PDF for some well known classes of non-Hermitian random matrices --- the complex Ginibre ensemble for example --- can be interpreted as the Boltzmann factor for one-component plasma systems in two-dimensional domains. We address this theme in a systematic fashion, identifying the plasma system for the Ginibre ensemble of non-Hermitian Gaussian random matrices $G$, the spherical ensemble of the product of an inverse Ginibre matrix and a Ginibre matrix $G_1^{-1} G_2$, and the ensemble formed by truncating unitary matrices, as well as for products of such matrices. We do this when each has either real, complex or real quaternion elements. One consequence of this analogy is that the leading form of the eigenvalue density follows as a corollary. Another is that the eigenvalue correlations must obey sum rules known to characterise the plasma system, and this leads us to a exhibit an integral identity satisfied by the two-particle correlation for real quaternion matrices in the neighbourhood of the real axis. Further random matrix ensembles investigated from this viewpoint are self dual non-Hermitian matrices, in which a previous study has related to the one-component plasma system in a disk at inverse temperature $β= 4$, and the ensemble formed by the single row and column of quaternion elements from a member of the circular symplectic ensemble.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymptotics of finite system Lyapunov exponents for some random matrix ensembles

For products $P_N$ of $N$ random matrices of size $d \times d$, there is a natural notion of finite $N$ Lyapunov exponents $\{μ_i\}_{i=1}^d$. In the case of standard Gaussian random matrices with real, complex or real quaternion elements, and extended to the general variance case for $μ_1$, methods known for the computation of $\lim_{N \to \infty} \langle μ_i \rangle$ are used to compute the large $N$ form of the variances of the exponents. Analogous calculations are performed in the case that the matrices making up $P_N$ are products of sub-blocks of random unitary matrices with Haar measure. Furthermore, we make some remarks relating to the coincidence of the Lyapunov exponents and the stability exponents relating to the eigenvalues of $P_N$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Diffusion processes and the asymptotic bulk gap probability for the real Ginibre ensemble

It is known that the bulk scaling limit of the real eigenvalues for the real Ginibre ensemble is equal in distribution to the rescaled $t \to \infty$ limit of the annihilation process $A + A \to \emptyset$. Furthermore, deleting each particle at random in the rescaled $t \to \infty$ limit of the coalescence process $A + A \to A$, a process equal in distribution to the annihilation process results. We use these inter-relationships to deduce from the existing literature the asymptotic small and large distance form of the gap probability for the real Ginibre ensemble. In particular, the leading form of the latter is shown to be equal to $\exp(-(ζ(3/2)/(2 \sqrt{2 π}))s)$, where $s$ denotes the gap size and $ζ(z)$ denotes the Riemann zeta function. It is shown how this can be rigorously established using an asymptotic formula for matrix Fredholm operators. A determinant formula is derived for the gap probability in the finite $N$ case, and this is used to illustrate the asymptotic formulas against numerical computations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite size corrections in random matrix theory and Odlyzko's data set for the Riemann zeros

Odlyzko has computed a data set listing more than $10^9$ successive Riemann zeros, starting at a zero number beyond $10^{23}$. The data set relates to random matrix theory since, according to the Montgomery-Odlyzko law, the statistical properties of the large Riemann zeros agree with the statistical properties of the eigenvalues of large random Hermitian matrices. Moreover, Keating and Snaith, and then Bogomolny and collaborators, have used $N \times N$ random unitary matrices to analyse deviations from this law. We contribute to this line of study in two ways. First, we point out that a natural process to apply to the data set is to thin it by deleting each member independently with some specified probability, and we proceed to compute empirical two-point correlation functions and nearest neighbour spacings in this setting. Second, we show how to characterise the order $1/N^2$ correction term to the spacing distribution for random unitary matrices in terms of a second order differential equation with coefficients that are Painlevé transcendents, and where the thinning parameter appears only in the boundary condition. This equation can be solved numerically using a power series method. Comparison with the Riemann zero data shows accurate agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Singular Values and Evenness Symmetry in Random Matrix Theory

Complex Hermitian random matrices with a unitary symmetry can be distinguished by a weight function. When this is even, it is a known result that the distribution of the singular values can be decomposed as the superposition of two independent eigenvalue sequences distributed according to particular matrix ensembles with chiral unitary symmetry. We give decompositions of the distribution of singular values, and the decimation of the singular values --- whereby only even, or odd, labels are observed --- for real symmetric random matrices with an orthogonal symmetry, and even weight. This requires further specifying the functional form of the weight to one of three types --- Gauss, symmetric Jacobi or Cauchy. Inter-relations between gap probabilities with orthogonal and unitary symmetry follow as a corollary. The Gauss case has appeared in a recent work of Bornemann and La Croix. The Cauchy case, when appropriately specialised and upon stereographic projection, gives decompositions for the analogue of the singular values for the circular unitary and circular orthogonal ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Singular values for products of complex Ginibre matrices with a source: hard edge limit and phase transition

The singular values squared of the random matrix product $Y = G_r G_{r-1} \cdots G_1 (G_0 + A)$, where each $G_j$ is a rectangular standard complex Gaussian matrix while $A$ is non-random, are shown to be a determinantal point process with correlation kernel given by a double contour integral. When all but finitely many eigenvalues of $A^*A$ are equal to $bN$, the kernel is shown to admit a well-defined hard edge scaling, in which case a critical value is established and a phase transition phenomenon is observed. More specifically, the limiting kernel in the subcritical regime of $0<b<1$ is independent of $b$, and is in fact the same as that known for the case $b=0$ due to Kuijlaars and Zhang. The critical regime of $b=1$ allows for a double scaling limit by choosing $b = (1-τ/\sqrt{N})^{-1}$, and for this the critical kernel and outlier phenomenon are established. In the simplest case $r=0$, which is closely related to non-intersecting squared Bessel paths, a distribution corresponding to the finite shifted mean LUE is proven to be the scaling limit in the supercritical regime of $b>1$ with two distinct scaling rates. Similar results also hold true for the random matrix product $T_r T_{r-1} \cdots T_1 (G_0 + A)$, with each $T_j$ being a truncated unitary matrix.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Eigenvalue statistics for product complex Wishart matrices

The eigenvalue statistics for complex $N \times N$ Wishart matrices $X_{r,s}^\dagger X_{r,s}$, where $ X_{r,s}$ is equal to the product of $r$ complex Gaussian matrices, and the inverse of $s$ complex Gaussian matrices, are considered. In the case $r=s$ the exact form of the global density is computed. The averaged characteristic polynomial for the corresponding generalized eigenvalue problem is calculated in terms of a particular generalized hypergeometric function ${}_{s+1} F_r$. For finite $N$ the eigenvalue probability density function is computed, and is shown to be an example of a biorthogonal ensemble. A double contour integral form of the corresponding correlation kernel is derived, which allows the hard edge scaled limit to be computed. The limiting kernel is given in terms of certain Meijer G-functions, and is identical to that found in the recent work of Kuijlaars and Zhang in the case $s=0$. Properties of the kernel and corresponding correlation functions are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equilibrium problems for Raney densities

The Raney numbers are a class of combinatorial numbers generalising the Fuss--Catalan numbers. They are indexed by a pair of positive real numbers $(p,r)$ with $p>1$ and $0 < r \le p$, and form the moments of a probability density function. For certain $(p,r)$ the latter has the interpretation as the density of squared singular values for certain random matrix ensembles, and in this context equilibrium problems characterising the Raney densities for $(p,r) = (θ+1,1)$ and $(θ/2+1,1/2)$ have recently been proposed. Using two different techniques --- one based on the Wiener--Hopf method for the solution of integral equations and the other on an analysis of the algebraic equation satisfied by the Green's function --- we establish the validity of the equilibrium problems for general $θ> 0$ and similarly use both methods to identify the equilibrium problem for $(p,r) = (θ/q+1,1/q)$, $θ> 0$ and $q \in \mathbb Z^+$. The Wiener--Hopf method is used to extend the latter to parameters $(p,r) = (θ/q + 1, m+ 1/q)$ for $m$ a non-negative integer, and also to identify the equilibrium problem for a family of densities with moments given by certain binomial coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local Central Limit Theorem for Determinantal Point Processes

We prove a local central limit theorem (LCLT) for the number of points $N(J)$ in a region $J$ in $\mathbb R^d$ specified by a determinantal point process with an Hermitian kernel. The only assumption is that the variance of $N(J)$ tends to infinity as $|J| \to \infty$. This extends a previous result giving a weaker central limit theorem (CLT) for these systems. Our result relies on the fact that the Lee-Yang zeros of the generating function for $\{E(k;J)\}$ --- the probabilities of there being exactly $k$ points in $J$ --- all lie on the negative real $z$-axis. In particular, the result applies to the scaled bulk eigenvalue distribution for the Gaussian Unitary Ensemble (GUE) and that of the Ginibre ensemble. For the GUE we can also treat the properly scaled edge eigenvalue distribution. Using identities between gap probabilities, the LCLT can be extended to bulk eigenvalues of the Gaussian Symplectic Ensemble (GSE). A LCLT is also established for the probability density function of the $k$-th largest eigenvalue at the soft edge, and of the spacing between $k$-th neigbors in the bulk.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probability densities and distributions for spiked and general variance Wishart $β$-ensembles

A Wishart matrix is said to be spiked when the underlying covariance matrix has a single eigenvalue $b$ different from unity. As $b$ increases through $b=2$, a gap forms from the largest eigenvalue to the rest of the spectrum, and with $b-2$ of order $N^{-1/3}$ the scaled largest eigenvalues form a well defined parameter dependent state. Recent works by Bloemendal and Virág [BV], and Mo, have quantified this parameter dependent state for real Wishart matrices from different viewpoints, and the former authors have done similarly for the spiked Wishart $β$-ensemble. The latter is defined in terms of certain random bidiagonal matrices. We use a recursive structure to give an alternative construction of the spiked and more generally the general variance Wishart $β$-ensemble, and we give the exact form of the joint eigenvalue PDF for the two matrices in the recurrence. In the case of real quaternion Wishart matrices ($β= 4$) the latter is recognised as having appeared in earlier studies on symmetrized last passage percolation, allowing the exact form of the scaled distribution of the largest eigenvalue to be given. This extends and simplifies earlier work of Wang, and is an alternative derivation to a result in [BV]. We also use the construction of the spiked Wishart $β$-ensemble from [BV] to give a simple derivation of the explicit form of the eigenvalue PDF.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Raney distributions and random matrix theory

Recent works have shown that the family of probability distributions with moments given by the Fuss-Catalan numbers permit a simple parameterized form for their density. We extend this result to the Raney distribution which by definition has its moments given by a generalization of the Fuss-Catalan numbers. Such computations begin with an algebraic equation satisfied by the Stieltjes transform, which we show can be derived from the linear differential equation satisfied by the characteristic polynomial of random matrix realizations of the Raney distribution. For the Fuss-Catalan distribution, an equilibrium problem characterizing the density is identified. The Stieltjes transform for the limiting spectral density of the singular values squared of the matrix product formed from $q$ inverse standard Gaussian matrices, and $s$ standard Gaussian matrices, is shown to satisfy a variant of the algebraic equation relating to the Raney distribution. Supported on $(0,\infty)$, we show that it too permits a simple functional form upon the introduction of an appropriate choice of parameterisation. As an application, the leading asymptotic form of the density as the endpoints of the support are approached is computed, and is shown to have some universal features.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The $q$-Dixon--Anderson integral and multi-dimensional $_1ψ_1$ summations

The Dixon--Anderson integral is a multi-dimensional integral evaluation fundamental to the theory of the Selberg integral. The $_1ψ_1$ summation is a bilateral generalization of the $q$-binomial theorem. It is shown that a $q$-generalization of the Dixon--Anderson integral, due to Evans, and multi-dimensional generalizations of the $_1ψ_1$ summation, due to Milne and Gustafson, can be viewed as having a common origin in the theory of $q$-difference equations as expounded by Aomoto. Each is shown to be determined by a $q$-difference equation of rank one, and a certain asymptotic behavior. In calculating the latter, essential use is made of the concepts of truncation, regularization and connection formulae.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Golden-Thompson inequality --- historical aspects and random matrix applications

The Golden-Thompson inequality, ${\rm Tr} \, (e^{A + B}) \le {\rm Tr} \, (e^A e^B)$ for $A,B$ Hermitian matrices, appeared in independent works by Golden and Thompson published in 1965. Both of these were motivated by considerations in statistical mechanics. In recent years the Golden-Thompson inequality has found applications to random matrix theory. In this survey article we detail some historical aspects relating to Thompson's work, giving in particular an hitherto unpublished proof due to Dyson, and correspondence with Pólya. We show too how the $2 \times 2$ case relates to hyperbolic geometry, and how the original inequality holds true with the trace operation replaced by any unitarily invariant norm. In relation to the random matrix applications, we review its use in the derivation of concentration type lemmas for sums of random matrices due to Ahlswede-Winter, and Oliveira, generalizing various classical results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Asymptotics of spacing distributions 50 years later

In 1962 Dyson used a physically based, macroscopic argument to deduce the first two terms of the large spacing asymptotic expansion of the gap probability for the bulk state of random matrix ensembles with symmetry parameter β. In the ensuing years, the question of asymptotic expansions of spacing distributions in random matrix theory has shown itself to have a rich mathematical content. As well as presenting the main known formulas, we give an account of the mathematical methods used for their proofs, and provide some new formulas. We also provide a high precision numerical computation of one of the spacing probabilities to illustrate the accuracy of the corresponding asymptotics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Probability of all eigenvalues real for products of standard Gaussian matrices

With $\{X_i\}$ independent $N \times N$ standard Gaussian random matrices, the probability $p_{N,N}^{P_m}$ that all eigenvalues are real for the matrix product $P_m = X_m X_{m-1} \cdots X_1$ is expressed in terms of an $N/2 \times N/2$ ($N$ even) and $(N+1)/2 \times (N+1)/2$ ($N$ odd) determinant. The entries of the determinant are certain Meijer $G$-functions. In the case $m=2$ high precision computation indicates that the entries are rational multiples of $π^2$, with the denominator a power of 2, and that to leading order in $N$ $p_{N,N}^{P_m}$ decays as $(π/4)^{N^2/2}$. We are able to show that for general $m$ and large $N$, $p_{N,N}^{P_m} \sim b_m^{N^2}$ with an explicit $b_m$. An analytic demonstration that $p_{N,N}^{P_m} \to 1$ as $m \to \infty$ is given.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ramanujan's $_1ψ_1$ summation theorem --- perspective, announcement of bilateral $q$-Dixon--Anderson and $q$-Selberg integral extensions, and context

The Ramanujan $_1ψ_1$ summation theorem in studied from the perspective of $q$-Jackson integrals, $q$-difference equations and connection formulas. This is an approach which has previously been shown to yield Bailey's very-well-poised $_6ψ_6$ summation. Bilateral Jackson integral generalizations of the Dixon--Anderson and Selberg integrals relating to the type $A$ root system are identified as natural candidates for multidimensional generalizations of the Ramanujan $_1ψ_1$ summation theorem. New results of this type are announced, and furthermore they are put into context by reviewing from previous literature explicit product formulas for Jackson integrals relating to other roots systems obtained from the same perspective.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Skew orthogonal polynomials for the real and quaternion real Ginibre ensembles and generalizations

There are some distinguished ensembles of non-Hermitian random matrices for which the joint PDF can be written down explicitly, is unchanged by rotations, and furthermore which have the property that the eigenvalues form a Pfaffian point process. For these ensembles, in which the elements of the matrices are either real, or real quaternion, the kernel of the Pfaffian is completely determined by certain skew orthogonal polynomials, which permit an expression in terms of averages over the characteristic polynomial, and the characteristic polynomial multiplied by the trace. We use Schur polynomial theory, knowledge of the value of a Schur polynomial averaged against real, and real quaternion Gaussian matrices, and the Selberg integral to evaluate these averages.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotic forms for hard and soft edge general $β$ conditional gap probabilities

An infinite log-gas formalism, due to Dyson, and independently Fogler and Shklovskii, is applied to the computation of conditioned gap probabilities at the hard and soft edges of random matrix $β$-ensembles. The conditioning is that there are $n$ eigenvalues in the gap, with $n \ll |t|$, $t$ denoting the end point of the gap. It is found that the entropy term in the formalism must be replaced by a term involving the potential drop to obtain results consistent with known asymptotic expansions in the case $n=0$. With this modification made for general $n$, the derived expansions - which are for the logarithm of the gap probabilities - are conjectured to be correct up to and including terms O$(\log|t|)$. They are shown to satisfy various consistency conditions, including an asymptotic duality formula relating $β$ to $4/β$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotics of spacing distributions at the hard edge for $β$-ensembles

In a previous work [J. Math. Phys. {\bf 35} (1994), 2539--2551], generalized hypergeometric functions have been used to a give a rigorous derivation of the large $s$ asymptotic form of the general $β> 0$ gap probability $E_β^{\rm hard}(0;(0,s);βa/2)$, provided both $βa /2 \in \mathbb Z_\ge 0$ and $2/β\in \mathbb Z^+$. It shown how the details of this method can be extended to remove the requirement that $2/β\in \mathbb Z^+$. Furthermore, a large deviation formula for the gap probability $E_β(n;(0,x);{\rm ME}_{β,N}(λ^{a β/2} e^{βN λ/2}))$ is deduced by writing it in terms of the charateristic function of a certain linear statistic. By scaling $x = s/(4N)^2$ and taking $N \to \infty$, this is shown to reproduce a recent conjectured formula for $E_β^{\rm hard}(n;(0,s);βa/2)$, $βa /2 \in \mathbb Z_{\ge 0}$, and moreover to give a prediction without the latter restriction. This extended formula, which for the constant term involves the Barnes double gamma function, is shown to satisfy an asymptotic functional equation relating the gap probability with parameters $(β,n,a)$, to a gap probability with parameters $(4/β,n',a')$, where $n'=β(n+1)/2-1$, $a'=β(a-2)/2+2$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Expanded Vandermonde powers and sum rules for the two-dimensional one-component plasma

The two-dimensional one-component plasma (2dOCP) is a system of $N$ mobile particles of the same charge $q$ on a surface with a neutralising background. The Boltzmann factor of the 2dOCP at temperature $T$ can be expressed as a Vandermonde determinant to the power $Γ=q^{2}/(k_B T)$. Recent advances in the theory of symmetric and anti-symmetric Jack polymonials provide an efficient way to expand this power of the Vandermonde in their monomial basis, allowing the computation of several thermodynamic and structural properties of the 2dOCP for $N$ values up to 14 and $Γ$ equal to 4, 6 and 8. In this work, we explore two applications of this formalism to study the moments of the pair correlation function of the 2dOCP on a sphere, and the distribution of radial linear statistics of the 2dOCP in the plane.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Large deviation eigenvalue density for the soft edge Laguerre and Jacobi $β$-ensembles

We analyze the eigenvalue density for the Laguerre and Jacobi $β$-ensembles in the cases that the corresponding exponents are extensive. In particular, we obtain the asymptotic expansion up to terms $o(1)$, in the large deviation regime outside the limiting interval of support. As found in recent studies of the large deviation density for the Gaussian $β$-ensemble, and Laguerre $β$-ensemble with fixed exponent, there is a scaling from this asymptotic expansion to the right tail asymptotics for the distribution of the largest eigenvalue at the soft edge.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lyapunov exponents for products of complex Gaussian random matrices

The exact value of the Lyapunov exponents for the random matrix product $P_N = A_N A_{N-1}...A_1$ with each $A_i = Σ^{1/2} G_i^{\rm c}$, where $Σ$ is a fixed $d \times d$ positive definite matrix and $G_i^{\rm c}$ a $d \times d$ complex Gaussian matrix with entries standard complex normals, are calculated. Also obtained is an exact expression for the sum of the Lyapunov exponents in both the complex and real cases, and the Lyapunov exponents for diffusing complex matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Painleve II in random matrix theory and related fields

We review some occurrences of Painlevé II transcendents in the study of two-dimensional Yang-Mills theory, fluctuation formulas for growth models, and as distribution functions within random matrix theory. We first discuss settings in which the parameter $α$ in the Painlevé equation is zero, and the boundary condition is that of the Hasting-MacLeod solution. As well as expressions involving the Painlevé transcendent itself, one encounters the sigma form of the Painlevé II equation, and Lax pair equations in which the Painlevé transcendent occurs as coefficients. We then consider settings which give rise to general $α$ Painlevé II transcendents. In a particular random matrix setting, new results for the corresponding boundary conditions in the cases $α= \pm 1/2$, 1 and 2 are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reunion probability of N vicious walkers: typical and large fluctuations for large N

We consider three different models of N non-intersecting Brownian motions on a line segment [0,L] with absorbing (model A), periodic (model B) and reflecting (model C) boundary conditions. In these three cases we study a properly normalized reunion probability, which, in model A, can also be interpreted as the maximal height of N non-intersecting Brownian excursions on the unit time interval. We provide a detailed derivation of the exact formula for these reunion probabilities for finite N using a Fermionic path integral technique. We then analyse the asymptotic behavior of this reunion probability for large N using two complementary techniques: (i) a saddle point analysis of the underlying Coulomb gas and (ii) orthogonal polynomial method. These two methods are complementary in the sense that they work in two different regimes, respectively for L\ll O(\sqrt{N}) and L\geq O(\sqrt{N}). A striking feature of the large N limit of the reunion probability in the three models is that it exhibits a third-order phase transition when the system size L crosses a critical value L=L_c(N)\sim \sqrt{N}. This transition is akin to the Douglas-Kazakov transition in two-dimensional continuum Yang-Mills theory. While the central part of the reunion probability, for L \sim L_c(N), is described in terms of the Tracy-Widom distributions (associated to GOE and GUE depending on the model), the emphasis of the present study is on the large deviations of these reunion probabilities, both in the right [L \gg L_c(N)] and the left [L \ll L_c(N)] tails. In particular, for model B, we find that the matching between the different regimes corresponding to typical L \sim L_c(N) and atypical fluctuations in the right tail L \gg L_c(N) is rather unconventional, compared to the usual behavior found for the distribution of the largest eigenvalue of GUE random matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The averaged characteristic polynomial for the Gaussian and chiral Gaussian ensembles with a source

In classical random matrix theory the Gaussian and chiral Gaussian random matrix models with a source are realized as shifted mean Gaussian, and chiral Gaussian, random matrices with real $(β= 1)$, complex ($β= 2)$ and real quaternion $(β= 4$) elements. We use the Dyson Brownian motion model to give a meaning for general $β> 0$. In the Gaussian case a further construction valid for $β> 0$ is given, as the eigenvalue PDF of a recursively defined random matrix ensemble. In the case of real or complex elements, a combinatorial argument is used to compute the averaged characteristic polynomial. The resulting functional forms are shown to be a special cases of duality formulas due to Desrosiers. New derivations of the general case of Desrosiers' dualities are given. A soft edge scaling limit of the averaged characteristic polynomial is identified, and an explicit evaluation in terms of so-called incomplete Airy functions is obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fluctuation universality for a class of directed solid-on-solid models

Our interest is in a class of directed solid-on-solid models, which may be regarded as continuum versions of boxed plane partitions. In the case that the heights are chosen from a uniform distribution, the joint PDF of the heights is the same as that for the positions in a finitized bead process recently introduced by the authors and Nordenstam. We use knowledge of the correlation functions for the latter to show that upon a certain scaling the fluctuations of the heights along the back row of the solid-on-solid model are given by the Airy process from random matrix theory, as is the case for boxed plane partitions. Moreover, we show that this limiting distribution remains true if instead of the uniform distribution, the heights are sampled from a general absolutely continuous distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-intersecting Brownian walkers and Yang-Mills theory on the sphere

We study a system of N non-intersecting Brownian motions on a line segment [0,L] with periodic, absorbing and reflecting boundary conditions. We show that the normalized reunion probabilities of these Brownian motions in the three models can be mapped to the partition function of two-dimensional continuum Yang-Mills theory on a sphere respectively with gauge groups U(N), Sp(2N) and SO(2N). Consequently, we show that in each of these Brownian motion models, as one varies the system size L, a third order phase transition occurs at a critical value L=L_c(N)\sim \sqrt{N} in the large N limit. Close to the critical point, the reunion probability, properly centered and scaled, is identical to the Tracy-Widom distribution describing the probability distribution of the largest eigenvalue of a random matrix. For the periodic case we obtain the Tracy-Widom distribution corresponding to the GUE random matrices, while for the absorbing and reflecting cases we get the Tracy-Widom distribution corresponding to GOE random matrices. In the absorbing case, the reunion probability is also identified as the maximal height of N non-intersecting Brownian excursions ("watermelons" with a wall) whose distribution in the asymptotic scaling limit is then described by GOE Tracy-Widom law. In addition, large deviation formulas for the maximum height are also computed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

One-component plasma on a spherical annulus and a random matrix ensemble

The two-dimensional one-component plasma at the special coupling β= 2 is known to be exactly solvable, for its free energy and all of its correlations, on a variety of surfaces and with various boundary conditions. Here we study this system confined to a spherical annulus with soft wall boundary conditions, paying special attention to the resulting asymptotic forms from the viewpoint of expected general properties of the two-dimensional plasma. Our study is motivated by the realization of the Boltzmann factor for the plasma system with β= 2, after stereographic projection from the sphere to the complex plane, by a certain random matrix ensemble constructed out of complex Gaussian and Haar distributed unitary matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pfaffian point process for the Gaussian real generalised eigenvalue problem

The generalised eigenvalues for a pair of $N\times N$ matrices $(X_1,X_2)$ are defined as the solutions of the equation $\det (X_1-λX_2)=0$, or equivalently, for $X_2$ invertible, as the eigenvalues of $X_2^{-1}X_1$. We consider Gaussian real matrices $X_1,X_2$, for which the generalised eigenvalues have the rotational invariance of the half-sphere, or after a fractional linear transformation, the rotational invariance of the unit disk. In these latter variables we calculate the joint eigenvalue probability density function, the probability $p_{N,k}$ of finding $k$ real eigenvalues, the densities of real and complex eigenvalues (the latter being related to an average over characteristic polynomials), and give an explicit Pfaffian formula for the higher correlation functions $ρ_{(k_1,k_2)}$. A limit theorem for $p_{N,k}$ is proved, and the scaled form of $ρ_{(k_1,k_2)}$ is shown to be identical to the analogous limit for the correlations of the eigenvalues of real Gaussian matrices. We show that these correlations satisfy sum rules characteristic of the underlying two-component Coulomb gas.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spectral density asymptotics for Gaussian and Laguerre $β$-ensembles in the exponentially small region

The first two terms in the large $N$ asymptotic expansion of the $β$ moment of the characteristic polynomial for the Gaussian and Laguerre $β$-ensembles are calculated. This is used to compute the asymptotic expansion of the spectral density in these ensembles, in the exponentially small region outside the leading support, up to terms $o(1)$ . The leading form of the right tail of the distribution of the largest eigenvalue is given by the density in this regime. It is demonstrated that there is a scaling from this, to the right tail asymptotics for the distribution of the largest eigenvalue at the soft edge.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Finitization of the Bead Process

The bead process is the particle system defined on parallel lines, with underlying measure giving constant weight to all configurations in which particles on neighbouring lines interlace, and zero weight otherwise. Motivated by the statistical mechanical model of the tiling of an $abc$-hexagon by three species of rhombi, a finitized version of the bead process is defined. The corresponding joint distribution can be realized as an eigenvalue probability density function for a sequence of random matrices. The finitized bead process is determinantal, and we give the correlation kernel in terms of Jacobi polynomials. Two scaling limits are considered: a global limit in which the spacing between lines goes to zero, and a certain bulk scaling limit. In the global limit the shape of the support of the particles is determined, while in the bulk scaling limit the bead process kernel of Boutillier is reclaimed, after approriate identification of the anisotropy parameter therein.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Fuchsian matrix differential equation for Selberg correlation integrals

We characterize averages of $\prod_{l=1}^N|x - t_l|^{α- 1}$ with respect to the Selberg density, further contrained so that $t_l \in [0,x]$ $(l=1,...,q)$ and $t_l \in [x,1]$ $(l=q+1,...,N)$, in terms of a basis of solutions of a particular Fuchsian matrix differential equation. By making use of the Dotsenko-Fateev integrals, the explicit form of the connection matrix from the Frobenius type power series basis to this basis is calculated, thus allowing us to explicitly compute coefficients in the power series expansion of the averages. From these we are able to compute power series for the marginal distributions of the $t_j$ $(j=1,...,N)$. In the case $q=0$ and $α< 1$ we compute the explicit leading order term in the $x \to 0$ asymptotic expansion, which is of interest to the study of an effect known as singularity dominated strong fluctuations. In the case $q=0$ and $α\in \mathbb Z^+$, and with the absolute values removed, the average is a polynomial, and we demonstrate that its zeros are highly structured.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A generalized plasma and interpolation between classical random matrix ensembles

The eigenvalue probability density functions of the classical random matrix ensembles have a well known analogy with the one component log-gas at the special couplings β= 1,2 and 4. It has been known for some time that there is an exactly solvable two-component log-potential plasma which interpolates between the β=1 and 4 circular ensemble, and an exactly solvable two-component generalized plasma which interpolates between β= 2 and 4 circular ensemble. We extend known exact results relating to the latter --- for the free energy and one and two-point correlations --- by giving the general (k_1+k_2)-point correlation function in a Pfaffian form. Crucial to our working is an identity which expresses the Vandermonde determinant in terms of a Pfaffian. The exact evaluation of the general correlation is used to exhibit a perfect screening sum rule.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Connection preserving deformations and $q$-semi-classical orthogonal polynomials

We present a framework for the study of $q$-difference equations satisfied by $q$-semi-classical orthogonal systems. As an example, we identify the $q$-difference equation satisfied by a deformed version of the little $q$-Jacobi polynomials as a gauge transformation of a special case of the associated linear problem for $q$-$\mathrm{P}_{\mathrm{VI}}$. We obtain a parameterization of the associated linear problem in terms of orthogonal polynomial variables and find the relation between this parameterization and that of Jimbo and Sakai.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Difference system for Selberg correlation integrals

The Selberg correlation integrals are averages of the products $\prod_{s=1}^m\prod_{l=1}^n (x_s - z_l)^{μ_s}$ with respect to the Selberg density. Our interest is in the case $m=1$, $μ_1 = μ$, when this corresponds to the $μ$-th moment of the corresponding characteristic polynomial. We give the explicit form of a $(n+1) \times (n+1)$ matrix linear difference system in the variable $μ$ which determines the average, and we give the Gauss decomposition of the corresponding $(n+1) \times (n+1)$ matrix. For $μ$ a positive integer the difference system can be used to efficiently compute the power series defined by this average.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interlaced particle systems and tilings of the Aztec diamond

Motivated by the problem of domino tilings of the Aztec diamond, a weighted particle system is defined on $N$ lines, with line $j$ containing $j$ particles. The particles are restricted to lattice points from 0 to $N$, and particles on successive lines are subject to an interlacing constraint. It is shown that marginal distributions for this particle system can be computed exactly. This in turn is used to give unified derivations of a number of fundamental properties of the tiling problem, for example the evaluation of the number of distinct configurations and the relation to the GUE minor process. An interlaced particle system associated with the domino tiling of a certain half Aztec diamond is similarly defined and analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Jack polynomial fractional quantum Hall states and their generalizations

In the the study of fractional quantum Hall states, a certain clustering condition involving up to four integers has been identified. We give a simple proof that particular Jack polynomials with $α= - (r-1)/(k+1)$, $(r-1)$ and $(k+1)$ relatively prime, and with partition given in terms of its frequencies by $[n_00^{(r-1)s}k 0 ^{r-1}k 0 ^{r-1}k \cdots 0 ^{r-1} m]$ satisfy this clustering condition. Our proof makes essential use of the fact that these Jack polynomials are translationally invariant. We also consider nonsymmetric Jack polynomials, symmetric and nonsymmetric generalized Hermite and Laguerre polynomials, and Macdonald polynomials from the viewpoint of the clustering.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The limiting Kac random polynomial and truncated random orthogonal matrices

An exact calculation of the eigenvalue statistics of truncated random Haar distributed real orthogonal matrices has recently been carried out by Khoruzhenko, Sommers and Zyczkowski. We further develop this calculation, and use it to deduce a Pfaffian form of the correlations for the zeros of the limiting Kac random polynomial. This contrasts with the forms known from previous studies of the real zeros (a multidimensional Gaussian integral with the integrand multiplied by the absolute values of the variables) and the complex zeros (a Hafnian).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Derivation of an eigenvalue probability density function relating to the Poincare disk

A result of Zyczkowski and Sommers [J.Phys.A, 33, 2045--2057 (2000)] gives the eigenvalue probability density function for the top N x N sub-block of a Haar distributed matrix from U(N+n). In the case n \ge N, we rederive this result, starting from knowledge of the distribution of the sub-blocks, introducing the Schur decomposition, and integrating over all variables except the eigenvalues. The integration is done by identifying a recursive structure which reduces the dimension. This approach is inspired by an analogous approach which has been recently applied to determine the eigenvalue probability density function for random matrices A^{-1} B, where A and B are random matrices with entries standard complex normals. We relate the eigenvalue distribution of the sub-blocks to a many body quantum state, and to the one-component plasma, on the pseudosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Matrix averages relating to the Ginibre ensembles

The theory of zonal polynomials is used to compute the average of a Schur polynomial of argument $AX$, where $A$ is a fixed matrix and $X$ is from the real Ginibre ensemble. This generalizes a recent result of Sommers and Khorozhenko [J. Phys. A {\bf 42} (2009), 222002], and furthermore allows analogous results to be obtained for the complex and real quaternion Ginibre ensembles. As applications, the positive integer moments of the general variance Ginibre ensembles are computed in terms of generalized hypergeometric functions, these are written in terms of averages over matrices of the same size as the moment to give duality formulas, and the averages of the power sums of the eigenvalues are expressed as finite sums of zonal polynomials.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Tridiagonal realization of the anti-symmetric Gaussian $β$-ensemble

The Householder reduction of a member of the anti-symmetric Gaussian unitary ensemble gives an anti-symmetric tridiagonal matrix with all independent elements. The random variables permit the introduction of a positive parameter $β$, and the eigenvalue probability density function of the corresponding random matrices can be computed explicitly, as can the distribution of $\{q_i\}$, the first components of the eigenvectors. Three proofs are given. One involves an inductive construction based on bordering of a family of random matrices which are shown to have the same distributions as the anti-symmetric tridiagonal matrices. This proof uses the Dixon-Anderson integral from Selberg integral theory. A second proof involves the explicit computation of the Jacobian for the change of variables between real anti-symmetric tridiagonal matrices, its eigenvalues and $\{q_i\}$. The third proof maps matrices from the anti-symmetric Gaussian $β$-ensemble to those realizing particular examples of the Laguerre $β$-ensemble. In addition to these proofs, we note some simple properties of the shooting eigenvector and associated Prüfer phases of the random matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Determinantal Correlations for Classical Projection Processes

Recent applications in queuing theory and statistical mechanics have isolated the process formed by the eigenvalues of successive minors of the GUE. Analogous eigenvalue processes, formed in general from the eigenvalues of nested sequences of matrices resulting from random corank 1 projections of classical random matrix ensembles, are identified for the LUE and JUE. The correlations for all these processes can be computed in a unified way. The resulting expressions can then be analyzed in various scaling limits. At the soft edge, with the rank of the minors differing by an amount proportional to $N^{2/3}$, the scaled correlations coincide with those known from the soft edge scaling of the Dyson Brownian motion model.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Skew orthogonal polynomials and the partly symmetric real Ginibre ensemble

The partly symmetric real Ginibre ensemble consists of matrices formed as linear combinations of real symmetric and real anti-symmetric Gaussian random matrices. Such matrices typically have both real and complex eigenvalues. For a fixed number of real eigenvalues, an earlier work has given the explicit form of the joint eigenvalue probability density function. We use this to derive a Pfaffian formula for the corresponding summed up generalized partition function. This Pfaffian formula allows the probability that there are exactly $k$ eigenvalues to be written as a determinant with explicit entries. It can be used too to give the explicit form of the correlation functions, provided certain skew orthogonal polynomials are computed. This task is accomplished in terms of Hermite polynomials, and allows us to proceed to analyze various scaling limits of the correlations, including that in which the matrices are only weakly non-symmetric.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Anti-Symmetric GUE Minor Process

Our study is initiated by a multi-component particle system underlying the tiling of a half hexagon by three species of rhombi. In this particle system species $j$ consists of $\lfloor j/2 \rfloor$ particles which are interlaced with neigbouring species. The joint probability density function (PDF) for this particle system is obtained, and is shown in a suitable scaling limit to coincide with the joint eigenvalue PDF for the process formed by the successive minors of anti-symmetric GUE matrices, which in turn we compute from first principles. The correlations for this process are determinantal and we give an explicit formula for the corresponding correlation kernel in terms of Hermite polynomials. Scaling limits of the latter are computed, giving rise to the Airy kernel, extended Airy kernel and bead kernel at the soft edge and in the bulk, as well as a new kernel at the hard edge.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A random matrix decimation procedure relating $β= 2/(r+1)$ to $β= 2(r+1)$

Classical random matrix ensembles with orthogonal symmetry have the property that the joint distribution of every second eigenvalue is equal to that of a classical random matrix ensemble with symplectic symmetry. These results are shown to be the case $r=1$ of a family of inter-relations between eigenvalue probability density functions for generalizations of the classical random matrix ensembles referred to as $β$-ensembles. The inter-relations give that the joint distribution of every $(r+1)$-st eigenvalue in certain $β$-ensembles with $β= 2/(r+1)$ is equal to that of another $β$-ensemble with $β= 2(r+1)$. The proof requires generalizing a conditional probability density function due to Dixon and Anderson.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Eigenvalue statistics of the real Ginibre ensemble

The real Ginibre ensemble consists of random $N \times N$ matrices formed from i.i.d. standard Gaussian entries. By using the method of skew orthogonal polynomials, the general $n$-point correlations for the real eigenvalues, and for the complex eigenvalues, are given as $n \times n$ Pfaffians with explicit entries. A computationally tractable formula for the cumulative probability density of the largest real eigenvalue is presented. This is relevant to May's stability analysis of biological webs.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Symmetrized models of last passage percolation and non-intersecting lattice paths

It has been shown that the last passage time in certain symmetrized models of directed percolation can be written in terms of averages over random matrices from the classical groups $U(l)$, $Sp(2l)$ and $O(l)$. We present a theory of such results based on non-intersecting lattice paths, and integration techniques familiar from the theory of random matrices. Detailed derivations of probabilities relating to two further symmetrizations are also given.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A note on biorthogonal ensembles

We consider ensembles of random matrices, known as biorthogonal ensembles, whose eigenvalue probability density function can be written as a product of two determinants. These systems are closely related to multiple orthogonal functions. It is known that the eigenvalue correlation functions of such ensembles can be written as a determinant of a kernel function. We show that the kernel is itself an average of a single ratio of characteristic polynomials. In the same vein, we prove that the type I multiple polynomials can be expressed as an average of the inverse of a characteristic polynomial. We finally introduce a new biorthogonal matrix ensemble, namely the chiral unitary perturbed by a source term.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Hermite and Laguerre $β$-ensembles: asymptotic corrections to the eigenvalue density

We consider Hermite and Laguerre $β$-ensembles of large $N\times N$ random matrices. For all $β$ even, corrections to the limiting global density are obtained, and the limiting density at the soft edge is evaluated. We use the saddle point method on multidimensional integral representations of the density which are based on special realizations of the generalized (multivariate) classical orthogonal polynomials. The corrections to the bulk density are oscillatory terms that depends on $β$. At the edges, the density can be expressed as a multiple integral of the Konstevich type which constitutes a $β$-deformation of the Airy function. This allows us to obtain the main contribution to the soft edge density when the spectral parameter tends to $\pm\infty$.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Relationships between $τ$-function and Fredholm determinant expressions for gap probabilities in random matrix theory

The gap probabilities at the hard and soft edges of scaled random matrix ensembles with orthogonal symmetry are known in terms of $τ$-functions. Extending recent work relating to the soft edge, it is shown that these $τ$-functions, and their generalizations to contain a generating function parameter, can be expressed as Fredholm determinants. These same Fredholm determinants also occur in exact expressions for gap probabilities in scaled random matrix ensembles with unitary and symplectic symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Dynamical Correlations for Circular Ensembles of Random Matrices

Circular Brownian motion models of random matrices were introduced by Dyson and describe the parametric eigenparameter correlations of unitary random matrices. For symmetric unitary, self-dual quaternion unitary and an analogue of antisymmetric hermitian matrix initial conditions, Brownian dynamics toward the unitary symmetry is analyzed. The dynamical correlation functions of arbitrary number of Brownian particles at arbitrary number of times are shown to be written in the forms of quaternion determinants, similarly as in the case of hermitian random matrix models.

0 citations0 reviewsNo code linkedOpen access