Search

In this work, we develop two methods to obtain a factorisation as a product of bidiagonal matrices for the Hessenberg recurrence matrix of a system of multiple orthogonal polynomials. One method is based on the Gauss-Borel factorisation of the moment matrix of the system and of its Darboux transformations, and the other uses the connection of multiple orthogonal polynomials with production matrices and branched continued fractions. As a case study, we present an explicit bidiagonal factorisation for the Hessenberg recurrence matrix of the Jacobi-Pi\~neiro polynomials, using both methods, and of the multiple Laguerre polynomials of first kind, as a limiting case., Comment: 22 pages, 1 figure

General Geronimus transformations, defined by regular matrix polynomials that are neither required to be monic nor restricted by the rank of their leading coefficients, are applied through both right and left multiplication to a rectangular matrix of measures associated with mixed multiple orthogonal polynomials. These transformations produce Christoffel-type formulas that establish relationships between the perturbed and original polynomials. Moreover, it is proven that the existence of Geronimus-perturbed orthogonality is equivalent to the non-cancellation of certain $\tau$-determinants. The effect of these transformations on the Markov-Stieltjes matrix functions is also determined. As a case study, we examine the Jacobi-Pi\~neiro orthogonal polynomials with three weights., Comment: 40 pages

Mixed orthogonal Laurent polynomials on the unit circle of CMV type are constructed utilizing a matrix of moments and its Gauss--Borel factorization and employing a multiple extension of the CMV ordering. A systematic analysis of the associated multiple orthogonality and biorthogonality relations, and an examination of the degrees of the Laurent polynomials is given. Recurrence relations, expressed in terms of banded matrices, are found. These recurrence relations lay the groundwork for corresponding Christoffel-Darboux kernels and relations, as well as for elucidating the ABC theorem. The paper also develops the theory of diagonal Christoffel and Geronimus perturbations of the matrix of measures. Christoffel formulas are found for both perturbations., Comment: 48 pages

This paper demonstrates how to explicitly construct a bidiagonal factorization of the banded recurrence matrix that appears in mixed multiple orthogonality on the step-line in terms of the coeffcients of the mixed multiple orthogonal polynomials. The construction is based on the \(LU\) factorization of the moment matrix and Christoffel transformations applied to the matrix of measures and the associated mixed multiple orthogonal polynomials., Comment: 14 pages

This work explores classical discrete multiple orthogonal polynomials, including Hahn, Meixner of the first and second kinds, Kravchuk, and Charlier polynomials, with an arbitrary number of weights. Explicit expressions for the recursion coefficients of Hahn multiple orthogonal polynomials are derived. By leveraging the multiple Askey scheme and the recently discovered explicit hypergeometric representation of type I multiple Hahn polynomials, corresponding explicit hypergeometric representations are provided for the type I polynomials and recursion coefficients of all the aforementioned descendants within the Askey scheme. Additionally, integral representations for these families within the Hahn class in the Askey scheme are presented. The multiple Askey scheme is further completed by providing the corresponding limits for the weights, polynomials, and recurrence coefficients., Comment: 37 pages, 5 figures

This paper addresses two primary objectives in the realm of classical multiple orthogonal polynomials with an arbitrary number of weights. Firstly, it establishes new and explicit hypergeometric expressions for type I Hahn multiple orthogonal polynomials. Secondly, applying the residue theorem and the Mellin transform, the paper derives contour integral representations for several families of orthogonal polynomials. Specifically, it presents contour integral formulas for both type I and type II multiple orthogonal polynomials in the Laguerre of the first kind, Jacobi-Pi\~neiro, and Hahn families. The evaluation of these integrals leads to explicit hypergeometric representations., Comment: 21 pages

In this paper, we extend our investigation into semiclassical multiple discrete orthogonal polynomials by considering an arbitrary number of weights. We derive multiple versions of the Toda equations and the Laguerre-Freud equations for the multiple generalized Charlier and multiple generalized Meixner II families., Comment: 23 pages

Performing both right and left multiplication operations using general regular matrix polynomials, which need not be monic and may possess leading coefficients of arbitrary rank, on a rectangular matrix of measures associated with mixed multiple orthogonal polynomials, reveals corresponding Christoffel formulas. These formulas express the perturbed mixed multiple orthogonal polynomials in relation to the original ones. Utilizing the divisibility theorem for matrix polynomials, we establish a criterion for the existence of perturbed orthogonality, expressed through the non-cancellation of certain $\tau$ determinants., Comment: 32 pages. Some minor typos corrected in the second version

This paper serves as an introduction to banded totally positive matrices, exploring various characterizations and associated properties. A significant result within is the demonstration that the collection of such matrices forms a semigroup, notably including a subset permitting positive bidiagonal factorization. Moreover, the paper applies this concept to investigate step line normality concerning the degrees of associated recursion polynomials. It presents a spectral Favard theorem, ensuring the existence of measures, thereby guaranteeing that these recursion polynomials represent mixed multiple orthogonal polynomials that maintain normality on the step line indices.

This paper delves into classical multiple orthogonal polynomials with an arbitrary number of weights, including Jacobi-Pi\~neiro, Laguerre of both first and second kinds, as well as multiple orthogonal Hermite polynomials. Novel explicit expressions for nearest-neighbor recurrence coefficients, as well as the step line case, are provided for all these polynomial families. Furthermore, new explicit expressions for type I multiple orthogonal polynomials are derived for Laguerre of the second kind and also for multiple Hermite polynomials.

Additive perturbations, specifically, matrix Uvarov transformations for matrix orthogonal polynomials, are under consideration. Christoffel-Uvarov formulas are deduced for the perturbed biorthogonal families, along with their matrix norms. These formulations are expressed in terms of the spectral jets of the Christoffel-Darboux kernels., Comment: arXiv admin note: substantial text overlap with arXiv:1605.04617

For a general number $p\geq 2$ of measures, we provide explicit expressions for the Jacobi-Pi\~neiro and Laguerre of the first kind multiple orthogonal polynomials of type I, presented in terms of multiple hypergeometric functions.

This paper explores a factorization using bidiagonal matrices of the recurrence matrix of Hahn multiple orthogonal polynomials. The factorization is expressed in terms of ratios involving the generalized hypergeometric function ${}_3F_2$ and is proven using recently discovered contiguous relations. Moreover, employing the multiple Askey scheme, a bidiagonal factorization is derived for the Hahn descendants, including Jacobi-Pi\~neiro, multiple Meixner (kinds I and II), multiple Laguerre (kinds I and II), multiple Kravchuk, and multiple Charlier, all represented in terms of hypergeometric functions. For the cases of multiple Hahn, Jacobi-Pi\~neiro, Meixner of kind II, and Laguerre of kind I, where there exists a region where the recurrence matrix is nonnegative, subregions are identified where the bidiagonal factorization becomes a positive bidiagonal factorization., Comment: 14 pages, 2 figures

This paper investigates stochastic finite matrices and the corresponding finite Markov chains constructed using recurrence matrices for general families of orthogonal polynomials and multiple orthogonal polynomials. The paper explores the spectral theory of transition matrices, utilizing both orthogonal and multiple orthogonal polynomials. Several properties are derived, including classes, periodicity, recurrence, stationary states, ergodicity, expected recurrence times, time-reversed chains, and reversibility. Furthermore, the paper uncovers factorization in terms of pure birth and pure death processes. The case study focuses on hypergeometric orthogonal polynomials, where all the computations can be carried out effectively. Particularly within the Askey scheme, all descendants under Hahn (excluding Bessel), such as Hahn, Jacobi, Meixner, Kravchuk, Laguerre, Charlier, and Hermite, present interesting examples of recurrent reversible birth and death finite Markov chains. Additionally, the paper considers multiple orthogonal polynomials, including multiple Hahn, Jacobi-Pi\~neiro, Laguerre of the first kind, and Meixner of the second kind, along with their hypergeometric representations and derives the corresponding recurrent finite Markov chains and time-reversed chains., Comment: 40 pages, 3 figures

Banded bounded matrices, which represent non normal operators, of oscillatory type that admit a positive bidiagonal factorization are considered. To motivate the relevance of the oscillatory character the Favard theorem for Jacobi matrices is revisited and it is shown that after an adequate shift of the Jacobi matrix one gets an oscillatory matrix. In this work we present a spectral theorem for this type of operators and show how the theory of multiple orthogonal polynomials apply., Comment: arXiv admin note: substantial text overlap with arXiv:2203.13578

In this paper, the authors investigate the case of discrete multiple orthogonal polynomials with two weights on the step line, which satisfy Pearson equations. The discrete multiple orthogonal polynomials in question are expressed in terms of tau functions, which are double Wronskians of generalized hypergeometric series. The shifts in the spectral parameter for type II and type I multiple orthogonal polynomials are described using banded matrices. It is demonstrated that these polynomials offer solutions to multicomponent integrable extensions of the nonlinear Toda equations. Additionally, the paper characterizes extensions of the Nijhoff-Capel totally discrete Toda equations. The hypergeometric $\tau$-functions are shown to provide solutions to these integrable nonlinear equations. Furthermore, the authors explore Laguerre-Freud equations, nonlinear equations for the recursion coefficients, with a particular focus on the multiple Charlier, generalized multiple Charlier, multiple Meixner II, and generalized multiple Meixner II cases., Comment: 31 pages

In this work we show how to get advantage from the Riemann--Hilbert analysis in order to obtain information about the matrix orthogonal polynomials and functions of second kind associated with a weight matrix. We deduce properties for the recurrence relation coefficients from differential properties of the weight matrix. We take the matrix polynomials of Hermite, Laguerre and Jacobi type as a case study.

In this work we survey on connections of Markov chains and the theory of multiple orthogonality. Here we mainly concentrate on give a procedure to generate stochastic tetra diagonal Hessenberg matrices, coming from some specific families of multiple orthogonal, such as the ones of Jacobi--Pi\~neiro and Hypergeometric Lima--Loureiro. We show that associated with a positive tetra diagonal nonnegative bounded Hessenberg matrix we can construct two stochastic tetra diagonal ones. These two stochastic tridiagonal nonnegative Hessenberg matrices are shown to be, enlightened by the Poincar\'e theorem, limit transpose of each other.

Explicit expressions for the Hahn multiple polynomials of type I, in terms of Kamp\'e de F\'eriet hypergeometric series, are given. Orthogonal and biorthogonal relations are proven. Then, part of the Askey scheme for multiple orthogonal polynomials type I is completed. In particular, explicit expressions in terms of generalized hypergeometric series and Kamp\'e de F\'eriet hypergeometric series, are given for the multiple orthogonal polynomials of type I for the Jacobi-Pi\~neiro, Meixner I, Meixner II, Kravchuk, Laguerre I, Laguerre II and Charlier families., Comment: 27 pages. Minor typos corrected

Spectral and factorization properties of oscillatory matrices leads to a spectral Favard theorem for bounded banded matrices, that admit a positive bidiagonal factorization, in terms of sequences of mixed multiple orthogonal polynomials with respect to a set positive Lebesgue-Stieltjes measures. A mixed multiple Gauss quadrature formula with corresponding degrees of precision is given., Comment: 33 pages

The Cholesky factorization of the moment matrix is considered for discrete orthogonal polynomials of hypergeometrical type. We derive the Laguerre-Freud equations when the first moments of the weights are given by the ${}_1F_2$, ${}_2F_2$ and ${}_3F_2$ generalized hypergeometrical functions., Comment: 19 pages. arXiv admin note: substantial text overlap with arXiv:2107.02177, arXiv:2107.01747

Recently a spectral Favard theorem for bounded banded lower Hessenberg matrices that admit a positive bidiagonal factorization was presented. In this paper conditions, in terms of continued fractions, for an oscillatory tetradiagonal Hessenberg matrix to have such positive bidiagonal factorization are found. Oscillatory tetradiagonal Toeplitz matrices are taken as a case study of matrix that admits a positive bidiagonal factorization. Moreover, it is proved that oscillatory banded Hessenberg matrices are organized in rays, with the origin of the ray not having the positive bidiagonal factorization and all the interior points of the ray having such positive bidiagonal factorization., Comment: 20 pages This is the second part of the splitting of the paper arXiv:2203.13578 into three

Recently a spectral Favard theorem for bounded banded lower Hessenberg matrices that admit a positive bidiagonal factorization was presented. These type of matrices are oscillatory. In this paper the Lima-Loureiro hypergeometric multiple orthogonal polynomials and the Jacobi-Pi\~neiro multiple orthogonal polynomials are discussed at the light of this bidiagonal factorization for tetradiagonal matrices. The Darboux transformations of tetradiagonal Hessenberg matrices is studied and Christoffel formulas for the elements of the bidiagonal factorization are given, i.e., the bidiagonal factorization is given in terms of the recursion polynomials evaluated at the origin., Comment: This is the third part of the splitting of the paper arXiv:2203.13578 into three. 15 pages and 1 figure

We consider matrix orthogonal polynomials related to Jacobi type matrices of weights that can be defined in terms of a given matrix Pearson equation. Stating a Riemann-Hilbert problem we can derive first and second order differential relations that these matrix orthogonal polynomials and the second kind functions associated to them verify. For the corresponding matrix recurrence coefficients, non-Abelian extensions of a family of discrete Painlev\'e d-PIV equations are obtained for the three term recurrence relation coefficients.

A spectral Favard theorem for bounded banded lower Hessenberg matrices that admit a positive bidiagonal factorization is found. The large knowledge on the spectral and factorization properties of oscillatory matrices leads to this spectral Favard theorem in terms of sequences of multiple orthogonal polynomials of types I and II with respect to a set of positive Lebesgue-Stieltjes~measures. Also a multiple Gauss quadrature is proven and corresponding degrees of precision are found. This spectral Favard theorem is applied to Markov chains with $(p+2)$-diagonal transition matrices, i.e. beyond birth and death, that admit a positive stochastic bidiagonal factorization. In the finite case, the Karlin-McGregor spectral representation is given. It is shown that the Markov chains are recurrent and explicit expressions in terms of the orthogonal polynomials for the stationary distributions are given. Similar results are obtained for the countable infinite Markov chain. Now the Markov chain is not necessarily recurrent, and it is characterized in terms of the first measure. Ergodicity of the Markov chain is discussed in terms of the existence of a mass at $1$, which is an eigenvalue corresponding to the right and left eigenvectors., Comment: 38 pages. This a very much improved version of the initial one. We have divided the original paper in three parts. Two more parts to follow: arXiv:2210.10728, Positive bidiagonal factorization of tetradiagonal Hessenberg matrices arXiv:2210.10727, Bidiagonal factorization of tetradiagonal matrices and Darboux transformations

Contiguous hypergeometric relations for semiclassical discrete orthogonal polynomials are described as Christoffel and Geronimus transformations. Using the Christoffel-Geronimus-Uvarov formulas quasi-determinatal expressions for the shifted semiclassical discrete orthogonal polynomials are obtained.

The Cholesky factorization of the moment matrix is considered for the generalized Charlier, generalized Meixner, and Gauss hypergeometric discrete orthogonal polynomials. For the generalized Charlier, we present an alternative derivation of the Laguerre-Freud relations found by Smet and Van Assche. Third-order and second-order nonlinear ordinary differential equations are found for the recursion coefficient $\gamma_n$, that happen to be forms of the Painlev\'e $\text{deg-P}_{\text V}$ in disguise. Laguerre-Freud relations are also found for the generalized Meixner case, which are compared with those of Smet and Van Assche. Finally, the Gauss hypergeometric discrete orthogonal polynomials, also known as generalized Hahn of type I, are also studied. Laguerre-Freud equations are found, and the differences with the equations found by Dominici and by Filipuk and Van Assche are provided., Comment: Some minor typos corrected In the new version the bibliography has been updated and the text improved. Typo in Theorem 8 corrected. arXiv admin note: substantial text overlap with arXiv:2107.01747

The Cholesky factorization of the moment matrix is applied to discrete orthogonal polynomials on the homogeneous lattice. In particular, semiclassical discrete orthogonal polynomials, which are built in terms of a discrete Pearson equation, are studied. The Laguerre-Freud structure semi-infinite matrix that models the shifts by $\pm 1$ in the independent variable of the set of orthogonal polynomials is introduced. In the semiclassical case it is proven that this Laguerre-Freud matrix is banded. From the well known fact that moments of the semiclassical weights are logarithmic derivatives of generalized hypergeometric functions, it is shown how the contiguous relations for these hypergeometric functions translate as symmetries for the corresponding moment matrix. It is found that the 3D Nijhoff-Capel discrete Toda lattice describes the corresponding contiguous shifts for the squared norms of the orthogonal polynomials. The continuous Toda for these semiclassical discrete orthogonal polynomials is discussed and the compatibility equations are derived. It also shown that the Kadomtesev-Petvishvilii equation is connected to an adequate deformed semiclassical discrete weight, but in this case the deformation do not satisfy a Pearson equation.

The recently found hypergeometric multiple orthogonal polynomials on the step-line by Lima and Loureiro are shown to be random walk polynomials. It is proven that the corresponding Jacobi matrix and its transpose, which are nonnegative matrices and describe higher recurrence relations, can be normalized to two stochastic matrices, dual to each other. Using the Christoffel-Darboux formula on the step-line and the Poincar\'e theory for non-homogeneous recurrence relations it is proven that both stochastic matrices are related by transposition in the large $n$ limit. These random walks are beyond birth and death, as they describe a chain in where transitions to the two previous states are allowed, or in the dual to the two next states.The corresponding Karlin-McGregor representation formula is given for these new Markov chains. The regions of hypergeometric parameters where the Markov chains are recurrent or transient are given. Stochastic factorizations, in terms of pure birth and of pure death factors, for the corresponding Markov matrices of types I and II, are provided.Twelve uniform Jacobi matrices and the corresponding random walks, related to a Jacobi matrix of Toeplitz type, and theirs stochastic or semi-stochastic matrices (with sinks and sources), that describe Markov chains beyond birth and death, are found and studied. One of these uniform stochastic cases, which is a recurrent random walk, is the only hypergeometric multiple random walk having a uniform stochastic factorization. The corresponding weights, Jacobi and Markov transition matrices and sequences of type II multiple orthogonal polynomials are provided. Chain of Christoffel transformations connecting the stochastic uniform tuples between them, and the semi-stochastic uniform tuples, between them, are presented., Comment: In the revision we added explicit expressions for the type I multiple orthogonal polynomials associated with the uniform tuples

Multiple orthogonal polynomials with respect to two weights on the step-line are considered. A connection between different dual spectral matrices, one banded (recursion matrix) and one Hessenberg, respectively, and the Gauss-Borel factorization of the moment matrix is given. It is shown a hidden freedom exhibited by the spectral system related to the multiple orthogonal polynomials. Pearson equations are discussed, a Laguerre-Freud matrix is considered, and differential equations for type I and II multiple orthogonal polynomials, as well as for the corresponding linear forms are given. The Jacobi-Pi\~neiro multiple orthogonal polynomials of type I and type II are used as an illustrating case and the corresponding differential relations are presented. A permuting Christoffel transformation is discussed, finding the connection between the different families of multiple orthogonal polynomials. The Jacobi-Pi\~neiro case provides a convenient illustration of these symmetries, giving linear relations between different polynomials with shifted and permuted parameters. We also present the general theory for the perturbation of each weight by a different polynomial or rational function aka called Christoffel and Geronimus transformations. The connections formulas between the type II multiple orthogonal polynomials, the type I linear forms, as well as the vector Stieltjes-Markov vector functions is also presented. We illustrate these findings by analyzing the special case of modification by an even polynomial., Comment: Revised version, completely new section on general Christoffel and Geronimus for multiple orthogonal polynomials on the stepline

Given a non-negative Jacobi matrix describing higher order recurrence relations for multiple orthogonal polynomials of type~II and corresponding linear forms of type I, a general strategy for constructing a pair of stochastic matrices, dual to each other, is provided. The corresponding Markov chains (or 1D random walks) allow, in one transition, to reach for the N-th previous states, to remain in the state or reach for the immediately next state. The dual Markov chains allow, in one transition, to reach for the N-th next states, to remain in the state or reach for immediately previous state. The connection between both dual Markov chains is discussed at the light of the Poincar\'e's theorem on ratio asymptotics for homogeneous linear recurrence relations and the Christoffel-Darboux formula within the sequence of multiple orthogonal polynomials and linear forms of type I. The Karlin-McGregor representation formula is extended to both dual random walks, and applied to the discussion of the corresponding generating functions and first-passage distributions. Recurrent or transient character of the Markov chain is discussed. Steady state and some conjectures on its existence and the relation with mass points are also given. The Jacobi-Pi\~neiro multiple orthogonal polynomials are taken as a case study of the described results. For the first time in the literature, an explicit formula for the type~I Jacobi--Pi\~neiro polynomials is determined. The region of parameters where the Markov chains are recurrent or transient is given, and it is conjectured that when recurrent, the Markov chains are null recurrent and, consequently, the expected return times are infinity. Examples of recurrent and transient Jacobi--Pi\~neiro random walks are constructed explicitly.

The Gauss-Borel or $LU$ factorization of Gram matrices of bilinear forms is the pivotal element in the discussion of the theory of biorthogonal polynomials. The construction of biorthogonal families of polynomials and its second kind functions, of the spectral matrices modeling the multiplication by the independent variable $x$, the Christoffel-Darboux kernel and its projection properties, are discussed from this point of view. Then, the Hankel case is presented and different properties, specific of this case, as the three terms relations, Heine formulas, Gauss quadrature and the Christoffel-Darboux formula are given. The classical orthogonal polynomial of Hermite, Laguerre and Jacobi type are discussed and characterized within this scheme. Finally, it is shown who this approach is instrumental in the derivation of Christoffel formulas for general Christoffel and Geronimus perturbations of the bilinear forms., Comment: 23 pages, 2 figures

In this paper the Riemann-Hilbert problem, with jump supported on a appropriate curve on the complex plane with a finite endpoint at the origin, is used for the study of corresponding matrix biorthogonal polynomials associated with Laguerre type matrices of weights ---which are constructed in terms of a given matrix Pearson equation. First and second order differential systems for the fundamental matrix, solution of the mentioned Riemann-Hilbert problem are derived. An explicit and general example is presented to illustrate the theoretical results of the work. Related matrix eigenvalue problems for second order matrix differential operators and non-Abelian extensions of a family of discrete Painlev\'e IV equations are discussed., Comment: arXiv admin note: text overlap with arXiv:1807.07119

In this paper we use the Riemann-Hilbert problem, with jumps supported on appropriate curves in the complex plane, for matrix biorthogonal polynomials and apply it to find Sylvester systems of differential equations for the orthogonal polynomials and its second kind functions as well. For this aim, Sylvester type differential Pearson equations for the matrix of weights are shown to be instrumental. Several applications are given, in order of increasing complexity. First, a general discussion of non-Abelian Hermite biorthogonal polynomials in the real line, understood as those whose matrix of weights is a solution of a Sylvester type Pearson equation with coefficients first order matrix polynomials, is given. All these is applied to the discussion of possible scenarios leading to eigenvalue problems for second order linear differential operators with matrix eigenvalues. Nonlinear matrix difference equations are discussed next. Firstly, for the general Hermite situation a general non linear relation (non trivial because the non commutativity features of the setting) for the recursion coefficients is gotten. In the next case of higher difficulty, degree two polynomials are allowed in the Pearson equation, but the discussion is simplified by considering only a left Pearson equation. In the case, the support of the measure is on an appropriate branch of an hyperbola. The recursion coefficients are shown to fulfill a non-Abelian extension of the alternate discrete Painlev\'e I equation. Finally, a discussion is given for the case of degree three polynomials as coefficients in the left Pearson equation characterizing the matrix of weights. However, for simplicity only odd polynomials are allowed. In this case, a new and more general matrix extension of the discrete Painlev\'e I equation is found., Comment: 38 pages, 1 figure

We introduce a large class of Sobolev bi-orthogonal polynomial sequences arising from a $LU$-factorizable moment matrix and associated with a suitable measure matrix that characterizes the Sobolev bilinear form. A theory of deformations of Sobolev bilinear forms is also proposed. We consider both polynomial deformations and a class of transformations related to the action of linear operators on the entries of a given bilinear form. Transformation formulae among new and old polynomial sequences are determined. Finally, integrable hierarchies of evolution equations arising from the factorization of a time deformation of the moment matrix are presented., Comment: 52 pages

This paper is a continuation of the recent paper "CMV biorthogonal Laurent polynomials: Christoffel formulas for Christoffel and Geronimus transformations" by the same authors. The behavior of quasidefinite sesquilinear forms for Laurent polynomials in the complex plane, characterized by bivariate linear functionals, and corresponding CMV biorthogonal Laurent polynomial families --including Sobolev and discrete Sobolev orthogonalities--- under two type of Geronimus--Uvarov transformations is studied. Either the linear functionals are multiplied by a Laurent polynomial and divided by the complex conjugate of a Laurent polynomial, with the addition of appropriate masses (linear functionals supported on the zeros of the perturbing Laurent polynomial in the denominator) or vice-versa, multiplied by the complex conjugate of a Laurent polynomial and divided by a Laurent polynomial. The connection formulas for the CMV biorthogonal Laurent polynomials, their norms, and Christoffel--Darboux kernels are given. For prepared Laurent polynomials these connection formulas lead to quasideterminantal (quotient of determinants) Christoffel formulas expressing an arbitrary degree perturbed biorthogonal Laurent polynomial in terms of $2N_{\mathfrak C}+2N_\Gamma$ unperturbed biorthogonal Laurent polynomials, their second kind functions or Christoffel--Darboux kernel and its mixed versions. When the linear functionals are supported on the unit circle, a particularly relevant role is played by the reciprocal polynomial, and the Christoffel formulas provide now with two possible ways of expressing the same perturbed quantities in terms of the original ones, one using only the nonperturbed biorthogonal family of Laurent polynomials, and the other using the Christoffel--Darboux kernel and its mixed versions.

Quasidefinite sesquilinear forms for Laurent polynomials in the complex plane and corresponding CMV biorthogonal Laurent polynomial families are studied. Bivariate linear functionals encompass large families of orthogonalities like Sobolev and discrete Sobolev types. Two possible Christoffel transformations of these linear functionals are discussed. Either the linear functionals are multiplied by a Laurent polynomial, or are multiplied by the complex conjugate of a Laurent polynomial. For the Geronimus transformation, the linear functional is perturbed in two possible manners as well, by a division by a Laurent polynomial or by a complex conjugate of a Laurent polynomial, in both cases the addition of appropriate masses (linear functionals supported on the zeros of the perturbing Laurent polynomial) is considered. The connection formulas for the CMV biorthogonal Laurent polynomials, its norms, and Christoffel-Darboux kernels, in all the four cases, are given. For the Geronimus transformation, the connection formulas for the second kind functions and mixed Christoffel-Darboux kernels are also given in the two possible cases. For prepared Laurent polynomials, i.e. of the form $L(z)=L_nz^n+\cdots+L_{-n}z^{-n}$, $L_nL_{-n}\neq 0$, these connection formulas lead to quasideterminantal (quotient of determinants) Christoffel formulas for all the four transformations, expressing an arbitrary degree perturbed biorthogonal Laurent polynomial in terms of $2n$ unperturbed biorthogonal Laurent polynomials, their second kind functions or Christoffel-Darboux kernels and its mixed versions. The unit circle case, given its exceptional properties, is discussed in more detail.