Dual Cubic String and Novikov Peakon Dynamics

• Dual Cubic String is a third-order spectral framework reformulated as a 3×3 first-order system with discrete mass measures on a finite interval.
• It underpins Novikov peakon dynamics via isospectral deformations, yielding explicit invariants and inverse maps through determinants and Pfaffians.
• The framework also connects dual scattering problems and interpolates between Toda-type integrable lattices through Miura-type reductions and bilinear relations.

Searching arXiv for papers on dual cubic string, Novikov peakons, and related cubic-string inverse scattering. Dual cubic string denotes a boundary-value and inverse-spectral framework that is formally paired with the cubic string and is linked, in the cited literature, to the Novikov equation, pure peakon dynamics, and integrable lattice hierarchies. In the most explicit formulation, the dual cubic string is posed on the interval $-1<\tilde y<1$ with a discrete positive mass measure and a first-order $3\times 3$ differential system whose spectral data are encoded by Weyl functions with partial-fraction expansions (Chang, 5 Sep 2025). A distinct but related usage appears in inverse scattering for a cubic string with a step-shaped potential, where the “dual” problem refers to scattering of waves coming from $-\infty$ rather than $+\infty$ (Zolotarev, 8 Sep 2025). Taken together, these works place the dual cubic string at the intersection of spectral theory for third-order string-type operators, peakon isospectral deformations, determinant/Pfaffian solution theory, and dual scattering formalisms (Chang, 5 Sep 2025, Zolotarev, 8 Sep 2025).

1. Spectral boundary-value formulation

The dual cubic-string boundary problem is presented in terms of a spectral parameter $z\in\mathbb C$ and an unknown vector

$$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$

with a weight measure on $-1<\tilde y<1$ of the form

$$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$

(Chang, 5 Sep 2025). The governing first-order system is

$$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$

with boundary conditions

$$\tilde\phi_2(-1)=\tilde\phi_3(-1)=0,\qquad\tilde\phi_3(1)=0$$

(Chang, 5 Sep 2025). The same source characterizes this as a first-order reformulation of a third-order string.

A discrete transfer description is given interval by interval. On each gap $3\times 3$0,

$3\times 3$1

while across a mass point $3\times 3$2,

$3\times 3$3

(Chang, 5 Sep 2025). This formulation makes the discrete geometry explicit: the problem is encoded by interval lengths $3\times 3$4 and point masses $3\times 3$5.

The associated Weyl functions are

$3\times 3$6

and they admit the partial-fraction expansions

$3\times 3$7

(Chang, 5 Sep 2025). In this presentation, the data $3\times 3$8 are the spectral invariants. This establishes the dual cubic string as an inverse-spectral object: the geometric variables $3\times 3$9 and the spectral variables $-\infty$0 are two coordinate systems for the same problem.

2. Relation to the Novikov equation and pure peakons

The cited literature states that the Novikov equation “can be formally regarded as linked to the dual cubic string,” and develops this connection through isospectral deformation (Chang, 5 Sep 2025). A compatible time flow is imposed by

$-\infty$1

in such a way that the spectrum $-\infty$2 remains fixed and

$-\infty$3

(Chang, 5 Sep 2025). The same source identifies this flow exactly with the Novikov peakon ODE system.

The PDE-side Lax representation is given as

$-\infty$4

and, with the peakon ansatz

$-\infty$5

a Liouville transform $-\infty$6 recovers the dual-string system (Chang, 5 Sep 2025). This places the dual cubic string within an isospectral representation of the Novikov dynamics.

In the pure-peakon sector, the ansatz is

$-\infty$7

with ODEs

$-\infty$8

$-\infty$9

(Chang, 5 Sep 2025). These ODEs are stated to be exactly the condition $+\infty$0 under the inverse-spectral map $+\infty$1.

A broader structural claim in the same work is that there is a bijective relationship between the DP and Novikov pure peakon trajectories, and also a one-to-one correspondence between the corresponding discrete cubic and dual cubic boundary value problems (Chang, 5 Sep 2025). This suggests that the dual cubic string is not merely an auxiliary spectral device for Novikov peakons but part of a paired architecture relating two third-order peakon integrable systems.

3. Spectral invariants and constants of motion

For the Novikov peakon flow associated with the dual cubic string, the polynomial

$+\infty$2

has coefficients $+\infty$3 that remain constant under the flow (Chang, 5 Sep 2025). These coefficients therefore serve as constants of motion expressed directly in the dual-string spectral polynomial.

The same source gives a closed-form formula: $+\infty$4 with $+\infty$5 so that $+\infty$6 (Chang, 5 Sep 2025). The paper also gives an equivalent characterization: each $+\infty$7 equals the sum of all $+\infty$8 minors of the $+\infty$9 matrix $z\in\mathbb C$0 (Chang, 5 Sep 2025).

This part of the theory is significant because it converts abstract spectral conservation into explicit algebraic invariants in peakon coordinates. The paper further emphasizes that these are “not previously known, explicit expressions for the constants of motion in the Novikov peakon dynamical system” (Chang, 5 Sep 2025). A plausible implication is that the dual cubic string provides a more directly computable invariant framework for Novikov peakons than formulations expressed solely at the PDE level.

4. Discrete–continuous inverse correspondence

A central result is a one-to-one correspondence between the discrete dual string and the continuous peakon formulation (Chang, 5 Sep 2025). The forward spectral map sends the discrete data $z\in\mathbb C$1 to the spectral variables $z\in\mathbb C$2 via the Weyl-function expansions. The inverse map is then written using bimoment determinants and Pfaffians derived from the discrete spectral measure

$z\in\mathbb C$3

(Chang, 5 Sep 2025).

The determinants are

$z\in\mathbb C$4

and the Pfaffians are

$z\in\mathbb C$5

with

$z\in\mathbb C$6

(Chang, 5 Sep 2025).

From these quantities, the inverse map is expressed explicitly by

$z\in\mathbb C$7

$z\in\mathbb C$8

$z\in\mathbb C$9

(Chang, 5 Sep 2025).

These formulas are not merely reconstruction identities; they place the dual cubic string in a determinant/Pfaffian framework. This suggests a structural affinity with BKP/CKP-type tau-function methods, although such a classification is not explicitly stated in the cited data. What is explicit is that determinants $$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$0 and Pfaffians $$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$1 provide a unified language for spectral data, inverse maps, and associated lattice flows (Chang, 5 Sep 2025).

5. Interpolating lattice and Toda-type connections

The same determinant/Pfaffian formalism yields a new integrable lattice connected to the dual cubic string (Chang, 5 Sep 2025). Writing

$$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$2

the paper derives the bilinear relations

$$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$3

$$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$4

$$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$5

(Chang, 5 Sep 2025).

New field variables are then introduced: $$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$6 leading to the four-field system

$$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$7

$$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$8

(Chang, 5 Sep 2025). Eliminating $$\tilde\Phi(\tilde y;z)=(\tilde\phi_1,\tilde\phi_2,\tilde\phi_3)^T,$$9 gives the two-field lattice

$-1<\tilde y<1$0

(Chang, 5 Sep 2025).

Two Miura-type reductions are then described. Under

$-1<\tilde y<1$1

the system becomes the classical B–Toda lattice

$-1<\tilde y<1$2

while with

$-1<\tilde y<1$3

one recovers the C–Toda lattice

$-1<\tilde y<1$4

(Chang, 5 Sep 2025). The paper therefore describes the two-field lattice as interpolating between the B–Toda and C–Toda flows.

Within the stated framework, this interpolation is one of the strongest indications that the dual cubic string is embedded in a wider hierarchy of integrable structures. The determinant/Pfaffian construction unifies spectral data, peakon dynamics, and Toda-type lattices rather than treating them as separate phenomena (Chang, 5 Sep 2025).

6. Dual scattering for step-shaped cubic strings

A different but related notion of duality appears in the inverse scattering problem for a cubic string having the shape of a step (Zolotarev, 8 Sep 2025). There, the governing equation is

$-1<\tilde y<1$5

with a real weight satisfying

$-1<\tilde y<1$6

and

$-1<\tilde y<1$7

for some $-1<\tilde y<1$8 (Zolotarev, 8 Sep 2025). The paper proves that in this setting there are two scattering problems: the direct one, for waves from $-1<\tilde y<1$9, and its dual, for waves from $$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$0 (Zolotarev, 8 Sep 2025).

For the dual problem, one introduces Jost solutions $$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$1, $$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$2, normalized at $$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$3 by

$$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$4

where

$$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$5

(Zolotarev, 8 Sep 2025). The physical expansion

$$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$6

is rewritten as

$$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$7

with

$$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$8

(Zolotarev, 8 Sep 2025). The corresponding unitarity relation is

$$\tilde g(\tilde y)=\sum_{k=1}^N\tilde g_k\,\delta(\tilde y-\tilde y_k),\quad -1=\tilde y_0<\tilde y_1<\cdots<\tilde y_N<\tilde y_{N+1}=1,\;\tilde g_k>0$$9

(Zolotarev, 8 Sep 2025).

The reconstruction scheme is based on a matrix Riemann–Hilbert problem on a three-ray star and yields a system of linear singular integral equations for unknown continuous functions $$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$0, $$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$1 and discrete residues $$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$2, $$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$3 (Zolotarev, 8 Sep 2025). The scattering data entering the problem are exactly

$$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$4

(Zolotarev, 8 Sep 2025). From the solution, one reconstructs

$$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$5

in the sector $$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$6, uses the asymptotic

$$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$7

where

$$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$8

and finally recovers the potential from

$$\frac{d}{d\tilde y}\tilde\Phi =\begin{pmatrix} 0 & \tilde g(\tilde y) & 0\ 0 & 0 & \tilde g(\tilde y)\ -\,z & 0 & 0 \end{pmatrix}\tilde\Phi,$$9

(Zolotarev, 8 Sep 2025).

This usage of “dual” differs from the discrete dual cubic-string boundary problem of the Novikov literature. In the step-potential setting, duality refers to the opposite incidence direction in scattering. The shared terminology nonetheless points to a broader pattern: cubic-string theory naturally supports paired formulations distinguished either by spectral/peakon correspondence or by left/right scattering orientation.

7. Conceptual scope and related structures

Across the cited works, the dual cubic string occupies two closely neighboring but nonidentical conceptual roles. In one role, it is a discrete boundary-value problem whose Weyl data generate the Novikov peakon isospectral flow and a determinant/Pfaffian inverse theory (Chang, 5 Sep 2025). In the other, it refers to the dual scattering problem for a step-shaped cubic string, describing waves incident from $$\tilde\phi_2(-1)=\tilde\phi_3(-1)=0,\qquad\tilde\phi_3(1)=0$$0 and reconstructed through singular integral equations and a star-contour Riemann–Hilbert formalism (Zolotarev, 8 Sep 2025).

The literature explicitly contrasts the dual cubic string with the cubic string through the statement that the DP equation can be viewed as an isospectral deformation of the boundary value problem for the cubic string, while the Novikov equation can be formally regarded as linked to the dual cubic string (Chang, 5 Sep 2025). It also states that there is a one-to-one correspondence between the corresponding discrete cubic and dual cubic boundary value problems (Chang, 5 Sep 2025). This suggests that “dual cubic string” is best understood not as a single isolated model but as one member of a paired third-order spectral framework.

A common misconception would be to identify the dual cubic string exclusively with either Novikov peakons or dual scattering from $$\tilde\phi_2(-1)=\tilde\phi_3(-1)=0,\qquad\tilde\phi_3(1)=0$$1. The cited sources do not support such a reduction. Instead, they document two technical contexts in which duality arises: one algebraic and isospectral, the other scattering-theoretic (Chang, 5 Sep 2025, Zolotarev, 8 Sep 2025). Another plausible implication is that the determinant/Pfaffian machinery developed for the discrete dual cubic string may inform future treatments of more general dual scattering problems, but that extension is not stated explicitly in the cited data.

In the present arXiv literature, the dual cubic string is therefore characterized by three main features: a third-order spectral structure recast as a $$\tilde\phi_2(-1)=\tilde\phi_3(-1)=0,\qquad\tilde\phi_3(1)=0$$2 first-order system, an exact correspondence with Novikov pure-peakon isospectral evolution, and a broader duality principle in cubic-string scattering theory (Chang, 5 Sep 2025, Zolotarev, 8 Sep 2025).