{
 "artifacts": null,
 "category": "dynamics",
 "description": "Replace an unconstrained recurrent or state-space transition with an energy-metric symmetric part plus an exactly energy-skew part. The symmetric component controls contraction or damping, while the skew component performs transport and mixing without changing the quadratic energy, enabling deeper networks and larger explicit integration steps.",
 "formulas_latex": [
  "$$A_{0}(V)\\partial_{t}V+A_{1}(V)\\partial_{x}V+A_{2}\\partial_{x}^{2}V=S_{\\phi}(V),\\qquad A_{0}(V)=\\begin{pmatrix}1\u00260\u00260\\\\0\u0026\\frac{1}{\\eta}\u00260\\\\0\u00260\u00261\\end{pmatrix},\\quad A_{1}(V)=\\begin{pmatrix}0\u00261\u00260\\\\1\u0026\\frac{2\\omega}{\\eta^{2}}\u00260\\\\0\u00260\u00260\\end{pmatrix},\\quad A_{2}=\\begin{pmatrix}0\u00260\u00260\\\\0\u00260\u0026-1\\\\0\u00261\u00260\\end{pmatrix}.$$",
  "$$\\mathfrak{E}(t)=\\sum_{i+j\\leq 2}\\left(A_{0}(V)\\partial_{t}^{i}\\partial_{x}^{j}V,\\partial_{t}^{i}\\partial_{x}^{j}V\\right)_{L^{2}}.$$",
  "$$B_{1}=\\underbrace{\\begin{pmatrix}\\frac{2m_{1}}{h^{2}}\u0026\\frac{m_{2}}{2h^{2}}\\\\\\frac{m_{2}}{2h^{2}}\u0026\\frac{m_{1}}{h^{2}}\\end{pmatrix}}_{B_{\\mathrm{1,sym}}}+\\frac{m_{2}}{2h^{2}}J,\\qquad J=\\begin{pmatrix}0\u00261\\\\-1\u00260\\end{pmatrix}.$$",
  "$$\\dot z=-(G^{-1}S+G^{-1}K)z+G^{-1}b,\\quad G=G^\\top\\succ0,\\ S=S^\\top\\succeq0,\\ K=-K^\\top,\\quad \\frac{d}{dt}\\frac12z^\\top Gz=-z^\\top Sz+z^\\top Kz=-z^\\top Sz.$$"
 ],
 "id": 38,
 "implementation": "Integrate this at the transition or mixing operator of a recurrent network, linear state-space model, or graph message-passing block. Let the hidden state be \\(z\\in\\mathbb R^d\\). Parameterize \\(G=LL^\\top+\\delta I\\), with lower-triangular \\(L\\) and \\(\\delta\u003e0\\); parameterize \\(S=RR^\\top\\); and parameterize \\(K=Q-Q^\\top\\). For an input-dependent transition, produce \\(b_t=W_xx_t\\) and optionally let diagonal gates modulate \\(R,Q\\) while preserving these constraints. Use an explicit midpoint or Euler step: (1) compute \\(v=G^{-1}(Sz+Kz-b_t)\\) by Cholesky solve; (2) update \\(z\\leftarrow z-\\Delta t\\,v\\); (3) apply the output projection and nonlinear readout, but do not insert an unconstrained activation inside the conservative transition. For the simplest MVP, use fixed diagonal \\(G\\), low-rank \\(K=UV^\\top-VU^\\top\\), and diagonal \\(S\\), making the transition cheap. The paper's mathematics supplies the exact symmetry/skew constraints and the energy identity; estimate the empirical energy ratio \\(E_{t+1}/E_t\\), where \\(E_t=z_t^\\top Gz_t/2\\), and tune \\(\\Delta t\\) using the largest observed stable step. Compare against a standard GRU, diagonal SSM, and unconstrained residual state-space layer on sequential MNIST and a small long-range-copy task using equal hidden size and FLOPs. Success signals are non-exploding hidden norms, stable training at a larger step size, lower gradient variance, and improved accuracy at sequence lengths where the unconstrained baseline becomes unstable. Ablate \\(K\\), \\(S\\), and learned \\(G\\) separately to verify that gains come from the decomposition rather than parameter count.",
 "math_summary": "The paper starts from the shallow-water system \\(\\partial_t h+\\operatorname{Div}(hu)=0\\) and \\(\\partial_tu+u\\cdot\\nabla u+\\nabla h-\\nabla\\Delta h=0\\), where \\(h\\) is height and \\(u\\) is velocity. Its generalized one-dimensional form is \\(A_0(V)\\partial_tV+A_1(V)\\partial_xV+A_2\\partial_x^2V=S_\\phi(V)\\), with \\(V=(\\eta-1,\\omega,\\phi)^\\top\\), \\(A_0(V)=\\operatorname{diag}(1,1/\\eta,1)\\) positive definite when \\(\\eta\u003e0\\), \\(A_1(V)=\\begin{pmatrix}0\u00261\u00260\\\\1\u00262\\omega/\\eta^2\u00260\\\\0\u00260\u00260\\end{pmatrix}\\) symmetric, and \\(A_2=\\begin{pmatrix}0\u00260\u00260\\\\0\u00260\u0026-1\\\\0\u00261\u00260\\end{pmatrix}\\) skew-symmetric. The associated energy is \\(\\mathfrak E(t)=\\sum_{i+j\\le2}(A_0(V)\\partial_t^i\\partial_x^jV,\\partial_t^i\\partial_x^jV)_{L^2}\\). In two dimensions, the nonsymmetric velocity matrices are decomposed as \\(B_k=B_{k,\\mathrm{sym}}+B_{k,\\mathrm{skw}}\\), isolating the skew contribution; for example \\(B_1=\\begin{pmatrix}2m_1/h^2\u0026m_2/h^2\\\\m_2/h^2\u0026m_1/h^2\\end{pmatrix}\\) is rewritten as \\(B_{1,\\mathrm{sym}}+\\frac{m_2}{2h^2}J\\), where \\(J=\\begin{pmatrix}0\u00261\\\\-1\u00260\\end{pmatrix}\\). The ML adaptation uses a positive-definite learned metric \\(G\\), a symmetric dissipative matrix \\(S=S^\\top\\succeq0\\), and a skew matrix \\(K=-K^\\top\\), giving \\(\\dot z=-(G^{-1}S+G^{-1}K)z+G^{-1}b\\); the \\(G\\)-energy \\(E(z)=z^\\top Gz/2\\) changes as \\(\\dot E=-z^\\top Sz+z^\\top Kz=-z^\\top Sz\\), because \\(z^\\top Kz=0\\).",
 "math_tags": [
  "dynamical-systems",
  "linear-algebra",
  "pde"
 ],
 "ml_areas": [
  "ssm",
  "rnn",
  "graph-nn"
 ],
 "paper": {
  "arxiv_id": "2608.21704",
  "arxiv_url": "https://arxiv.org/abs/2608.21704",
  "summary_what_math_gives_to_ml": "The paper exposes a structure-preserving decomposition for coupled dynamics: a positive symmetrizer defines the energy metric, while the principal transport and surface-tension operators are arranged so that their contributions are energy-neutral rather than amplifying. In two dimensions, the nonsymmetric velocity coupling is explicitly split into symmetric and skew parts, and the specific vorticity \\(\\theta=\\operatorname{curl}u/h\\) is exactly advected. These constructions can be transferred to recurrent, state-space, graph, or sequence-mixing layers by parameterizing the interaction operator as a sum of energy-symmetric damping and energy-skew transport. The resulting layers should permit substantially larger stable step sizes and reduce exploding or vanishing activations without relying only on spectral normalization or small residual scales.",
  "title": "Local Well-Posedness of the Inviscid Shallow Water Equations with Surface Tension",
  "year": "2026"
 },
 "ratings": {
  "difficulty": 5,
  "novelty": 6,
  "usefulness": 7
 },
 "solves": [
  "stability",
  "accuracy",
  "scalability"
 ],
 "title": "Symmetrized Energy-Skew Neural Dynamics",
 "url": "https://synthcore.org/idea/38/symmetrized-energy-skew-neural-dynamics",
 "verification": {
  "peer_reviewed": false,
  "status": "unverified",
  "status_label": "Unverified",
  "verdict_source": "deterministic test code (paired-seed permutation statistics)",
  "verification_axes": {
   "benchmark_mechanism": {
    "confirmed": null,
    "tested": false
   },
   "practical_benchmark": {
    "beats_baseline": null,
    "tested": false
   },
   "toy_mechanism_gate": {
    "confirmed": null,
    "tested": false
   }
  }
 }
}