TD-to-PDE Continuation Training
Implementation & benchmark of arXiv:2608.17117 — Physics-informed Reinforcement Learning for Stochastic Reach-Avoid Analysis
Source paper: Physics-informed Reinforcement Learning for Stochastic Reach-Avoid Analysis arXiv:2608.17117 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Train a neural value or latent-dynamics model with temporal-difference targets before enforcing a stiff differential-equation residual, and ramp the physics weight only after the critic has become predictive. For a stochastic dynamical model, the residual is computed using the infinitesimal generator, while terminal, safe, and failure boundary conditions are imposed through separate penalties.
Formulas
Mathematical statement
For the stochastic system $\mathrm{d}X_t=f(X_t,U_t)\mathrm{d}t+\sigma(X_t,U_t)\mathrm{d}W_t$, define $a(x,u)=\sigma(x,u)\sigma(x,u)^\top$ and the generator $\mathcal{L}^u V(x)=f(x,u)^\top\nabla V(x)+\frac{1}{2}\operatorname{tr}(a(x,u)\nabla^2V(x))$. A candidate reach-avoid value $V_\theta(x)$ should satisfy a stationary or discounted Bellman-HJB equation of the form $\mathcal{H}[V_\theta](x)=0$, where $\mathcal{H}[V](x)=\max_{u\in\mathbb{U}}\{\mathcal{L}^uV(x)-\beta V(x)+r(x,u)\}$ for control maximization, discount rate $\beta\geq0$, and reward/source term $r$. Given transitions $(x_t,u_t,x_{t+1},r_t)$, the TD error is $\delta_t=r_t+\gamma V_\theta(x_{t+1})-V_\theta(x_t)$, with discount factor $\gamma$. Optimize $L(\theta;\lambda)=L_{TD}(\theta)+\lambda_{PDE}L_{PDE}(\theta)+\lambda_B L_B(\theta)$, where $L_{PDE}=\mathbb{E}_{x\sim\mu}|\mathcal{H}[V_\theta](x)|^2$ and $L_B$ penalizes violations of prescribed boundary or terminal values $g(x)$. The mechanism is continuation in $\lambda_{PDE}$: begin at $\lambda_{PDE}=0$ and increase it only after TD error and value variance fall below specified thresholds.
Implementation notes
Integrate this at the critic-training step of an actor-critic agent, or at the loss level of a neural operator or world model that has known drift and diffusion. Use an MLP or transformer critic $V_\theta(x)$ and obtain first- and second-order derivatives with automatic differentiation. For each minibatch, sample environment transitions for the TD term and independently sample collocation states from replay-buffer states, uniform domain points, or boundary-focused points for the PDE and boundary terms. Pseudocode: compute $\delta=r+\gamma V_\theta(x')-V_\theta(x)$ and $L_{TD}=\operatorname{mean}(\delta^2)$; evaluate $\nabla V_\theta$ and the Hessian at collocation states; form $\mathcal{L}^uV_\theta$; maximize over a discretized action set or evaluate at the current actor action $u=\pi_\phi(x)$; compute $L=L_{TD}+\lambda_{PDE}L_{PDE}+\lambda_BL_B$; update $\theta$; and update the actor separately using $\nabla_\phi V_\theta(x)$. Increase $\lambda_{PDE}$ by factor $q$ only when an exponential moving average of TD loss is below $\tau_{TD}$ and critic-output variance is below $\tau_V$. The SDE generator and PDE form are supplied by the known model; TD loss, output variance, residual estimates, and boundary error are measured empirically. First test on a 1D or 2D stochastic navigation task with an analytic or finite-difference reach-avoid solution. Compare pure PDE-PINN, pure TD critic, simultaneous fixed-weight training, and continuation training using identical networks and collocation points. The quantitative prediction is a sharp reduction in TD loss before strong PDE-residual reduction: staged training should avoid the early-loss plateau or divergence caused by a large fixed $\lambda_{PDE}$. Sweep the initial PDE weight and measure the largest stable value; continuation should tolerate a larger final $\lambda_{PDE}$ and achieve boundary error at least as low as fixed-weight training without increasing TD error by more than 20%.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Not reproduced ✗
Stage 2 · Practical benchmark result: Significantly worse ✗
Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works
Stage 1 — Mechanism check agent confidence 8/10
Built a 1-D OU diffusion critic experiment with automatic-differentiation generator residuals, TD loss, boundary penalty, fixed-PDE, TD-only, and gated continuation schedules. The generator check passed exactly, and weak-Euler error decreased from 0.00796 to 0.000922 as dt fell by 8x, close to the predicted linear ratio 0.125 (observed 0.116). The continuation gate showed a threshold transition, and activated continuation achieved lower PDE loss than fixed weighting, although fixed weighting had slightly better value RMSE and TD loss; this is a promising mechanism signal rather than a decisive overall win.
- Agent confidence
- 8/10
- Baseline
- TD-only: TD loss 0.000911, PDE loss 0.537, value RMSE 0.375. Fixed PDE weight 10: TD loss 0.004876, PDE loss 0.001084, value RMSE 0.02934.
- Idea
- Continuation at TD threshold 0.002: first ramp at step 5, final lambda 10, TD loss 0.005029, PDE loss 0.0000661, value RMSE 0.03167. Generator autodiff error 0; weak-generator error ratio 0.116 observed versus 0.125 predicted.
Limitations: Only a small manufactured 1-D diffusion and one random seed were tested; no actor learning, action maximization, reach-avoid boundaries, stochastic-policy evaluation, hyperparameter robustness, wall-clock comparison, or larger architectures were evaluated.
How to run: python3 run_experiment.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.17117)
Significantly WORSE than the tuned baseline (+340.5%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Poisson_boundary
- Model
- smooth_mlp_shared
- Paired seeds
- 8
- Baseline mean
- 0.0197
- Idea mean
- 0.0868
- Effect (Δ)
- 0.0671 (+340.5%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±42.5%
- Mechanism
- Not confirmed ✗
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
Only the registered poisson_boundary PDE track and one shared smooth MLP architecture were tested. Genuine stochastic TD transitions, actor-critic training, action maximization, alternative PDE weights beyond the tested continuation sweep, and larger models were not evaluated.
How to run: python3 run_bench.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Stage 2 — Benchmark latest Failed ✗
Agent confidence: 10/10
Built and evaluated TD-to-PDE continuation on the registered poisson_boundary PDE track using a shared smooth MLP. Baseline test MSE was 0.0197138, while the best continuation setting achieved 0.0868471; paired delta was +0.0671333 with permutation p=0.0081, so the idea was significantly worse. The trained-model mechanism signature was also not confirmed because PDE residual increased from 0.0057077 early to 0.7555173 late.
- Baseline
- poisson_boundary, smooth_mlp_shared, tuned lr=0.01; 8-seed mean test MSE 0.019713770831003785
- Idea
- poisson_boundary, smooth_mlp_shared, continuation lr=0.003 and PDE weight 1.0; 8-seed mean test MSE 0.0868471022695303
Benchmark result
Significantly WORSE than the tuned baseline (+340.5%, p=0.0081, wins 0 / 8 paired seeds).
- Benchmark
- Poisson_boundary
- Model
- smooth_mlp_shared
- Paired seeds
- 8
- Baseline mean
- 0.0197
- Idea mean
- 0.0868
- Effect (Δ)
- 0.0671 (+340.5%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 0.0081 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±42.5%
- Mechanism
- Not confirmed ✗
- Practical effect
- Hurts
- Baseline tuning
- swept over 3 configs
Only the registered poisson_boundary PDE track and one shared smooth MLP architecture were tested. Genuine stochastic TD transitions, actor-critic training, action maximization, alternative PDE weights beyond the tested continuation sweep, and larger models were not evaluated.
How to run: python3 run_bench.py
Artifacts
- bench_report.json 7.3 KB View Raw
- poisson_track.py 0.8 KB View Raw
- report.md 1.5 KB View
- report_bench_2026-09-01T162428.md 4.0 KB View
- results.json 1.3 KB View Raw
- run_bench.py 6.7 KB View Raw
- run_experiment.py 4.6 KB View Raw