Joint-Particle Distributional Critic
Implementation & benchmark of arXiv:2608.22765 — Learning to Control Coupled-Dynamics Environments with Joint Markov Decision Processes
Source paper: Learning to Control Coupled-Dynamics Environments with Joint Markov Decision Processes arXiv:2608.22765 ⓘ · analyzed Aug 29, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Replace independent per-action distributional value heads with a critic whose shared latent particle produces a vector of return samples for all actions simultaneously. Train the predicted joint return vector against a Bellman target vector formed from coupled counterfactual reward-transition samples, using a sliced Wasserstein loss. The greedy action is selected by the mean of the corresponding marginal particles, while shared particles retain cross-action dependence for learning and uncertainty estimation.
Formulas
Mathematical statement
Definition 1 defines the discounted return as $Z^{\pi}(s,a):=\sum_{t=0}^{\infty}\gamma^{t}R_t$, where $s\in\mathcal{S}$ is a state, $a\in\mathcal{A}$ is an action, $R_t$ is the reward at time $t$, and $\gamma\in(0,1)$ is the discount factor. In a JMDP, the one-step variables $(R_a,S'_a)_{a\in\mathcal{A}}$ are sampled jointly from one coupled kernel, so different counterfactual actions share exogenous randomness. The adapted joint Bellman target for particle $k$ is $Y_k=(Y_{a,k})_{a\in\mathcal{A}}$ with $Y_{a,k}=R_{a,k}+\gamma f_{\bar\theta}(S'_{a,k},\pi_{\bar\theta}(S'_{a,k}),\xi'_k)_a$, where $\bar\theta$ are target-network parameters, $\xi'_k$ is a shared latent particle, and $\pi_{\bar\theta}(s)=\arg\max_b K^{-1}\sum_{j=1}^{K}f_{\bar\theta}(s,\xi_j)_b$ is the mean-greedy policy. The critic predicts $X_k=f_\theta(s,\xi_k)\in\mathbb{R}^{A}$, where $A=|\mathcal{A}|$ and $X_{a,k}$ is a sample from the return distribution for action $a$. Use the sliced Wasserstein discrepancy $\operatorname{SW}_1(X,Y)=\mathbb{E}_{u\sim\operatorname{Unif}(\mathbb{S}^{A-1})}[W_1(\{u^\top X_k\}_{k=1}^{K},\{u^\top Y_k\}_{k=1}^{K})]$, approximated with random directions and sorted projected samples. The paper's convergence result states that iterates of the distributional Bellman optimality operator converge in Wasserstein distance to the optimal joint return law when the induced marginal MDP has a unique optimal policy; first and second moments may converge under weaker tie conditions. The implementation uses the shared coupling and Wasserstein geometry to preserve cross-action return dependence.
Implementation notes
(1) Integration point: implement this as a distributional RL critic replacing a scalar Q-network or independent quantile heads. Given a state embedding $h(s)$ and a latent particle $\xi_k\sim\mathcal N(0,I)$, an MLP outputs an $A$-dimensional vector $X_k=f_\theta(s,\xi_k)$, one return sample for every action. Use the same particle index across all action coordinates. The replay buffer must store coupled transitions $(R_a,S'_a)_{a=1}^{A}$ generated by one shared simulator random seed; ordinary single-action environments cannot test the joint advantage directly. (2) Pseudocode: sample a minibatch of coupled transitions and latent particles; compute $X_k=f_\theta(s,\xi_k)$; compute $\mu_b=K^{-1}\sum_k f_{\bar\theta}(s'_a,\xi_k)_b$ and choose $b^*(s'_a)=\arg\max_b\mu_b$; construct $Y_{a,k}=R_a+\gamma f_{\bar\theta}(S'_a,\xi'_k)_{b^*(S'_a)}$ using shared $\xi'_k$; draw $M$ random unit vectors $u_m$; project both $X_k$ and $Y_k$ onto every direction, sort the $K$ projections, and minimize the mean absolute sorted difference in the sliced-Wasserstein loss. Update $\theta$ with Adam and periodically copy it to $\bar\theta$. (3) Mathematical versus empirical quantities: the return recursion, shared-vector target, and Wasserstein objective are mathematical components; $K$, $M$, target-update period, and $\gamma$ are hyperparameters. Estimate action means by particle averages. Log empirical cross-action covariance $\widehat{\operatorname{Cov}}(X_{\cdot,a},X_{\cdot,b})$ and compare it with covariance from held-out coupled rollouts. (4) First cheap experiment: create a stochastic two-action bandit or gridworld in which both counterfactual rewards depend on a controllable shared Bernoulli or Gaussian exogenous variable. Compare scalar DQN, independent quantile regression, and the joint critic at equal network size and environment interactions. Then test a small vectorized control task with common-randomness simulator seeds. Measure held-out joint-return Wasserstein error, mean return, action-selection accuracy when means are close, and cross-action covariance calibration. Success means lower joint-law error or faster return improvement at equal interaction count; if shared-particle training gives no improvement over independently resampled transitions, reject the added architecture.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Not tested
Stage 2 · Practical benchmark result: Not run
Methodology: Toy-system gate first; the benchmark stage runs only after a pass. How verification works
Stage 1 — Mechanism check agent confidence 8/10
Built a shared-latent joint-particle critic with sliced-Wasserstein training and an independent marginal-head baseline on a coupled two-action Gaussian terminal task. The math check passed: identical particle sets had SW1=0, and translation error was 2.98e-8. The joint critic achieved lower held-out joint SW1 (0.1413 vs 0.2709) and reproduced cross-action covariance (0.750 vs true 0.706), while the independent baseline produced incorrect covariance (-0.092), demonstrating the intended dependence-preservation effect.
- Agent confidence
- 8/10
- Baseline
- Independent marginal heads: held-out SW1=0.270889, mean absolute marginal mean error=0.0684397, cross-action covariance=-0.09184 vs true 0.70614.
- Idea
- Joint-particle critic: held-out SW1=0.141323, mean absolute marginal mean error=0.0668253, cross-action covariance=0.75017 vs true 0.70614. Training runtime was 15.17 seconds on CUDA.
Limitations: This is a terminal coupled Gaussian bandit, so gamma=0 and no multi-step Bellman backup or learned mean-greedy control policy was evaluated. Results use one fixed seed and finite particles/directions; the independent baseline is designed to match marginals but does not receive coupled-pair supervision.
How to run: python3 joint_particle_experiment.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Super-audit (legacy) latest Failed ✗
Agent confidence: 9/10
Implemented and ran a 10-seed super-check with a harder nonlinear, multi-step coupled MDP, Bellman targets, target networks, mean-greedy action selection, independent-head and shared-particle ablations, and the proposed sliced-Wasserstein joint critic. The math check passed, but the joint critic beat the independent baseline on covariance error in only 6/10 seeds (paired p=0.414796), while having worse SW1 and marginal mean error on average. It fails the required 8/10 and p<0.05 criterion, so the previous toy improvement does not robustly survive.
- Baseline
- Independent heads: SW1 1.235583 +/- 0.041576; covariance error 0.063689 +/- 0.049522; marginal mean error 0.275467 +/- 0.163172. Per-seed covariance errors: [0.071819, 0.082057, 0.119519, 0.055204, 0.005183, 0.012329, 0.153291, 0.018086, 0.025947, 0.093460].
- Data
- On-the-fly nonlinear 3-D coupled Gaussian-exogenous synthetic MDP; horizon 5, gamma 0.90; 10 seeds.
- Idea
- Joint critic: SW1 1.266271 +/- 0.028910; covariance error 0.043537 +/- 0.046468; marginal mean error 0.334466 +/- 0.119019. Per-seed covariance errors: [0.000684, 0.020638, 0.007064, 0.150824, 0.030302, 0.064120, 0.057331, 0.081086, 0.012157, 0.011163]. Covariance wins: 6/10; paired p=0.414796; runtime 227.42 seconds on CUDA. Math: identical SW1=0, translation absolute error=2.98e-08.
- Shared ablation
- Shared-particle marginal heads: SW1 1.219561 +/- 0.054550; covariance error 0.051546 +/- 0.050420; marginal mean error 0.265675 +/- 0.162898.
The harder benchmark remains synthetic rather than a real environment. Only one small MLP architecture, short horizon, and training budget were tested; no hyperparameter sweep, large-scale RL benchmark, deployment speed study, or uncertainty calibration study was performed.
How to run: python3 supercheck.py
Artifacts
- joint_particle_experiment.py 6.7 KB View Raw
- make_report.py 3.0 KB View Raw
- report.md 1.5 KB View
- report_super_2026-08-30T104107.md 3.1 KB View
- results.json 1.4 KB View Raw
- super_results.json 19.8 KB View Raw
- supercheck.py 7.2 KB View Raw
- supercheck_report.md 2.2 KB View