Spectrally identifiable phaseless recurrent layer
Implementation & benchmark of arXiv:2607.26705 — Dynamical phase retrieval for Schr{ö}dinger evolution on finite graphs
Source paper: Dynamical phase retrieval for Schr{ö}dinger evolution on finite graphs arXiv:2607.26705 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Build a complex-valued recurrent or graph-neural layer whose hidden state evolves under a fixed graph Schrödinger operator and is exposed to the downstream network only through coordinate magnitudes at several times. Choose the diagonal potential so that the spectrum has unique unordered pair sums, the squared-eigenvector matrix is invertible, and every eigenvector pair overlaps in at least one observed coordinate; the resulting temporal intensity code is theoretically injective up to one unavoidable global phase.
Formulas
Mathematical statement
The paper studies the real symmetric operator $H_Q=\Delta_G+Q$, where $\Delta_G$ is the graph Laplacian and $Q=\operatorname{diag}(q_1,\ldots,q_n)$ is a real diagonal potential. For an initial complex state $u_0\in\mathbb{C}^n$, the observable is $y_j(t)=|e^{-itH_Q}u_0(j)|$ for vertices $j$ and times $t\in[0,T]$. Writing $H_Q\phi_k=\lambda_k\phi_k$ with orthonormal real eigenvectors $\phi_k$, the sufficient criterion uses Condition $B_2$: all values $\lambda_j+\lambda_k$ for $1\leq j\leq k\leq n$ are pairwise distinct; Condition $(I)$: the squared-eigenvector matrix $S\in\mathbb{R}^{n\times n}$ with $S_{jk}=\phi_k(j)^2$ is invertible; and Condition $(S)$: every pair of eigenvectors has a vertex $j$ with $\phi_k(j)\phi_\ell(j)\neq0$. Under these conditions, the magnitude trajectory determines $u_0$ up to multiplication by a scalar $e^{i\theta}$.
Implementation notes
1. Integration point: replace the hidden-state transition in a small complex-valued RNN or graph neural network with a fixed graph Schrödinger encoder. Let $h\in\mathbb{C}^n$ be the incoming hidden state, choose a connected graph Laplacian $\Delta_G$, and compute $K$ intensity observations $z_r=|U_rh|$ with $U_r=\exp(-i\tau_rH_Q)$ and $\tau_r=r\Delta t$. Concatenate $z_0,\ldots,z_{K-1}$ and feed them to a real MLP or classifier. 2. Pseudocode: initialize $H=\Delta_G+\operatorname{diag}(q)$; precompute $U=\operatorname{matrix\_exponential}(-i\Delta tH)$; set $x=h$; for $r=0,\ldots,K-1$, append $|x|$ to a list and set $x=Ux$; concatenate the list and apply the readout. Use $K=4$--$16$ and normalize each magnitude vector by its RMS. 3. During initialization, eigendecompose $H=\Phi\Lambda\Phi^T$, calculate $S=\Phi\odot\Phi$, reject potentials with small $\sigma_{\min}(S)$, reject near-collisions among all pair sums $\lambda_j+\lambda_k$, and reject eigenvector pairs with no overlapping support. In floating point, use thresholds such as $\sigma_{\min}(S)>10^{-2}$ and minimum pair-sum separation above $10^{-3}$ after spectral normalization. The uniqueness theorem is imported from the paper; the time grid and conditioning thresholds are empirical. 4. First experiment: compare a 2-layer complex graph-RNN using this encoder against a same-parameter real GRU and a generic magnitude-only encoder on synthetic graph-sequence classification. Generate random complex initial states and labels invariant to global phase, train at equal parameter count and FLOPs, and test robustness to random global phase and additive sensor noise. Success is higher accuracy and lower variance across global phases, ideally matching the full complex-state baseline while exposing only magnitude observations.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Confirmed ✓
Stage 2 · Practical benchmark result: Significantly worse ✗
Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works
Stage 1 — Mechanism check agent confidence 8/10
Built a fixed complex graph Schrödinger encoder that exposes only coordinate magnitudes over multiple times, with spectral diagnostics and a reproducible toy classification benchmark. For n=6, the Jacobian rank rose from 6 at K=1 to the predicted 11=2n−1 at K=2 and remained 11 thereafter; global-phase perturbation changed observations by only 6.7e-16. Increasing random potential scale from 0.02 to 1.0 increased median pair-sum separation from 9.6e-5 to 3.1e-2 and sigma_min(S) from 0.0026 to 0.106, while multi-time classification reached 0.878 versus 0.603 for a single-time magnitude baseline. The mechanism manifested, but this does not establish superiority over a trained GRU or full recurrent architecture.
- Agent confidence
- 8/10
- Baseline
- Single-time coordinate magnitude with standardized logistic readout: clean accuracy 0.603
- Idea
- Six-time fixed spectral phaseless encoder with standardized logistic readout: clean accuracy 0.878; noisy test accuracy 0.858; local Jacobian rank 11/11 after K>=2; global-phase max error 6.66e-16; selected operator sigma_min(S)=0.115 and pair-sum gap=0.00945
Limitations: The experiment used a path graph and a linear logistic readout, not a trained complex RNN, GRU, or graph-neural sequence model. It checked local differential identifiability rather than proving or numerically exhaustively testing global injectivity, and did not perform FLOP-matched recurrent training or broad graph/noise/time-grid sweeps.
How to run: python3 spectral_phaseless_experiment.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2607.26705)
Significantly WORSE than the tuned baseline (+4490.7%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0004
- Idea mean
- 0.0181
- Effect (Δ)
- 0.0177 (+4490.7%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±1565.8%
- Mechanism
- Confirmed ✓
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
The idea model has more parameters than the GRU baseline (24961 vs 13313). It mean-pools the input window before spectral encoding rather than evolving a spectral hidden state at every input step. Larger datasets, alternative graphs, sensor noise, longer training, and a K sweep were not tested.
How to run: python3 stage2_bench.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Stage 2 — Benchmark latest Failed ✗
Agent confidence: 10/10
The spectrally identifiable phaseless recurrent layer was implemented and evaluated on the registered dynamics track. Its phase-invariance mechanism was confirmed, but it was significantly worse on task MSE than the tuned GRU baseline: 0.0180801 versus 0.000393844, paired delta +0.0176862 with permutation p=0.0081.
- Baseline
- Mean test MSE 0.0003938444 over 8 seeds; best learning rate 0.01.
- Idea
- Mean test MSE 0.01808006 over 8 seeds; best tested learning rate 0.01; 0/8 paired wins.
Benchmark result
Significantly WORSE than the tuned baseline (+4490.7%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0004
- Idea mean
- 0.0181
- Effect (Δ)
- 0.0177 (+4490.7%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±1565.8%
- Mechanism
- Confirmed ✓
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
The idea model has more parameters than the GRU baseline (24961 vs 13313). It mean-pools the input window before spectral encoding rather than evolving a spectral hidden state at every input step. Larger datasets, alternative graphs, sensor noise, longer training, and a K sweep were not tested.
How to run: python3 stage2_bench.py
Artifacts
- bench_report.json 3.5 KB View Raw
- prediction_report.json 2.1 KB View Raw
- report.md 1.8 KB View
- report_bench_2026-08-31T182012.md 5.0 KB View
- results.json 3.6 KB View Raw
- spectral_phaseless_experiment.py 7.1 KB View Raw
- stage2_bench.py 4.0 KB View Raw