A Residence-Time Approach for Determining Position-Dependent Diffusivities from Biased Molecular Simulations

Michael von Domaros; Praveen Ranganath Prabhakar; Rinto Thomas

A Residence-Time Approach for Determining Position-Dependent Diffusivities from Biased Molecular Simulations

T0 review · 3 major / 2 minor · reviewed 2026-07-13 · grok-4.5

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

Pith's one-line read Residence-time statistics extract local diffusivities from biased MD after free-energy gradient compensation.

desk verdict Useful methods paper for D(z) from biased MD; abstract-only so confidence is low, but the checks and the lag-time caveat look honest. read the letter →

arxiv 2604.01940 v3 pith:DS2PIBAY submitted 2026-04-02 cond-mat.soft physics.bio-phphysics.chem-ph

Rinto Thomas , Praveen Ranganath Prabhakar , Michael von Domaros This is my paper · ORCID

classification cond-mat.softphysics.bio-phphysics.chem-ph

keywords position-dependentdiffusivityresidence-timeapproachbiasedmolecularsimulationspotentialofmeanforcemembranepermeationfirst-exitstatisticspropagatorreconstruction

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Position-dependent diffusivities are the key missing parameters when one reduces molecular transport in heterogeneous environments to a simple one-dimensional stochastic description. This paper shows that those local diffusivities can be recovered from first-exit (residence-time) statistics measured in biased simulations, once the mean free-energy gradient has been compensated. Applied to oxygen crossing a hexadecane/water slab, water crossing a POPC bilayer, and water plus volatiles crossing a model skin barrier, the method yields PMF–diffusivity pairs whose propagators match unbiased molecular-dynamics propagators over substantial lag-time windows and recover known bulk values in the slab. The same tests also expose cases in which no single lag-time-independent diffusivity profile can describe every timescale, clarifying both the reach and the limit of the effective one-dimensional picture.

What carries the argument

The residence-time approach (RTA): first-exit times measured inside spatially localized windows after free-energy-gradient compensation are inverted for a local diffusivity that is then paired with the PMF to form a reduced stochastic model.

What would settle it

In a heterogeneous system where independent bulk or lag-time-resolved diffusivities are known, check whether the RTA profile recovers those references and whether the PMF–diffusivity pair’s propagator matches unbiased MD propagators at every lag time of interest; systematic mismatch at long or short lags falsifies the claim of a sufficient lag-independent profile.

Watch

Extended reading notes

Core claim

After the mean free-energy gradient is removed, first-exit statistics collected in biased simulations furnish effective position-dependent diffusivities that, together with the potential of mean force, reconstruct the propagator dynamics of the corresponding unbiased system over appreciable lag-time ranges and recover independent bulk reference values.

Load-bearing premise

After the mean free-energy gradient is compensated, a single lag-time-independent local diffusivity extracted from first-exit statistics is enough to reconstruct the full propagator of the unbiased dynamics.

Editorial extensions

If this is right

PMF–diffusivity pairs obtained by RTA can be used to build reduced stochastic models whose short-to-intermediate lag propagators match unbiased MD.
Bulk reference diffusivities in slab geometries can be recovered directly from biased first-exit statistics.
When a single lag-independent profile fails, the same diagnostics flag the breakdown of the one-dimensional Markovian description.
The approach supplies a practical route for transport coefficients in membranes and other soft heterogeneous media without requiring long unbiased trajectories.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

RTA could be combined with enhanced-sampling free-energy methods to map diffusivity profiles across entire reaction coordinates rather than only one-dimensional membranes.
Cases where no lag-independent profile works point to memory or orthogonal degrees of freedom that may require multi-dimensional or non-Markovian extensions of the same residence-time idea.
The method’s success in recovering bulk values suggests it can serve as an internal consistency check when experimental or independent computational diffusivities are available.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

This is a methods paper that claims a practical residence-time route to position-dependent diffusivities from biased MD: compensate the mean free-energy gradient, then pull local D(z) from first-exit statistics. The punchline is that the resulting PMF–D pairs match bulk references in a hexadecane/water slab and rebuild unbiased MD propagators over useful lag times in membrane systems, while the authors themselves note that a single lag-time-independent D(z) sometimes fails.

What looks new is the packaging: a bias-compatible first-exit estimator aimed at the soft-matter / biophysical workflow people already run (umbrella sampling and friends). Related free-energy and diffusivity estimation from biased trajectories exists, so this is incremental rather than a new theory of diffusion. Credit where due: the validation targets are the right ones—external bulk numbers and independent unbiased propagators—not self-consistency theater. The abstract’s own caveat that no single D(z) always captures all timescales is a real strength; it delimits the modeling premise instead of hiding it.

Soft spots, in proportion: we only have the abstract. No equations, no binning protocol, no lag-time dependence plots, no code or data. Free parameters (region width, lag windows, bias/FE hyperparameters) are inevitable and uninspectable here. Circularity risk is moderate-low because the free-energy profile used for gradient compensation usually comes from the same biased runs, but the external propagator and bulk checks reduce that. The load-bearing assumption is exactly what the authors flag: after mean-gradient compensation, first-exit stats give an effective D(z) that is lag-independent enough to reconstruct the propagator. When that fails, the method still reports useful effective profiles; it does not magically fix non-Markovian or multi-timescale physics.

Who it is for: people who build reduced stochastic models of membrane and interface transport and need D(z) from biased trajectories. Worth a serious referee if the full derivation, lag-time analysis, and numerical details hold up. I would send it to peer review rather than desk-reject; I would not cite from the abstract alone. Bring to reading group only if someone is actively estimating position-dependent diffusivities this quarter.

Referee Report

3 major / 2 minor

Summary. The manuscript proposes a residence-time approach (RTA) that extracts position-dependent diffusivities from first-exit statistics measured in biased molecular simulations after compensation of the mean free-energy gradient. Applications include oxygen diffusion across a hexadecane/water slab, water permeation across a POPC bilayer, and transport of water and volatile organic compounds through a model skin-barrier membrane. In the slab, RTA diffusivities are reported to agree with independent bulk references; in the membrane systems, PMF–diffusivity pairs constructed via RTA are reported to reproduce unbiased MD propagators over substantial lag-time ranges. The abstract also states that, in some cases, no single lag-time-independent diffusivity profile captures the dynamics across all timescales.

Significance. Position-dependent diffusivities are central inputs to reduced stochastic models of molecular transport in heterogeneous soft-matter environments. A practical estimator that can be applied to biased simulations and that is externally validated against bulk references and unbiased propagators would be useful for membrane permeation and related problems. The abstract’s external checks (bulk slab references; reconstruction of unbiased MD propagators) are the right kind of validation, and the explicit acknowledgment that lag-time-independent D(z) is not always sufficient is a presentational strength if the domain of validity is carefully mapped. Assessment of whether those strengths are realized is limited here to the abstract alone.

major comments (3)

Abstract: The load-bearing modeling premise is that, after mean free-energy-gradient compensation, first-exit (residence-time) statistics yield an effective position-dependent diffusivity that is lag-time-independent and, paired with the PMF, sufficient to reconstruct unbiased propagator dynamics. The abstract itself reports cases where no single lag-time-independent profile captures all timescales. That caveat must be turned into a concrete delimitation of when the Markovian RTA description holds and when multi-timescale or non-Markovian effects require a more general treatment; without that delimitation the practical claim remains under-specified.
Abstract (validation claims): Agreement with bulk references in the hexadecane/water slab and reconstruction of unbiased MD propagators in the membrane systems are appropriate external checks, but the abstract gives no quantitative metrics (lag-time windows of successful reconstruction, error measures, bin-width or region-width sensitivity). Those quantities are essential to judge whether the reported agreement is robust or confined to particular observation windows.
Abstract (method statement): The procedure depends on free parameters that are not quantified in the abstract—spatial bin/region width for first-exit statistics, lag-time/observation windows for propagator comparison, and bias and free-energy estimation hyperparameters. Systematic sensitivity of the reported D(z) profiles and of propagator agreement to these choices is load-bearing for the claim that RTA is a reliable practical route; it must be documented.

minor comments (2)

Abstract: The phrase “substantial lag-time ranges” is vague; specific numerical ranges (and the systems to which they apply) should be stated so that the scope of the propagator reconstruction claim is clear.
Abstract: The precise operational definition of the compensated dynamics and of the residence/first-exit time (including how the mean free-energy gradient is estimated and removed) should be introduced early and used consistently; the abstract alone leaves that definition implicit.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; abstract-only method paper with external validation against bulk references and independent unbiased MD.

full rationale

Only the abstract is available, so no equations, self-citations, uniqueness theorems, or fitted-parameter-to-prediction reductions can be exhibited. The claimed method extracts effective D(z) from first-exit (residence-time) statistics after free-energy-gradient compensation in biased simulations, then validates by (i) matching independently determined bulk reference diffusivities in a hexadecane/water slab and (ii) reconstructing unbiased MD propagators over lag-time ranges in membrane systems. Those benchmarks are external to the estimator itself; agreement is an empirical test, not a definitional identity. The abstract itself flags the modeling premise's limits (no single lag-time-independent D(z) always works), which is an honesty caveat rather than circularity. Absent any quoteable self-definitional step, fitted-input-called-prediction, load-bearing self-citation chain, or renaming of a known result, the circularity score is 0. The reader's moderate-low residual risk (PMF from same trajectories) is a shared-data correlation concern, not a by-construction reduction under the enumerated patterns.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

Abstract-only: free parameters and axioms are inferred from standard practice for this class of method. Central claim rests on overdamped 1D diffusion with position-dependent D after mean free-energy compensation, plus numerical choices (bins, lag times, bias) that are not quantified here. No new physical entities are introduced.

free parameters (3)

spatial bin / region width for first-exit statistics
Local residence-time estimators require a finite spatial window; width is a numerical choice that affects estimated D and is not fixed by the abstract.
lag-time / observation window for propagator comparison
Abstract reports agreement over 'substantial lag-time ranges' and failures of lag-time-independent D; the chosen windows are free analysis parameters.
bias and free-energy estimation hyperparameters
Biased simulations and PMF reconstruction involve force constants, collective-variable definitions, and reweighting choices that enter the gradient compensation step.

assumptions (3)

domain assumption After compensating the mean free-energy gradient, short-time first-exit statistics are governed by a local diffusivity D(z) in an overdamped 1D diffusion picture.
Core modeling premise of RTA; abstract applies it and notes cases where no single lag-time-independent D(z) works.
domain assumption Unbiased MD propagators are the external ground truth for validating reduced PMF–D models.
Used as the membrane-system validation standard in the abstract.
domain assumption Bulk reference diffusivities in homogeneous phases are reliable independent benchmarks for the slab system.
Abstract claims agreement with independently determined bulk values.

how reviews work

0 comments

Cite this review

Pith. "Pith review of A Residence-Time Approach for Determining Position-Dependent Diffusivities from Biased Molecular Simulations." pith.science (2026). https://pith.science/paper/DS2PIBAY

@misc{pith2026260401940,
  author       = {Pith},
  title        = {Pith review of: A Residence-Time Approach for Determining Position-Dependent Diffusivities from Biased Molecular Simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/DS2PIBAY}},
  note         = {Machine review of arXiv:2604.01940}
}

read the original abstract

Position-dependent diffusivities are central parameters in reduced stochastic descriptions of molecular transport in heterogeneous environments, but their reliable estimation from molecular dynamics simulations remains challenging. We present a residence-time approach (RTA) that extracts local diffusivities from first-exit statistics measured in biased simulations after compensation of the mean free-energy gradient. We apply the method to oxygen diffusion across a hexadecane/water slab, water permeation across a POPC lipid bilayer, and transport of water and volatile organic compounds through a model skin-barrier membrane. In the slab system, RTA diffusivities agree with independently determined bulk reference values. In the membrane systems, propagator predictions based on RTA-derived diffusivities reproduce unbiased molecular dynamics propagators over substantial lag-time ranges, while also revealing that, in some cases, no single lag-time-independent diffusivity profile captures the dynamics across all timescales. These results support residence-time statistics as a practical route for determining effective position-dependent diffusivities from biased molecular simulations.

Discussion (0). Continue with ORCID to comment.

Pith tools

Reviewed July 13, 2026 · model on record in the stance chip above.