Grounding LLMs in Scientific Discovery via Embodied Actions
Organizations: Tsinghua University, Beijing, China.
Abstract
Large Language Models (LLMs) have shown significant potential in scientific discovery but struggle to bridge the gap between theoretical reasoning and verifiable physical simulation. Existing solutions operate in a passive "execute-then-response" loop and thus lacks runtime perception, obscuring agents to transient anomalies (e.g., numerical instability or diverging oscillations). To address this limitation, we propose EmbodiedAct, a framework that transforms established scientific software into active embodied agents by grounding LLMs in embodied actions with a tight perception-execution loop. We instantiate EmbodiedAct within MATLAB and evaluate it on complex engineering design and scientific modeling tasks. Extensive experiments show that EmbodiedAct significantly outperforms existing baselines, achieving SOTA performance by ensuring satisfactory reliability and stability in long-horizon simulations and enhanced accuracy in scientific modeling.
Figures & tables
| Models | Overall | |||||||||
| Core/Ext | Core/Ext | Core/Ext | Core/Ext | Core/Ext | Core/Ext | Core/Ext | Core/Ext | Core/Ext | Core/Ext | |
| Generative Models | ||||||||||
| Qwen3-8B-VL-Instruct | 18.0 / 18.0 | 14.2 / 18.3 | 80.0 / 24.9 | 13.3 / 22.2 | 45.4 / 45.4 | 13.3 / 13.3 | 15.0 / 40.1 | – / 44.3 | 26.7 / 27.5 | 21.7 / 25.9 |
| Qwen3-235B-Instruct | 33.0 / 33.0 | 28.5 / 30.5 | 80.0 / 28.9 | 13.3 / 32.6 | 68.8 / 68.8 | 13.3 / 13.3 | 53.3 / 55.2 | – / 56.7 | 78.9 / 53.1 | 39.6 / 42.4 |
| GPT-5-mini | 38.3 / 38.3 | 56.5 / 36.3 | 80.0 / 36.4 | 20.0 / 26.7 | 25.0 / 25.0 | 40.0 / 40.0 | 76.7 / 82.6 | – / 85.2 | 81.1 / 41.0 | 42.0 / 46.9 |
| Claude-Sonnet-4.5 | 38.2 / 38.2 | 25.7 / 23.4 | 80.0 / 40.0 | 13.3 / 32.3 | 68.3 / 68.3 | 33.3 / 33.3 | 40.0 / 63.3 | – / 65.9 | 72.2 / 50.1 | 42.1 / 44.5 |
| Models | Chemistry | Physics | Math | Avg. | ||||||||||
| atkins | chemmc | quan | matter | avg. | fund | class | thermo | avg. | diff | stat | calc | avg. | ||
| Generative Models | ||||||||||||||
| Qwen3-8B-VL-Instruct | 12.50 | 25.00 | 12.50 | 0.00 | 12.50 | 12.50 | 9.52 | 10.00 | 10.25 | 6.25 | 10.00 | 20.00 | 11.11 | 11.21 |
| Gemini-3-Flash | 12.50 | 25.00 | 25.00 | 12.50 | 18.75 | 12.50 | 19.05 | 20.00 | 17.95 | 25.00 | 40.00 | 70.00 | 41.67 | 26.17 |
| Qwen3-235B-Instruct | 12.50 | 62.50 | 25.00 | 0.00 | 25.00 | 25.00 | 14.29 | 30.00 | 20.52 | 18.75 | 60.00 | 50.00 | 38.89 | 28.04 |
| GPT-5-mini | 25.00 | 50.00 | 50.00 | 25.00 | 37.50 | 37.50 | 19.05 | 50.00 | 30.77 | 31.25 | 70.00 | 80.00 | 55.56 | 41.12 |
| Models | GPT-5.2 | GPT-5-mini | ||
| EngDesign | SciBench-107 | EngDesign | SciBench-107 | |
| Generative Models | 48.0 | 44.9 | 42.0 | 41.1 |
| CodeAct | 55.4 | 49.5 | 52.9 | 48.6 |
| EmbodiedAct | 70.6 | 54.2 | 59.3 | 49.5 |
| w/o | 65.0 | 51.0 | 55.0 | 48.0 |
| w/o | 56.5 | 50.0 | 48.5 | 46.5 |
Appendix figures & tables14 assets
Supplementary material from the paper’s appendix.
Appendix
| Module | Output | Key Prompt Constraints |
| Hierarchical decomposition; toolbox-grounded steps; first-principles derivation; validation checkpoints | ||
| Executable primitives only; intermediate variable logging; ASCII-only syntax; scope-safe closures | ||
| N, E, W | Tri-state classification (Normal/Error/Warning); streaming observation parsing; Hot-Fix trigger | |
| Feedback t , | LLM-as-judge verification; local adjustment vs. global re-planning; physics-informed error analysis |
| Control Variable | CodeAct | EmbodiedAct |
| LLM backbone | Identical (GPT-5.2 / Claude-4.5 / Qwen3) | |
| Max interaction turns | 5 | 5 |
| Problem input | Identical statements, no ground-truth | |
| Task-specific tuning | None | None |
| MATLAB toolbox access | Equivalent (Python bridge / native) | |
| Differentiator | Single-prompt ReAct | Four-module architecture |
| Textbook | Domain | Original | Subset |
| class | Classical Mechanics | 56 | 21 |
| diff | Differential Equations | 50 | 16 |
| thermo | Thermodynamics | 66 | 10 |
| stat | Statistics | 72 | 10 |
| calculus | Calculus | 42 | 10 |
| quan | Quantum Mechanics | 33 | 10 |
| Subject | Direct | ODE/ | Root- | Conserv. | Numer. | Other | Total |
| Formula | IVP | Finding | Laws | Methods | |||
| Classical Mech. | 8 | 5 | 3 | 3 | 1 | 1 | 21 |
| Diff. Equations | 4 | 7 | 2 | 0 | 2 | 1 | 16 |
| Thermodynamics | 6 | 0 | 1 | 1 | 1 | 1 | 10 |
| Statistics | 7 | 0 | 0 | 0 | 2 | 1 | 10 |
| Calculus | 5 | 0 | 1 | 0 | 3 | 1 | 10 |
| Field | Type | Description |
| id | string | Unique message identifier (UUID) |
| type | enum | Message type (see Table 9 ) |
| payload | object | Type-specific data payload |
| timestamp | float | Unix timestamp of message creation |
| session_id | string | Session identifier for routing |
| operation_id | string | Links related messages in a flow |
| Category | Message Type | Direction |
| Operation Lifecycle | operation_request | |
| operation_ack | ||
| operation_start | ||
| operation_progress | ||
| operation_complete | ||
| operation_failed |
| Category | Case | Mechanism | Outcome |
| 1. Runtime Error Recovery | 1.1 | I/O mismatch Local Adjustment | ✓ |
| 1.2 | Undefined function Local Adjustment | ✓ | |
| 1.3 | Errors fixed, semantic extraction wrong | × | |
| 2. Intermediate Monitoring | 2.1 | 5% threshold verification | ✓ |
| 2.2 | Grid search + reflective refinement | ✓ | |
| 2.3 | Script structure errors persisted | × |
| Method | Total Tokens | Ratio | Avg Score | Pass Rate |
| Baseline (Direct) | 104,889 | 1.0 | 53.3 | 33.3% |
| CodeAct | 657,141 | 6.3 | 44.0 | 15.6% |
| EmbodiedAct | 1,115,740 | 10.6 | 61.6 | 51.1% |
| Module | Symbol | Input | Output | Total | % |
| Strategic Planner | 124,764 | 9,264 | 134,028 | 12% | |
| Primitive Generator | 601,858 | 44,778 | 646,636 | 58% | |
| Runtime Perception | 187,146 | 13,894 | 201,040 | 18% | |
| Reflective Decision | 123,772 | 10,264 | 134,036 | 12% | |
| Total | – | 1,037,540 | 78,200 | 1,115,740 | 100% |