
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.

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.

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 HarishChandra/ItzyksonZuber
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.

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
multitime generalizations of the microscopic results.

The MuttalibBorodin 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.

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_kx_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 rankone modification of
a Hermitian matrix, and the extension of a Hermitian matrix by adding an extra
row and column with complex Gaussians.

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 Gkernels found for products of complex
Ginibre random matrices.

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 starlike
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 RiemannHilbert
problem of size $(d+1)\times (d+1)$. The vector equilibrium problem and the
RiemannHilbert problem that we investigate are generalizations of those
studied recently by BleherKuijlaars in the case $d=2$.

The normal matrix model with a cubic potential is illdefined 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 RiemannHilbert problem. It follows that the zeros of
the orthogonal polynomials tend to $\Sigma_1$, with the exception of at most
three spurious zeros.

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 RiemannHilbert problem for orthogonal
polynomials, suitably modified to cover the present situation, and the
DeiftZhou 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$.

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 M1
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.

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 Gfunctions. 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 CauchyLaguerre twomatrix model, which indicates that these kernels
represent a new universality class in random matrix theory.

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.

We consider the normal matrix model with a cubic potential. The model is
illdefined, and in order to reguralize it, Elbau and Felder introduced a model
with a cutoff and corresponding system of orthogonal polynomials with respect
to a varying exponential weight on the cutoff 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 cutoff,
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 RiemannHilbert (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 DeiftZhou 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.

We study multiple orthogonal polynomials of MeixnerPollaczek 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 socalled sixvertex model with domain wall boundary conditions,
in a particular regime known as the free fermion line. We show how the multiple
MeixnerPollaczek polynomials arise in an inhomogeneous version of this model.

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 MehlerHeine type
asymptotic formula for the polynomials in this regime. The method used to
analyze the problem is the steepest descent technique for RiemannHilbert
problems. A key point in the analysis is the construction of a new local
parametrix.

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
RiemannHilbert 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.

We study a model of $n$ onedimensional nonintersecting 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 timespace 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 RiemannHilbert problem of size $4\times 4$. We
prove solvability of the RiemannHilbert problem and establish a remarkable
connection with the HastingsMcLeod 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 nonintersecting 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.

We consider the two sequences of biorthogonal polynomials (p_{k,n})_k and
(q_{k,n})_k related to the Hermitian twomatrix 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 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 RiemannHilbert 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 nonintersecting
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.

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 RiemannHilbert 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 taplane, where either the constraint becomes
active, or the two intervals in the support come together (or both).

We consider n nonintersecting 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 RiemannHilbert 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 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.

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$ matrixvalued RiemannHilbert 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 onecut 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.

We study polynomials that are orthogonal with respect to the modified
Laguerre weight $z^{n + \nu} e^{Nz} (z1)^{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 $(z1)^{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
RiemannHilbert problem associated with orthogonal polynomials.