REVIEW 3 major objections 4 minor 2 cited by
Generative Flow Networks: Theory and Applications to Structure Learning
T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The thesis establishes that sampling from an intractable discrete distribution reduces to finding edge flows that satisfy a conservation law, and that this is enough to approximate Bayesian posteriors over causal graph structures.
desk verdict A rigorous thesis-length synthesis of GFlowNets; the theory holds up, but the practical claims still depend on the usual amortized-inference gap and the genuinely new content is thin. 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 a Markovian flow network on a pointed DAG: a non-negative edge flow $F(s\to s')$ with a single source $s_0$ and a single sink $\bot$. The flow matching condition enforces conservation at every non-initial state; the boundary condition $F(x\to\bot)=R(x)$ fixes the terminating flows. The identity that carries the argument is the induced forward policy $P_F^\star(s'|s)=F(s\to s')/F(s)$, which, by Proposition 3.1.8, samples terminating states exactly with probability $F(x\to\bot)/Z^\star$. Detailed balance, trajectory balance, and sub-trajectory balance are alternative re-expressions of the same conservation law at transition, full-trajectory, and partial-trajectory granularity, each usable as a training loss.
What would settle it
On a small, fully enumerable state space (for example, all DAGs on four nodes), compute the exact terminating-state distribution $P_F^\top(x)$ by dynamic programming, train a GFlowNet with a flow matching loss, and test whether the learned distribution approaches $R(x)/Z$ as training time and capacity grow; a persistent gap would falsify the claim that flow matching suffices for sampling.
Extended reading notes
Core claim
The central claim is that a non-negative reward function $R$ on a finite set of compositional objects can be realized as the terminating flow of a Markovian flow network. Setting $F(x\to\bot)=R(x)$ on each terminating edge and choosing edge flows satisfying the flow matching condition $\sum_{s\in\mathrm{Pa}(s')}F(s\to s')=\sum_{s''\in\mathrm{Ch}(s')}F(s'\to s'')$ at every non-initial state yields a forward policy $P_F(s'|s)=F(s\to s')/F(s)$ whose terminating-state distribution is exactly $R(x)/\sum_x R(x)$ (Proposition 3.1.8). The thesis further claims that this conservation law can be turned into a family of learning objectives — flow matching, detailed balance, trajectory balance, and sub-trajectory balance — and that training a neural-network policy against these losses yields an amortized sampler that approximates the target distribution, with convergence guarantees for exact flows. In structure learning, this yields DAG-GFlowNet for the marginal posterior over graphs and JSP-GFN for the joint posterior over graph and parameters, and the Markov-chain view extends the construction to continuous state spaces.
Load-bearing premise
The practical claim assumes that a parametric flow trained with flow-matching losses can be driven close enough to exact conservation over a huge discrete space; the theorems are exact only for exactly consistent flows and do not bound the error from finite-sample, function-approximation training.
Editorial extensions
If this is right
- A GFlowNet trained with a flow matching loss can sample DAGs from the Bayesian posterior $P(G|D)$ rather than from a path-count-biased distribution, because the conservation law cancels trajectory multiplicities.
- Off-policy training with a behavior policy can be used to improve exploration without changing the target distribution, and the convergence guarantees hold for the exact flow matching objectives.
- GFlowNets and maximum-entropy reinforcement learning are equivalent after a reward correction, so the known multi-path bias of MaxEnt RL disappears when the reward is corrected as in Chapter 5.
- The Markov-chain perspective extends GFlowNets to continuous and general state spaces, giving generalized detailed balance and trajectory balance conditions.
- A single GFlowNet can approximate the joint posterior over the structure and the parameters of a Bayesian network, with DAG-GFlowNet as a special case.
Reading between the lines
- A reader could test the practical limit directly: on small enumerable state spaces, compute $P_F^\top(x)$ exactly by dynamic programming and measure how the flow-matching residual shrinks with capacity; this would separate training error from representational bias.
- The same conservation-law sampler suggests an amortized proposal for MCMC or a reusable sampler for other discrete combinatorial objects, a connection the thesis only touches on.
- Because the exact theorems are for exact flows, the practical quality of a GFlowNet may be predicted by the degree to which learned edge flows satisfy conservation at intermediate states; monitoring this residual is an implicit diagnostic.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The thesis develops the theory of Generative Flow Networks (GFlowNets) as a class of sequential, amortized generative models for sampling from unnormalized distributions over discrete and compositional objects, and applies them to Bayesian structure learning. Part I introduces flow networks, defines Markovian flows, and proves that exact flows satisfying the flow matching condition (3.2.1), detailed balance (3.2.21), trajectory balance (3.2.31), or sub-trajectory balance (3.2.49) conditions are equivalent characterizations of a unique Markovian flow; Proposition 3.1.8 shows that the terminating distribution of such a flow is the normalized terminating edge flow. The thesis also connects GFlowNets to maximum entropy reinforcement learning and variational inference, and extends them to general state spaces. Part II constructs GFlowNets over DAGs, approximates the Bayesian posterior over structures (DAG-GFlowNet) and over structures and parameters jointly (JSP-GFN), and validates the approximations against exact posteriors computed by dynamic programming on small graphs and against MCMC baselines.
Significance. If the central claims hold, GFlowNets offer a principled and scalable alternative to MCMC for posterior inference over discrete structures, with clear connections to variational inference and max-entropy RL. The manuscript's strengths include the careful proofs of the flow conditions, the explicit empirical check against exact posterior enumeration for small graphs (Section 7.5), the integration of interventional and observational data (Section 7.7), and the candid limitation sections (7.4.3, 7.4.4, 8.3.5). The main caveat is that the exact-flow theory does not by itself justify the trained parametric approximation; the degree to which this gap is quantified or appropriately labeled determines whether the application claims are fully supported.
major comments (3)
- [§4.2.1; Prop. 3.1.8; Thm. 3.2.1] The central practical claim is that a trained GFlowNet samples approximately from the target reward distribution. Proposition 3.1.8 and Theorem 3.2.1 establish this only for exact Markovian flows satisfying the flow matching condition (3.2.1) and the boundary reward F(x→⊥)=R(x). The training losses in Section 4.2.1 minimize sampled residuals; their global optimum is zero at best on the support of the training distribution, and nothing in the presentation ensures that the parametric family contains a valid flow or that the training distribution covers the state space. The limitation sections 7.4.3, 7.4.4, and 8.3.5 acknowledge this gap, but provide no quantitative bound. Please either supply an error bound (for example, total variation distance in terms of the expected flow-matching residual and a coverage or exploration term) or explicitly reframe the practical claims as heuristic approximations.
- [§4.2.3] The chapter introduction promises 'novel convergence guarantees of the approximate distribution towards the target' and Section 4.2.3 is titled 'Convergence guarantees'. The precise assumptions of these guarantees should be stated clearly. If they require exact minimization of the flow-matching loss over a hypothesis class that contains a valid flow, or access to the full training distribution, then they do not cover the finite-capacity, finite-sample training used in Chapters 7 and 8, and the application chapters remain conditional on an unquantified approximation error.
- [§7.5, §7.6.3] The exact-posterior comparison in Section 7.5 is limited to small graphs, while the medium-scale experiments (d=20) rely on metrics whose limitations are themselves acknowledged in Section 7.6.3. The thesis should either provide a quantitative calibration of those metrics (for example, by comparing expected SHD and AUROC against exact posteriors on graphs small enough for enumeration) or temper the scaling claims accordingly, so that the reader can distinguish validated accuracy from heuristic evaluation.
minor comments (4)
- [Résumé, p. iii] The French text reads 'lesréseaux à flots génératifs'; there is a missing space between 'les' and 'réseaux'.
- [Figure 2.2 caption] The caption uses 'quicunx'; the correct spelling is 'quincunx'.
- [Glossary, p. xiii] The glossary entry for 'P⊗T F' has inconsistent sub- and superscript rendering; please ensure it matches the notation used in Appendix A.
- [Conclusion] Since Sections 7.4.3, 7.4.4, and 8.3.5 each discuss limitations in different places, a consolidated limitations paragraph in the Conclusion that relates these caveats to the theoretical guarantees of Part I would help the reader assess the scope of the contributions.
Circularity Check
No circularity: the flow-matching theorem is proved in-text, the reward is an input, and the acknowledged approximation gap is a limitation, not a definitional reduction.
full rationale
Walking the derivation chain: the target distribution P*(x)=exp(-E(x))/Z is an input reward, not a quantity derived from the fitted model. Proposition 3.1.8 shows that a valid flow with terminating edge flow F(x->bottom)=R(x) induces the terminating distribution R(x)/Z; the proof uses Lemmas 3.1.6 and Proposition 3.1.4, not the target distribution. Theorem 3.2.1 proves that the flow matching condition (3.2.1) characterizes the edge flows of a unique Markovian flow, with a full proof in the text rather than an imported uniqueness claim. The training objectives in Section 4.2.1 minimize flow-matching residuals whose global zero is exactly the flow matching condition, and the convergence guarantees in Section 4.2.3 concern that exact minimization. The manuscript explicitly acknowledges in Sections 7.4.3, 7.4.4, and 8.3.5 that finite-sample, function-approximation, and training-distribution coverage gaps are not bounded; this is an honest limitation, not a circular step, because the theory is stated as a sufficiency result for exact flows. Empirical validation in Section 7.5.1 compares against exact posteriors computed by dynamic programming, an external benchmark, so the evaluation does not reuse GFlowNet outputs as ground truth. Although the thesis builds on several papers by the same author, the load-bearing theorems are restated and proved in the manuscript, so the self-citations are not load-bearing arguments. No fitted parameter is renamed as a prediction, no known result is merely relabeled, and no ansatz is smuggled in via citation. The derivation is self-contained relative to its stated assumptions.
Assumptions & free parameters
free parameters (2)
- temperature alpha in MaxEnt RL equivalence =
1 (set by assumption)
- prior hyperparameters for BDe/BGe scores =
not specified in thesis
assumptions (5)
- domain assumption Faithfulness: the data-generating distribution is faithful to the true DAG (Section 1.4.1).
- standard math The state space is a finite pointed DAG with a unique source and sink (Definitions 2.2.1-2.2.3).
- ad hoc to paper Existence of a strictly positive Markovian flow matching the boundary reward (implicit in Theorem 3.2.1).
- domain assumption In Bayesian structure learning, the marginal likelihood (BGe/BDe score) is computable and used as the unnormalized posterior (Section 7.1, Appendix B).
- domain assumption For continuous-state extensions, the generalized Markov chain is Harris recurrent with an accessible atom (Section 6.2).
invented entities (2)
-
Flow (trajectory flow, state flow, edge flow)
independent evidence
-
Generative Flow Network (GFlowNet) as a model class
independent evidence
Cite this review
Pith. "Pith review of Generative Flow Networks: Theory and Applications to Structure Learning." pith.science (2026). https://pith.science/paper/53XSQA62
@misc{pith2026250105498,
author = {Pith},
title = {Pith review of: Generative Flow Networks: Theory and Applications to Structure Learning},
year = {2026},
howpublished = {\url{https://pith.science/paper/53XSQA62}},
note = {Machine review of arXiv:2501.05498}
}
read the original abstract
Without any assumptions about data generation, multiple causal models may explain our observations equally well. To avoid selecting a single arbitrary model that could result in unsafe decisions if it does not match reality, it is therefore essential to maintain a notion of epistemic uncertainty about our possible candidates. This thesis studies the problem of structure learning from a Bayesian perspective, approximating the posterior distribution over the structure of a causal model, represented as a directed acyclic graph (DAG), given data. It introduces Generative Flow Networks (GFlowNets), a novel class of probabilistic models designed for modeling distributions over discrete and compositional objects such as graphs. They treat generation as a sequential decision making problem, constructing samples of a target distribution defined up to a normalization constant piece by piece. In the first part of this thesis, we present the mathematical foundations of GFlowNets, their connections to existing domains of machine learning and statistics such as variational inference and reinforcement learning, and their extensions beyond discrete problems. In the second part of this thesis, we show how GFlowNets can approximate the posterior distribution over DAG structures of causal Bayesian Networks, along with the parameters of its causal mechanisms, given observational and experimental data.
Figures
Figures from the paper (93 more)
Forward citations
Cited by 2 Pith papers
-
Revisiting Non-Acyclic GFlowNets in Discrete Environments
In cyclic discrete environments, GFlowNet flows are expected visit counts, and training a non-acyclic GFlowNet with the smallest expected trajectory length is equivalent to minimizing total flow.
-
Efficient Diversity-Preserving Diffusion Alignment via Gradient-Informed GFlowNets
A gradient-informed GFlowNet objective, residual nabla-DB, finetunes diffusion models to sample according to a reward while preserving diversity and prior knowledge.
Reference graph
Works this paper leans on
-
[1]
Raj Agrawal, Chandler Squires, Karren Yang, Karthikeyan Shanmugam, and Caroline Uhler (2019). ABCD-Strategy: Budgeted Experimental Design for Targeted Causal Structure Discovery. International Conference on Artificial Intelligence and Statistics (AISTATS)(cited on pp. 193, 200). Hirotugu Akaike (1974). A New Look at the Statistical Model Identification.IE...
work page 2019
-
[3]
This concludes the proof, showing that invariance also holds ats0
Therefore ∑ s∈S F (s)PF (s0|s) = ∑ s∈S ∞∑ n=0 F (s0)Es0 [ 1(n<σ s0)1(Xn =s) ] PF (s0|s) (A.2.5) =F (s0) ∑ s∈S ∞∑ n=1 Es0 [ 1(σs0 =n)1(Xn−1 =s) ] (A.2.6) =F (s0) ∞∑ n=1 Es0 [ 1(σs0 =n) ∑ s∈S 1(Xn−1 =s) /bracehtipupleft/bracehtipdownright/bracehtipdownleft/bracehtipupright = 1 ] (A.2.7) =F (s0) ∞∑ n=1 Es0 [ 1(σs0 =n) ] =F (s0), (A.2.8) where we substitutedF...
-
[2019]
Proposition A.1.2 (Data processing inequality)
and recalled for completeness, plays an important role in proving Proposition 4.3.1, which establishes a bound on thef-divergence between the terminating state probability distribution of a GFlowNet and the target Gibbs distribution, in terms of the f-divergence between distributions over complete trajectories (Malkin et al., 2022). Proposition A.1.2 (Dat...
work page 2022
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.