On the Christoffel–Darboux formula for generalized matrix orthogonal polynomials

We obtain an extension of the Christoffel–Darboux formula for matrix orthogonal polynomials with a generalized Hankel symmetry, including the Adler–van Moerbeke generalized orthogonal polynomials.     

Perturbations on the subdiagonals of Toeplitz matrices

Let L be an Hermitian linear functional defined on the linear space of Laurent polynomials. It is very well known that the Gram matrix of the associated bilinear functional in the linear space of polynomials is a Toeplitz matrix. In this contribution we analyze some linear spectral transforms of L such that the corresponding Toeplitz matrix is the result of the addition of a constant in a subdiagonal of the initial Toeplitz matrix. We focus our attention in the analysis of the quasi-definite character of the perturbed linear functional as well as in the explicit expressions of the new monic orthogonal polynomial sequence in terms of the first one.     

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.     

Szegő–Lobatto quadrature rules

Gauss-type quadrature rules with one or two prescribed nodes are well known and are commonly referred to as Gauss–Radau and Gauss–Lobatto quadrature rules, respectively. Efficient algorithms are available for their computation. Szeg? quadrature rules are analogs of Gauss quadrature rules for the integration of periodic functions; they integrate exactly trigonometric polynomials of as high degree as possible. Szeg? quadrature rules have a free parameter, which can be used to prescribe one node. This paper discusses an analog of Gauss–Lobatto rules, i.e., Szeg? quadrature rules with two prescribed nodes. We refer to these rules as Szeg?–Lobatto rules. Their properties as well as numerical methods for their computation are discussed.     

Riemann‐Hilbert problem and matrix discrete Painlevé II systems

Matrix Szeg? biorthogonal polynomials for quasi‐definite matrices of Hölder continuous weights are studied. A Riemann‐Hilbert problem is uniquely solved in terms of the matrix Szeg? polynomials and its Cauchy transforms. The Riemann‐Hilbert problem is given as an appropriate framework for the discussion of the Szeg? matrix and the associated Szeg? recursion relations for the matrix orthogonal polynomials and its Cauchy transforms. Pearson‐type differential systems characterizing the matrix of weights are studied. These are linear systems of ordinary differential equations that are required to have trivial monodromy. Linear ordinary differential equations for the matrix Szeg? polynomials and its Cauchy transforms are derived. It is shown how these Pearson systems lead to nonlinear difference equations for the Verblunsky matrices and two examples, of Fuchsian and non‐Fuchsian type, are considered. For both cases, a new matrix version of the discrete Painlevé II equation for the Verblunsky matrices is found. Reductions of these matrix discrete Painlevé II systems presenting locality are discussed.     

Direct and inverse polynomial perturbations of hermitian linear functionals

This paper is devoted to the study of direct and inverse (Laurent) polynomial modifications of moment functionals on the unit circle, i.e., associated with hermitian Toeplitz matrices. We present a new approach which allows us to study polynomial modifications of arbitrary degree.The main objective is the characterization of the quasi-definiteness of the functionals involved in the problem in terms of a difference equation relating the corresponding Schur parameters. The results are presented in the general framework of (non-necessarily quasi-definite) hermitian functionals, so that the maximum number of orthogonal polynomials is characterized by the number of consistent steps of an algorithm based on the referred recurrence for the Schur parameters.The non-uniqueness of the inverse problem makes it more interesting than the direct one. Due to this reason, special attention is paid to the inverse modification, showing that different approaches are possible depending on the data about the polynomial modification at hand. These different approaches are translated as different kinds of initial conditions for the related inverse algorithm.Some concrete applications to the study of orthogonal polynomials on the unit circle show the effectiveness of this new approach: an exhaustive and instructive analysis of the functionals coming from a general inverse polynomial perturbation of degree one for the Lebesgue measure; the classification of those pairs of orthogonal polynomials connected by a certain type of linear relation with constant polynomial coefficients; and the determination of those orthogonal polynomials whose associated ones are related to a degree one polynomial modification of the original orthogonality functional.     

Solution of a multiple Nevanlinna–Pick problem for Carathéodory functions via orthogonal rational functions in the matrix case

We consider an interpolation problem of Nevanlinna–Pick type for matrix‐valued Carathéodory functions, where the values of the functions and its derivatives up to certain orders are given at finitely many points of the open unit disk. For the non‐degenerate case, i.e., in the particular situation that a specific block matrix (which is formed by the given data in the problem) is positive Hermitian, the solution set of this problem is described in terms of orthogonal rational matrix‐valued functions. These rational matrix functions play here a similar role as Szegő's orthogonal polynomials on the unit circle in the classical case of the trigonometric moment problem. In particular, we present and use a connection between Szegő and Schur parameters for orthogonal rational matrix‐valued functions which in the primary situation of orthogonal polynomials was found by Geronimus. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Distinguished Solutions of Truncated Matricial Hamburger Moment Problems

We study two slightly different versions of the truncated matricial Hamburger moment problem. A central topic is the construction and investigation of distinguished solutions of both moment problems under consideration. These solutions turn out to be nonnegative Hermitian q × q Borel measures on the real axis which are concentrated on a finite number of points. These points and the corresponding masses will be explicitly described in terms of the given data. Furthermore, we investigate a particular class of sequences (s_j)^∞_{j = 0} of complex q × q matrices for which the corresponding infinite matricial Hamburger moment problem has a unique solution. Our approach is mainly algebraic. It is based on the use of particular matrix polynomials constructed from a nonnegative Hermitian block Hankel matrix. These matrix polynomials are immediate generalizations of the monic orthogonal matrix polynomials associated with a positive Hermitian block Hankel matrix. We generalize a classical theorem due to Kronecker on infinite Hankel matrices of finite rank to block Hankel matrices and discuss its consequences for the nonnegative Hermitian case.     

Matrix measures on the unit circle, moment spaces, orthogonal polynomials and the Geronimus relations

We study the moment space corresponding to matrix measures on the unit circle. Moment points are characterized by non-negative definiteness of block Toeplitz matrices. This characterization is used to derive an explicit representation of orthogonal polynomials with respect to matrix measures on the unit circle and to present a geometric definition of canonical moments. It is demonstrated that these geometrically defined quantities coincide with the Verblunsky coefficients, which appear in the Szegö recursions for the matrix orthogonal polynomials. Finally, we provide an alternative proof of the Geronimus relations which is based on a simple relation between canonical moments of matrix measures on the interval [−1, 1] and the Verblunsky coefficients corresponding to matrix measures on the unit circle.     

Green’s matrices

Our goal is to identify and understand matrices A that share essential properties of the unitary Hessenberg matrices M that are fundamental for Szegö’s orthogonal polynomials. Those properties include: (i) Recurrence relations connect characteristic polynomials {r_k(x)} of principal minors of A. (ii) A is determined by generators (parameters generalizing reflection coefficients of unitary Hessenberg theory). (iii) Polynomials {r_k(x)} correspond not only to A but also to a certain “CMV-like” five-diagonal matrix. (iv) The five-diagonal matrix factors into a product BC of block diagonal matrices with 2 × 2 blocks. (v) Submatrices above and below the main diagonal of A have rank 1. (vi) A is a multiplication operator in the appropriate basis of Laurent polynomials. (vii) Eigenvectors of A can be expressed in terms of those polynomials.Conditions (v) connects our analysis to the study of quasi-separable matrices. But the factorization requirement (iv) narrows it to the subclass of “Green’s matrices” that share Properties (i)-(vii).The key tool is “twist transformations” that provide 2ⁿ matrices all sharing characteristic polynomials of principal minors with A. One such twist transformation connects unitary Hessenberg to CMV. Another twist transformation explains findings of Fiedler who noticed that companion matrices give examples outside the unitary Hessenberg framework. We mention briefly the further example of a Daubechies wavelet matrix. Infinite matrices are included.     

Random block matrices generalizing the classical Jacobi and Laguerre ensembles

In this paper we consider random block matrices which generalize the classical Laguerre ensemble and the Jacobi ensemble. We show that the random eigenvalues of the matrices can be uniformly approximated by the zeros of matrix orthogonal polynomials and obtain a rate for the maximum difference between the eigenvalues and the zeros. This relation between the random block matrices and matrix orthogonal polynomials allows a derivation of the asymptotic spectral distribution of the matrices.     

Polynomial perturbations of bilinear functionals and Hessenberg matrices

This paper deals with symmetric and non-symmetric polynomial perturbations of symmetric quasi-definite bilinear functionals. We establish a relation between the Hessenberg matrices associated with the initial and the perturbed functionals using LU and QR factorizations. Moreover we give an explicit algebraic relation between the sequences of orthogonal polynomials associated with both functionals.     

Biorthogonal Laurent Polynomials,Toplitz Determinants,Minimal Toda Orbits and Isomonodromic Tau Functions

We consider the class of biorthogonal polynomials that are used to solve the inverse spectral problem associated to elementary co-adjoint orbits of the Borel group of upper triangular matrices; these orbits are the phase space of generalized integrable lattices of Toda type. Such polynomials naturally interpolate between the theory of orthogonal polynomials on the line and orthogonal polynomials on the unit circle and tie together the theory of Toda, relativistic Toda, Ablowitz-Ladik and Volterra lattices. We establish corresponding Christoffel-Darboux formulae. For all these classes of polynomials a 2 × 2 system of Differential-Difference-Deformation equations is analyzed in the most general setting of pseudo-measures with arbitrary rational logarithmic derivative. They provide particular classes of isomonodromic deformations of rational connections on the Riemann sphere. The corresponding isomonodromic tau function is explicitly related to the shifted Toplitz determinants of the moments of the pseudo-measure. In particular, the results imply that any (shifted) Toplitz (Hankel) determinant of a symbol (measure) with arbitrary rational logarithmic derivative is an isomonodromic tau function.     

On differential properties for bivariate orthogonal polynomials

In this paper, we consider bivariate orthogonal polynomials associated with a quasi-definite moment functional which satisfies a Pearson-type partial differential equation. For these polynomials differential properties are obtained. In particular, we deduce some structure and orthogonality relations for the successive partial derivatives of the polynomials.      

Szegő's theorem on Parreau–Widom sets

In this paper, we generalize Szeg?'s theorem for orthogonal polynomials on the real line to infinite gap sets of Parreau–Widom type. This notion includes Cantor sets of positive measure. The Szeg? condition involves the equilibrium measure which in turn is absolutely continuous. Our approach builds on a canonical factorization of the M-function and the covering space formalism of Sodin–Yuditskii.     

On five-diagonal Toeplitz matrices and orthogonal polynomials on the unit circle

This paper deals with modifications of the Lebesgue moment functional by trigonometric polynomials of degree 2 and their associated orthogonal polynomials on the unit circle. We use techniques of five-diagonal matrix factorization and matrix polynomials to study the existence of such orthogonal polynomials.Dedicated to Prof. Luigi Gatteschi on his 70th birthdayThis research was partially supported by Diputación General de Aragón under grant P CB-12/91.     

Biorthogonal Rational Functions and the Generalized Eigenvalue Problem

We present some general results concerning so-called biorthogonal polynomials of R_II type introduced by M. Ismail and D. Masson. These polynomials give rise to a pair of rational functions which are biorthogonal with respect to a linear functional. It is shown that these rational functions naturally appear as eigenvectors of the generalized eigenvalue problem for two arbitrary tri-diagonal matrices. We study spectral transformations of these functions leading to a rational modification of the linear functional. An analogue of the Christoffel–Darboux formula is obtained.     

Orthogonal Newton polynomials

The problem is to determine all nonnegative measures on the Borel subsets of the complex plane with respect to which all polynomials are square integrable and with respect to which the Newton polynomials form an orthogonal set. The Newton polynomials do not belong to any classical scheme of orthogonal polynomials. The discovery that a plane measure exists with respect to which they form an orthogonal set was only recently made by T. L. Kriete and D. Trutt [Amer. J. Math.93 (1971), 215–225]. A general structure theory for such measures is now obtained under hypotheses suggested by the expansion theory of Cesàro operators.     

Bi-orthogonal systems on the unit circle, regular semi-classical weights and integrable systems — II

We derive the Christoffel–Geronimus–Uvarov transformations of a system of bi-orthogonal polynomials and associated functions on the unit circle, that is to say the modification of the system corresponding to a rational modification of the weight function. In the specialisation of the weight function to the regular semi-classical case with an arbitrary number of regular singularities {z₁,…,z_M} the bi-orthogonal system is known to be monodromy preserving with respect to deformations of the singular points. If the zeros and poles of the Christoffel–Geronimus–Uvarov factors coincide with the singularities then we have the Schlesinger transformations of this isomonodromic system. Compatibility of the Schlesinger transformations with the other structures of the system — the recurrence relations, the spectral derivatives and deformation derivatives is explicitly deduced. Various forms of Hirota–Miwa equations are derived for the τ-functions or equivalently Toeplitz determinants of the system.