Van Vleck spectra of high-order Heun operators:\ finite-band universality and exterior asymptotics

Abstract: We study high-order analogues of the classical Heun operator of Fuchs index one, [ \dq=\sum_{i=1}^k Q_i(z)\frac{d^i}{dz^i}, \qquad °Q_i\le i+1, \qquad °Q_k=k+1. ] For a fixed degree $n$ we consider the linear Van Vleck polynomials $V$ for which $\dq+V$ has a polynomial solution of degree $n$, and we form the spectral polynomial $Sp_n$ whose zeros are the zeros of these Van Vleck polynomials. The main result is a finite-band determinant representation and the resulting universality theorem: after normalization, all fixed power sums of the zeros of $Sp_n$ have limits given by explicit constant-term formulae depending only on the leading coefficient $Q_k$. The lower coefficients of $\dq$ enter only lower order correction terms. Combining this with the localization theorem for Van Vleck roots, we strengthen the usual germ-at-infinity conclusion to locally uniform convergence of the normalized Cauchy transforms and logarithmic potentials on the whole exterior of the convex hull of the zeros of $Q_k$. We also prove a determinacy criterion: if the spectral roots are asymptotically confined to a compact set with empty interior and connected complement, then the finite-band moments determine the actual weak limit. In particular, when the zeros of $Q_k$ are collinear the root-counting measures of $Sp_n$ converge weakly to a probability measure supported on the corresponding segment; this limit is independent of all lower coefficients of $\dq$. Finally, we prove holonomicity of the exterior Cauchy transform and derive Picard--Fuchs equations for the WKB periods, with an explicit third-order equation in the first non-classical case $k=3$. The paper ends with a precise mother-body conjecture for the genuinely complex case, clearly separated from the unconditional results.

• The paper establishes a universality theorem showing that the asymptotic distribution of Van Vleck roots in high-order Heun operators depends solely on the leading polynomial coefficient.
• It introduces a finite-band matrix formulation and moment analysis to rigorously describe the exterior Cauchy transform and logarithmic potential of spectral measures.
• A novel WKB approach and the conjectured spectral tree framework provide fresh insights into the mother-body phenomenon and interior spectral structure in complex configurations.

Van Vleck Spectra for High-Order Heun Operators: Universality, Potentials, and Asymptotics

Introduction

This work addresses the spectral geometry of polynomial Van Vleck parameters for high-order analogues of the classical Heun operator—linear differential operators of Fuchs index one. The central object of inquiry is the spectral polynomial $Sp_n$, whose roots (Van Vleck parameters) parameterize polynomial eigenfunctions (Stieltjes polynomials) for a class of operators generalizing the Lamé equation. The analysis establishes a comprehensive universality theorem for the asymptotic distribution of these roots and elucidates their exterior potential-theoretic structure, with particular attention to the distinction between unconditional, averaged results and conjectural, tree-like fine structure in the generic complex case.

Main Theoretical Framework

Given a high-order differential operator

$$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$

the focus is on those at most linear polynomials $V$ for which $\mathcal{D}_Q + V$ admits a polynomial solution of fixed degree $n$. The associated Van Vleck spectral polynomial is then

$Sp_n(z) = \prod_{j=1}^{n+1} (z - z_{n,j}),$

with $\{z_{n,j}\}$ the spectral parameters of interest. These roots are analyzed via root-counting measures $\mu_n$.

A deterministic, finite-band matrix formulation arises for $Sp_n$ in terms of truncations of the action of $\mathcal{D}_Q-\lambda_n(z-t)$ on degree-$$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$0 polynomials, yielding substantial technical leverage, especially for asymptotic limits $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$1.

Finite-Band Universality and Exterior Potentials

A central result is the finite-band moment theorem, which characterizes the normalized power sums of spectral roots: $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$2 where $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$3 is the constant term in $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$4 and $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$5 are coefficients of the leading polynomial $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$6.

Crucially, all power sum moments in the limit depend solely on $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$7, and the influence of non-leading coefficients enters only at lower correction orders, establishing strong spectral universality. This moment structure determines the exterior Cauchy transform and logarithmic potential of the asymptotic measures.

Figure 1: The roots of the spectral polynomial $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$8 for the fourth-order operator with a generic complex quintic $$\mathcal{D}_Q = \sum_{i=1}^k Q_i(z) \frac{d^i}{dz^i}, \quad \deg Q_i \le i+1, \quad \deg Q_k = k+1,$$9.

The finiteness-band determinant representation is exploited to derive locally uniform convergence of the normalized Cauchy transforms $V$0 and logarithmic potentials $V$1 to analytic and harmonic functions, respectively, throughout the exterior of the convex hull $V$2 of zeros of $V$3. Thus, the exterior field is independent of lower-order perturbations.

Figure 2: The roots of the spectral polynomial $V$4 for a second-order operator with a cubic $V$5. The Van Vleck roots cluster inside the convex hull of the cubic's roots.

Moment Determinacy and Rigidity Results

The determinacy result is sharp: if the roots of $V$6 asymptotically localize on a compact set $V$7 of empty interior and connected complement, then the normalized moments uniquely specify the weak limit of $V$8. In particular, if the roots of $V$9 are collinear (including the real-rooted case), the limit is a probability measure supported on the segment spanned by those roots, and is again dictated exclusively by $\mathcal{D}_Q + V$0.

Asymptotics, Holonomicity, and Picard-Fuchs Equations

The paper establishes that the limiting exterior Cauchy transform $\mathcal{D}_Q + V$1 is holonomic in $\mathcal{D}_Q + V$2, via creative telescoping on the algebraic integrand arising from the finite-band description. Explicit Picard-Fuchs equations are derived for the so-called WKB periods in the spectral parameter $\mathcal{D}_Q + V$3. For $\mathcal{D}_Q + V$4 (the first non-classical case), the differential equation is: $\mathcal{D}_Q + V$5 admitting explicit comparison with finite-band Cauchy transforms and conjectural spectral trees.

Figure 3: Numerical approximations of four auxiliary averaged measures for $\mathcal{D}_Q + V$6 (colored sets), alongside high-precision roots of $\mathcal{D}_Q + V$7 and the zeros of $\mathcal{D}_Q + V$8.

Mother-Body Phenomenon and the Spectral Tree Conjecture

While the averaged potential and exterior moments are robustly determined, the internal structure of the limiting measure is more subtle in the general complex case. For complex $\mathcal{D}_Q + V$9 with simple roots, numerical evidence demonstrates that the Van Vleck roots cluster on a thin, tree-like support connecting the roots of $n$0, rather than being supported on a two-dimensional region.

The Spectral Tree/Mother-Body Conjecture posits that the limiting measure $n$1 is the positive mother body of the exterior field generated by the finite-band average—i.e., it is a probability measure on a finite planar tree, uniquely determined by its potential outside $n$2 and having support inclusive of all roots of $n$3. This is a sharp and technical refinement of the universality result: the support is in general much smaller than the two-dimensional support of the averaged finite-band measure.

Figure-Based Illustration of the Asymptotics and Mother-Body Tree

Numerical calculations show (Figures 1–3 above) that, for high-degree $n$4, the Van Vleck roots aggregate on loci forming a tree inside $n$5, whose leaves are precisely the roots of $n$6. These figures contrast the spectral root set with the supports of the frozen averaged measures, highlighting the skeletonization (that is, one-dimensional reduction) induced by the genuine spectral measure.

Figure 4: The union of roots of 861 quadratic Van Vleck polynomials for classical Lamé equation, visualizing the rich distributional structure in the non-determinantal case.

WKB Approach and Structural Conjectures

A WKB (semiclassical) analysis is advanced for the geometric structure of the spectral tree. The conjecture is that the support of the spectral measure is characterized by loci where the real parts of certain WKB periods vanish, subject to a positivity (S-property) constraint and quantization condition, leading to the selection of a minimal carrier tree. The approach blends asymptotic analysis, period calculations, and algebraic geometry of differentials over auxiliary Riemann surfaces.

Implications and Directions for Future Work

This framework solidifies a new universality principle for polynomial spectral problems associated with high-order Heun operators: the asymptotic exterior behavior of spectral roots is determined solely by the leading term, up to all finite moments. The nuanced issue of the “interior” (support structure and uniqueness of the mother body) for complex configurations is partially conjectural, with strong evidence for selection of tree-like supports. These results impact both the spectral theory of non-self-adjoint operators and logarithmic potential theory in the complex plane.

The research leaves open rigorous confirmation of uniqueness and positivity for the mother body in the most general settings, as well as a deeper connection between WKB quantization/periods and finite-band moment formulas, particularly in the presence of multiple active cycles and nontrivial support topologies. These questions are compelling from both mathematical and mathematical-physics perspectives.

Conclusion

This paper establishes a detailed, technically robust theory of the asymptotic root distribution for the Van Vleck spectra associated with high-order Heun operators. It demonstrates a strong universality in moments and exterior potential, provides explicit algebraic tools for analysis, and articulates a precise conjectural framework (supported by numerics) for the internal (support) structure of the limiting measures. The interplay of potential theory, spectral asymptotics, and algebraic geometry highlighted here has significant theoretical leverage and sets the stage for further analytic and computational advances in the spectral theory of differential operators (2607.02700).