LLM-ODE: Data-driven Discovery of Dynamical Systems with Large Language Models
About
Discovering the governing equations of dynamical systems is a central problem across many scientific disciplines. As experimental data become increasingly available, automated equation discovery methods offer a promising data-driven approach to accelerate scientific discovery. Among these methods, genetic programming (GP) has been widely adopted due to its flexibility and interpretability. However, GP-based approaches often suffer from inefficient exploration of the symbolic search space, leading to slow convergence and suboptimal solutions. To address these limitations, we propose LLM-ODE, a large language model-aided model discovery framework that guides symbolic evolution using patterns extracted from elite candidate equations. By leveraging the generative prior of large language models, LLM-ODE produces more informed search trajectories while preserving the exploratory strengths of evolutionary algorithms. Empirical results on 91 dynamical systems show that LLM-ODE variants consistently outperform classical GP methods in terms of search efficiency and Pareto-front quality. Overall, our results demonstrate that LLM-ODE improves both efficiency and accuracy over traditional GP-based discovery and offers greater scalability to higher-dimensional systems compared to linear and Transformer-only model discovery methods.
Related benchmarks
| Task | Dataset | Result | Rank | |
|---|---|---|---|---|
| Symbolic Discovery | ODEBench Generalization 1.0 | NMSE0.0047 | 12 | |
| Symbolic Discovery | ODEBench Out-of-Distribution 1.0 | NMSE0.0243 | 12 | |
| Symbolic Discovery | ODEBench Reconstruction 1.0 | NMSE4.12e-5 | 12 | |
| Symbolic Discovery | ODEBase | Reconstruction NMSE7.46e-5 | 12 | |
| Symbolic Discovery | ODEBench 1.0 | Complexity25.2 | 12 | |
| Dynamical System Discovery | Dynamical Systems D=1, n=23 | Success Rate (Tol 10^-1)21 | 6 | |
| Dynamical System Discovery | Dynamical Systems D=2, n=28 | Success Rate (Tol 10^-1)18 | 6 | |
| Dynamical System Discovery | Dynamical Systems D=4, n=18 | Performance at Tolerance 1e-16 | 6 | |
| Dynamical System Discovery | Dynamical Systems D=3, n=22 | Performance ($10^{-1}$ Tolerance)8 | 6 |