Researcher profile

Arno B. J. Kuijlaars

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

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
7topics
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

15 published item(s)

preprint2022arXiv

Critical measures on higher genus Riemann surfaces

Critical measures in the complex plane are saddle points for the logarithmic energy with external field. Their local and global structure was described by Martinez-Finkelshtein and Rakhmanov. In this paper we start the development of a theory of critical measures on higher genus Riemann surfaces, where the logarithmic energy is replaced by the energy with respect to a bipolar Green's kernel. We study a max-min problem for the bipolar Green's energy with external fields Re V where dV is a meromorphic differential. Under reasonable assumptions the max-min problem has a solution and we show that the corresponding equilibrium measure is a critical measure in the external field. In a special genus one situation we are able to show that the critical measure is supported on maximal trajectories of a meromorphic quadratic differential. We are motivated by applications to random lozenge tilings of a hexagon with periodic weightings. Correlations in these models are expressible in terms of matrix valued orthogonal polynomials. The matrix orthogonality is interpreted as (partial) scalar orthogonality on a Riemann surface. The theory of critical measures will be useful for the asymptotic analysis of a corresponding Riemann-Hilbert problem as we outline in the paper.

preprint2021arXiv

Matrix-valued orthogonal polynomials related to hexagon tilings

In this paper, we study a class of matrix-valued orthogonal polynomials (MVOPs) that are related to 2-periodic lozenge tilings of a hexagon. The general model depends on many parameters. In the cases of constant and $2$-periodic parameter values we show that the MVOP can be expressed in terms of scalar polynomials with non-Hermitian orthogonality on a closed contour in the complex plane. The 2-periodic hexagon tiling model with a constant parameter has a phase transition in the large size limit. This is reflected in the asymptotic behavior of the MVOP as the degree tends to infinity. The connection with the scalar orthogonal polynomials allows us to find the limiting behavior of the zeros of the determinant of the MVOP. The zeros tend to a curve $\widetildeΣ_0$ in the complex plane that has a self-intersection. The zeros of the individual entries of the MVOP show a different behavior and we find the limiting zero distribution of the upper right entry under a geometric condition on the curve $\widetildeΣ_0$ that we were unable to prove, but that is convincingly supported by numerical evidence.

preprint2020arXiv

A periodic hexagon tiling model and non-Hermitian orthogonal polynomials

We study a one-parameter family of probability measures on lozenge tilings of large regular hexagons that interpolates between the uniform measure on all possible tilings and a particular fully frozen tiling. The description of the asymptotic behavior can be separated into two regimes: the low and the high temperature regime. Our main results are the computations of the disordered regions in both regimes and the limiting densities of the different lozenges there. For low temperatures, the disorded region consists of two disjoint ellipses. In the high temperature regime the two ellipses merge into a single simply connected region. At the transition from the low to the high temperature a tacnode appears. The key to our asymptotic study is a recent approach introduced by Duits and Kuijlaars providing a double integral representation for the correlation kernel. One of the factors in the integrand is the Christoffel-Darboux kernel associated to polynomials that satisfy non-Hermitian orthogonality relations with respect to a complex-valued weight on a contour in the complex plane. We compute the asymptotic behavior of these orthogonal polynomials and the Christoffel-Darboux kernel by means of a Riemann-Hilbert analysis. After substituting the resulting asymptotic formulas into the double integral we prove our main results by classical steepest descent arguments.

preprint2020arXiv

A vector equilibrium problem for symmetrically located point charges on a sphere

We study the equilibrium measure on the two dimensional sphere in the presence of an external field generated by r+1 equal point charges that are symmetrically located around the north pole. The support of the equilibrium measure is known as the droplet. The droplet has a motherbody which we characterize by means of a vector equilibrium problem (VEP) for r measures in the complex plane. The model undergoes two transitions which is reflected in the support of the first component of the minimizer of the VEP, namely the support can be a finite interval containing 0, the union of two intervals, or the full half-line. The two interval case corresponds to a droplet with two disjoint components, and it is analyzed by means of a genus one Riemann surface.

preprint2019arXiv

A numerical method for oscillatory integrals with coalescing saddle points

The value of a highly oscillatory integral is typically determined asymptotically by the behaviour of the integrand near a small number of critical points. These include the endpoints of the integration domain and the so-called stationary points or saddle points -- roots of the derivative of the phase of the integrand -- where the integrand is locally non-oscillatory. Modern methods for highly oscillatory quadrature exhibit numerical issues when two such saddle points coalesce. On the other hand, integrals with coalescing saddle points are a classical topic in asymptotic analysis, where they give rise to uniform asymptotic expansions in terms of the Airy function. In this paper we construct Gaussian quadrature rules that remain uniformly accurate when two saddle points coalesce. These rules are based on orthogonal polynomials in the complex plane. We analyze these polynomials, prove their existence for even degrees, and describe an accurate and efficient numerical scheme for the evaluation of oscillatory integrals with coalescing saddle points.

preprint2019arXiv

The two periodic Aztec diamond and matrix valued orthogonal polynomials

We analyze domino tilings of the two-periodic Aztec diamond by means of matrix valued orthogonal polynomials that we obtain from a reformulation of the Aztec diamond as a non-intersecting path model with periodic transition matrices. In a more general framework we express the correlation kernel for the underlying determinantal point process as a double contour integral that contains the reproducing kernel of matrix valued orthogonal polynomials. We use the Riemann-Hilbert problem to simplify this formula for the case of the two-periodic Aztec diamond. In the large size limit we recover the three phases of the model known as solid, liquid and gas. We describe fine asymptotics for the gas phase and at the cusp points of the liquid-gas boundary, thereby complementing and extending results of Chhita and Johansson.

preprint2010arXiv

Critical behavior of non-intersecting Brownian motions at a tacnode

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é II equation. We show that this Painlevé 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.

preprint2010arXiv

Random matrix model with external source and a constrained vector equilibrium problem

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

preprint2009arXiv

Multiple orthogonal polynomial ensembles

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.

preprint2008arXiv

Painleve IV asymptotics for orthogonal polynomials with respect to a modified Laguerre weight

We study polynomials that are orthogonal with respect to the modified Laguerre weight $z^{-n + ν} e^{-Nz} (z-1)^{2b}$ in the limit where $n, N \to \infty$ with $N/n \to 1$ and $ν$ 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 $Ψ$-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.

preprint2007arXiv

Universality in unitary random matrix ensembles when the soft edge meets the hard edge

Unitary random matrix ensembles Z_{n,N}^{-1} (\det M)^alpha exp(-N Tr V(M)) dM defined on positive definite matrices M, where alpha > -1 and V is real analytic, have a hard edge at 0. The equilibrium measure associated with V typically vanishes like a square root at soft edges of the spectrum. For the case that the equilibrium measure vanishes like a square root at 0, we determine the scaling limits of the eigenvalue correlation kernel near 0 in the limit when n, N tend to infinity such that n/N - 1 = O(n^{-2/3}). For each value of alpha > -1 we find a one-parameter family of limiting kernels that we describe in terms of the Hastings-McLeod solution of the Painleve II equation with parameter alpha + 1/2.

preprint2006arXiv

Large n limit of Gaussian random matrices with external source, Part III: Double scaling limit

We consider the double scaling limit in the random matrix ensemble with an external source $\frac{1}{Z_n} e^{-n \Tr({1/2}M^2 -AM)} dM$ defined on $n\times n$ Hermitian matrices, where $A$ is a diagonal matrix with two eigenvalues $\pm a$ of equal multiplicities. The value $a=1$ is critical since the eigenvalues of $M$ accumulate as $n \to \infty$ on two intervals for $a > 1$ and on one interval for $0 < a < 1$. These two cases were treated in Parts I and II, where we showed that the local eigenvalue correlations have the universal limiting behavior known from unitary random matrix ensembles. For the critical case $a=1$ new limiting behavior occurs which is described in terms of Pearcey integrals, as shown by Brézin and Hikami, and Tracy and Widom. We establish this result by applying the Deift/Zhou steepest descent method to a $3 \times 3$-matrix valued Riemann-Hilbert problem which involves the construction of a local parametrix out of Pearcey integrals. We resolve the main technical issue of matching the local Pearcey parametrix with a global outside parametrix by modifying an underlying Riemann surface.

preprint2005arXiv

Universality of the double scaling limit in random matrix models

We study unitary random matrix ensembles in the critical case where the limiting mean eigenvalue density vanishes quadratically at an interior point of the support. We establish universality of the limits of the eigenvalue correlation kernel at such a critical point in a double scaling limit. The limiting kernels are constructed out of functions associated with the second Painleve equation. This extends a result of Bleher and Its for the special case of a critical quartic potential. The two main tools we use are equilibrium measures and Riemann-Hilbert problems. In our treatment of equilibrium measures we allow a negative density near the critical point, which enables us to treat all cases simultaneously. The asymptotic analysis of the Riemann-Hilbert problem is done with the Deift/Zhou steepest descent analysis. For the construction of a local parametrix at the critical point we introduce a modification of the approach of Baik, Deift, and Johansson so that we are able to satisfy the required jump properties exactly.

preprint2004arXiv

Large $n$ limit of Gaussian random matrices with external source, part II

We continue the study of the Hermitian random matrix ensemble with external source $\frac{1}{Z_n} e^{-n \Tr({1/2}M^2 -AM)} dM$ where $A$ has two distinct eigenvalues $\pm a$ of equal multiplicity. This model exhibits a phase transition for the value $a=1$, since the eigenvalues of $M$ accumulate on two intervals for $a > 1$, and on one interval for $0 < a < 1$. The case $a > 1$ was treated in part I, where it was proved that local eigenvalue correlations have the universal limiting behavior which is known for unitarily invariant random matrices, that is, limiting eigenvalue correlations are expressed in terms of the sine kernel in the bulk of the spectrum, and in terms of the Airy kernel at the edge. In this paper we establish the same results for the case $0 < a < 1$. As in part I we apply the Deift/Zhou steepest descent analysis to a $3 \times 3$-matrix Riemann-Hilbert problem. Due to the different structure of an underlying Riemann surface, the analysis includes an additional step involving a global opening of lenses, which is a new phenomenon in the steepest descent analysis of Riemann-Hilbert problems.