Adaptive Control of Scientific Instruments with Nonstationary Responses through Multimodal Agents
Organizations: Beijing ENN Fusion Energy Science and Technology Co., Ltd., Beijing 101111, China · Beijing Key Laboratory of High Magnetic Field Spherical Torus Fusion Energy, Beijing 101111, China · Hebei Key Laboratory of Compact Fusion, Langfang 065000, China
Abstract
Controlling scientific instruments requires parameters to be adjusted progressively in response to feedback from successive operations. When unobserved operating conditions evolve with operational history, responses become nonstationary, and a static mapping from actions to responses becomes inadequate for decisions across shots. Shot-to-shot tuning of neutral beam injection (NBI) ion sources in magnetic confinement fusion devices is representative of this class of problems: the conditions of electrode surfaces and vacuum are difficult to observe directly, discharge outcomes vary between shots, and experiments are costly and carry operational risks. We propose a physically grounded multimodal agent with three independently auditable components. Deterministic algorithms extract symbolic physical quantities and visual evidence from multichannel waveforms; a memory with two layers organizes trajectories from the current task, experience from earlier tasks, and procedural knowledge, from which a multimodal large language model (MLLM) generates candidate actions; and an independent action feasibility check enforces action constraints before execution. In a strictly chronological replay of 455 shots from four EXL-50U platforms, physical symbols alone outperformed waveform images alone by 0.21-0.32 in macro-F1 for quality classification across all four evaluated models. Retrieved context increased the mean agreement with historical action directions by 22.4-42.7 percentage points, whereas the contribution of visual evidence varied across models. Under supervised operation on five ion sources, the system achieved a beam-extraction rate of 86.80% across 409 executed shots. These results demonstrate the engineering feasibility of the framework for supervised sequential tuning on a physical scientific instrument.
Figures & tables
| Method | Handles latent- state drift | Requires a high- fidelity simulator | Uses unstructured knowledge | Local–global feature integration | Action-constraint mechanism | Physical-system closed loop |
| Ion-source/NBI static surrogate models and GA/PSO [ 20 , 22 , 23 , 24 , 25 ] | Usually no | No | No | Scalar features; no explicit integration | Static bounds and constraints can be encoded | No reported NBI tuning closed loop |
| Simulator-driven plasma RL [ 17 ] | Indirectly, through state observations | Yes | No | Local magnetic measurements and global shape objectives | Action bounds and device interlocks | Yes (plasma control) |
| Standard/safe Bayesian optimization [ 27 , 28 ] | Not addressed for latent-state drift | No | No | Point observations and a global surrogate | Conditional guarantees in safe variants | Applicable to online experiments |
| LLM agents for scientific experiments [ 31 , 32 ] | Through context or memory | No | Yes | Task-relevant inputs; no explicit integration | Can use semantic guardrails or external validators [ 38 ] | Yes (chemical experiments) |
| This work | History-conditioned dual-layer memory | No | Yes | Quantitative details and overall waveform morphology | Deterministic action feasibility check | Supervised shot-to-shot execution |
| Class | Number of shots |
|---|---|
| Full breakdown | 67 |
| Insufficient pulse width | 113 |
| Long-pulse high current | 80 |
| Long-pulse low current | 81 |
| High-frequency oscillation | 51 |
| Fault noise | 63 |
| Configuration | Symbolic grounding | Visual evidence | Experience and knowledge retrieval | Feasibility check | Human feedback |
|---|---|---|---|---|---|
| Baseline | |||||
| Physical-only | |||||
| +Physical | |||||
| +RAG | |||||
| Full |
| Platform | Shots | Current MAE/A | Main-width MAE/ms | Flat-top-count accuracy | Flat-top event- count F 1 | Breakdown-count accuracy | Breakdown event- count F 1 |
|---|---|---|---|---|---|---|---|
| Test Bench 1 | 127 | 0.16 | 4.46 | 100.00% | 100.00% | 100.00% | 100.00% |
| Test Bench 2 | 116 | 0.64 | 14.35 | 95.69% | 96.08% | 100.00% | 100.00% |
| Experimental Bench 1 | 110 | 0.70 | 13.24 | 98.18% | 97.44% | 99.09% | 97.56% |
| Experimental Bench 2 | 102 | 1.34 | 10.35 | 94.12% | 96.20% | 100.00% | 100.00% |
| Pooled | 455 | 0.68 | 10.42 | 97.14% | 97.61% | 99.78% | 99.66% |
| Model | Configuration | Quality-classification macro-F 1 | Current MAE/A | Pulse-width MAE/ms | Flat-top-count MAE | Mean directional agreement |
|---|---|---|---|---|---|---|
| Qwen3.5-Plus | Baseline | 0.47 | 3.67 | 138.19 | 0.27 | 36.7% |
| Qwen3.5-Plus | Physical-only | 0.79 | 1.40 | 28.37 | 0.09 | 54.4% |
| Qwen3.5-Plus | +Physical | 0.73 | 1.54 | 42.23 | 0.09 | 55.8% |
| Qwen3.5-Plus | +RAG | 0.71 | 1.24 | 19.36 | 0.09 | 78.2% |
| Qwen3.5-Plus | Full | 0.68 | 1.31 | 24.33 | 0.08 | 86.0% |
| MiniMax3 | Baseline | 0.57 | 3.91 | 107.80 | 0.24 | 35.6% |
| Ion source | Executed shots | With beam extraction | Without beam extraction | Beam-extraction rate |
|---|---|---|---|---|
| Test Bench 2 | 49 | 45 | 4 | 91.84% |
| Source 1, Experimental Bench 2 | 64 | 53 | 11 | 82.81% |
| Source 2, Experimental Bench 2 | 83 | 69 | 14 | 83.13% |
| Source 3, Experimental Bench 2 | 103 | 91 | 12 | 88.35% |
| Source 4, Experimental Bench 2 | 110 | 97 | 13 | 88.18% |
| All | 409 | 355 | 54 | 86.80% |
| Metric | Trial 1 | Trial 2 | Trial 3 | Total |
|---|---|---|---|---|
| Shot range | 11394–11402 | 11884–11891 | 11948–11953 | — |
| Number of shots | 9 | 8 | 6 | 23 |
| Shots with a valid beam interval | 8 | 7 | 6 | 21 |
| Valid-beam-interval rate | 88.89% | 87.50% | 100.00% | 91.30% |
| Median acceleration-voltage range/kV | 26.39–33.36 | 21.06–44.54 | 22.43–38.44 | 21.06–44.54 |
| Median acceleration-current range/A | 16.60–26.88 | 10.31–23.46 | 10.89–28.56 | 10.31–28.56 |
Appendix figures & tables6 assets
Supplementary material from the paper’s appendix.
Appendix
| Model | Configuration | Current/A | Pulse width/ms | Flat-top count | Quality classification | ||||
|---|---|---|---|---|---|---|---|---|---|
| MAE | RMSE | MAE | RMSE | MAE | RMSE | Accuracy | F 1 | ||
| Qwen3.5-Plus | Baseline | 3.67 | 14.41 | 138.19 | 328.19 | 0.27 | 0.66 | 0.51 | 0.47 |
| Qwen3.5-Plus | Physical-only | 1.40 | 5.13 | 28.37 | 86.58 | 0.09 | 0.41 | 0.79 | 0.79 |
| Qwen3.5-Plus | +Physical | 1.54 | 7.08 | 42.23 | 121.89 | 0.09 | 0.34 | 0.73 | 0.73 |
| Qwen3.5-Plus | +RAG | 1.24 | 4.92 | 19.36 | 67.51 | 0.09 | 0.37 | 0.74 | 0.71 |
| Qwen3.5-Plus | Full | 1.31 | 5.05 | 24.33 | 77.34 | 0.08 | 0.36 | 0.84 | 0.68 |
| Model | Configuration | Mean | |||||
|---|---|---|---|---|---|---|---|
| Qwen3.5-Plus | Baseline | 55.6% | 47.9% | 26.6% | 18.7% | 34.9% | 36.7% |
| Qwen3.5-Plus | Physical-only | 66.6% | 69.2% | 59.6% | 39.8% | 36.9% | 54.4% |
| Qwen3.5-Plus | +Physical | 66.6% | 76.9% | 56.9% | 40.9% | 37.8% | 55.8% |
| Qwen3.5-Plus | +RAG | 87.9% | 82.2% | 72.3% | 75.4% | 73.2% | 78.2% |
| Qwen3.5-Plus | Full | 89.7% | 84.4% | 85.1% | 81.3% | 89.7% | 86.0% |
| MiniMax3 | Baseline | 49.2% | 48.8% | 23.7% | 19.1% | 37.4% | 35.6% |
| Ion source | Manual tuning | Agent tuning | Difference (percentage points) | ||
|---|---|---|---|---|---|
| Beam extracted/ executed | Extraction rate | Beam extracted/ executed | Extraction rate | ||
| ID 5 (Source 1) | 377/436 | 86.47% | 53/64 | 82.81% | |
| ID 6 (Source 2) | 337/419 | 80.43% | 69/83 | 83.13% | |
| ID 7 (Source 3) | 361/419 | 86.16% | 91/103 | 88.35% | |
| ID 8 (Source 4) | 346/421 | 82.19% | 97/110 | 88.18% | |
| Total | 1421/1695 | 83.83% | 310/360 | 86.11% | |
| Symbol | Physical quantity | Unit | PLC permissible range | Maximum per-shot change |
|---|---|---|---|---|
| Acceleration voltage | kV | – | Adjustment step: 1–5 kV | |
| Filament voltage | V | – | Adjustment step: 0.1–0.5 V | |
| Arc voltage | V | – | Adjustment step: 1–5 V | |
| Gas flow | sccm | – | Adjustment step: 10–50 sccm | |
| Suppression voltage | V | – | Constrained during the check to 4%–10% of |
| Symbol | Meaning | Unit | Value and determination method |
|---|---|---|---|
| Median-filter window width | samples | ||
| Measurement-spike rejection bound | dimensionless | ; is compared with the global standard deviation, rather than a local baseline | |
| Fraction trimmed from each end for the trimmed mean | % | ||
| Sliding-window width for the trimmed mean | samples | ||
| Absolute-amplitude threshold for candidate pulses | A | ; the acceleration-current threshold is identical across all four platforms | |
| Relative coefficient for candidate pulses | dimensionless | ; relative to the maximum pulse mean |
Adaptive Control
Adaptive control lets systems update models, policies, or controller parameters online to maintain performance as dynamics and disturbances change.
See all 148 papers