Abstract:We establish a new perturbation theory for orthogonal polynomials using a Riemann--Hilbert approach and consider applications in numerical linear algebra and random matrix theory. This new approach shows that the orthogonal polynomials with respect to two measures can be effectively compared using the difference of their Stieltjes transforms on a suitably chosen contour. Moreover, when two measures are close and satisfy some regularity conditions, we use the theta functions of a hyperelliptic Riemann surface to derive explicit and accurate expansion formulae for the perturbed orthogonal polynomials. In contrast to other approaches, a key strength of the methodology is that estimates can remain valid as the degree of the polynomial grows. The results are applied to analyze several numerical algorithms from linear algebra, including the Lanczos tridiagonalization procedure, the Cholesky factorization and the conjugate gradient algorithm. As a case study, we investigate these algorithms applied to a general spiked sample covariance matrix model by considering the eigenvector empirical spectral distribution and its limits. For the first time, we give precise estimates on the output of the algorithms, applied to this wide class of random matrices, as the number of iterations diverges. In this setting, beyond the first order expansion, we also derive a new mesoscopic central limit theorem for the associated orthogonal polynomials and other quantities relevant to numerical algorithms.

On a Class of Elliptic Orthogonal Polynomials and their Integrability

On a class of elliptic orthogonal polynomials and their integrability

Determinantal approach to multiple orthogonal polynomials, and the corresponding integrable equations

Superelliptic Affine Lie algebras and orthogonal polynomials

Semiclassical orthogonal polynomials on the unit circle: A Riemann-Hilbert perspective

Orthogonal Polynomials on Planar Cubic Curves

A Riemann--Hilbert Approach to Skew-Orthogonal Polynomials of Symplectic Type

Hypergeometric orthogonal polynomials of Jacobi type

Integral and hypergeometric representations for multiple orthogonal polynomials

A Riemann-Hilbert Approach to Skew-Orthogonal Polynomials of Symplectic Type

Elliptic Curves and a New Construction of Integrable Systems

On Certain Wronskians of Multiple Orthogonal Polynomials.

On the -Classical Orthogonal Polynomials

Skew-Orthogonal Polynomials and Pfaff Lattice Hierarchy Associated With an Elliptic Curve

Classical discrete multiple orthogonal polynomials: hypergeometric and integral representations

Introduction to the theory of elliptic hypergeometric integrals

Orthogonal Trigonometric Polynomials: Riemann-Hilbert Analysis and Relations with OPUC

Connecting Exceptional Orthogonal Polynomials of Different Kind

Orthogonal trigonometric polynomials from Riemann-Hilbert view

A Riemann--Hilbert approach to the perturbation theory for orthogonal polynomials: Applications to numerical linear algebra and random matrix theory

Applications of multiple orthogonal polynomials with hypergeometric moment generating functions