GFlowNet Foundations

Abstract

Generative Flow Networks (GFlowNets) have been introduced as a method to sample a diverse set of candidates in an active learning context, with a training objective that makes them approximately sample in proportion to a given reward function. In this paper, we show a number of additional theoretical properties of GFlowNets, including a new local and efficient training objective called detailed balance for the analogy with MCMC. GFlowNets can be used to estimate joint probability distributions and the corresponding marginal distributions where some variables are unspecified and, of particular interest, can represent distributions over composite objects like sets and graphs. GFlowNets amortize the work typically done by computationally expensive MCMC methods in a single but trained generative pass. They could also be used to estimate partition functions and free energies, conditional probabilities of supersets (supergraphs) given a subset (subgraph), as well as marginal distributions over all supersets (supergraphs) of a given set (graph). We introduce variations enabling the estimation of entropy and mutual information, sampling from a Pareto frontier, connections to reward-maximizing policies, and extensions to stochastic environments, continuous actions and modular energy functions.

1 Introduction

Building upon the introduction of Generative Flow Networks (GFlowNets) by Bengio et al. (2021), we provide here an in-depth formal foundation and expansion of the set of theoretical results in ways that may be of interest for the active learning scenario of Bengio et al. (2021) but also much more broadly.

1.1 What is a GFlowNet ?

GFlowNets have properties which make them well-suited to perform amortized probabilistic inference in general, whether for sampling or for marginalizing. Sampling takes place at training time while run-time sampling or computations of marginalized quantities can be done in a single pass through a sequence of constructive stochastic steps. This makes GFlowNets an interesting alternative to Monte-Carlo Markov chains (MCMC) and related to amortized variational inference (Malkin et al., 2023).

Because sampling of a compositional object $s$ can be achieved through a sequence of stochastic steps, very rich multimodal distributions $P_T(s)$ over such objects can be represented, and the offline training objectives make it possible to explore and discover modes of the distribution of interest. The key property of GFlowNets is that their sampling policy is trained to make the probability $P_T(s)$ of sampling an object $s$ approximately proportional to the value $R(s)$ of a given reward function applied to that object. We also talk of an energy function $ℰ(s)=-\log R(s)$, i.e., the reward function is non-negative and corresponds to an unnormalized probability. Whereas one typically trains a generative model from a dataset of positive examples, a GFlowNet is trained to match the given energy or reward function and convert it into a sampler. We view that sampler as a generative policy because the composite object $s$ is constructed through a sequence of smaller stochastic steps (see Fig. 1), often corresponding to constructively composing different elements of $s$, like the edges of a graph.

This conversion of an energy function or unnormalized probability function to a sampler is similar to what MCMC methods achieve but once trained, GFlowNets will generate a sample in one shot instead of generating a long sequence of samples whose distribution would gradually approach the desired one. GFlowNets thus avoid the lengthy stochastic search in the space of such objects and the associated mode-mixing intractability challenge of MCMC methods (Jasra et al., 2005; Bengio et al., 2013; Pompe et al., 2020). Multiple iid samples can be obtained from the GFlowNet by calling the sampler multiple times. GFlowNets exchange that intractability of sampling with MCMC for the challenge of amortized training of the generative policy. The latter problem would be equally intractable if the modes of the reward function did not have a inherent (but not necessarily known) structure over which the learner could generalize, i.e., the learner had almost no chance to correctly guess where to find new modes based on (i.e., training on) those it had already visited.

Figure 1: A diagram of how a GFlowNet iteratively constructs an object. We adopt notation that is common in the reinforcement learning literature: $s_t$ represents the state of the partially constructed object (in this case, a graph) at time $t$, $a_t$ represents the action taken by the GFlowNet at time $t$ to transition to state $s_{t+1}=T(s_t,a_t)$. In this diagram, the GFlowNet takes a 3-node graph as input and determines an action to take. The action, combined with the environment transition function $T(s_t,a_t)$, determines $s_{t+1}$: a four-node graph. This process repeats until an exit action is sampled and the sample is complete.

The energy function or reward function (exponential of minus energy) is evaluated only at the end of the sequential construction process for objects $s$, in what we call a terminating state. Every such constructive sequence starts in the single initial state $s_0$ and ends in a terminal state. As illustrated in Figure 2, we can visualize the set of all trajectories starting from $s_0$ and ending in a terminal state $s$. The term ”flow” in ”generative flow networks” refers to unnormalized probabilities that can be learned by GFlowNet learning procedures. The flow in an intermediate state $s$ is a weighted sum of the non-negative rewards of the terminating states reachable from $s$. Those weights are such as to avoid double-counting: if we were to inject a fixed flow of liquid in $s_0$ and dispatch that liquid in each child of any state $s$ proportionally to the GFlowNet policy for choosing a child of $s$, we would obtain the flow at each state and the flow at terminating states would match the reward function at those states. As shown in greater detail here and for the first time in the first GFlowNet paper (Bengio et al., 2021), this can be achieved with a flow constraint at each state: the sum of incoming flows must match the sum of outgoing flows.

Figure 2: Illustration of the structure of a Generative Flow Network (GFlowNet), as a pointed DAG over states $s$, with particles flowing along edges to represent the flow function. Any object sampled by the GFlowNet policy can be obtained by starting from initial state $s_0$ and then at each step choosing a child with probability proportional to the GFlowNet policy’s transition probability. This process stops when a terminating action is chosen from a terminating state $s$ (yielding a terminal state $s_f$), at which point a reward $R(s)$ is obtained. The figure shows a tiny GFlowNet and the possible trajectories from $s_0$ to any of the terminal states. It illustrates that in general a state can be reached through several trajectories. GFlowNet algorithms learn a policy such that the probability of sampling terminating state $s$ is proportional to $R(s)$. It tries to learn a flow function $F(s)$ and $F(s\to s^{\prime })$ over all states (including intermediate states) $s$ and transitions $s\to s^{\prime }$ with $F(s)=R(s)$ at terminal states and $F(s_0)$ being the sum of rewards over all terminal states. A sufficient property to achieve this is that at each state the sum of incoming flows equals the sum of outgoing flows.

1.2 Contributions of this paper

In this paper, an important contribution is the notion of conditional GFlowNet, which enables estimation of intractable sums corresponding to marginalization over many steps of object construction, and can thus be used to compute free energies over different types of joint distributions, perhaps most interestingly over sets and graphs. This marginalization also enables estimation of entropies, conditional entropies  and mutual information. GFlowNets can thus be generalized to estimate multiple flows corresponding to modeling a rich outcome (rather than a scalar reward function) .

We refer the reader to Bengio et al. (2021) and Sec. 7 for a discussion of related approaches and differences with common generative models and reinforcement learning (RL) methods. In an RL context, two interesting properties of GFlowNets already noted in that paper are that they (1) can be trained in an offline manner with trajectories sampled from a distribution different from the one represented by the GFlowNet and (2) they match the reward function in probability rather than try to find a configuration which maximizes rewards or returns. The latter property is particularly interesting in the context of exploration, to ensure the configurations sampled from the generative policy are both interesting and diverse. It is also interesting to transform GFlowNets into amortized probabilistic inference machines: if we choose the reward function to be a prior (over some random variable) times a likelihood (how well some data is fit given that choice of random variable value), then the GFlowNet policy learns to sample from the corresponding Bayesian posterior (which is proportional to prior times likelihood). The ability of GFlowNets to generate a diverse set of samples then corresponds to the ability to sample from the modes of the target distribution.

An important source of inspiration for GFlowNets is the way information propagates in temporal-difference RL methods (Sutton and Barto, 2018). Both rely on a principle of coherence for credit assignment which may only be achieved asymptotically when training converges. While exact gradient calculation may be intractable, because the number of paths in state space to consider is exponentially large, both methods rely on local coherence between different components and a training objective that states that if all the learned components are coherent with each other locally, then we obtain a system that estimates the quantities of interest globally. Examples include estimation of expected discounted returns in temporal-difference methods and probability measures with GFlowNets.

This paper extends the theory of the original GFlowNet construction (Bengio et al., 2021) in several directions, including a new local training objective called detailed balance (for the analogy with the detailed balance condition of Monte-Carlo Markov chains) which avoids forming explicit sums required by the previously proposed flow matching loss, as well as formulations enabling the calculation of marginal probabilities (or free energies) for subsets of variables, more generally for subsets of larger sets, or subgraphs, their application to estimating entropy and mutual information, and the introduction of an unsupervised form of GFlowNets (the reward function is not needed while training, only observations of outcomes) enabling sampling from a Pareto frontier, for example. Although basic GFlowNets are more similar to bandits (in that a reward is only provided at the end of a sequence of actions), they can be extended to take into account intermediate rewards and thus a notion of return, and sample according to these returns. The original formulation of GFlowNets is also limited to discrete and deterministic environments, while this paper suggests how these two limitations could be lifted. Finally, whereas the basic formulation of GFlowNets assumes a given reward or energy function, this paper considers how the energy function could be jointly learned with the GFlowNet, opening the door to novel energy-based modeling methodologies and a modular structure for both the energy function and the GFlowNet.

1.3 GFlowNets in other works

In addition to the theory presented in this paper, Malkin et al. (2023) and Zimmermann et al. (2022) prove some partial equivalences between GFlowNets and hierarchical variational methods, providing yet more theoretical evidence for the efficacy of GFlowNets in learning to sample proportionally to a given reward function. These works also provide evidence for the superiority of GFlowNets in off-policy settings.

GFlowNets have found a wide array of applications due to the associated diversity of generated samples. In contexts where a cheap proxy for the true reward function exists, GFlowNets have been used to surface samples under which to query the proxy before more expensive evaluation under the true reward function. In these settings, the diversity of samples generated by GFlowNets can be used for robustness to proxy misspecification and to incorporate epistemic uncertainty. For example, Zhang et al. (2023) use GFlowNets to produce sample schedules for operations in a computation graph, where evaluating the runtimes of sample schedules via a proxy is fast but evaluating the same schedules on target hardware is expensive. In active learning problems, Jain et al. (2022, 2023) use GFlowNet sampling as a subroutine inside an active learning loop as a substitute for Bayesian Optimization or RL-based methods. Jain et al. (2022) apply GFlowNets to search for novel anti-microbial peptides, discover DNA sequences that have high binding activity with human transcription factors, and to find proteins with high fluorescence. Additionally, Jain et al. (2023) develops preference-conditional GFlowNets, where a preference weight vector is used to scalarize multiple objective functions into a single reward. The authors apply their techniques to various molecule and DNA sequence generation tasks and find that their methods are able to find different Pareto-optimal samples along the Pareto frontier.

GFlowNets have found applications in several other machine learning problems. For example, Zhang et al. (2022) simultaneously train an energy-based model and a GFlowNet; the energy function is trained with samples from a GFlowNet, which, in turn, uses the energy function to form its reward. Their method results in a generative model for binary vectors in high dimensions, e.g., binarized digits. Deleu et al. (2022) use a GFlowNet for structure learning; the GFlowNet produces samples that approximates the true posterior over causal graphs given a dataset. Their method works on both observational and interventional data, and compares favorably to MCMC- and variational inference-based methods. Hu et al. (2023) find maximum-likelihood estimates of latent variable models with discrete compositional latents by jointly training a GFlowNet to approximately sample from the generally intractable posterior in the E-step of the expectation-maximization (EM) algorithm.

2 Flow Networks and Markovian Flows

2.1 Some elements of graph theory

In this section, we recall some basic definitions and properties of graphs, which are the basis of flow networks and GFlowNets.

Definition 1.

A directed graph is a tuple $G=(𝒮,𝔸)$, where $𝒮$ is a finite set of states, and $𝔸$ a subset of $𝒮\times 𝒮$ representing directed edges. Elements of $𝔸$ are denoted $s\to s^{\prime }$ and called edges or transitions.
A trajectory in such a graph is a sequence $\tau =(s_1,\mathrm{\dots },s_n)$ of elements of $𝒮$ such that every transition $s_t\to s_{t+1}\in 𝔸$ and $n>1$. We denote $s\in \tau$ to mean that $s$ is in the trajectory $\tau$, i.e., $\exists t\in \{1,\mathrm{\dots },n\}{s}_t=s$, and similarly $s\to s^{\prime }\in \tau$ to mean that $\exists t\in \{1,\mathrm{\dots },n-1\}{s}_t=s,s_{t+1}=s^{\prime }$. For convenience, we also use the notation $\tau =s_1\to \mathrm{\dots }\to s_n$. The length of a trajectory is the number of edges in it (the length of $\tau =(s_1,\mathrm{\dots },s_n)$ is thus $n-1$).
A directed acyclic graph (DAG) is a directed graph in which there is no trajectory $\tau =(s_1,\mathrm{\dots },s_n)$ satisfying $s_n=s_1$.

Given a DAG $G=(𝒮,𝔸)$, and two states $s,s^{\prime }\in 𝒮$, if there exists a trajectory in $G$ starting in $s$ and ending in $s^{\prime }$, then we write . The binary relationship “” defines a strict partial order (i.e. it is irreflexive, asymmetric and transitive). We write $s\le s^{\prime }$ if or $s=s^{\prime }$. The binary relation “$\le$” is a (non-strict) partial order (i.e. it is reflexive, antisymmetric and transitive).

If there is no order relation between $s$ and $s^{\prime }$, we write $s\lessgtr s^{\prime }$.

Definition 2.

Given a DAG $G=(𝒮,𝔸)$, the parent set of a state $s\in 𝒮$, which we denote $Par(s)$, contains all of the direct parents of $s$ in $G$, i.e., $Par(s)=\{s^{\prime }\in 𝒮:s^{\prime }\to s\in 𝔸\}$; similarly, the child set $Child(s)$ contains all of the direct children of $s$ in $G$, i.e., $Child(s)=\{s^{\prime }\in 𝒮:s\to s^{\prime }\in 𝔸\}$.

Definition 3.

Given a DAG $G=(𝒮,𝔸)$. $G$ is called a pointed DAG if there exist two states $s_0,s_f\in 𝒮$ that satisfy:

$s_0$ is called the source state or initial state. $s_f$ is called the sink state or final state. Because “” is a strict partial order, these two states are unique.

A complete trajectory in such a DAG is any trajectory starting in $s_0$ and ending in $s_f$. We denote such a trajectory as $\tau =(s_0,s_1,\mathrm{\dots },s_n,s_{n+1}=s_f)$.

We denote by $𝒯$ the set of all complete trajectories in $G$, and by ${𝒯}^{partial}$ the set of (possibly incomplete) trajectories in $G$.

A state $s\in 𝒮$ is called a terminating state if it is a parent of the sink state, i.e. $s\to s_f\in 𝔸$. The transition $s\to s_f$ is called a terminating edge. We denote by:

• ${𝔸}^{-f}=\{s\to s^{\prime }\in 𝔸,s^{\prime }\ne s_f\}$, the set of non-terminating edges in $G$,

• ${𝔸}^f=\{s\to s^{\prime }\in 𝔸,s^{\prime }=s_f\}=𝔸\setminus {𝔸}^{-f}$, the set of terminating edges in $G$,

• ${𝒮}^f=\{s\in 𝒮,s\to s_f\in {𝔸}^f\}=Par(s_f)$, the set of terminating states in $G$.

In Fig. 3, we visualize the concepts introduced in the previous definitions.

Figure 3: Example of a pointed DAG $G$ illustrating the notions of initial state ($s_0$), final or sink state ($s_f$), terminating states in ${𝒮}^f$, with a transition to $s_f$ called a terminating edge, in ${𝔸}^f$. A terminating state may have other children different from the sink state (e.g., the terminating state $s_7$).

Note that the constraint of a single source state and single sink state is only a mathematical convenience since a bijection exists between general DAGs and those with this constraint (by the addition of a unique source/sink state connected to all the other source/sink states).

Definition 4.

Let $G$ be a pointed DAG with source state $s_0$ and sink state $s_f$. A forward (resp. backward) probability function consistent with $G$ is any non-negative function ${\widehat{P}}_F$ (resp. ${\widehat{P}}_B$) defined on $𝔸$ that satisfies $\forall s\in 𝒮\setminus \{s_f\},{\sum }_{s^{\prime }\in Child(s)}{\widehat{P}}_F(s^{\prime }\mid s)=1$ (resp. $\forall s\in 𝒮\setminus \{s_0\},{\sum }_{s^{\prime }\in Par(s)}{\widehat{P}}_B(s^{\prime }\mid s)=1$).

With pointed DAGs, consistent forward and backward probability functions, that are probabilities over states, can be used to define probabilities over trajectories, i.e. probability measures on some subsets of ${𝒯}^{partial}$. The following lemma shows how to construct such factorized probability measures:

Lemma 5.

Let $G=(𝒮,𝔸)$ be a pointed DAG, and consider a forward probability function $\widehat{P_F}$, and a backward probability function $\widehat{P_B}$ both consistent with $G$. For any state $s\in 𝒮\setminus \{s_f\}$, we denote by ${𝒯}_{s,f}\subseteq {𝒯}^{partial}$ the set of trajectories in $G$ starting in $s$ and ending in $s_f$; and for any state $s\in 𝒮\setminus \{s_0\}$, we denote by ${𝒯}_{0,s}\subseteq {𝒯}^{partial}$ the set of trajectories in $G$ starting in $s_0$ and ending in $s$.

Consider the extensions of ${\widehat{P}}_F$ and ${\widehat{P}}_B$ on ${𝒯}^{partial}$ defined by:

$\forall \tau =(s_1,\mathrm{\dots },s_n)\in {𝒯}^{partial}{\widehat{P}}_F(\tau ):={\displaystyle \prod _{t=1}^{n-1}}{\widehat{P}}_F(s_{t+1}\mid s_t)$ (1)

$\forall \tau =(s_1,\mathrm{\dots },s_n)\in {𝒯}^{partial}{\widehat{P}}_B(\tau ):={\displaystyle \prod _{t=1}^{n-1}}{\widehat{P}}_B(s_t\mid s_{t+1})$ (2)

We have the following:

$\forall s\in 𝒮\setminus \{s_f\}{\displaystyle \sum _{\tau \in {𝒯}_{s,f}}}{\widehat{P}}_F(\tau )=1$ (3)

$\forall s^{\prime }\in 𝒮\setminus \{s_0\}{\displaystyle \sum _{\tau \in {𝒯}_{0,s^{\prime }}}}{\widehat{P}}_B(\tau )=1$ (4)

Proof For convenience, we will use ${𝒯}_{s\to s^{\prime },s_f}$ to denote the set of trajectories starting with $s\to s^{\prime }$ and ending in $s_f$, and ${𝒯}_{0,s\to s^{\prime }}$ to denote the set of trajectories starting in $s_0$ and ending with $s\to s^{\prime }$. This allows to write:

$\forall s\ne s_f{𝒯}_{s,f}={\displaystyle \bigcup _{s^{\prime }\in Child(s)}}{𝒯}_{s\to s^{\prime },s_f},\{{𝒯}_{s\to s^{\prime },s_f},s^{\prime }\in Child(s)\}\text{pairwise disjoint},$

$\forall s^{\prime }\ne s_0{𝒯}_{0,s^{\prime }}={\displaystyle \bigcup _{s\in Par(s^{\prime })}}{𝒯}_{0,s\to s^{\prime }},\{{𝒯}_{0,s\to s^{\prime }},s\in Par(s^{\prime })\}\text{pairwise disjoint}.$

Additionally, for any $s\ne s_f$, we denote by $d_{s,f}$ the maximum trajectory length in ${𝒯}_{s,f}$; and for any $s^{\prime }\ne s_0$, we denote by $d_{0,s^{\prime }}$ the maximum trajectory length in ${𝒯}_{0,s}$.

We will prove Eq. 3 by strong induction on $d_{s,f}$ and Eq. 4 by strong induction on $d_{0,s^{\prime }}$.

Base cases: If $d_{s,f}=1$ and $d_{0,s^{\prime }}=1$, then ${𝒯}_{s,f}=\{(s\to s_f)\}$ and ${𝒯}_{0,s^{\prime }}=\{(s_0\to s^{\prime })\}$. Hence, ${\sum }_{\tau \in {𝒯}_{s,f}}{\widehat{P}}_F(\tau )={\widehat{P}}_F(s\to s_f)={\widehat{P}}_F(s_f\mid s)=1$ given that $s_f$ is the only child of $s$ (otherwise $d_{s,f}$ cannot be $1$), and ${\sum }_{\tau \in {𝒯}_{0,s^{\prime }}}{\widehat{P}}_B(\tau )={\widehat{P}}_B(s_0\mid s^{\prime })=1$ given that $s_0$ is the only parent of $s^{\prime }$ (otherwise $d_{0,s^{\prime }}$ cannot be $1$).

Induction steps: Consider $s\ne s_f$ such that $d_{s,f}>1$ and $s^{\prime }\ne s_0$ such that $d_{0,s^{\prime }}>1$. Because of the disjoint unions written above, we have:

${\displaystyle \sum _{\tau \in {𝒯}_{s,f}}}{\widehat{P}}_F(\tau )={\displaystyle \sum _{\tilde{s}\in Child(s)}}{\displaystyle \sum _{\tau \in {𝒯}_{s\to \tilde{s},f}}}{\widehat{P}}_F(\tau )={\displaystyle \sum _{\tilde{s}\in Child(s)}}{\widehat{P}}_F(\tilde{s}\mid s){\displaystyle \sum _{\tau \in {𝒯}_{\tilde{s},f}}}{\widehat{P}}_F(\tau )=1,$

${\displaystyle \sum _{\tau \in {𝒯}_{0,s^{\prime }}}}{\widehat{P}}_B(\tau )={\displaystyle \sum _{{\tilde{s}}^{\prime }\in Par(s^{\prime })}}{\displaystyle \sum _{\tau \in {𝒯}_{0,{\tilde{s}}^{\prime }\to s}}}{\widehat{P}}_B(\tau )={\displaystyle \sum _{{\tilde{s}}^{\prime }\in Par(s^{\prime })}}{\widehat{P}}_B({\tilde{s}}^{\prime }\mid s^{\prime }){\displaystyle \sum _{\tau \in {𝒯}_{0,{\tilde{s}}^{\prime }}}}{\widehat{P}}_B(\tau )=1,$

where we used the induction hypotheses in the third equality of each line.

2.2 Trajectories and Flows

We augment pointed DAGs it with a function $F$ called a flow. An analogy which helps to picture flows is a stream of particles flowing through a network where each particle starts at $s_0$ and flowing through some trajectory terminating in $s_f$. The flow $F(\tau )$ associated with each complete trajectory $\tau$ contains the number of particles sharing the same path $\tau$.

Definition 6.

Given a pointed DAG, a trajectory flow (or “flow”) is any non-negative function $F:𝒯\mapsto {ℝ}^{+}$ defined on the set of complete trajectories $𝒯$. $F$ induces a measure over the $\sigma$-algebra $\mathrm{\Sigma }=2^{𝒯}$, the power set on the set of complete trajectories $𝒯$. In particular, for every subset $A\subseteq 𝒯$, we have

$$F(A)=\sum _{\tau \in A}F(\tau ).$$ (5)

The pair $(G,F)$ is called a flow network.

This definition ensures that $(𝒯,2^{𝒯},F)$ is a measure space. We abuse the notation here, using $F$ to denote both a function of complete trajectories, and its corresponding measure over $(𝒯,2^{𝒯})$. A special case is when the event $A$ is the singleton trajectory $\{\tau \}$, where we just write its measure as $F(\tau )$. We also abuse the notation to define the flow through either a particular state $s$, or through a particular edge $s\to s^{\prime }$ in the following way.

Definition 7.

The flow through a state (or state flow) $F:𝒮\mapsto {ℝ}^{+}$ corresponds to the measure of the set of complete trajectories going through a particular state:

$$F(s):=F(\{\tau \in 𝒯:s\in \tau \})=\sum _{\tau \in 𝒯:s\in \tau }F(\tau ).$$ (6)

Similarly, the flow through an edge (or edge flow) $F:𝔸\mapsto {ℝ}^{+}$ corresponds to the measure of the set of complete trajectories going through a particular edge:

$$F(s\to s^{\prime }):=F(\{\tau \in 𝒯:s\to s^{\prime }\in \tau \})=\sum _{\tau \in 𝒯:s\to s^{\prime }\in \tau }F(\tau ).$$ (7)

Note that with this definition, we have $F(s\to s^{\prime })=0$ if $s\to s^{\prime }\notin 𝔸$ is not an edge in the pointed DAG (since $F(\mathrm{\varnothing })=0$). We call the flow of a terminating transition $F(s\to s_f)$ a terminating flow. The following proposition relates the state flows and the edge flows:

Proposition 8.

Given a flow network $(G,F)$. The state flows and edge flows satisfy:

$\forall s\in 𝒮\setminus \{s_f\}F(s)={\displaystyle \sum _{s^{\prime }\in Child(s)}}F(s\to s^{\prime })$ (8)

$\forall s^{\prime }\in 𝒮\setminus \{s_0\}F(s^{\prime })={\displaystyle \sum _{s\in Par(s^{\prime })}}F(s\to s^{\prime })$ (9)

Proof Given $s\ne s_f$, the set of complete trajectories going through $s$ is the (disjoint) union of the sets of trajectories going through $s\to s^{\prime }$, for all $s^{\prime }\in Child(s)$:

$$\{\tau \in 𝒯:s\in \tau \}=\bigcup _{s^{\prime }\in Child(s)}\{\tau \in 𝒯:s\to s^{\prime }\in \tau \}.$$

Therefore, it follows that:

$$F(s)=\sum _{\tau :s\in \tau }F(\tau )=\sum _{s^{\prime }\in Child(s)}\sum _{\tau :s\to s^{\prime }\in \tau }F(\tau )=\sum _{s^{\prime }\in Child(s)}F(s\to s^{\prime })$$

Similarly, Eq. 9 follows by writing the set of complete trajectories going though $s^{\prime }\ne s_0$ as the (disjoint) union of the sets of trajectories going through $s\to s^{\prime }$ for all $s\in Par(s^{\prime })$.

2.3 Flow Induced Probability Measures

Definition 9.

Given a flow network $(G,F)$, the total flow $Z$ is the measure of the whole set $𝒯$, corresponding to the sum of the flows of all the complete trajectories:

$$Z:=F(𝒯)=\sum _{\tau \in 𝒯}F(\tau ).$$ (10)

Proposition 10.

The flow through the initial state equals the flow through the final state equals the total flow $Z$.

Proof Since $\forall \tau \in 𝒯,s_0,s_f\in \tau$, applying Eq. 6 to $s_0$ and $s_f$ yields

$F(s_0)$ $={\displaystyle \sum _{\tau \in 𝒯}}F(\tau )=Z,$ (11)

$F(s_f)$ $={\displaystyle \sum _{\tau \in 𝒯}}F(\tau )=Z.$ (12)

Intuitively, Prop. 10 justifies the use of the term “flow”, introduced by Bengio et al. (2021), by analogy with a stream of particles flowing from the initial state to the final states.

We use the letter $Z$ in Def. 9, often used to denote the partition function in probabilistic models and statistical mechanics, because it is a normalizing constant which can turn the measure space $(𝒯,2^{𝒯},F)$ defined above into the probability space $(𝒯,2^{𝒯},P)$:

Definition 11.

Given a flow network $(G,F)$, the flow probability is the probability measure $P$ over the measurable space $(𝒯,2^{𝒯})$ associated with $F$:

$$\forall A\subseteq 𝒯P(A):=\frac{F(A)}{F(𝒯)}=\frac{F(A)}{Z}.$$ (13)

For two events $A,B\subseteq 𝒯$, the conditional probability $P(A\mid B)$ thus satisfies:

$P(A\mid B)$ $:={\displaystyle \frac{F(A\cap B)}{F(B)}}.$ (14)

Similar to the flow $F$, we abuse the notation $P$ to define the probability of going through a state:

$$\forall s\in 𝒮P(s):=\frac{F(s)}{Z},$$ (15)

and similarly for the probability of going through an edge. Note that $P(s)$ does not correspond to a distribution over states, in the sense that ${\sum }_{s\in 𝒮}P(s)\ne 1$; in particular, it is easy to see that $P(s_0)=1$ (in other words, the probability of a trajectory passing through the initial state $s_0$ is $1$). Additionally, for a trajectory $\tau \in 𝒯$, we also use the abuse of notation $P(\tau )$ instead of $P(\{\tau \})$ to denote the probability of going through a specific trajectory $\tau$.

Definition 12.

Given a flow network $(G,F)$, the forward transition probability operator $P_F$ is a function on $𝒮\times 𝒮$, that is a special case of the conditional probabilities induced by $F$ (Eq. 14):

$$\forall s\to s^{\prime }\in 𝔸P_F(s^{\prime }\mid s):=P(s\to s^{\prime }\mid s)=\frac{F(s\to s^{\prime })}{F(s)}.$$ (16)

Similarly, the backwards transition probability is the operator defined by:

$$\forall s\to s^{\prime }\in 𝔸P_B(s\mid s^{\prime }):=P(s\to s^{\prime }\mid s^{\prime })=\frac{F(s\to s^{\prime })}{F(s^{\prime })}.$$ (17)

Note how $P_F$ and $P_B$ are consistent with $G$ (in the sense of Def. 4), as a consequence of Prop. 8.

Because flows define probabilities over states and edges, they can be used to define probability distributions over the terminating states of a graph (denoted by ${𝒮}^f=Par(s_f)$) as follows:

Definition 13.

Given a flow network $(G,F)$, the terminating state probability $P_T$ is the probability over terminating states ${𝒮}^f$ under the flow probability $P$:

$$\forall s\in {𝒮}^f{P}_T(s):=P(s\to s_f)=\frac{F(s\to s_f)}{Z}$$ (18)

Contrary to the probability $P(s)$ of going through a state $s$, the terminating state probability $P_T$ is a well-defined distribution over the terminating states $s\in {𝒮}^f$, in the following sense:

Proposition 14.

The terminating state probability $P_T$ is a well-defined distribution over the terminating states $s\in {𝒮}^f$, in that $P_T(s)\ge 0$ for all $s\in {𝒮}^f$, and

$$\sum _{s\in {𝒮}^f}{P}_T(s)=1.$$

Proof Since the flow $F(s\to s_f)$ is non-negative, it is easy to see that $P_T(s)\ge 0$. Moreover, using the definition of ${𝒮}^f=Par(s_f)$, Prop. 8 (relating the edge flows and the state flows), and Prop. 10 ($F(s_f)=Z$), we have

$$\sum _{s\in {𝒮}^f}{P}_T(s)=\frac{1}{Z}\sum _{s\in {𝒮}^f}F(s\to s_f)=\frac{1}{Z}\sum _{s\in Par(s_f)}F(s\to s_f)=\frac{F(s_f)}{Z}=1.$$

The terminating state probability is particularly important in the context of estimating flow networks (see Sec. 3), as it shows that a flow network $(G,F)$ induces a probability distribution over terminating states which is proportional to the terminating flows $F(s\to s_f)$, the normalization constant $Z$ being given by initial flow $F(s_0)$.

2.4 Markovian Flows

Defining a flow requires the specification of $|𝒯|$ non-negative values (one for every trajectory $\tau \in 𝒯$), which is generally exponential in the number of graph edges. Markovian flows however have the remarkable property that they can be defined with much fewer “numbers”, given that trajectory flows factorize according to $G$.

Definition 15.

Let $(G,F)$ be a flow network, with flow probability measure $P$. $F$ is called a Markovian flow (or equivalently $(G,F)$ a Markovian flow network) if, for any state $s\ne s_0$, outgoing edge $s\to s^{\prime }$, and for any trajectory $\tau =(s_0,s_1,\mathrm{\dots },s_n=s)\in {𝒯}^{partial}$ starting in $s_0$ and ending in $s$:

$$P(s\to s^{\prime }\mid \tau )=P(s\to s^{\prime }\mid s)=P_F(s^{\prime }\mid s).$$ (19)

Note that the Markovian property does not hold for all of the flows as defined in the previous sections (e.g. Fig. 4). Intuitively, a flow can be considered non-Markovian if a particle in the “flow stream” can remember its past history; if not, its future behavior can only depend on its current state and the flow must be Markovian. In this work, we will primarily be concerned with Markovian flows, though later we will re-introduce a form of memory via state-conditional flows that allow each flow “particle” to remember parts of its history. The following proposition shows that Markovian flows have the property that the flows at (or the probabilities of) complete trajectories factorize according the the graph, and that it is a sufficient condition for defining Markovian flows.

Proposition 16.

Let $(G,F)$ be a flow network, and $P$ the corresponding flow probability. The following three statements are equivalent:

1. $F$ is a Markovian flow

2. There exists a unique probability function ${\widehat{P}}_F$ consistent with $G$ such that for all complete trajectories $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)$:

   $$P(\tau )=\prod _{t=1}^{n+1}{\widehat{P}}_F(s_t\mid s_{t-1}).$$ (20)

   Moreover, the probability function ${\widehat{P}}_F$ is exactly the forward transition probability associated with the flow probability $P$: ${\widehat{P}}_F=P_F$.

3. There exists a unique probability function ${\widehat{P}}_B$ consistent with $G$ such that for all complete trajectories $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)$:

   $$P(\tau )=\prod _{t=1}^{n+1}{\widehat{P}}_B(s_{t-1}\mid s_t).$$ (21)

   Moreover, the probability function ${\widehat{P}}_B$ is exactly the backwards transition probability associated with the flow probability $P$: ${\widehat{P}}_B=P_B$.

Proof Recall from Lemma 5 the notations ${𝒯}_{0,s}$ to denote the set of partial trajectories from $s_0$ to $s$, and ${𝒯}_{s^{\prime },f}$ to denote the set of partial trajectories from $s^{\prime }$ to $s_f$. We will prove the equivalences $1\iff 2$ and $1\iff 3$.

• $1\Rightarrow 2$: Suppose that $F$ is a Markovian flow. Then using the laws of probability, the Markov property in Eq. 19, and $P(s_0)=1$, for some complete trajectory $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)$:

  $P(\tau )$ $=P(s_0\to s_1\to \mathrm{\dots }\to s_{n+1})=P(s_0\to s_1){\displaystyle \prod _{t=1}^n}P(s_t\to s_{t+1}\mid s_0\to \mathrm{\dots }\to s_t)$

  $=P(s_0\to s_1){\displaystyle \prod _{t=1}^n}P(s_t\to s_{t+1}\mid s_t)$

  $=P(s_0)P_F(s_1\mid s_0){\displaystyle \prod _{t=1}^n}{P}_F(s_{t+1}\mid s_t)$

  $={\displaystyle \prod _{t=1}^{n+1}}{P}_F(s_t\mid s_{t-1}),$

  where the second line uses to Markov property, and the third line uses the definition of the forward transition probability $P_F$. $P_F$ thus satisfies Eq. 20 for all complete trajectories.

  To show uniqueness of $P_F$, assume Eq. 20 is satisfied by some ${\widehat{P}}_F$ for all complete trajectories. By definition of the forward transition probability:

  $$P_F(s^{\prime }\mid s):=P(s\to s^{\prime }\mid s)=\frac{P(s\to s^{\prime })}{P(s)}.$$

  Any complete trajectory $\tau$ going through a state $s$ can be (uniquely) decomposed into a partial trajectory ${\tau }^{\prime }\in {𝒯}_{0,s}$ from $s_0$ to $s$, and a partial trajectory ${\tau }^{\prime \prime }\in {𝒯}_{s,f}$ from $s$ to $s_f$. Using the definition of $P(s)$, we have:

  $P(s)$ $={\displaystyle \sum _{\tau :s\in \tau }}P(\tau )={\displaystyle \sum _{\tau :s\in \tau }}{\displaystyle \prod _{(s_t\to s_{t+1})\in \tau }}{\widehat{P}}_F(s_{t+1}\mid s_t)$

  $=\left[{\displaystyle \sum _{{\tau }^{\prime }\in {𝒯}_{0,s}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime }}}{\widehat{P}}_F(s_{t+1}\mid s_t)\right]\underset{=\mathrm{\hspace{0.17em}1}\text{(}\text{<a class="ltx\_ref" href="\#Thmtheorem5" style="font-size:70\%;" title="Lemma 5. ‣ 2.1 Some elements of graph theory ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Lemma</span> <span class="ltx\_text ltx\_ref\_tag">5</span></a>}\text{)}}{\underbrace{\left[{\displaystyle \sum _{{\tau }^{\prime \prime }\in {𝒯}_{s,f}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime \prime }}}{\widehat{P}}_F(s_{t+1}\mid s_t)\right]}}$

  $={\displaystyle \sum _{{\tau }^{\prime }\in {𝒯}_{0,s}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime }}}{\widehat{P}}_F(s_{t+1}\mid s_t).$

  Similarly, any complete trajectory going through $s\to s^{\prime }$ can be (uniquely) decomposed into a partial trajectory ${\tau }^{\prime }\in {𝒯}_{0,s}$ from $s_0$ to $s$, and a partial trajectory ${\tau }^{\prime \prime }\in {𝒯}_{s^{\prime },f}$ from $s^{\prime }$ to $s_f$. Again, using the definition of $P(s\to s^{\prime })$:

  $P(s\to s^{\prime })$ $={\displaystyle \sum _{\tau :(s\to s^{\prime })\in \tau }}P(\tau )={\displaystyle \sum _{\tau :(s\to s^{\prime })\in \tau }}{\displaystyle \prod _{(s_t\to s_{t+1})\in \tau }}{\widehat{P}}_F(s_{t+1}\mid s_t)$

  $=\underset{=P(s)}{\underbrace{\left[{\displaystyle \sum _{{\tau }^{\prime }\in {𝒯}_{0,s}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime }}}{\widehat{P}}_F(s_{t+1}\mid s_t)\right]}}{\widehat{P}}_F(s^{\prime }\mid s)\underset{=\mathrm{\hspace{0.17em}1}\text{(}\text{<a class="ltx\_ref" href="\#Thmtheorem5" style="font-size:70\%;" title="Lemma 5. ‣ 2.1 Some elements of graph theory ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Lemma</span> <span class="ltx\_text ltx\_ref\_tag">5</span></a>}\text{)}}{\underbrace{\left[{\displaystyle \sum _{{\tau }^{\prime \prime }\in {𝒯}_{s^{\prime },f}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime \prime }}}{\widehat{P}}_F(s_{t+1}\mid s_t)\right]}}$

  $=P(s){\widehat{P}}_F(s^{\prime }\mid s).$

  Combining the two results above, we get:

  $$P_F(s^{\prime }\mid s)=\frac{P(s\to s^{\prime })}{P(s)}={\widehat{P}}_F(s^{\prime }\mid s).$$

• $2\Rightarrow 1$: Suppose that there exists a probability function ${\widehat{P}}_F$ consistent with $G$ such that for some complete trajectory $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)$

  $$P(\tau )=\prod _{t=1}^{n+1}{\widehat{P}}_F(s_t\mid s_{t-1}).$$

  For the same reasons as those used to justify the uniqueness in the $1\Rightarrow 2$ proof, ${\widehat{P}}_F$ is necessarily equal to the forward transition probability $P_F$, associated with $P$.

  We now want to show that the flow $F$ associated with $P$ is Markovian, by showing the Markov property from Eq. 19. Let ${\tau }^{\prime }\in {𝒯}_{0,s}$ be any partial trajectory from $s_0$ to $s$; using the definition of conditional probability:

  $$P(s\to s^{\prime }\mid {\tau }^{\prime })=\frac{P(s_0\to \mathrm{\dots }\to s\to s^{\prime })}{P(s_0\to \mathrm{\dots }\to s)}.$$

  Following the same idea as above, we will now rewrite $P(s_0\to \mathrm{\dots }\to s)$, as a sum over complete trajectories that share the same prefix trajectory ${\tau }^{\prime }$. Any such complete trajectory $\tau$ can be (uniquely) decomposed into this common prefix ${\tau }^{\prime }$, and a partial trajectory ${\tau }^{\prime \prime }\in {𝒯}_{s,f}$ from $s$ to $s_f$.

  $P(s_0\to \mathrm{\dots }\to s)$ $={\displaystyle \sum _{\tau :{\tau }^{\prime }\subseteq \tau }}P(\tau )={\displaystyle \sum _{\tau :{\tau }^{\prime }\subseteq \tau }}{\displaystyle \prod _{(s_t\to s_{t+1})\in \tau }}{P}_F(s_{t+1}\mid s_t)$

  $=\left[{\displaystyle \prod _{s_{t-1}\to s_t\in {\tau }^{\prime }}}{P}_F(s_t\mid s_{t-1})\right]\underset{=\mathrm{\hspace{0.17em}1}\text{(}\text{<a class="ltx\_ref" href="\#Thmtheorem5" style="font-size:70\%;" title="Lemma 5. ‣ 2.1 Some elements of graph theory ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Lemma</span> <span class="ltx\_text ltx\_ref\_tag">5</span></a>}\text{)}}{\underbrace{\left[{\displaystyle \sum _{{\tau }^{\prime \prime }\in {𝒯}_{s,f}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime \prime }}}{P}_F(s_{t+1}\mid s_t)\right]}}$

  $={\displaystyle \prod _{s_{t-1}\to s_t\in {\tau }^{\prime }}}{P}_F(s_t\mid s_{t-1}).$

  Similarly, any complete trajectory $\tau$ that share the same prefix trajectory $(s_0,\mathrm{\dots },s,s^{\prime })$ can be (uniquely) decomposed into this common prefix, and a partial trajectory ${\tau }^{\prime \prime }\in {𝒯}_{s^{\prime },f}$ from $s^{\prime }$ to $s_f$, leading to:

  $P(s_0\to \mathrm{\dots }\to s\to s^{\prime })=P(s_0\to \mathrm{\dots }\to s)P_F(s^{\prime }\mid s)$

  Combining the two results above, we can conclude that $P$ satisfies the Markov property, and therefore that the flow $F$ is Markovian:

  $$P(s^{\prime }\to s\mid {\tau }^{\prime })=\frac{P(s_0\to \mathrm{\dots }\to s\to s^{\prime })}{P(s_0\to \mathrm{\dots }\to s)}=P_F(s^{\prime }\mid s)=P(s^{\prime }\to s\mid s)$$

• $\{1,2\}\Rightarrow 3$: Suppose that $F$ is a Markovian flow. We have shown above that this is equivalent to $P$ being decomposed into a product of forward transition probabilities $P_F$. For some complete trajectory $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)$:

  $$P(\tau )=\prod _{t=1}^{n+1}{P}_F(s_t\mid s_{t-1})=\prod _{t=1}^{n+1}\frac{P(s_{t-1}\to s_t)}{P(s_{t-1})}=\prod _{t=1}^{n+1}\frac{P(s_{t-1}\to s_t)}{P(s_t)}=\prod _{t=1}^{n+1}{P}_B(s_{t-1}\mid s_t),$$

  where the third equality uses the fact that $P(s_0)=P(s_f)=1$, and using the definition of the backwards transition probability $P_B$. The proof of uniqueness of $P_B$ is similar to that of $P_F$ in $1\Rightarrow 2$, and uses:

  $P(s\to s^{\prime })$ $={\displaystyle \sum _{\tau :(s\to s^{\prime })\in \tau }}P(\tau )={\displaystyle \sum _{\tau :(s\to s^{\prime })\in \tau }}{\displaystyle \prod _{(s_t\to s_{t+1})\in \tau }}{\widehat{P}}_B(s_t\mid s_{t+1})$

  $=\underset{=\mathrm{\hspace{0.17em}1}\text{(}\text{<a class="ltx\_ref" href="\#Thmtheorem5" style="font-size:70\%;" title="Lemma 5. ‣ 2.1 Some elements of graph theory ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Lemma</span> <span class="ltx\_text ltx\_ref\_tag">5</span></a>}\text{)}}{\underbrace{\left[{\displaystyle \sum _{{\tau }^{\prime }\in {𝒯}_{0,s}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime }}}{\widehat{P}}_B(s_t\mid s_{t+1})\right]}}{\widehat{P}}_B(s\mid s^{\prime })\underset{=P(s^{\prime })}{\underbrace{\left[{\displaystyle \sum _{{\tau }^{\prime \prime }\in {𝒯}_{s^{\prime },f}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in {\tau }^{\prime \prime }}}{\widehat{P}}_B(s_t\mid s_{t+1})\right]}}$

  $=P(s^{\prime }){\widehat{P}}_B(s\mid s^{\prime }),$

• $3\Rightarrow 1$: Similar to the proof of $2\Rightarrow 1$, ${\widehat{P}}_B$ is necessarily equal to the backwards transition probability $P_B$ associated with $P$. Additionally, $P_B$ is related to the forward transition probability $P_F$:

  $$P(s\to s^{\prime })=P_B(s\mid s^{\prime })P(s^{\prime })=P_F(s^{\prime }\mid s)P(s).$$

  We can therefore write the decomposition of $P$ in terms of $P_F$, instead of $P_B$. For some complete trajectory $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)$:

  $P(\tau )$ $={\displaystyle \prod _{t=1}^{n+1}}{P}_B(s_{t-1}\mid s_t)={\displaystyle \prod _{t=1}^{n+1}}{\displaystyle \frac{P(s_{t-1})}{P(s_t)}}{P}_F(s_{t+1}\mid s_t)={\displaystyle \frac{P(s_0)}{P(s_f)}}{\displaystyle \prod _{t=1}^{n+1}}{P}_F(s_{t+1}\mid s_t)$

  $={\displaystyle \prod _{t=1}^{n+1}}{P}_F(s_{t+1}\mid s_t),$

  where we used the fact that $P(s_0)=P(s_f)=1$. Using “$2\Rightarrow 1$”, we can conclude that $F$ is a Markovian flow.

The decomposition of Eq. 20 shows how Markovian flows can be used to draw terminating states from the terminating state probability $P_T$ (Eq. 18). Namely, we have the following result:

Corollary 17.

Let $(G,F)$ be a Markovian flow network, and $P_F$ the corresponding forward transition probability. Consider the procedure starting from $s=s_0$, and iteratively drawing one sample from $P_F(.\mid s)$ until reaching $s_f$. Then the probability of the procedure terminating in a state $s$ is $P_T(s)$.

Proof First, note that the procedure terminates with probability 1, given that $G$ is acyclic.

For the procedure to terminate in a state $s$, it means that the trajectory $\tau \in 𝒯$ implicitly constructed during the procedure contains the edge $s\to s_f$. The probability of the procedure terminating in $s$ is thus:

${\displaystyle \sum _{\tau \in 𝒯:s\to s_f\in \tau }}\underset{P(\tau )\text{, according to}\text{<a class="ltx\_ref" href="\#S2.E20" style="font-size:70\%;" title="In Item 2 ‣ Proposition 16. ‣ 2.4 Markovian Flows ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Eq.</span> <span class="ltx\_text ltx\_ref\_tag">20</span></a>}}{\underbrace{{\displaystyle \prod _{s^{\prime }\to s^{\prime \prime }\in \tau }}{P}_F(s^{\prime \prime }\mid s^{\prime })}}=P(s\to s_f)=P_T(s)$

The following proposition shows that, as a consequence of the Prop. 16, we obtain three different parametrizations of Markovian flows.

Proposition 18.

Given a pointed DAG $G=(𝒮,𝔸)$, a Markovian flow on $G$ is completely and uniquely specified by one of the following:

1. the combination of the total flow $\widehat{Z}$ and the forward transition probabilities ${\widehat{P}}_F(s^{\prime }\mid s)$ for all edges $s\to s^{\prime }\in 𝔸$,

2. the combination of the total flow $\widehat{Z}$ and the backward transition probabilities ${\widehat{P}}_B(s\mid s^{\prime })$ for all edges $s\to s^{\prime }\in 𝔸$.

3. the combination of the terminating flows $\widehat{F}(s\to s_f)$ for all terminating edges $s\to s_f\in {𝔸}^f$ and the backwards transition probabilities ${\widehat{P}}_B(s\mid s^{\prime })$ for all non-terminating edges $s\to s^{\prime }\in {𝔸}^{-f}$,

Proof In the first two settings, we define a flow function $F:𝒯\to {ℝ}^{+}$, at a trajectory $\tau =(s_0,s_1,\mathrm{\dots },s_n,s_{n+1}=s_f)$ as:

1. $F(\tau ):=\widehat{Z}{\prod }_{t=1}^{n+1}{\widehat{P}}_F(s_t\mid s_{t-1})$,

2. $F(\tau ):=\widehat{Z}{\prod }_{t=1}^{n+1}{\widehat{P}}_B(s_{t-1}\mid s_t)$

We need to prove that it is the only Markovian flow that can be defined for both settings. The proof for the third setting will follow from that of the second setting.

First setting:

First, we need to show that the total flow $Z$ associated with the flow function $F$ (Eq. 10) matches $\widehat{Z}$. This is a consequence of Lemma 5:

$$Z=\sum _{\tau \in 𝒯}F(\tau )=\widehat{Z}\underset{=1\text{, according to}\text{<a class="ltx\_ref" href="\#Thmtheorem5" style="font-size:70\%;" title="Lemma 5. ‣ 2.1 Some elements of graph theory ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Lemma</span> <span class="ltx\_text ltx\_ref\_tag">5</span></a>}}{\underbrace{\sum _{\tau =(s_0,s_1,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯}\prod _{t=1}^{n+1}{\widehat{P}}_F(s_t\mid s_{t-1})}}=\widehat{Z}$$

Then, we need to show that the forward transition probability function $P_F$ associated with $F$ (Eq. 16) matches ${\widehat{P}}_F$, and that the flow $F$ is Markovian. To this end, note that the corresponding flow probability $P$ satisfies Eq. 20. Thus, as a consequence of Prop. 16, $F$ is a Markovian flow, and its forward transition probability function is ${\widehat{P}}_F$.

As a last requirement, we need to show that if a Markovian flow $F^{\prime }$ has a partition function $Z^{\prime }=\widehat{Z}$ and a forward transition probability function $P_F^{\prime }={\widehat{P}}_F$, then it is necessarily equal to $F$. This is a direct consequence of Prop. 16, given that for any $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯$:

$$F^{\prime }(\tau )=Z^{\prime }\prod _{t=1}^{n+1}{P}_F^{\prime }(s_t\mid s_{t-1})=\widehat{Z}\prod _{t=1}^{n+1}{\widehat{P}}_F(s_t\mid s_{t-1})=F(\tau )$$

Second setting:

First, we show that as a consequence of Lemma 5, the total flow $Z$ associated with $F$ matches $\widehat{Z}$:

$$Z=\sum _{\tau \in 𝒯}F(\tau )=\widehat{Z}\underset{=1\text{, according to}\text{<a class="ltx\_ref" href="\#Thmtheorem5" style="font-size:70\%;" title="Lemma 5. ‣ 2.1 Some elements of graph theory ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Lemma</span> <span class="ltx\_text ltx\_ref\_tag">5</span></a>}}{\underbrace{\sum _{\tau =(s_0,s_1,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯}\prod _{t=1}^{n+1}{\widehat{P}}_B(s_{-1}\mid s_t)}}=\widehat{Z}$$

Second, we note that the flow probability $P$ associated with $F$ satisfies Eq. 21. Thus, as a consequence of Prop. 16, $F$ is a Markovian flow, and its backward transition probability function is ${\widehat{P}}_B$.

Finally, if a Markovian flow $F^{\prime }$ has a partition function $Z^{\prime }=\widehat{Z}$ and a backward transition probability function $P_B^{\prime }={\widehat{P}}_B$, then following Prop. 16, $\forall \tau \in 𝒯,F^{\prime }(\tau )=F(\tau )$.

Third setting:

From the terminating flows $\widehat{F}(s\to s_f)$ and the backwards transition probabilities ${\widehat{P}}_B(s\mid s^{\prime })$ for non-terminating edges, we can uniquely define a total flow $\widehat{Z}$, and extend ${\widehat{P}}_B$ to all edges as follows:

$\widehat{Z}$ $:={\displaystyle \sum _{s\in Par(s_f)}}\widehat{F}(s\to s_f)$

${\widehat{P}}_B(s\mid s^{\prime })$ $:=\{\begin{array}{cc}{\widehat{P}}_B(s\mid s^{\prime })\text{if}{s}^{\prime }\ne s_f\hfill & \\ \frac{\widehat{F}(s\to s_f)}{\widehat{Z}}\text{otherwise.}\hfill & \end{array}$

This takes us back to the second setting, for which we have already proven that with $\widehat{Z}$ and ${\widehat{P}}_B$ defined for all edges, a Markovian flow is uniquely defined.

2.5 Flow Matching Conditions

In Prop. 18, we saw how forward and backward probability functions can be used to uniquely define a Markovian flow. We will show in the next proposition how non-negative functions of states and edges can be used to define a Markovian flow. Such functions cannot be unconstrained (as ${\widehat{P}}_F$ and $\widehat{Z}$ in Prop. 18 e.g.), as we have seen in Prop. 8.

Proposition 19.

Let $G=(𝒮,𝔸)$ be a pointed DAG. Consider a non-negative function $\widehat{F}$ taking as input either a state $s\in 𝒮$ or a transition $s\to s^{\prime }\in 𝔸$. Then $\widehat{F}$ corresponds to a flow if and only if the flow matching conditions:

$\forall s^{\prime }>s_0,\widehat{F}(s^{\prime })$ $={\displaystyle \sum _{s\in Par(s^{\prime })}}\widehat{F}(s\to s^{\prime })$

$={\displaystyle \sum _{s^{\prime \prime }\in Child(s^{\prime })}}\widehat{F}(s^{\prime }\to s^{\prime \prime })$ (22)

are satisfied. More specifically, $\widehat{F}$ uniquely defines a Markovian flow $F$ matching $\widehat{F}$ on states and transitions:

$$\forall \tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯F(\tau )=\frac{{\prod }_{t=1}^{n+1}\widehat{F}(s_{t-1}\to s_t)}{{\prod }_{t=1}^n\widehat{F}(s_t)}.$$ (23)

Proof Necessity is a direct consequence of Prop. 8. Let’s show sufficiency. Let ${\widehat{P}}_F$ be the forward probability function defined by:

$\forall s\to s^{\prime }\in 𝔸{\widehat{P}}_F(s^{\prime }\mid s):={\displaystyle \frac{\widehat{F}(s\to s^{\prime })}{\widehat{F}(s)}}.$

${\widehat{P}}_F$ is consistent with $G$ given that $\widehat{F}$ satisfies the flow matching conditions (Eq. 22). Let $\widehat{Z}=\widehat{F}(s_0)$. According to Prop. 18, there exists a unique Markovian flow $F$ with forward transition probability function $P_F={\widehat{P}}_F$ and partition function $Z=\widehat{Z}$, and such that for a trajectory $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯$:

$$\forall \tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯F(\tau )=\widehat{Z}\prod _{t=1}^{n+1}{\widehat{P}}_F(s_t\mid s_{t-1})=\frac{{\prod }_{t=1}^{n+1}\widehat{F}(s_{t-1}\to s_t)}{{\prod }_{t=1}^n\widehat{F}(s_t)}.$$ (24)

Additionally, similar to the proof of Prop. 16, we can write for any state $s^{\prime }\ne s_0$:

$F(s^{\prime })$ $=\widehat{Z}{\displaystyle \sum _{\tau \in {𝒯}_{0,s^{\prime }}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in \tau }}{\widehat{P}}_F(s_{t+1}\mid s_t)$

$=\widehat{Z}{\displaystyle \frac{\widehat{F}(s^{\prime })}{\widehat{F}(s_0)}}\underset{=1\text{, according to}\text{<a class="ltx\_ref" href="\#Thmtheorem5" style="font-size:70\%;" title="Lemma 5. ‣ 2.1 Some elements of graph theory ‣ 2 Flow Networks and Markovian Flows ‣ GFlowNet Foundations"><span class="ltx\_text ltx\_ref\_tag">Lemma</span> <span class="ltx\_text ltx\_ref\_tag">5</span></a>}}{\underbrace{{\displaystyle \sum _{\tau \in {𝒯}_{0,s^{\prime }}}}{\displaystyle \prod _{(s_t\to s_{t+1})\in \tau }}{\widehat{P}}_B(s_t\mid s_{t+1})}}$

$=\widehat{F}(s^{\prime }),$

where ${\widehat{P}}_B(s^{\prime }\mid s):=\frac{\widehat{F}(s\to s^{\prime })}{\widehat{F}(s^{\prime })}$ defines a backward probability function consistent with $G$. And because $\forall s\to s^{\prime }\in 𝔸P_F(s^{\prime }\mid s)={\widehat{P}}_F(s^{\prime }\mid s)$, it follows that $\forall s\to s^{\prime }\in 𝔸F(s\to s^{\prime })=\widehat{F}(s\to s^{\prime })$.

To show uniqueness, let’s consider a Markovian flow $F^{\prime }$ that matches $\widehat{F}$ on states and edges. Following Prop. 16, for any trajectory $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯$

$F^{\prime }(\tau )=\widehat{Z}{\displaystyle \prod _{t=1}^{n+1}}{\widehat{P}}_F(s_t\mid s_{t-1})={\displaystyle \frac{{\prod }_{t=1}^{n+1}\widehat{F}(s_{t-1}\to s_t)}{{\prod }_{t=1}^n\widehat{F}(s_t)}}=F(\tau ).$

Note how Eq. 22 can be used to recursively define the flow in all the states if $Z$ is given and either the forward or the backwards transition probabilities are given. Either way, we would start from the flow at one of the extreme states $s_0$ or $s_f$ and then distribute it recursively through the directed acyclic graph of the flow network, either going forward or going backward. A setting of particular interest, that will be central in Sec. 3, is when we are given all the terminal flows $F(s\to s_f)$, and we would like to deduce a state flow function $F(s)$ and a forward transition probability function $P_F(s^{\prime }\mid s)$ for the rest of the flow network.

Next, we will see how to parametrize Markovian flows using forward and backward probability functions consistent with the DAG. Unlike the condition in Prop. 19, the new condition does not involve a sum over transitions, which could be problematic if each state can have a large number of successors or if the state-space is continuous. Interestingly, the resulting condition is analogous to the detailed balance condition of Monte-Carlo Markov chains.

Definition 20.

Given a pointed DAG $G=(𝒮,𝔸)$, a forward transition probability function ${\widehat{P}}_F$ and a backward transition probability function ${\widehat{P}}_B$ consistent with $G$, ${\widehat{P}}_F$ and ${\widehat{P}}_B$ are compatible if there exists an edge flow function $\widehat{F}:𝔸\to {ℝ}^{+}$ such that

$\forall s\to s^{\prime }\in 𝔸{\widehat{P}}_F(s^{\prime }\mid s)$ $={\displaystyle \frac{\widehat{F}(s\to s^{\prime })}{{\sum }_{s^{\prime }\in Child(s)}\widehat{F}(s\to s^{\prime })}},{\widehat{P}}_B(s\mid s^{\prime })$ $={\displaystyle \frac{\widehat{F}(s\to s^{\prime })}{{\sum }_{s^{\prime \prime }\in Par(s^{\prime })}\widehat{F}(s^{\prime \prime }\to s^{\prime })}}$ (25)

Proposition 21.

Let $G=(𝒮,𝔸)$ be a pointed DAG. Consider a non-negative function $\widehat{F}$ over states, a forward transition probability function ${\widehat{P}}_F$ and a backwards transition probability function ${\widehat{P}}_B$ consistent with $G$. Then, $\widehat{F},{\widehat{P}}_B,$ and $\widehat{P_F}$ jointly correspond to a flow if and only if the detailed balance conditions holds:

$$\forall s\to s^{\prime }\in 𝔸\widehat{F}(s){\widehat{P}}_F(s^{\prime }\mid s)=\widehat{F}(s^{\prime }){\widehat{P}}_B(s\mid s^{\prime }).$$ (26)

More specifically, $\widehat{F},{\widehat{P}}_F,$ and ${\widehat{P}}_B$ uniquely define a Markovian flow $F$ matching $\widehat{F}$ on states, and with transition probabilities matching ${\widehat{P}}_F$ and ${\widehat{P}}_B$. Furthermore, when this condition is satisfied, the forward and backward transition probability functions ${\widehat{P}}_F$ and ${\widehat{P}}_B$ are compatible.

Proof For necessity, consider a flow $F$, with state flow function denoted $F$, and forward and backward transitions $P_F$ and $P_B$. It is clear from the definition of $P_F$ and $P_B$ (Def. 12) that Eq. 26 holds. We prove the sufficiency of the condition by first defining the edge flow

$$\forall s\to s^{\prime }\in 𝔸\widehat{F}(s\to s^{\prime }):=\widehat{F}(s){\widehat{P}}_F(s^{\prime }\mid s).$$ (27)

We then sum both sides of Eq. 26 over $s$, yielding

$$\forall s^{\prime }>s0\sum _{s\in Par(s^{\prime })}\widehat{F}(s){\widehat{P}}_F(s^{\prime }\mid s)=\widehat{F}(s^{\prime })\sum _{s\in Par(s^{\prime })}{\widehat{P}}_B(s\mid s^{\prime })=\widehat{F}(s^{\prime })$$ (28)

where we used the fact that ${\widehat{P}}_B$ is a normalized probability distribution. Combining this with Eq. 27, we get

$$\forall s^{\prime }>s_0\widehat{F}(s^{\prime })=\sum _{s\in Par(s^{\prime })}\widehat{F}(s\to s^{\prime })$$ (29)

which is the first equality of the flow-matching condition (Eq. 22) of Prop. 19. We can obtain the second equality by first using the normalization of $\widehat{P}$, and then using our definition of the edge flow (Eq. 27):

$\forall s^{\prime }>s_0\widehat{F}(s^{\prime })$ $=\widehat{F}(s^{\prime }){\displaystyle \sum _{s^{\prime \prime }\in Child(s^{\prime })}}{\widehat{P}}_F(s^{\prime \prime }\mid s^{\prime })$

$={\displaystyle \sum _{s^{\prime \prime }\in Child(s^{\prime })}}\widehat{F}(s^{\prime }){\widehat{P}}_F(s^{\prime \prime }\mid s^{\prime })$

$={\displaystyle \sum _{s^{\prime \prime }\in Child(s^{\prime })}}\widehat{F}(s^{\prime }\to s^{\prime \prime }).$ (30)

Following Prop. 19, there exists a unique Markovian flow $F$ with state and edge flows given by $\widehat{F}$. Using Eq. 27 and Eq. 26, it follows that $F$ has transition probabilities ${\widehat{P}}_F$ and ${\widehat{P}}_B$ as required. The uniqueness is also a consequence of Eq. 27. This proves sufficiency.

To show that ${\widehat{P}}_F$ and ${\widehat{P}}_B$ are compatible (Def. 20), we first combine Eq. 27 and Eq. 30 (with relabeling of variables) to obtain

$$\forall s\to s^{\prime }\in 𝔸{\widehat{P}}_F(s^{\prime }\mid s)=\frac{\widehat{F}(s\to s^{\prime })}{{\sum }_{s^{\prime }\in Child(s)}\widehat{F}(s\to s^{\prime })},$$

we then isolate ${\widehat{P}}_B$ in Eq. 26, yielding

$$\forall s\to s^{\prime }\in 𝔸{\widehat{P}}_B(s\mid s^{\prime })=\frac{\widehat{F}(s)}{\widehat{F}(s^{\prime })}{\widehat{P}}_F(s^{\prime }\mid s)=\frac{\widehat{F}(s\to s^{\prime })}{\widehat{F}(s^{\prime })}=\frac{\widehat{F}(s\to s^{\prime })}{{\sum }_{s^{\prime \prime }\in Par(s^{\prime })}\widehat{F}(s^{\prime \prime }\to s^{\prime })},$$

We thus get Eq. 25 of Def. 20, as desired.

At first glance, it may seem that when ${\widehat{P}}_B$ is unconstrained, the detailed balance condition can trivially be achieved by setting

$$\forall s\to s^{\prime }\in 𝔸{\widehat{P}}_B(s\mid s^{\prime })=\frac{{\widehat{P}}_F(s^{\prime }\mid s)\widehat{F}(s)}{\widehat{F}(s^{\prime })}$$ (31)

However, because we also have the constraint ${\sum }_{s\in Par(s^{\prime })}{\widehat{P}}_B(s\mid s^{\prime })=1$, then Eq. 31 can only be satisfied if the flows are consistent with the forward transition:

$$\sum _{s\in Par(s^{\prime })}{\widehat{P}}_F(s^{\prime }\mid s)\widehat{F}(s)=\widehat{F}(s^{\prime }).$$

2.6 Backwards Transitions can be Chosen Freely

Consider the setting in which we are given terminating flows to be matched, i.e. where the goal is to find a flow function with the right terminating flows. This is the setting introduced in Bengio et al. (2021), and that will be studied in Sec. 3. In this case, Prop. 18 tells us that in order to fully determine the forward transition probabilities and the state or state-action flows, it is not sufficient in general to specify only the terminating flows; it is also necessary to specify the backwards transition probabilities on the edges other than the terminal ones (the latter being given by the terminating flows).

What this means is that the terminating flows do not specify the flow completely, e.g., because many different paths can land in the same terminating state. The preference over such different ways to achieve the same final outcome is specified by the backwards transition probability $P_B$ (except for $P_B(s\mid s_f)$ which is a function of the terminating flows and $Z$). For example, we may want to give equal weight to all parents of a node $s$, or we may prefer shorter paths, which can be achieved if we keep track in the state $s$ of the length of the shortest path to the node $s$, or we may let a learner discover a $P_B$ that makes learning $P_F$ or $F$ easier.

2.7 Equivalence Between Flows

In the previous sections, we have seen that Markovian flows have the property that trajectory flows or probabilities factorize according to the DAG, and we have seen different ways of characterizing Markovian flows. In Sec. 3, we show how to approximate Markovian flows in order to define probability measures over terminating states. In this section, through an equivalence relation between trajectory flows, we justify the focus on Markovian flows. Given a pointed DAG $G=(𝒮,𝔸)$, we denote by:

• $ℱ(G)$: the set of flows on $G$, i.e. the set of functions from $𝒯$, the set of complete trajectories in $G$, to ${ℝ}^{+}$,

• ${ℱ}_{Markov}(G)$: the set of flows in $ℱ(G)$ that are Markovian.

Definition 22.

Let $G=(𝒮,𝔸)$ be a pointed DAG, and $F_1,F_2\in ℱ(G)$ two trajectory flow functions. We say that $F_1$ and $F_2$ are equivalent if they coincide on edge-flows, i.e.:

$$\forall s\to s^{\prime }\in 𝔸F_1(s\to s^{\prime })=F_2(s\to s^{\prime })$$

Fig. 4 shows four flow functions in a simple pointed DAG that are pairwise equivalent.

$\tau$	$F_1(\tau )$	$F_2(\tau )$	$F_3(\tau )$	$F_4(\tau )$
$s_0,s_2,s_f$	$1$	$4/5$	$1$	$6/5$
$s_0,s_1,s_2,s_f$	$1$	$6/5$	$1$	$4/5$
$s_0,s_2,s_3,s_f$	$1$	$6/5$	2	$9/5$
$s_0,s_1,s_2,s_3,s_f$	$2$	$9/5$	$1$	$6/5$

Figure 4: Equivalent flows and Markovian flows. Flows $F_1$ and $F_2$ are equivalent. $F_3$ and $F_4$ are equivalent, but not equivalent to $F_1$ and $F_2$. $F_2$ and $F_4$ are Markovian. $F_1$ and $F_3$ are not Markovian. $F_1,F_2,F_3$ and $F_4$ coincide on the terminating flows i.e. at $s_2\to s_f$ and $s_3\to s_f$.

This defines an equivalence relation (i.e., a relation that is reflexive, symmetric, and transitive). Hence, each flow $F$ belongs to an equivalence class, and the set of flows $ℱ(G)$ can be partitioned into equivalence classes. Note that if two flows are equivalent, then the corresponding state flow functions also coincide (as a direct consequence of Prop. 8).

Proposition 23.

Given a pointed DAG $G$. If two flow function $F_1,F_2\in {ℱ}_{Markov}(G)$ are equivalent, then they are equal. Additionally, for any flow function $F^{\prime }\in ℱ(G)$, there exists a unique Markovian flow function $F\in {ℱ}_{Markov}(G)$ such that $F$ and $F^{\prime }$ are equivalent.

Proof Because $F_1$ and $F_2$ are Markovian, then for any trajectory $\tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)$:

$F_1(\tau )$ $={\displaystyle \frac{{\prod }_{t=1}^{n+1}{F}_1(s_{t-1}\to s_t)}{{\prod }_{t=1}^n{F}_1(s_t)}}$

$={\displaystyle \frac{{\prod }_{t=1}^{n+1}{F}_2(s_{t-1}\to s_t)}{{\prod }_{t=1}^n{F}_2(s_t)}}$

$=F_2(\tau ),$

where we combined the definition of equivalent flows and Prop. 16.

Given a flow function $F^{\prime }$, because its state and edge flow functions satisfy the flow matching conditions (as a consequence of Prop. 8), then according to Prop. 19, the flow $F$ defined by:

$$\forall \tau :=(s_0,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯F(\tau )=\frac{{\prod }_{t=1}^{n+1}{F}^{\prime }(s_{t-1}\to s_t)}{{\prod }_{t=1}^n{F}^{\prime }(s_t)}$$

is Markovian, and coincides with $F^{\prime }$ on state and edge flows. Combining this with the statement above, we conclude that $F$ is the unique Markovian flow that is equivalent to $F^{\prime }$.
The previous proposition shows that in each equivalence class stands out a particular flow function, that has a property the other flows in the same equivalence class don’t have: it is Markovian.

A consequence of this is that, if we care essentially about state and edge flows, instead of dealing with the full set of flows $ℱ(G)$, it suffices to restrict any flow learning problem to the set of Markovian flows ${ℱ}_{Markov}(G)$. The advantage of this restriction is that defining a flow requires the specification of $F(\tau )$ for all trajectories $\tau \in 𝒯$, whereas defining a Markovian flow requires the specification of $F(s\to s^{\prime })$ for all edges $s\to s^{\prime }\in 𝔸$, which is generally exponentially smaller than $𝒯$ (note that the edge flows still need to satisfy the flow-matching conditions in Prop. 19). Thus, in order to approximate or learn a flow function that satisfies some conditions on its edge or state values, it suffices to approximate or learn a Markovian flow, by learning the edge flow function, which is a much smaller object than the actual flow function.

3 GFlowNets: Learning a Flow

With the theoretical preliminaries established in Sec. 1 and Sec. 2, we now consider the general class of problems introduced by Bengio et al. (2021) where some constraints or preferences over flows are given. Our goal is to find functions such as the state flow function $F(s)$ or the transition probability function $P(s\to s^{\prime }\mid s)$ that best match these desiderata using corresponding estimators $\widehat{F}(s)$ and $\widehat{P}(s\to s^{\prime }\mid s)$ which may not correspond to a proper flow. Such learning machines are called Generative Flow Networks (or GFlowNets for short). We focus on scenarios where we are given a target reward function $R:{𝒮}^f\to {ℝ}^{+}$, and aim at estimating flows $F$ that satisfy:

$$\forall s\in {𝒮}^{f}F(s\to s_f)=R(s)$$ (32)

Because of the equivalences that exist in the set of flows, then without loss of generality, we choose GFlowNets to approximate Markovian flows only. We are thus interested in the following set of flows:

$${ℱ}_{Markov}(G,R)=\{F\in {ℱ}_{Markov}(G),\forall s\in {𝒮}^{f}F(s\to s_f)=R(s)\}$$ (33)

For now, we informally define a GFlowNet as an estimator of a Markovian flow function $F\in {ℱ}_{Markov}(G,R)$. We provide a more formal definition later-on.

With an estimator $\widehat{F}$ of such a Markovian flow $F$, we can define an approximate forward transition probability function ${\widehat{P}}_F$, as in Prop. 16, in order to draw trajectories $\tau \in 𝒯$ (the set of complete trajectories in $G$) by iteratively sampling each state given the previous one, starting at $s_0$ and then with $s_{t+1}\sim {\widehat{P}}_F(.\mid s_t)$ until we reach the sink state $s_{n+1}=s_f$ for some $n$.

Next, we will clarify how such an estimator can be obtained.

3.1 GFlowNets as an Alternative to MCMC Sampling

The main established methods to approximately sample from the distribution associated with an energy function $ℰ$ are Monte-Carlo Markov chain (MCMC) methods, which require significant computation (running a potentially very long Markov chain) to obtain samples. Instead, the GFlowNet approach amortizes upfront computation to train a generator that yields very efficient computation (a single configuration is constructed, no chain needed) for each new sample. For example, Bengio et al. (2021) build a GFlowNet that constructs a molecule via a small sequence of actions, each of which adds an atom or a molecular substructure to an existing molecule represented by a graph, starting from an empty graph. Only one such configuration needs to be considered, in contrast with MCMC methods, which require potentially very long chains of such configurations, and suffer from the challenge of mode-mixing (Jasra et al., 2005; Bengio et al., 2013; Pompe et al., 2020), which can take time exponentially long in the distance between modes. In GFlowNets, this computational challenge is avoided but the computational demand is converted to that of training the GFlowNet. To see how this can be extremely beneficial, consider having already constructed some configurations $x$ and obtained their unnormalized probability or reward $R(x)$. With these pairs $(x,R(x))$, a machine learning system could potentially generalize about the value of $R$ elsewhere, and if it is a generative model, sample new $x$’s in places of large $R(x)$. Hence, if there is an underlying statistical structure in how the modes of $R$ are related to each other, a generative learner that generalizes could guess the presence of modes it has not visited yet, taking advantage of the patterns it has already uncovered from the $(x,R(x))$ pairs it has seen. On the other hand, if there is no structure (the modes are randomly placed), then we should not expect GFlowNets to do significantly better than MCMC because training becomes intractable in high-dimensional spaces (since it requires visiting every area of the configuration space to ascertain its reward).

3.2 GFlowNets and flow-matching losses

We have seen in Sec. 2.4 and Sec. 2.5 different ways of parametrizing a flow. For example, with a partition function and forward transition probabilities, or with edge flows that satisfy the flow matching conditions. Because there are many ways to parametrize GFlowNets, we start with an abstract formulation for them, where $o\in 𝒪$ represents a parameter configuration (e.g., resulting from or while training of a GFlowNet), $\mathrm{\Pi }(o)$ gives the corresponding probability measure over trajectories $\tau \in 𝒯$, and $\mathscr{H}$ maps a Markovian flow $F$ to its parametrization $o$. In the following definition, we show what conditions should be satisfied in order for such a parametrization to be valid.

Definition 24.

Given a pointed DAG $G=(𝒮,𝔸)$, with an initial and sink states $s_0$ and $s_f$ respectively, and a target reward function $R:{𝒮}^f\to {ℝ}^{+}$, we say that the triplet $(𝒪,\mathrm{\Pi },\mathscr{H})$ is a flow parametrization of $(G,R)$ if:

1. $𝒪$ is a non-empty set,

2. $\mathrm{\Pi }$ is a function mapping each object $o\in 𝒪$ to an element $\mathrm{\Pi }(o)\in \mathrm{\Delta }(𝒯)$, the set of probability distributions on $𝒯$,

3. $\mathscr{H}$ is an injective functional from ${ℱ}_{Markov}(G,R)$ to $𝒪$,

4. For any $F\in {ℱ}_{Markov}(G,R)$, $\mathrm{\Pi }(\mathscr{H}(F))$ is the probability measure associated with the flow $F$ (Def. 11).

To each object $o\in 𝒪$, the distribution $\mathrm{\Pi }(o)$ implicitly defines a terminating state probability measure:

$$\forall s\in {𝒮}^f{P}_T(s):=\sum _{\tau \in 𝒯:s\to s_f\in \tau }\mathrm{\Pi }(o)(\tau ),$$ (34)

where the dependence on $o$ in $P_T$ is omitted for clarity.

The intuition behind the introduction of $(𝒪,\mathrm{\Pi },\mathscr{H})$ is that we can define a probability measure over $𝒯$ for each object $o\in 𝒪$, but only some of these objects correspond to a Markovian flow with the right terminating flows. For such objects $o$ (i.e. those that can be written as $o=\mathscr{H}(F)$ for some flow $F\in {ℱ}_{Markov}(G,R)$), the probability measure $P_T$ corresponds to the distribution of interest, according to Def. 13, i.e.:

$$\forall s\in {𝒮}^f{P}_T(s)\propto R(s)$$

GFlowNets thus provide a solution to the generally intractable problem of sampling from a target reward function $R$, or its associated energy function:

$$\forall s\in {𝒮}^fℰ(s):=-\log R(s)$$ (35)

Directly approximating flows $F\in {ℱ}_{Markov}(G,R)$ is a hard problem, whereas with some sets $𝒪$, searching for an object $o\in \mathscr{H}({ℱ}_{Markov}(G,R))\subseteq 𝒪$ is a simpler problem that can be tackled with function approximation techniques.

Note that not the set $𝒪$ cannot be arbitrary, as there needs to be a way to define an injective function from ${ℱ}_{Markov}(G,R)$ to $𝒪$. Below, for a given DAG $G$, we show three examples clarifying the abstract concept of parametrization:

Example 1.

Edge-flow parametrization: Consider ${𝒪}_{edge}=ℱ({𝔸}^{-f},{ℝ}^{+})$, the set of functions from ${𝔸}^{-f}$ to ${ℝ}^{+}$, and the functionals ${\mathscr{H}}_{edge}:{ℱ}_{Markov}(G,R)\to {𝒪}_{edge}$ and ${\mathrm{\Pi }}_{edge}:{𝒪}_{edge}\to \mathrm{\Delta }(𝒯)$ defined by:

${\mathscr{H}}_{edge}(F):(s\to s^{\prime })\in {𝔸}^{-f}\mapsto F(s\to s^{\prime }),$

$\forall \tau =(s_0,\mathrm{\dots },s_n=s_f)\in 𝒯{\mathrm{\Pi }}_{edge}(\widehat{F})(\tau )\propto {\displaystyle \prod _{t=1}^n}{P}_{\widehat{F}}(s_t\mid s_{t-1}),$

where

$$P_{\widehat{F}}(s^{\prime }\mid s)=\{\begin{array}{cc}\frac{\widehat{F}(s\to s^{\prime })}{{\sum }_{s^{\prime \prime }\ne s_f}\widehat{F}(s\to s^{\prime \prime })+R(s)}\text{if}{s}^{\prime }\ne s_f\hfill & \\ \frac{R(s)}{{\sum }_{s^{\prime \prime }\ne s_f}\widehat{F}(s\to s^{\prime \prime })+R(s)}\text{if}{s}^{\prime }=s_f\hfill & \end{array}$$ (36)

The injectivity of ${\mathscr{H}}_{edge}$ follows directly from Prop. 23 (two Markovian flows that coincide on both their terminating and non-terminating edge flow values are equal). And for any Markovian flow $F\in {ℱ}_{Markov}(G,R)$, ${\mathrm{\Pi }}_{edge}({\mathscr{H}}_{edge}(F))$ equals the probability measure associated with $F$, as is shown in Prop. 16.

$({𝒪}_{edge},{\mathrm{\Pi }}_{edge},{\mathscr{H}}_{edge})$ is thus a valid flow parametrization of $(G,R)$.

Example 2.

Forward transition probability parametrization: Consider the set ${𝒪}_{PF}={𝒪}_1\times {𝒪}_2$, where ${𝒪}_1=ℱ(𝒮\setminus \{s_f\},{ℝ}^{+})$ is the set of function from $𝒮\setminus \{s_f\}$ to ${ℝ}^{+}$ and ${𝒪}_2$ is the set of forward probability functions ${\widehat{P}}_F$ consistent with $G$ , and the functionals ${\mathscr{H}}_{PF}:{ℱ}_{Markov}(G,R)\to {𝒪}_{PF}$ and ${\mathrm{\Pi }}_{PF}:{𝒪}_{PF}\to \mathrm{\Delta }(𝒯)$ defined by:

${\mathscr{H}}_{PF}(F)=(s\in 𝒮\setminus \{s_f\}\mapsto F(s),(s\to s^{\prime })\in 𝔸\mapsto P_F(s^{\prime }\mid s)),$

$\forall \tau =(s_0,\mathrm{\dots },s_n=s_f)\in 𝒯{\mathrm{\Pi }}_{PF}(\widehat{F},{\widehat{P}}_F)(\tau )\propto {\displaystyle \prod _{t=1}^n}{\widehat{P}}_F(s_t\mid s_{t-1}),$

where $P_F$ is the forward transition probability function associated with $F$ (Eq. 16). To verify that ${\mathscr{H}}_{PF}$ is injective, consider $F_1,F_2\in {ℱ}_{Markov}(G,R)$ such that ${\mathscr{H}}_{PF}(F_1)={\mathscr{H}}_{PF}(F_2)$. It means that $\forall s\in {𝒮}^f$, $F_1(s)=F_2(s)$, and $\forall s\to s^{\prime }\in 𝔸$, $\frac{F_1(s\to s^{\prime })}{F_1(s)}=\frac{F_2(s\to s^{\prime })}{F_2(s)}$. It follows that $\forall s\to s^{\prime }\in 𝔸$, $F_1(s\to s^{\prime })=F_2(s\to s^{\prime })$. Which, according to Prop. 23, means that $F_1=F_2$. And for any Markovian flow $F\in {ℱ}_{Markov}(G,R)$, ${\mathrm{\Pi }}_{PF}({\mathscr{H}}_{PF}(F))$ equals the probability measure associated with $F$, as is shown in Prop. 16.

$({𝒪}_{PF},{\mathrm{\Pi }}_{PF},{\mathscr{H}}_{PF})$ is thus a valid flow parametrization of $(G,R)$.

Example 3.

Transition probabilities parametrization: Similar to Ex. 2, we can parametrize a Markovian flow using the state-flow function and both its forward and backward transition probabilities, i.e. with ${𝒪}_{PFB}={𝒪}_{PF}\times {𝒪}_3$, ${\mathscr{H}}_{PFB}$, and ${\mathrm{\Pi }}_{PFB}$ defined as:

${\mathscr{H}}_{PFB}(F)=({\mathscr{H}}_{PF}(F),(s\to s^{\prime })\in {𝔸}^{-f}\mapsto P_B(s\mid s^{\prime }),),$

$\forall \tau =(s_0,\mathrm{\dots },s_n=s_f)\in 𝒯{\mathrm{\Pi }}_{PFB}(\widehat{F},{\widehat{P}}_F,{\widehat{P}}_B)(\tau )\propto {\displaystyle \prod _{t=1}^n}{\widehat{P}}_F(s_t\mid s_{t-1}),$

where $P_B$ is the function defined by Eq. 17. and ${𝒪}_3$ is the set of backward probability functions ${\widehat{P}}_B$ consistent with $G$. The injectivity of ${\mathscr{H}}_{PFB}$ is a direct consequence of that of ${\mathscr{H}}_{PF}$. And for any Markovian flow $F$, ${\mathrm{\Pi }}_{PFB}({\mathscr{H}}_{PFB}(F))$ equals the probability measure associated with $F$, as is shown in Prop.3.

$({𝒪}_{PFB},{\mathrm{\Pi }}_{PFB},{\mathscr{H}}_{PFB})$ is thus a valid flow parametrization of $(G,R)$.

We now have all the ingredients to formally define a GFlowNet:

Definition 25.

A GFlowNet is a tuple $(G,R,𝒪,\mathrm{\Pi },\mathscr{H})$, where:

• $G=(𝒮,𝔸)$ is a pointed DAG with initial state $s_0$ and sink state $s_f$,

• $R:{𝒮}^f\to {ℝ}^{+}$ a target reward function,

• $(𝒪,\mathrm{\Pi },\mathscr{H})$ a flow parametrization of $(G,R)$.

Each object $o\in 𝒪$ is called a GFlowNet configuration. When it is clear from context, we will use the term GFlowNet to refer to both $(G,R,𝒪,\mathrm{\Pi },\mathscr{H})$ and a particular configuration $o$; similar to how the term “Neural Network” refers to both the class of functions that can be represented with a particular architecture, and to a particular element of that class / weight configuration.

If $o\in \mathscr{H}({ℱ}_{Markov}(G,R))$, then the corresponding terminating state probability measure (Eq. 34) is proportional to the target reward $R$.

Once we have a GFlowNet $(G,R,𝒪,\mathrm{\Pi },\mathscr{H})$, we still need a way to find objects $o\in \mathscr{H}({ℱ}_{Markov}(G,R))\subseteq 𝒪$. To this end, it suffices to design a loss function $ℒ$ on $𝒪$ that equals zero on objects $o\in \mathscr{H}({ℱ}_{Markov}(G,R))$ and only on those objects. If our loss function $ℒ$ is chosen to be non-negative, then an approximation of the target distribution (on ${𝒮}^f$) is obtained by approximating the minimum of the function $ℒ$. This provides a recipe for casting the search problem of interest to a minimization problem, as we typically do in machine learning. Such loss functions can be easily designed for the natural parametrizations we considered in Ex. 1, Ex. 2, and Ex. 3, as we will illustrate below.

Definition 26.

Let $(G,R,𝒪,\mathrm{\Pi },\mathscr{H})$ be a GFlowNet. A flow-matching loss is any function $ℒ:𝒪\to {ℝ}^{+}$ such that:

$\forall o\in 𝒪ℒ(o)=0\iff \exists F\in {ℱ}_{Markov}(G,R)o=\mathscr{H}(F)$ (37)

We say that $ℒ$ is edge-decomposable, if there exists a function $L:𝒪\times 𝔸\to {ℝ}^{+}$ such that:

$$\forall o\in 𝒪ℒ(o)=\sum _{s\to s^{\prime }\in 𝔸}L(o,s\to s^{\prime }),$$

We say that $ℒ$ is state-decomposable, f there exists a function $L:𝒪\times 𝒮\to {ℝ}^{+}$ such that:

$$\forall o\in 𝒪ℒ(o)=\sum _{s\in 𝒮}L(o,s),$$

We say that $ℒ$ is trajectory-decomposable if there exists a function $L:𝒪\times 𝒯\to {ℝ}^{+}$ such that:

$$\forall o\in 𝒪ℒ(o)=\sum _{\tau \in 𝒯}L(o,\tau )$$

As mentioned above, with a such a loss function, our search problems can be written as minimization problems of the form

$$\underset{o\in 𝒪}{\min }ℒ(o),$$ (38)

which can be tackled with gradient-based learning if the function $ℒ$ is differentiable. Note that with an edge-decomposable flow-matching loss, the minimization problem in Eq. 38 is equivalent to:

$$\underset{o\in 𝒪}{\min }{𝔼}_{(s\to s^{\prime })\sim {\pi }_T}[L(o,s\to s^{\prime })],$$ (39)

where ${\pi }_T$ is any full support probability distribution on $𝔸$, i.e. a probability distribution such that $\forall s\to s^{\prime }\in 𝔸{\pi }_T(s\to s^{\prime })>0$. A similar statement can be made for state-decomposable or trajectory-decomposable flow-matching losses.

Example 4.

Consider the edge-flow parametrization $({𝒪}_{edge},{\mathrm{\Pi }}_{edge},{\mathscr{H}}_{edge})$, and the function $L_{FM}:{𝒪}_{edge}\times 𝒮\to {ℝ}^{+}$ defined for each $\widehat{F}\in {𝒪}_{edge}$ and $s^{\prime }\in 𝒮$ as

$L_{FM}(\widehat{F},s^{\prime })=\{\begin{array}{cc}{\left(\log \left(\frac{\delta +{\sum }_{s\in Par(s^{\prime })}\widehat{F}(s\to s^{\prime })}{\delta +R(s^{\prime })+{\sum }_{s^{\prime \prime }\in Child(s^{\prime })\setminus \{s_f\}}\widehat{F}(s^{\prime }\to s^{\prime \prime })}\right)\right)}^2\text{if}{s}^{\prime }\ne s_f,\hfill & \\ 0\text{otherwise}\hfill & \end{array}$

where $\delta \ge 0$ is a hyper-parameter. The function ${ℒ}_{FM}$ mapping each $\widehat{F}\in {𝒪}_{edge}$ to

$${ℒ}_{FM}(\widehat{F})=\sum _{s\in 𝒮}{L}_{FM}(\widehat{F},s)$$ (40)

is a flow-matching loss, that is (by definition) state-decomposable.

To see this, let $\widehat{F}\in {𝒪}_{edge}$ such that ${ℒ}_{FM}(\widehat{F})=0$, and extend it to terminating edge:

$$\forall s\in {𝒮}^f\widehat{F}(s\to s_f):=R(s)$$

Now that $\widehat{F}$ is defined for all edges in $G$, we can write that

$$\forall s^{\prime }\in 𝒮\sum _{s\in Par(s^{\prime })}\widehat{F}(s\to s^{\prime })=\sum _{s^{\prime \prime }\in Child(s)}\widehat{F}(s^{\prime }\to s^{\prime \prime }).$$

Which, according to Prop. 19, means that there exists a Markovian flow $F\in {ℱ}_{Markov}(G,R)$ such that ${\mathscr{H}}_{edge}(F)=\widehat{F}$. The converse

$$\forall F\in {ℱ}_{Markov}(G,R){ℒ}_{FM}({\mathscr{H}}_{edge}(F))=0$$

is a trivial consequence of Prop. 19.

This is the loss function proposed in Bengio et al. (2021). $\delta$ allows to reduce the importance given to small flows (those smaller than $\delta$), and the usage of the square of the log-ratio is justified as a way to ensure that states with large flows do not contribute to the gradients of ${ℒ}_{FM}$ much more than states with small flows.

Example 5.

Detailed-balance loss: Consider the transition probabilities parametrization $({𝒪}_{PFB},{\mathrm{\Pi }}_{PFB},{\mathscr{H}}_{PFB})$, and the function $L_{DB}:{𝒪}_{PFB}\times 𝔸\to {ℝ}^{+}$ defined for each $(\widehat{F},{\widehat{P}}_F,{\widehat{P}}_B)\in {𝒪}_{PFB}$ and $s\to s^{\prime }\in 𝔸$ as

$L_{DB}(\widehat{F},{\widehat{P}}_F,{\widehat{P}}_B,s\to s^{\prime })=\{\begin{array}{cc}{\left(\log \left(\frac{\delta +\widehat{F}(s){\widehat{P}}_F(s^{\prime }\mid s)}{\delta +\widehat{F}(s^{\prime }){\widehat{P}}_B(s\mid s^{\prime })}\right)\right)}^2\text{if}{s}^{\prime }\ne s_f,\hfill & \\ {\left(\log \left(\frac{\delta +\widehat{F}(s){\widehat{P}}_F(s^{\prime }\mid s)}{\delta +R(s)}\right)\right)}^2\text{otherwise},\hfill & \end{array}$

where $\delta \ge 0$ is a hyper-parameter. The function ${ℒ}_{DB}$ mapping each $(\widehat{F},{\widehat{P}}_F,{\widehat{P}}_B)\in {𝒪}_{PFB}$ to

$${ℒ}_{DB}(\widehat{F},\widehat{P},{\widehat{P}}_B))=\sum {}_{s\to s^{\prime }\in 𝔸}L{}_{DB}(\widehat{F},\widehat{P},{\widehat{P}}_B,s\to s^{\prime })$$

is a flow-matching loss that is (by definition) edge-decomposable. The proof of this statement is similar to the one of the example above, using Prop. 21.

According to Sec. 2.6, the reward function does not completely specify the flow. Thus, the detailed-balance loss of Ex. 5 can be used with the $({𝒪}_{PF},{\mathrm{\Pi }}_{PF},{\mathscr{H}}_{PF})$ parametrization, using any function ${\widehat{P}}_B\in {𝒪}_3$ as input to the detailed-balance loss.

Example 6.

Trajectory-balance loss: This loss has been introduced in Malkin et al. (2022) for the parametrization $({𝒪}_{TB},{\mathrm{\Pi }}_{TB},{\mathscr{H}}_{TB})$, where ${𝒪}_{TB}={𝒪}_1\times {𝒪}_2\times {𝒪}_3$, with ${𝒪}_1={ℝ}^{+}$ parametrizes the partition function $\widehat{Z}$, and ${𝒪}_2$ and ${𝒪}_3$ introduced in Ex. 2 and Ex. 3 (the set of forward and backward probabilities consistent with $G$). ${\mathscr{H}}_{TB}$ maps a Markovian flow in ${ℱ}_{Markov}(G,R)$ to the corresponding triplet $(Z,P_F,P_B)$, and ${\mathrm{\Pi }}_{TB}$ maps a parametrization $(\widehat{Z},{\widehat{P}}_F,{\widehat{P}}_B)$ to a probability over trajectories defined by ${\widehat{P}}_F$ as in Ex. 2. Prop. 18 justifies the validity of this parametrization. The loss ${ℒ}_{TB}$ maps each $(\widehat{Z},{\widehat{P}}_F,{\widehat{P}}_B)\in {𝒪}_{TB}$ to:

$${ℒ}_{TB}(\widehat{Z},{\widehat{P}}_F,{\widehat{P}}_B)=\sum _{\tau \in 𝒯}{L}_{TB}(\widehat{Z},{\widehat{P}}_F,{\widehat{P}}_B,\tau ),$$

where

$$\forall \tau =(s_0,\mathrm{\dots },s_{n+1}=s_f)\in 𝒯L_{TB}(\widehat{Z},{\widehat{P}}_F,{\widehat{P}}_B,\tau )={\left(\log \frac{\widehat{Z}{\prod }_{t=1}^{n+1}{\widehat{P}}_F(s_t\mid s_{t-1})}{R(s_n){\prod }_{t=1}^n{\widehat{P}}_B(s_{t-1}\mid s_t)}\right)}^2.$$ (41)

Malkin et al. (2022) prove that ${ℒ}_{TB}$ is a flow-matching loss and call it trajectory balance. It is trajectory-decomposable by definition.

Training by stochastic gradient descent:

In the examples of the previous section, given a GFlowNet $(G,R,𝒪,\mathrm{\Pi },\mathscr{H})$ and a flow-matching loss $ℒ$, objects $o\in 𝒪$ are themselves functions or combinations of functions, and we can thus parametrize $𝒪$ with function approximators such as Neural Networks. However, most of the times, the evaluation (let alone the minimization) of $ℒ(o)$ is intractable, given that even with a full support distribution, only a subset of edges (or states or trajectories) can be visited in finite time. In practice, with an edge-decomposable loss e.g., we resort to a stochastic gradient, such as

$${\nabla }_{o}L(o,s\to s^{\prime }),s\to s^{\prime }\sim {\pi }_o$$ (42)

for edge-decomposable losses, or

$${\nabla }_{o}L(o,\tau ),\tau \sim {\pi }_o$$ (43)

for trajectory-decomposable losses, where ${\pi }_o$, called the training distribution, is a distribution over edges or trajectories that can be associated with $\mathrm{\Pi }(o)$, corresponding to the online setting in RL, or defined in other ways, corresponding to the behavior policy in offline RL, see Sec. 3.3.3 below.

3.3 Extensions

In this section, we discuss possible relaxations to the GFlowNet training paradigm introduced thus far.

3.3.1 Introducing Time Stamps to Allow Cycles

Note that the state-space of a GFlowNet can easily be modified to accommodate an underlying state space for which the transitions do not form a DAG, e.g., to allow cycles. Let $𝒮$ be such an underlying state-space. Define the augmented state space ${𝒮}^{\prime }=𝒮\times \mathbb{N}$, where $\mathbb{N}=\{0,1,2,\mathrm{\dots }\}$ is the set of natural numbers, and $s_t^{\prime }=(s_t,t)$ is the augmented state, where $t$ is the position of the state $s_t$ in the trajectory. With this augmented state space, we automatically avoid cycles. Furthermore, we may design or train the backwards transition probabilities $P_B(s_t^{\prime }\mid s_{t+1}^{\prime }=(s_{t+1},t+1))$ to create a preference for shorter paths towards $s_{t+1}$, as discussed in Sec. 2.6. Note that we can further generalize this setup by replacing $\mathbb{N}$ with any totally ordered indexing set; the augmented state space will still have an associated DAG. The ordering “” in the original state-space is lifted to the augmented state-space: if and only if and .

3.3.2 Stochastic Rewards

We also consider the setting in which the given reward is stochastic rather than being a deterministic function of the state, yielding training procedures based on stochastic gradient descent. For example, with the trajectory balance loss of Eq. 41, if $R(s)$ is stochastic (even when given $s$), we can think of what is being really optimized is the squared loss with $\log R(s)$ replaced by its expectation (given $s$). This is a straightforward consequence of minimizing the expected value of a squared error loss (as for example in neural networks trained with a squared error loss and a stochastic target output, where the neural network effectively tries to estimate the expected value of that target).

3.3.3 GFlowNets can be trained offline

As discussed in Sec. 3.2, we do not need to train a GFlowNet using samples from its own trajectory distribution $\widehat{P}=\mathrm{\Pi }(o)$. Those training trajectories can be drawn from any training distribution ${\pi }_T$ with full support, as already shown by Bengio et al. (2021). It means that a GFlowNet can be trained offline, as in offline reinforcement learning (Ernst et al., 2005; Riedmiller, 2005; Lange et al., 2012).

It should also be noted that with a proper adaptive choice of ${\pi }_T$, and assuming that computing $R$ is cheaper or comparable in cost to running the GFlowNet on a trajectory, it should be more efficient to continuously draw new training samples from ${\pi }_T$ than to rehearse the same trajectories multiple times. An exception would be rehearsing the trajectories leading to high rewards if these are rare.

How should one choose the training distribution ${\pi }_T$? It needs to cover the support of $R$ but if it were uniform it would be very wasteful and if it were equal to the current GFlowNet policy $\pi$ it might not have sufficient effective support and thus miss modes of $R$, i.e., regions where $R(x)$ is substantially greater than 0 but $R(x)\gg P_T(x)$. Hence the training distribution should be sampled from an exploratory policy that visits places that have not been visited yet and may have a high reward. High epistemic uncertainty around the current policy would make sense and the literature on acquisition functions for Bayesian optimization (Srinivas et al., 2010) may be a good guide. More generally, this means the training distribution should be adaptive. For example, ${\pi }_T$ could be the policy of a second GFlowNet trained mostly to match a different reward function that is high when the losses observed by the main GFlowNet are large. It would also be good to regularly visit those trajectories corresponding to known large $R$, i.e., according to samples from $\pi$, to make sure those are not forgotten, even temporarily.

3.4 Exploiting Data as Known Terminating States

In some applications we may have access to a dataset of $(s,R(s))$ pairs and we would like to use them in a purely offline way to train a GFlowNet, or we may want to combine such data with queries of the reward function $R$ to train the GFlowNet. For example, the dataset may contain examples of some of the high-reward terminating states $s$ which would be difficult to obtain by sampling from a randomly initialized GFlowNet. How can we compute a gradient update for the GFlowNet parameters using such $(s,R(s))$ pairs?

If we choose to parametrize the backwards transition probabilities $P_B$ (which is necessary for implementing the detailed balance loss), then we can just sample a trajectory $\tau$ leading to $s$ using $P_B$ and use these trajectories to update the flows and forward transition probabilities along the traversed transitions. However, this alone is not guaranteed to produce the correct GFlowNet sampling distribution because the empirical distribution over training trajectories $\tau$ defined as above does not have full support. Suppose for example that the dataset only contains high-reward terminating states with $R(s)=1$. The GFlowNet could then just sample trajectories uniformly (which would be wrong, we would like the probability of most states not in the training set to be very small). On the other hand, if we combine the distribution of trajectories leading to terminal transitions in the dataset with a training distribution whose support covers all possible trajectories, then the offline property of GFlowNet guarantees that we can recover a flow-matching model.

4 Conditional Flows and Free energies

A remarkable property of flow networks is that we can recover the normalizing constant $Z$ from the initial state flow $F(s_0)$ (Prop. 10). $Z$ also gives us the partition function associated with a given terminal reward function $R$ specifying the terminating flows.

What about internal states $s$ with ? If we had something like a normalizing constant for only the terminating flows achievable from $s$, we would be able to obtain a form of marginalization given state $s$, i.e., a conditional probability for terminating states $s^{\prime }\ge s$, given $s$. Naturally, one could ask: does the flow $F(s)$ through state $s$ give us that kind of marginalization over only the downstream terminating flows? Unfortunately in general, the answer to this question is no, as illustrated in Fig. 5: in this example $F(s_2)=4$, whereas the sum of terminating flows achievable from $s_2$ is $6$ (the terminating states reachable from $s_2$ are $\{s_5,s_6,s_7\}$). The discrepancy is caused by the flow through $(s_0,s_1,s_5)$ that contributes to the terminating flow $F(s_5\to s_f)$, but not to $F(s_2)$ since there is no order relation between $s_1$ and $s_2$.

(a)

(b)

(c)

Figure 5: Example of a state-conditional flow network. (a) The original (Markovian) flow network. (b) The subgraph of states reachable from $s_2$; there is a flow through $(s_0,s_1,s_5)$ that contributed to $F(s_5\to s_f)$, but not to $F(s_2)$, showing that $F(s_2)$ does not marginalize the rewards of its descendant. (c) State-conditional flow network $F_{s_2}$, which differs from the original flow $F$ on the subgraph, but satisfies the desired marginalization property.

In Sec. 5.3, we show how GFlowNets applied to sampling sets of random variables can be used to estimate the marginal probability for the values given to a subset of the variables. It requires computing the kind of intractable sum discussed above (over the rewards associated with all the descendants of a state $s$, with $s$ corresponding to such a subset of variables and a descendant to a full specification of all the variables). That motivates the following definition:

Definition 27.

Given a pointed DAG $G=(𝒮,𝔸)$, the corresponding partial order denoted by $\ge$, and a function $ℰ:𝒮\to ℝ$, called the energy function, we define the free energy $ℱ(s)$ of a state $s$ as:

$$e^{-ℱ(s)}:=\sum _{s^{\prime }:s^{\prime }\ge s}{e}^{-ℰ(s^{\prime })}.$$ (44)

Free energies are generic formulations for the marginalization operation (i.e. summing over a large number of terms) associated with energy functions, and we find their estimation to open the door to interesting applications where expensive MCMC methods would typically be the main approach otherwise.

4.1 Conditional flow networks

In Sec. 2.2, we defined a flow network as a DAG, augmented with some function $F$ over the set of complete trajectories $𝒯$. We can extend this notion of flow networks by conditioning each component on some information $x$. In general, this conditioning variable can represent any conditioning information, either external to the flow network (but influencing the terminating flows), or internal (e.g., $x$ can be a property of complete trajectories over another flow network, like passing through a particular state).

Definition 28.

Let $𝒳$ be a set of conditioning variables. We consider a family of DAGs $G_x=({𝒮}_x,{𝒜}_x)$ indexed by $x\in 𝒳$, along with a family of initial and terminal states denoted by $(s_0\mid x)\in {𝒮}_x$ and $(s_f\mid x)\in {𝒮}_x$ respectively. For each DAG $G_x$, we denote by ${𝒯}_x$ the set of complete trajectories in $G_x$, and we denote by $𝒯$ their union:

$$𝒯=\bigcup _{x\in 𝒳}{𝒯}_x.$$

A conditional flow network is the specification of $𝒳$, the family $\{G_x,x\in 𝒳\}$, along with a conditional flow function $F$, i.e. a function $F:𝒳\times 𝒯\to {ℝ}^{+}$ such that $F(x,\tau )=0$ if $\tau \notin {𝒯}_x$. For clarity, we will denote, for each $x\in 𝒳$, by $F_x$ the function mapping each $\tau \in {𝒯}_x$ to $F(x,\tau )$. Similar to Sec. 2.2, $F_x$ induces a measure of the $\sigma$-algebra $2^{{𝒯}_x}$ for each $x$.

Conditional flow networks effectively represent a family of flow networks, indexed by the value of $x$. Since conditional flow networks are defined using the same components as an unconditional flow network, they inherit from all the properties of flow networks for all DAGs $G_x$ and flow functions $F_x$. In particular, we can directly extend the notion of probability distribution over flows, state and edge flows, forward and backward transition probabilities (Sec. 2.3), of Markovian flows (Sec. 2.4), and any flow matching condition (Sec. 2.5) to conditional flows; the only difference is that now every term explicitly depends of the conditioning variable $x$.

In Sec. 4.2 and Sec. 4.3, we will elaborate two important examples of conditional flow networks: flow networks conditioned on external information that changes the reward $R(s\mid x)$, and state-conditional flow networks that depend on internal information, i.e., previously visited states.

4.2 Reward-conditional flow networks

Definition 29.

Let $𝒳$ be a set of conditioning variables. Consider a flow network given by a pointed DAG $G=(𝒮,𝔸)$ and a flow function $F$. Consider a family $ℛ$ of non-negative functions of $𝒮$: $\{R_x:{𝒮}^f\to {ℝ}^{+},x\in 𝒳\}$. A reward-conditional flow network compatible with the family $ℛ$ is a conditional flow network (Def. 28), with $G_x=G$ for every $x\in 𝒳$, such that the edge-flow functions induced by the conditional flow function $F$ satisfy:

$\forall x\in 𝒳\forall s\in {𝒮}^f{F}_x(s\to s_f)=R_x(s).$

We will use the notations $R_x(s)$ and $R(s\mid x)$ interchangeably.

Note that the definition above implies that all the DAGs of a reward-conditional flow network are identical, and only the terminating flows differ amongst the members of the family.

Example 7.

We will see in Sec. 4.4 that we can estimate a conditional flow network using a GFlowNet (Sec. 3), given a reward function $R(s\mid x)$. In an Energy-Based Model, the model $P_{\theta }(s)$ is associated with a given energy function ${ℰ}_{\theta }(s)$, parametrized by $\theta$, with

$$P_{\theta }(s)=\frac{\exp (-{ℰ}_{\theta }(s))}{Z(\theta )}.$$

This model can be parametrized using a reward-conditional flow network, conditioned on $\theta$ with the reward function $R(s\mid \theta )=\exp (-{ℰ}_{\theta }(s))$. We show in Sec. 4.5 how to use such a conditional flow-network to learn an Energy-Based Model.

4.3 State-conditional flow networks

Definition 30.

Consider a flow network given by a DAG $G=(𝒮,𝔸)$ and a flow function $F$. For each state $s\in 𝒮$, let $G_s$ be the subgraph of $G$ containing all the states $s^{\prime }$ such that $s^{\prime }\ge s$. A state-conditional flow network is given by the family $\{G_s,s\in 𝒮\}$, along with a conditional flow function $F:𝒮\times 𝒯\to {ℝ}^{+}$, where $𝒯={\bigcup }_{s\in 𝒮}{𝒯}_s$, and ${𝒯}_s$ the set of complete trajectories in ${𝒢}_s$, that satisfies:

$$F_s(s^{\prime }\to s_f)=F(s^{\prime }\to s_f).$$ (45)

Note that in the definition above, we abused the notation $F$ to refer to both flow functions and edge flow functions, but also used $F_s$ to refer to the conditional flow function (or the corresponding edge flow function) $\tau \mapsto F(s,\tau )$. Unlike the reward-conditional flow networks defined in Sec. 4.2, the structure of the DAG in a state-conditional flow network depends on the anchor state $s$. In particular, this means that the initial state $(s_0\mid s)=s$ changes, but the final state $(s_f\mid s)=s_f$ remains unchanged, for any state $s$.

Since the definition of a state-conditional flow network depends on an original flow network, we must ensure that this definition is indeed correct, i.e. that such a state-conditional flow network that satisfies the conditions in Eq. 45 exists.

Proposition 31.

For any flow network given by a DAG $G=(𝒮,𝔸)$ and a flow F, we can define a state-conditional flow network as per Def. 30.

Proof Let $s\in 𝒮$ be a state. Since the structure of the DAG $G_s$ is clearly well-defined, we just need to show that there exists a flow function $F_s:{𝒯}_s\to {ℝ}^{+}$ that satisfies Eq. 45. If such a function exists for every $s\in 𝒮$, then it would suffice to define the conditional flow function $F:𝒮\times 𝒯\to {ℝ}^{+}$ as:

$$F(s,\tau )=\{\begin{array}{cc}{F}_s(\tau )\text{if}\tau \in {𝒯}_s\hfill & \\ 0\text{otherwise}.\hfill & \end{array}$$

Let $A_{s^{\prime }\mid s}$ be the set of complete trajectories in ${𝒯}_s$ terminating in $s^{\prime }\ge s$; the condition in Eq. 45 then reads:

$$F_s(s^{\prime }\to s_f)=F_s(A_{s^{\prime }\mid s})=\sum _{\tau \in A_{s^{\prime }\mid s}}{F}_s(\tau )=F(s^{\prime }\to s_f).$$ (46)

Note that in Eq. 46, $F(s^{\prime }\to s_f)$ is a given quantity because the flow $F$ is known. Since the sets of trajectories $\{A_{s^{\prime }\mid s},,s^{\prime }\ge s\}$ form a partition of all the complete trajectories ${𝒯}_s$, Eq. 46 is a system of linear equations, whose unknowns are $F_s(\tau )$ for all $\tau \in {𝒯}_s$, where each equation involves separate sets of unknowns. Therefore there exists at least a solution $F_s(\tau )$ of this system.

We can construct such a solution in the following way. For some $\tau \in {𝒯}_s$, we can first start by selecting the complete trajectories $\overline{\tau }\in 𝒯$ that contain $\tau$:

$$C_{\tau }=\{\overline{\tau }\in 𝒯:\tau \subseteq \overline{\tau }\}.$$

The key difference between the DAG $G$ and the subgraph $G_s$ though is that $G$ may contain trajectories that terminate in some $s^{\prime }\ge s$ but do not pass through $s$, and those are therefore not covered by the trajectories of $G_s$. Let $U_{s^{\prime }\mid s}$ be the set of complete trajectories of $G$ defined as

$$U_{s^{\prime }\mid s}=\{\overline{\tau }\in 𝒯:\exists s^{\prime \prime }>s,s^{\prime \prime }\in \overline{\tau },s^{\prime }\in \overline{\tau },s\notin \overline{\tau }\}.$$

For all $\tau \in {𝒯}_s$ such that $\tau$ terminates in some $s^{\prime }\ge s$, we can therefore construct the flow $F_s(\tau )$ as

$$F_s(\tau ):=F(C_{\tau })+\frac{1}{n}F(U_{s^{\prime }\mid s}),$$

where $n=|A_{s^{\prime }\mid s}|$ is the number of trajectories ${\tau }^{\prime }\in {𝒯}_s$ that terminate in $s^{\prime }$. It is easy to verify that $F_s(\tau )$ is a solution of Eq. 46.

While we saw in Prop. 10 that the initial flow $F(s_0)$ was equal to the partition function, the initial state-conditional flow also benefit from a marginalizing property, and is now related to the free energy at $s$.

Proposition 32.

Given a state-conditional flow network $(G,F)$ as in Def. 30, for any state $s$, the initial flow of the state-conditional flow network corresponds to marginalizing the terminating flows $F(s^{\prime }\to s_f)$ for $s^{\prime }\ge s$:

$$F_s(s_0\mid s)=F_s(s)=\sum _{s^{\prime }:s^{\prime }\ge s}F(s^{\prime }\to s_f)=\exp (-ℱ(s)),$$

where $ℱ(s)$ is the free energy associated to the energy function $ℰ(s^{\prime })=-\log F(s^{\prime }\to s_f)$.

Proof This is a direct consequence of Prop. 10, applied to the state-conditional flow function ${ℱ}_s$, along with Def. 27.

$F_s(s_0\mid s)$ $={\displaystyle \sum _{\tau \in {𝒯}_s}}{F}_s(\tau )={\displaystyle \sum _{s^{\prime }:s^{\prime }\ge s}}{F}_s(s^{\prime }\to s_f)$

$={\displaystyle \sum _{s^{\prime }:s^{\prime }\ge s}}F(s^{\prime }\to s_f)={\displaystyle \sum _{s^{\prime }:s^{\prime }\ge s}}{e}^{-ℰ(s^{\prime })}$

Note that the definition of state-conditional flow networks is consistent with our original definition of (unconditional) flow networks in Section 2.1, in the sense that the original flow network is a valid state-conditional flow networks anchored at the initial state.

Another quantity of interest that state-conditional flow networks allow us to evaluate, is the probability of terminating a trajectory in a state $s^{\prime }$ if all terminating edge flows were diverted towards an earlier state :

Corollary 33.

Consider a flow network given by a DAG $G=(𝒮,𝔸)$ and a flow $F$, from which we define any state-conditional flow network, as per Def. 30. Given a state $s$, the flow function $F_s$ induces a probability distribution over $\{s^{\prime \prime }\in {𝒮}^f:s^{\prime \prime }\ge s\}\subseteq {𝒮}^f$, that we denote by $P_T(.\mid s)$.

Under this measure, the probability of terminating a trajectory in $G_s$ in a state $s^{\prime }$ (i.e. the last edge of the trajectory is $s^{\prime }\to s_f$) is:

$$P_T(s^{\prime }\mid s)={\mathrm{𝟏}}_{s^{\prime }\ge s}{e}^{-ℰ(s^{\prime })+ℱ(s)},$$ (47)

where $ℰ$ is the energy function mapping each state $s^{\prime }$ that is parent of $s_f$ to $-\log F(s^{\prime }\to s)$, and $ℱ$ is the corresponding free energy function.

Proof Because $F_s$ is a flow function, Def. 11 and Prop. 10 tell us that:

Combining this with Prop. 32, and Eq. 45, we obtain for $s^{\prime }\ge s$:

$P_T(s^{\prime }\mid s)$ $={\mathrm{𝟏}}_{s^{\prime }\ge s}{\displaystyle \frac{F(s^{\prime }\to s_f)}{e^{-ℱ(s)}}}$

$={\mathrm{𝟏}}_{s^{\prime }\ge s}{\displaystyle \frac{e^{-ℰ(s^{\prime })}}{e^{-ℱ(s)}}}$

4.4 Conditional GFlowNets

Similar to the way we used a GFlowNet to estimate the flow of a flow network, we can also use a (conditional) GFlowNet in order to estimate a conditional flow network, with given target reward functions. A conditional GFlowNet follows the construction presented in Section 3, with the exception that all quantities to be learned now depend on the conditioning variable $x\in 𝒳$ (e.g., $x$ is an additional input of the neural network).

All parametrizations and losses presented in Sec. 3.2 could in principle be used to train a conditional GFlowNet, regardless of the conditioning set. Below we discuss yet another loss, first presented in Deleu et al. (2022), that could be used to train both GFlowNets and conditional GFlowNets.

Example 8.

Given a family of DAGs $G_x$ and reward functions $R_x$ indexed by $x\in 𝒳$, where each state $s\in G_x$ is terminating (i.e. is a parent of $s_f$), and following Exs. 2 and 3, we consider a parametrization given by the forward and backward transition probabilities ${𝒪}_P^x={𝒪}_2^x\times {𝒪}_3^x$, where ${𝒪}_2^x$ (resp. ${𝒪}_3^x$) is the set of forward (resp. backward) probability functions ${\widehat{P}}_F$ (resp. ${\widehat{P}}_B$) consistent with $G_x$ for every $x\in 𝒳$, and $({\mathrm{\Pi }}_P^x,{\mathscr{H}}_P^x)$ defined as in Exs. 3 and 3. Each $({𝒪}_P^x,{\mathrm{\Pi }}_P^x,{\mathscr{H}}_P^x)$ is a flow parametrization of $(G_x,R_x)$, which can be trained with an edge-decomposable flow-matching loss, as proved in Deleu et al. (2022), and defined for every $s\to s^{\prime }\in {𝔸}^{-f}$:

${ℒ}_{DB}({\widehat{P}}_F,{\widehat{P}}_B,s\to s^{\prime },x)={\left(\log {\displaystyle \frac{R_x(s^{\prime })P_B(s\mid s^{\prime },x)P_F(s_f\mid s,x)}{R_x(s)P_F(s^{\prime }\mid s,x)P_F(s_f\mid s^{\prime },x)}}\right)}^2$

4.5 Training Energy-Based Models with a GFlowNet

A GFlowNet can be trained to convert an energy function into an approximate corresponding sampler. Thus, it can be used as an alternative to MCMC sampling (Sec. 3.1). Consider the model $P_{\theta }(s)$ associated with a given parametrized energy function ${ℰ}_{\theta }(s)$ with parameters $\theta$: $P_{\theta }(s)=\frac{e^{-{ℰ}_{\theta }(s)}}{Z}$. Sampling from $P_{\theta }(s)$ could be approximated by sampling from the terminating probability distribution $P_T(s)$ of a GFlowNet trained with target terminal reward $R(s)=e^{-{ℰ}_{\theta }(s)}$ (see Eq. 34). In practice, ${\widehat{P}}_T$ would be an estimator for the true $P_{\theta }$ because the GFlowNet training objective is not zeroed (insufficient capacity or finite training time). The GFlowNet samples drawn according to ${\widehat{P}}_T$ could then be used to obtain a stochastic gradient estimator for the negative log-likelihood of observed data $x$ with respect to parameters $\theta$ of an energy function ${ℰ}_{\theta }$:

$$\frac{\partial -\log P_{\theta }(x)}{\partial \theta }=\frac{\partial {ℰ}_{\theta }(x)}{\partial \theta }-\sum _s{P}_{\theta }(s)\frac{\partial {ℰ}_{\theta }(s)}{\partial \theta }.$$ (48)

An approximate stochastic estimator of the second term could thus be obtained by sampling one or more terminating states $s\sim {\widehat{P}}_T(s)$, i.e., from the trained GFlowNet’s sampler. Furthermore, if the GFlowNet’s loss is $0$, i.e. ${\widehat{P}}_T=P_{\theta }$, the gradient estimator would be unbiased.

One could thus potentially jointly train an energy function ${ℰ}_{\theta }$ and a corresponding GFlowNet by alternating updates of $\theta$ using the above equation (with sampling from $P_{\theta }$ replaced by sampling from ${\widehat{P}}_T$) and updates of the GFlowNet using the updated energy function for the target terminal reward.

If we fix $\widehat{F}(s\to s_f)=R(s)$ by construction (which we can do if the reward function is deterministic), then we can parametrize the energy function with the same neural network that computes the flow, since $ℰ(s)=-\log R(s)=-\log \widehat{F}(s\to s_f)$. Hence the same parameters are used for the energy function and for the GFlowNet, which is appealing.

The above strategy for learning jointly an energy function and how to sample from it could be generalized to learning conditional distributions by using a conditional GFlowNet instead. Let $x$ be an observed random variable and $h$ be a hidden variable, with the GFlowNet generating the pair $(x,h)$ in two sub-trajectories: either first generate $x$ and then generate $h$ given $x$, or first generate $h$ and then generate $x$ given $h$. This can be achieved by introducing a 6-valued component $u$ in the state to make sure that both $h$ and $x$ are generated before exiting into $s_f$, with the following values and constraints:

$s=s_0$ $\Rightarrow u=0$ (49)

$s=s_f$ $\Rightarrow u=5$ (50)

$(u_t\to u_{t+1})$ $\in \{0\to 1,0\to 2,1\to 1,2\to 2,1\to 3,2\to 4,3\to 3,4\to 4,3\to 5,4\to 5\}$ (51)

where $u=1$ indicates that $x$ is being generated (before $h$), $u=2$ indicates that $h$ is being generated (before $x$), $u=3$ indicates that $h$ is being generated (conditioned on $x$), and $u=4$ indicates that $x$ is being generated (given $h$). The GFlowNet cannot reach the final state $s_f$ until both $x$ and $h$ have been generated. The conditional GFlowNet can thus approximately sample $P_T(x)$, $P_T(h\mid x)$, $P_T(h)$, $P_T(x\mid h)$ as well as $P_T(x,h)$. If we only want to sample $x$ (or only $h$), we allow exiting as soon as it is generated (resp. $h$ is generated). See Sec. 5.3 for a more general discussion on how to represent, estimate and sample marginal distributions.

Let us denote by $P_{\theta }(x,h)\propto e^{-{ℰ}_{\theta }(x,h)}$ the joint distribution over $(x,h)$ associated with the energy function, i.e., with $F(s\to s_f)$. When $x$ is observed but $h$ is not, $\theta$ could thus be updated by approximating the marginal log-likelihood gradient

$$\frac{\partial -\log P_{\theta }(x)}{\partial \theta }=\sum _h{P}_T(h\mid x)\frac{\partial {ℰ}_{\theta }(x,h)}{\partial \theta }-\sum _{x^{\prime },h}{P}_T(x^{\prime },h)\frac{\partial {ℰ}_{\theta }(x^{\prime },h)}{\partial \theta }$$ (52)

using samples from the estimated terminal sampling probabilities ${\widehat{P}}_T$ of a trained GFlowNet to approximate in a stochastic gradient way the above sums (using one or a batch of samples).

Note how we now have outer loop updates (of the energy function, i.e., the reward function) from actual data, and an inner loop updates (of the GFlowNet) using the energy function as a driving target for the GFlowNet. How many inner loop updates are necessary for such a scheme to work is an interesting open question but most likely depends on the form of the underlying data generating distribution. If the work on GANs (Goodfellow et al., 2014) is a good analogy, a good strategy may be to interleave updates of the energy function (as minus the log-terminal flow of a GFlowNet) based on a batch of data, and updates of the GFlowNet as a sampler based on both these samples (trajectories can be sampled backwards from a terminating state $s$ using $P_B$) and forward samples from the tempered training policy ${\pi }_T$ defined by the forward transition probabilities of the GFlowNet.

4.6 Active Learning with a GFlowNet

An interesting variant on the above scheme is one where the GFlowNet sampler is used not just to produce negative examples for the energy function but also to actively explore the environment. Jain et al. (2022) use an active learning scheme where the GFlowNet is used to sample candidates $x$ for which we expect the reward $R(x)$ to be generally large (since the GFlowNet approximately samples proportionally to $R(x)$). The challenge is that evaluating the true reward $R$ for any $x$ is computationally expensive and can potentially be noisy (for example, a biological assay to measure the binding energy of a drug to a given target protein). Thus, instead of using the true reward directly, the authors introduce a proxy $\widehat{f}$ (which approximates the true reward function $f$), which is used to train the GFlowNet. This would lead to a setup similar to Sec. 4.5, with an inner loop where a GFlowNet is trained to match the proxy $\widehat{f}$, and an outer-loop where the proxy $\widehat{f}$ is learned in a supervised fashion using $(x,y)$ pairs, where $x$ is proposed by the GFlowNet, and $y$ is the corresponding true reward from the environment (for example, outcome of a biological of chemical assay). It is important to note here that the GFlowNet and the proxy are intricately linked since the coverage of proxy $\widehat{f}$ over the domain of $x$ relies on diverse candidates from the GFlowNet. And similarly, since the GFlowNet matches a reward distribution defined by the proxy reward function $\widehat{f}$, it also depends on the quality of the true reward function $f$.

This setup can be further extended by incorporating information about how novel a given candidate is, or how much epistemic uncertainty, $u(x,f)$, there is in the prediction of $\widehat{f}$. We can use the acquisition function heuristics (like Upper Confidence Bound (UCB) or Expected Improvement (EI)) from Bayesian optimization (Močkus, 1975; Srinivas et al., 2010) to combine the predicted usefulness $\widehat{f}(x)$ of configuration $x$ with an estimate of the epistemic uncertainty around that prediction. Using this as the reward can allow the GFlowNet to explore areas where the predicted usefulness is high ($\widehat{f}(x)$ is large) and at the same time explore areas where there is more information to be gathered about useful configurations of $x$. The uncertainty over the predictions of $\widehat{f}$ with the appropriate acquisition function can provide more control over the exploratory behaviour of GFlowNets.

As discussed by Bengio et al. (2021) when comparing GFlowNets with return-maximizing reinforcement learning methods, an interesting property of sufficiently trained GFlowNets is that they will sample from all the modes of the reward function, which is particularly desirable in a setting where exploration is required, as in active learning. The experiments in the paper also demonstrate this advantage experimentally in terms of the diversity of the solutions sampled by the GFlowNet compared with PPO, an RL method that had previously been used for generating molecular graphs and that tends to focus on a single mode of the reward function.

4.7 Estimating Entropies, Conditional Entropies and Mutual Information

Definition 34.

Given a reward function $R$ with $\forall s$, we define the entropic reward function $R^{\prime }$ associated with $R$ as:

$$R^{\prime }(s)=-R(s)\log R(s).$$ (53)

In brief, in this section, we show that we can estimate entropies by training two GFlowNets: one that estimates flows as usual for a target terminal reward function $R(s)$, and one that estimates flows for the corresponding entropic reward function. We show below that we obtain an estimator of entropy by looking up the flow in the initial state, and if we do this exercise with conditional flows, we get conditional entropy. Once we have the conditional entropy, we can also estimate the mutual information.

Proposition 35.

Consider a flow network $(G,F)$ such that the terminating flows match a given reward function $R$, i.e. $\forall s\in {𝒮}^f,F(s\to s_f)=R(s)$, with for all $s$, and a second flow network $(G,F^{\prime })$ with the same pointed DAG, but with a flow function for which the terminating flows match the entropic reward function $R^{\prime }$ (Eq. 53), then the entropy $H[S]$ associated with the terminating state random variable $S\in {𝒮}^f$ with distribution $P_T(S=s)=\frac{R(s)}{Z}$ (Eq. 18) is

$$H[S]:=-\sum _s{P}_T(s)\log P_T(s)=\frac{F^{\prime }(s_0)}{F(s_0)}+\log F(s_0).$$ (54)

Proof First apply the definition of $P_T(s)$, then Eq. 11 on both flows:

$-{\displaystyle \sum _s}{P}_T(s)\log P_T(s)$ $=-{\displaystyle \sum _s}{\displaystyle \frac{R(s)}{F(s_0)}}(\log R(s)-\log F(s_0))$

$={\displaystyle \frac{\left(-{\sum }_{s}R(s)\log R(s)\right)+\left(\log F(s_0){\sum }_{s}R(s)\right)}{F(s_0)}}$

$={\displaystyle \frac{F^{\prime }(s_0)}{F(s_0)}}+\log F(s_0).$

Note that we need to make sure that the rewards $R^{\prime }(s)$ (and thus the flows) are positive.

Proposition 36.

Given a set $𝒳$ of conditioning variables, consider a conditional flow network defined by a conditional flow function $F$, for which the terminating flows match a target reward $R$ family (conditioned on $x\in 𝒳$) that satisfies for all $s$, and a second conditional flow network defined by a conditional flow function $F^{\prime }$, for which the terminating flows match the entropic reward functions $R_x^{\prime }$ (Eq. 53), then the conditional entropy $H[S\mid x]$ of random terminating states $S\in {𝒮}^f$ consistent with condition $x$ is given by

$$H[S\mid x]=\frac{F^{\prime }(s_0\mid x)}{F(s_0\mid x)}+\log F(s_0\mid x).$$ (55)

In particular, for a state-conditional GFlowNet ($𝒳=𝒮$ is the state space of the DAG), we obtain

$$H[S\mid s]=\frac{F^{\prime }(s\mid s)}{F(s\mid s)}+\log F(s\mid s).$$ (56)

More generally, the mutual information MI$(S;X)$ between the random draw of a terminating state $S=s$ according to $P_T(s\mid x)$ and the conditioning random variable $X$ is

$$\mathrm{MI}(S;X)=H[S]-E_X[H[S\mid X]]=\frac{F^{\prime }(s_0)}{F(s_0)}+\log F(s_0)-E_X\left[\frac{F^{\prime }(s_0\mid X)}{F(s_0\mid X)}+\log F(s_0\mid X)\right]$$ (57)

where $F(s)$ and $F^{\prime }(s)$ indicate the unconditional flows (trained with no condition $x$ given) while $F(s\mid x)$ and $F^{\prime }(s\mid x)$ are their conditioned counterparts.

Proof The proof of Eq. 55 follows from the fact that each $(G_x,F_x)$ is a flow network, to which we can apply Prop. 35. Eq. 57 is a direct consequence of the definition of the Mutual Information, Eq. 54 and Eq. 55.

If we have a sampling mechanism for $P(X)$, we can thus approximate the expectation in Eq. 57 by a Monte-Carlo average with draws from $P(X)$.

5 GFlowNets on Sets, Graphs, and to Marginalize Joint Distributions

5.1 Set GFlowNets

We first define an action space for constructing sets and we view the GFlowNet as a means to generate a random set $S$ and to estimate quantities like probabilities, conditional probabilities or marginal probabilities for realizations of this random variable. The elements of those sets are taken from a larger “universe” set $𝒰$.

Definition 37.

Given a “universe” set $𝒰$, consider the pointed DAG $G=(𝒮,𝔸)$, where $𝒮:=2^{𝒰}\cup \{s_f\}$ is the set of all subsets of $𝒰$ with an additional state $s_f$, $s_0=\{\}$ is the empty set, and for any two subsets $s,s^{\prime }$ of $𝒰$, $s\to s^{\prime }\in 𝔸\iff \exists a\in 𝒰\setminus s,s^{\prime }=s\cup \{a\}$; meaning that each transition in the DAG corresponds to adding one element of $𝒰$ to the current subset. Additionally all subsets are connected to $s_f$, i.e. $\forall s\in 𝒮,s\to s_f\in 𝔸$. A set flow network is a flow network on this graph $G$, and a set GFlowNet is an estimator of such a flow network, as defined in Sec. 3. The target terminal reward function $R:s\mapsto F(s\to s_f)$ satisfies:

(58)

A set flow network defines a terminating probability distribution $P_T$ on states (see Def. 13 and Eq. 44), with

$$P_T(s)=e^{-ℰ(s)+ℱ(s_0)}=\frac{F(s\to s_f)}{F(s_0)}$$ (59)

where $ℰ$ represents the energy $-\log R$. Similarly, Cor. 33 provides us with a formula for conditional probabilities of a given superset $s^{\prime }$ of a given set $s$ under $P_T$,

$$P_T(s^{\prime }\mid s^{\prime }\supseteq s):=e^{-ℰ(s^{\prime })+ℱ(s)}=\frac{F(s^{\prime }\to s_f)}{F(s\mid s)}$$ (60)

where $ℱ$ indicates free energy (see Def. 27).

Remember that with a GFlowNet with state and edge flow estimator $\widehat{F}$, it is not guaranteed that $\widehat{F}(s)=R(s)$ for all states $s\in 𝒮\setminus \{s_f\}$, so we could estimate probabilities with

$${\widehat{P}}_T(s)=\frac{\widehat{F}(s\to s_f)}{\widehat{F}(s_0)}$$ (61)

or alternatively

$${\widehat{P}}_T(s)=\frac{R(s)}{\widehat{F}(s_0)}.$$ (62)

Similarly, we can estimate conditional superset probabilities with Eq. 60 or with

$${\widehat{P}}_T(s^{\prime }\mid s^{\prime }\supseteq s)=\frac{R(s^{\prime })}{\widehat{F}(s\mid s)},$$ (63)

none of which are guaranteed to exactly sum to 1.

We can also compute the marginal probability over all supersets of a given set $s$, as shown below.

Proposition 38.

Following the notations of Def. 37, let $𝔖(s)=\{s^{\prime }\supseteq s\}$ be the set of all supersets of a set $s$. The probability of drawing any element from $𝔖(s)$ given a set flow network is

$$P_T(𝔖(s))=\sum _{s^{\prime }\supseteq s}{P}_T(s^{\prime })=\frac{e^{-ℱ(s)}}{Z}=\frac{F(s\mid s)}{F(s_0)}.$$ (64)

Proof We can rewrite the sum as follows, first applying the definition of $P_T$ (Eq. 18), and Prop. 32:

${\displaystyle \sum _{s^{\prime }\supseteq s}}{P}_T(s^{\prime })$ $={\displaystyle \sum _{s^{\prime }\ge s}}{\displaystyle \frac{F(s^{\prime }\to s_f)}{F(s_0)}}$

$={\displaystyle \frac{F(s\mid s)}{F(s_0)}}$

$={\displaystyle \frac{e^{-ℱ(s)}}{Z}}$

where we notice that for states that are sets, the order relationship $s\le s^{\prime }$ is equivalent to the subset relationship $s\subseteq s^{\prime }$.

To summarize, a GFlowNet which is trained to match a given energy function (i.e. derived rewards) over sets can be used to represent that distribution, sample from it, estimate the probability of a set under it, estimate the partition function, search for the lowest energy set, sample a conditional distribution over supersets of a given set, estimate that conditional distribution for a given pair of set and superset, compute the marginal probability of a subset (i.e., summing over the probabilities of the supersets), and compute the entropy of the set distribution or of the conditional distribution of supersets of a set. For example, using a GFlowNet to model distributions over sets has been successfully used in Malik et al. (2023), where the goal is to iteratively select an informative batch of points for efficient active learning. The authors additionally show that the distribution over sets defined by a mutual information-based reward (Sec. 4.7) can be approximated to satisfying levels using a neural network as a function approximator for the GFlowNet policy.

5.2 GFlowNet on Graphs

A graph is a special kind of set in which there are two kinds of elements: nodes and edges, with edges being pairs of node indices. Graphs may also have content attached to nodes and/or edges. The set operations described in the previous section can thus be specialized accordingly. Some actions (i.e. edges in the GFlowNet DAG) could insert a node while other actions could insert an edge. The set of allowable actions can be limited, for example to make sure the graph has a single connected component, or to ensure acyclicity. Like for sets in general, one cannot have an action which adds a node or edge which is already in the set. Deleu et al. (2022) used a GFlowNet over DAGs in order to learn an approximation of the posterior distribution over the graphical structure of a Bayesian Network. While the GFlowNet was learned by minimizing some loss derived from flow-matching conditions as in Sec. 3.2, they showed that the resulting distribution is an accurate approximation of the target posterior. Since graphs are sets, all the GFlowNet operations on sets can be applied on graphs.

5.3 Marginalizing over Missing Variables

The ability of GFlowNets to capture probability distributions over sets can be applied to modeling the joint distribution over random variables, to calculating marginal probabilities over given subsets of variable values, and to sampling or computing probabilities for any conditional (e.g., for a subset of variables given another subset of variables).

Let $X=(X_1,X_2,\mathrm{\dots },X_n)$ be a composite random variable with $n$ element random variables $X_i$, $1\le i\le n$, each with possible values $x_i\in {𝒳}_i$ (not necessarily numbers). If we are given an energy function or a terminal reward function $R(x)$ to score any instance $X=x$, we can train a particular kind of set GFlowNet for which the set elements are pairs $(i,x_i)$ and at most one element in the set with index $i$, for all $i$. The only allowed terminating transitions are when the set has exactly size $n$ and every size-$n$ set $s$ terminates on the next transition.

Note how that GFlowNet can sample an $X$ in any possible order, if $P_B$ allows that order. Given an existing set of $(i,x_i)$ pairs (represented by a state $s={\{(i,x_i)\}}_i$), we can estimate the marginal probability of that subset of variables (implicitly summing over all the missing ones, see Eq. 64). We can sample the other variables by setting the state at $s$ and continuing to sample from the GFlowNet’s policy (the learned forward transition probability function). We can sample from a chosen subset $S^{\prime }$ of the other variables by constraining that policy to only add elements which are in $S^{\prime }$. In addition, we can do all the other things that are feasible on set GFlowNets, such as estimating the partition function, sampling in an order which prefers the early subsequences with the largest marginal probability, searching for the most probable configuration of variables, or estimating the entropy of the distribution.

5.4 Modular Energy Function Decomposition

Let us see how we can apply the graph GFlowNet framework to a special kind of graph: a factor graph (Kschischang et al., 2001) with reusable factors. This will yield a distribution $P_T(g)$ over graphs $g$, each of which is associated with an energy function value $ℰ(g)$ (and associated reward $R(g)$). Energy-based models are convenient because they can decompose a joint probability into independent pieces (possibly corresponding to independent mechanisms, Schölkopf et al., 2012; Goyal et al., 2019; Goyal and Bengio, 2020), each corresponding to a factor of a factor graph. In our case, we would like a shared set of factors $𝔽$ to be reusable across many factor graphs $g$. The factor graph will provide an energy and a probability over a set of random variables $𝒱$. Let the graph $g={\{(F^i,v^i)\}}_i$ be written as a set of pieces $(F^i,v^i)$, where $F^i\in 𝔽$ is the index of a factor with energy function term ${ℰ}_{F^i}$, selected from the pool $𝔽$ of possible factors, and where $v^i=(v_1,v_2,\mathrm{\dots })$ is a list of realizations of the random variables $V_j\leftarrow v_j$, where $V_j\in 𝒱$ is a node of the factor graph. That list defines the edges of the factor graph connecting variable $V_j$ with the $j$-th argument of ${ℰ}_{F^i}$. Let us denote ${ℰ}_{F^i}(v^i)$ the value of this energy function term ${ℰ}_{F^i}$ applied to those values $v^i$, i.e.,

$${ℰ}_{F^i}(v^i)={ℰ}_{F^i}(v_1,v_2,\mathrm{\dots }).$$ (65)

The total energy function of such a graph can then be decomposed as follows

$$ℰ(g)=\sum _i{ℰ}_{F^i}(v^i).$$ (66)

What is interesting with this construction is that the graph GFlowNet can now sample a graph $g$, possibly given some conditioning observations $x$: see Sec. 4.5 on how GFlowNets can be trained jointly with an energy function, including the case where only some random variables are observed. Hence, given some observed variables (not necessarily always the same), the graph GFlowNet can sample a latent factor graph containing and connecting (with energy function terms) both observed and latent random variables, and whose structure defines an energy function over values of the joint observed and latent variables.

Not only can we use the compositional nature of the objects generated by a GFlowNet to decompose the total energy into reusable energy terms corresponding to ideally independent mechanisms, but we can also decompose the GFlowNet itself into modules associated with each mechanism. The action space of this graph GFlowNet is fairly complex, with each action corresponding with the addition of a latent variable $V_k$ or the addition of a graph piece $(F^i,v^i)$. Such an action is taken in the context of the state of the GFlowNet, which is a partially constructed graph (arising from the previous actions). The GFlowNet and its associated energy function parameters are thus decomposed into modules. Each module knows how to compute an energy function ${ℰ}_{F_i}$ and how to score and sample competitively (against the other modules) a new graph piece (to insert the corresponding factor in the graph).

Consider some observed variables (a subset of the $V_j$’s with their values $v_j$), collectively denoted $x$. Consider a graph $g$ among those compatible with $x$ (i.e. with some nodes corresponding to $V_j=v_j$ for the observed variables) and denote $h$ the specification of the part of $g$ not already provided by $x$. We can think about latent variable $h$ as the explanation for the observed $x$. Note how marginalizing over all the possible $h$, we can compute the free energy of $x$. The principles of Sec. 4.5 can be applied to train such an energy-based GFlowNet. It also makes sense to represent a prior over graph structures in the energy function. For example, we may prefer sparse factors (with few arguments), and we may introduce soft or hard constraints having to do with a notion of type that is commonly used in computer programming and in natural language. Each random variable in the graph can have as one of its attributes a type, and each factor energy function argument can expect a type. Energy function terms can be added to construct this prior by favouring graph pieces $(F_i,v^i={(v_k)}_k)$ in which the type of variable $V_k$ (of which $v_k$ is a realization) matches the type expected of the $k$-th argument of $F_i$. This is very similar to attention mechanisms (Bahdanau et al., 2014; Vaswani et al., 2017), which can be seen to match a query (an expected type) with a key (the type associated with an element). For instance, Pan et al. (2023) propose the Forward-Looking GFlowNet framework that exploits the modularity of the energy function, as in Eq. 66 and obtain a better approximation of the target distribution than those resulting from regular losses.

6 Continuous or Hybrid Actions and States

All of the mathematical developments above have used sums over states or actions, with the idea that these would be elements of a discrete space. However, for the most part one can replace these sums by integrals in case the states or actions are either continuous or hybrid (with some discrete components and some continuous components). Beyond this, we discuss below what the presence of continuous-valued actions and states changes to the GFlowNet framework.

Although there are explicit sums respectively over successors and predecessors which come up in Eq. 40, such sums are also hiding in the detailed balance constraint of Eq. 26. Indeed, these sums are implicit as part of the normalizing constant in the conditional density of the next state or previous state in $P_F(s_{t+1}\mid s_t)$ and $P_B(s_t\mid s_{t+1})$. We consider below ideas to deal with this challenge.

6.1 Integrable Normalization Constants

We first note that if we can handle a continuous state, we can also handle a hybrid state, as follows. Let the state be decomposed into

$$s=(s^i,s^x)$$ (67)

where $s^i$ is discrete and $s^x$ is continuous. Then we can decompose any of the transition conditionals as follows:

$$P_F(s_{t+1}\mid s_t)=P(s_{t+1}^x\mid s_{t+1}^i,s_t)P(s_{t+1}^i\mid s_t).$$ (68)

We note that this is formally equivalent to decomposing the transition into two transitions, first to perform the discrete choice into the next state, and second to perform the continuous choice into the next state (given the discrete choice). Having continuous-valued inputs to a neural net is no problem. The challenge is to represent continuous densities on the output, with the need to both being able to compute the density of a particular value (say $P(s_{t+1}^x\mid s_{t+1}^i,s_t)$) and to be able to sample from it. Computing categorical probabilities and sampling from a conditional categorical is standard fare, so we only discuss the continuous conditional. One possibility is to parametrize $s_{t+1}^x\mid s_{t+1}^i,s_t$ with a density for which the normalization constant is a known tractable integral, like the Gaussian. However, that may limit capacity too much, and may prevent a good minimization of the detailed balance or flow-matching loss. One workaround is to augment the discrete part of the state $s_i$ with extra dimensions corresponding to “cluster IDs”, i.e., partition the continuous density into a mixture. We know that with enough mixture components, we can arbitrarily well approximate densities from a very large family. Other approaches include modeling the conditional density with an autogressive or normalizing flow model (Rezende and Mohamed, 2015, with a different meaning of the word flow), or, like in denoising diffusion models (Sohl-Dickstein et al., 2015), decomposing sampling of $s^x$ into several resampling steps, transforming its ditribution from a simple one to complex one.

To guarantee that the detailed balance constraint can be exactly satisfied, we could go further and think about parametrizing the edge flow $F((s_i,s_x)\to (s_i^{\prime },s_x^{\prime }))$, and note that this is the natural parametrization if we use the node-based flow-matching loss. For example, keeping with the Gaussian example, we would now have a joint Gaussian energy in the vector $(s_x,s_x^{\prime })$ for each feasible discrete component indexed by $(s_i,s_i^{\prime })$. Note that in practical applications of GFlowNets, not all the transitions that satisfy the order relationship are generally allowed in the GFlowNet’s underlying DAG. For example, with set GFlowNets, the only allowed actions add one element to the set (not an arbitrary number of elements). These constraints on the action space mean that the number of legal $(s_i,s_i^{\prime })$ pairs is manageable and correspond to the number of discrete actions. The overall action is therefore seen as having a discrete part (choosing $s_i^{\prime }$ given $s_i$) and a continuous part (choosing $s_x^{\prime }$ given $s_i^{\prime }$, $s_i$ and $s_x$). With such a joint flow formulation, the forward and backward conditional densities can be computed exactly and be compatible with each other.

6.2 GFlowNets in GFlowNets

Another way to implement an edge flow involving continuous variables is to use a lower-level $⟨$GFlowNet, energy function$⟩$ pair to represent its flow, conditional probabilities and sample from them. Remember that such a pair can be trained following the approach discussed in Sec. 4.5. Instead of a joint Gaussian for $(s_x,s_x^{\prime })$ given $(s_i,s_i^{\prime })$ we could have a smaller-scale GFlowNet and energy function (representing an edge flow in the outer GFlowNet) to handle a whole family of transitions of a particular type in a larger-scale outer GFlowNet. Imagine that we have a fairly arbitrary energy function for such a transition, with parameters that we will learn. Then we can also train a GFlowNet to sample in either direction (either from $s_t$ to $s_{t+1}$ or from $s_{t+1}$ to $s_t$) and to evaluate the corresponding normalizing constants (and hence, the corresponding conditional probabilities). The discrete aspect of the state and of its transition may correspond to a family of transitions (e.g., insert a particular type of node in a graph GFlowNet), and a separate $⟨$GFlowNet, energy function$⟩$ module may be specialized and trained to handle such transitions.

While this paper was under review, Lahlou et al. (2023) extended the theory of GFlowNets to more general state spaces, including continuous ones, and experimentally validated that the usual losses can be used to effectively perform inference in continuous domains.

7 Related Work

There are several classes of related literature that concern the problem of generating a diversity of samples, given some energy or reward signal, in particular:

• generative models (in particular deep learning ones),

• RL methods that maximize reward with some form of exploratory behavior or smoothness prior,

• MCMC methods that solve the problem of sampling from $p(x)\propto f(x)$ in principle,

• evolutionary methods, that can leverage group diversity through iterations over a population of solutions.

In what follows we discuss these and offer insights into similarities and differences between GFlowNets and these approaches. Note that the literature related to this problem is much larger than we can reference here, and extends to many other subfields of ML, such as GANs (Kumar et al., 2019), VAEs (Kingma and Welling, 2013; Kusner et al., 2017), and normalizing flows (Dinh et al., 2014, 2016; Rezende and Mohamed, 2015). Yet another related type of approach are the Bayesian optimization methods (Močkus, 1975; Srinivas et al., 2010), which have also been used for searching in the space of molecules (Griffiths and Hernández-Lobato, 2017). The main relation with Bayesian optimization methods is that GFlowNets are generative and can thus complement Bayesian optimization methods which scan a tractable list of candidates. When the search space is too large to be able to separately compute a Bayesian optimization acquisition function score on every candidate, using a generative model is appealing. In addition, GFlowNets are used to explore the modes of the distribution rather than to search for the single most dominant mode. This difference is similar to that with classical RL methods, discussed further below.

7.1 Contrast with Generative Models

The main difference between GFlowNets and established deep generative models like VAEs or GANs is that whereas the latter are trained by being provided a finite set of examples sampled from the distribution of interest, a GFlowNet is normally trained by being provided an energy function or a reward function.

This reward function tells us not just about the samples that are likely under the distribution of interest (which we can think of as positive examples) but also about those that are unlikely (which we can think of as negative examples) and also about those in-between (whose reward is not large but is not zero either). If we think of the maximum likelihood training objective in those terms, it is like a reward function that gives a high reward to every training example (seen as a positive example, where the probability should be high) and a zero reward everywhere else. However, other reward functions are possible, as seen in the application of GFlowNets to the discovery of new molecules (Bengio et al., 2021), where the reward is not binary and increases monotonically as a function of the value of a desirable property of the candidate molecule.

Note however that the difference with other generative modeling approaches blurs when we include the learning of the energy function along with the learning of the GFlowNet sampler, as outlined in Sec. 4.5. In that case, the pair comprising the trainable GFlowNet sampler and the trainable energy function achieves a similar objective as a trainable generative model. Note that GFlowNets have been designed for generating discrete variable-size compositional structures (like sets or graphs), for both latent and observed variables, whereas GANs, VAEs or normalizing flows start from the point of view of modeling real-valued fixed-size vectors using real-valued fixed-size latent variables.

An interesting difference between GFlowNets and most generative model training frameworks (typically some variation on maximum likelihood) is in the very nature of the training objective for GFlowNets, which came about in the context of active learning scenarios. Whereas the GFlowNet training pairs $(s,R(s))$ can come from any distribution over $s$ (any full-support training policy ${\pi }_T$), which does not have to be stationary (and indeed will generally not be, in an active learning setting), the maximum likelihood framework is very sensitive to changes in the distribution of the data it sees. This is connected to the “offline learning” property of the flow matching objective (Sec. 3.3.3, among others).

7.2 Contrast with Regularized Reinforcement Learning

The flow-matching loss of GFlowNets (Bengio et al., 2021) arose from the inspiration of the temporal-difference training (Sutton and Barto, 2018) objectives associated with the Bellman equation. The flow-matching equations are analogous to the Bellman equation in the sense that the training objective is local (in time and states), credit assignment propagates through a bootstrap process and tries to fix the parametrization so that these equations are satisfied, knowing that if they were (everywhere), we would obtain the desired properties. However, these desired properties are different, as elaborated in the next paragraph. The context in which GFlowNets were developed is also different from the typical way of thinking about agents learning in some environment: we can think of the deterministic environments of GFlowNets as involving internal actions typically needed by a cognitive agent that needs to perform some kind of inference through a sequence of steps (predict or sample some things given other things), i.e., through actions internal to the agent and controlling its computation. This is in contrast with the origins of RL, focused on the actions of an agent in an external and unknown stochastic environment. GFlowNets were introduced as a tool for learning an internal policy, similar to the use of attention in modern deep learning, where we know the effect of actions, and the composition of these actions defines an inference machinery for that agent.

Classical RL (Sutton and Barto, 2018) control methods work by maximizing return in Markov Decision Processes (MDPs); their focus is on finding the policy ${\pi }^{\ast }\in {\text{argmax}}_{\pi }{V}^{\pi }(s)\forall s$ maximizing the expected return $V^{\pi }(s)$, which happens to provably be achieved with a deterministic policy (Sutton and Barto, 2018), even in stochastic MDPs. In a deterministic MDP, of interest here, this means that training an RL agent is a search for the most rewarding trajectory, or in the case of terminal-reward-only MDPs (again of interest here), the most rewarding terminating state.

Another perspective, that emerged out of both the probabilistic inference literature (Toussaint and Storkey, 2006) and the bandits literature (Auer et al., 2002), is concerned with finding policies of the form $\pi (a|s)\propto f(s,a)$. It turns out that maximizing both return and entropy (Ziebart et al., 2008) of policies in a control setting yield policies such that

$$p(\tau )=\left[p(s_0)\prod _{t=0}^{T-1}P(s_{t+1}|s_t,a_t)\right]\exp \left(\eta \sum _{t=0}^{T-1}R(s_t,a_t)\right)$$ (69)

where $\tau =(s_0,a_0,s_1,a_1,\mathrm{\dots },s_T)$ and $\eta$ can be seen as a temperature parameter. This result can also be found under the control-as-inference framework (Haarnoja et al., 2017; Levine, 2018). In deterministic MDPs with terminal rewards and no discounting of future rewards, this simplifies to $p(\tau )\propto \exp (\eta \rho (\tau ))$, where $\rho$ is the return.

In recent literature, this entropy maximization (MaxEnt) is often interpreted as a regularization scheme (Nachum et al., 2017), entropy being used either as an intrinsic reward signal or as an explicit regularization objective to be maximized. Another way to understand this scheme is to imagine ourselves in an adversarial bandit setting (Auer et al., 2002) where each arm corresponds to a unique trajectory, drawn with probability $\propto \exp (\rho (\tau ))$.

An important distinction to make between MaxEnt RL and GFlowNets is that, in the general case they do not find the same result. A GFlowNet learns a policy such that $P_T(s)\propto R(s)$, whereas MaxEnt RL (with appropriately chosen temperature and $R$) learns a policy such that $P_T(s)\propto n(s)R(s)$, where $n(s)$ is the number of paths in the DAG of all trajectories that lead to $s$ (a proof is provided in Bengio et al., 2021). An equivalence only exists if the DAG minus $s_f$ is a tree rooted at $s_0$, which has been found to be useful (Buesing et al., 2019). What this overweighting by a factor $n(s)$ means practically is that states corresponding to longer sequences (which typically will have exponentially more paths to them) will tend to be sampled much more often (typically exponentially more often) than states corresponding to shorter sequences. Clearly, this breaks the objective of sampling terminating states in proportion to their reward and provides a strong motivation for considering GFlowNets instead.

Another perspective on maximizing entropy in RL is that one can also maximize entropy on the states’ stationary distribution $d^{\pi }$ (Ziebart et al., 2008), rather than the policy. In fact, one can show that the objective of training a policy such that $P_T(s)\propto R(s)$ is equivalent to training a policy that maximizes $r(s,a)=\log R(s,a)-\log d^{\pi }(s,a)$. Unfortunately, computing stationary distributions, although possible (Nachum et al., 2019; Wen et al., 2020), is not always tractable nor precise enough for purposes of reward regularization.

7.3 Contrast with Monte-Carlo Markov Chain methods

MCMC has a long and rich history (Metropolis et al., 1953; Hastings, 1970; Andrieu et al., 2003), and is particularly relevant to the present work, since it is also a principled class of methods towards sampling from $P_T(s)\propto R(s)$. MCMC-based methods have already found some amount of success with learned deep neural networks used to drive sampling (Grathwohl et al., 2021; Dai et al., 2020; Xie et al., 2021; Nash and Durkan, 2019; Seff et al., 2019).

An important drawback of MCMC is its reliance on iterative sampling (forming the Markov chain, one configuration at a time, each of which is like a terminating state of a GFlowNet): a new state configuration is obtained at each step of the chain by making a small stochastic change to the configuration in the previous step. Although these methods guarantee that asymptotically (in the length of the chain) we obtain samples drawn from the correct distribution, there is an important set of distributions for which finite chains are unlikely to provide enough diversity of the modes of the distribution.

This is known as the mode-mixing problem (Jasra et al., 2005; Bengio et al., 2013; Pompe et al., 2020): the chances of going from one mode to a neighboring one may become exponentially small (and thus require exponentially long chains) if the modes are separated by a long sequence of low-probability configurations. This can be alleviated by burning more computation (sampling longer chains) but becomes exponentially unsustainable with increased mode separation. The issue can also be reduced by introducing random sampling (e.g., drawing multiple chains) and simulated annealing (Andrieu et al., 2003) to facilitate jumping between modes. However, this becomes less effective in high dimensions and when the modes occupy a tiny volume (which can become an exponentially small fraction of the total space as its dimension increases) since random sampling is unlikely to land in the neighborhood of a mode.

In contrast, GFlowNets belong to the family of amortized sampling methods (which includes VAEs, Kingma and Welling, 2013), where we train a machine learning system to produce samples: we have exchanged the complexity of sampling through long chains for the complexity of training the sampler. The potential advantage of such amortized samplers is when the distribution of interest has generalizable structure: when it is possible to guess reasonably well where high-probability samples can be found, based on the knowledge of a set of known high-probability samples (the training set). This is what makes deep generative models work in the first place and thus suggests that in such high-dimensional settings where modes occupy tiny volumes (as per the manifold hypothesis, Cayton, 2005; Narayanan and Mitter, 2010; Rifai et al., 2011), one can capitalize on the already observed $(x,R(x))$ pairs (where $x$ is an already visited configuration and $R(x)$ its reward) to “jump” from known modes to yet unvisited ones, even if these are far from the ones already visited.

How well this will work then depends on the ability to generalize of the learner, i.e., on the strength and appropriateness of its inductive biases, as usual in machine learning. In the case where there is no structure at all (and thus no possibility to generalize when learning about the distribution), there is no reason to expect that amortized ML methods will fare better than MCMC. But if there is structure, then the exponential cost of mixing between modes could go away. There is plenty of evidence that ML methods can do a good job in such high-dimensional spaces (like the space of natural images) and this suggests that GFlowNets and other amortized sampling methods would be worth considering where ML generally works well. Molecular graph generation experiments (Bengio et al., 2021) comparing GFlowNets and MCMC methods appear to confirm this.

Another factor to consider (independent of the mode mixing issue) is the amortization of the computational costs: GFlowNets pay a large price upfront to train the network and then a small price (sampling once from $P_T$) to generate each new sample. Instead, MCMC has no upfront cost but pays a lot for each independent sample. Hence, if we want to only sample once, MCMC may be more efficient, whereas if we want to generate a lot of samples, amortized methods may be advantageous. One can imagine settings where GFlowNets and MCMC could be combined to achieve some of the advantages of both approaches.

Evolutionary Methods

Evolutionary methods work similarly to MCMC methods, via an iterative process of stochastic local search, and populations of candidates are found that maximize one or many objectives (Brown et al., 2004; Salimans et al., 2017; Jensen, 2019; Swersky et al., 2020). From such a perspective, they have similar advantages and disadvantages. One practical advantage of these methods is that natural diversity is easily obtainable via group metrics and subpopulation selection (Mouret and Doncieux, 2012). This is not something that is explicitly tackled by GFlowNet, which instead relies on i.i.d. sampling and giving non-zero probability to suboptimal samples as a diversity mechanism.

Sequential Monte-Carlo

Sequential Monte Carlo (SMC, Naesseth et al., 2019; Arulampalam et al., 2002) methods are a class of methods aimed at solving inference problems. Similar to GFlowNets, SMC samplers are trained to sample from a distribution given by its unnormalized probability density, and require forward and backward kernels, as in GFlowNets. Unlike GFlowNets though, they require specifying intermediate targets ${\gamma }_t(z_t)$. In addition, with GFlowNets, the reward normally only comes at the end of the sequence, unlike the per-time-step likelihoods used to reweight particles in SMC. GFlowNets can be applied in settings which do not fit the typical particle filter setting, such as those whose intermediate states do not correspond to valid elements of the sample space. The trajectories do not necessarily represent a sequence of latent variables associated with corresponding observations, as in filtering tasks. The GFlowNet trajectory distributions are only defined by the terminal state reward function.

8 Conclusions and Open Questions

This paper extends and deepens the mathematical framework and mathematical properties of GFlowNets (Bengio et al., 2021). It connects the notion of flow in GFlowNets with that of measure over trajectories and introduces a novel training objective (the detailed balance loss) which makes it possible to choose a parametrization separating the backward policy $P_B$ which controls preferences over the order in which things are done from the constraints imposed by the target reward function.

An important contribution of this paper is the mathematical framework for marginalization or free energy estimation using GFlowNets. It relies on the simple idea of conditioning the GFlowNet so as to push the ability to estimate a partition function already introduced by Bengio et al. (2021) to a much more general setting. This makes it possible in principle to estimate intractable sums of rewards over the terminating states reachable by an arbitrary state, opening the door to marginalization over supergraphs of graphs, supersets of sets, and supersets of (variable,value) pairs. In turn, this provides formulae for estimating entropies, conditional entropies and mutual information.

Appendix A introduces alternatives to the flow-matching objective which may bypass the slow “bootstrapping” propagation of credit information from the end of the action sequence to its beginnings, with so-called direct credit assignment (which has some similarity to policy gradient). This question also motivated the trajectory balance objective introduced by Malkin et al. (2022) and further work in this direction is warranted.

In an attempt to better discern the links between GFlowNets and the more common forms of RL, Appendix C considers the cases where rewards are provided not just at the end but possibly after every action, and where the environment may be stochastic. It also shows how a greedy policy that maximises returns could be obtained from a trained GFlowNet. Another link with RL is introduced in Appendix D by considering something similar to a value function in GFlowNets, i.e., the expected downstream reward from a given state $s$. In a similar spirit, Appendix E generalizes GFlowNets to the case where rewards can be earned along the trajectory, thus introducing a notion of return, but keeping the objective of sampling in proportion to the return.

Appendix F generalizes GFlowNets in another direction, considering the possibility of having a family of flows being learned simultaneously over the same graph. This opens the door to a distributional generalization of GFlowNets (by analogy with distributional RL) as well a notion of unsupervised GFlowNet (where the reward function can be defined after the GFlowNet has been trained in an unsupervised fashion) and learning to sample from the Pareto front defined by a set of reward functions.

Many open questions obviously remain, from the extension to continuous actions and states to hierarchical versions of GFlowNets with abstract actions and integrating the energy function in the GFlowNet parametrization itself, enabling an interesting form of modularization and knowledge decomposition. Importantly, many of the mathematical formulations presented in this paper will require empirical validation to ascertain their usefulness, improve these ideas, turn them into impactful algorithms and explore a potentially very broad range of interesting applications, from replacing MCMC or being combined with MCMC in some settings, to probabilistic reasoning to further applications in active learning for scientific discovery.

Appendix A Direct Credit Assignment in GFlowNets

Similarly to temporal-difference methods, which are based on minimizing the mismatch with respect to the Bellman equations, the flow-matching and detailed-balance losses will take many updates (and sampling many trajectories) to propagate a reward mismatch on a terminating state into the transition probabilities inside the flow network. This would be particularly acute for longer trajectories and prompts the question of alternative more direct training objectives. This section attempts to answer the following question: “Given a training trajectory $\tau$, are there more direct ways of assigning credit to the earlier transitions in the trajectory?”. Malkin et al. (2022) provide an alternative answer by introducing the Trajectory Balance loss (Ex. 6).

We can view the process of sampling a trajectory with a GFlowNet as analogous to sampling a sequence of states in a stochastic recurrent neural network. What makes things complicated is that such a neural network (i) does not directly output a prediction to be matched by some target, and (ii) that the state may be discrete (or a combination of discrete and continuous components).

Regarding (i), we recall that the flow in a transition $s\to s^{\prime }$ (involved in the training objective) is an intractable sum over all the possible trajectories leading to $s$. However, we may be able to obtain a stochastic gradient that, in average over such trajectories, estimates the desired quantity. For this, we exploit the properties of flows to obtain a stochastic gradient estimator, derived through the next three propositions.

Notation for derivatives through the flow-matching constraint: below we sometimes need to distinguish total derivatives (noted $\frac{dy}{dx}$) that take into account the indirect effects due to the flow matching constraint from other derivatives (noted $\frac{\partial y}{\partial x}$) and capturing only direct gradients. Either notation can be used when the constraint does not change the result. Computing indirect gradients through implicit dependencies is an active area of research and commonly utilizes the Implicit Function Theorem, e.g., for implicit layers in fixed-point iteration layers and Deep Equilibrium Models (Bai et al., 2019).

Proposition 39.

Consider the effect of a slight change in the log of the flow at under the flow-matching constraint: it yields a change in the log of the flow at $s^{\prime }$, following the conditional probability $P(s|s^{\prime })$:

$$\frac{d\log F(s^{\prime })}{d\log F(s)}=P(s|s^{\prime })$$ (70)

where $P$ is the distribution on events over trajectories.

Proof We are going to consider a partition of the complete trajectories going through $s^{\prime }$ into those that also go through $s$ (set ${𝒯}_{s^{\prime }}^s$) and those that don’t (set ${𝒯}_{s^{\prime }}^{-s}$). By definition we can write:

$P(s^{\prime })=P({𝒯}_{s^{\prime }}^s)+P({𝒯}_{s^{\prime }}^{-s})$

Because ${𝒯}_{s^{\prime }}^s=\{\tau \in 𝒯:s\in \tau ,s^{\prime }\in \tau \}$, we can write $P({𝒯}_{s^{\prime }}^s)=P(s)P(s^{\prime }|s)$. Hence:

$F(s^{\prime })=F(s)P(s^{\prime }|s)+F({𝒯}_{s^{\prime }}^{-s}).$

Additionally, because the flow in $s$ does not influence trajectories in ${𝒯}_{s^{\prime }}^{-s}$, then $\frac{dF({𝒯}_{s^{\prime }}^{-s})}{dF(s)}=0$; which leads to:

${\displaystyle \frac{dF(s^{\prime })}{dF(s)}}$ $=P(s^{\prime }|s)$

${\displaystyle \frac{dF(s^{\prime })}{d\log F(s)}}$ $=F(s)P(s^{\prime }|s)$

$=F(s\cap s^{\prime })$

$=P(s|s^{\prime })F(s^{\prime })$

${\displaystyle \frac{d\log F(s^{\prime })}{d\log F(s)}}$ $=P(s|s^{\prime }).$

Proposition 40.

Consider the effect of a slight change in the log of the flow on the edge $s\to s^{\prime }$: we obtain a change in the log of the flow at $s^{\prime }$ following the backward conditional probability $P_B(s|s^{\prime })$:

$$\frac{d\log F(s^{\prime })}{d\log F(s\to s^{\prime })}=P_B(s|s^{\prime })$$ (71)

where $P$ is the distribution on events over trajectories.

Proof We first use the chain rule and properties of the derivatives of the log, and then we use from the flow matching constraint that $F(s^{\prime })={\sum }_{s\in Par(s^{\prime })}F(s\to s^{\prime })\Rightarrow \frac{dF(s^{\prime })}{dF(s\to s^{\prime })}=1$:

${\displaystyle \frac{d\log F(s^{\prime })}{d\log F(s\to s^{\prime })}}$ $={\displaystyle \frac{1}{F(s^{\prime })}}{\displaystyle \frac{dF(s^{\prime })}{dF(s\to s^{\prime })}}{\displaystyle \frac{dF(s\to s^{\prime })}{d\log F(s\to s^{\prime })}}$

$={\displaystyle \frac{F(s\to s^{\prime })}{F(s^{\prime })}}{\displaystyle \frac{dF(s^{\prime })}{dF(s\to s^{\prime })}}$

$={\displaystyle \frac{F(s\to s^{\prime })}{F(s^{\prime })}}$

$=P_B(s|s^{\prime }).$

Intuitively, Eq. 70 and Eq. 71 tell us how a perturbation in the flow at one place would result in a change elsewhere in order to maintain the flow match everywhere, using the flow-matching conditions to propagate infinitesimal changes in flow backwards and forwards. Hence, they only become true as we approach the limit of matched flows, and in practice (with an imperfectly trained GFlowNet) the corresponding expressions will be biased. However, we can exploit them to estimate long-range equilibrium gradients and obtain an estimator of credit assignment across a long trajectory in Prop. 41 below.

In the GFlowNet setting, suppose we parametrize the edge flow estimator $F_{\theta }(s\to s^{\prime })$ via parameters $\theta$. In order to understand the effect of a change in $\theta$ on our loss function $ℒ$, we must compute the total derivative $\frac{dL}{d\theta }$ by summing the direct and indirect gradients. In our context, the direct gradient $\frac{\partial L}{\partial \theta }$ is due to the explicit change in loss from changing $\theta$ (not taking the flow-matching constraint into account) and the indirect gradient includes the induced changes in the flow due to the constraint. With this, we are in a position to formalize unbiased estimators for the total derivative $\frac{dL}{d\theta }$ in Prop. 41:

Proposition 41.

Let $ℒ$ be a flow-matching loss computed along a trajectory $\tau =(s_0,s_1,\mathrm{\dots },s_f)$ sampled according to the GFlowNet trajectory distribution $P(\tau )$. Let $ℒ={\sum }_{i}L(s_i)$ decompose the total loss into per-state losses $L(s_i)$ along the trajectory. Let $\theta$ parametrize the edge flow estimator $F_{\theta }(s\to s^{\prime })$. Then, in the limit of the flows becoming matched,

$$G_1≔\sum _i\frac{\partial L(s_i)}{\partial \theta }+\frac{\partial L(s_i)}{\partial \log F(s_i)}\sum _{t=1}^{i-1}\frac{\partial \log F_{\theta }(s_{t-1}\to s_t)}{\partial \theta }$$ (72)

is an unbiased estimator of the total derivative $\frac{dℒ}{d\theta }$, as is

$G_2≔{\displaystyle \sum _i}{\displaystyle \frac{\partial L(s_i)}{\partial \theta }}+{\displaystyle \frac{\partial L(s_i)}{\partial \log F(s_i)}}{\displaystyle \sum _{t=1}^{i-1}}{\displaystyle \sum _{s\in Par(s_t)}}{P}_B(s|s_t){\displaystyle \frac{\partial \log F_{\theta }(s\to s_t)}{\partial \theta }}$ (73)

as is the convex combination

$G=\lambda G_1+(1-\lambda )G_2$ (74)

for any $0\le \lambda \le 1$. Intuitively, $\frac{\partial L(s_i)}{\partial \theta }$ indicates the gradient directly through the occurrences of $F_{\theta }$ in $L(s_i)$, and $F(s_i)$ indicates either the forward or backward flow sum present in $L(s_i)$ to obtain the flow through $s_i$.

Proof We will use the partial derivative notation when not considering the indirect influence due to the matching flow constraint and the total derivative notation when considering it, to apply Prop. 39. Consider the expected value under the GFlowNet’s trajectory distribution of the $G$ in Eq. 72, use the previous proposition (Eq. 70) and the chain rule:

$E[G_1]$ $=E\left[{\displaystyle \sum _i}{\displaystyle \frac{\partial L(s_i)}{\partial \theta }}+{\displaystyle \sum _{t=1}^{i-1}}{\displaystyle \frac{\partial L(s_i)}{\partial \log F(s_i)}}{\displaystyle \frac{\partial \log F_{\theta }(s_{t-1}\to s_t)}{\partial \theta }}\right]$

$={\displaystyle \sum _{s_i}}P(s_i)({\displaystyle \frac{\partial L(s_i)}{\partial \theta }}$

$={\displaystyle \sum _{s_i}}P(s_i)({\displaystyle \frac{\partial L(s_i)}{\partial \theta }}$

$=E\left[{\displaystyle \frac{dℒ}{d\theta }}\right]$

$={\displaystyle \frac{dE[ℒ]}{d\theta }}$

The above demonstration shows that $G$ in Eq. 74 is asymptotically (as the flows become matched) unbiased when $\lambda =1$ because we recover the $G_1$ of Eq. 72. The same proof technique can then be used for $G_2$ which uses transitions $s\to s_t$ sampled from $P_B(s|s_t)$ instead of the trajectory transitions $s_{t-1}\to s_t$ and we obtain that the estimator $G$ in Eq. 74 is asymptotically unbiased when $\lambda =0$:

$E[G_2]$ $=E\left[{\displaystyle \sum _i}{\displaystyle \sum _{t=1}^{i-1}}{\displaystyle \frac{\partial L(s_i)}{\partial \theta }}+{\displaystyle \frac{\partial L(s_i)}{\partial \log F(s_i)}}{\displaystyle \sum _{s\in Par(s^{\prime })}}{P}_B(s|s^{\prime }){\displaystyle \frac{\partial \log F_{\theta }(s\to s_t)}{\partial \theta }}\right]$

${\displaystyle \sum _{s\in Par(s^{\prime })}}{P}_B(s|s^{\prime }){\displaystyle \frac{\partial L(s_i)}{\partial \log F(s_i)}}{\displaystyle \frac{\partial \log F_{\theta }(s\to s^{\prime })}{\partial \theta }})$

$={\displaystyle \sum _{s_i}}P(s_i)({\displaystyle \frac{\partial L(s_i)}{\partial \theta }}+$

$={\displaystyle \sum _{s_i}}P(s_i)({\displaystyle \frac{\partial L(s_i)}{\partial \theta }}+$

$={\displaystyle \frac{dE[ℒ]}{d\theta }}$

where the last identity follows the fourth line in the proof for $G_1$. Finally a convex combination of two unbiased estimators is unbiased, so we obtain that $G$ in Eq. 74 is asymptotically unbiased for any $0\le \lambda \le 1$.

This surprising result says that something very close to policy gradient actually provides an asymptotically (i.e., when flows are matched) unbiased gradient on the parameters of the edge flow, in expectation. Note that it only works exactly in an online setting, i.e., when the trajectory is sampled according to the learner’s current policy. Otherwise, the gradient estimator may be biased (it would be biased anyways in practice because the flows are never perfectly matched). However, if instead of sampling trajectories $\tau$ from the GFlowNet transition probabilities $P_F(s_{t+1}|s_t)$ we sample them from a training distribution $\tilde{P}$ with transition probabilities $\tilde{P}(s_{t+1}|s_t)$ we can calculate the importance weights (by the ratio $P(\tau )/\tilde{P}(\tau )$) and correct the estimator accordingly. Since the training distribution $\tilde{P}$ should be broader and have a full support, the importance ratio cannot explode but there could still be the usual numerical problems with the variance of such importance-weighted estimators.

We now consider the setting in which the sampling policy is only slightly different from $\widehat{P}$ , which is typically the case because we want the sampling policy to be broader and more exploratory, and because we may be using delayed data, e.g., with a replay buffer. This slight difference may induce a bias but it might still be advantageous to use the above gradient estimator. Note how it does not come in conflict with the gradient of the flow matching loss (which is the first term in $G$). The expected advantage of using $G$ is that it may initially speed up training by directly providing updates to earlier transitions of a complete trajectory. However, analogous to the trade-off between temporal-difference methods and policy-gradient methods, this may come at the price of higher variance.

This estimator is unbiased when the flows are matched and when the trajectory is sampled according to the GFlowNet’s distribution, but it also makes a lot of intuitive sense: if the estimated flow at $s_i$ is too small (in the eye of $L_i$) one can clearly push that flow up by increasing the probability of a transition on a path leading to $s_i$. Even if we consider a slightly different trajectory sampling distribution, so long as it leads to $s_i$ we would expect that increasing its probability would increase the probability of ending up in $s_i$ (see Prop. 39).

If the state has a continuous component, we could also increase the probability of ending up in $s_i$ by choosing more often a more probable path to $s_i$. This could be calculated by backpropagating through the state transitions (with some form of backpropagation through time). However, if the transitions are not fully known or are not differentiable, this approach may be more challenging, and is related to similar questions raised with credit assignment in reinforcement learning with continuous states.

Finally, keep in mind that the more direct credit assignment terms in $G$ have to be combined with the local terms $\frac{\partial L_i}{\partial \theta }$ which make sure that the flows becomes better matched, since flow-matching is a necessary condition for $G$ to be unbiased.

Appendix B Conditional GFlowNets

Definition 42.

Consider a set of conditioning information $𝒳$, a family of DAGs $𝒢=\{G_x=({𝒮}_x,{𝔸}_x),x\in 𝒳\}$, a family of target reward functions $ℛ=\{R_x:{𝒮}_x^f\to {ℝ}^{+},x\in 𝒳\}$, and a flow parametrization $({𝒪}_x,{\mathrm{\Pi }}_x,{\mathscr{H}}_x)$ of $(G_x,R_x)$ for every $x\in 𝒳$. The three functions $𝒪:x\in 𝒳\mapsto {𝒪}_x$, $\mathrm{\Pi }:x\in 𝒳\mapsto {\mathrm{\Pi }}_x$, and $\mathscr{H}:x\in 𝒳\mapsto {\mathscr{H}}_x$ parametrize the family of DAGs and target reward functions, and we say that $(𝒪,\mathrm{\Pi },\mathscr{H})$ form a conditional flow parametrization of $(𝒳,𝒢,ℛ)$. The tuple $(𝒳,𝒢,ℛ,𝒪,\mathrm{\Pi },\mathscr{H})$ is called a conditional GFlowNet.

For clarity, for an object $o\in 𝒪$, we will use $o_x$ and $o(x)$ interchangeably.

As with the unconditional case, conditional GFlowNets provide a way to sample from different target reward functions $R_x$ simultaneously. Given a configuration $o\in 𝒪$ of the conditional GFlowNet and a condition $x\in 𝒳$, ${\pi }_x:={\mathrm{\Pi }}_x(o_x)$ is a distribution over ${𝒯}_x$, the set of complete trajectories in $G_x$, which implicitly defines a terminating state probability measure in $G_x$:

$$\forall x\in 𝒳\forall s\in {𝒮}_x^f{P}_T(s\mid x):=\sum _{\tau \in {𝒯}_x:s\to s_f\mid x\in \tau }{\pi }_x(\tau ),$$ (75)

where the dependence on $o$ in $P_T$ is omitted for clarity.

Conditional GFlowNets cast the problem of sampling from target reward functions to a search problem: searching for objects $o\in 𝒪$ such that $o_x\in {\mathscr{H}}_x({ℱ}_{Markov}(G_x,R_x))\subseteq {𝒪}_x$. For such objects $o\in 𝒪$, the terminating state probability measures correspond to the distribution of interest, i.e.:

$$\forall x\in 𝒳\forall s\in {𝒮}_x^f{P}_T(s\mid x)\propto R_x(s).$$

Similar to the unconditional case, we need to design a conditional loss function $ℒ$ on $𝒪$ that equals zero on such objects $o$ and only on those objects.

Definition 43.

Let $(𝒳,𝒢,ℛ,𝒪,\mathrm{\Pi },\mathscr{H})$ be a conditional GFlowNet. A conditional flow-matching loss is any function $ℒ:𝒪\to {ℝ}^{+}$ such that:

$$\forall o\in 𝒪ℒ(o)=0\iff \forall x\in 𝒳\exists F_x\in {ℱ}_{Markov}(G_x,R_x)o_x={\mathscr{H}}_x(F_x).$$ (76)

We say that $ℒ$ is condition-decomposable, if there are functions ${ℒ}_x:{𝒪}_x\to {ℝ}^{+}$ such that:

$$\forall o\in 𝒪ℒ(o)=\sum _{x\in 𝒳}{ℒ}_x(o_x)$$

We can obtain conditional flow-matching losses that are condition-decomposable starting from any family of flow-matching losses ${ℒ}_x$ for $(G_x,R_x,{𝒪}_x,{\mathrm{\Pi }}_x,{\mathscr{H}}_x)$ (Def. 26). In particular, we can choose these losses to be state-decomposable, edge-decomposable or trajectory-decomposable.

For instance, if each ${ℒ}_x={\sum }_{s\in {𝒮}_x}{L}_x(.,s)$ is state-decomposable, then the simultaneous sampling problem is cast to the following minimization problem:

$$\underset{o\in 𝒪}{\min }{𝔼}_{(x,s)\sim {\pi }_T}[L_x(o,s)],$$ (77)

where ${\pi }_T$ is any conditional full support distribution on $𝒳\times {\bigcup }_{x\in 𝒳}{𝒮}_x$, i.e. a probability distribution that satisfies:

$$\forall x\in 𝒳\forall s\in \bigcup _{x\in 𝒳}{𝒮}_x{\pi }_T(x,s)>0\iff s\in {𝒮}_x.$$

Such a conditional full support distribution can be obtained starting from any distribution ${\pi }_{𝒳}$ with full support on $𝒳$ and distributions ${\pi }_x$ with full support on ${𝒮}_x$ for any $x\in 𝒳$, as:

$${\pi }_T(x,s)=\{\begin{array}{cc}{\pi }_{𝒳}(x){\pi }_x(s)\text{if}s\in {𝒮}_x\hfill & \\ 0\text{otherwise}.\hfill & \end{array}$$

Example 9.

Consider the set:

$$𝒪=\{\widehat{F}:𝒳\times \bigcup _{x\in 𝒳}{𝔸}_x^{-f}\to {ℝ}^{+},\widehat{F}(s\to s^{\prime }\mid x)=0\text{if}s\to s^{\prime }\notin {𝔸}_x\}.$$

For each $x\in 𝒳$ the function ${\widehat{F}}_x:=\widehat{F}(.\mid x)$ is an element of ${𝒪}_{edge,x}$(Ex. 1), i.e. it is a function from ${𝔸}_x^{-f}$ to ${ℝ}^{+}$. Meaning that the set $𝒪$ can be seen as a function mapping each $x\in 𝒳$ to $F_x\in {𝒪}_{edge,x}$. Denoting by $\mathscr{H}:x\in 𝒳\to {\mathscr{H}}_{edge,x}$ and $\mathrm{\Pi }:x\in 𝒳\to {\mathrm{\Pi }}_{edge,x}$, we obtain a valid conditional flow parametrization $(𝒪,\mathscr{H},\mathrm{\Pi })$.

Instead of learning each function ${\widehat{F}}_x$ separately, this parametrization enables learning functions $\widehat{F}\in 𝒪$ of both the condition $x$ and the non-terminating edge $s\to s^{\prime }$, thus exploiting the generalization capabilities of machine learning algorithms not only on edges, but also on conditions.

Consider the functions $L_x$ defined for each $x\in 𝒳$ as:

$$L_x({\widehat{F}}_x,s^{\prime })=\{\begin{array}{cc}{\left(\log \left(\frac{\delta +{\sum }_{s\in Par(s^{\prime })}\widehat{F}(s\to s^{\prime }\mid x)}{\delta +R(s^{\prime }\mid x)+{\sum }_{s^{\prime \prime }\in Child(s^{\prime })\setminus \{s_f\mid x\}}\widehat{F}(s^{\prime }\to s^{\prime \prime }\mid x)}\right)\right)}^2\text{if}{s}^{\prime }\ne s_f,\hfill & \\ 0\text{otherwise}\hfill & \end{array}$$

where $\delta \ge 0$ is a hyperparameter. The function $ℒ$ mapping each $\widehat{F}\in 𝒪$ to

$$ℒ(\widehat{F})=\sum _{x\in 𝒳}\sum _{s\in {𝒮}_x}{L}_x({\widehat{F}}_x,s^{\prime }),$$

is a conditional flow-matching loss, that is both condition-decomposable and state-decomposable.

Appendix C Policies in Deterministic and Stochastic Environments

Until now, we have focused on a deterministic environment where state changes can perfectly be calculated from a given action. This makes sense when the actions are cognitive actions, internal to an agent, e.g., sequentially constructing a candidate solution to a problem (an explanation, a plan, an inferred guess, etc). What about the scenario where the actions are external and affect the real world? The outcome are likely to be only imperfectly predictable. To address this scenario, we will now extend the GFlowNet framework to learn a policy $\pi$ for an agent in an environment that could be deterministic or stochastic.

Definition 44.

A policy $\pi :𝒜\times 𝒮\mapsto ℝ$ is a probability distribution $\pi (a|s)$ over actions $a\in 𝒜$ for each state $s$. To denote the fact that the action space may be restricted based on $s$, we write $𝒜(s)$ for the valid actions in state $s$.

To denote the introduction of actions in the GFlowNet framework, we will decompose transitions in two steps: first an action $a_t$ is sampled according to a policy $\pi$ from state $s_t$, and then the environment transforms this (in a possibly stochastic way) into a new state $s_{t+1}$.

Definition 45.

We generalize the notion of state as follows: even states are of the form $s\in 𝒮$ while odd states are of the form $(s,a)\in 𝒮\times 𝒜$. The policy $\pi$ governs the transition from an even state to a compatible next odd state with $a\in 𝒜(s)$, while the environment $P(s_t\to s_{t+1}|s_t,a_t)$ governs the transition from an odd state to the next even state.

As a result of the above definition, the even-to-even transition is summarized by

$$P_F(s_{t+1}|s_t)=\sum _{a_t}P(s_t\to s_{t+1}|s_t,a_t)\pi (a_t|s_t).$$ (78)

Note that the detailed balance condition, which involves a backward transition $P_B$, will also be decomposed in two parts: (1) for inverting the even-to-odd transition,

$$P_B(s_t|(s_t,a_t))=1$$ (79)

by definition, and (2) for inverting the odd-to-even transition, we have to actually represent (and learn)

$$P_B((s_t,a_t)|s_{t+1}).$$ (80)

This conditional distribution incorporates the preference we may have over different paths leading to the same state while consistent with the environment $P(s_t\to s_{t+1}|s_t,a_t)$. The normalization constraint on ${\widehat{P}}_B$ can guarantee flow-matching via detailed balance, as argued around Eq. 31.

C.1 Known Deterministic Environments

A deterministic environment is perfectly controllable: we can choose the action that leads to the most desired next state, among the valid actions from the previous state. In the case where the environment is deterministic, we can directly apply the results of Sec. 3, as follows. At each time step $t$ and from state $s_t$, the agent picks an allowed action $a_t\in 𝒜(s_t)$ according to a policy $\pi (a_t|s_t)$. The set of allowed actions should coincide with those actions for which $\pi (a_t|s_t)>0$. Since the environment is deterministic and known, there is a deterministic function $T:𝒮\times 𝒜\mapsto 𝒮$ which gives us the next state $s_{t+1}=T(s_t,a_t)$. In that case, we can ignore the even/odd state distinction and identify the learnable policy $\pi$ with the learnable transition probability function $P_F(s_{t+1}|s_t)$ of GFlowNets, as follows.

Proposition 46.

In a deterministic environment and a GFlowNet agent with policy $\pi (a_t|s_t)$ and state transitions given by $s_{t+1}=T(s_t,a_t)$, the transition probability $P_F(s_{t+1}|s_t)$ is given by

$$P_F(s_{t+1}|s_t)=\sum _{a:T(s_t,a)=s_{t+1}}\pi (a|s_t).$$ (81)

Hence if only one action $a_t$ can transition from $s_t$ to $s_{t+1}=T(s_t,a_t)$, then

$$P_F(s_{t+1}|s_t)=\pi (a_t|s_t).$$ (82)

Proof The result is obtained by marginalizing over $a$:

$P_F(s_{t+1}|s_t)$ $={\displaystyle \sum _a}P(s_t\to s_{t+1},a_t|s_t)$

$={\displaystyle \sum _a}P(s_t\to s_{t+1}|s_t,a_t)\pi (a|s_t)$

$={\displaystyle \sum _{a:s_{t+1}=T(s_t,a)}}\pi (a|s_t)$

with $P(s_t\to s_{t+1}|s_t,a_t)=1_{s_{t+1}=T(s_t,a)}$.

The case with a single possible action to obtain the transition is obtained because the sum contains only one term.

Proposition 47.

In a deterministic environment with $s_{t+1}=T(s_t,a_t)$, a backwards transition probability function can be derived from a backwards policy ${\pi }_B$,

$$P_B(s_t|s_{t+1})=\sum _{a:T(s_t,a)=s_{t+1}}{\pi }_B(a|s_{t+1})$$ (83)

and in the case where a single action $a_t$ explains each transition $s_{t+1}=T(s_t,a_t)$,

$$P_B(s_t|s_{t+1})={\pi }_B(a_t|s_{t+1})$$ (84)

Proof The proof goes along exactly the same lines as for Prop. 46.

C.2 Unknown Deterministic Environments

If the environment is deterministic but unknown, we have to learn the transition function $T$ and we should also learn its inverse $T^{-1}$ which recovers the previous state given the next state and the action:

$$T^{-1}:𝒮\times 𝒜\mapsto 𝒮\mathrm{s}.\mathrm{t}.T^{-1}(T(s,a),a)=s.$$ (85)

Unfortunately, if the state and action spaces are discrete and in high dimension, learning $T$ and $T^{-1}$ in a way that generalizes to unseen transitions may be difficult and might be more easily achievable via a continuous relaxation. The methods for stochastic environments could be used for this purpose.

C.3 Stochastic Environments

The setting of stochastic environments is less straightforward but more general. We will decompose the transition as per Eq. 78 but not assume that $P(s_t\to s_{t+1}|s_t,a_t)$ is a dirac. The first thing to note is that we can still obtain a Markovian flow, but that we are not guaranteed to find a policy which matches the desired terminal reward function.

Proposition 48.

In a stochastic environment with environment transitions $P(s_t\to s_{t+1}|s_t,a_t)$, any policy $\pi (a|s)$ can yield a Markovian flow and it may not be possible to perfectly achieve desired flows $\widehat{F}(s\to s_f)=R(s)$.

Proof We obtain a flow by satisfying the flow-matching or detailed balance equations for both even and odd steps, which can always be done for the following reason. From the even states, we can define an edge flow

$$\widehat{F}(s_t\to (s_t,a_t))=\widehat{F}(s_t)\pi (a_t|s_t)$$

and the backwards transition is ${\pi }_B(s_t|(s_t,a_t))=1$. This leads to the intermediate state flow

$$\widehat{F}((s_t,a_t))=\widehat{F}(s_t\to (s_t,a_t))$$

since there is only one edge into $(s_t,a_t)$, the one starting at $s_t$ and taking action $a_t$. From the odd states, we have the edge flow

$\widehat{F}((s_t,a_t)\to s_{t+1})=$ $\widehat{F}((s_t,a_t))P(s_t\to s_{t+1}|s_t,a_t)$

$=$ $\widehat{F}(s_t)\pi (a_t|s_t)P(s_t\to s_{t+1}|s_t,a_t)$

with $P(s_t\to s_{t+1}|s_t,a_t)$ representing the environment, and we obtain the even state flow with the usual formula (Eq. 22)

$$\widehat{F}(s_{t+1})=\sum _{(s_t,a_t)}\widehat{F}((s_t,a_t)\to s_{t+1})=\sum _{(s_t,a_t)}\widehat{F}(s_t)\pi (a_t|s_t)P(s_t\to s_{t+1}|s_t,a_t).$$

If unknown, the environment transitions $P(s_t\to s_{t+1}|s_t,a_t)$ can be estimated in the usual supervised way by observing the triplets $(s_t,a_t,s_{t+1})$ and estimating transition probabilities that approximately maximize the empirical log-likelihood of these observations. However, whereas in a stationary environment the transitions $P(s_t\to s_{t+1}|s_t,a_t)$ do not depend on the policy, the backwards transitions ${\widehat{P}}_B((s_t,a_t)|s_{t+1})$ depend on the forward environment transition and on the state flows, i.e., on the policy. With enough training time and capacity, the forward and backward transitions can be made compatible, but as usual with GFlowNets, in a realistic settings the flow matching equations will not be perfectly achieved.

If there is enough capacity and training time (training to completion), we thus obtain a flow. Then, defining the transition probabilities by the sequential sampling of transitions from the even and odd steps above, we obtain a Markovian flow (Prop. 16).

To show that the desired terminal flows are not necessarily achievable, it is sufficient to identify a counter-example. Consider a terminal reward $R(s)>0$ while the environment transitions into $s$ have zero probability. In that case, no matter how we choose our policy, we cannot put the desired flow into state $s$.

Keep also in mind that in practice, even in a completely controllable environment, we will not be guaranteed to find a flow that matches the target terminal reward function simply because of finite capacity and finite training time for the GFlowNet.

Whereas with a deterministic environment for the GFlowNet, one can freely choose $P_B$ for non-terminal edges, it is not so for stochastic environments, as argued below. On even-to-odd transitions, $P_B(s|s_t,a_t)=1_{s=s_t}$ by construction. On odd-to-even transitions $(s_t,a_t)\to s_{t+1}$ the problem is that the forward transition is not a free parameter (it corresponds to the environment’s $P(s_t\to s_{t+1}|s_t,a_t)$).

Counterexample 49.

With a GFlowNet with a fixed stochastic environment, it may not be possible for $P_B(s_t,a_t|s_{t+1})$ to be chosen freely while also matching the flows and the terminal rewards.

For a counterexample, consider the setting in Fig. 6. Suppose $R(s^{\prime \prime \prime })\ne 0$ and the environment-provided transition $T(s^{\prime \prime \prime }|s^{\prime \prime },a^{\prime \prime })=0$. Then, in order to match the terminal reward $R(s^{\prime \prime \prime })$, we must require that $P_B(s^{\prime },a^{\prime }|s^{\prime \prime \prime })\ne 0$, which means that $P_B$ cannot be chosen freely.

Figure 6: Consider a simple counter-example with only two paths from $s$ to $s^{\prime \prime \prime }$, with a given $R(s^{\prime \prime \prime })>0$. One path goes through $s^{\prime }$ and the other through $s^{\prime \prime }$. From $s^{\prime }$ the only feasible action $a^{\prime }$ leads to $s^{\prime \prime \prime }$ and similarly from $s^{\prime \prime }$ with action $a^{\prime \prime }$. However, it may be that the environment probability $P(s^{\prime \prime \prime }|(s^{\prime \prime },a^{\prime \prime }))=0$, constraining $P_B((s^{\prime \prime },a^{\prime \prime })|s^{\prime \prime \prime })=0$. Therefore it may not possible to choose the backward transitions freely while matching the flows and terminal rewards.

Appendix D Expected Downstream Reward and Reward-Maximizing Policy

We have already introduced the probability distribution $P_T(s)$ and conditional probabilities $P_T(s^{\prime }|s\le s^{\prime })$ over terminating states (the states visited just before exiting into $s_f$). More generally, one can consider any distribution $P_{\pi }(s)$ over terminating states arising from some arbitrary choice of GFlowNet policy $\pi$ and compute the expected reward under this distribution:

Definition 50.

The expected reward after visiting state $s$ of a flow with terminal reward function $R$, under some distribution over terminating states $P$, is

$$V_{P_{\pi }}(s):=E_{P_{\pi }(S)}[R(S)|S\ge s]=\sum _{s^{\prime }\ge s}R(s^{\prime })P_{\pi }(s^{\prime }|s\le s^{\prime }).$$ (86)

Proposition 51.

When the probability distribution over terminating states is $P_T$ given by the flow (see Def. 13), the expected reward under $P_T$ is

$$V_{P_T}(s)=\frac{{\sum }_{s^{\prime }\ge s}R{(s^{\prime })}^2}{{\sum }_{s^{\prime }\ge s}R(s^{\prime })}.$$ (87)

Proof We apply the definition of conditional $P_T$ (Cor. 33) to Eq. 86 and obtain the result.

While we have a simple expression of the expected reward under $P_T$, the expected reward is defined more broadly for the distribution $P_{\pi }$ arising from any policy $\pi$. In particular for any policy $\pi (a|s)$, we can also define the expected reward $V_{P_{\pi }}$ under the distribution $P_{\pi }$ over terminating states induced by $\pi$. Expected rewards play a role similar to the state and state-action value functions in reinforcement learning, and as a consequence they also satisfy an equivalent of the policy improvement theorem when intermediate rewards are $0$ and the discount factor $\gamma =1$:

Proposition 52.

Let $P_{\pi }$ be a distribution over terminating states arising from a policy $\pi$, and $\overline{\pi }$ a greedy policy under the expected reward $V_{P_{\pi }}$, i.e.,

$\overline{\pi }(a|s)=0\mathrm{unless}$

$V_{P_{\pi }}((s,a))\ge V_{P_{\pi }}((s,a^{\prime }))\forall a^{\prime }.$ (88)

Then for all $s$

$$V_{P_{\overline{\pi }}}(s)\ge V_{P_{\pi }}(s).$$ (89)

That is, the expected reward under the probability induced by $\overline{\pi }$ is no worse than the one induced by $\pi$.

Proof Let us denote by $\overline{\pi }(s)$ the action deterministically chosen by greedy policy $\overline{\pi }$ from $s$, and $s_n$ the stochastically sampled terminating state. Then:

$V_{P_{\pi }}(s)$ $\le V_{P_{\pi }}((s,\overline{\pi }(s)))$

$=E_{\overline{\pi }}[V_{P_{\pi }}(s_{t+1})|s_t=s]$

$\le E_{\overline{\pi }}[V_{P_{\pi }}(s_{t+1},\overline{\pi }(s_{t+1}))|s_t=s]$

$=E_{\overline{\pi }}[{𝔼}_{\overline{\pi }}[V_{P_{\pi }}(s_{t+2})|s_t=s|s_t=s]$

$=E_{\overline{\pi }}[V_{P_{\pi }}(s_{t+2})|s_t=s]$

$\mathrm{\dots }$

$\le E_{\overline{\pi }}[V_{P_{\pi }}(s_n)|s_t=s]$

$=E_{\overline{\pi }}[R(s_n)|s_t=s]$

$=V_{P_{\overline{\pi }}}(s)$

where we have used the fact that, for all $s^{\prime }$, $V_{P_{\pi }}(s^{\prime })\le V_{P_{\pi }}(s^{\prime },\overline{\pi }(s^{\prime }))$ since $\overline{\pi }$ is a greedy policy.

An immediate consequence is the following:

Corollary 53.

There exists a policy ${\pi }^{\ast }(a|s)$ that maximizes the expected reward for all states $s$, namely the greedy policy of Prop. 52 associated with the GFlowNet’s policy $\pi$ yielding terminal distribution $P_T(s)=R(s)/Z$.

How do can we estimate the expected reward under $P_T$? We just need to train another flow (or set of heads for a GFlowNet, see Appendix F) with $R^2$ as the reward function.

Proposition 54.

Consider two flows $F$ and $F^{\prime }$, one matching terminal reward function $R$ and the other matching terminal reward function $R^2$. Then the expected reward under $P_T$ (the distribution over terminating states defined by the flow $F$) is

$$V_{P_T}(s)=\frac{F^{\prime }(s|s)}{F(s|s)}.$$ (90)

Proof We start from Eq. 87 of the above corollary and notice that the numerator is the self-flow for $F^{\prime }$ while the denominator is the self-flow for $F$ (see Eq. 44).

D.1 Preference for High-Reward Early Trajectory

We have seen in Sec. 2.6 that by imposing a particular preference on $P_B$, one can make the GFlowNet sampling mechanism prefer to construct states in some orders more than others, e.g., one could prefer to start with states with larger expected reward (over their potential continuations), using Def. 50. It suffices to define $P_B(s_t,a_t|s_{t+1})$ so that it puts more probability mass on state-action pairs $(s_t,a_t)$ with larger $V((s_t,a_t))$.

Appendix E Intermediate Rewards and Trajectory Returns

Up to now and in the GFlowNet paper (Bengio et al., 2021), we have considered terminal rewards as events happening only once per trajectory, at its end. Consider instead an agent experiencing a complete trajectory $\tau$ and declare its return to be the sum of some intermediate environment rewards associated with all the transitions into the sink node from each of the visited states.

Definition 55.

The trajectory return $\rho (\tau )$ associated with a partial trajectory $\tau =(s_i,s_{i+1}\mathrm{\dots },s_n,s_f)$ is defined as

$\rho (s_i,s_{i+1},\mathrm{\dots },s_n,s_f)$ $={\displaystyle \sum _{t=i}^n}R(s_t)={\displaystyle \sum _{t=i}^n}F(s_t\to s_f)$ (91)

$\rho (s_f)$ $=0$ (92)

and the expected future return $\overline{\rho }(s_t)$ associated with a state $s_t$ is defined as

$\overline{\rho }(s_t)$ $=E[\rho (s_t,s_{t+1},s_{t+2},\mathrm{\dots },s_n,s_f)|s_t]$

$={\displaystyle \sum _{s_{t+1},\mathrm{\dots },s_n}}P(s_{t+1},\mathrm{\dots },s_n,s_f|s_t)\rho (s_t,s_{t+1},\mathrm{\dots },s_f)$ (93)

where the expectation is defined under the flow’s probability measure over trajectories (conditioned on the trajectory going through $s_t$).

Proposition 56.

The expected future return $\overline{\rho }(s)$ achievable from trajectories starting at $s$ satisfies the following recursion:

$$\overline{\rho }(s)=R(s)+\sum _{s^{\prime }\in Child(s)}{P}_F(s^{\prime }|s)\overline{\rho }(s^{\prime }).$$ (94)

Proof From the definition of return (Def. 55), we obtain the following:

$\overline{\rho }(s_t)$ $={\displaystyle \sum _{s_{t+1},s_{t+2},\mathrm{\dots }}}P(s_{t+1},\mathrm{\dots },s_f|s_t)(R(s_t)+\rho (s_{t+1},s_{t+2},\mathrm{\dots },s_n))$

$=R(s_t)+{\displaystyle \sum _{s_{t+1}}}{P}_F(s_{t+1}|s_t){\displaystyle \sum _{s_{t+2},\mathrm{\dots },s_n}}P(s_{t+2},\mathrm{\dots },s_n|s_{t+1})\rho (s_{t+1},\mathrm{\dots },s_n)$

$=R(s_t)+{\displaystyle \sum _{s_{t+1}}}{P}_F(s_{t+1}|s_t)E[\rho (s_{t+1},s_{t+2},\mathrm{\dots })|s_{t+1}]$

$=R(s_t)+{\displaystyle \sum _{s_{t+1}}}{P}_F(s_{t+1}|s_t)\overline{\rho }(s_{t+1}).$

One reason why the above recursion is interesting is that it corresponds to the Bellman equation (Sutton and Barto, 2018) for the value function (which is the expected downstream return) in the case of no discounting (with an episodic setting). It is interesting to compare it with one of the equations we obtain for the state flow (Eq. 22):

$F(s_t)$ $={\displaystyle \sum _{s_{t+1}}}F(s_{t+1})P_B(s_t|s_{t+1})$

$=R(s_t)+{\displaystyle \sum _{s_{t+1}\ne s_f}}{P}_B(s_t|s_{t+1})F(s_{t+1})$

The two recursions are different: one uses the forward transition to propagate values backward, while the other uses the backward transitions to propagate flows.

Definition 57.

Let us denote $r(s)$ the possibly stochastic environment reward, provided when an agent visits state $s$ (and generally distinct from the GFlowNet terminal reward at $s$), and consider the environment reward accumulated in the partial trajectory $\tau =(s_0,s_1,\mathrm{\dots },s_n)$ leading to state $s_n=s$. Let us call the GFlowNet state return-augmented if the state $s$ includes the accumulated reward $\nu (s)$, i.e., there exists a function $\nu (s)={\sum }_{t=0}^{n}r(s_t)$. Let us call a GFlowNet with a return-augmented state a return-augmented GFlowNet.

We want to keep track of accumulated intermediate rewards in the GFlowNet state to compute the terminal reward from the state, and thus train GFlowNets to sample in proportion to the accumulated reward:

Proposition 58.

Suppose $G$ is a return-augmented GFlowNet with target terminal reward function $R(s)$ equal to the accumulated environment reward $\nu (s)$. Furthermore, suppose $G$ is trained to completion. Then sampling from $G$ produces accumulated reward $\nu$ with probability proportional to $\nu$.

Proof Since $G$ is a GFlowNet trained to completion with terminal reward function $R(s)=\nu (s)$, we know that sampling from $G$ samples terminal accumulated rewards $\nu (s)$ with probability proportional to $\nu$.

Note that such a GFlowNet can only be trained offline, i.e., using trajectories which have terminated and for which we have observed the return.

Note also that if we did not augment the state with the accumulated reward, then the GFlowNet terminal reward would not be a function of the state. Having a return-augmented state also makes it possible to handle stochastic environment rewards in the GFlowNet framework.

Appendix F Multi-Flows, Distributional GFlowNets, Unsupervised GFlowNets and Pareto GFlowNets

Consider an environment with stochastic rewards. As with Def. 57, we could augment the state to include the random accumulated rewards, thus (by Prop. 58) making the GFlowNet sample trajectories with returns $\rho$ occurring with probability proportional to $\rho$. However, similarly to Distributional RL (Bellemare et al., 2017), it could be interesting to generalize GFlowNets to capture not just the expected value of achievable terminal rewards but also other statistics of its distribution. More generally, we can think of this like a family of GFlowNets, each of which models in its flow a particular future environmental outcome of interest. With the particle analogy of GFlowNets, it would be as if the particles had a colour or label (just like the frequency of each photon in an group of photons travelling together in a beam of light) and that we separately account for the flows associated with all the possible label types.

If the number of outcomes (the number of possible labels) is small, this could be implemented with different output heads of the GFlowNet (e.g., one output for the flow associated with each label). When a trajectory associated with a particular label outcome is observed, the corresponding heads get gradients. A more powerful and general implementation puts the outcome event as an input of the GFlowNet, thus amounting to training a conditional GFlowNet (see Sec. 4.4), and formalized below.

Definition 59.

Let us define the outcome $y=f(s)$ as a known function $f(s)$ of the state $s$. An outcome can be whether an environment reward takes a particular value, or it can be a vector of important features of $s$ which are sufficient to determine many possible environment reward functions (in particular, $f$ can be the identity function). Let us call the conditional GFlowNet taking $y$ as conditioning input, with conditional flows $F(A|y)$ for events $A$ over the trajectories consistent with reward function

$$R_y(s)=1_{f(s)=y}$$ (95)

an outcome-conditioned GFlowNet with outcome function $f$.

We will limit ourselves here to a discrete set of outcomes for simplicity but expect that this approach can be generalized to a continuous set. Note how, if the outcome-conditioned GFlowNet is trained to completion, it makes it possible to only sample terminating states $s$ yielding the chosen outcome $y$. In principle, this allows sampling objects guaranteed to have a high reward (under the given reward function). In practice, a GFlowNet will never be perfectly trained to completion, and we should think of such an outcome-conditioned GFlowNet similarly to a goal-conditioned policy in RL (Ghosh et al., 2018) or the reward-conditioned upside-down RL (Schmidhuber, 2019). An interesting question for future work is to extend these outome-conditioned GFlowNets to the case of stochastic rewards or stochastic environments.

Definition 60.

A distributional GFlowNet is an outcome-conditioned GFlowNet taking as conditioning input the value of the environment return associated with complete trajectories. This can be achieved by making the GFlowNet return-augmented, so that the return can be read from the terminating state of the trajectory.

Training an outcome-conditioned GFlowNet can only be done offline because the conditioning input (e.g., the final return) may only be known after the trajectory has been sampled. A reasonable contrastive training procedure could thus proceed as follows:

1. Sample a trajectory ${\tau }^{+}$ according to an unconditional training policy ${\pi }_T$.

2. Obtain the outcome $y^{+}=f(s^{+})$ from the terminating state $s^{+}$ (occurring just before the sink state $s_f$ in ${\tau }^{+}$).

3. Update the conditional GFlowNet with ${\tau }^{+}$ and target terminating reward $R(s^{+}|y^{+})=1_{f(s^{+})=y^{+}}=1$.

4. Sample a trajectory ${\tau }^{-}$ according to the conditional GFlowNet policy with condition $y^{+}$.

5. Obtain the actual outcome $y^{-}=f(s^{-})$ for the terminating state $s^{-}$ (occurring just before the sink state $s_f$) in ${\tau }^{-}$. If the GFlowNet was perfectly trained, we should have $y^{+}=y^{-}$ but otherwise, especially if the number of possible outcomes is large, this becomes unlikely.

6. Update the conditional GFlowNet using a flow-matching loss with trajectory ${\tau }^{-}$ and target terminating reward $R(s^{-}|y)=1_{y^{-}=y^{+}}$ (likely to be 0 if there are many possible values for $y$).

An interesting question for future work is to consider a smoother reward function instead of the sharp but sparse reward $R(s|y)=1_{f(s)=y}$ as conditional reward, in order to make training easier.

F.1 Defining a reward function a posteriori

The reward function may not be known a priori or it may be known only up to some unknown constants (e.g., defining a Pareto front) or we may wish to generalize GFlowNets so they can be trained or pre-trained in a more unsupervised way, with the specific generative task only specified afterwards. The following important proposition allows us to do just that, and convert an outcome-conditioned GFlowNet into one that samples according to a given reward function, without having to retrain the network.

Proposition 61.

Consider an outcome-conditioned GFlowNet trained to completion with respect to the possible outcomes $y=f(s)$ over terminating states $s$ and a terminal reward function $R(s)=r(f(s))$ given a posteriori (possibly after training the GFlowNet) as a function $r(y)$ of the outcome $y=f(s)$. Then a GFlowNet with flow $F_{r\circ f}(A)$ over events $A$ which matches target terminal reward function $R=r\circ f$ can be obtained from the flow $F(A|y)$ of the outcome-conditioned GFlowNet via

$$F_{r\circ f}(A)=\sum _{y}r(y)F(A|y).$$ (96)

For example, the GFlowNet policy ${\pi }_{r\circ f}(a|s)$ for terminal reward $R=r\circ f$ can be obtained from

$${\pi }_{r\circ f}(a|s)=\frac{{\sum }_{y}r(y)F((s,a)|y)}{{\sum }_{y}r(y)F(s|y)}$$ (97)

or

$${\pi }_{r\circ f}(a|s)=\frac{{\sum }_{y}r(y)F(s|y)\pi (a|s,y)}{{\sum }_{y}r(y)F(s|y)}$$ (98)

where $F(s|y)$, $F((s,a)|y)$ and $\pi (a|s,y)$ are the outcome-conditioned state flow, state-action flow and action policy respectively.

Proof We first clarify what $F_{r\circ f}(A)$ means:

$$F_{r\circ f}(A)=\sum _{s}r(f(s))P(A|s\to s_f)=\sum _{s}r(f(s))P_B(A|s\to s_f)$$ (99)

where $P(A|s\to s_f)$ is the probability of event $A$ (e.g., a particular state or transition) among the trajectories that end in terminal transition $s\to s_f$, and it can be determined entirely by $P_B$ (considering all the backward paths starting at $s$ and going back to the particular state or transition $A$). This means that $P(A|s\to s_f)$ does not depend of the choice of reward function, so Eq. 99 is valid for any $r$.

Clearly, we see that Eq. 96 works in the extreme case where $r(y^{\prime })=1_{y=y^{\prime }}$ is an indicator function at some specific $y$ value since the sum in Eq. 96 reduces to $F(A|y)$ which corresponds to reward function $R_y$ as per Eq. 95. This yields

$$F_{1_y}(A)=F(A|y)=\sum _s{1}_{f(s)=y}P(A|s\to s_f)$$ (100)

where $1_y$ denotes the function that, given $s$, returns $1_{y=s}$.

To complete the proof let us start with the right-hand side of Eq. 96 and insert the above definition of $F(A|y)$ (Eq. 100), then swap the sums and use the indicator function to cancel the sum over $y$, and finally apply the definition of $F_{r\circ f}(A)$ in Eq. 99:

${\displaystyle \sum _y}r(y)F(A|y)$ $={\displaystyle \sum _y}r(y){\displaystyle \sum _s}{1}_{f(s)=y}P(A|s\to s_f)$

$={\displaystyle \sum _s}{\displaystyle \sum _y}r(y)1_{f(s)=y}P(A|s\to s_f)$

$={\displaystyle \sum _s}r(f(s))P(A|s\to s_f)$

$=F_{r\circ f}(A)$ (101)

recovering the left-hand side of Eq. 96 as desired.

This makes it possible to predict probabilities and perform sampling actions for all the possible outcomes $y$ arising in different states (in the extreme where we know nothing about possible reward functions and the outcome is $y=s$), and then convert that GFlowNet on the fly to one specialized to a given terminal reward function $R=r\circ f$. However, we note that it requires more computation for each action at run-time: we have to perform these sums (possibly via Monte-Carlo integration) over the outcome space, and there may be a computational time versus accuracy trade-off in the resulting decisions (based on how many Monte-Carlo samples are used to approximate the above sums).

F.2 Pareto GFlowNets

A related application of these ideas concerns Pareto optimization, where we are not sure about the correct reward function up to a few coefficients forming a convex combination of underlying objectives.

Definition 62.

The Pareto additive terminal reward functions can be written as

$$R_{\omega }(s)=\sum _i{\omega }_i{f}_i(s)$$ (102)

where $\omega \in \{\omega \in W\subset {ℝ}^d:{\omega }_i\ge 0,{\sum }_i{\omega }_i=1\}$ are convex weights and the outcomes of interest are $d$ objectives $y_i=f_i(s)$, and $W$ is a discrete set of convex weights.

Definition 63.

The Pareto multiplicative terminal reward functions can be written as

$$R_{\omega }(s)=e^{-{\sum }_i{\omega }_i{e}_i(s)}$$ (103)

where $\omega \in \{\omega \in W\subset {ℝ}^d:{\omega }_i\ge 0,{\sum }_i{\omega }_i=1\}$ are convex weights and the outcomes of interest are $d$ objectives $y_i=f_i(s)$, and $W$ is a discrete set of convex weights.

In these cases, we can train a conditional GFlowNet with $\omega$ as conditioning input and $R_{\omega }(s)=R(s|\omega )$ as conditional terminal reward function. At run-time, we can scan the set of $\omega$’s in order to obtain different policies or predicted probabilities or free energies. The above can easily be generalized to a non-convex and non-linear combination of the objectives, so long as the combined objective is parametrized by $\omega$. Note the similarity between this idea (and more generally outcome-conditioned GFlowNets) and the earlier work by Dosovitskiy and Djolonga (2019). The same idea of conditioning by a form of specification of the loss can be applied to GFlowNets to obtain a family of GFlowNets, one for each variant of the loss function.

A useful application of a Pareto GFlowNet with such reward functions is to draw samples from the Pareto frontier. Once the Pareto GFlowNet is trained, we can draw samples from the Pareto frontier by first sampling the convex weights $\omega$ and then sampling trajectories. This can be useful in multi-objective optimization or sampling, where we want to draw a diversity of solutions corresponding to different trade-off points of the various objectives.

We could also train an outcome-driven GFlowNet by providing the vector of objective values $y$ as input and exploit prior knowledge about the objectives. For example, we may believe that different objectives can be modeled independently of each other and that $F(s|y)$ (or similarly for $F((s,a)|y)$) can be written as a basis expansion $F(s|y)={\sum }_{i=1}^d{\sum }_{j=1}^N{\varphi }_j(y_i)F_{i,j}(s)$ with $N$ bases ${\varphi }_j$ used to represent the objective $y_i$. This could be advantageous from a generalization point of view if learning about the different objectives should be disentangled from one another (e.g., one is stationary and the other is not).