Dynamic Bayesian Optimization Framework for Instruction Tuning in Partial Differential Equation Discovery
Junqi Qu, Yan Zhang, Shangqian Gao, Shibo Li
TL;DR
The paper addresses instruction brittleness in LLM-based equation discovery by treating prompt engineering as a dynamic, discrete optimization problem. NeuroSym-BO uses a closed-loop three-agent system where a Bayesian Optimization agent selects among a bank of 100 reasoning strategies to assemble adaptive prompts for the LLM, guided by a numerical evaluator that fits coefficients and balances accuracy with parsimony. Empirically, the approach yields higher recovery rates and more parsimonious symbolic forms across five PDE benchmarks (e.g., Fisher approaching $R^2=0.9999$) compared to a fixed-prompt baseline, demonstrating improved sample efficiency and robustness. This framework significantly advances automated scientific discovery by integrating adaptive prompt design with rigorous numerical feedback, with potential extensions to higher-dimensional or more complex dynamical systems.
Abstract
Large Language Models (LLMs) show promise for equation discovery, yet their outputs are highly sensitive to prompt phrasing, a phenomenon we term instruction brittleness. Static prompts cannot adapt to the evolving state of a multi-step generation process, causing models to plateau at suboptimal solutions. To address this, we propose NeuroSymBO, which reframes prompt engineering as a sequential decision problem. Our method maintains a discrete library of reasoning strategies and uses Bayesian Optimization to select the optimal instruction at each step based on numerical feedback. Experiments on PDE discovery benchmarks show that adaptive instruction selection significantly outperforms fixed prompts, achieving higher recovery rates with more parsimonious solutions.
