Wronskian Evolutions in Integrable Systems

• Wronskian Evolutions are defined as explicit determinantal representations that generate τ-functions, linking linear differential systems with nonlinear Hirota bilinear structures.
• They unify continuous, discrete, and supersymmetric models by providing conserved integrals of motion and explicit solution methods in integrable hierarchies.
• The approach extends to algebraic-geometric and combinatorial contexts, ensuring invariance through Plücker relations and facilitating explicit formulations for exceptional polynomials and soliton solutions.

A Wronskian evolution refers to the process by which functions, equations, or integrable systems evolve with time or discrete variables in such a manner that their underlying solutions, invariants, or τ-functions admit explicit determinantal representations in terms of Wronskians. These Wronskians encapsulate multi-component linear structures, provide explicit integrals of motion, and enforce bilinear integrability in a variety of settings: from classical ODE/PDE, through integrable soliton hierarchies, to algebraic and combinatorial constructions such as the KP hierarchy and exceptional orthogonal polynomials. The Wronskian evolution is underpinned by the interplay between linearity and nonlinear determinantal constraints, Plücker relations, and the algebraic geometry of flags and Grassmannians.

1. Wronskian Evolutions in Integrable Hierarchies

The Wronskian technique constructs rapidly verifiable explicit solutions (“τ-functions”) for soliton equations and their integrable discretizations. For instance, solutions to the KdV, modified KdV, and the “good” Boussinesq equation can be expressed as

$$\tau(x,t) = \operatorname{Wr}(\varphi_1, \dots, \varphi_N) = \det\left[\partial_x^{j-1}\varphi_i(x,t)\right]_{i,j=1}^N\,,$$

with the field variable $u$ supplied through a logarithmic derivative:

$u = 2\,\partial_x^2 \ln \tau(x,t).$

The entries $\varphi_i$ satisfy compatible linear differential or difference systems—a Linear Differential Equation System (LDES). Evolution in time (e.g., $\varphi_t = B\,\varphi$) naturally induces an evolution of $\tau$ consistent with the integrable hierarchy's Hirota bilinear identities, and explicit time evolution is constructed via exponentials or Jordan chains (Zhang, 2019, Dai et al., 2023, Zhang et al., 2012).

This structure persists for the modified KdV (with an auxiliary matrix implementing complex conjugation constraints in the entry equations) and for higher symmetries or reductions (e.g., lattice KdV, AKNS hierarchies) (Zhang et al., 2012, Li et al., 2014). Plücker relations for determinants guarantee that the $\tau$ function satisfies the required bilinear equations, unifying the analytic evolution (via Lax pairs) with the algebraic/determinantal structure (Zhang, 2019, Dai et al., 2023).

2. Wronskian Evolutions and Invariants in ODE Systems

For linear ODE and matrix systems, the evolution of the Wronskian determinant encodes key integrals of motion. Given a linear ODE system $d\Phi/dt = A(t)\Phi$, the Wronskian $W(t) = \det \Phi(t)$ evolves as

$\frac{dW}{dt} = W(t) \cdot \operatorname{Tr}A(t).$

If the trace vanishes (as for the moment hierarchy derived from the Jacobi equation along an Euler–Lagrange extremal), the Wronskian is conserved—providing a constant of motion along classical trajectories. This applies to the evolution of variational moments and underpins conservation laws in both ODE and quantum analogues ($$\langle\delta x^2\rangle\langle\delta p^2\rangle - \langle\delta x \delta p\rangle^2$$ constant) (Koshcheev, 15 Aug 2025).

For Stieltjes-type, time-scale, or difference equations, suitably generalized Wronskians (with respect to nontrivial “derivators”) obey analogous Abel-type evolution identities, certifying linear independence and enabling explicit solution construction for problems with discontinuities or time-scale structure (Fernández et al., 2022).

3. Wronskian Evolutions: Algebraic and Geometric Aspects

Wronskian evolutions bridge the analytic and geometric/arithmetic worlds. The KP hierarchy and its tau functions are parametrized by points in infinite-dimensional Grassmannians, with initial data supplied by flag varieties via the Mukhin–Varchenko–Schechtman (MVS) Wronskian map. Given a flag $F$ (represented by $g \in \operatorname{GL}_N$), Wronskian polynomials

$$y_i(x;g) = \operatorname{Wr}(b_1(\cdot;g), \dots, b_i(\cdot;g))$$

(with $b_i(x;g)$ polynomial bases formed from $g$) specify all initial data for the KP flows. The full integrable evolution is recovered by exponentiating with the KP times:

$$p_i(x;t) = \exp\left(\sum_{k \ge 1} t_k \partial_x^k\right) b_i(x;g),\qquad \tau(t) = \operatorname{Wr}(p_1,\dots,p_n).$$

Under KP flows, this reproduces all tau-functions, and the Wronskian image is a contraction of the universal Plücker (Grassmannian) coordinates, encoding both the algebraic and time-evolution structure of integrable PDEs (Gorbounov et al., 2020).

4. Combinatorics and Equivalence of Wronskian Evolutions

Wronskian evolutions exhibit deep combinatorial identities, particularly in the context of rational Darboux transformations and exceptional orthogonal polynomials. In the Hermite case, repeated state-adding and state-deleting transformations (encoded in Maya diagrams and partitions) generate a class of so-called pseudo-Wronskian determinants. All Wronskians corresponding to shifts of the origin in the Maya diagram are proportional, the minimal order representative realized at the maximal Durfee rectangle (inside corners of Ferrers diagrams)—this structure fully governs exceptional Hermite and Painlevé IV solutions (Gomez-Ullate et al., 2016).

More generally, in Laguerre and Jacobi cases, pseudo-Wronskian determinants are indexed by pairs of Maya diagrams, yielding even richer equivalence relations due to additional symmetry operations (shape invariance, discrete symmetries). Every pseudo-Wronskian in the same equivalence class corresponds to the same rational extension of the potential, with explicit proportionality relations and parameter shifts (Gomez-Ullate et al., 2018).

5. Discrete, Semi-Discrete, and Toda-Type Wronskian Evolutions

Wronskian evolutions admit discrete and semi-discrete analogues (Casoratians, mixed Wronskian–Casoratian forms). In such settings (e.g., the constrained discrete KP hierarchy and Ablowitz–Ladik lattice), evolution is driven by discrete flows, gauge/dressing transformations, and Sato-type difference equations. Wronskian determinants

$$\tau_\Delta(n,t) = \det [\Delta^{i-1} f_j(n,t)]_{i,j=1}^m$$

evolve through explicit shifts (e.g., $\Delta^k$), with the rank-restriction criteria ensuring that the constructed $\tau$-functions generate true integrable soliton solutions (Y-type, dark-bright, etc.) (Li et al., 2014).

These structures connect with 2D Toda systems. The evolution of principal minors of Wronskians (or their discrete analogues) is governed by classical Toda equations:

$$\partial_x\partial_y \ln w_n = \frac{w_{n+1}w_{n-1}}{w_n^2},$$

valid in both continuous, semi-discrete, and lattice settings. These recurrences yield the full system of integrals and general solutions for Toda chains and their reductions (Demskoi et al., 2014).

6. Wronskian Evolutions in Supersymmetric Quantum Mechanics

In $k$-confluent SUSY quantum mechanics, Wronskian determinantal formulas provide explicit solutions to higher-order intertwining of Schrödinger operators. In the confluent limit (all factorization energies coincident), Jordan chain solutions $\{u_j(x)\}$ and their parametric derivatives constitute the set of entries in the $k$-fold Wronskian:

$$V_k(x) = V_0(x) - 2\,\partial_x^2 \ln W(u_1, \dots, u_k).$$

Such constructions enable the classification and explicit realization of non-singular SUSY partner potentials, bypassing earlier integral-based constructions (Bermudez, 2015).

7. Summary of Structural Roles and Unifying Principles

Wronskian evolutions furnish:

• Explicit $\tau$-function solutions for integrable PDEs/P$\Delta$Es and their reductions, with time/discrete evolution governed by LDES for the entries.
• Algebraic invariants serving as integrals of motion, conserved under linear or variational flow (ODE, classical/quantum settings).
• The geometric bridges between flag varieties, Grassmannians, and the algebraic modulation of integrable system initial data.
• The combinatorial underpinning of equivalence classes for pseudo-Wronskians, directly controlling the structure and minimality of exceptional function spaces.
• Closed-form solutions to higher-order soliton, SUSY, and Toda/equivalent systems, unifying disparate integrable models.

Table: Representative Wronskian Evolutions Across Different Domains

Domain	Wronskian Type / Role	Reference
Soliton PDE hierarchies	$\tau=\operatorname{Wr}(\varphi_1, ..., \varphi_N)$, evolution via LDES and Hirota bilinear equations	(Zhang, 2019, Dai et al., 2023, Zhang et al., 2012)
Variational ODEs	$\det \Phi(t)$ conserved for traceless ODE systems	(Koshcheev, 15 Aug 2025)
Discrete / Semi-discrete integrable systems	Casoratian, discrete Wronskian determinants as $\tau$-functions, Toda recurrences	(Li et al., 2014, Demskoi et al., 2014)
KP hierarchy / algebraic geometry	Wronskian map from flags to tau-functions	(Gorbounov et al., 2020)
SUSY quantum mechanics	$k$-fold confluent Wronskians for SUSY partner potentials	(Bermudez, 2015)
Orthogonal polynomials, combinatorics	Pseudo-Wronskians, equivalence via partitions, Durfee rectangles	(Gomez-Ullate et al., 2016, Gomez-Ullate et al., 2018)

Wronskian evolutions thus provide a mathematically robust framework for encoding, evolving, and classifying explicit solutions, invariants, and symmetries across continuous and discrete integrable systems, variational mechanics, and algebraic combinatorics.