REVIEW 3 major objections 4 minor 15 references
A closed-loop evolutionary algorithm that uses measured training outcomes to guide an LLM reduced PINN solution error by 95.38% on a 1D wave equation.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-01 22:56 UTC pith:BIAHTLJV
load-bearing objection A transparent feasibility study of EA-guided LLM search for PINN architecture, with the expected limitation: no matched baselines, so the feedback loop's causal role is unproven. the 3 major comments →
Evolutionary Algorithm-Guided LLMs for Physics-Informed Neural Network Design
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that an LLM, when embedded as a conditional mutation/crossover operator inside a population-based evolutionary loop, can accumulate empirical experience and discover PINN configurations that outperform its own initial proposals. The discovery rests on a lineage-resolved experiment: the best final configuration was assembled from separate validated changes (residual connections and depth five) combined via cross-generation memory, then refined through targeted mutations of width and collocation density. The paper also reports that the lowest-MSE model had a PDE residual of 8.12 against 0.109 for a less accurate model, showing that MSE-only selection does not guarantee phy
What carries the argument
The central object is the AlgorithmSpec, a structured declaration of the full PINN algorithm (network, loss, sampling, optimizer, constraints). The architecture maintains an evaluated population, computes effective fingerprints that exclude metadata so that renames don't masquerade as new designs, selects elite and diverse parents, applies parent-conditioned mutation and distinct-parent crossover, and writes a cross-generation memory of measured successes and failures into the next LLM prompt. Every proposed configuration is executed directly under an exact optimizer-step budget with no silent fallback, making each proposal auditable.
Load-bearing premise
The paper assumes the measured improvement comes from the evolutionary feedback loop, but it never tests an LLM without feedback or random search under matched compute; if the LLM's prior knowledge or chance alone produced the 95.38% result, the central closed-loop claim would fail.
What would settle it
Run the same Wave-C search with an otherwise identical LLM that receives no evaluated feedback (e.g., fresh prompts with no lineage or memory) and a random-search baseline over the same AlgorithmSpec space, both under the same 600,000-optimizer-step budget. If either achieves a comparable or larger MSE reduction, the claim that the evolutionary feedback loop is responsible for the improvement would be refuted. Alternatively, a factorial ablation that disables the cross-generation memory should show a significant drop in the best-MSE reduction if memory is load-bearing.
If this is right
- If the approach transfers, LLM-guided search can be made to accumulate measured experience for PINN design rather than relying on one-shot proposals.
- Structured specifications and effective fingerprinting allow recombination of validated sub-designs across generations.
- Exact-budget execution without fallback makes each LLM hypothesis directly comparable and auditable.
- The observed objective mismatch motivates physics-aware multi-objective evolution, e.g., constraining PDE residual or maintaining a Pareto front.
- The lineage analysis offers a template for attributing improvements to specific design changes.
Where Pith is reading between the lines
- The mechanism may extend to other PDE families and inverse problems, but the paper's single-equation, two-seed scope leaves transfer unshown; testing on multiple benchmarks with more seeds would be a natural next step.
- The objective mismatch suggests a concrete repair: evolve over a Pareto front of MSE and residual, or add a physics-residual constraint to the selection step; the paper explicitly points toward this.
- The deduplication of effective specifications is a general idea: any LLM search over structured algorithms can benefit from fingerprinting that ignores cosmetic metadata.
- The wide gap between seeds (2.97% vs 95.38% reduction) implies high sensitivity to initial population and stochasticity; users of such a system should run multiple restarts or seed diversity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces an evolutionary architecture that guides an LLM to propose complete PINN configurations (AlgorithmSpecs) across generations, using measured training outcomes to inform parent selection, mutation/crossover, fingerprint deduplication, elite survival, and cross-generation memory. All proposed configurations are executed with an exact optimizer-step budget. On a one-dimensional multiscale wave equation, two independent ten-generation runs are reported; the best configuration appears in the final generation, with MSE reductions of 2.97% and 95.38% relative to the best initial-generation proposal. The paper also reports an observed mismatch between low solution MSE and high PDE residual, and it explicitly acknowledges limitations in experimental breadth and missing baselines.
Significance. If the closed-loop attribution holds, the paper would be a useful demonstration that LLM-driven structured search can improve PINN configurations beyond one-shot proposals. The strengths are real: the execution contract is auditable (no generated Python, exact step budgets, no silent fallback), the artifacts are thoroughly documented, and the paper is unusually transparent about its limitations, including an implementation issue in the strict-budget audit. The lineage-resolved analysis in Table 2 is compelling as a case study. However, the current evidence supports only a narrow feasibility claim, not the broader causal claim that the evolutionary feedback loop, rather than LLM prior knowledge or unguided sampling, produced the reported improvement.
major comments (3)
- [§6, §5.2] The central causal claim is not established. The abstract and contributions state that the stronger run's 95.38% reduction was 'achieved through a sequence of interpretable evolutionary changes' and that feedback drove improvement, but Section 6 explicitly concedes there is no comparison against Bayesian optimization, random search, PBT, OPRO, or PINNsAgent under matched compute, and no ablation removing reflection memory. The lineage in Table 2 is observational. Without a matched no-feedback control (e.g., the same LLM sampler with per-generation context reset, or random search over the same AlgorithmSpec space under the same 10,000-step budget), the improvement could stem from LLM pretrained PINN knowledge or from sampling variance. This is load-bearing for the paper's main contribution, so the causal wording must be softened or the missing control must be supplied.
- [§4.2, Appendix C, Table 1] The statistical basis is too thin for the headline claim. Two seeds give no variance estimate, and the '95.38%' result is a single favorable run; the other seed improves only 2.97%. The paper acknowledges this, but the abstract and contributions still lead with the stronger run. Additionally, the advertised '60 successful PINNs' is not accurate as stated: Appendix C reports that seed 1's online search produced 28 successful models, and two post-hoc control configurations were later trained and assigned to missing generations to complete the 30-model accounting. The abstract and Section 4.2 should either explicitly state that the 60-PINN total includes two audit controls, or report the 58-model total from the actual evolutionary search. The controls do not change the best-so-far trajectory, but the current wording is misleading.
- [§5.4] The objective-mismatch observation is interesting but its interpretation needs a concrete check. The claim that 'low solution error can coexist with a high PDE residual' rests on comparing MSE-optimal models across two runs. It would be more convincing if the PDE residual were reported on the same evaluation grid or with the same sampling density for both models, and if the residual magnitude were compared against a baseline scale (e.g., the residual of the analytic solution interpolated by a standard PINN). As written, the residual value 8.12 is an absolute number whose significance depends on the residual definition and normalization. This is not fatal, but it is a load-bearing supporting observation for the proposed multi-objective extension.
minor comments (4)
- [Table 1] The 'Mean' row is unclear: it averages initial and final MSE values across runs, but the 'Reduction' column shows a value with a negative sign and no clear meaning. Either remove the row or explain the aggregation (e.g., geometric mean of reductions).
- [Figure 2 caption] The caption mentions that the candidate scatter includes two post-audit contract controls, but the main text does not identify which points correspond to those controls. Please mark them explicitly so readers can exclude them when assessing the evolutionary trajectory.
- [§3.3, Eq. (4)] The pre-scoring function q(a) uses hand-chosen weights and the LLM's self-reported predicted improvement Δ̂. This is a heuristic that could introduce selection bias; the paper should state explicitly that these weights are not tuned and that q(a) has no verified correlation with post-training quality, as the axiom ledger in related discussions suggests.
- [§4.2] The search protocol says 'generation zero creates unconstrained proposals' but does not specify how the initial population of 10 is formed (e.g., all from the same prompt, or with varying contexts). Please clarify, because the initial-population baseline is used for the reduction percentages.
Circularity Check
No significant circularity: the paper reports empirical evolutionary search results with an external MSE fitness; limitations concern causal attribution, not construction.
full rationale
This paper contains no first-principles derivation whose conclusion is equivalent to its inputs. The core object J(a) (Eq. 2) is an externally evaluated MSE on a fixed grid, and the reported improvements (2.97% and 95.38%) are measured training outcomes, not quantities reconstructed from the search inputs. The pre-training score q(a) (Eq. 4) is a hand-specified selection heuristic with fixed weights; it is not fitted to the reported final MSE, so the final result is not forced by construction. The lineage narrative is supported by recorded parent identifiers and cross-generation context rather than by a mathematical identity. The paper itself states in Section 6 that it does not compare against Bayesian optimization, random search, PBT, OPRO, or PINNsAgent under matched compute and that the operator analysis is observational, not a factorial ablation; this is an honest limitation about causal attribution and experimental breadth, not evidence of circularity. Appendix C discloses that two post-hoc control configurations were trained to complete accounting, and explicitly states they did not change any best-so-far value or lineage, so they do not circularly define the reported result. There are no load-bearing self-citations and no known result is renamed or repackaged as a new prediction. The only substantive caveat is that the central claim that the feedback loop caused the improvement is not isolated from LLM prior knowledge or unguided search; that is a correctness/experimental-design concern, not a circularity concern, and the paper appropriately bounds its claims as feasibility evidence.
Axiom & Free-Parameter Ledger
free parameters (4)
- q(a) scoring weights =
w_Δ=0.40, w_ν=0.30, w_τ=0.15, w_e=0.10, w_r=0.05
- Elite preservation size k
- Mutation change count (1–3 fields)
- Population/generation protocol =
10 per gen, 3 successful, 10,000 steps
axioms (5)
- domain assumption DeepSeek-chat reliably produces syntactically valid AlgorithmSpec objects and follows the structured prompt.
- domain assumption The analytic solution of Wave-C (Eq. 9) is the correct ground truth for MSE; this benchmark is representative enough to draw conclusions about PINN design.
- domain assumption Canonicalizing only effective fields via SHA256 (Eq. 3) preserves algorithmic identity; metadata like rationale and lineage can be discarded.
- ad hoc to paper The candidate pre-scoring function q(a) correlates with post-training quality.
- domain assumption 10,000 optimizer steps are sufficient for MSE to be a stable indicator of configuration quality.
read the original abstract
Physics-informed neural networks (PINNs) are unusually sensitive to interacting choices of architecture, activation, loss weighting, collocation, optimization, and constraint enforcement. Large language models (LLMs) can propose these choices, but independent recommendations do not accumulate experience from previously trained PINNs. We propose a closed-loop evolutionary algorithm that guides an LLM to generate complete, executable PINN configurations across generations, using measured training outcomes to determine subsequent search decisions. The algorithm maintains an evaluated population and lineage, applies parent-conditioned mutation and crossover, preserves elite and diverse solutions, rejects effective duplicates, and converts parent-relative successes and failures into the next-generation context supplied to the LLM. Every proposed configuration is executed directly under an exact optimizer-step budget. On a one-dimensional multiscale wave equation, two independent ten-generation runs trained 60 PINNs for 600,000 optimizer steps. In both runs, the best configuration appeared in the final generation, with best mean-squared error reduced by 2.97\% and 95.38\% relative to the initial population. The stronger run validated residual connections and increased depth on separate branches, combined them in a later generation, and then refined width and collocation density. It also revealed that low solution error can coexist with a high PDE residual. These results demonstrate the feasibility of evolutionary-algorithm-guided LLMs for PINN design on a controlled PDE while motivating broader, physics-aware evaluation.
Figures
Reference graph
Works this paper leans on
-
[1]
Journal of Computational Physics , volume=
Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations , author=. Journal of Computational Physics , volume=. 2019 , doi=
2019
-
[2]
Nature Reviews Physics , volume=
Physics-informed machine learning , author=. Nature Reviews Physics , volume=. 2021 , doi=
2021
-
[3]
SIAM Journal on Scientific Computing , volume=
Understanding and mitigating gradient flow pathologies in physics-informed neural networks , author=. SIAM Journal on Scientific Computing , volume=. 2021 , doi=
2021
-
[4]
When and why
Wang, Sifan and Yu, Xinling and Perdikaris, Paris , journal=. When and why. 2022 , doi=
2022
-
[5]
Advances in Neural Information Processing Systems , volume=
Characterizing possible failure modes in physics-informed neural networks , author=. Advances in Neural Information Processing Systems , volume=. 2021 , url=
2021
-
[6]
Computer Methods in Applied Mechanics and Engineering , volume=
A comprehensive study of non-adaptive and residual-based adaptive sampling for physics-informed neural networks , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2023 , doi=
2023
-
[7]
2021 , doi=
Lu, Lu and Meng, Xuhui and Mao, Zhiping and Karniadakis, George Em , journal=. 2021 , doi=
2021
-
[8]
Computer Methods in Applied Mechanics and Engineering , volume=
Exact imposition of boundary conditions with distance functions in physics-informed deep neural networks , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2022 , doi=
2022
-
[9]
Advances in Neural Information Processing Systems , volume=
Implicit neural representations with periodic activation functions , author=. Advances in Neural Information Processing Systems , volume=. 2020 , url=
2020
-
[10]
2024 , doi=
Hao, Zhongkai and Yao, Jiachen and Su, Chang and Su, Hang and Wang, Ziao and Lu, Fanzhi and Xia, Zeyu and Zhang, Yichi and Liu, Songming and Lu, Lu and Zhu, Jun , booktitle=. 2024 , doi=
2024
-
[11]
2025 , url=
Wuwu, Qingpo and Gao, Chonghan and Chen, Tianyu and Huang, Yihang and Zhang, Yuekai and Wang, Jianing and Li, Jianxin and Zhou, Haoyi and Zhang, Shanghang , booktitle=. 2025 , url=
2025
-
[12]
arXiv preprint arXiv:2401.02051 , year=
Evolution of Heuristics: Towards efficient automatic algorithm design using large language model , author=. arXiv preprint arXiv:2401.02051 , year=
-
[13]
International Conference on Learning Representations , year=
Large language models as optimizers , author=. International Conference on Learning Representations , year=
-
[14]
Nature , volume=
Mathematical discoveries from program search with large language models , author=. Nature , volume=. 2024 , doi=
2024
-
[15]
IEEE Transactions on Neural Networks , volume=
Artificial neural networks for solving ordinary and partial differential equations , author=. IEEE Transactions on Neural Networks , volume=. 1998 , doi=
1998
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.