FES-FM: Free Energy Surface Sampling via Reduced Flow Matching

Tiejun Li; Zichen Liu

REVIEW 4 major objections 4 minor 1 cited by

FES-FM samples the free energy surface directly by learning a transport map in collective-variable space, so generation only requires integrating a low-dimensional ODE.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 15:04 UTC pith:WE3VNBB7

load-bearing objection A plausible reduced flow-matching sampler for free energy surfaces, but the core optimisation claim is unproven and the abstract overpromises. the 4 major comments →

arxiv 2605.00337 v2 pith:WE3VNBB7 submitted 2026-05-01 cs.LG

FES-FM: Free Energy Surface Sampling via Reduced Flow Matching

Zichen Liu , Tiejun Li This is my paper

classification cs.LG

keywords free energy surfaceflow matchingcollective variablesBoltzmann samplingreduced dynamicstransport equationJarzynski equalityHessian-informed prior

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tries to establish that the free energy surface of a high-dimensional Boltzmann system can be sampled by a learned transport map living entirely in the low-dimensional space of collective variables (CVs), without simulating the full configuration space at generation time. The core object is the reduced transport equation (20), which relates the velocity of the reduced dynamics to conditional expectations—the local mean force and the time derivative of the interpolated potential—along the level sets of the CV map. To train this velocity, the paper pairs a squared-loss regression of the local mean force with a residual of the transport equation, estimating both expectations through non-equilibrium trajectories and Jarzynski weights. The claim is that this reduces sampling to integrating a short, low-dimensional ODE, which is far cheaper than projecting full-space trajectories, and that the reconstructed free energy surface is obtained directly from those generated CV samples.

Core claim

The paper's central claim is that there exists a time-dependent velocity field u(y,t) in CV space whose ODE pushes the projected prior distribution onto the target CV density ρ(y), and that this u can be trained without ever computing the high-dimensional target density. The training objective combines a neural velocity u, an auxiliary mean-force network v, and a scalar network c to satisfy the reduced transport equation (20), in which the coefficients are conditional expectations over the CV level set Σy. Combined with the Hessian-informed harmonic prior of Section 5, which gives an E(3)-invariant and physically meaningful starting distribution for many-particle systems, the method claims t

What carries the argument

The reduced transport equation (20) — E_μ[∂tU] + u·E_μ[D] − ∇·u + c = 0, where D is the local mean force (17), E_μ denotes expectation under the time-dependent Boltzmann measure conditioned on the CV level set ξ(x)=y, and c approximates ∂t log Z(t) — is the mathematical engine. The training loss (26) writes the equation as a pointwise residual plus a mean-force regression, and estimates the required expectations by running the coupled non-equilibrium dynamics (8)–(9) and applying the Jarzynski reweighting identity (10). The Hessian-informed harmonic prior (28) is the companion object for many-particle systems: it restricts the prior to the non-trivial Hessian modes at a local minimum, making

Load-bearing premise

The load-bearing premise is that jointly minimizing the mean-force regression ℓ1 and the transport residual ℓ2 lands on the intersection where both are satisfied, so the learned velocity actually transports the prior to the target CV distribution; the paper states this without proof, and a spurious joint optimum would silently break the sampler.

What would settle it

On a low-dimensional system where the conditional expectations in (20) can be computed exactly by numerical quadrature on the CV level sets, train FES-FM and compare the learned u and v to the exact values. If, as λ grows, the residuals do not both vanish or the generated marginal fails to converge to the target ρ(y), the joint-minimization reasoning is false. A second check: feed the trained u into the Jarzynski reweighting estimator; the theory predicts the path weights become constant (zero variance), and a substantial residual weight spread would indicate the transport equation was not sol

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

Sample generation after training is a d-dimensional Euler integration, so cost per sample depends on the CV dimension, not the full configuration dimension n.
The free energy surface can be reconstructed directly by histogramming generated CV samples, without projecting or reweighting full-space trajectories.
The Jarzynski reweighting in training provides an unbiased estimate of the conditional expectations needed for the transport equation, and the method inherits variance reduction from the warm-up drift.
The E(3)-invariant Hessian prior should make the approach compatible with molecular systems whose potential has a well-defined local minimum and near-harmonic fluctuations.
On the tested potentials (Müller-Brown, double-well in 50–200 dimensions, two many-particle systems), the method matches or improves accuracy while cutting generation time by roughly an order of magnitude.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Because the reduced dynamics are formulated purely in terms of the CV map, the method could be combined with learned, data-driven collective variables to form an end-to-end free energy sampler; the paper does not explore this coupling.
The training stage still relies on NETS-style high-dimensional non-equilibrium simulation; removing that dependence, as the authors mention in Appendix E, would be the step that makes the method scale to realistic force fields.
A direct test of the joint-minimization claim: in a system with analytically known conditional mean force, the trained v should match it; if it does not, the ℓ1+λℓ2 objective is being trapped in a spurious minimum.
The coarea-formula derivation assumes ∇ξ has full rank everywhere in the region of interest; systems where the CV becomes degenerate (e.g., at a saddle) would need a regularized or smoothed CV map.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

A plausible reduced flow-matching sampler for free energy surfaces, but the core optimisation claim is unproven and the abstract overpromises.

read the letter

First, what this paper actually does: it takes the NETS Boltzmann sampler, derives the transport equation in collective-variable space, and constructs a training loss that estimates the conditional mean force via Jarzynski reweighting. That reduction is natural and the derivation up to equation (20) is clean. The Hessian-informed harmonic prior for many-particle systems is a genuinely useful idea, and the invariance proofs are solid. For someone working on generative samplers for free energy surfaces, this is a relevant and readable paper.

The problems are specific. The abstract advertises an alanine-dipeptide benchmark; no such experiment appears in the paper. That is a textual inconsistency that should have been caught before submission.

The bigger issue is §4.3. The paper asserts that jointly minimising ℓ1+λℓ2 over (u,c,v) yields a point satisfying both (22) and (24), and hence the reduced transport equation (20). That does not follow. It is true that in unrestricted function space the intersection of the two solution sets is non-empty. But the parameterised family is not unrestricted, the loss is nonconvex, and u and v enter multiplicatively in ℓ2. The global minimiser of the combined loss may trade off the two terms and satisfy neither equation pointwise. The paper gives no universal-approximation argument, no optimality guarantee, and no numerical check that the trained u actually solves (20). This matters because (20) is exactly what makes the reduced ODE transport the prior to the target distribution. The experiments show good histograms, but they do not close the gap.

The speed claim also needs care. FES-FM only avoids full-space simulation at generation time. Training still runs the full-space NETS dynamics to estimate the Jarzynski weights, and the paper admits this in the limitations. The abstract's 'drastically reduces computational costs' is misleading when only generation time is reported.

The math leading to the loss is sound, the loss design is sensible, and the experiments, while all toy systems, do show consistent accuracy. The Hessian prior alone is a citable contribution.

I would send this to peer review, but with a clear expectation of major revision. The authors need to fix the abstract, add the missing benchmark or remove it, prove or empirically verify the joint-minimiser claim, and report training time alongside generation time. Until then, treat the theoretical guarantee as conditional.

For a reading group, it is a useful discussion piece. I would cite the Hessian prior and the reduced-loss idea in related work.

Referee Report

4 major / 4 minor

Summary. FES-FM is a reduced flow-matching method for sampling the distribution ρ(y) of collective variables ξ(X) under a Boltzmann target p(x). The authors derive a reduced transport equation (20) for a CV-space velocity u, express the conditional force terms via the local mean force D(x,t), and define a combined loss (26) in which the expensive conditional expectations are reweighted by the NETS/Jarzynski identity (10). Training has a warm-up stage for the full-space drift bθ0 followed by main training of uθ1,cθ1,vθ1; generation solves the d-dimensional ODE (13). The paper also introduces a Hessian-informed harmonic prior with E(3)-invariance. Experiments on Müller-Brown, double-well potentials in n=50–200, and two many-particle systems report lower generation time than projecting full-space NETS samples (NETS-P), with comparable or better 1-Wasserstein error.

Significance. If the central transport claim were rigorously established, FES-FM would be a genuinely useful method: generation cost is decoupled from the ambient dimension, the loss reduction in (26) is elegant and avoids explicitly sampling the CV-conditioned manifolds, and the Hessian-informed prior addresses an important practical issue for molecular systems. The paper makes good use of standard identities (coarea formula, local mean force, Jarzynski reweighting), and the experiments, while limited to synthetic systems, provide initial evidence. The main weaknesses are that the mathematical justification of the training objective is incomplete and one advertised benchmark is absent.

major comments (4)

[§4.3, Eqs. (21)–(27)] The assertion that minimizing ℓ1+λℓ2 yields uθ1,cθ1,vθ1 satisfying both (22) and (24) is not established. This requires (i) existence of a point in the intersection of the two solution sets and (ii) attainment of zero loss by the chosen neural parameterization. The paper only states that the intersection is non-empty; it does not prove that the reduced Liouville equation (14) [equivalently (20)] admits a solution u, nor that the global minimizer of (26) attains the zero-loss point. Without this, Algorithm 2 is not guaranteed to transport the prior to ρ(y). Please add an existence argument (e.g., construct u from the conditional expectation of the full-space velocity b under μΣy,t) and either prove or empirically verify that the trained u satisfies (20) on held-out (y,t) points.
[Abstract and §6] The abstract states that the method is evaluated 'including alanine dipeptide in implicit solvent as a molecular benchmark,' but no alanine dipeptide experiment appears in Section 6 or the appendices. The experiments cover only Müller-Brown, high-dimensional double-well potentials, and two synthetic many-particle systems. This is a factual mismatch between the claimed evaluation and the reported results. The claim should be removed or the experiment added.
[§1 and Table 1] The paper claims that FES-FM 'drastically reduces computational costs,' but the 'Time' reported in Tables 1–2 is only the generation time (solving ODE (13)). Training (Algorithm 1) still requires full-space simulations of (8)–(9) in every epoch, and training wall-clock time is not reported. The abstract's unqualified cost claim is therefore misleading. Please qualify the claim as 'generation cost' and report end-to-end training plus generation time for a fair comparison.
[§5, Eq. (28)] The Hessian-informed prior U0(x) is defined through the alignment argmin R*0(x)=argmin_{R∈O(3)}∥(I⊗R)x−x0∥. Theorem 2 proves O(3)-invariance but not the regularity needed for the SDE (8), whose drift uses ∇U(x,t)=(1−t)∇U0(x)+t∇Utarget(x). The Kabsch alignment is generally only piecewise smooth, and at points where the argmin is non-unique U0 may be non-differentiable. The paper should either prove that this set has measure zero and the SDE is well-posed, or replace U0 by a smooth surrogate.

minor comments (4)

[§4.3, Eq. (27)] The notation ∇·u(ξ(x),t) should be explicitly defined as divergence with respect to the y argument of u; otherwise it can be confused with a full-space divergence.
[§2.1] The phrase 'generates samples for optimization' is unclear; please rephrase to indicate what is optimized.
[§6.1 / Fig. 1] The caption says 'the contour lines colored by the viridis colormap' but the figure is not colorblind-safe; consider using a labeled contour or a different colormap.
[Appendix D / Fig. D.5] The caption of Fig. D.5(a) describes FES-FM as using Npre training epochs; clarify that Npre is the number of warm-up epochs for bθ0, not for uθ1.

Circularity Check

0 steps flagged

No significant circularity; core reduced-transport derivation is self-contained, and same-author citations are non-load-bearing.

full rationale

The derivation chain is not circular. The target CV density ρ(y) is fixed by U and ξ through Eq. (1); the reduced transport equation (20) is obtained by projecting the full Liouville equation via the coarea formula and identities (16)-(18). The training loss (26) is the expectation of the squared residuals of (22) and (24) under νt, with Eνt estimated by the Jarzynski reweighting identity (10), which holds exactly for any drift b̂; no ground-truth CV samples or precomputed free-energy values are used as labels. Thus the learned u would be a genuine solution of (20), not a renamed fit to the target. The joint-minimizer claim in §4.3 is an unproven optimization/representability assumption (finite capacity, nonconvexity), which is a correctness risk, not a circular reduction: it does not presuppose the target distribution or the trained velocity. The only same-author citations (Liu et al. 2025b,c, Appendix C.3) are offered as alternative E(3)-equivariant architecture choices and are not load-bearing for the main derivation; no uniqueness theorem or ansatz is imported from the authors' prior work. Hence no circular step is exhibited, and the score is 0.

Axiom & Free-Parameter Ledger

4 free parameters · 7 axioms · 0 invented entities

The method's mathematics relies on standard identities (coarea formula, Liouville equation, Jarzynski equality) and on domain assumptions about the CV map and the potential. The distinctively paper-specific assumptions are the unproved joint-minimizer claim, the well-posedness of the reduced ODE and the Hessian prior, and the smoothness of the alignment-based prior. No new physical entities are introduced.

free parameters (4)

λ (loss weight in L1) = 1.0
Trade-off between mean-force regression and transport-equation residual; fixed for all experiments (Appendix C).
ϵ (diffusion coefficient in SDE (8)) = 0.2 (Müller-Brown), 0.1 (DW-50/100/200D), 0.02 (R2-3P, R3-4P)
Hand-chosen per dataset; ablation Figure D.5c shows performance is sensitive to it.
Npre (warm-up epochs) = 2000 for all experiments
Controls accuracy of bθ0 and hence variance of Jarzynski reweighting; ablation shows performance varies with it.
Prior hyperparameters for Müller-Brown = mean (-0.7, 0.7), std 0.3
Hand-chosen Gaussian prior; no sensitivity analysis provided.

axioms (7)

standard math Coarea formula and Liouville equation for ODE marginals (Eq. 1, Eq. 14)
Used to express the CV density and the reduced transport equation.
standard math Jarzynski equality / NETS reweighting identity (Eq. 10-11)
Taken from Albergo & Vanden-Eijnden 2025; used to estimate expectations under νt in the training loss.
domain assumption ξ ∈ C^1 and rank(∇ξ)=d in the region of interest (Section 1)
Required for the coarea formula and Moore-Penrose inverse expressions.
domain assumption U is twice continuously differentiable, O(3)-invariant, and translation-invariant; H(x0) has rank n-6 (Theorem 1)
Needed for the Hessian-informed prior and its null-space characterization.
ad hoc to paper The joint minimizer of ℓ1 + λℓ2 lies in the intersection of solution sets of (22) and (24) (Section 4.3)
Unproved optimization claim; if false, the trained velocity does not satisfy the transport equation.
ad hoc to paper The reduced ODE (13) with velocity satisfying (20) has a unique solution whose marginal is ρ(y,t)
Well-posedness of the learned flow; not discussed in the paper.
ad hoc to paper The Hessian-informed prior U0(x) of Eq. (28) is smooth enough for the SDE (8) despite the alignment argmin
The Kabsch-alignment argmin may be non-smooth or non-unique; the paper does not analyze regularity.

pith-pipeline@v1.3.0-alltime-deepseek · 16610 in / 15474 out tokens · 147477 ms · 2026-08-02T15:04:23.904064+00:00 · methodology

0 comments

read the original abstract

Sampling the distribution of collective variables (CVs) and estimating the associated free energy surface are crucial problems in statistical physics, as they underpin a better understanding of chemical reactions and conformational transitions. Traditional methods usually rely on simulations in high-dimensional configuration space and project the resulting configurations onto the CV space. To improve sampling speed, we propose FES-FM, a reduced flow matching (FM) method for free energy surface (FES) sampling. We train a dynamical transport map in the CV space, thereby enabling direct sampling of CV distributions and reconstruction of the corresponding free energy surface. For many-particle systems, we construct a prior distribution based on the Hessian at a local minimum of the potential, which ensures both rotation-translation invariance and physically meaningful configurations. We evaluate the proposed method across a variety of potential functions and collective variables, including alanine dipeptide in implicit solvent as a molecular benchmark. Comparative experiments demonstrate that our approach significantly improves sampling speed while maintaining accuracy.

Figures

Figures reproduced from arXiv: 2605.00337 by Tiejun Li, Zichen Liu.

**Figure 2.** Figure 2: Results of the many-particle systems. The red histograms in all subfigures denote the ground-truth [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

**Figure 3.** Figure 3: Comparison of sampling workflows between FES-FM and NETS-P. Starting from the prior [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗

**Figure 4.** Figure 4: Illustration of multi-particle systems. (a) A three-particle system in a 2D plane, where the CV is [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗

Review history (2 revisions) →

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

iSMART: An Iterative Sampling-and-Regression Technique for Solving Martingale-Based PDEs
math.NA 2026-07 conditional novelty 6.0

Solving martingale-based PDEs becomes a sequence of least-squares regressions on SDE path samples, avoiding nested expectations and adversarial training.

Reference graph

Works this paper leans on

15 extracted references · 11 linked inside Pith · cited by 1 Pith paper

[3]

arXiv preprint arXiv:2211.01364

An optimal control perspective on diffusion-based generative modeling. arXiv preprint arXiv:2211.01364 . Bortoli, V.D., Hutchinson, M.J., Wirnsberger, P., Doucet, A.,

Pith/arXiv arXiv
[4]

ArXiv abs/2402.08667

Target score matching. ArXiv abs/2402.08667. Carbone, D., Hua, M., Coste, S., Vanden-Eijnden, E.,

Pith/arXiv arXiv
[5]

arXiv preprint arXiv:2410.03282

Neural sampling from boltzmann densities: Fisher-rao curves in the wasserstein geometry. arXiv preprint arXiv:2410.03282 . Chen, J., Richter, L., Berner, J., Blessing, D., Neumann, G., Anandkumar, A.,

Pith/arXiv arXiv
[6]

arXiv preprint arXiv:2412.07081

Sequential controlled langevin diffusions. arXiv preprint arXiv:2412.07081 . Ciccotti, G., Lelievre, T., Vanden-Eijnden, E.,

Pith/arXiv arXiv
[7]

arXiv preprint arXiv:2209.03003

Flow straight and fast: Learning to generate and transfer data with rectified flow. arXiv preprint arXiv:2209.03003 . Liu, Z., Zhang, W., Li, T., 2025b. Improving the euclidean diffusion generation of manifold data by mitigating score function singularity, in: The Thirty-ninth Annual Conference on Neural Information Processing Systems. Liu, Z., Zhang, W.,...

Pith/arXiv arXiv
[9]

ArXiv abs/2402.06320

Particle denoising diffusion sampler. ArXiv abs/2402.06320. Plainer, M., Wu, H., Klein, L., Günnemann, S., Noé, F.,

Pith/arXiv arXiv
[11]

arXiv preprint arXiv:2407.07873

Dynamical measure transport and neural pde solvers for sampling. arXiv preprint arXiv:2407.07873 . Tian, Y., Panda, N., Lin, Y.T.,

Pith/arXiv arXiv
[12]

Xu, Y., Wang, Y., Luo, S., Gao, K., He, T., Liu, C., He, D.,

Iterated energy-based flow matching for sampling from boltzmann densities.arXiv:2408.16249. Xu, Y., Wang, Y., Luo, S., Gao, K., He, T., Liu, C., He, D.,

Pith/arXiv arXiv
[1955]

Mathematical Proceedings of the Cambridge Philosophical Society 51, 406–413

A generalized inverse for matrices. Mathematical Proceedings of the Cambridge Philosophical Society 51, 406–413. doi:10.1017/S0305004100030401. Phillips, A., Dau, H.D., Hutchinson, M.J., Bortoli, V.D., Deligiannidis, G., Doucet, A.,

doi:10.1017/s0305004100030401
[2018]

arXiv preprint arXiv:1809.10188

Monge-amp\ere flow for generative modeling. arXiv preprint arXiv:1809.10188 . Zhang, Q., Chen, Y.,

Pith/arXiv arXiv
[2021]

Forthethree-particlesystemin R2, theparametersin (31)arechosenas α1 = 5000/49, α2 = 5000/49, α3 = 50, r1 = 2 , r2 = 2 , r3 = 2 .4, r4 = 3 .1

and those based on alignment with respect to a given reference configuration (Liu et al., 2025b,c). Forthethree-particlesystemin R2, theparametersin (31)arechosenas α1 = 5000/49, α2 = 5000/49, α3 = 50, r1 = 2 , r2 = 2 , r3 = 2 .4, r4 = 3 .1. For the four-particle system in R3, the parameters in (32) are chosen as α1 = α2 = α3 = α4 = α5 = 5000 /49, α6 = 20...

2000
[2022]

arXiv preprint arXiv:2209.15571

Building normalizing flows with stochastic inter- polants. arXiv preprint arXiv:2209.15571 . Albergo, M.S., Vanden-Eijnden, E.,

Pith/arXiv arXiv
[2024]

ArXiv abs/2402.06121

Iterated denoising energy matching for sampling from boltzmann densities. ArXiv abs/2402.06121. Albergo, M.S., Vanden-Eijnden, E.,

Pith/arXiv arXiv
[2025]

ArXiv abs/2506.17139

Consistent sampling and simulation: Molecular dynamics with energy-based diffusion models. ArXiv abs/2506.17139. Raissi, M., Perdikaris, P., Karniadakis, G.E.,

arXiv
[2048]

The 1-Wasserstein distance is computed using 10,000 points sampled from the generated distribution and 10,000 points sampled from the ground-truth distribution

ODE(4) and ODE (13) are solved via the Euler method with a step size of 0.001. The 1-Wasserstein distance is computed using 10,000 points sampled from the generated distribution and 10,000 points sampled from the ground-truth distribution. Remaining hyperparameters are summarized in Table C.4. Detailed experimental setups are provided in the subsections b...

2000

[1] [3]

arXiv preprint arXiv:2211.01364

An optimal control perspective on diffusion-based generative modeling. arXiv preprint arXiv:2211.01364 . Bortoli, V.D., Hutchinson, M.J., Wirnsberger, P., Doucet, A.,

Pith/arXiv arXiv

[2] [4]

ArXiv abs/2402.08667

Target score matching. ArXiv abs/2402.08667. Carbone, D., Hua, M., Coste, S., Vanden-Eijnden, E.,

Pith/arXiv arXiv

[3] [5]

arXiv preprint arXiv:2410.03282

Neural sampling from boltzmann densities: Fisher-rao curves in the wasserstein geometry. arXiv preprint arXiv:2410.03282 . Chen, J., Richter, L., Berner, J., Blessing, D., Neumann, G., Anandkumar, A.,

Pith/arXiv arXiv

[4] [6]

arXiv preprint arXiv:2412.07081

Sequential controlled langevin diffusions. arXiv preprint arXiv:2412.07081 . Ciccotti, G., Lelievre, T., Vanden-Eijnden, E.,

Pith/arXiv arXiv

[5] [7]

arXiv preprint arXiv:2209.03003

Flow straight and fast: Learning to generate and transfer data with rectified flow. arXiv preprint arXiv:2209.03003 . Liu, Z., Zhang, W., Li, T., 2025b. Improving the euclidean diffusion generation of manifold data by mitigating score function singularity, in: The Thirty-ninth Annual Conference on Neural Information Processing Systems. Liu, Z., Zhang, W.,...

Pith/arXiv arXiv

[6] [9]

ArXiv abs/2402.06320

Particle denoising diffusion sampler. ArXiv abs/2402.06320. Plainer, M., Wu, H., Klein, L., Günnemann, S., Noé, F.,

Pith/arXiv arXiv

[7] [11]

arXiv preprint arXiv:2407.07873

Dynamical measure transport and neural pde solvers for sampling. arXiv preprint arXiv:2407.07873 . Tian, Y., Panda, N., Lin, Y.T.,

Pith/arXiv arXiv

[8] [12]

Xu, Y., Wang, Y., Luo, S., Gao, K., He, T., Liu, C., He, D.,

Iterated energy-based flow matching for sampling from boltzmann densities.arXiv:2408.16249. Xu, Y., Wang, Y., Luo, S., Gao, K., He, T., Liu, C., He, D.,

Pith/arXiv arXiv

[9] [1955]

Mathematical Proceedings of the Cambridge Philosophical Society 51, 406–413

A generalized inverse for matrices. Mathematical Proceedings of the Cambridge Philosophical Society 51, 406–413. doi:10.1017/S0305004100030401. Phillips, A., Dau, H.D., Hutchinson, M.J., Bortoli, V.D., Deligiannidis, G., Doucet, A.,

doi:10.1017/s0305004100030401

[10] [2018]

arXiv preprint arXiv:1809.10188

Monge-amp\ere flow for generative modeling. arXiv preprint arXiv:1809.10188 . Zhang, Q., Chen, Y.,

Pith/arXiv arXiv

[11] [2021]

Forthethree-particlesystemin R2, theparametersin (31)arechosenas α1 = 5000/49, α2 = 5000/49, α3 = 50, r1 = 2 , r2 = 2 , r3 = 2 .4, r4 = 3 .1

and those based on alignment with respect to a given reference configuration (Liu et al., 2025b,c). Forthethree-particlesystemin R2, theparametersin (31)arechosenas α1 = 5000/49, α2 = 5000/49, α3 = 50, r1 = 2 , r2 = 2 , r3 = 2 .4, r4 = 3 .1. For the four-particle system in R3, the parameters in (32) are chosen as α1 = α2 = α3 = α4 = α5 = 5000 /49, α6 = 20...

2000

[12] [2022]

arXiv preprint arXiv:2209.15571

Building normalizing flows with stochastic inter- polants. arXiv preprint arXiv:2209.15571 . Albergo, M.S., Vanden-Eijnden, E.,

Pith/arXiv arXiv

[13] [2024]

ArXiv abs/2402.06121

Iterated denoising energy matching for sampling from boltzmann densities. ArXiv abs/2402.06121. Albergo, M.S., Vanden-Eijnden, E.,

Pith/arXiv arXiv

[14] [2025]

ArXiv abs/2506.17139

Consistent sampling and simulation: Molecular dynamics with energy-based diffusion models. ArXiv abs/2506.17139. Raissi, M., Perdikaris, P., Karniadakis, G.E.,

arXiv

[15] [2048]

The 1-Wasserstein distance is computed using 10,000 points sampled from the generated distribution and 10,000 points sampled from the ground-truth distribution

ODE(4) and ODE (13) are solved via the Euler method with a step size of 0.001. The 1-Wasserstein distance is computed using 10,000 points sampled from the generated distribution and 10,000 points sampled from the ground-truth distribution. Remaining hyperparameters are summarized in Table C.4. Detailed experimental setups are provided in the subsections b...

2000