Multivariate Orthogonal Laurent Polynomials and Integrable Systems

Abstract

An ordering for Laurent polynomials in the algebraic torus $({\mathbb{C}}^{\ast }{)}^D$, inspired by the Cantero–Moral–Vel´azquez approach to orthogonal Laurent polynomials in the unit circle, leads to the construction of a moment matrix for a given Borel measure in the unit torus ${\mathbb{T}}^D$. The Gauss–Borel factorization of this moment matrix allows for the construction of multivariate biorthogonal Laurent polynomials in the unit torus, which can be expressed as last quasi-determinants of bordered truncations of the moment matrix. The associated second-kind functions are expressed in terms of the Fourier series of the given measure. Persymmetries and partial persymmetries of the moment matrix are studied and Cauchy integral representations of the second-kind functions are found, as well as Plemelj-type formulae. Spectral matrices give string equations for the moment matrix, which model the three-term relations as well as the Christoffel–Darboux formulae.

Christoffel-type perturbations of the measure given by the multiplication by Laurent polynomials are studied. Sample matrices on poised sets of nodes, which belong to the algebraic hypersurface of the perturbing Laurent polynomial, are used to find a Christoffel formula that expresses the perturbed orthogonal Laurent polynomials in terms of a last quasi-determinant of a bordered sample matrix constructed in terms of the original orthogonal Laurent polynomials. Poised sets exist only for prepared Laurent polynomials, which are analyzed from the perspective of Newton polytopes and tropical geometry. Then, an algebraic geometrical characterization of prepared Laurent polynomial perturbation and poised sets is given; full-column-rankness of the corresponding multivariate Laurent–Vandermonde matrices and a product of different prime prepared Laurent polynomials leads to such sets. Some examples are constructed in terms of perturbations of the Lebesgue–Haar measure.

Discrete and continuous deformations of the measure lead to a Toda-type integrable hierarchy, being the corresponding flows described through Lax and Zakharov–Shabat equations; bilinear equations and vertex operators are found. Varying size matrix nonlinear partial difference and differential equations of two-dimensional Toda lattice type are shown to be solved by matrix coefficients of the multivariate orthogonal polynomials. The discrete flows are connected with a Gauss–Borel factorization of the Jacobi-type matrices and its quasi-determinants allow for expressions for the multivariate orthogonal polynomials in terms of shifted quasi-tau matrices, which generalize those that relate the Baker functions with ratios of Miwa shifted $\tau$-functions in the one-dimensional scenario. It is shown that the discrete and continuous flows are deeply connected and determine nonlinear partial difference-differential equations that involve only one site in the integrable lattice behaving as a Kadomtsev–Petviashvili-type system.