REVIEW 3 major objections 5 minor 15 references
Thermodynamic Interpolation: A generative approach to molecular thermodynamics and kinetics
T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Thermodynamic Interpolation learns a differentiable map between Boltzmann ensembles at different temperatures, and the paper claims this map generates accurate ensembles, free energies, and approximate kinetic rates at temperatures not…
desk verdict A promising temperature-interpolating flow method, but the printed importance weight in Eq. (20) has both energy terms negative, which breaks the unbiased sampling construction as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the learned diffeomorphic map $f_{AB}$ between Boltzmann distributions $\mu_A$ and $\mu_B$ at temperatures $T_A$ and $T_B$, implemented as a continuous normalizing flow with an SE(3)-equivariant ChiroPaiNN velocity field $b^{(\theta)}(t,x;T_A,T_B)$ trained by the stochastic interpolant objective. The two-sided interpolant $x(t) = (1-t)x_A + t x_B + \gamma(t) z$ defines a stochastic path between paired samples, and minimizing the regression loss (eq. 19) yields a velocity whose integrated ODE transports $x_A \sim \mu_A$ to $x_B \sim \rho_B$. The same flow supplies the log-Jacobian determinant entering the importance weights $w_a$ (eq. 20), which correct the surrogate samples to unbiased estimates at the target temperature. For kinetics, the machinery is the gEDMD estimator of the Koopman generator for overdamped Langevin dynamics, whose eigenvalues give relaxation rates, applied to the TI-generated equilibrium samples.
What would settle it
Run long unbiased underdamped molecular dynamics for N-methylformamide at 300 K, build a Markov state model, and compare its slowest relaxation rate with the gEDMD rate obtained from TI-generated samples at the same temperature; a substantial mismatch would localize the failure in the Brownian-generator assumption rather than the learned map.
Extended reading notes
Core claim
The central claim is that a continuous normalizing flow can learn a velocity field, trained by stochastic interpolants, that maps a Boltzmann ensemble at one temperature directly onto another, and that this map generalizes to target temperatures outside the training set. On the model systems studied (a one-dimensional asymmetric double well and two molecules from the QM9 set, N-methylformamide and 3-propan-2-ylhex-1-yne), the ambient TI map produces torsion-angle histograms and TICA projections that match reference low-temperature MD, and Kish effective sample sizes remain high enough that free-energy differences from as few as five transformed samples match what direct reweighting needs roughly twenty times more samples to achieve. The latent TI variant, a temperature-conditioned Boltzmann generator, samples multiple temperatures on demand from a shared latent Gaussian. The authors also show that gEDMD, applied to TI-generated equilibrium samples, yields rates for slow metastable exchanges whose temperature trends agree with rates estimated from MD, provided the underlying dynamics are taken to be overdamped Langevin (Brownian) dynamics.
Load-bearing premise
The kinetic rate estimates assume the molecular dynamics follows overdamped Langevin (Brownian) dynamics, so the generator used in gEDMD is an approximation for real, underdamped molecules; if that assumption fails, the predicted rates are not the true transition rates.
Editorial extensions
If this is right
- Within the trained temperature range, one aTI model can produce ensembles at arbitrary target temperatures, replacing multiple dedicated low-temperature simulations.
- Free-energy differences between temperatures can be estimated from very few transformed samples, which is valuable when sampling is expensive.
- TI-generated equilibrium samples feed directly into gEDMD, giving temperature-dependent kinetic rates without time-correlated trajectories.
- The latent TI variant shares statistical information across temperatures, so a single model can generate on-demand samples at several temperatures with improved accuracy compared with separately trained models.
Reading between the lines
- The same learned-map architecture should extend to thermodynamic variables other than temperature, such as pressure or chemical potential, since the flow is conditioned on the state parameters; this is a natural but untested extrapolation of the paper's approach.
- Because the paper finds gEDMD rates are robust to mapping bias even without reweighting, kinetic studies might skip importance weights entirely, accepting a small systematic error for large computational savings; the paper hints at this but does not recommend it.
- The high sample efficiency at very small sample counts suggests TI could be applied to expensive potential-energy surfaces, such as quantum-mechanical calculations, where generating many MD frames is prohibitive.
- Multi-temperature training may reduce the amount of simulation data needed per temperature, mirroring the information sharing seen in multi-ensemble Markov models; the paper observes this effect but does not quantify data requirements.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript introduces Thermodynamic Interpolation (TI), a generative method for sampling molecular configurations across thermodynamic states. Two variants are presented: ambient TI, which learns a direct diffeomorphic map between Boltzmann ensembles at different temperatures in configuration space, and latent TI, which maps from a normal latent distribution to a temperature-conditioned ensemble. The velocity fields are trained with stochastic interpolants and SE(3)-equivariant message passing networks. The authors validate TI on a 1D asymmetric double well and two small molecules (N-methylformamide and 3-propan-2-ylhex-1-yne), reporting high effective sample sizes, accurate free energy differences via targeted free energy perturbation, and temperature-dependent kinetic rates obtained from generator EDMD. The paper claims interpolation and extrapolation in temperature, including temperatures not seen during training.
Significance. If fully correct, the TI framework would be a useful addition to the Boltzmann generator toolbox: it shares information across temperatures, avoids replica ladders, and connects generative sampling to TFEP and gEDMD. The paper's strengths are its grounding in the stochastic interpolant framework, the use of standard TFEP relations, and held-out temperature validation against independent MD simulations and numerical integration. Experimental comparisons are plausible, and the authors are transparent about some limitations (Jacobian cost, Brownian generator assumption). However, the central importance-weight formula for ambient TI contains a sign error that invalidates the unbiasedness argument as written, so the quantitative claims must be re-evaluated. The kinetics section is also explicitly model-dependent and should be framed as Brownian dynamics rates rather than true molecular rates.
major comments (3)
- [Methods, Ambient TI, Eq. (20)] The importance weight in Eq. (20) is printed with two negative energy terms. A change-of-variables derivation from xA ~ mu_A to y = f_AB(xA) gives a target/source density ratio proportional to exp(-E(y)/(kTB) + E(xA)/(kTA)) times |det J|, so the first energy term must be positive. In the identity limit TA=TB and f_AB=id, the correct weight is 1 up to the usual partition-function constant, whereas Eq. (20) yields exp(-2E(x)/(kT)), which after normalization preferentially weights low-energy conformations and is not the target Boltzmann measure. Since Eq. (20) is the basis for the reported aTI ESS values and for unbiased observable estimates from ambient TI, the central sampling claim is not supported by the manuscript as written. The formula must be corrected and the numerical results re-examined.
- [Methods, gEDMD; Results, Generator Estimation] The kinetic-rate estimates are derived under the overdamped Langevin generator in Eqs. (24)-(25), and the text states "assuming the generator is Brownian." For the molecular systems, the reference MD is underdamped, so the estimated rates are not true molecular transition rates; the paper's own MSM comparison for the ADW system shows that overdamped rates bound underdamped rates from below. The claims about "kinetic exchange rates" for molecules should be explicitly restricted to the Brownian model, and the comparison to MD-derived rates obtained with the same misspecified estimator does not validate true kinetics. This is a caveated but load-bearing limitation for the kinetics section.
- [Results, Figures 3 and 4; IQR filtering] The authors state that they use an IQR-filtering strategy to "avoid high variance on the free energy estimators" (Figures 3B-E and 4B-C). If outliers are removed based on the values of the estimator or the importance weights, this can bias the reported free energies and effective sample sizes. The Supporting Information must specify the filtering rule, clarify whether it depends on the estimator, and demonstrate (e.g., by reporting unfiltered or alternatively robust estimates) that the reported results are not an artifact of the filter.
minor comments (5)
- [Abstract and Introduction] The term "Thermodynamic Interpolation" is introduced twice, and the sentence "We introduce TI... We introduce TI flavors" is redundant; please tighten the wording.
- [Results, first paragraph] The molecule is referred to as both "3p2y1y" (abstract and Figure 5) and "3p2y1p" in the main text; please make the abbreviation consistent.
- [Methods, Eq. (21)] The normalized temperature \hat T is defined with an overline T that is not explicitly defined; please define all symbols in the equation.
- [Conclusion] The statement that exact computation of free energies "is not tractable in the current architectures" for larger systems is a strong limitation; please quantify the scaling for the tested systems or add a reference to the stated linear Jacobian cost.
- [Manuscript formatting] The "TOC Graphic" appears twice on the final pages, which appears to be a duplication in the manuscript.
Circularity Check
No significant circularity: the held-out-temperature predictions are genuine external benchmarks, and the cited methods are independent tools rather than self-referential premises.
full rationale
The paper's central claims are the learned ambient/latent TI maps can generate samples at temperatures held out from training, that free energy differences can be estimated from those maps, and that kinetics can be approximated from equilibrium samples. These claims are validated against independent references: for the asymmetric double well, free energies are compared to numerical integration of the Boltzmann distribution, and for the molecular systems, generated histograms and effective sample sizes are compared to reference MD simulations. The aTI models are explicitly trained with the target temperature excluded from the training set, so the reported results are genuine predictions rather than fitted quantities. The importance weights in Eq. (20) define the unbiased estimator, but the empirical ESS values and free energy estimates are benchmarked externally rather than being forced by construction. The paper does contain a separate internal reproducibility issue: Eq. (20) as printed has a sign error in the source-energy term, so at TA = TB with the identity map it yields weights proportional to exp(-2E/kT) instead of unity; however, this is a correctness/reproducibility defect, not a circularity in the derivation chain. Self-citations appear for the ChiroPaiNN architecture and for gEDMD/RFF estimators, but these are methodological tools from prior work by overlapping authors, not uniqueness theorems or ansatz assumptions used to force the present conclusions; they are externally established and machine-implementable. The overdamped Langevin generator assumption for kinetics is explicitly stated and caveated by the authors, and rates from underdamped simulations are also discussed. No step in the paper reduces, by definition or by fitted input, to the quantity it claims to predict. The derivation chain is therefore self-contained against external benchmarks, and the appropriate circularity score is 0.
Assumptions & free parameters
free parameters (4)
- l0
- gamma(t) noise schedule
- gEDMD RFF bandwidth sigma and feature count p =
sigma=0.6, p=50
- IQR filter threshold
assumptions (5)
- standard math The canonical Boltzmann distribution mu(x) = Z^{-1} exp(-E(x)/(kT)) is the target distribution for each thermodynamic state.
- domain assumption There exists a smooth invertible map (diffeomorphism) between Boltzmann distributions at different temperatures that the CNF can approximate.
- domain assumption Molecular dynamics can be approximated as overdamped Langevin dynamics (eq 24) so that the gEDMD generator (eq 25) yields kinetic rates.
- domain assumption The learned models are trained on simulation data from a subset of temperatures and generalize to held-out temperatures.
- domain assumption SE(3)-equivariance is sufficient for molecular systems because chirality is conserved in classical MD data.
Cite this review
Pith. "Pith review of Thermodynamic Interpolation: A generative approach to molecular thermodynamics and kinetics." pith.science (2026). https://pith.science/paper/6G2YL7W3
@misc{pith2026241110075,
author = {Pith},
title = {Pith review of: Thermodynamic Interpolation: A generative approach to molecular thermodynamics and kinetics},
year = {2026},
howpublished = {\url{https://pith.science/paper/6G2YL7W3}},
note = {Machine review of arXiv:2411.10075}
}
read the original abstract
Using normalizing flows and reweighting, Boltzmann Generators enable equilibrium sampling from a Boltzmann distribution, defined by an energy function and thermodynamic state. In this work, we introduce Thermodynamic Interpolation (TI), which allows for generating sampling statistics in a temperature-controllable way. We introduce TI flavors that work directly in the ambient configurational space, mapping between different thermodynamic states or through a latent, normally distributed reference state. Our ambient-space approach allows for the specification of arbitrary target temperatures, ensuring generalizability within the temperature range of the training set and demonstrating the potential for extrapolation beyond it. We validate the effectiveness of TI on model systems that exhibit metastability and non-trivial temperature dependencies. Finally, we demonstrate how to combine TI-based sampling to estimate free energy differences through various free energy perturbation methods and provide corresponding approximated kinetic rates estimated through generator extended dynamic mode decomposition (gEDMD).
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
The Journal of Physical Chemistry B 2004, 108, 6571–6581
Theory. The Journal of Physical Chemistry B 2004, 108, 6571–6581. (20) Buchete, N.-V.; Hummer, G. Coarse Master Equations for Peptide Folding Dynamics. The Journal of Physical Chemistry B 2008, 112, 6057–6069. (21) Husic, B. E.; Pande, V. S. Markov State Models: From an Art to a Science. Journal of the American Chemical Society 2018, 140, 2386–2396. (22) ...
work page 2004
-
[6]
Boltzmann generators: Sampling equilibrium states of many-body systems with deep learning
30 (45) No´ e, F.; Olsson, S.; K¨ ohler, J.; Wu, H. Boltzmann generators: Sampling equilibrium states of many-body systems with deep learning. Science 2019, 365, eaaw1147. (46) Rezende, D. J.; Mohamed, S. Variational Inference with Normalizing Flows
work page 2019
-
[10]
Neural Thermodynamic Integration: Free Energies from Energy-based Diffusion Models
(54) Rizzi, A.; Carloni, P.; Parrinello, M. Targeted Free Energy Perturbation Revisited: Accurate Free Energies from Mapped Reference Potentials. The Journal of Physical Chemistry Letters 2021, 12, 9449–9454, PMID: 34555284. 31 (55) Wirnsberger, P.; Ballard, A. J.; Papamakarios, G.; Abercrombie, S.; Racani` ere, S.; Pritzel, A.; Jimenez Rezende, D.; Blund...
work page Pith review arXiv 2021
-
[14]
(72) Williams, M. O.; Kevrekidis, I. G.; Rowley, C. W. A Data–Driven Approximation of the Koopman Operator: Extending Dynamic Mode Decomposition. Journal of Nonlinear Science 2015, 25, 1307–1346. (73) Klus, S.; N¨ uske, F.; Koltai, P.; Wu, H.; Kevrekidis, I.; Sch¨ utte, C.; No´ e, F. Data-Driven Model Reduction and Transfer Operator Approximation. Journal...
work page 2015
-
[233]
(48) Diez, J. V.; Atance, S. R.; Engkvist, O.; Olsson, S. Generation of conformational en- sembles of small molecules via surrogate model-assisted molecular dynamics. Machine Learning: Science and Technology 2024, 5, 025010. (49) K¨ ohler, J.; Klein, L.; No´ e, F. Equivariant flows: exact likelihood generative learning for symmetric densities. Internation...
arXiv 2024
-
[342]
Variational characterization of free energy: Theory and algorithms
(65) Hartmann, C.; Richter, L.; Sch¨ utte, C.; Zhang, W. Variational characterization of free energy: Theory and algorithms. Entropy 2017, 19,
work page 2017
-
[626]
Density estimation using Real NVP
32 (66) Dinh, L.; Sohl-Dickstein, J.; Bengio, S. Density estimation using Real NVP. 2017; https://arxiv.org/abs/1605.08803. (67) Chen, R. T.; Rubanova, Y.; Bettencourt, J.; Duvenaud, D. K. Neural ordinary differ- ential equations. Advances in neural information processing systems 2018, 31 . (68) Thomas, N.; Smidt, T.; Kearnes, S.; Yang, L.; Li, L.; Kohlho...
arXiv 2017
-
[1965]
(75) Molgedey, L.; Schuster, H. G. Separation of a mixture of independent signals using time delayed correlations. Phys. Rev. Lett. 1994, 72, 3634–3637. 33 (76) Ziehe, A.; M¨ uller, K.-R.ICANN 98 ; Springer London, 1998; p 675–680. (77) P´ erez-Hern´ andez, G.; Paul, F.; Giorgino, T.; De Fabritiis, G.; No´ e, F. Identification of slow molecular order para...
work page 1994
Show all 15 references
-
[1998]
D.; 27 Sch¨ utte, C.; No´ e, F
(18) Prinz, J.-H.; Wu, H.; Sarich, M.; Keller, B.; Senne, M.; Held, M.; Chodera, J. D.; 27 Sch¨ utte, C.; No´ e, F. Markov models of molecular kinetics: Generation and validation. The Journal of Chemical Physics 2011, 134, 174105. (19) Swope, W. C.; Pitera, J. W.; Suits, F. De...
2011
-
[2002]
M.; Swendsen, R
(3) Ferrenberg, A. M.; Swendsen, R. H. Optimized Monte Carlo data analysis. Physical Review Letters 1989, 63, 1195–1198. (4) Shirts, M. R.; Chodera, J. D. Statistically optimal analysis of samples from multiple equilibrium states. The Journal of Chemical Physics 2008, 129 . (5...
1989
-
[2016]
G.; Vanden-Eijnden, E
(47) Tabak, E. G.; Vanden-Eijnden, E. DENSITY ESTIMATION BY DUAL ASCENT OF THE LOG-LIKELIHOOD. Communications in Mathematical Sciences 2010, 8, 217–
2010
-
[2021]
Implicit Transfer Operator Learning: Multi- ple Time-Resolution Models for Molecular Dynamics
(70) Schreiner, M.; Winther, O.; Olsson, S. Implicit Transfer Operator Learning: Multi- ple Time-Resolution Models for Molecular Dynamics. Advances in Neural Information Processing Systems. 2023; pp 36449–36462. (71) Schebek, M.; Invernizzi, M.; No´ e, F.; Rogal, J. Efficient ...
2023
-
[2022]
(52) Zheng, S. et al. Predicting equilibrium distributions for molecular systems with deep learning. Nature Machine Intelligence 2024, 6, 558–567. (53) Herron, L.; Mondal, K.; Schneekloth, J. S.; Tiwary, P. Inferring phase transitions and critical exponents from limited observ...
2024
-
[2023]
V AMPnets for deep learning of molecular kinetics
(31) Mardt, A.; Pasquali, L.; Wu, H.; No´ e, F. V AMPnets for deep learning of molecular kinetics. Nature Communications 2018, 9 . (32) Chen, W.; Sidky, H.; Ferguson, A. L. Nonlinear discovery of slow molecular modes using state-free reversible V AMPnets.The Journal of Chemica...
2018 arXiv
-
[2024]
(41) Klein, L.; Foong, A. Y. K.; Fjelde, T. E.; Mlodozeniec, B.; Brockschmidt, M.; Nowozin, S.; No´ e, F.; Tomioka, R. Timewarp: Transferable Acceleration of Molec- ular Dynamics by Learning Time-Coarsened Dynamics. 2023; https://arxiv.org/ abs/2302.01170. (42) Hsu, T.; Sadigh...
2023 arXiv
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.