Non-commutative Multiple Orthogonal Polynomials

• Non-commutative multiple orthogonal polynomials are polynomial systems with coefficients in a free division ring that satisfy multiple orthogonality conditions via quasideterminantal constructions.
• They employ Gelfand–Retakh quasideterminants to replace classical determinants, yielding explicit formulas for polynomials and tau-function normalizations linked to integrable Toda dynamics.
• Generalized nearest-neighbor recurrences, duality relations, and Hirota-type equations interconnect the polynomial structure with multidimensional discrete-time Toda systems.

Non-commutative multiple orthogonal polynomials are polynomial systems in a central indeterminate $x$ with coefficients in a non-commutative algebraic setting—specifically, in the recent formal theory, a free division ring generated by moments or bi-moments—subject to multiple orthogonality conditions indexed by a multi-index $s=(s_1,\dots,s_r)$. In this setting, determinant formulas are replaced by quasideterminants, left and right module structures become essential, and the normalization factors of the monic polynomials play the role of non-commutative tau-functions satisfying Hirota-type relations. The subject lies at the intersection of multiple orthogonality, matrix orthogonal polynomials, bi-orthogonality, and integrable systems, with a formal quasideterminantal treatment developed in "Non-commutative multiple bi-orthogonal polynomials: formal approach and integrability" (Doliwa, 2 Oct 2025), and with a broader Padé–Toda context surveyed in "Hermite-Padé approximation, multiple orthogonal polynomials, and multidimensional Toda equations" (Doliwa, 2023).

1. Algebraic setting and definition

The formal theory takes place over a free division ring $\mathcal{R}$ generated by moments or bi-moments, with polynomials in a commutative indeterminate $x$, so the ambient polynomial ring is $\mathcal{R}[x]$ (Doliwa, 2 Oct 2025). For the specialization to non-commutative multiple orthogonal polynomials (NCMOP), one fixes $r$ sequences of moments $(\nu_j^{(k)})_{j\in\mathbb{N}_0}$, $k=1,\dots,r$, together with an involutive anti-automorphism $(\cdot)^*$ satisfying $(\nu_j^{(k)})^*=\nu_j^{(k)}$, extended to $s=(s_1,\dots,s_r)$0 coefficientwise (Doliwa, 2 Oct 2025).

The basic orthogonality data are sesquilinear forms

$s=(s_1,\dots,s_r)$1

Given a multi-index $s=(s_1,\dots,s_r)$2, $s=(s_1,\dots,s_r)$3, the monic NCMOP $s=(s_1,\dots,s_r)$4 are defined as polynomials of degree $s=(s_1,\dots,s_r)$5 satisfying

$s=(s_1,\dots,s_r)$6

for each $s=(s_1,\dots,s_r)$7, together with normalization

$s=(s_1,\dots,s_r)$8

These conditions are the non-commutative analogue of type II multiple orthogonality, but with ordering controlled by the anti-automorphism and by left/right actions (Doliwa, 2 Oct 2025).

This formal NCMOP theory appears as a specialization of a broader theory of non-commutative multiple bi-orthogonal polynomials. In that larger setting one starts from $s=(s_1,\dots,s_r)$9 bi-moment arrays $\mathcal{R}$0 and bilinear forms

$\mathcal{R}$1

The reduction to NCMOP is the Hankel condition

$\mathcal{R}$2

which converts the general bi-moment problem into a moment problem of multiple-orthogonality type (Doliwa, 2 Oct 2025).

A common source of confusion is the phrase “non-commutative.” In the NCMOP setting it refers to the coefficients, moments, and normalization data living in a non-commutative division ring. This differs from constructions in which ordinary scalar polynomials are generated by automorphisms of non-commutative operator algebras; that second meaning becomes relevant later.

2. Quasideterminantal construction

The formal machinery is built from Gelfand–Retakh quasideterminants. For a square matrix $\mathcal{R}$3, the $\mathcal{R}$4-quasideterminant is defined by

$\mathcal{R}$5

whenever the relevant inverse exists, where $\mathcal{R}$6 is obtained by deleting row $\mathcal{R}$7 and column $\mathcal{R}$8, and $\mathcal{R}$9, $x$0 are the truncated row and column. In the commutative case,

$x$1

so quasideterminants recover the usual determinantal quotient only after commutativity is imposed (Doliwa, 2 Oct 2025).

Under the Hankel reduction, the monic NCMOP admit an explicit block Hankel quasideterminantal formula: $x$2 The normalization factors are given by the analogous quasideterminants

$x$3

These $x$4 are the normalization functions and also the tau-functions of the integrable structure (Doliwa, 2 Oct 2025).

The simplest nontrivial example occurs for $x$5 and $x$6, when

$x$7

with $x$8 the upper-left $x$9 block. The corresponding normalizations are

$\mathcal{R}[x]$0

provided the indicated inverses exist (Doliwa, 2 Oct 2025).

Existence is therefore inseparable from non-commutative invertibility. The required quasideterminants exist only when the relevant Schur complements are invertible, and the paper formulates this in terms of normal multi-indices and perfect systems: $\mathcal{R}[x]$1 is normal if the linear system for the non-leading coefficients of $\mathcal{R}[x]$2 has a unique solution, and the system is perfect if all multi-indices are normal (Doliwa, 2 Oct 2025).

3. Orthogonality, duality, and recurrence structure

The broader non-commutative multiple bi-orthogonal theory contains two monic polynomial families, $\mathcal{R}[x]$3 and $\mathcal{R}[x]$4, with left and right orthogonality conditions. For the $\mathcal{R}[x]$5-family,

$\mathcal{R}[x]$6

while for the $\mathcal{R}[x]$7-family,

$\mathcal{R}[x]$8

Their normalizations are $\mathcal{R}[x]$9 and $r$0. Under the Hankel reduction to NCMOP, these dual families are related by

$r$1

so the multiple-orthogonality problem inherits a canonical transposed companion (Doliwa, 2 Oct 2025).

The NCMOP satisfy a generalized nearest-neighbor recurrence: $r$2 This is the non-commutative multiple analogue of the three-term relation, but it is indexed by the full multi-index lattice and contains ordered right coefficients. The coefficients are expressed by the normalization data: $r$3 and

$r$4

where $r$5 is the corresponding moment-shifted quasideterminant (Doliwa, 2 Oct 2025).

These recurrence coefficients satisfy compatibility relations that replace the commutative zero-curvature identities by ordered identities: $r$6

$r$7

and

$r$8

The formal significance of these formulas is that orthogonality, normalization, and recurrence are not independent layers: each is encoded by the same quasideterminantal data (Doliwa, 2 Oct 2025).

4. Hirota equations, Lax systems, and multidimensional Toda dynamics

The normalization factors $r$9 function as non-commutative tau-functions. For $(\nu_j^{(k)})_{j\in\mathbb{N}_0}$0, they satisfy a potential form of the non-commutative Hirota system: $(\nu_j^{(k)})_{j\in\mathbb{N}_0}$1 and, for distinct $(\nu_j^{(k)})_{j\in\mathbb{N}_0}$2,

$(\nu_j^{(k)})_{j\in\mathbb{N}_0}$3

The polynomials themselves solve the corresponding linear problem

$(\nu_j^{(k)})_{j\in\mathbb{N}_0}$4

so the orthogonal polynomial system appears as the wave function of the non-commutative Hirota structure (Doliwa, 2 Oct 2025).

A discrete-time variable is introduced by shifting moments,

$(\nu_j^{(k)})_{j\in\mathbb{N}_0}$5

which leads to time-dependent block Hankel matrices and a Lax-type pair: $(\nu_j^{(k)})_{j\in\mathbb{N}_0}$6 and

$(\nu_j^{(k)})_{j\in\mathbb{N}_0}$7

Compatibility yields a non-commutative multidimensional discrete-time Toda system: $(\nu_j^{(k)})_{j\in\mathbb{N}_0}$8

$(\nu_j^{(k)})_{j\in\mathbb{N}_0}$9

and

$k=1,\dots,r$0

The same discrete-time formalism also implies mixed space-time Hirota-type relations for $k=1,\dots,r$1 (Doliwa, 2 Oct 2025).

The broader significance of these identities is clarified by the 2023 survey, which reviews the connection between Hermite–Padé approximation, multiple orthogonal polynomials, and multidimensional Toda equations, and explicitly points to a non-commutative extension in which determinants are replaced by quasideterminants and the resulting Paszkowski-type constraints become non-commutative Hirota–Miwa equations (Doliwa, 2023). In that sense, NCMOP occupy the non-commutative counterpart of a well-established Padé–Toda correspondence.

5. Relation to multiple orthogonality, matrix orthogonality, and bi-orthogonality

The formal theory was introduced not merely as a non-commutative variant of ordinary MOPs, but as a simultaneous generalization of multiple orthogonality, matrix orthogonal polynomials, and bi-orthogonality (Doliwa, 2 Oct 2025). This point matters structurally.

First, the multi-index $k=1,\dots,r$2 and the $k=1,\dots,r$3 distinct forms $k=1,\dots,r$4 retain the multiple-orthogonality aspect. Second, the non-commutative coefficient ring accommodates matrix- and operator-valued specializations. Third, the presence of two monic families $k=1,\dots,r$5 and $k=1,\dots,r$6, with left and right orthogonality, places the construction in a genuinely bi-orthogonal framework. The NCMOP reduction is therefore not the replacement of scalar moments by matrices in an otherwise unchanged theory; it is a specialization of a larger left/right formalism.

The 2023 survey formulates this broader perspective in terms of Hermite–Padé approximation and block moment structures. In the commutative case, type II MOPs are characterized by

$k=1,\dots,r$7

and their moment determinants produce multidimensional Toda equations (Doliwa, 2023). The non-commutative extension modifies every component of this picture: moments become operator-valued, determinants become quasideterminants, recurrence coefficients become ordered quantities, and compatibility conditions acquire non-Abelian form.

A plausible implication is that the phrase “non-commutative multiple orthogonal polynomials” should be reserved for settings in which the orthogonality data themselves are non-commutative, not merely for matrix rephrasings of scalar multiple orthogonality. The formal theory of (Doliwa, 2 Oct 2025) makes that distinction explicit through left/right bilinear forms, anti-automorphisms, and quasideterminantal tau-functions.

6. Operator-algebraic constructions, bispectrality, and a second meaning of non-commutativity

An earlier and complementary line of work constructs vector orthogonal polynomials and certain type II multiple orthogonal polynomials from automorphisms of non-commutative algebras, especially the Weyl algebra and its difference analogues (Horozov, 2016). There the central objects are not non-commuting moments but non-commuting operators. One fixes

$k=1,\dots,r$8

or, in the discrete case,

$k=1,\dots,r$9

defines an automorphism

$(\cdot)^*$0

and constructs polynomials by dressing the basic bispectral wave function: $(\cdot)^*$1 The transformed operator

$(\cdot)^*$2

acts diagonally,

$(\cdot)^*$3

while the transformed anti-isomorphism produces a fixed-length finite-term recurrence in $(\cdot)^*$4 (Horozov, 2016).

In that framework, vector orthogonal polynomials are $(\cdot)^*$5-orthogonal in the sense of van Iseghem and Maroni, equivalently characterized by a $(\cdot)^*$6-term recurrence, and they can be regarded—after regrouping orthogonality conditions—as certain type II multiple orthogonal polynomials. The paper also shows that continuous and discrete families are related by a Mellin-type transform $(\cdot)^*$7, which maps $(\cdot)^*$8 to $(\cdot)^*$9 and intertwines the continuous and discrete constructions (Horozov, 2016).

This is not the same theory as NCMOP over a division ring. The distinction is substantial. In (Horozov, 2016), “non-commutative” refers to the algebraic machinery used to generate scalar polynomial families with Bochner’s property and bispectrality. In (Doliwa, 2 Oct 2025), it refers to the coefficient and moment algebra of the polynomials themselves. The two perspectives are nevertheless adjacent: both replace positivity-based orthogonality by algebraic functionals, both privilege recurrence and spectral equations, and both connect orthogonal polynomial systems to integrable structures.

7. Scope, formal character, and open directions

The current formal theory is explicitly quasideterminantal and integrable-systems oriented. It assumes a free division ring, a central indeterminate, and the existence of the necessary inverses for the quasideterminants. Normality of multi-indices and perfect systems are therefore structural assumptions rather than secondary regularity conditions (Doliwa, 2 Oct 2025).

A direct non-commutative Riemann–Hilbert characterization is not developed in the formal NCMOP work, and the 2023 survey likewise treats operator-valued jump relations and ordering issues as beyond its scope (Doliwa, 2023). The emphasis is instead on block moment constructions, linear systems, Hirota identities, and Toda-type compatibility. This suggests that the current theory should be read as a formal algebraic foundation rather than as a measure-theoretic or analytic classification.

Related work indicates several adjacent directions. The 2023 survey points to non-commutative Hermite–Padé approximation via quasideterminants and to non-commutative multidimensional discrete Toda equations (Doliwa, 2023). The operator-algebraic VOP framework identifies extensions via Darboux transformations, commuting algebras for Krall-type families, Toda-type integrable-system connections, and bispectral algebras beyond $(\nu_j^{(k)})^*=\nu_j^{(k)}$0 (Horozov, 2016). Taken together, these developments suggest that non-commutative multiple orthogonality is becoming a meeting point for quasideterminant calculus, multicomponent recurrences, Hermite–Padé theory, and non-Abelian integrable hierarchies.

Within that landscape, the defining contribution of the recent formal theory is precise: it provides explicit Hankel-type quasideterminantal formulas for monic non-commutative multiple orthogonal polynomials, identifies their normalization functions as non-commutative tau-functions, derives their generalized nearest-neighbor recurrences, and embeds the whole construction into a non-commutative multidimensional discrete-time Toda system (Doliwa, 2 Oct 2025).