BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Peter J. Forrester

Peter J. Forrester contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math-phmath.MPmath.PRComputationcond-mat.stat-mechmath.CAmath.COphysics.data-an

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 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