Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Estimating Committor Functions via Deep Adaptive Sampling on Rare Transition Paths

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read DASTR breaks the transition-data bottleneck in neural committor estimation by sampling from the gradient-weighted Boltzmann density, and its autoencoder variant makes the scheme usable on molecules.

desk verdict Solid adaptive sampling scheme with convincing synthetic gains, but the latent-space variant's stated target distribution is unjustified and needs work before the molecular claims stand. read the letter →

arxiv 2501.15522 v1 pith:XLSJOCBQ submitted 2025-01-26 stat.ML cs.LGq-bio.QM

classification stat.MLcs.LGq-bio.QM
keywords committorfunctiondeepadaptivesamplingrareeventtransitionpathKRnetautoencodercollectivevariablesnormalizingflows
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

The paper addresses a bottleneck in neural committor estimation: transition-state configurations are rarely visited by straightforward molecular dynamics, so training sets are starved of exactly the points that matter. It proposes DASTR (Deep Adaptive Sampling for TRansition paths), an iterative scheme in which the current committor network's gradient magnitude times the Boltzmann weight, $|\nabla q_\theta(x)|^2 e^{-\beta V(x)}$, is treated as an unnormalized density and approximated by a normalizing-flow model whose new samples are added to the training set. In the 10- and 20-dimensional benchmark problems the resulting estimators are reported to be up to an order of magnitude more accurate than sampling from biased dynamics at equal sample counts. For the alanine dipeptide system, running DASTR on autoencoder-learned latent variables yields molecular configurations that pass a 150 kJ/mol energy filter over 97% of the time and cuts molecular sampling time to under 4% of the umbrella-sampling alternative. If these claims hold, DASTR gives a general route to high-quality transition data for committor networks without hand-picked collective variables or repeated SDE simulation.

What carries the argument

The load-bearing object is the target density $p_{V,q}(x)=|\nabla q_\theta(x)|^2 e^{-\beta V(x)}/C$ (or with a bias potential added), which assigns more weight where the committor changes fastest, that is, in the transition-state region. The second object is a normalizing-flow density model (the paper uses KRnet) that approximates this target by minimizing cross-entropy with importance weights and generates new training points through the inverse flow. The third is a trained autoencoder whose latent variables serve as automatically selected collective variables, so that flow samples can be decoded into molecular geometries and then filtered by an energy threshold. Together these pieces turn the variational loss into an adaptive importance-sampling loop: train $q_\theta$, locate where its gradient is large, generate fresh samples there, retrain.

What would settle it

Replace the latent-space identification with the exact change-of-variables formula, multiplying the target density by the encoder Jacobian determinant, and compare the decoded samples and committor error on alanine dipeptide against the reported results; if the corrected sampler does not match or improve the accuracy, the stated mechanism is not what produces the gains. A cheaper check is to compare histograms of $|\nabla q_\theta|^2 e^{-\beta V}$ at decoded versus directly sampled configurations.

Watch

Extended reading notes

Core claim

The committor $q(x)$ is the probability that a trajectory starting at $x$ reaches $B$ before $A$; it solves a high-dimensional linear PDE. DASTR's central move is to define an adaptive target $p_{V,q}(x) \propto |\nabla q_\theta(x)|^2 e^{-\beta V(x)}$, optionally with a bias potential that lowers the barrier, approximate this target by a normalizing-flow model, draw new training points from it, retrain $q_\theta$, and repeat. Because the gradient of the committor concentrates near the transition state, the generated points land preferentially on rare transition paths. For realistic molecular systems the same loop is run on low-dimensional latent collective variables: an autoencoder maps atomic coordinates to latent variables, the flow samples in latent space, and a decoder reconstructs molecular configurations, which are then filtered by potential energy. The paper reports that this latent-space DASTR produces physically valid alanine dipeptide configurations with high probability while reducing molecular sampling time by more than 95% relative to DASTR combined with umbrella sampling.

Load-bearing premise

The load-bearing premise is Eq. (11), which treats the sampling density on the autoencoder's latent variables as equal to the density on atomic coordinates without correcting for the fact that many atomic configurations can map to the same latent point; if that identification fails, the autoencoder version of DASTR samples from a distribution that is not the intended importance target.

Editorial extensions

If this is right

  • At fixed sample size, DASTR reduces the $L^2$ error of the estimated committor by roughly an order of magnitude on the rugged Mueller test, with faster error decay as the sample budget grows.
  • The autoencoder variant removes the need to hand-pick collective variables and avoids SDE-based umbrella sampling for generating transition-state molecules, making the pipeline more general for molecular systems.
  • In the alanine dipeptide test, committor values on the 0.5 level set produced by DASTR cluster around 0.5, whereas a network trained on metadynamics samples alone spreads widely, indicating better resolution of the transition-state ensemble.
  • In the latent-space setting, over 97% of decoded molecular configurations pass the 150 kJ/mol energy filter, while direct flow sampling on heavy-atom coordinates yields almost no physically valid configurations, so the latent space is what makes the construction physical.

Reading between the lines

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

  • The paper's identification $p_{V,q}(s(x))=p_{V,q}(x)$ in Eq. (11) drops the encoder's Jacobian and many-to-one structure; a corrected target that accounts for the change of variables is a natural extension and could be tested on the same benchmarks.
  • Because the sampling target is driven by the current network's gradient, DASTR is a form of self-training; a stopping rule or convergence analysis based on successive changes of the flow target would clarify when the iterative loop has extracted the available transition information.
  • The energy threshold is applied as a post-processing filter rather than being built into the sampling density; including a truncated or reweighted target might raise the fraction of usable samples further.
  • The same 'gradient as density' recipe applies to any variational loss with a nonnegative integrand, so a similar adaptive collocation principle should transfer to other high-dimensional PDE solvers beyond committor estimation.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper proposes DASTR, an iterative deep adaptive sampling method for training neural-network approximations of the committor function. Starting from a standard variational formulation, DASTR defines a sampling density proportional to |∇q_θ|^2 e^{-βV}, approximates it with a normalizing flow (KRnet), and uses newly generated samples to refine the training set. For high-dimensional molecular systems the algorithm is applied in a reduced latent space, either with hand-picked collective variables combined with umbrella sampling or with autoencoder-learned latent variables decoded back to atomic coordinates. Numerical experiments cover a 10D rugged Mueller potential, a 20D Brownian-motion problem, and alanine dipeptide (d=66), reporting improved L2 errors relative to SDE, uniform, and artificial-temperature baselines, as well as a 97% molecular-validity rate and more than 95% sampling-time reduction for the autoencoder variant.

Significance. If the claims hold, DASTR addresses a real bottleneck in neural committor estimation: the scarcity of transition-region data. The paper has clear strengths: the variational loss is standard and correctly derived; the synthetic benchmarks use external reference solutions (FEniCS and an analytic formula); implementation details in Appendix B are unusually complete; and the KRnet sampling procedure is fully specified. The comparison of sampling-time in Table 3 is a useful practical metric. The main technical risk is the latent-space density identification in Section 3.4, which is load-bearing for the alanine-dipeptide claim. The evaluation on the 1/2-isosurface is informative but not sufficient to establish global accuracy. The absence of comparisons with the closest adaptive committor methods [11,14] limits the strength of the claimed improvement. These issues are fixable, so the correct disposition is major revision rather than rejection.

major comments (3)
  1. [Section 3.4, Eq. (11)] The identity p_{V,q}(s(x)) = p_{V,q}(x) is not a valid density transformation. Since the encoder s: R^d -> R^m is many-to-one with m < d, the density induced by p_{V,q} on the latent coordinates is the pushforward ρ(s) = ∫_{x: s(x)=s} p_{V,q}(x) (det(J_s(x) J_s(x)^T))^{-1/2} dσ(x), not p_{V,q}(x) evaluated at a single preimage. As a result, Eq. (14) and Algorithm 2 train the latent KRnet on an inconsistent objective, and the decoded samples are not, in general, samples from a distribution proportional to |∇q|^2 e^{-βV}. The reported alanine-dipeptide results may still hold empirically, but the stated mechanism for the autoencoder variant is not justified by the equations as written. Please derive the correct pushforward and the associated importance weights, or explicitly reformulate the latent-space objective as an approximation and validate it, for example by comparing weighted moments or effective sample sizes against the intended target.
  2. [Section 4.3, Figures 11 and 15] The realistic-system evaluation rests entirely on the 1/2-isosurface of q_θ: test points are selected by the condition |q_θ(x) - 0.5| ≤ 5 × 10^{-5}, and their committor values are then estimated by short SDE simulations. This is a necessary condition but not a sufficient one: a systematically biased q_θ can have a level set on which the empirical PB concentrates near 0.5 while the function is inaccurate elsewhere, and the filtering step uses q_θ itself to select the test points. To support the claim that the latent autoencoder variant provides a good approximation of the committor, the authors should add a more global check, such as evaluating PB on multiple isosurfaces (for example 0.3 and 0.7) with points selected independently of q_θ, or comparing against a reference obtained by an independent method.
  3. [Section 1.1 and Section 4] The closest adaptive committor methods — active importance sampling for variational objectives [11] and deep adaptive committor sampling [14], and arguably the committor-to-committor method [13] — are cited but never used as baselines. The numerical comparisons are limited to SDE sampling, uniform sampling, artificial temperature, and umbrella/metadynamics. Since the paper's central claim concerns the advantage of adaptive sampling over non-adaptive and potential-modification strategies, at least one direct numerical comparison with [14] or [11] is needed to calibrate the significance of the reported gains. Without such a comparison, the improvement over 'baseline methods' is not benchmarked against the state of the art.
minor comments (4)
  1. [Throughout] Typographical errors such as 'commttor', 'standaredd-dimensional', and 'committer' should be corrected.
  2. [Eq. (12)] The proposal density is written as p_IS(s(x_i)) = e^{-βV_modified(x_i)}, which is not a normalized density on the collective variable s. If unnormalized weights are intentional because constant normalization factors cancel in the cross-entropy minimization, this should be stated explicitly.
  3. [Table 2] In the DASTR row with 20 neurons and |S| = 5 × 10^3, the reported mean is 0.0266 with standard deviation 0.1396, so the standard deviation is larger than the mean. This indicates high variance across runs and should be discussed or addressed with more seeds.
  4. [Section 4.3.2] The support or bounding box of the KRnet in the latent space is not specified; since the decoder may map out-of-support latent points to unphysical states, a sentence on how the latent range is determined would improve reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: DASTR's accuracy is checked against external references; adaptive use of q_theta is training feedback, not a relabeled fit.

full rationale

The paper's central derivation is self-contained: DASTR defines a sampling density p_V,q(x) = |nabla q_theta(x)|^2 e^{-beta V(x)} / C that depends on the current neural approximation q_theta, uses a normalizing flow (KRnet) to draw fresh collocation points, and retrains q_theta. This is an adaptive importance-sampling loop rather than a prediction forced by construction, because the reported accuracies are measured against independent references: a finite-element solution for the 10D rugged Mueller potential, the analytical solution for 20D Brownian motion, and SDE-estimated first-hitting probabilities on the 1/2-isosurface for alanine dipeptide. The isosurface validation does select points using q_theta, but the empirical hitting probabilities come from new stochastic simulations, so a wrong q_theta would produce histograms not centered at 0.5. Self-citations to KRnet and to the authors' earlier DAS-PINNs papers describe standard tools and a general adaptive-sampling framework; no uniqueness theorem or unverified prior result is load-bearing, and the text explicitly notes other exact-likelihood generative models could replace KRnet. The main caveat is not circularity: Eq. (11), p_V,q(s(x)) = p_V,q(x), equates a density on atomic coordinates with a density on latent collective variables without a Jacobian or fiber-integration factor, so the latent-space target in Eq. (14) is not the true pushforward measure. This is an omitted mathematical justification for the autoencoder variant, and it should be weighed as a correctness risk, but the empirical claims (97% validity, timing reduction, isosurface histograms) are not equivalent to the inputs of the method. Hence no circular step is exhibited, and the circularity score is 0.

Assumptions & free parameters 4 free parameters · 5 assumptions · 1 invented entities

The core mathematical guarantee is the standard equivalence between the committor PDE and the variational problem; everything else is heuristic. The sampling target p_{V,q} is stipulated, the latent-space density identification is unproven, and several threshold parameters are hand-set. These items are load-bearing for the realistic-system claims but are not backed by independent evidence.

free parameters (4)
  • Energy threshold for generated configuration filtering = 150 kJ/mol
    Section 4.3.2 sets this threshold to discard generated molecular configurations. It is chosen with reference to the 115.5 kJ/mol maximum from umbrella sampling, and it directly determines the reported valid-configuration proportions.
  • Boundary penalty weight λ = 10 or 1000
    In Eq. (4), λ is hand-set per experiment (10 for Mueller and alanine dipeptide, 1000 for Brownian motion) and controls how strongly the boundary conditions are enforced during training.
  • Umbrella force constant k_us = 200 to 400 kJ/mol
    Appendix B.3 increases k_us during training. It determines how tightly generated samples are pulled to the target collective variables in the umbrella sampling refinement.
  • Metadynamics deposition parameters = w = 1.0 kJ/mol, σ = 0.1 rad, period = 1000 steps
    These parameters define the initial training distribution for the alanine dipeptide experiments and the bias potential Vbias used in Eq. (7).
assumptions (5)
  • standard math The committor function satisfies the Fokker-Planck PDE (2), and the variational problem (3) is equivalent to it.
    Appendix A derives the variational principle from the PDE by integration by parts and the Neumann boundary condition, assuming sufficient regularity of q.
  • ad hoc to paper The sampling target p_{V,q} ∝ |∇q|² e^{-βV} is a suitable importance-sampling distribution for the variational loss.
    The paper stipulates this in Eq. (6) without a variance-reduction or convergence guarantee. It is the core algorithmic heuristic and is not derived from first principles.
  • ad hoc to paper In the latent space, p_{V,q}(s(x)) = p_{V,q}(x) defines a valid sampling density on collective variables.
    Eq. (11) identifies a density on atomic coordinates with a density on latent variables without accounting for the encoder Jacobian or the many-to-one nature of the map. This is load-bearing for the realistic-system claims.
  • domain assumption Autoencoder-decoded latent samples produce physically reasonable molecular configurations after an energy filter.
    Section 3.4.2 relies on the autoencoder having learned the manifold of valid heavy-atom coordinates, and the 150 kJ/mol energy filter is an additional ad hoc correction applied to generated samples.
  • domain assumption Metadynamics samples provide a sufficient initial training set covering the transition region.
    Appendix B.3 uses metadynamics with 5000 Gaussian terms to build S0. If this exploration misses relevant transition paths, the adaptive stages start from an incomplete dataset.
invented entities (1)
  • Latent collective variables learned by the autoencoder
    purpose: Low-dimensional sampling space for KRnet so that generated molecular configurations remain physically valid without hand-picked collective variables or SDE simulations.
    The latent variables have no independent physical meaning or external verification outside the paper's own decoder and energy filter. Their validity is demonstrated only through the paper's internal experiments, and the density identification in Eq. (11) is unproven.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Estimating Committor Functions via Deep Adaptive Sampling on Rare Transition Paths." pith.science (2026). https://pith.science/paper/XLSJOCBQ

@misc{pith2026250115522,
  author       = {Pith},
  title        = {Pith review of: Estimating Committor Functions via Deep Adaptive Sampling on Rare Transition Paths},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XLSJOCBQ}},
  note         = {Machine review of arXiv:2501.15522}
}
read the original abstract

The committor functions are central to investigating rare but important events in molecular simulations. It is known that computing the committor function suffers from the curse of dimensionality. Recently, using neural networks to estimate the committor function has gained attention due to its potential for high-dimensional problems. Training neural networks to approximate the committor function needs to sample transition data from straightforward simulations of rare events, which is very inefficient. The scarcity of transition data makes it challenging to approximate the committor function. To address this problem, we propose an efficient framework to generate data points in the transition state region that helps train neural networks to approximate the committor function. We design a Deep Adaptive Sampling method for TRansition paths (DASTR), where deep generative models are employed to generate samples to capture the information of transitions effectively. In particular, we treat a non-negative function in the integrand of the loss functional as an unnormalized probability density function and approximate it with the deep generative model. The new samples from the deep generative model are located in the transition state region and fewer samples are located in the other region. This distribution provides effective samples for approximating the committor function and significantly improves the accuracy. We demonstrate the effectiveness of the proposed method through both simulations and realistic examples.

Figures

Figures reproduced from arXiv: 2501.15522 by the authors.

Figure 1
Figure 1. The schematic of DASTR for computing the committor function. Training a deep neural network qθ to approximate the high-dimensional committor function must use a high-quality dataset (i.e. data points from the transition area). Typically, the data points from Langevin dynamics are not in the transition state region since the transition between two metastable states is rare and difficult to sample. The proposed DASTR … view at source ↗
Figure 2
Figure 2. Molecular configurations generated by two different settings in DASTR: (a) the inputs of KRnet are [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. The schematic of adaptive sampling in the latent space. We first train an autoencoder to obtain the latent variables as the collective variables (CVs), and then use KRnet to approximate the distribution of the CVs. After training KRnet, we use a random sample z0 from the standard Gaussian distribution to generate a new sample of latent CVs. We can feed this new sample of latent CVs into the decoder to obtain a new s… view at source ↗
Figures from the paper (12 more)
Figure 4
Figure 4. Figure 4: DASTR, samples for the 10-dimensional rugged Mueller potential problem. The red line denotes the [PITH_FULL_IMAGE:figures/full_fig_p013_4.png]
Figure 5
Figure 5. Figure 5: Solutions, 10-dimensional rugged Mueller potential test problem. Figure 5a shows the error behavior of different methods. In Figure 5b-5d, we compare the reference solution qref obtained by the finite element method, the DASTR solution given by 4 × 105 samples and the …
Figure 6
Figure 6. Figure 6: The error w.r.t. sample size |S|. 4.2. Standard Brownian Motion In this test problem, we consider the committor function under the standard Brownian motion [46, 47]. For a stochastic process (Xt)t≥0 ∈ R d , which is a standard Brownian motion starting at x ∈ R d , that…
Figure 7
Figure 7. Figure 7: Solutions evaluated along a curve and the behavior of relative errors, [PITH_FULL_IMAGE:figures/full_fig_p015_7.png]
Figure 8
Figure 8. Figure 8: Histogram of the norm of samples, 20-dimensional test problem [PITH_FULL_IMAGE:figures/full_fig_p016_8.png]
Figure 9
Figure 9. Figure 9: The two metastable states and two transition states of the alanine dipeptide. [PITH_FULL_IMAGE:figures/full_fig_p017_9.png]
Figure 10
Figure 10. Figure 10: Samples during training for the alanine dipeptide test problem. We use DASTR to generate target [PITH_FULL_IMAGE:figures/full_fig_p018_10.png]
Figure 11
Figure 11. Figure 11: The alanine dipeptide test problem: the histograms of the committor function values on the [PITH_FULL_IMAGE:figures/full_fig_p019_11.png]
Figure 12
Figure 12. Figure 12: Visualization of the latent collective variables, the two metastable states [PITH_FULL_IMAGE:figures/full_fig_p020_12.png]
Figure 13
Figure 13. Figure 13: The proportions of valid molecular configurations for two different settings in DASTR are shown. [PITH_FULL_IMAGE:figures/full_fig_p021_13.png]
Figure 14
Figure 14. Figure 14: Samples during training for different latent dimensions, the alanine dipeptide test problem. The [PITH_FULL_IMAGE:figures/full_fig_p022_14.png]
Figure 15
Figure 15. Figure 15: Conducting DASTR in the latent space for the alanine dipeptide test problem: the histograms of [PITH_FULL_IMAGE:figures/full_fig_p022_15.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Commitment Before Realization: When Classifier-Free Guidance Becomes Unnecessary in Masked Diffusion Language Models

    cs.CL 2026-08 conditional novelty 6.0 of 10

    A prompt-specific commitment horizon, identified by comparing guided versus base-only continuations, marks an early point where classifier-free guidance can be removed with little loss in constraint success.

Reference graph

Works this paper leans on

55 extracted references · 45 canonical work pages · cited by 1 Pith paper

  1. [11]

    G. M. Rotskoff, A. R. Mitchell, E. Vanden-Eijnden, Active importance sampling for vari- ational objectives dominated by rare events: Consequences for optimization and general- ization, in: Mathematical and Scientific Machine Learning, PMLR, 2022, pp. 757–780. 23

  2. [14]

    B. Lin, W. Ren, Deep learning method for computing committor functions with adaptive sampling, arXiv preprint arXiv:2404.06206 (2024)

  3. [13]

    P. Kang, E. Trizio, M. Parrinello, Computing the committor with the committor to study the transition state ensemble, Nature Computational Science (2024) 1–10

  4. [1]

    Okuyama-Yoshida, M

    N. Okuyama-Yoshida, M. Nagaoka, T. Yamabe, Transition-state optimization on free en- ergy surface: Toward solution chemical reaction ergodography, International Journal of Quantum Chemistry 70 (1) (1998) 95–103

  5. [2]

    W. E, E. Vanden-Eijnden, Towards a theory of transition paths, Journal of Statistical Physics 123 (3) (2006) 503–523

  6. [3]

    Berteotti, A

    A. Berteotti, A. Cavalli, D. Branduardi, F. L. Gervasio, M. Recanatini, M. Parrinello, Pro- tein conformational transitions: the closure mechanism of a kinase explored by atomistic simulations, Journal of the American Chemical Society 131 (1) (2009) 244–250

  7. [4]

    W. E, E. Vanden-Eijnden, Transition-path theory and path-finding algorithms for the study of rare events., Annual Review of Physical Chemistry 61 (2010) 391–420

  8. [5]

    R. Lai, J. Lu, Point cloud discretization of Fokker–Planck operators for committor func- tions, Multiscale Modeling & Simulation 16 (2) (2018) 710–726

Show all 55 references
  1. [6]

    Q. Li, B. Lin, W. Ren, Computing committor functions for the study of rare events using deep learning, The Journal of Chemical Physics 151 (5) (2019) 054112

  2. [7]

    Y. Chen, J. Hoskins, Y. Khoo, M. Lindsey, Committor functions via tensor networks, Journal of Computational Physics 472 (2023) 111646

  3. [8]

    Y. Khoo, J. Lu, L. Ying, Solving for high-dimensional committor functions using artificial neural networks, Research in the Mathematical Sciences 6 (1) (2019) 1–13

  4. [9]

    H. Li, Y. Khoo, Y. Ren, L. Ying, A semigroup method for high dimensional committor functions based on neural network, in: Mathematical and Scientific Machine Learning, PMLR, 2022, pp. 598–618

  5. [10]

    H. Li, Y. Khoo, Y. Ren, L. Ying, Solving for high dimensional committor functions using neural network with online approximation to derivatives, arXiv preprint arXiv:2012.06727 (2020)

  6. [12]

    M. R. Hasyim, C. H. Batton, K. K. Mandadapu, Supervised learning and the finite- temperature string method for computing committor functions and reaction rates, The Journal of Chemical Physics 157 (18) (2022)

  7. [15]

    W. Gao, C. Wang, Active learning based sampling for high-dimensional nonlinear partial differential equations, Journal of Computational Physics 475 (2023) 111848

  8. [16]

    K. Tang, X. Wan, C. Yang, DAS-PINNs: A deep adaptive sampling method for solv- ing high-dimensional partial differential equations, Journal of Computational Physics 476 (2023) 111868

  9. [17]

    X. Wang, K. Tang, J. Zhai, X. Wan, C. Yang, Deep Adaptive Sampling for Surrogate Modeling Without Labeled Data, Journal of Scientific Computing 101 (3) (2024) 77.doi: 10.1007/s10915-024-02711-1

  10. [18]

    K. Tang, J. Zhai, X. Wan, C. Yang, Adversarial adaptive sampling: Unify PINN and op- timal transport for the approximation of PDEs, in: The Twelfth International Conference on Learning Representations, 2024

  11. [19]

    Z. Gao, L. Yan, T. Zhou, Failure-informed adaptive sampling for pinns, SIAM Journal on Scientific Computing 45 (4) (2023) A1971–A1994

  12. [20]

    Y. Jiao, D. Li, X. Lu, J. Z. Yang, C. Yuan, A Gaussian mixture distribution-based adap- tive sampling method for physics-informed neural networks, Engineering Applications of Artificial Intelligence 135 (2024) 108770

  13. [21]

    Czibula, A.-I

    G. Czibula, A.-I. Albu, M. I. Bocicor, C. Chira, Autoppi: An ensemble of deep autoen- coders for protein–protein interaction prediction, Entropy 23 (6) (2021) 643

  14. [22]

    F. F. Alam, T. Rahman, A. Shehu, Learning reduced latent representations of protein structure data, in: Proceedings of the 10th ACM International Conference on Bioinfor- matics, Computational Biology and Health Informatics, 2019, pp. 592–597

  15. [23]

    Hawkins-Hooker, F

    A. Hawkins-Hooker, F. Depardieu, S. Baur, G. Couairon, A. Chen, D. Bikard, Generating functional protein variants with variational autoencoders, PLoS Computational Biology 17 (2) (2021) e1008736

  16. [24]

    Sirignano, K

    J. Sirignano, K. Spiliopoulos, DGM: A deep learning algorithm for solving partial differ- ential equations, Journal of Computational Physics 375 (2018) 1339–1364

  17. [25]

    Raissi, P

    M. Raissi, P. Perdikaris, G. E. Karniadakis, Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations, Journal of Computational Physics 378 (2019) 686–707. 24

  18. [26]

    G. E. Karniadakis, I. G. Kevrekidis, L. Lu, P. Perdikaris, S. Wang, L. Yang, Physics- informed machine learning, Nature Reviews Physics 3 (6) (2021) 422–440

  19. [27]

    W. E, B. Yu, The deep Ritz method: A deep learning-based numerical algorithm for solving variational problems, Communications in Mathematics and Statistics 6 (1) (2018) 1–12

  20. [28]

    Y. Liao, P. Ming, Deep Nitsche Method: Deep Ritz Method with Essential Boundary Conditions, Communications in Computational Physics 29 (5) (2021) 1365–1384

  21. [29]

    Y. Lu, J. Lu, M. Wang, A priori generalization analysis of the deep Ritz method for solving high dimensional elliptic partial differential equations, in: Conference on Learning Theory, PMLR, 2021, pp. 3196–3241

  22. [30]

    J.Kästner, Umbrellasampling, WileyInterdisciplinaryReviews: ComputationalMolecular Science 1 (6) (2011) 932–942

  23. [31]

    Bussi, A

    G. Bussi, A. Laio, Using metadynamics to explore complex free-energy landscapes, Nature Reviews Physics 2 (4) (2020) 200–212

  24. [32]

    A.Barducci, G.Bussi, M.Parrinello, Well-temperedmetadynamics: asmoothlyconverging and tunable free-energy method, Physical Review Letters 100 (2) (2008) 020603

  25. [33]

    L. Dinh, J. Sohl-Dickstein, S. Bengio, Density estimation using real NVP, arXiv preprint arXiv:1605.08803 (2016)

  26. [34]

    D. P. Kingma, P. Dhariwal, Glow: Generative flow with invertible 1x1 convolutions, in: Advances in Neural Information Processing Systems, 2018, pp. 10215–10224

  27. [35]

    T. Q. Chen, Y. Rubanova, J. Bettencourt, D. K. Duvenaud, Neural ordinary differential equations, in: Advances in Neural Information Processing Systems, 2018, pp. 6571–6583

  28. [36]

    Y. Song, J. Sohl-Dickstein, D. P. Kingma, A. Kumar, S. Ermon, B. Poole, Score-based generative modeling through stochastic differential equations, in: International Conference on Learning Representations, 2021

  29. [37]

    K.Tang, X.Wan, Q.Liao, Deepdensityestimationviainvertibleblock-triangularmapping, Theoretical & Applied Mechanics Letters 10 (2020) 143–148

  30. [38]

    X. Wan, S. Wei, VAE-KRnet and its applications to variational Bayes, Communications in Computational Physics 31 (4) (2022) 1049–1082

  31. [39]

    K. Tang, X. Wan, Q. Liao, Adaptive deep density approximation for Fokker-Planck equa- tions, Journal of Computational Physics 457 (2022) 111080

  32. [40]

    De Boer, D

    P.-T. De Boer, D. P. Kroese, S. Mannor, R. Y. Rubinstein, A tutorial on the cross-entropy method, Annals of Operations Research 134 (1) (2005) 19–67

  33. [41]

    R. Y. Rubinstein, D. P. Kroese, The cross-entropy method: a unified approach to combi- natorial optimization, Monte-Carlo simulation and machine learning, Springer Science & Business Media, 2013. 25

  34. [42]

    Fiorin, M

    G. Fiorin, M. L. Klein, J. Hénin, Using collective variables to drive molecular dynamics simulations, Molecular Physics 111 (22-23) (2013) 3345–3362

  35. [43]

    Schrödinger, The PyMOL Molecular Graphics System, Version 1.8., (No Title) (2015)

    L. Schrödinger, The PyMOL Molecular Graphics System, Version 1.8., (No Title) (2015)

  36. [44]

    Alnæs, J

    M. Alnæs, J. Blechta, J. Hake, A. Johansson, B. Kehlet, A. Logg, C. Richardson, J. Ring, M.E.Rognes, G.N.Wells, TheFEniCSProjectVersion1.5, ArchiveofNumericalSoftware 3 (100) (2015)

  37. [45]

    Logg, K.-A

    A. Logg, K.-A. Mardal, G. Wells, Automated solution of differential equations by the finite element method: The FEniCS book, Vol. 84, Springer Science & Business Media, 2012

  38. [46]

    Hartmann, O

    C. Hartmann, O. Kebiri, L. Neureither, L. Richter, Variational approach to rare event simulation using least-squares regression, Chaos 29 (6) (2019)

  39. [47]

    Nüsken, L

    N. Nüsken, L. Richter, Interpolating between BSDEs and PINNs: Deep learning for elliptic and parabolic boundary value problems, Journal of Machine Learning 2 (1) (2023) 31–64

  40. [48]

    Vershynin, High-dimensional probability: An introduction with applications in data science, Vol

    R. Vershynin, High-dimensional probability: An introduction with applications in data science, Vol. 47, Cambridge University Press, 2018

  41. [49]

    Ma, High-dimensional data analysis with low-dimensional models: Principles, Computation, and Applications, Cambridge University Press, 2022

    J.Wright, Y. Ma, High-dimensional data analysis with low-dimensional models: Principles, Computation, and Applications, Cambridge University Press, 2022

  42. [50]

    S. Jo, T. Kim, V. G. Iyer, W. Im, Charmm-gui: a web-based graphical user interface for charmm, Journal of Computational Chemistry 29 (11) (2008) 1859–1865

  43. [51]

    B. R. Brooks, C. L. Brooks III, A. D. Mackerell Jr, L. Nilsson, R. J. Petrella, B. Roux, Y. Won, G. Archontis, C. Bartels, S. Boresch, et al., Charmm: the biomolecular simulation program, Journal of Computational Chemistry 30 (10) (2009) 1545–1614

  44. [52]

    J. Lee, X. Cheng, S. Jo, A. D. MacKerell, J. B. Klauda, W. Im, Charmm-gui input generator for namd, gromacs, amber, openmm, and charmm/openmm simulations using the charmm36 additive force field, Biophysical Journal 110 (3) (2016) 641a

  45. [53]

    McInnes, J

    L. McInnes, J. Healy, J. Melville, Umap: Uniform manifold approximation and projection for dimension reduction, arXiv preprint arXiv:1802.03426 (2018)

  46. [54]

    L. Zeng, X. Wan, T. Zhou, Bounded KRnet and its applications to density estimation and approximation, arXiv:2305.09063 (2023)

  47. [55]

    Eastman, J

    P. Eastman, J. Swails, J. D. Chodera, R. T. McGibbon, Y. Zhao, K. A. Beauchamp, L.-P. Wang, A. C. Simmonett, M. P. Harrigan, C. D. Stern, et al., Openmm 7: Rapid development of high performance algorithms for molecular dynamics, PLoS Computational Biology 13 (7) (2017) e100565...

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.