Researcher profile

A. B. J. Kuijlaars

A. B. J. Kuijlaars contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
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

22 published item(s)

preprint2018arXiv

The local universality of Muttalib-Borodin biorthogonal ensembles with parameter $θ= \frac{1}{2}$

The Muttalib-Borodin biorthogonal ensemble is a probability density function for $n$ particles on the positive real line that depends on a parameter $θ$ and an external field $V$. For $θ=\frac{1}{2}$ we find the large $n$ behavior of the associated correlation kernel with only few restrictions on $V$. The idea is to relate the ensemble to a type II multiple orthogonal polynomial ensemble that can in turn be related to a $3\times 3$ Riemann-Hilbert problem which we then solve with the Deift-Zhou steepest descent method. The main ingredient is the construction of the local parametrix at the origin, with the help of Meijer G-functions, and its matching condition with a global parametrix. We will present a new iterative technique to obtain the matching condition, which we expect to be applicable in more general situations as well.

preprint2012arXiv

Critical behavior in Angelesco ensembles

We consider Angelesco ensembles with respect to two modified Jacobi weights on touching intervals [a,0] and [0,1], for a < 0. As a \to -1 the particles around 0 experience a phase transition. This transition is studied in a double scaling limit, where we let the number of particles of the ensemble tend to infinity while the parameter a tends to -1 at a rate of order n^{-1/2}. The correlation kernel converges, in this regime, to a new kind of universal kernel, the Angelesco kernel K^{Ang}. The result follows from the Deift/Zhou steepest descent analysis, applied to the Riemann-Hilbert problem for multiple orthogonal polynomials.

preprint2011arXiv

Universality

Universality of eigenvalue spacings is one of the basic characteristics of random matrices. We give the precise meaning of universality and discuss the standard universality classes (sine, Airy, Bessel) and their appearance in unitary, orthogonal, and symplectic ensembles. The Riemann-Hilbert problem for orthogonal polynomials is one possible tool to derive universality in unitary random matrix ensembles. An overview is presented of the Deift/Zhou steepest descent analysis of the Riemann-Hilbert problem in the one-cut regular case. Non-standard universality classes that arise at singular points in the spectrum are discussed at the end.

preprint2010arXiv

Asymptotic zero distribution of complex orthogonal polynomials associated with Gaussian quadrature

In this paper we study the asymptotic behavior of a family of polynomials which are orthogonal with respect to an exponential weight on certain contours of the complex plane. The zeros of these polynomials are the nodes for complex Gaussian quadrature of an oscillatory integral on the real axis with a high order stationary point, and their limit distribution is also analyzed. We show that the zeros accumulate along a contour in the complex plane that has the S-property in an external field. In addition, the strong asymptotics of the orthogonal polynomials is obtained by applying the nonlinear Deift--Zhou steepest descent method to the corresponding Riemann--Hilbert problem.

preprint2010arXiv

Non-intersecting squared Bessel paths: critical time and double scaling limit

We consider the double scaling limit for a model of $n$ non-intersecting squared Bessel processes in the confluent case: all paths start at time $t=0$ at the same positive value $x=a$, remain positive, and are conditioned to end at time $t=1$ at $x=0$. After appropriate rescaling, the paths fill a region in the $tx$--plane as $n\to \infty$ that intersects the hard edge at $x=0$ at a critical time $t=t^{*}$. In a previous paper (arXiv:0712.1333), the scaling limits for the positions of the paths at time $t\neq t^{*}$ were shown to be the usual scaling limits from random matrix theory. Here, we describe the limit as $n\to \infty$ of the correlation kernel at critical time $t^{*}$ and in the double scaling regime. We derive an integral representation for the limit kernel which bears some connections with the Pearcey kernel. The analysis is based on the study of a $3\times 3$ matrix valued Riemann-Hilbert problem by the Deift-Zhou steepest descent method. The main ingredient is the construction of a local parametrix at the origin, out of the solutions of a particular third-order linear differential equation, and its matching with a global parametrix.

preprint2008arXiv

Asymptotics for a special solution of the thirty fourth Painleve equation

In a previous paper we studied the double scaling limit of unitary random matrix ensembles of the form Z_{n,N}^{-1} |\det M|^{2α} e^{-N \Tr V(M)} dM with α> -1/2. The factor |\det M|^{2α} induces critical eigenvalue behavior near the origin. Under the assumption that the limiting mean eigenvalue density associated with V is regular, and that the origin is a right endpoint of its support, we computed the limiting eigenvalue correlation kernel in the double scaling limit as n, N \to \infty such that n^{2/3}(n/N-1) = O(1) by using the Deift-Zhou steepest descent method for the Riemann-Hilbert problem for polynomials on the line orthogonal with respect to the weight |x|^{2α} e^{-NV(x)}. Our main attention was on the construction of a local parametrix near the origin by means of the ψ-functions associated with a distinguished solution u_α of the Painleve XXXIV equation. This solution is related to a particular solution of the Painleve II equation, which however is different from the usual Hastings-McLeod solution. In this paper we compute the asymptotic behavior of u_α(s) as s \to \pm \infty. We conjecture that this asymptotics characterizes u_α and we present supporting arguments based on the asymptotic analysis of a one-parameter family of solutions of the Painleve XXXIV equation which includes u_α. We identify this family as the family of tronquee solutions of the thirty fourth Painleve equation.

preprint2007arXiv

Asymptotics for Hermite-Pade rational approximants for two analytic functions with separated pairs of branch points (case of genus 0)

We investigate the asymptotic behavior for type II Hermite-Pade approximation to two functions, where each function has two branch points and the pairs of branch points are separated. We give a classification of the cases such that the limiting counting measures for the poles of the Hermite-Pade approximants are described by an algebraic function of order 3 and genus 0. This situation gives rise to a vector-potential equilibrium problem for three measures and the poles of the common denominator are asymptotically distributed like one of these measures. We also work out the strong asymptotics for the corresponding Hermite-Pade approximants by using a 3x3 Riemann-Hilbert problem that characterizes this Hermite-Pade approximation problem.

preprint2007arXiv

Critical edge behavior in unitary random matrix ensembles and the thirty fourth Painleve transcendent

We describe a new universality class for unitary invariant random matrix ensembles. It arises in the double scaling limit of ensembles of random $n \times n$ Hermitian matrices $Z_{n,N}^{-1} |\det M|^{2α} e^{-N \Tr V(M)} dM$ with $α> -1/2$, where the factor $|\det M|^{2α}$ induces critical eigenvalue behavior near the origin. Under the assumption that the limiting mean eigenvalue density associated with $V$ is regular, and that the origin is a right endpoint of its support, we compute the limiting eigenvalue correlation kernel in the double scaling limit as $n, N \to \infty$ such that $n^{2/3}(n/N-1) = O(1)$. We use the Deift-Zhou steepest descent method for the Riemann-Hilbert problem for polynomials on the line orthogonal with respect to the weight $|x|^{2α} e^{-NV(x)}$. Our main attention is on the construction of a local parametrix near the origin by means of the $ψ$-functions associated with a distinguished solution of the Painleve XXXIV equation. This solution is related to a particular solution of the Painleve II equation, which however is different from the usual Hastings-McLeod solution.

preprint2007arXiv

Locating the zeros of partial sums of exp(z) with Riemann-Hilbert methods

In this paper we derive uniform asymptotic expansions for the partial sums of the exponential series. We indicate how this information will be used in a later publication to obtain full and explicitly computable asymptotic expansions with error bounds for all zeros of the Taylor polynomials $p_{n-1}(z) = \sum_{k=0}^{n-1} z^k/ k!$. Our proof is based on a representation of $p_{n-1}(nz)$ in terms of an integral of the form $\int_γ \frac{e^{nϕ(s)}}{s-z}ds$. We demonstrate how to derive uniform expansions for such integrals using a Riemann-Hilbert approach. A comparison with classical steepest descent analysis shows the advantages of the Riemann-Hilbert analysis in particular for points $z$ that are close to the critical points of $ϕ$.

preprint2007arXiv

The asymptotic behaviour of recurrence coefficients for orthogonal polynomials with varying exponential weights

We consider orthogonal polynomials $\{p_{n,N}(x)\}_{n=0}^{\infty}$ on the real line with respect to a weight $w(x)=e^{-NV(x)}$ and in particular the asymptotic behaviour of the coefficients $a_{n,N}$ and $b_{n,N}$ in the three term recurrence $x π_{n,N}(x) = π_{n+1,N}(x) + b_{n,N} π_{n,N}(x) + a_{n,N} π_{n-1,N}(x)$. For one-cut regular $V$ we show, using the Deift-Zhou method of steepest descent for Riemann-Hilbert problems, that the diagonal recurrence coefficients $a_{n,n}$ and $b_{n,n}$ have asymptotic expansions as $n \to \infty$ in powers of $1/n^2$ and powers of $1/n$, respectively.

preprint2005arXiv

Multi-critical unitary random matrix ensembles and the general Painleve II equation

We study unitary random matrix ensembles of the form $Z_{n,N}^{-1} |\det M|^{2α} e^{-N \Tr V(M)}dM$, where $α>-1/2$ and $V$ is such that the limiting mean eigenvalue density for $n,N\to\infty$ and $n/N\to 1$ vanishes quadratically at the origin. In order to compute the double scaling limits of the eigenvalue correlation kernel near the origin, we use the Deift/Zhou steepest descent method applied to the Riemann-Hilbert problem for orthogonal polynomials on the real line with respect to the weight $|x|^{2α}e^{-NV(x)}$. Here the main focus is on the construction of a local parametrix near the origin with $ψ$-functions associated with a special solution $q_α$ of the Painlevé II equation $q&#39;&#39;=sq+2q^3-α$. We show that $q_α$ has no real poles for $α> -1/2$, by proving the solvability of the corresponding Riemann-Hilbert problem. We also show that the asymptotics of the recurrence coefficients of the orthogonal polynomials can be expressed in terms of $q_α$ in the double scaling limit.

preprint2005arXiv

Multiple orthogonal polynomials of mixed type and non-intersecting Brownian motions

We present a generalization of multiple orthogonal polynomials of type I and type II, which we call multiple orthogonal polynomials of mixed type. Some basic properties are formulated, and a Riemann-Hilbert problem for the multiple orthogonal polynomials of mixed type is given. We derive a Christoffel-Darboux formula for these polynomials using the solution of the Riemann-Hilbert problem. The main motivation for studying these polynomials comes from a model of non-intersecting one-dimensional Brownian motions with a given number of starting points and endpoints. The correlation kernel for the positions of the Brownian paths at any intermediate time coincides with the Christoffel-Darboux kernel for the multiple orthogonal polynomials of mixed type with respect to Gaussian weights.

preprint2005arXiv

Type II Hermite-Padé approximation to the exponential function

We obtain strong and uniform asymptotics in every domain of the complex plane for the scaled polynomials $a (3nz)$, $b (3nz)$, and $c (3nz)$ where $a$, $b$, and $c$ are the type II Hermite-Padé approximants to the exponential function of respective degrees $2n+2$, $2n$ and $2n$, defined by $a (z)e^{-z}-b (z)=Ø(z^{3n+2})$ and $a (z)e^{z}-c (z)=Ø(z^{3n+2})$ as $z\to 0$. Our analysis relies on a characterization of these polynomials in terms of a $3\times 3$ matrix Riemann-Hilbert problem which, as a consequence of the famous Mahler relations, corresponds by a simple transformation to a similar Riemann-Hilbert problem for type I Hermite-Padé approximants. Due to this relation, the study that was performed in previous work, based on the Deift-Zhou steepest descent method for Riemann-Hilbert problems, can be reused to establish our present results.

preprint2004arXiv

A Christoffel-Darboux formula for multiple orthogonal polynomials

Bleher and Kuijlaars recently showed that the eigenvalue correlations from matrix ensembles with external source can be expressed by means of a kernel built out of special multiple orthogonal polynomials. We derive a Christoffel-Darboux formula for this kernel for general multiple orthogonal polynomials. In addition, we show that the formula can be written in terms of the solution of the Riemann-Hilbert problem for multiple orthogonal polynomials, which will be useful for asymptotic analysis.

preprint2003arXiv

Orthogonality of Jacobi polynomials with general parameters

In this paper we study the orthogonality conditions satisfied by Jacobi polynomials $P_n^{(α,β)}$ when the parameters $α$ and $β$ are not necessarily $>-1$. We establish orthogonality on a generic closed contour on a Riemann surface. Depending on the parameters, this leads to either full orthogonality conditions on a single contour in the plane, or to multiple orthogonality conditions on a number of contours in the plane. In all cases we show that the orthogonality conditions characterize the Jacobi polynomial $P_n^{(α, β)}$ of degree $n$ up to a constant factor.

preprint2003arXiv

Quadratic Hermite-Pade approximation to the exponential function: a Riemann-Hilbert approach

We investigate the asymptotic behavior of the polynomials p, q, r of degrees n in type I Hermite-Pade approximation to the exponential function, defined by p(z)e^{-z}+q(z)+r(z)e^{z} = O(z^{3n+2}) as z -> 0. These polynomials are characterized by a Riemann-Hilbert problem for a 3x3 matrix valued function. We use the Deift-Zhou steepest descent method for Riemann-Hilbert problems to obtain strong uniform asymptotics for the scaled polynomials p(3nz), q(3nz), and r(3nz) in every domain in the complex plane. An important role is played by a three-sheeted Riemann surface and certain measures and functions derived from it. Our work complements recent results of Herbert Stahl.

preprint2003arXiv

Random matrices with external source and multiple orthogonal polynomials

We show that the average characteristic polynomial P_n(z) = E [\det(zI-M)] of the random Hermitian matrix ensemble Z_n^{-1} \exp(-Tr(V(M)-AM))dM is characterized by multiple orthogonality conditions that depend on the eigenvalues of the external source A. For each eigenvalue a_j of A, there is a weight and P_n has n_j orthogonality conditions with respect to this weight, if n_j is the multiplicity of a_j. The eigenvalue correlation functions have determinantal form, as shown by Zinn-Justin. Here we give a different expression for the kernel. We derive a Christoffel-Darboux formula in case A has two distinct eigenvalues, which leads to a compact formula in terms of a Riemann-Hilbert problem that is satisfied by multiple orthogonal polynomials.

preprint2003arXiv

The Riemann-Hilbert approach to strong asymptotics for orthogonal polynomials on [-1,1]

We consider polynomials that are orthogonal on $[-1,1]$ with respect to a modified Jacobi weight $(1-x)^α(1+x)^βh(x)$, with $α,β>-1$ and $h$ real analytic and stricly positive on $[-1,1]$. We obtain full asymptotic expansions for the monic and orthonormal polynomials outside the interval $[-1,1]$, for the recurrence coefficients and for the leading coefficients of the orthonormal polynomials. We also deduce asymptotic behavior for the Hankel determinants. For the asymptotic analysis we use the steepest descent technique for Riemann--Hilbert problems developed by Deift and Zhou, and applied to orthogonal polynomials on the real line by Deift, Kriecherbauer, McLaughlin, Venakides, and Zhou. In the steepest descent method we will use the Szegő function associated with the weight and for the local analysis around the endpoints $\pm 1$ we use Bessel functions of appropriate order, whereas Deift et al. use Airy functions.

preprint2002arXiv

Asymptotic zero behavior of Laguerre polynomials with negative parameter

We consider Laguerre polynomials $L_n^{(α_n)}(nz)$ with varying negative parameters $α_n$, such that the limit $A = -\lim_n α_n/n$ exists and belongs to $(0,1)$. For $A > 1$, it is known that the zeros accumulate along an open contour in the complex plane. For every $A \in (0,1)$, we describe a one-parameter family of possible limit sets of the zeros. Under the condition that the limit $r= - \lim_n \frac{1}{n} \log \dist(α_n, \mathbb Z)$ exists, we show that the zeros accumulate on $Γ_r \cup [β_1,β_2]$ with $β_1$ and $β_2$ only depending on $A$. For $r \in [0,\infty)$, $Γ_r$ is a closed loop encircling the origin, which for $r = +\infty$, reduces to the origin. This shows a great sensitivity of the zeros to $α_n$&#39;s proximity to the integers. We use a Riemann-Hilbert formulation for the Laguerre polynomials, together with the steepest descent method of Deift and Zhou to obtain asymptotics for the polynomials, from which the zero behavior follows.

preprint2002arXiv

Direct and inverse spectral transform for the relativistic Toda lattice and the connection with Laurent orthogonal polynomials

We introduce a spectral transform for the finite relativistic Toda lattice (RTL) in generalized form. In the nonrelativistic case, Moser constructed a spectral transform from the spectral theory of symmetric Jacobi matrices. Here we use a non-symmetric generalized eigenvalue problem for a pair of bidiagonal matrices (L,M) to define the spectral transform for the RTL. The inverse spectral transform is described in terms of a terminating T-fraction. The generalized eigenvalues are constants of motion and the auxiliary spectral data have explicit time evolution. Using the connection with the theory of Laurent orthogonal polynomials, we study the long-time behaviour of the RTL. As in the case of the Toda lattice the matrix entries have asymptotic limits. We show that L tends to an upper Hessenberg matrix with the generalized eigenvalues sorted on the diagonal, while M tends to the identity matrix.

preprint2002arXiv

Riemann--Hilbert analysis for Laguerre polynomials with large negative parameter

We study the asymptotic behavior of Laguerre polynomials $L_n^{(α_n)}(nz)$ as $n \to \infty$, where $α_n$ is a sequence of negative parameters such that $-α_n/n$ tends to a limit $A > 1$ as $n \to \infty$. These polynomials satisfy a non-hermitian orthogonality on certain contours in the complex plane. This fact allows the formulation of a Riemann--Hilbert problem whose solution is given in terms of these Laguerre polynomials. The asymptotic analysis of the Riemann--Hilbert problem is carried out by the steepest descent method of Deift and Zhou, in the same spirit as done by Deift et al. for the case of orthogonal polynomials on the real line. A main feature of the present paper is the choice of the correct contour.