Abstract:The joint moments of the derivatives of the characteristic polynomial of a random unitary matrix, and also a variant of the characteristic polynomial that is real on the unit circle, in the large matrix size limit, have been studied intensively in the past twenty five years, partly in relation to conjectural connections to the Riemann zeta-function and Hardy's function. We completely settle the most general version of the problem of convergence of these joint moments, after they are suitably rescaled, for an arbitrary number of derivatives and with arbitrary positive real exponents. Our approach relies on a hidden, higher-order exchangeable structure, that of an exchangeable array. Using these probabilistic techniques, we then give a combinatorial formula for the leading order coefficient in the asymptotics of the joint moments, when the power on the characteristic polynomial itself is a positive real number and the exponents of the derivatives are integers, in terms of a finite number of finite-dimensional integrals which are explicitly computable. Finally, we develop a method, based on a class of Hankel determinants shifted by partitions, that allows us to give an exact representation of all these joint moments, for finite matrix size, in terms of derivatives of Painlevé V transcendents, and then for the leading order coefficient in the large-matrix limit in terms of derivatives of solutions of the $\sigma$-Painlevé III' equation. Equivalently, we can represent all the joint moments of power sum linear statistics of a certain determinantal point process behind this problem in terms of derivatives of $\sigma$-Painlevé III' transcendents. This gives an efficient way to compute all these quantities explicitly. Our methods can be used to obtain analogous results for a number of other models sharing the same features.

Construction of Effective Randomized Projective Estimates for Solutions of Integral Equations Based on Legendre Polynomials

Conditional Optimization of Algorithms for Estimating Distributions of Solutions to Stochastic Differential Equations

Randomized Iterative Methods for Linear Systems

Optimal Lagrange Interpolation Projectors and Legendre Polynomials

The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers

A Bregman-Kaczmarz method for nonlinear systems of equations

Recursive construction of biorthogonal polynomials for handling polynomial regression

Improved Laguerre Spectral Methods with Less Round-off Errors and Better Stability

Randomized Matrix Approximation to Enhance Regularized Projection Schemes in Inverse Problems

Robust Projection Method for Fredholm Integral Equation of the Second Kind

Exchangeable arrays and integrable systems for characteristic polynomials of random matrices

Regular polynomial interpolation and approximation of global solutions of linear partial differential equations

Rational approximation of Euler's constant using multiple orthogonal polynomials

Efficient Numerical Integration in Reproducing Kernel Hilbert Spaces via Leverage Scores Sampling

Maximum Likelihood Estimation of Nonlinear Rational Expectations Models by Orthogonal Polynomial Projection Methods

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

Randomized Methods for Linear Constraints: Convergence Rates and Conditioning

Numerical solution of a weakly singular integral equation by the method of orthogonal polynomials

Projection and modified projection methods for nonlinear Hammerstein integral equations on the real line using Hermite polynomials

A Riemann--Hilbert approach to computing the inverse spectral map for measures supported on disjoint intervals

Optimal Randomized First-Order Methods for Least-Squares Problems