# Lie-Poisson Hamiltonian latent block

- ID: 1255
- Canonical URL: https://synthcore.org/idea/1255/lie-poisson-hamiltonian-latent-block
- API JSON: https://synthcore.org/api/idea/1255.json
- API Markdown: https://synthcore.org/api/idea/1255.md
- Verification status: failed_benchmark
- Source: [arXiv:2607.21813](https://arxiv.org/abs/2607.21813)
- Category: dynamics
- Solves: stability, accuracy
- ML areas: ssm, rnn, world-model
- Math tags: dynamical-systems, geometry, differential-geometry, linear-algebra
- Ratings: usefulness 7/10; difficulty 5/10; novelty 6/10

## Idea description

Replace an unconstrained latent ODE or recurrent update with Hamiltonian dynamics on a product of Euclidean coordinates and a Lie-algebra momentum. The momentum dynamics contain the explicit coadjoint term generated by the Lie-group structure constants, allowing the model to represent rotational or frame-dependent memory without learning this antisymmetric coupling from data.

## Mathematical statement

Let q ∈ ℝⁿ be latent position, p ∈ ℝⁿ its momentum, and Hθ(q,p) a learned scalar Hamiltonian. Let cᵃᵢⱼ be Lie-algebra structure constants defined by [bᵢ,bⱼ] = cᵏᵢⱼ bₖ. The paper's Poisson bivector is 𝒥 = Dqᵢ ⊗ Dpⁱ − Dpⁱ ⊗ Dqᵢ + pₐ cᵃᵢⱼ Dpʲ ⊗ Dpⁱ, giving ż = 𝒥(z)∇Hθ(z). In coordinates, q̇ᵢ = ∂Hθ/∂pᵢ and ṗᵢ = −∂Hθ/∂qᵢ + pₐ cᵃᵢⱼ ∂Hθ/∂pⱼ. For SO(3), cᵃᵢⱼ = εₐᵢⱼ, so the additional term is a cross product. If Hθ is autonomous, dHθ/dt = ∇Hθᵀ𝒥∇Hθ = 0 because 𝒥 is antisymmetric.

## Key formulas

- 𝒥 = Dqᵢ ⊗ Dpⁱ − Dpⁱ ⊗ Dqᵢ + pₐ cᵃᵢⱼ Dpʲ ⊗ Dpⁱ
- q̇ᵢ = ∂Hθ/∂pᵢ,    ṗᵢ = −∂Hθ/∂qᵢ + pₐ cᵃᵢⱼ ∂Hθ/∂pⱼ
- dHθ/dt = ∇Hθᵀ𝒥∇Hθ = 0,    𝒥ᵀ = −𝒥
- ż = 𝒥(z)∇Hθ(z),    z = (q,p)

## Implementation notes

Use this block as the transition function inside a latent state-space model, Hamiltonian RNN, or world model. Encode an observation into q₀ and p₀, evaluate a neural scalar Hamiltonian Hθ(q,p) with an MLP, and integrate the resulting ODE for one or more prediction steps. For rotational latents, use n = 3 and hard-code cᵃᵢⱼ = εₐᵢⱼ; for a general Lie algebra, store c[a,i,j] and enforce antisymmetry in i and j. Pseudocode is: compute H; obtain qdot = gradient of H with respect to p; obtain pdot = minus the gradient of H with respect to q; add pdot[i] += p[a] c[a,i,j] qdot[j]; then advance z with implicit midpoint or a Lie-group-aware symplectic method. For implicit midpoint solve z_next = z + Δt 𝒥((z + z_next)/2) gradient H((z + z_next)/2) using three to five Newton iterations. The structure constants are computed from the chosen group, gradients use automatic differentiation, and energy drift is measured empirically. First test a synthetic rigid-body or rotating-pendulum sequence against a GRU and a standard Hamiltonian neural ODE at equal parameter count, using 100 to 500 step rollouts. The prediction is bounded oscillatory energy error of order Δt² with midpoint integration, whereas Euler has secular drift of order TΔt. Removing the coadjoint term should increase orientation phase error at nonzero momentum, while both variants should agree when p = 0.

## Verification

- Status: failed_benchmark
- Mechanism evidence: yes
- Mechanism confirmed: no
- Practical verdict: harms
- Verdict: Built a self-contained SO(3) Lie-Poisson rigid-body experiment with explicit coadjoint coupling, implicit midpoint integration, Euler baseline, and no-coadjoint ablation. The core predictions were confirmed: antisymmetry and energy conservation held to machine precision, the coadjoint effect vanished at p=0, momentum scaling was 2.06 versus predicted 2, and nonlinear midpoint state error showed approximately second-order timestep scaling.

### Mechanism check

- Verdict: Built a self-contained SO(3) Lie-Poisson rigid-body experiment with explicit coadjoint coupling, implicit midpoint integration, Euler baseline, and no-coadjoint ablation. The core predictions were confirmed: antisymmetry and energy conservation held to machine precision, the coadjoint effect vanished at p=0, momentum scaling was 2.06 versus predicted 2, and nonlinear midpoint state error showed approximately second-order timestep scaling.
- Confidence: 8/10
- Limitations: This was a toy analytic rigid-body test, not a learned latent neural block. No GRU, neural ODE, encoder, parameter matching, dataset benchmark, GPU benchmark, or wall-clock/FLOP comparison was tested.

### Practical benchmark

- Paired seeds: 8
- Baseline mean: 0
- Idea mean: 0
- p-value: 0.0081
- Paired wins: 0/8
- Benchmark verdict: idea worse (significant)

## Artifacts

- [bench_report.json](https://synthcore.org/code/503/bench_report.json)
- [bench_run.py](https://synthcore.org/code/503/bench_run.py)
- [experiment.py](https://synthcore.org/code/503/experiment.py)
- [nonlinear_check.py](https://synthcore.org/code/503/nonlinear_check.py)
- [nonlinear_results.json](https://synthcore.org/code/503/nonlinear_results.json)
- [report.md](https://synthcore.org/code/503/report.md)
- [report_bench_2026-08-31T162049.md](https://synthcore.org/code/503/report_bench_2026-08-31T162049.md)
- [results.json](https://synthcore.org/code/503/results.json)
- [Download all files as ZIP](https://synthcore.org/download/503)

## Disclaimer

AI-generated research hypothesis, automatically tested. Not peer-reviewed.