King Differential Equation Overview

• King differential equation is a context-sensitive term that describes a 14-point lattice relation, a nonlocal modulation law in critical heat flow, and a second-order ODE in shifted-Gaussian spectral theory.
• In discrete integrable systems, the King–Schief equation arises from the lattice BKP/AKP hierarchies, revealing equivalence with the dual AKP equation and inheriting soliton and multidimensional consistency properties.
• In analytical frameworks, it governs the blow-up scale in four-dimensional semilinear heat flow and defines the Sturm–Liouville ODE for King functions, establishing self-adjoint spectral properties via unitary equivalence with radial Schrödinger operators.

Searching arXiv for the cited works and closely related usage of “King differential equation.” The expression “King differential equation” is not a single canonical term with a uniform meaning across mathematical physics. In the arXiv literature, it appears in at least three distinct but internally precise senses: as the King–Schief 14-point lattice equation in discrete integrable systems, as the effective modulation equation governing the blow-up scale in the Fila–King analysis of the four-dimensional energy-critical semilinear heat equation, and as the explicit second-order linear ODE satisfied by the radial King function $e^{-r^2} i_l(k r)$ arising from shifted Gaussian kernels in spherical harmonic expansions (Kamp et al., 2019, Wei et al., 2022, Wang et al., 5 Jun 2026).

1. Terminological scope and principal usages

The literature represented here assigns the name to mathematically different objects, each tied to a specific research program rather than to a universal historical definition.

Context	Object called “King differential equation”	Defining setting
Discrete integrable systems	King–Schief equation	14-point lattice relation on $\mathbb Z^3$
Critical parabolic PDE	Effective King modulation equation	Nonlocal evolution law for the blow-up scale
Shifted Gaussian spectral theory	King equation	ODE for $e^{-r^2} i_l(k r)$

In the discrete setting, “King” refers to A.D. King, coauthor with W.K. Schief of the geometric work on Bianchi hypercubes from which the 14-point lattice equation arises. In the parabolic setting, “King” refers to Fila & King, whose matched asymptotic analysis predicted a logarithmic infinite-time blow-up law in dimension four. In the shifted-Gaussian setting, “King function” names the radial kernels in the laboratory-frame spherical harmonic expansion of a shifted Gaussian, and the corresponding “King differential equation” is the ODE they satisfy (Kamp et al., 2019, Wei et al., 2022, Wang et al., 5 Jun 2026).

A common misconception is that the phrase necessarily denotes a continuum differential equation. That is incorrect in the integrable-systems usage, where the central object is explicitly a difference equation on a lattice, and only in the 2026 shifted-Gaussian work does the term denote an actual second-order ODE (Kamp et al., 2019, Wang et al., 5 Jun 2026).

2. The King–Schief equation in discrete integrable systems

In van der Kamp, Zhang, and Quispel’s treatment, the relevant object is the King–Schief equation, a scalar lattice equation for $\tau(k,l,m)$ on the three-dimensional integer lattice with shifts

$$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$

Its full form is a 14-point equation involving 14 distinct lattice sites and parameters $A,B,C,D$. The same construction has a 12-point reduction when $D=0$, yielding an AKP-related equation (Kamp et al., 2019).

The paper states that the 14-point King–Schief equation was obtained by King and Schief as a consequence of the lattice BKP equation, while the 12-point reduction is tied to the lattice AKP equation. In this framework, the King–Schief equation is not introduced as an ad hoc stencil but as a reduced relation emerging from the bilinear BKP/AKP hierarchy. The geometrical origin is in King and Schief’s description via “Bianchi hypercubes” and a unification of Hirota’s AKP and Miwa’s BKP equations (Kamp et al., 2019).

The central structural fact is its equivalence to the dual AKP equation. Van der Kamp et al. show that if $\tau(k,l,m)$ solves the 14-point King–Schief equation, then

$\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$

solves the dual AKP equation under the parameter identification

$$a_1=A^2,\qquad a_2=B^2,\qquad a_3=C^2,\qquad a_4=-D^2.$$

Thus the coordinate transformation

$\mathbb Z^3$0

carries the King–Schief lattice frame into the dual AKP lattice frame (Kamp et al., 2019).

This equivalence has two immediate consequences. First, the full King–Schief 14-point equation corresponds to the dual AKP equation with all parameters present, while the 12-point reduction corresponds to the case where one parameter vanishes. Second, the equivalence transfers integrability properties: since the King–Schief equation is a consequence of BKP or AKP, the dual AKP equation inherits the relevant soliton structure, multidimensional consistency, and Bäcklund-type features through this identification (Kamp et al., 2019).

The paper’s main claim is that this connection establishes the integrability of the dual AKP equation and proves the existence of the conjectured $\mathbb Z^3$1-soliton solution for all values of the parameters. Because the King–Schief and dual AKP equations are equivalent, the same conclusion applies to the King–Schief equation itself (Kamp et al., 2019).

3. The Fila–King modulation equation in four-dimensional critical heat flow

In the PDE literature represented by Wei, Zhang, and Zhou, “King differential equation” refers not to the primary PDE but to the effective evolution law for the slowly varying scale parameter in the four-dimensional energy-critical semilinear heat equation

$\mathbb Z^3$2

For the general equation

$\mathbb Z^3$3

the energy-critical exponent is $\mathbb Z^3$4, so $\mathbb Z^3$5 when $\mathbb Z^3$6. Accordingly, $\mathbb Z^3$7 is the energy-critical case in $\mathbb Z^3$8 (Wei et al., 2022).

Fila and King studied the critical equation with radially symmetric positive initial data having prescribed power decay at infinity. In dimension four, their matched asymptotic analysis predicted global unbounded solutions exhibiting infinite-time blow-up with logarithmic rate, and the paper identifies the associated effective ODE, more precisely a nonlocal integro-differential equation, for the scale parameter $\mathbb Z^3$9 as the object informally referred to as the “King differential equation” (Wei et al., 2022).

Wei, Zhang, and Zhou rigorously realize this scenario. They construct a positive global solution with

$e^{-r^2} i_l(k r)$0

and with the asymptotic blow-up form

$e^{-r^2} i_l(k r)$1

where the elliptic bubble is

$e^{-r^2} i_l(k r)$2

and the modulation parameters satisfy

$e^{-r^2} i_l(k r)$3

Since $e^{-r^2} i_l(k r)$4, this yields

$e^{-r^2} i_l(k r)$5

(Wei et al., 2022).

The modulation law arises by projecting the error of the bubble ansatz onto the kernel of the linearized operator

$e^{-r^2} i_l(k r)$6

The translational modes are

$e^{-r^2} i_l(k r)$7

and the scaling mode is

$e^{-r^2} i_l(k r)$8

Orthogonality against $e^{-r^2} i_l(k r)$9 produces a nonlocal equation for $\tau(k,l,m)$0, written in the paper in the canonical form

$\tau(k,l,m)$1

with $\tau(k,l,m)$2 and $\tau(k,l,m)$3, leading to the asymptotic law $\tau(k,l,m)$4 (Wei et al., 2022).

The same paper proves stability of the infinite-time blow-up regime. In the nonradial setting it assumes perturbations with

$\tau(k,l,m)$5

while in the radial setting it states that the conjectured threshold $\tau(k,l,m)$6 suffices and $\tau(k,l,m)$7. The authors further remark that extending the nonradial stability result to the full range $\tau(k,l,m)$8 should be possible, although it is not carried out in detail (Wei et al., 2022).

4. King functions and the ODE arising from shifted Gaussian kernels

A third and fully literal usage appears in the study of shifted Gaussian distributions. There, the starting point is the shifted Maxwellian

$\tau(k,l,m)$9

whose laboratory-frame spherical harmonic expansion produces radial coefficients called King functions. For spherical harmonic degree $$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$0, the radial kernel is

$$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$1

where $$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$2 is the modified spherical Bessel function (Wang et al., 5 Jun 2026).

Introducing the dimensionless variables

$$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$3

the essential radial profile becomes

$$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$4

This is the real-parameter King function. The corresponding imaginary branch is defined by $$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$5, $$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$6, which yields

$$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$7

with $$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$8 the spherical Bessel function (Wang et al., 5 Jun 2026).

The same work clarifies the relation between these kernels and the co-moving Laguerre hierarchy. In the co-moving frame, the radial eigenfunctions are

$$\tilde{\tau}=\tau(k+1,l,m),\qquad \hat{\tau}=\tau(k,l+1,m),\qquad \dot{\tau}=\tau(k,l,m+1).$$9

and the King–Laguerre expansion shows that each King function is an infinite superposition of Laguerre modes: $A,B,C,D$0 where $A,B,C,D$1 is the unshifted Maxwellian (Wang et al., 5 Jun 2026).

In this setting, the King differential equation is the ODE satisfied by

$A,B,C,D$2

Using the standard Bessel equation for $A,B,C,D$3, the paper derives

$A,B,C,D$4

Here the term “King equation” is exact rather than heuristic: it is a second-order linear differential equation with spectral parameter $A,B,C,D$5 (Wang et al., 5 Jun 2026).

5. Sturm–Liouville structure, self-adjoint realization, and spectrum

The shifted-Gaussian King equation admits a natural Sturm–Liouville formulation. Defining

$A,B,C,D$6

the differential expression is

$A,B,C,D$7

and the King equation becomes the eigenvalue problem

$A,B,C,D$8

(Wang et al., 5 Jun 2026).

The natural Hilbert space is

$A,B,C,D$9

On this space, the maximal domain is

$D=0$0

The associated operator $D=0$1 is defined by $D=0$2, with no extra boundary condition for $D=0$3, while for $D=0$4 the paper imposes

$D=0$5

With these domains, $D=0$6 is self-adjoint on $D=0$7 (Wang et al., 5 Jun 2026).

The decisive structural result is a unitary equivalence with the free radial Schrödinger operator. The unitary map

$D=0$8

sends $D=0$9 onto $\tau(k,l,m)$0, and the transformed operator is

$\tau(k,l,m)$1

Hence

$\tau(k,l,m)$2

The King operator is therefore unitarily equivalent to the standard half-line free radial Schrödinger operator with angular momentum $\tau(k,l,m)$3 (Wang et al., 5 Jun 2026).

This yields the full spectral description: $\tau(k,l,m)$4 The generalized eigenfunctions are given by the imaginary-parameter King functions

$\tau(k,l,m)$5

with generalized orthogonality

$\tau(k,l,m)$6

Accordingly, every $\tau(k,l,m)$7 has the spectral representation

$\tau(k,l,m)$8

with

$\tau(k,l,m)$9

and Parseval identity

$\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$0

(Wang et al., 5 Jun 2026).

The paper then distinguishes the spectral role of the imaginary branch from the approximation-theoretic role of the real branch. In the radial space

$\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$1

the family of real-parameter King functions

$\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$2

forms a dense non-orthogonal system. Since these real-parameter King functions correspond to $\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$3, they lie in the resolvent set rather than the spectrum of $\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$4, yet they still provide a dense approximation dictionary for King mixture representations (Wang et al., 5 Jun 2026).

6. Comparative interpretation and recurring themes

Across these three usages, the phrase “King differential equation” designates objects that are formally different: a rational 14-point lattice equation, a nonlocal modulation equation for a blow-up scale, and a linear Sturm–Liouville ODE. The shared name is therefore historical and contextual, not classificatory in the strict sense (Kamp et al., 2019, Wei et al., 2022, Wang et al., 5 Jun 2026).

Several structural parallels nevertheless recur. Each usage isolates a reduced equation governing a distinguished mode or kernel: the King–Schief equation condenses information from BKP/AKP bilinear identities into a 14-point relation; the Fila–King modulation equation governs the slow evolution of the concentrating bubble scale $\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$5; and the shifted-Gaussian King equation characterizes the radial kernel $\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$6 underlying the spherical harmonic expansion of a shifted Maxwellian. This suggests a family resemblance centered on reduced governing equations extracted from richer ambient structures, although that interpretation is synthetic rather than terminologically fixed (Kamp et al., 2019, Wei et al., 2022, Wang et al., 5 Jun 2026).

Another frequent source of confusion is the status of “integrability” or “spectrality.” In the discrete case, integrability is meant in the discrete-integrable sense: existence of $\tilde{\tau}(k,l,m):=\tau(l+m,k+m,k+l)$7-soliton solutions, relation to BKP/AKP, multidimensional consistency, and Bäcklund transformations. In the four-dimensional heat-flow case, the central issue is not integrability but singularity formation, specifically stable infinite-time blow-up with logarithmic rate. In the shifted-Gaussian case, the key notion is self-adjoint spectral theory, culminating in unitary equivalence with a free radial Schrödinger operator and density of the real-parameter branch in a natural radial Hilbert space (Kamp et al., 2019, Wei et al., 2022, Wang et al., 5 Jun 2026).

Taken together, these works show that “King differential equation” is best understood as a context-sensitive label. In discrete integrable systems it points to the King–Schief lattice relation; in critical heat-flow analysis it denotes the effective Fila–King scale law; and in shifted-Gaussian analysis it denotes the exact ODE for the King radial profile. Any precise use of the term therefore requires specifying the surrounding framework, the dependent variable, and the operator or hierarchy from which the equation is derived.