Szegő-Type Recurrence Relations

• Szegő-type recurrence relations are algebraic identities that define orthogonal polynomials on the unit circle and extend to multivariate and multiple orthogonality settings.
• They incorporate key features such as Verblunsky coefficients, compatibility conditions, and determinantal representations, which are crucial for spectral analysis.
• Applications span spectral theory, moment problems, and approximation theory, providing explicit computational methods and structural insights.

Szegő-type recurrence relations are algebraic identities governing sequences of orthogonal polynomials on the unit circle and their generalizations to multiple orthogonality and higher dimensions. These recurrences encode spectral information about underlying measures or moment functionals and are fundamental in analysis, approximation theory, and spectral theory. The Szegő recurrence’s key features—Verblunsky coefficients, direct and inverse recurrences, compatibility relations, and determinantal representations—extend naturally to multiple orthogonality and multivariate settings, facilitating explicit computations and structural insights.

1. Classical Szegő Recurrence and Verblunsky Coefficients

The Szegő recurrence for orthogonal polynomials on the unit circle (OPUC) is established as follows. Given a nontrivial probability measure $\sigma$ on the unit circle $T$, the sequence of monic OPUC $\{\Phi_n\}$ and their reversed polynomials $\Phi_n^*$ are determined by

$$\Phi_{n+1}(z) = z\,\Phi_n(z) - \overline{\alpha_n}\,\Phi_n^*(z), \qquad \Phi_{n+1}^*(z) = \Phi_n^*(z) - \alpha_n z \Phi_n(z),$$

where $\Phi_n^*(z) = z^n \overline{\Phi_n(1/\bar z)}$ and $\alpha_n$ are the Verblunsky coefficients, with $|\alpha_n| < 1$ (Castillo et al., 2015). These coefficients uniquely parameterize the measure, and their square summability is equivalent to Szegő’s theorem: $$\int_0^{2\pi} \log w(e^{i\theta})\,\frac{d\theta}{2\pi} = \sum_{j=1}^\infty \log(1 - |\alpha_j|^2),$$ where $w$ is the density of the absolutely continuous part of $\sigma$ (Gibson, 2022).

2. Multiple Orthogonality and Szegő-type Recurrences

Szegő-type recurrence relations have been generalized to multiple orthogonal polynomials on the unit circle (MOPUC), involving simultaneous orthogonality with respect to $r$ probability measures $\{\mu_j\}_{j=1}^r$ (Vaktnäs et al., 2024, Kozhan et al., 8 Jan 2026). For multi-indices $\vec n = (n_1, ..., n_r)\in\mathbb{Z}_+^r$, the type II MOPUC $\Phi_{\vec n}(z)$ and its reversed companion $\Phi^*_{\vec n}(z)$ satisfy:

Direct recurrence: $$\Phi_{\vec n+\mathbf e_k}(z) = z\,\Phi_{\vec n}(z) - \overline{\alpha_{\vec n,k}}\,\Phi^*_{\vec n}(z),\qquad k=1,\dots,r$$ where each $\alpha_{\vec n,k}$ is a generalized Verblunsky coefficient obeying $|\alpha_{\vec n,k}| < 1$. Off-diagonal coefficients $\rho_{\vec n,j}$ emerge via expansions in lower-degree bases (Vaktnäs et al., 2024).

Inverse recurrence: $$\Phi^*_{\vec n}(z) = \Phi^*_{\vec n-\mathbf e_k}(z) + \beta_{\vec n} z\,\Phi_{\vec n-\mathbf e_k}(z),\qquad n_k > 0$$ with $\beta_{\vec n} = \overline{\alpha_{\vec n,k}}$ for the single-measure case (Vaktnäs et al., 2024).

Multi-parameter families of recurrences connect shifts in $\bm{n}$ and $\bm{m}$ in the Laurent polynomial setting, and analogous relations hold for type I and type I$^*$ families (Kozhan et al., 8 Jan 2026). Compatibility (“zero curvature”) conditions among recurrence coefficients enforce algebraic consistency: $$\alpha_{\vec n} \beta_{\vec n} + \sum_{j=1}^r \rho_{\vec n,j} = 1$$ and more intricate partial-difference systems for multiple indices (Vaktnäs et al., 2024, Kozhan et al., 8 Jan 2026).

3. Multivariate and Higher-dimensional Szegő Machinery

The Szegő recurrence and associated Verblunsky-type theorems have multivariate extensions, established on $\mathbb{T}^d$ for $d\geq 1$ (Gibson, 2022). Schur functions $f:\mathbb{D}^d\to\mathbb{D}$ yield polynomials $\Phi_n$ and $\Phi_n^*$ via transfer matrices $M_n$, parametrized by Schur parameters $r_n$ and variable-allocation indices $v(n)$: $$\begin{pmatrix}\Phi_{n+1}(z)\\Phi_{n+1}^*(z)\end{pmatrix} = \frac{1}{\rho_{n+1}}\begin{pmatrix} z_{v(n+1)} & -\overline{r_{n+1}}\ - r_{n+1} z_{v(n+1)} & 1 \end{pmatrix}\begin{pmatrix}\Phi_n(z)\\Phi_n^*(z)\end{pmatrix}$$ with $\rho_{n+1} = \sqrt{1 - |r_{n+1}|^2}$ (Gibson, 2022). For $d=1$, this reduces to the classical Szegő case.

Szegő–Verblunsky-type theorems generalize, asserting that the convergence of the Szegő integral on $\mathbb{T}^d$ is equivalent to square summability of $\{r_j\}$, and analogous almost-periodic results hold for systems on irrational rotation flows.

4. Determinantal Representations and Christoffel–Darboux Formulas

Heine-type determinantal formulas provide explicit construction for MOPUC and their recurrence coefficients in terms of the moment matrix $T_{\bm{n};\bm{m}}$ (Kozhan et al., 8 Jan 2026): $\Phi_{\bm{n};\bm{m}}(z) = \frac{1}{\det T_{\bm{n};\bm{m}}}\det \begin{pmatrix} T_{\bm{n};\bm{m}} & \vdots & [z^{-|\bm{m}|},\dots,z^{|\bm{n}|}] \ \hdashline 0 & \dots & 0 \end{pmatrix}$ and similarly for $\Phi^*_{\bm{n};\bm{m}}(z)$. All recurrence coefficients admit closed-form determinants involving matrix shifts.

Christoffel–Darboux (CD) formulas extend to the multiple setting. For a path of normal indices $(\bm{n}_0;\bm{m}) \to \ldots \to (\bm{n}_N;\bm{m})$,

$$(\xi-z)\sum_{k=0}^{N-1}\Phi_{\bm{n}_k;\bm{m}}(z)\,\bm{\Xi}_{\bm{n}_{k+1};\bm{m}}(\xi) = \Phi^*_{\bm{n}_N;\bm{m}}(z)\,\bm{\Xi}^*_{\bm{n}_N;\bm{m}}(\xi) - z\xi \sum_{j=1}^r \rho_{\bm{n}_N;\bm{m},j} \Phi_{\bm{n}_N-\bm{e}_j;\bm{m}}(z)\,\bm{\Xi}_{\bm{n}_N+\bm{e}_j;\bm{m}}(\xi)$$

with a dual identity for the reversed families (Kozhan et al., 8 Jan 2026, Vaktnäs et al., 2024).

5. Szegő Mapping, Geronimus Relations, and Real-Line Analogues

The Szegő mapping explicitly relates OPUC and multiple orthogonal polynomials on the real line (MOPRL) via the Joukowski transformation: $$x = \frac{1}{2}(z + z^{-1}),\qquad d\sigma(\theta) = \tfrac{1}{2} d\mu(\cos\theta)$$ The mapping is formalized by

$M[x^k] = L[w^k + w^{-k}],\qquad L[z^k] = L[z^{-k}]$

(Castillo et al., 2015, Kozhan et al., 8 Jan 2026). Geronimus relations connect recurrence coefficients: $$d_{n+1} = \tfrac{1}{4}(1 - \alpha_{2n-1})(1 - \alpha_{2n}^2)(1 + \alpha_{2n+1}),\qquad b_{n+1} = \tfrac{1}{2}\bigl[\alpha_{2n}(1 - \alpha_{2n-1}) - \alpha_{2n-2}(1 + \alpha_{2n-1})\bigr]$$ Transformations and associated/anti-associated polynomial constructions follow via index shifts in these formulas. The measure, polynomial, and transform connections yield new recurrences, CD identities, and new expressions for associated perturbations of the OPUC/OPRL data.

6. Applications, Analogies, and Future Directions

Szegő-type recurrences unify the analytic and algebraic machinery for orthogonal polynomials on the unit circle, their multiple and Laurent generalizations, and higher-dimensional analogues. The structural analogy with nearest-neighbour recurrences for MOPRL extends to block-determinant identities, compatibility conditions, and explicit spectral parameters (Vaktnäs et al., 2024, Kozhan et al., 8 Jan 2026).

Applications range from spectral theory to moment problems, trace formulas for Schrödinger equations in impedance form, and Hermite–Padé approximation. The multivariate and almost-periodic Szegő–Verblunsky theorems facilitate analogues of classical results and the analysis of new function spaces (Gibson, 2022).

The ongoing development of Szegő-type recurrences in multiple, Laurent, and higher-dimensional settings continues to refine the interplay between orthogonal polynomial theory, measure transformations, and spectral analysis.