---
title: King Differential Equation Overview
url: https://www.emergentmind.com/topics/king-differential-equation
type: topic
---

# King Differential Equation Overview

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 [1912.02299] [2210.04352] [2606.12455].

## 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 [1912.02299] [2210.04352] [2606.12455].

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 [1912.02299] [2606.12455].

## 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 [1912.02299].

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 [1912.02299].

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
\[
(k,l,m)\mapsto (l+m,k+m,k+l)
\]
carries the King–Schief lattice frame into the dual AKP lattice frame [1912.02299].

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 [1912.02299].

The paper’s main claim is that this connection **establishes the integrability of the dual AKP equation** and proves the existence of the conjectured **\(N\)-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 [1912.02299].

## 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
\[
\begin{cases}
u_t=\Delta u+u^3 & \text{in } \mathbb R^4\times (0,\infty),\\
u(x,0)=u_0(x) & \text{in } \mathbb R^4.
\end{cases}
\]
For the general equation
\[
u_t=\Delta u+|u|^{p-1}u\quad \text{in }\mathbb R^n,
\]
the energy-critical exponent is \(p=(n+2)/(n-2)\), so \(p=3\) when \(n=4\). Accordingly, \(u_t=\Delta u+u^3\) is the **energy-critical** case in \(\mathbb R^4\) [2210.04352].

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 \(\mu(t)\) as the object informally referred to as the “King differential equation” [2210.04352].

Wei, Zhang, and Zhou rigorously realize this scenario. They construct a positive global solution with
\[
\|u(\cdot,t)\|_{L^\infty(\mathbb R^4)}\sim \ln t\qquad \text{as } t\to\infty,
\]
and with the asymptotic blow-up form
\[
u(x,t)=\eta\!\left(\frac{x-\xi(t)}{\sqrt{t}}\right)\mu(t)^{-1}w\!\left(\frac{x-\xi(t)}{\mu(t)}\right)+\text{smaller terms},
\]
where the elliptic bubble is
\[
w(y)=\frac{2\sqrt{2}}{1+|y|^2},
\]
and the modulation parameters satisfy
\[
\mu(t)=\frac{1+o(1)}{\ln t},\qquad \xi(t)=O(t^{-1})\qquad \text{as } t\to\infty.
\]
Since \(w(0)=2\sqrt{2}\), this yields
\[
\|u(\cdot,t)\|_{L^\infty}\sim \mu(t)^{-1}\sim \ln t
\]
[2210.04352].

The modulation law arises by projecting the error of the bubble ansatz onto the kernel of the linearized operator
\[
L=\Delta+3w^2.
\]
The translational modes are
\[
Z_i(y)=\partial_{y_i}w(y),\qquad i=1,\dots,4,
\]
and the scaling mode is
\[
Z_5(y)=y\cdot \nabla w+2w.
\]
Orthogonality against \(Z_5\) produces a nonlocal equation for \(\mu(t)\), written in the paper in the canonical form
\[
p_1'(t)+B_{\nu_1}(t)p_1(t)=f_{\nu_1}[p_1](t),
\]
with \(p(t)=\mu(t)\) and \(B_{\nu_1}(t)\sim \frac1t\frac{(1-\nu_1)\ln t+2}{\ln t}\), leading to the asymptotic law \(p(t)\sim (\ln t)^{-1}\) [2210.04352].

The same paper proves **stability** of the infinite-time blow-up regime. In the nonradial setting it assumes perturbations with
\[
|g_0(x)|\le C_g |x|^{-\ell},\qquad \ell>3,
\]
while in the radial setting it states that the conjectured threshold \(\ell>2\) suffices and \(\xi[g_0]\equiv 0\). The authors further remark that extending the nonradial stability result to the full range \(\ell>2\) should be possible, although it is not carried out in detail [2210.04352].

## 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
\[
M_{\mathbf u,\varsigma}(\mathbf v)
=
\frac{1}{\pi^{3/2}\varsigma^{3}}
\exp\!\Bigl(-\frac{|\mathbf v-\mathbf u|^{2}}{\varsigma^{2}}\Bigr),
\]
whose laboratory-frame spherical harmonic expansion produces radial coefficients called **King functions**. For spherical harmonic degree \(l\), the radial kernel is
\[
\mathcal K_l(v;u,\varsigma)
=
\frac{1}{\pi^{3/2}\varsigma^{3}(2l+1)}
\exp\!\Bigl(-\frac{u^2+v^2}{\varsigma^2}\Bigr)
\, i_l\!\Bigl(\frac{2uv}{\varsigma^2}\Bigr),
\qquad u>0,
\]
where \(i_l\) is the modified spherical Bessel function [2606.12455].

Introducing the dimensionless variables
\[
r:=\frac{v}{\varsigma},\qquad k:=\frac{2u}{\varsigma},
\]
the essential radial profile becomes
\[
\Psi_{l,k}(r):=e^{-r^2} i_l(k r).
\]
This is the **real-parameter King function**. The corresponding **imaginary branch** is defined by \(k=i\kappa\), \(\kappa\ge0\), which yields
\[
\Phi_{l,\kappa}(r):=e^{-r^2} j_l(\kappa r),
\]
with \(j_l\) the spherical Bessel function [2606.12455].

The same work clarifies the relation between these kernels and the **co-moving Laguerre hierarchy**. In the co-moving frame, the radial eigenfunctions are
\[
R_{p,l}^{\mathrm{Lag}}(\widehat w)=\widehat w^l L_p^{(l+1/2)}(\widehat w^2),
\qquad p\in\mathbb N_0,
\]
and the King–Laguerre expansion shows that each King function is an infinite superposition of Laguerre modes:
\[
\mathcal K_l(r;k,\varsigma)
=
M_{0,\varsigma}(r)\,
\frac{\sqrt{\pi}}{2}(2l+1)
\Bigl(\frac{k}{2}\Bigr)^l
\sum_{p=0}^\infty
\frac{(-1)^p (k/2)^{2p}}{\Gamma(p+l+3/2)}
\,L_p^{(l+1/2)}(r^2),
\]
where \(M_{0,\varsigma}(r)=\pi^{-3/2}\varsigma^{-3}e^{-r^2}\) is the unshifted Maxwellian [2606.12455].

In this setting, the **King differential equation** is the ODE satisfied by
\[
f(r)=e^{-r^2} i_l(k r).
\]
Using the standard Bessel equation for \(i_l\), the paper derives
\[
f''(r)
+
\Bigl(\frac{2}{r}+4r\Bigr)f'(r)
+
\Bigl[
-\frac{l(l+1)}{r^2}+4r^2+6-k^2
\Bigr]f(r)=0,
\qquad r>0.
\]
Here the term “King equation” is exact rather than heuristic: it is a second-order linear differential equation with spectral parameter \(k^2\) [2606.12455].

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

The shifted-Gaussian King equation admits a natural Sturm–Liouville formulation. Defining
\[
\omega(r)=r^2 e^{2r^2},
\qquad
V_l(r)=\frac{l(l+1)}{r^2}-(4r^2+6),
\]
the differential expression is
\[
\tau_l f(r)
=
-\frac1{\omega(r)}
\frac{d}{dr}\Bigl(\omega(r)f'(r)\Bigr)
+
V_l(r)f(r),
\]
and the King equation becomes the eigenvalue problem
\[
\tau_l f = k^2 f
\]
[2606.12455].

The natural Hilbert space is
\[
H=L^2\big((0,\infty);\,r^2 e^{2r^2}\,dr\big).
\]
On this space, the maximal domain is
\[
\operatorname{Dom}_{\max}(\tau_l)
=
\Bigl\{
f\in H\ \big|\ 
f,\ \omega f'\in AC_{\mathrm{loc}}(0,\infty),\ \tau_l f\in H
\Bigr\}.
\]
The associated operator \(\mathcal L_l\) is defined by \(\mathcal L_l f:=\tau_l f\), with no extra boundary condition for \(l\ge1\), while for \(l=0\) the paper imposes
\[
\lim_{r\to 0^+} f(r)=0.
\]
With these domains, \(\mathcal L_l\) is self-adjoint on \(H\) [2606.12455].

The decisive structural result is a unitary equivalence with the free radial Schrödinger operator. The unitary map
\[
(\mathcal U f)(r)=r e^{r^2} f(r)
\]
sends \(H\) onto \(L^2(0,\infty)\), and the transformed operator is
\[
\mathcal U \tau_l \mathcal U^{-1}
=
h_l:=-\frac{d^2}{dr^2}+\frac{l(l+1)}{r^2}.
\]
Hence
\[
\mathcal L_l=\mathcal U^{-1} h_l\,\mathcal U.
\]
The King operator is therefore unitarily equivalent to the standard half-line free radial Schrödinger operator with angular momentum \(l\) [2606.12455].

This yields the full spectral description:
\[
\sigma(\mathcal L_l)=\sigma_{\mathrm{ac}}(\mathcal L_l)=[0,\infty),
\qquad
\sigma_{\mathrm{p}}(\mathcal L_l)=\sigma_{\mathrm{sc}}(\mathcal L_l)=\emptyset.
\]
The generalized eigenfunctions are given by the imaginary-parameter King functions
\[
\Phi_{l,\kappa}(r)=e^{-r^2} j_l(\kappa r),
\qquad \kappa\ge0,
\]
with generalized orthogonality
\[
\bigl\langle \Phi_{l,\kappa},\Phi_{l,\kappa'}\bigr\rangle_H
=
\frac{\pi}{2\kappa^2}\,\delta(\kappa-\kappa').
\]
Accordingly, every \(f\in H\) has the spectral representation
\[
f(r)=\int_0^\infty \widetilde f(\kappa)\,\Phi_{l,\kappa}(r)\,d\kappa,
\]
with
\[
\widetilde f(\kappa)
=
\frac{2\kappa^2}{\pi}
\int_0^\infty
f(r)\overline{\Phi_{l,\kappa}(r)}\,r^2 e^{2r^2}\,dr,
\]
and Parseval identity
\[
\|f\|_H^2
=
\int_0^\infty
|\widetilde f(\kappa)|^2\,
\frac{\pi}{2\kappa^2}\,d\kappa
\]
[2606.12455].

The paper then distinguishes the spectral role of the imaginary branch from the approximation-theoretic role of the real branch. In the radial space
\[
H_{\rm rad}
=
L^2\big((0,\infty);\,r^2 e^{-r^2}\,dr\big),
\]
the family of real-parameter King functions
\[
\mathscr K_l^{\mathbb R}
=
\{\Psi_{l,k}(r)=e^{-r^2}i_l(k r)\mid k>0\}
\]
forms a **dense non-orthogonal system**. Since these real-parameter King functions correspond to \(\lambda=-k^2<0\), they lie in the **resolvent set** rather than the spectrum of \(\mathcal L_l\), yet they still provide a dense approximation dictionary for **King mixture representations** [2606.12455].

## 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 [1912.02299] [2210.04352] [2606.12455].

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 \(\mu(t)\); and the shifted-Gaussian King equation characterizes the radial kernel \(e^{-r^2} i_l(k r)\) 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 [1912.02299] [2210.04352] [2606.12455].

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 \(N\)-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 [1912.02299] [2210.04352] [2606.12455].

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.

Source: https://www.emergentmind.com/topics/king-differential-equation