BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Arno B. J. Kuijlaars

Arno B. J. Kuijlaars appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPmath.CAmath.PRmath.CVcond-mat.stat-mechmath.FAmath.NANumerical Analysis

Catalog footprint

What is connected

28works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The normal matrix model with a monomial potential, a vector equilibrium problem, and multiple orthogonal polynomials on a star

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 $(μ_{1},\ldots,μ_{d})$ supported on star-like sets in the plane. We show that in the subcritical regime, the first component $μ_{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 $μ_{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$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Singular values of products of random matrices and polynomial ensembles

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large Deviations for a Non-Centered Wishart Matrix

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Singular values of products of Ginibre random matrices, multiple orthogonal polynomials and hard edge scaling limits

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Weakly Admissible Vector Equilibrium Problems

We establish lower semi-continuity and strict convexity of the energy functionals for a large class of vector equilibrium problems in logarithmic potential theory. This in particular implies the existence and uniqueness of a minimizer for such vector equilibrium problems. Our work extends earlier results in that we allow unbounded supports without having strongly confining external fields. To deal with the possible noncompactness of supports, we map the vector equilibrium problem onto the Riemann sphere and our results follow from a study of vector equilibrium problems on compacts in higher dimensions. Our results cover a number of cases that were recently considered in random matrix theory and for which the existence of a minimizer was not clearly established yet.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Double scaling limit for modified Jacobi-Angelesco polynomials

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multiple Meixner-Pollaczek polynomials and the six-vertex model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-intersecting squared Bessel paths with one positive starting and ending point

We consider a model of $n$ non-intersecting squared Bessel processes with one starting point $a>0$ at time t=0 and one ending point $b>0$ at time $t=T$. After proper scaling, the paths fill out a region in the $tx$-plane. Depending on the value of the product $ab$ the region may come to the hard edge at 0, or not. We formulate a vector equilibrium problem for this model, which is defined for three measures, with upper constraints on the first and third measures and an external field on the second measure. It is shown that the limiting mean distribution of the paths at time $t$ is given by the second component of the vector that minimizes this vector equilibrium problem. The proof is based on a steepest descent analysis for a $4 \times 4$ matrix valued Riemann-Hilbert problem which characterizes the correlation kernel of the paths at time $t$. We also discuss the precise locations of the phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Orthogonal polynomials in the normal matrix model with a cubic potential

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 $(μ_1,μ_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 $(μ_1,μ_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 $μ_1$, the first component of the equilibrium measure. We also obtain analytical results for the measure $μ_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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multiple orthogonal polynomials in random matrix theory

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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Hermitian two matrix model with an even quartic potential

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access