Complex Phase-Corrected State Integrator
Implementation & benchmark of arXiv:2608.22315 — Sharp CFL stability and temporal-dispersion optimization of symmetric splitting schemes for time-domain Maxwell equations
Source paper: Sharp CFL stability and temporal-dispersion optimization of symmetric splitting schemes for time-domain Maxwell equations arXiv:2608.22315 ⓘ · analyzed Aug 29, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Use the complex-conjugate palindromic coefficient that cancels the leading temporal phase defect of oscillatory modes. Implement complex arithmetic directly or use an exactly equivalent doubled-real state, then project the final state to its real component for real-valued prediction tasks.
Formulas
Mathematical statement
For the same palindromic family, \(g_2=a(1-2a)/2\). The leading temporal phase defect cancels when \(g_2=1/12\), but every real coefficient satisfies \(g_2\leq1/16\). Solving the cancellation equation gives \(a=1/4\pm i\sqrt{3}/12\). The paper reports the resulting stability threshold \(s_*=6\sqrt{3}/(7\sqrt d)\), and distinguishes fourth-order phase accuracy for each fixed semidiscrete Fourier mode from second-order global field accuracy. For a real neural state, use a complex state \(z=u+iv\), apply the palindromic composition with the complex coefficient, and return \(u=\operatorname{Re}(z)\). A real-only implementation stores \((u,v)\) and applies the block-real representation of every complex linear map.
Implementation notes
Integrate this into an oscillatory SSM or Hamiltonian neural ODE for long-context sequence prediction, especially when hidden dynamics encode rotations, Fourier modes, or wave-like memory. Start with a linear split state consisting of h and e, and set a=0.25+i*sqrt(3)/12. In complex PyTorch or JAX, implement each shear with complex scalar multiplication. In real-only code, duplicate every hidden channel into real and imaginary parts and implement a complex matrix W=A+iB as the real block matrix [[A,-B],[B,A]]. Use the same coefficient throughout the palindromic composition; do not learn unrelated real and imaginary paths. Estimate the coupling norm with power iteration and restrict the step to delta_t*L_hat <= 0.9*6*sqrt(3)/(7*sqrt(d)). For a cheap falsification experiment, train a small sinusoid-forecasting SSM and a copy-memory model with sequence lengths 1000 and 10000. Compare Euler, the real a=1/4 splitter, and the complex splitter at matched numbers of operator calls and matched parameter counts where possible. Measure phase error after 1000 and 10000 steps, spectral radius, hidden-state norm, gradient norms, training loss, and validation loss. The expected signal is lower accumulated phase drift and better long-horizon prediction at equal step size. Reject the method if doubled-real memory and compute erase the accuracy gain, if the real projection causes optimization instability, or if it does not outperform the real stability-optimized splitter on oscillatory data.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.