REVIEW 6 minor 42 references
Deep neural networks can provably solve Bellman equations for Markov decision processes without the curse of dimensionality
T0 review · 0 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read This paper proves that Q-functions of discounted Markov decision processes can be approximated by leaky ReLU neural networks with polynomial parameter growth in the state-space dimension and in the reciprocal error, provided the reward…
desk verdict The central theorem survives the stress test: the 'every probability measure' claim is justified because w≡1 satisfies Corollary 5.4's Lyapunov condition, and the apparent counterexample violates the explicit first-moment hypothesis in Theorem 1.1. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the full-history recursive multilevel fixed-point (MLFP) approximation scheme of reference [5], coupled with an explicit ANN calculus. The scheme replaces the Bellman fixed point by a multilevel Monte Carlo recursion in which each level averages over $M^{n-l}$ independent nested transition samples and subtracts lower-level corrections, and the paper's stability estimates show that the error of the $n$-th level decays geometrically in $n$. The ANN side supplies the concrete representations: a two-layer leaky ReLU block realizes the finite maximum function, composition and parallelization combine networks for the reward and transition maps, and linear combinations sum the multilevel terms. The key meta-identity is that every MLFP iterate is itself a function of the same form—composed, summed, and maxed primitives—so the whole iteration compiles into one leaky ReLU network whose depth and width stay polynomial.
What would settle it
Take an MDP with an analytically known Q-function, such as a linear-quadratic regulator or a threshold policy, run the MLFP network construction at dimensions $d=10,100,1000$ and accuracies $\varepsilon=10^{-2},10^{-4},10^{-6}$, and record the smallest parameter count $P$ achieving the $L^2$ error; the central claim predicts $P \le c d^c \varepsilon^{-c}$ for a fixed $c$, so a growth rate visibly faster than any polynomial in $d$ and $1/\varepsilon$ would falsify it.
Extended reading notes
Core claim
The paper's central claim, stated as Theorem 1.1 and generalized in Theorem 5.1, is that the unique bounded solution $\mathscr{q}$ of the Bellman equation $(\mathscr{q}(x))(a) = (g_d(x))(a) + \delta \mathbb{E}[\max_{b \in A} (\mathscr{q}(t_d^{(a)}(x, \xi_d)))(b)]$ admits, for every $d \in \mathbb{N}$, every $\varepsilon \in (0,1)$, and every probability measure $\mu$, a leaky ReLU network $Q$ with $\left(\int_{\mathbb{R}^d} \|\mathscr{q}(x) - (\mathcal{R}(Q))(x)\|^2 \mu(dx)\right)^{1/2} \le \varepsilon$ and $P(Q) \le c d^c \varepsilon^{-c}$. The mechanism is inheritance: the Bellman equation is a fixed-point equation, and the multilevel fixed-point scheme approximates its solution using only compositions, sums, and finite maxima of the primitive functions $g_d$ and $t_d^{(a)}$; each of those operations is realized by a leaky ReLU network, so the network representation of the fixed point has the same activation type and polynomial complexity as the primitives. Uniqueness and existence of $\mathscr{q}$ are also established under the stated contraction and Lipschitz assumptions.
Load-bearing premise
The load-bearing premise is that the primitive data—the reward function $g_d$ and the random transition maps $t_d^{(a)}$—are themselves approximable by leaky ReLU networks with error at most $\varepsilon c d^c$ and at most $c d^c \varepsilon^{-c}$ parameters; the proof also takes the MLFP error estimate (reference [5, Proposition 3.8]) as an external input rather than proving it.
Editorial extensions
If this is right
- For every MDP in the assumed class, the Q-function has a leaky ReLU representation with at most $c d^c \varepsilon^{-c}$ parameters, so achieving error $\varepsilon$ remains feasible when the state-space dimension grows.
- Because the maximum over the finite action set is built into the network, the approximation applies directly to the optimal Q-function and not to a softened or regularized surrogate.
- The same $L^2$ guarantee holds for any probability measure $\mu$ on the state space, so the result is compatible with the measure-based error notions used in reinforcement learning.
- The construction is fully explicit: given networks for $g_d$ and $t_d^{(a)}$, one can write down the Q-network by composing the MLFP iteration, which means the theorem is a proof of existence rather than an abstract counting argument.
Reading between the lines
- Editorial extension: the same MLFP-plus-ANN argument should carry over to finite-horizon problems and to related fixed-point equations, because the paper's stability machinery only uses a contraction in Wasserstein distance and does not appear to need the discount factor in a special way.
- Editorial extension: a practical reading is that deep Q-learning architectures could be initialized or structured with the MLFP network when the transition model is known; the paper does not address training dynamics or sample complexity.
- Editorial extension: the theorem locates the actual bottleneck of high-dimensional reinforcement learning in the primitive functions—problems whose rewards or transitions oscillate too fast to be leaky-ReLU approximable at polynomial cost will remain hard regardless of the Bellman solver.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper proves that, for a class of discrete-time discounted Markov decision processes, the Q-function can be approximated in L^2 by leaky ReLU neural networks with parameter counts growing at most polynomially in the state dimension d and in 1/ε. The main theorem, Theorem 1.1, assumes that the reward and transition functions are themselves approximable pointwise by such networks with polynomial complexity, together with Lipschitz regularity and a global first-moment condition on the transitions. The proof builds on the full-history recursive MLFP scheme imported from Beck, Jentzen, Kleinberg, and Kruse [5], and combines it with a detailed ANN calculus that represents MLFP iterations as explicit networks. The paper states Theorem 1.1 as an immediate consequence of Corollary 5.4, which is formulated with a general weight function w_d.
Significance. If the result holds, it is a meaningful contribution: it transfers the curse-of-dimensionality-free approximation results previously obtained for MLP/DNN methods for PDEs to Bellman equations and Q-functions, a central object in reinforcement learning. The proof is detailed, with explicit parameter bounds, and the ANN constructions are concrete. The main theorem is conditional on strong structural assumptions—pointwise DNN approximability of the primitives and a uniform-in-state first moment of the transition—but these are stated explicitly and are not hidden. The paper relies on [5, Proposition 3.8] as an external input; since [5] is a published parameter-free MLFP error bound, this is an acceptable, though visible, dependence. The central derivation is coherent, and the stress-test concern about the passage from Corollary 5.4 to Theorem 1.1 is resolved by the uniform first-moment assumption already present in Theorem 1.1, with the choice w_d≡1.
minor comments (6)
- [§1, Theorem 1.1 and §5.3, Corollary 5.4] The claimed 'immediate consequence' is valid but not transparent: taking w_d ≡ 1 and λ=1 reduces hypotheses (5.59)-(5.64) to assumptions (1.3)-(1.4) of Theorem 1.1, with the uniform first-moment condition in Theorem 1.1 supplying (5.62). Because the stress-test concern turns on this reduction, I recommend adding one explicit sentence in §1 to spell out this verification.
- [§5.3, proof of Corollary 5.4] The definition 'q_d(x))(a) = (g_d(x))(a) + (u_d(x))(x)' is not well formed, since x is a vector in R^d and u_d(x) is indexed by actions; it should read '(u_d(x))(a)'.
- [Corollary 5.4, proof] The line 'let K = max{2|A|, c, c}' repeats c; one of the entries is presumably a different constant, and the notation should be cleaned up.
- [§1, Theorem 1.1 and §2.2, Lemma 2.6] There are typographical slips, including '∪d,∈N' in the introduction and the stray 'y' in the integral inequality (2.6); these should be corrected.
- [§5.2, Corollary 5.3 proof] Equation (5.71) writes 'Rϕ(Sd,ε)' where the ANN is S_{|A|,2}; the subscript is inconsistent with the construction in (5.69).
- [References] Reference [25] duplicates [24], and the name 'Tistsiklis' in [11] is misspelled; the reference list should be deduplicated and corrected.
Circularity Check
No significant circularity: the Q-function approximation is derived compositionally from assumed DNN approximability of the primitive data and from an external MLFP error bound.
full rationale
The derivation chain is Theorem 1.1 -> Corollary 5.4 -> Corollary 5.3 -> Theorem 5.1, with the MLFP error estimate imported from [5, Proposition 3.8]. The assumptions in (1.3) concern DNN approximations of the reward g_d and the transition maps t_d^{(a)}; the conclusion in (1.6) concerns DNN approximations of the Q-function q_d. These are distinct objects. No parameter is fitted to q_d: the approximating ANN Q is constructed explicitly, as a leaky-ReLU network realizing the MLFP iteration, by composing the given ANNs G_{d,epsilon}, t_{d,epsilon}^{(a)}, and the max-network built in Lemmas 4.1 and 4.2. Lemma 4.3 realizes the MLFP approximation as an ANN, and the error estimate is then obtained from the stability result in Proposition 2.9 and from [5, Proposition 3.8]. The cited MLFP result is a published, parameter-free theorem about Monte Carlo estimators for fixed-point equations, not about DNN approximation, and its assumptions do not include the target claim. The overlapping authorship is therefore self-citation, but it is not circular in the prohibited sense. The same is true for the ANN-calculus lemmas reused from [1], [18], and [19], which are standard composition and parallelization facts. The statement that Theorem 1.1 is an immediate consequence of Corollary 5.4 is also structurally coherent: Theorem 1.1 explicitly assumes sup_x E[||t_d^{(a)}(x, xi_d)||] < infinity, so taking w_d identically equal to 1, lambda = 1, and absorbing constants into c makes the weight hypotheses (5.62)-(5.64) satisfiable, and the measure restriction in Corollary 5.4, namely (int |w_d(x)|^2 mu(dx))^{1/2} <= c d^c, holds for every probability measure mu. Whether all details of that application are fully spelled out is a presentation or correctness matter, not circularity. No equation in the paper is defined in terms of the conclusion, and no fitted value is relabeled as a prediction.
Assumptions & free parameters
assumptions (5)
- domain assumption For each d and ε, the reward g_d and transition t_d are approximable by leaky ReLU ANNs with pointwise error at most ε c d^c and parameter count at most c d^c ε^{-c} (Theorem 1.1, Eq. (1.3)).
- domain assumption The stochastic kernels and transition kernels satisfy the moment and contraction conditions in (5.4)-(5.6), including sup_d λ_d L_d < 1 and sup_{d,ε} λ_d ℓ_{d,ε} < 1.
- domain assumption The full-history recursive multilevel fixed-point error bound from [5, Proposition 3.8] is valid as stated.
- standard math The ANN calculus in Section 3, including composition, parallelization, linear combinations, and the leaky ReLU maximum representation in Lemma 4.2, holds as stated.
- standard math Measurability and integrability of the random fields are assumed so that conditional expectations, Fubini's theorem, and coupling arguments apply.
Cite this review
Pith. "Pith review of Deep neural networks can provably solve Bellman equations for Markov decision processes without the curse of dimensionality." pith.science (2026). https://pith.science/paper/6WGDUODO
@misc{pith2026250622851,
author = {Pith},
title = {Pith review of: Deep neural networks can provably solve Bellman equations for Markov decision processes without the curse of dimensionality},
year = {2026},
howpublished = {\url{https://pith.science/paper/6WGDUODO}},
note = {Machine review of arXiv:2506.22851}
}
abstract
Discrete time stochastic optimal control problems and Markov decision processes (MDPs) are fundamental models for sequential decision-making under uncertainty and as such provide the mathematical framework underlying reinforcement learning theory. A central tool for solving MDPs is the Bellman equation and its solution, the so-called $Q$-function. In this article, we construct deep neural network (DNN) approximations for $Q$-functions associated to MDPs with infinite time horizon and finite control set $A$. More specifically, we show that if the the payoff function and the random transition dynamics of the MDP can be suitably approximated by DNNs with leaky rectified linear unit (ReLU) activation, then the solutions $Q_d\colon \mathbb R^d\to \mathbb R^{|A|}$, $d\in \mathbb{N}$, of the associated Bellman equations can also be approximated in the $L^2$-sense by DNNs with leaky ReLU activation whose numbers of parameters grow at most polynomially in both the dimension $d\in \mathbb{N}$ of the state space and the reciprocal $1/\varepsilon$ of the prescribed error $\varepsilon\in (0,1)$. Our proof relies on the recently introduced full-history recursive multilevel fixed-point (MLFP) approximation scheme.
Reference graph
Works this paper leans on
-
[5]
Beck, C., Jentzen, A., Kleinberg, K., and Kruse, T. Nonlinear Monte Carlo methods with polynomial runtime for Bellman equations of discrete time high-dimensional stochastic optimal control problems.Appl. Math. Optim. 91, 1 (2025), 1–42
work page 2025
-
[1]
Ackermann, J., Jentzen, A., Kruse, T., Kuckuck, B., and Padgett, J. L. Deep neural networks with ReLU, leaky ReLU, and softplus activation provably overcome the curse of dimensionality for kolmogorov partial differential equations with lipschitz nonlin- earities in thelp-sense. arXiv:2309.13722 (2023)
work page Pith review arXiv 2023
-
[2]
Ackermann, J., Jentzen, A., Kuckuck, B., and Padgett, J. L. Deep neu- ral networks with ReLU, leaky ReLU, and softplus activation provably overcome the curse of dimensionality for space-time solutions of semilinear partial differential equations. arXiv:2406.10876 (2024)
work page Pith review arXiv 2024
-
[3]
Beck, C., Gonon, L., and Jentzen, A. Overcoming the curse of dimensionality in the numerical approximation of high-dimensional semilinear elliptic partial differential equa- tions. Partial Differ. Equ. Appl. 5, 31 (2024)
work page 2024
-
[4]
Beck, C., Hornung, F., Hutzenthaler, M., Jentzen, A., and Kruse, T. Over- coming the curse of dimensionality in the numerical approximation of Allen-Cahn partial differential equations via truncated full-history recursive multilevel Picard approximations. J. Numer. Math. 28, 4 (2020), 197–222
work page 2020
-
[6]
Beck, C., Jentzen, A., and Kruse, T. NonlinearMonteCarlomethodswithpolynomial runtime for high-dimensional iterated nested expectations.arXiv:2009.13989 (2020)
arXiv 2020
-
[7]
Becker, S., Braunw arth, R., Hutzenthaler, M., Jentzen, A., and von Wurstemberger, P. Numerical simulations for full history recursive multilevel Picard approximations for systems of high-dimensional partial differential equations. Commun. Comput. Phys. 28, 5 (2020), 2109–2138
work page 2020
-
[8]
Bellman, R. Dynamic programming. Princeton Landmarks in Mathematics. Princeton University Press, Princeton, NJ, 2010. Reprint of the 1957 edition, With a new introduction by Stuart Dreyfus
work page 2010
Show all 42 references
-
[9]
Reinforcement learning and optimal control
Bertsekas, D. Reinforcement learning and optimal control. Athena Scientific, 2019
2019
-
[10]
Bertsekas, D., and Shreve, S. E. Stochastic optimal control: the discrete-time case, vol. 5. Athena Scientific, 1996. 47
1996
-
[11]
Bertsekas, D., and Tsitsiklis, J. N. Neuro-dynamic programming. Athena Scientific, 1996
1996
-
[12]
D., and W ang, Z
Cai, Q., Yang, Z., Lee, J. D., and W ang, Z. Neural temporal-difference learning converges to global optima.Advances in Neural Information Processing Systems 32(2019)
2019
-
[13]
Cioica-Licht, P., Hutzenthaler, M., and Werner, P. T. Deep neural networks overcome the curse of dimensionality in the numerical approximation of semilinear partial differential equations. arXiv:2205.14398 (2022)
2022 arXiv
-
[14]
E, W., Hutzenthaler, M., Jentzen, A., and Kruse, T. On multilevel Picard numer- ical approximations for high-dimensional nonlinear parabolic partial differential equations and high-dimensional nonlinear backward stochastic differential equations.J. Sci. Comput. 79, 3 (2019), 1534–1571
2019
-
[15]
Multilevel Picard iterations for solving smooth semilinear parabolic heat equations
E, W., Hutzenthaler, M., Jentzen, A., and Kruse, T. Multilevel Picard iterations for solving smooth semilinear parabolic heat equations. Partial Differ. Equ. Appl. 2, 6 (2021), 80
2021
-
[16]
A theoretical analysis of deep Q-learning
F an, J., W ang, Z., Xie, Y., and Yang, Z. A theoretical analysis of deep Q-learning. In Learning for dynamics and control(2020), PMLR, pp. 486–489
2020
-
[17]
B., Jentzen, A., and Welti, T
Giles, M. B., Jentzen, A., and Welti, T. Generalised multilevel Picard approxima- tions. arXiv:1911.03188 (2019)
2019 arXiv
-
[18]
Space-time error estimates for deep neural network approximations for differential equations.Adv
Grohs, P., Hornung, F., Jentzen, A., and Zimmermann, P. Space-time error estimates for deep neural network approximations for differential equations.Adv. Comput. Math. 49, 1 (2023), 4
2023
-
[19]
Deep neural network approximations for solutions of PDEs based on Monte Carlo algorithms.Partial Differ
Grohs, P., Jentzen, A., and Salimov a, D. Deep neural network approximations for solutions of PDEs based on Monte Carlo algorithms.Partial Differ. Equ. Appl. 3, 4 (2022), 45
2022
-
[20]
Overcoming the curse of dimen- sionality in the numerical approximation of parabolic partial differential equations with gradient-dependent nonlinearities
Hutzenthaler, M., Jentzen, A., and Kruse, T. Overcoming the curse of dimen- sionality in the numerical approximation of parabolic partial differential equations with gradient-dependent nonlinearities. Found. Comput. Math.(2021), 1–62
2021
-
[21]
Hutzenthaler, M., Jentzen, A., Kruse, T., and Nguyen, T. A. Multilevel Picard approximations for high-dimensional semilinear second-order PDEs with Lipschitz nonlin- earities. arXiv:2009.02484 (2020)
2020 arXiv
-
[22]
Hutzenthaler, M., Jentzen, A., Kruse, T., and Nguyen, T. A. A proof that rectified deep neural networks overcome the curse of dimensionality in the numerical ap- proximation of semilinear heat equations.Partial Differ. Equ. Appl. 1(2020), 1–34
2020
-
[23]
A., and von Wurstem- berger, P
Hutzenthaler, M., Jentzen, A., Kruse, T., Nguyen, T. A., and von Wurstem- berger, P. Overcoming the curse of dimensionality in the numerical approximation of semilinear parabolic partial differential equations. Proc. R. Soc. A 476, 2244 (2020), 20190630
2020
-
[24]
Overcoming the curse of dimensionality in the approximative pricing of financial derivatives with default risks
Hutzenthaler, M., Jentzen, A., and von Wurstemberger, P. Overcoming the curse of dimensionality in the approximative pricing of financial derivatives with default risks. Electron. J. Probab. 25(2020), Paper No. 101, 73. 48
2020
-
[25]
Overcoming the curse of dimensionality in the approximative pricing of financial derivatives with default risks
Hutzenthaler, M., Jentzen, A., and von Wurstemberger, P. Overcoming the curse of dimensionality in the approximative pricing of financial derivatives with default risks. Electron. J. Probab. 25(2020), Paper No. 101, 73
2020
-
[26]
Multilevel Picard approximations of high- dimensional semilinear parabolic differential equations with gradient-dependent nonlinear- ities
Hutzenthaler, M., and Kruse, T. Multilevel Picard approximations of high- dimensional semilinear parabolic differential equations with gradient-dependent nonlinear- ities. SIAM J. Numer. Anal. 58, 2 (2020), 929–961
2020
-
[27]
Mathematicalintroduction to deep learning: Methods, implementations, and theory, 2023
Jentzen, A., Kuckuck, B., and von Wurstemberger, P. Mathematicalintroduction to deep learning: Methods, implementations, and theory, 2023
2023
-
[28]
Champion-level drone racing using deep reinforcement learning.Nature 620, 7976 (2023), 982–987
Kaufmann, E., Bauersfeld, L., Loquercio, A., Müller, M., Koltun, V., and Scaramuzza, D. Champion-level drone racing using deep reinforcement learning.Nature 620, 7976 (2023), 982–987
2023
-
[29]
P., Hunt, J
Lillicrap, T. P., Hunt, J. J., Pritzel, A., Heess, N., Erez, T., Tassa, Y., Sil- ver, D., and Wierstra, D. Continuous control with deep reinforcement learning.4th International Conference on Learning Repre- sentations, ICLR 2016(2016)
2016
-
[30]
A., Veness, J., Bellemare, M
Mnih, V., Ka vukcuoglu, K., Sil ver, D., Rusu, A. A., Veness, J., Bellemare, M. G., Gra ves, A., Riedmiller, M., Fidjeland, A. K., Ostrovski, G., et al. Human-level control through deep reinforcement learning.nature 518, 7540 (2015), 529– 533
2015
-
[31]
Neufeld, A., and Nguyen, T. A. Multilevel Picard approximations and deep neural networks with ReLU, leaky ReLU, and softplus activation overcome the curse of dimension- ality when approximating semilinear parabolic partial differential equations in Lp-sense. arXiv:2409.20431 (2024)
2024
-
[32]
Neufeld, A., and Nguyen, T. A. Rectified deep neural networks overcome the curse of dimensionality in the numerical approximation of gradient-dependent semilinear heat equations. arXiv preprint arXiv:2403.09200(2024)
2024 arXiv
-
[33]
A., and Wu, S
Neufeld, A., Nguyen, T. A., and Wu, S. Deep ReLU neural networks overcome the curseofdimensionalitywhenapproximatingsemilinearpartialintegro-differentialequations. arXiv:2310.15581 (2023)
2023 arXiv
-
[34]
A., and Wu, S
Neufeld, A., Nguyen, T. A., and Wu, S. Multilevel Picard approximations overcome the curse of dimensionality in the numerical approximation of general semilinear PDEs with gradient-dependent nonlinearities. Journal of Complexity(2025), 101946
2025
-
[35]
Neufeld, A., and Wu, S. Multilevel Picard approximation algorithm for semilinear partial integro-differential equations and its complexity analysis.Stochastics and Partial Differential Equations: Analysis and Computations(2025), 1–59
2025
-
[36]
The curse of dimension and a universal method for numer- ical integration
Nov ak, E., and Ritter, K. The curse of dimension and a universal method for numer- ical integration. InMultivariate approximation and splines (Mannheim, 1996), vol. 125 of Internat. Ser. Numer. Math.Birkhäuser, Basel, 1997, pp. 177–187
1996
-
[37]
Tractability of multivariate problems
Nov ak, E., and Woźniakowski, H. Tractability of multivariate problems. Vol. 1: Linear information, vol. 6 ofEMS Tracts in Mathematics. European Mathematical Society (EMS), Zürich, 2008. 49
2008
-
[38]
Powell, W. B. Approximate Dynamic Programming: Solving the curses of dimensionality, vol. 703. John Wiley & Sons, 2007
2007
-
[39]
Puterman, M. L. Markov decision processes: discrete stochastic dynamic programming. John Wiley & Sons, 2014
2014
-
[40]
J., Guez, A., Sifre, L., V an Den Driess- che, G., Schrittwieser, J., Antonoglou, I., Panneershel v am, V., Lanctot, M., et al
Sil ver, D., Huang, A., Maddison, C. J., Guez, A., Sifre, L., V an Den Driess- che, G., Schrittwieser, J., Antonoglou, I., Panneershel v am, V., Lanctot, M., et al. Mastering the game of go with deep neural networks and tree search.nature 529, 7587 (2016), 484–489
2016
-
[41]
S., and Barto, A
Sutton, R. S., and Barto, A. G. Reinforcement learning: An introduction. MIT press, 2018
2018
-
[42]
A finite-time analysis of Q-learning with neural network function approximation
Xu, P., and Gu, Q. A finite-time analysis of Q-learning with neural network function approximation. InInternational conference on machine learning(2020), PMLR, pp. 10555– 10565. 50
2020
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.