The Dynamics of Discoverability: How Trajectories and Priors Shape Equation Recovery
Not on the map. Why?
Abstract (arXiv)
How the dynamical regime of the observed system affects equation discovery has mainly been investigated through comparisons across systems. However, such comparisons vary both the equations and the dynamics, confounding the effect of the regime with the difficulty of recovering the equations symbolically. We separate the two by varying the forcing of Lorenz-84, moving its fully observed post-transient trajectories through fixed-point, periodic, and chaotic regimes while preserving the equations' functional form. Within each regime, we separately vary the amount of data, the noise, and prior knowledge of which terms the equations contain. We then measure how well two complementary approaches, sparse regression over a fixed library of candidate functions (SINDy) and an evolutionary search over symbolic expression trees (PySR), recover the true equations' terms and coefficients. We find that recovery depends on whether the sampled states distinguish combinations of candidate functions: equations remain poorly recovered from fixed-point data even when the candidate set contains only the true terms. We link the effect of the dynamics on both algorithms to one object: the moment matrix of the candidate functions under the invariant measure of the regime. Small eigenvalues mark weakly distinguishable combinations of candidate functions: we show that more data, less noise, and more prior knowledge can mitigate the resulting recovery difficulties, while a zero eigenvalue makes distinct equations indistinguishable on the visited states. Hence, a more precise prior needs less informative data, with consequences for data collection, method design, and evaluation.
Key facts. Cited by 1 paper in the map (data through 7 Oct 2026). Builds on 16 of its 35 references that are in the map. First posted on arXiv 20 Jul 2026.
Builds on
16 of 35 references.
References · 16
- 2021Contemporary Symbolic Regression Methods and their Relative Performancecited by 209
- 2020SINDy-PI: A Robust Algorithm for Parallel Implicit Sparse Identification of Nonlinear Dynamicscited by 125
- 2024LLM-SR: Scientific Equation Discovery via Programming with Large Language Modelscited by 117
- 2021PySINDy: A comprehensive Python package for robust sparse system identificationcited by 88
- 2021Chaos as an interpretable benchmark for forecasting and data-driven modellingcited by 87
- 2025LLM-SRBench: A New Benchmark for Scientific Equation Discovery with Large Language Modelscited by 61
- 2023ODEFormer: Symbolic Regression of Dynamical Systems with Transformerscited by 58
- 2017Extracting Sparse High-Dimensional Dynamics from Limited Datacited by 50
- 2022Rethinking Symbolic Regression Datasets and Benchmarks for Scientific Discoverycited by 44
- 2020Sparse Identification of Nonlinear Dynamical Systems via Reweighted ℓ1-regularized Least Squarescited by 21
- 2023Benchmarking sparse system identification with low-dimensional chaoscited by 19
- 2021A toolkit for data-driven discovery of governing equations in high-noise regimescited by 15
4 more
- 2024Benchmarking symbolic regression constant optimization schemescited by 6
- 2025When is a System Discoverable from Data? Discovery Requires Chaoscited by 6
- 2024Invariant Measures in Time-Delay Coordinates for Unique Dynamical System Identificationcited by 5
- 2025Information theory and discriminative sampling for model discoverycited by 3
Built on by
Citing papers per year: 2026: 1.
1 paper.