• ### Universality for conditional measures of the sine point process(1703.02349)

March 21, 2019 math.CA, math.PR
The sine process is a rigid point process on the real line, which means that for almost all configurations $X$, the number of points in an interval $I = [-R,R]$ is determined by the points of $X$ outside of $I$. In addition, the points in $I$ are an orthogonal polynomial ensemble on $I$ with a weight function that is determined by the points in $X \setminus I$. We prove a universality result that in particular implies that the correlation kernel of the orthogonal polynomial ensemble tends to the sine kernel as the length $|I|=2R$ tends to infinity, thereby answering a question posed by A.I. Bufetov.
• ### Exceptional Laguerre polynomials(1708.03106)

Jan. 23, 2018 math.CA
The aim of this paper is to present the construction of exceptional Laguerre polynomials in a systematic way, and to provide new asymptotic results on the location of the zeros. To describe the exceptional Laguerre polynomials we associate them with two partitions. We find that the use of partitions is an elegant way to express these polynomials and we restate some of their known properties in terms of partitions. We discuss the asymptotic behavior of the regular zeros and the exceptional zeros of exceptional Laguerre polynomials as the degree tends to infinity.
• ### Spherical functions approach to sums of random Hermitian matrices(1611.08932)

Nov. 27, 2016 math-ph, math.MP, math.CA, math.PR, math.RT
We present an approach to sums of random Hermitian matrices via the theory of spherical functions for the Gelfand pair $(\mathrm{U}(n) \ltimes \mathrm{Herm}(n), \mathrm{U}(n))$. It is inspired by a similar approach of Kieburg and K\"osters for products of random matrices. The spherical functions have determinantal expressions because of the Harish-Chandra/Itzykson-Zuber integral formula. It leads to remarkably simple expressions for the spherical transform and its inverse. The spherical transform is applied to sums of unitarily invariant random matrices from polynomial ensembles and the subclass of polynomial ensembles of derivative type (in the additive sense), which turns out to be closed under addition. We finally present additional detailed calculations for the sum with a random matrix from a Laguerre Unitary Ensemble.
• ### Propagation of singular behavior for Gaussian perturbations of random matrices(1608.05870)

Aug. 20, 2016 math-ph, math.MP, math.CA, math.PR
We study the asymptotic behavior of the eigenvalues of Gaussian perturbations of large Hermitian random matrices for which the limiting eigenvalue density vanishes at a singular interior point or vanishes faster than a square root at a singular edge point. First, we show that the singular behavior propagates macroscopically for sufficiently small Gaussian perturbations, and we describe the macroscopic eigenvalue behavior for Gaussian perturbations of critical size. Secondly, for sufficiently small Gaussian perturbations of unitary invariant random matrices, we prove that the microscopic eigenvalue correlations near the singular point are described by the same limiting kernel as in the unperturbed case. We also interpret our results in terms of nonintersecting Brownian paths with random starting positions, and we establish multi-time generalizations of the microscopic results.
• ### A Vector Equilibrium Problem for Muttalib-Borodin Biorthogonal Ensembles(1604.03070)

July 5, 2016 math-ph, math.MP, math.CA, math.PR
The Muttalib-Borodin biorthogonal ensemble is a joint density function for $n$ particles on the positive real line that depends on a parameter $\theta$. There is an equilibrium problem that describes the large $n$ behavior. We show that for rational values of $\theta$ there is an equivalent vector equilibrium problem.
• ### Transformations of polynomial ensembles(1501.05506)

Jan. 30, 2015 math-ph, math.MP, math.PR
A polynomial ensemble is a probability density function for the position of $n$ real particles of the form $\frac{1}{Z_n} \, \prod_{j<k} (x_k-x_j) \, \det \left[ f_k (x_j) \right]_{j,k=1}^n$, for certain functions $f_1, \ldots, f_n$. Such ensembles appear frequently as the joint eigenvalue density of random matrices. We present a number of transformations that preserve the structure of a polynomial ensemble. These transformations include the restriction of a Hermitian matrix by removing one row and one column, a rank-one modification of a Hermitian matrix, and the extension of a Hermitian matrix by adding an extra row and column with complex Gaussians.
• ### Singular value statistics of matrix products with truncated unitary matrices(1501.03910)

Jan. 16, 2015 math-ph, math.MP, math.PR
We prove that the squared singular values of a fixed matrix multiplied with a truncation of a Haar distributed unitary matrix are distributed by a polynomial ensemble. This result is applied to a multiplication of a truncated unitary matrix with a random matrix. We show that the structure of polynomial ensembles and of certain Pfaffian ensembles is preserved. Furthermore we derive the joint singular value density of a product of truncated unitary matrices and its corresponding correlation kernel which can be written as a double contour integral. This leads to hard edge scaling limits that also include new finite rank perturbations of the Meijer G-kernels found for products of complex Ginibre random matrices.
• ### The normal matrix model with a monomial potential, a vector equilibrium problem, and multiple orthogonal polynomials on a star(1401.2419)

Jan. 16, 2015 math-ph, math.MP, math.CA
We investigate the asymptotic behavior of a family of multiple orthogonal polynomials that is naturally linked with the normal matrix model with a monomial potential of arbitrary degree $d+1$. The polynomials that we investigate are multiple orthogonal with respect to a system of $d$ analytic weights defined on a symmetric $(d+1)$-star centered at the origin. In the first part we analyze in detail a vector equilibrium problem involving a system of $d$ interacting measures $(\mu_{1},\ldots,\mu_{d})$ supported on star-like sets in the plane. We show that in the subcritical regime, the first component $\mu_{1}^{*}$ of the solution to this problem is the asymptotic zero distribution of the multiple orthogonal polynomials. It also characterizes the domain where the eigenvalues in the normal matrix model accumulate, in the sense that the Schwarz function associated with the boundary of this domain can be expressed explicitly in terms of $\mu_{1}^{*}$. The second part of the paper is devoted to the asymptotic analysis of the multiple orthogonal polynomials. The asymptotic results are obtained again in the subcritical regime, and they follow from the Deift/Zhou steepest descent analysis of a Riemann-Hilbert problem of size $(d+1)\times (d+1)$. The vector equilibrium problem and the Riemann-Hilbert problem that we investigate are generalizations of those studied recently by Bleher-Kuijlaars in the case $d=2$.
• ### The supercritical regime in the normal matrix model with cubic potential(1412.7597)

Dec. 27, 2014 math-ph, math.MP
The normal matrix model with a cubic potential is ill-defined and it develops a critical behavior in finite time. We follow the approach of Bleher and Kuijlaars to reformulate the model in terms of orthogonal polynomials with respect to a Hermitian form. This reformulation was shown to capture the essential features of the normal matrix model in the subcritical regime, namely that the zeros of the polynomials tend to a number of segments (the motherbody) inside a domain (the droplet) that attracts the eigenvalues in the normal matrix model. In the present paper we analyze the supercritical regime and we find that the large $n$ behavior is described by the evolution of a spectral curve satisfying the Boutroux condition. The Boutroux condition determines a system of contours $\Sigma_1$, consisting of the motherbody and whiskers sticking out of the domain. We find a second critical behavior at which the original motherbody shrinks to a point at the origin and only the whiskers remain. In the regime before the second criticality we also give strong asymptotics of the orthogonal polynomials by means of a steepest descent analysis of a $3 \times 3$ matrix valued Riemann-Hilbert problem. It follows that the zeros of the orthogonal polynomials tend to $\Sigma_1$, with the exception of at most three spurious zeros.
• ### Asymptotic behavior and zero distribution of polynomials orthogonal with respect to Bessel functions(1406.0969)

June 4, 2014 math.CA, math.CV
We consider polynomials $P_n$ orthogonal with respect to the weight $J_{\nu}$ on $[0,\infty)$, where $J_{\nu}$ is the Bessel function of order $\nu$. Asheim and Huybrechs considered these polynomials in connection with complex Gaussian quadrature for oscillatory integrals. They observed that the zeros are complex and accumulate as $n \to \infty$ near the vertical line $\textrm{Re}\, z = \frac{\nu \pi}{2}$. We prove this fact for the case $0 \leq \nu \leq 1/2$ from strong asymptotic formulas that we derive for the polynomials $P_n$ in the complex plane. Our main tool is the Riemann-Hilbert problem for orthogonal polynomials, suitably modified to cover the present situation, and the Deift-Zhou steepest descent method. A major part of the work is devoted to the construction of a local parametrix at the origin, for which we give an existence proof that only works for $\nu \leq 1/2$.
• ### Singular values of products of random matrices and polynomial ensembles(1404.5802)

April 23, 2014 math-ph, math.MP, math.CA, math.PR
Akemann, Ipsen, and Kieburg showed recently that the squared singular values of a product of M complex Ginibre matrices are distributed according to a determinantal point process. We introduce the notion of a polynomial ensemble and show how their result can be interpreted as a transformation of polynomial ensembles. We also show that the squared singular values of the product of M-1 complex Ginibre matrices with one truncated unitary matrix is a polynomial ensemble, and we derive a double integral representation for the correlation kernel associated with this ensemble. We use this to calculate the scaling limit at the hard edge, which turns out to be the same scaling limit as the one found by Kuijlaars and Zhang for the squared singular values of a product of M complex Ginibre matrices. Our final result is that these limiting kernels also appear as scaling limits for the biorthogonal ensembles of Borodin with parameter theta > 0, in case theta or 1/theta is an integer. This further supports the conjecture that these kernels have a universal character.
• ### Singular values of products of Ginibre random matrices, multiple orthogonal polynomials and hard edge scaling limits(1308.1003)

Dec. 23, 2013 math-ph, math.MP, math.CA, math.PR
Akemann, Ipsen and Kieburg recently showed that the squared singular values of products of M rectangular random matrices with independent complex Gaussian entries are distributed according to a determinantal point process with a correlation kernel that can be expressed in terms of Meijer G-functions. We show that this point process can be interpreted as a multiple orthogonal polynomial ensemble. We give integral representations for the relevant multiple orthogonal polynomials and a new double contour integral for the correlation kernel, which allows us to find its scaling limits at the origin (hard edge). The limiting kernels generalize the classical Bessel kernels. For M=2 they coincide with the scaling limits found by Bertola, Gekhtman, and Szmigielski in the Cauchy-Laguerre two-matrix model, which indicates that these kernels represent a new universality class in random matrix theory.
• ### Large Deviations for a Non-Centered Wishart Matrix(1204.6261)

March 12, 2013 math.CA, math.PR
We investigate an additive perturbation of a complex Wishart random matrix and prove that a large deviation principle holds for the spectral measures. The rate function is associated to a vector equilibrium problem coming from logarithmic potential theory, which in our case is a quadratic map involving the logarithmic energies, or Voiculescu's entropies, of two measures in the presence of an external field and an upper constraint. The proof is based on a two type particles Coulomb gas representation for the eigenvalue distribution, which gives a new insight on why such variational problems should describe the limiting spectral distribution. This representation is available because of a Nikishin structure satisfied by the weights of the multiple orthogonal polynomials hidden in the background.
• ### Orthogonal polynomials in the normal matrix model with a cubic potential(1106.6168)

June 30, 2011 math-ph, math.MP, math.CA
We consider the normal matrix model with a cubic potential. The model is ill-defined, and in order to reguralize it, Elbau and Felder introduced a model with a cut-off and corresponding system of orthogonal polynomials with respect to a varying exponential weight on the cut-off region on the complex plane. In the present paper we show how to define orthogonal polynomials on a specially chosen system of infinite contours on the complex plane, without any cut-off, which satisfy the same recurrence algebraic identity that is asymptotically valid for the orthogonal polynomials of Elbau and Felder. The main goal of this paper is to develop the Riemann-Hilbert (RH) approach to the orthogonal polynomials under consideration and to obtain their asymptotic behavior on the complex plane as the degree $n$ of the polynomial goes to infinity. As the first step in the RH approach, we introduce an auxiliary vector equilibrium problem for a pair of measures $(\mu_1,\mu_2)$ on the complex plane. We then formulate a $3\times 3$ matrix valued RH problem for the orthogonal polynomials in hand, and we apply the nonlinear steepest descent method of Deift-Zhou to the asymptotic analysis of the RH problem. The central steps in our study are a sequence of transformations of the RH problem, based on the equilibrium vector measure $(\mu_1,\mu_2)$, and the construction of a global parametrix. The main result of this paper is a derivation of the large $n$ asymptotics of the orthogonal polynomials on the whole complex plane. We prove that the distribution of zeros of the orthogonal polynomials converges to the measure $\mu_1$, the first component of the equilibrium measure. We also obtain analytical results for the measure $\mu_1$ relating it to the distribution of eigenvalues in the normal matrix model which is uniform in a domain bounded by a simple closed curve.
• ### Multiple Meixner-Pollaczek polynomials and the six-vertex model(1101.2982)

We study multiple orthogonal polynomials of Meixner-Pollaczek type with respect to a symmetric system of two orthogonality measures. Our main result is that the limiting distribution of the zeros of these polynomials is one component of the solution to a constrained vector equilibrium problem. We also provide a Rodrigues formula and closed expressions for the recurrence coefficients. The proof of the main result follows from a connection with the eigenvalues of block Toeplitz matrices, for which we provide some general results of independent interest. The motivation for this paper is the study of a model in statistical mechanics, the so-called six-vertex model with domain wall boundary conditions, in a particular regime known as the free fermion line. We show how the multiple Meixner-Pollaczek polynomials arise in an inhomogeneous version of this model.
• ### Double scaling limit for modified Jacobi-Angelesco polynomials(1102.1349)

Feb. 7, 2011 math.CA
We consider multiple orthogonal polynomials with respect to two modified Jacobi weights on touching intervals [a,0] and [0,1], with a < 0, and study a transition that occurs at a = -1. The transition is studied in a double scaling limit, where we let the degree n of the polynomial tend to infinity while the parameter a tends to -1 at a rate of O(n^{-1/2}). We obtain a Mehler-Heine type asymptotic formula for the polynomials in this regime. The method used to analyze the problem is the steepest descent technique for Riemann-Hilbert problems. A key point in the analysis is the construction of a new local parametrix.
• ### The Hermitian two matrix model with an even quartic potential(1010.4282)

Oct. 20, 2010 math-ph, math.MP, math.PR, math.CV
We consider the two matrix model with an even quartic potential W(y)=y^4/4+alpha y^2/2 and an even polynomial potential V(x). The main result of the paper is the formulation of a vector equilibrium problem for the limiting mean density for the eigenvalues of one of the matrices M_1. The vector equilibrium problem is defined for three measures, with external fields on the first and third measures and an upper constraint on the second measure. The proof is based on a steepest descent analysis of a 4 x 4 matrix valued Riemann-Hilbert problem that characterizes the correlation kernel for the eigenvalues of M_1. Our results generalize earlier results for the case alpha=0, where the external field on the third measure was not present.
• ### Critical behavior of non-intersecting Brownian motions at a tacnode(1009.2457)

Sept. 13, 2010 math.CA, math.PR
We study a model of $n$ one-dimensional non-intersecting Brownian motions with two prescribed starting points at time $t=0$ and two prescribed ending points at time $t=1$ in a critical regime where the paths fill two tangent ellipses in the time-space plane as $n \to \infty$. The limiting mean density for the positions of the Brownian paths at the time of tangency consists of two touching semicircles, possibly of different sizes. We show that in an appropriate double scaling limit, there is a new familiy of limiting determinantal point processes with integrable correlation kernels that are expressed in terms of a new Riemann-Hilbert problem of size $4\times 4$. We prove solvability of the Riemann-Hilbert problem and establish a remarkable connection with the Hastings-McLeod solution of the Painlev\'e II equation. We show that this Painlev\'e II transcendent also appears in the critical limits of the recurrence coefficients of the multiple Hermite polynomials that are associated with the non-intersecting Brownian motions. Universality suggests that the new limiting kernels apply to more general situations whenever a limiting mean density vanishes according to two touching square roots, which represents a new universality class.
• ### A vector equilibrium problem for the two-matrix model in the quartic/quadratic case(1007.3137)

July 19, 2010 math-ph, math.MP, math.CA
We consider the two sequences of biorthogonal polynomials (p_{k,n})_k and (q_{k,n})_k related to the Hermitian two-matrix model with potentials V(x) = x^2/2 and W(y) = y^4/4 + ty^2. From an asymptotic analysis of the coefficients in the recurrence relation satisfied by these polynomials, we obtain the limiting distribution of the zeros of the polynomials p_{n,n} as n tends to infinity. The limiting zero distribution is characterized as the first measure of the minimizer in a vector equilibrium problem involving three measures which for the case t=0 reduces to the vector equilibrium problem that was given recently by two of us. A novel feature is that for t < 0 an external field is active on the third measure which introduces a new type of critical behavior for a certain negative value of t. We also prove a general result about the interlacing of zeros of biorthogonal polynomials.
• ### Multiple orthogonal polynomials in random matrix theory(1004.0846)

April 6, 2010 math.CA, math.PR
Multiple orthogonal polynomials are a generalization of orthogonal polynomials in which the orthogonality is distributed among a number of orthogonality weights. They appear in random matrix theory in the form of special determinantal point processes that are called multiple orthogonal polynomial (MOP) ensembles. The correlation kernel in such an ensemble is expressed in terms of the solution of a Riemann-Hilbert problem, that is of size (r+1) x (r+1) in the case of r weights. A number of models give rise to a MOP ensemble, and we discuss recent results on models of non-intersecting Brownian motions, Hermitian random matrices with external source, and the two matrix model. A novel feature in the asymptotic analysis of the latter two models is a vector equilibrium problem for two or more measures, that describes the limiting mean eigenvalue density. The vector equilibrium problems involve both an external field and an upper constraint.
• ### Random matrix model with external source and a constrained vector equilibrium problem(1001.1238)

Jan. 8, 2010 math-ph, math.MP
We consider the random matrix model with external source, in case where the potential V(x) is an even polynomial and the external source has two eigenvalues a, -a of equal multiplicity. We show that the limiting mean eigenvalue distribution of this model can be characterized as the first component of a pair of measures (mu_1,mu_2) that solve a constrained vector equilibrium problem. The proof is based on the steepest descent analysis of the associated Riemann-Hilbert problem for multiple orthogonal polynomials. We illustrate our results in detail for the case of a quartic double well potential V(x) = x^4/4 - tx^2/2. We are able to determine the precise location of the phase transitions in the ta-plane, where either the constraint becomes active, or the two intervals in the support come together (or both).
• ### A graph-based equilibrium problem for the limiting distribution of non-intersecting Brownian motions at low temperature(0907.2310)

July 14, 2009 math.PR, math.CV
We consider n non-intersecting Brownian motion paths with p prescribed starting positions at time t=0 and q prescribed ending positions at time t=1. The positions of the paths at any intermediate time are a determinantal point process, which in the case p=1 is equivalent to the eigenvalue distribution of a random matrix from the Gaussian unitary ensemble with external source. For general p and q, we show that if a temperature parameter is sufficiently small, then the distribution of the Brownian paths is characterized in the large n limit by a vector equilibrium problem with an interaction matrix that is based on a bipartite planar graph. Our proof is based on a steepest descent analysis of an associated (p+q) by (p+q) matrix valued Riemann-Hilbert problem whose solution is built out of multiple orthogonal polynomials. A new feature of the steepest descent analysis is a systematic opening of a large number of global lenses.
• ### Multiple orthogonal polynomial ensembles(0902.1058)

Feb. 6, 2009 math.CA, math.PR
Multiple orthogonal polynomials are traditionally studied because of their connections to number theory and approximation theory. In recent years they were found to be connected to certain models in random matrix theory. In this paper we introduce the notion of a multiple orthogonal polynomial ensemble (MOP ensemble) and derive some of their basic properties. It is shown that Angelesco and Nikishin systems give rise to MOP ensembles and that the equilibrium problems that are associated with these systems have a natural interpretation in the context of MOP ensembles.
• ### Universality in the two matrix model: a Riemann-Hilbert steepest descent analysis(0807.4814)

July 30, 2008 math-ph, math.MP
The eigenvalue statistics of a pair $(M_1,M_2)$ of $n\times n$ Hermitian matrices taken random with respect to the measure $$\frac{1}{Z_n}\exp\big(-n\Tr (V(M_1)+W(M_2)-\tau M_1M_2)\big) {\rm d}M_1 {\rm d} M_2$$ can be described in terms of two families of biorthogonal polynomials. In this paper we give a steepest descent analysis of a $4 \times 4$ matrix-valued Riemann-Hilbert problem characterizing one of the families of biorthogonal polynomials in the special case $W(y)=y^4/4$ and $V$ an even polynomial. As a result we obtain the limiting behavior of the correlation kernel associated to the eigenvalues of $M_1$ (when averaged over $M_2$) in the global and local regime as $n\to \infty$ in the one-cut regular case. A special feature in the analysis is the introduction of a vector equilibrium problem involving both an external field and an upper constraint.
• ### Painleve IV asymptotics for orthogonal polynomials with respect to a modified Laguerre weight(0804.2564)

July 8, 2008 math.CA, math.CV
We study polynomials that are orthogonal with respect to the modified Laguerre weight $z^{-n + \nu} e^{-Nz} (z-1)^{2b}$ in the limit where $n, N \to \infty$ with $N/n \to 1$ and $\nu$ is a fixed number in $\mathbb{R} \setminus \mathbb{N}_0$. With the effect of the factor $(z-1)^{2b}$, the local parametrix near the critical point $z =1$ can be constructed in terms of $\Psi$-functions associated with the Painleve IV equation. We show that the asymptotics of the recurrence coefficients of orthogonal polynomials can be described in terms of specified solution of the Painleve IV equation in the double scaling limit. Our method is based on the Deift/Zhou steepest decent analysis of the Riemann-Hilbert problem associated with orthogonal polynomials.