---
title: Generalized Poisson Flow (GPFlow) Overview
url: https://www.emergentmind.com/topics/generalized-poisson-flow-gpflow
type: topic
---

# Generalized Poisson Flow (GPFlow) Overview

Generalized Poisson Flow (GPFlow) denotes a family of generative constructions organized around Poisson-type dynamics, but the term is used in two closely related senses in the recent literature. In the fixed-dimensional image-generation line, it refers to generalizations of Poisson Flow Generative Models (PFGM), where data are treated as electric charges in an augmented space and generation follows deterministic ODEs induced by a Poisson field [2209.11178]. In the variable-length protein-design line, GPFlow is a specific framework that learns the rate function of an inhomogeneous generalized Poisson process and couples length insertions to within-length continuous, discrete, or Riemannian dynamics, thereby modeling joint distributions over both length and content [2607.09039]. The common thread is the replacement of a purely Gaussian, fixed-dimensional forward corruption process by transport laws derived from Poisson equations, generalized Poisson processes, or their limiting relations.

## 1. Terminological scope and historical development

The modern point of departure is the 2022 formulation of PFGM, which proposes a generative model that maps a uniform distribution on a high-dimensional hemisphere into an arbitrary data distribution by interpreting data points as electrical charges on the \(z=0\) hyperplane in an augmented space and following electric field lines generated by the solution to a Poisson equation [2209.11178]. In that formulation, the additional coordinate \(z\) is not an auxiliary bookkeeping device: it anchors the backward ODE, and the sampling trajectory terminates when \(z\) reaches zero.

Subsequent work generalized this construction in two directions. First, PFGM++ extends the augmented space from one extra coordinate to \(D\) extra coordinates and treats \(D\) as a design parameter. The 2024 Poisson Flow Consistency Models paper explicitly frames this as a broader family of Poisson-type flows parameterized by \(D\), and shows that consistency distillation can be ported from diffusion probability-flow ODEs to PFGM++ by a change of variables and an updated noise distribution [2402.08159]. Second, the 2025 Poisson Flow Joint Model for multiphase CT treats \(D\) as a task-specific generation-path parameter and uses PFGM++ as the backbone for conditional joint modeling of routine-dose and low-dose multiphase CT data [2502.03079].

A common misconception is to treat all usages of “GPFlow” as simple renamings of PFGM++. The literature does not support that identification. In the image-generation line, “generalized” refers to extensions of electrostatic Poisson flow in augmented Euclidean space, especially through the augmentation dimension \(D\). In the 2026 protein-design line, “Generalized Poisson Flow” is the name of a variable-length probability-flow framework whose central object is an inhomogeneous generalized Poisson process over sequence length rather than an electrostatic field over a fixed-dimensional ambient space [2607.09039]. This suggests that GPFlow is best understood as a family resemblance rather than a single algorithm.

## 2. Electrostatic Poisson flow in augmented space

In the original PFGM construction, data live in \(\mathbb{R}^N\) and are embedded into an augmented space \(\mathbb{R}^{N+1}\) through
\[
\tilde{x} = (x,z),
\]
with the data distribution represented as a charge sheet on the hyperplane \(z=0\):
\[
\tilde{p}(\tilde{x}) = p(x)\,\delta(z).
\]
The corresponding Poisson equation in the augmented space is
\[
\nabla^2 \varphi(\tilde{x}) = -\tilde{p}(\tilde{x}),
\]
and the electric field is the negative gradient of the potential,
\[
E(\tilde{x}) = -\nabla \varphi(\tilde{x}) = \frac{1}{S_N(1)}\int \frac{\tilde{x}-\tilde{y}}{\|\tilde{x}-\tilde{y}\|^{N+1}}\,\tilde{p}(\tilde{y})\,d\tilde{y}.
\]
The model’s base distribution is a uniform distribution on the upper hemisphere \(S_N^+(r)\) of radius \(r\) in \(\mathbb{R}^{N+1}\), and a central theorem states that, under compact-support and regularity assumptions, backward flow from a uniform hemisphere recovers the data distribution in the limit \(r\to\infty\) [2209.11178].

The associated flow is deterministic. In augmented coordinates \(\tilde{x}(t)=(x(t),z(t))\), the forward and backward ODEs are
\[
\frac{d\tilde{x}}{dt} = E(\tilde{x}),
\qquad
\frac{d\tilde{x}}{dt} = -E(\tilde{x}).
\]
A key structural fact is rescalability: replacing \(E(x)\) by \(f(x)E(x)\) for any positive scalar field \(f(x)\) preserves the same trajectories up to time reparameterization. PFGM exploits this by learning only the normalized field direction,
\[
u(\tilde{x}) := -\sqrt{N}\,\frac{\hat{E}(\tilde{x})}{\|\hat{E}(\tilde{x})\|_2},
\]
rather than the full field magnitude. This makes the learned object a direction field whose integral curves coincide with those of the physical Poisson field.

Sampling is performed with a backward ODE anchored by the physically meaningful coordinate \(z\). After reparameterizing by \(z\) and then by \(\log z\), the model uses
\[
\frac{d(x,z)}{dt'} = \big(f_{x,\theta}(\tilde{x})\,f_{z,\theta}(\tilde{x})^{-1}z,\; z\big),
\]
so that \(z\) decays exponentially and the stopping criterion is simply \(z\approx 0\). This anchoring is central to the original claim that samples “hit the unaugmented data manifold when the \(z\) reaches zero” [2209.11178].

The significance of this formulation is empirical as well as geometric. On CIFAR-10, PFGM reports an Inception score of \(9.68\) and an FID score of \(2.35\), described as current state-of-the-art performance among normalizing flow models, while also performing on par with state-of-the-art SDE approaches and offering \(10\times\) to \(20\times\) acceleration on image generation tasks. The same paper also reports tolerance to estimation errors on a weaker network architecture and robustness to the step size in the Euler method [2209.11178].

## 3. PFGM++, augmentation dimension \(D\), and the diffusion limit

PFGM++ generalizes the augmented-space construction by embedding data into \(\mathbb{R}^{N+D}\) rather than \(\mathbb{R}^{N+1}\). With
\[
\tilde{\mathbf{x}} = (\mathbf{x},\mathbf{0}) \in \mathbb{R}^{N+D},
\qquad
\tilde{\mathbf{x}}_\sigma = (\mathbf{x}_\sigma,\mathbf{z}) \in \mathbb{R}^{N+D},
\]
the data are again treated as electric charges on the hyperplane \(z_1=\cdots=z_D=0\), and the electrostatic field is
\[
\mathbf{E}(\tilde{\mathbf{x}}_\sigma) =
\frac{1}{S_{N+D-1}(1)}
\int
\frac{\tilde{\mathbf{x}}_\sigma-\tilde{\mathbf{x}}}
{\|\tilde{\mathbf{x}}_\sigma-\tilde{\mathbf{x}}\|^{N+D}}
\,p(\mathbf{x})\,d\mathbf{x}.
\]
Because of rotational symmetry in the augmented coordinates, it is sufficient to track the data component \(\mathbf{x}_\sigma\) and the radius
\[
r=\|\mathbf{z}\|_2.
\]
The Poisson flow ODE then reduces to
\[
d\mathbf{x}_\sigma
=
\mathbf{E}(\tilde{\mathbf{x}}_\sigma)_{\mathbf{x}_\sigma}\;
E(\tilde{\mathbf{x}}_\sigma)_r^{-1}\;dr,
\]
which defines a bijective mapping between the data manifold at \(r=0\) and a noise hypersurface at \(r=r_{\max}\) [2402.08159].

Training is recast through a perturbation kernel rather than direct field evaluation:
\[
p_r(\mathbf{x}_\sigma|\mathbf{x})
\propto
\frac{1}{\big(\|\mathbf{x}_\sigma-\mathbf{x}\|_2^2+r^2\big)^{\frac{N+D}{2}}}.
\]
The objective is
\[
\mathbb{E}_{r\sim p(r)}
\mathbb{E}_{\mathbf{x}\sim p(\mathbf{x})}
\mathbb{E}_{\mathbf{x}_\sigma \sim p_r(\mathbf{x}_\sigma|\mathbf{x})}
\left[
\left\|
s_\phi(\tilde{\mathbf{x}}_\sigma)
-\frac{\mathbf{x}_\sigma-\mathbf{x}}{r/\sqrt{D}}
\right\|_2^2
\right],
\]
whose minimizer is
\[
s_\phi^*(\tilde{\mathbf{x}}_\sigma)
=
\sqrt{D}\,
\mathbf{E}(\tilde{\mathbf{x}}_\sigma)_{\mathbf{x}_\sigma}
E(\tilde{\mathbf{x}}_\sigma)_r^{-1}.
\]
Thus the network directly approximates the ODE vector field \(d\mathbf{x}/dr\) up to the factor \(\sqrt{D}\) [2402.08159].

The role of \(D\) is the decisive generalization lever. With the change of variables
\[
r=\sigma\sqrt{D},
\]
PFGM++ can be written in the same generic probability-flow ODE form used in diffusion models,
\[
d\mathbf{x} = \Phi(\mathbf{x},\sigma,\boldsymbol{\phi})\,dt,
\]
with \(\Phi_{\text{PFGM++}}(\mathbf{x},\sigma,\boldsymbol{\phi}) := s_\phi(\tilde{\mathbf{x}})\). The same line of work states that, with \(r=\sigma\sqrt{D}\) and \(D\to\infty\), the PFGM++ objective and sampling converge to EDM, while \(D=1\) recovers the original PFGM [2402.08159]. The 2025 CT paper reiterates this as the central unification result: PFGM++ becomes PFGM for \(D=1\) and the diffusion model for \(D\rightarrow\infty\) [2502.03079].

This parameterization gives GPFlow, in the fixed-dimensional sense, a continuum of behaviors rather than a binary choice between Poisson flow and diffusion. Small finite \(D\) yields non-Gaussian, heavy-tailed perturbations and different robustness properties; large \(D\) yields diffusion-like behavior. The literature repeatedly interprets this as a robustness–rigidity trade-off.

## 4. Consistency formulations and low-dose CT denoising

Poisson Flow Consistency Models (PFCM) combine PFGM++ with the consistency-model framework, yielding a family of Poisson-flow-based consistency models that support single-step or few-step sampling [2402.08159]. The key formal move is to generalize consistency distillation from diffusion probability-flow ODEs to Poisson-flow ODEs by modifying both the solver update and the perturbation distribution. In place of Gaussian perturbations, PFCM use the PFGM++ heavy-tailed kernel, which in radial form satisfies
\[
p_r(R)\propto \frac{R^{N-1}}{(R^2+r^2)^{\frac{N+D}{2}}}.
\]

The learned consistency function maps noisy states at scale \(\sigma\) to the near-clean endpoint:
\[
f_{\boldsymbol{\theta}}(\mathbf{x},\sigma)
=
c_{\text{skip}}(\sigma)\mathbf{x}
+
c_{\text{out}}(\sigma)F_{\boldsymbol{\theta}}(\mathbf{x},\sigma),
\]
with boundary condition \(f(\mathbf{x},\epsilon)\approx \mathbf{x}\). In the conditional CT formulation, the low-dose image \(\mathbf{y}\) is added as an extra input channel, and \(f_{\boldsymbol{\theta}}(\mathbf{x}_\sigma,\sigma,\mathbf{y})\) directly approximates the posterior mapping associated with \(p(\mathbf{x}|\mathbf{y})\). Distillation uses a reverse PF ODE step
\[
\breve{\mathbf{x}}_{\sigma_i}
=
\mathbf{x}_{\sigma_{i+1}}
+
(\sigma_i-\sigma_{i+1})\,s_{\boldsymbol{\phi}}(\tilde{\mathbf{x}}_{\sigma_{i+1}}),
\]
and a consistency matching loss with LPIPS distance and an EMA teacher.

The CT application is explicitly described as a “task-specific” sampler that “hijacks” the generative process by replacing an intermediate state with the low-dose CT image. The paper also states that this introduces a severe mismatch because the noise characteristics of low-dose CT images are different from those of intermediate states in the Poisson flow process. The central empirical claim is that small \(D\) mitigates this mismatch: PFCM inherit from PFGM++ the ability to trade off robustness for rigidity via the hyperparameter \(D\in(0,\infty)\), and the resulting low-\(D\) samplers are markedly more robust than standard consistency models in this conditional setting [2402.08159].

On the Mayo low-dose CT validation set, the reported mean \(\pm\) std metrics are: LDCT input, LPIPS \(0.152 \pm 0.02\), SSIM \(0.94 \pm 0.02\), PSNR \(41.5 \pm 1.6\); EDM at \(79\) NFE, LPIPS \(0.053 \pm 0.01\), SSIM \(0.97 \pm 0.01\), PSNR \(43.8 \pm 1.2\); CD at \(1\) NFE, LPIPS \(0.065 \pm 0.01\), SSIM \(0.94 \pm 0.01\), PSNR \(42.0 \pm 0.8\); PFGM++ with \(D=4096\) at \(79\) NFE, LPIPS \(0.055 \pm 0.01\), SSIM \(0.96 \pm 0.01\), PSNR \(43.6 \pm 1.2\); and PS-PFCM with \(D=4096\) at \(1\) NFE, LPIPS \(0.061 \pm 0.01\), SSIM \(0.96 \pm 0.01\), PSNR \(43.0 \pm 1.0\). The same study reports that the \(1\)-NFE sampler reduces CT validation processing from approximately \(4\) hours to approximately \(3\) minutes on a single A6000 GPU, and that unconditional PFCM on CIFAR-10 reaches FID \(\approx 3.86\) at \(1\) NFE [2402.08159].

## 5. Joint conditional generalized Poisson flow for multiphase CT

The Poisson Flow Joint Model (PFJM) is a conditional, joint instantiation of PFGM++ for low-dose multiphase contrast-enhanced CT, where the data variable is the joint routine-dose triplet
\[
y=[y^1,y^2,y^3]\in \mathbb{R}^{L\times W\times 3}
\]
and the conditioning signal is the corresponding low-dose triplet
\[
c=[c^1,c^2,c^3]\in \mathbb{R}^{L\times W\times 3}.
\]
The model treats the full triplet as a single datum in the generalized Poisson-flow framework, with effective data dimension \(N=L\times W\times 3\), and learns the conditional mapping associated with
\[
p(y^1,y^2,y^3\mid c^1,c^2,c^3)
\]
rather than modeling each phase independently [2502.03079].

Training follows the PFGM++ perturbation design. A noise scale \(\sigma\) is sampled, the augmented-space radius is set by
\[
r\leftarrow \sigma\sqrt{D},
\]
a scalar perturbation \(R\sim p_r(R)\) and a random direction \(v=u/\|u\|_2\) with \(u\sim\mathcal{N}(0,I)\) are drawn, and the perturbed joint image is
\[
[\hat{y}^1,\hat{y}^2,\hat{y}^3]
=
[y^1,y^2,y^3] + Rv.
\]
The supervised conditional loss is
\[
\mathcal{L}(\theta)
=
\lambda(\sigma)\,
\bigl\|
f_\theta([\hat{y}^1,\hat{y}^2,\hat{y}^3],\sigma,[c^1,c^2,c^3])
-
[y^1,y^2,y^3]
\bigr\|_2^2.
\]
This differs from the original PFGM++ objective by predicting the clean joint routine-dose image directly and by conditioning explicitly on the joint low-dose image.

Sampling adds a second conditioning mechanism at the path level. The initial state is set to the condition,
\[
[x^1,x^2,x^3]^{(0)} \gets [c^1,c^2,c^3],
\]
and at each ODE step the current state is pulled toward the condition by
\[
[x^1,x^2,x^3]^{(n)}
\leftarrow
(1-w)[x^1,x^2,x^3]^{(n)} + w[c^1,c^2,c^3].
\]
The update direction is then computed from the network output, followed by a trapezoidal correction. The paper describes this as stronger than merely conditioning the network, because conditioning enters both the denoiser/vector field and the generation path itself.

PFJM treats \(D\) as a task-specific generation-path parameter and evaluates \(D\in\{2,8,32,64,128,256,512,2048\}\). The reported optimum is \(D=128\), with performance improving rapidly from \(D=2\) to \(D=128\) and then degrading gradually for larger \(D\). Across all tested \(D\) and phases, PFJM outperforms PFGM++, with average gains of \(-7.58\) HU in MAE, \(+3.97\%\) in SSIM, and \(+2.40\) dB in PSNR [2502.03079].

The headline quantitative results for PFJM averaged over phases are MAE \(8.99\) HU, SSIM \(98.75\%\), and PSNR \(48.24\) dB. The reported FID over joint 3-channel images is \(4.30\), compared with \(6.34\) for PFGM++ with \(D=128\), \(10.87\) for DDIM with \(20\) steps, and \(27.23\) for PFGM. PFJM also uses \(10\) ODE steps, whereas the diffusion baselines use \(20\) and PFGM uses \(30\) [2502.03079]. In this sense, PFJM is a concrete demonstration that generalized Poisson flow can support joint, strongly conditional medical-image generation with fewer deterministic solver steps than the compared diffusion baselines.

## 6. Variable-length GPFlow for protein design

In the 2026 protein-design framework, GPFlow is no longer a fixed-dimensional electrostatic flow. It is a variable-length generative model over
\[
\mathcal{Y}=\bigcup_{k\ge 0}\{k\}\times S^k,
\]
where \(k\) is the current protein length and \(S\) is the component space for residues, such as \(\mathbb{R}^3\) for coordinates, a categorical alphabet for amino acids, \(\mathrm{SO}(3)\) for rotations, or \(\mathbb{T}\) for torsion angles. Length evolves through a continuous-time pure birth process \(X_t\) with rate function \(\lambda_t(k)\), and content evolves through within-length generators \(B_t\) appropriate to the modality [2607.09039].

For a fixed target length \(L\), GPFlow introduces a scheduler \(\kappa_t\) and defines
\[
h_t=
\begin{cases}
\frac{\dot{\kappa}_t}{1-\kappa_t}, & t<\tau,\\
0, & t\ge \tau,
\end{cases}
\qquad
\lambda_t(X_t\mid L)=(L-X_t)\,h_t,
\]
with completion time \(\tau=\inf\{t:\kappa_t=1\}\). Under this construction, \(X_t\sim\mathrm{Binomial}(L,\kappa_t)\) for \(t<\tau\), and \(X_t\to L\) almost surely as \(t\to\tau^{-}\). For an arbitrary target length distribution \(q(k)\), the marginal rate is obtained by posterior averaging,
\[
\lambda_t^*(k)
=
\mathbb{E}_{X_1\sim p(\cdot\mid X_t=k)}
\big[\lambda_t(k\mid X_1)\big].
\]
Theorem 3.2 states that the generalized Poisson process driven by \(\lambda_t^*(k)\) has marginals \(p_t\) matching the posterior mixture over conditional paths and satisfies \(p_1=q\).

A central technical contribution is the simulation-free likelihood objective. For a realization with event times \(\{t_k\}_{k=1}^L\), the negative log-likelihood of the generalized Poisson trajectory is
\[
\mathcal{L}_{\mathrm{GP}}
=
\int_0^1 \lambda_t(X_t)\,dt
-
\sum_{k=1}^{L}\log \lambda_{t_k}(X_{t_k^-}).
\]
After marginalizing event times analytically, the rate-learning objective becomes
\[
\mathcal{L}_{\mathrm{GP}}
=
\mathbb{E}_{t,X_1\sim q(X)}
\big[
\lambda_\theta(X_t,t)
-
h_t(X_1-X_t)\log \lambda_\theta(X_t,t)
\big].
\]
This trains the marginal rate network \(\lambda_\theta\) without simulating the event process during training.

The multimodal extension couples length insertions to within-length dynamics. For current state \(Y_t=(X_t,S_t)\), insertions are represented by an insertion variable \(m\) and map \(I(y,m)\), while the conditional jump measure is
\[
Q_t(A\mid y,Y_1)
=
\lambda_t(y\mid Y_1)
\int 1\{I(y,m)\in A\}\,\rho_t(dm\mid y,Y_1).
\]
The corresponding marginal objects are
\[
B_t^* f(y) := \mathbb{E}[B_t f(y)\mid Y_t=y],
\qquad
\lambda_t^*(y) := \mathbb{E}[\lambda_t(y\mid Y_1)\mid Y_t=y],
\]
\[
\rho_t^*(dm\mid y)
:=
\frac{\mathbb{E}[\lambda_t(y\mid Y_1)\rho_t(dm\mid y,Y_1)\mid Y_t=y]}
{\mathbb{E}[\lambda_t(y\mid Y_1)\mid Y_t=y]}.
\]
Training combines the rate loss with a reconstruction term for insertion sampling and a within-length flow-matching term,
\[
\mathcal{L}
=
\mathcal{L}_{\mathrm{GP}}
+
w_{\mathrm{rec}}\mathcal{L}_{\mathrm{rec}}
+
w_{\mathrm{FM}}\mathcal{L}_{\mathrm{FM}}.
\]

The associated theory is formulated at the level of general Markov generators. For a process with generator
\[
A_t f(x)
=
u_t(x)\cdot\nabla f(x)
+
\tfrac{1}{2}\mathrm{Tr}\big(\Sigma_t(x)\nabla^2 f(x)\big)
+
\int \big(f(y)-f(x)\big)\,Q_t(dy\mid x),
\]
Theorem 4.1 gives the bound
\[
D_{\mathrm{KL}}(p_T\Vert\tilde{p}_T)
\le
\int_0^T
\mathbb{E}_{x\sim p_t}[L_t(x)]\,dt,
\]
where
\[
L_t(x)
=
\frac{\|u_t(x)-\tilde{u}_t(x)\|^2}{2\sigma_t^2}
+
D_{\mathrm{KL}}\big(Q_t(\cdot\mid x)\Vert \tilde{Q}_t(\cdot\mid x)\big).
\]
For GPFlow, the jump-measure divergence decomposes as
\[
D_{\mathrm{KL}}(Q_t^*\Vert Q_{\theta,t})
=
D_{\mathrm{KL}}(\lambda_t^*\Vert\lambda_{\theta,t})
+
\lambda_t^*(y)\,D_{\mathrm{KL}}(\rho_t^*\Vert\rho_{\theta,t}),
\]
and the paper states that, at the global optimum where these terms vanish, the learned process recovers the joint multimodal distribution over variable-length proteins.

Empirically, GPFlow is evaluated on unconditional structure design, unconditional sequence design, structure-based motif scaffolding, sequence-based motif scaffolding, and peptide co-design. In unconditional structure design, it raises designability on PDB from \(77.6\%\) for Proteina (60M) to \(96.1\%\) for GPFlow (65M), and on AFDB from \(90.6\%\) to \(94.5\%\), while recovering the empirical length distribution with mean \(\approx 153.7\) versus empirical \(\approx 154\). In unconditional sequence design, it matches UniRef50 mean pLDDT more closely than DPLM, with a difference of \(3.17\) versus \(14.40\), while maintaining approximately \(17\%\) higher diversity than DPLM. In structure-based motif scaffolding, it ranks first on \(10\) of \(16\) tasks, with average unique successes \(117.4\) and total unique successes \(1878\). In sequence-based motif scaffolding, it passes \(10/17\) tasks versus \(6/17\) for DPLM and \(5/17\) for EvoDiff. In peptide co-design, despite not using native length conditioning, it improves over PepFlow in AAR (\(53.71\%\) versus \(46.35\%\)), RMSD (\(1.39\)Å versus \(1.99\)Å), designability (\(48.19\%\) versus \(37.27\%\)), and diversity (\(0.443\) versus \(0.390\)) [2607.09039].

In this protein-design usage, GPFlow refers specifically to probability flows over variable-dimensional state spaces in which length emerges as the realized count of insertion events. That distinguishes it sharply from the fixed-dimensional electrostatic construction of PFGM and PFGM++, even though both lines of work use Poisson-type objects to define generative trajectories.

Source: https://www.emergentmind.com/topics/generalized-poisson-flow-gpflow