When Does AI for PDEs Yield Scientific Evidence?
Organizations: School of Future Technology, South China University of Technology, China
Abstract
Existing AI-for-PDE benchmarks primarily assess models in terms of predictive or approximation accuracy. In physics research, however, AI outputs often serve as evidence for scientific claims. These two objectives are not equivalent: the former measures an output's agreement with a reference target or satisfaction of governing constraints; the latter asks whether, given a specified object of study, scientific claim, assumptions, and evidence standard, the output provides sufficient evidence for that claim. To bridge this gap, we extend a widely used PDE-simulation benchmark and a comprehensive benchmark for PDE inverse problems to enable, for the first time in AI for PDEs, evaluation of whether and to what extent model outputs support specified scientific claims. Our results show that numerical accuracy and evidential support can rank models differently, explain when and why they do so, and reveal that existing benchmarks can favor methods whose outputs provide weaker support for the scientific claims of interest. Together, we formalize, empirically demonstrate, and explain this evaluation--use mismatch in AI for PDEs.
Figures & tables
| Task | Studies | Evaluated output | Reference target | Tested constraint | Reported criterion |
| Forward | |||||
| [1.0pt/1.2pt] F1: Single-instance solution approximation | [ 3 , 11 ] | Solution field | Analytic or high-fidelity solution | — | Solution-field error |
| F2: Parametric operator learning | [ 12 , 13 ] | Solution for an unseen instance or geometry | Numerical or manufactured solution | — | Solution-field error |
| F3: Time stepping and rollout | [ 9 , 14 , 15 ] | Future trajectory | High-fidelity simulated trajectory | — | One-step or rollout error |
| F4: Probabilistic forward prediction | [ 16 , 17 ] | Predictive distribution | Held-out field or trajectory realizations | — | Probabilistic score, calibration or coverage, and predictive-mean error |
| Inverse | |||||
Appendix figures & tables10 assets
Supplementary material from the paper’s appendix.
Appendix
| Component | Submitted artifact | Independent unit | Registered / complete-primary contracts |
|---|---|---|---|
| PDEBench dynamics | Predicted future field or trajectory | Complete trajectory | 30 lineages 2 claims = 60 |
| PDEBench Darcy | Predicted solution field | Coefficient–solution pair | 5 lineages 2 claims = 10 |
| PDEInvBench scalar | Recovered scalar parameter | Complete trajectory group, including all windows sharing an initial condition | 15 registered cells 2 claims = 30; 9 cells 2 claims = 18 complete-primary |
| PDEInvBench Darcy | Recovered coefficient field | Complete Darcy field instance | 2 regimes 2 claims = 4 |
| PDEBench family | Regional quantity | Threshold/extent | Flux/transport | Event/time | Executed primary pair |
|---|---|---|---|---|---|
| 1D Advection | A | A | A | A | Regional + frozen threshold/event |
| 1D Burgers | A | A | A | A | Regional + frozen threshold/event |
| 1D Diffusion–Reaction | C | C | C | C | Regional + frozen threshold/event |
| 1D Diffusion–Sorption | A | A | C | A | Regional + frozen threshold/event |
| 2D Diffusion–Reaction | C | C | C | C | Regional + frozen threshold/event |
| 1D Compressible Navier–Stokes | C | C | C | C | Regional + frozen threshold/event |
| PDEInvBench target | Interval | Threshold/regime | Phase fraction | Regional contrast | Registered claim pair |
|---|---|---|---|---|---|
| Reaction–diffusion | A | A | N/A | N/A | Interval + directional threshold |
| Reaction–diffusion | A | C | N/A | N/A | Interval + audited directional threshold |
| Unforced Navier–Stokes viscosity | A | A | N/A | N/A | Interval + viscosity-regime threshold |
| Forced Navier–Stokes viscosity | A | A | N/A | N/A | Interval + viscosity-regime threshold |
| KdV | A | C | N/A | N/A | Interval + audited directional threshold |
| Darcy coefficient field | N/A | N/A | A | A | Regional high-phase fraction + ordered contrast |
| Support ID | Frozen scope | Compared methods |
|---|---|---|
| PDB-DYN-C6-30L | 30 dynamic lineages in ten PDE families | FNO, U-Net, DPOT-S-adapt, DPOT-M-adapt, VCNeF, OmniArch-B |
| PDB-PINN-C3-13L | 13 released-interface dynamic lineages | FNO, U-Net, PINN |
| PDB-DARCY-C2-5L | Five steady-Darcy lineages | FNO, U-Net |
| INV-SCALAR-C6-15R | Five registered scalar targets in 15 cells; complete primary: three targets in nine cells | FNO, ResNet, SC-FNO, CONFIDE, FNOPE, ten-start L-BFGS-B |
| INV-DARCY-C6-2R | Darcy ID and OOD-extreme | FNO, ResNet, DeepONet, iFNO, FNOPE, DGenNO |
| Audit | PDEBench | PDEInvBench | Outcome |
|---|---|---|---|
| Oracle attainability | Every contract has 96/96 supported focal-true and 96/96 refuted complement-true separated cases; simultaneous worst-case lower bound 92.1%; boundary-near unknown range 41/48–48/48 | Every contract has 96/96 in both separated directions; lower bound 92.8%; boundary-near unknown range 40/48–48/48 | All 104 contracts pass |
| Specificity | On 13,440 records per control, correct-resolution counts are 958, 722, 489, and 267; oracle–control separations are , , , and | On 6,528 records per control, counts are 558, 397, 291, and 169; separations are , , , and | Every simultaneous lower bound exceeds 0.75 |
| Calibration and direction | Coverage 94.6–96.8%; all-assigned wrong-direction 63/24,360 (0.259%); invalid 44/24,360 (0.181%) | Coverage 93.8–97.1%; all-assigned wrong-direction 38/6,976 (0.545%); invalid 34/6,976 (0.487%) | Near-nominal coverage; low aggregate error rates |
| Provenance | Zero prohibited exposures or evaluator-side solver calls; 11/11 nominal tracks pass | Zero prohibited exposures or calls; 12/12 nominal tracks pass | 23/23 nominal tracks pass |
| Benchmark | Method | (simultaneous 95% interval) | point rank | (simultaneous 95% interval) | point rank |
|---|---|---|---|---|---|
| PDEBench dynamics | OmniArch-B | 1 | 5 | ||
| DPOT-M-adapt | 2 | 2 | |||
| DPOT-S-adapt | 3 | 4 | |||
| VCNeF | 4 | 1 | |||
| FNO | 5 | 3 | |||
| U-Net | 6 | 6 |
| Record | -winner confidence set | -winner confidence set | Comparison |
|---|---|---|---|
| PDB-DYN-C6-RGNL | {OmniArch-B} | {VCNeF} | Disjoint |
| PDB-DYN-C6-EVNT | {OmniArch-B} | {VCNeF} | Disjoint |
| PDB-PINN-C3-RGNL | {FNO} | {FNO} | Match |
| PDB-PINN-C3-EVNT | {FNO} | {FNO, PINN} | Unresolved overlap |
| PDB-DARCY-C2-REGN | {U-Net} | {U-Net} | Match |
| PDB-DARCY-C2-FLUX | {U-Net} | {U-Net} | Match |
| Joint stratum | Whole groups | Recomputations | Disjoint selections | Matching selections | Unresolved |
|---|---|---|---|---|---|
| Adverse | 17/17 | 152 | 107 (70.4%) | 21 (13.8%) | 24 (15.8%) |
| Favorable | 16/17 | 148 | 13 (8.8%) | 119 (80.4%) | 16 (10.8%) |
| Intermediate | 17/17 | 84 | 35 (41.7%) | 31 (36.9%) | 18 (21.4%) |
| Intervention | Assigned | -equivalence band | ||
|---|---|---|---|---|
| Move equal-energy forward error into the claim region | 120 | |||
| Reverse forward error toward the adverse direction | 120 | |||
| Reverse inverse signed bias at matched error magnitude | 108 | |||
| Expand the inverse evidence half-width | 132 |
| Predictor | Multiclass log loss | Balanced accuracy |
|---|---|---|
| Mechanism features | ||
| Validation- profile only | ||
| Validation- profile task difficulty |
Scientific Claim Verification
Scientific claim verification checks whether research statements are supported, contradicted, or left uncertain by relevant evidence.
See all 54 papers