TY - GEN A1 - Zhang, Wei A1 - Schütte, Christof T1 - Reliable approximation of long relaxation timescales in molecular dynamics N2 - Many interesting rare events in molecular systems like ligand association, protein folding or con- formational changes happen on timescales that often are not accessible by direct numerical simulation. Therefore rare event approximation approaches like interface sampling, Markov state model building or advanced reaction coordinate based free energy estimation have attracted huge attention recently. In this article we analyze the reliability of such approaches: How precise is an estimate of long relaxation timescales of molecular systems resulting from various forms of rare event approximation methods? Our results give a theoretical answer to this question by relating it with the transfer operator approach to molecular dynamics. By doing so they also allow for understanding deep connections between the different approaches. T3 - ZIB-Report - 17-19 Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-63718 SN - 1438-0064 ER - TY - JOUR A1 - Andrés Arcones, Daniel A1 - Weiser, Martin A1 - Koutsourelakis, Phaedon-Stelios A1 - Unger, Jörg F. T1 - Embedded Model Form Uncertainty Quantification with Measurement Noise for Bayesian Model Calibration JF - Data-Centric Engineering N2 - A key factor in ensuring the accuracy of computer simulations that model physical systems is the proper calibration of their parameters based on real-world observations or experimental data. Inevitably, uncertainties arise, and Bayesian methods provide a robust framework for quantifying and propagating these uncertainties to model predictions. Nevertheless, Bayesian methods paired with inexact models usually produce predictions unable to represent the observed datapoints. Additionally, the quantified uncertainties of these overconfident models cannot be propagated to other Quantities of Interest (QoIs) reliably. A promising solution involves embedding a model inadequacy term in the inference parameters, allowing the quantified model form uncertainty to influence non-observed QoIs. This paper introduces a more interpretable framework for embedding the model inadequacy compared to existing methods. To overcome the limitations of current approaches, we adapt the existing likelihood models to properly account for noise in the measurements and propose two new formulations designed to address their shortcomings. Moreover, we evaluate the performance of this inadequacy-embedding approach in the presence of discrepancies between measurements and model predictions, including noise and outliers. Particular attention is given to how the uncertainty associated with the model inadequacy term propagates to the QoIs, enabling a more comprehensive statistical analysis of prediction’s reliability. Finally, the proposed approach is applied to estimate the uncertainty in the predicted heat flux from a transient thermal simulation using temperature bservations. Y1 - 2026 U6 - https://doi.org/10.1017/dce.2025.10035 VL - 7 ER - TY - JOUR A1 - Lelièvre, Tony A1 - Zhang, Wei T1 - Pathwise estimates for effective dynamics: the case of nonlinear vectorial reaction coordinates JF - Multiscale Modeling and Simulation N2 - Effective dynamics using conditional expectation was proposed in [F. Legoll and T. Lelièvre, Nonlinearity, 2010] to approximate the essential dynamics of high-dimensional diffusion processes along a given reaction coordinate. The approximation error of the effective dynamics when it is used to approximate the behavior of the original dynamics has been considered in recent years. As a continuation of the previous work [F. Legoll, T. Lelièvre, and S. Olla, Stoch. Process. Appl, 2017], in this paper we obtain pathwise estimates for effective dynamics when the reaction coordinate function is either nonlinear or vector-valued. Y1 - 2018 U6 - https://doi.org/10.1137/18M1186034 IS - 17 SP - 1019 EP - 1051 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Schütte, Christof A1 - Zhang, Wei T1 - Jarzynski's equality, fluctuation theorems, and variance reduction: Mathematical analysis and numerical algorithms JF - Journal of Statistical Physics N2 - In this paper, we study Jarzynski's equality and fluctuation theorems for diffusion processes. While some of the results considered in the current work are known in the (mainly physics) literature, we review and generalize these nonequilibrium theorems using mathematical arguments, therefore enabling further investigations in the mathematical community. On the numerical side, variance reduction approaches such as importance sampling method are studied in order to compute free energy differences based on Jarzynski's equality. Y1 - 2018 U6 - https://doi.org/10.1007/s10955-019-02286-4 VL - 175 IS - 6 SP - 1214 EP - 1261 ER - TY - JOUR A1 - Zhang, Wei T1 - Some new results on relative entropy production, time reversal, and optimal control of time-inhomogeneous diffusion processes JF - Journal of Mathematical Physics N2 - This paper studies time-inhomogeneous nonequilibrium diffusion processes, including both Brownian dynamics and Langevin dynamics. We derive upper bounds of the relative entropy production of the time-inhomogeneous process with respect to the transient invariant probability measures. We also study the time reversal of the reverse process in Crooks' fluctuation theorem. We show that the time reversal of the reverse process coincides with the optimally controlled forward process that leads to zero variance importance sampling estimator based on Jarzynski's equality. Y1 - 2021 U6 - https://doi.org/10.1063/5.0038740 VL - 62 IS - 4 ER - TY - JOUR A1 - Lelievre, Tony A1 - Stoltz, Gabriel A1 - Zhang, Wei T1 - Multiple projection MCMC algorithms on submanifolds JF - IMA Journal of Numerical Analysis N2 - We propose new Markov Chain Monte Carlo algorithms to sample probability distributions on submanifolds, which generalize previous methods by allowing the use of set-valued maps in the proposal step of the MCMC algorithms. The motivation for this generalization is that the numerical solvers used to project proposed moves to the submanifold of interest may find several solutions. We show that the new algorithms indeed sample the target probability measure correctly, thanks to some carefully enforced reversibility property. We demonstrate the interest of the new MCMC algorithms on illustrative numerical examples. Y1 - 2022 U6 - https://doi.org/10.1093/imanum/drac006 ER - TY - JOUR A1 - Zhang, Wei A1 - Li, Tiejun A1 - Schütte, Christof T1 - Solving eigenvalue PDEs of metastable diffusion processes using artificial neural networks JF - Journal of Computational Physics N2 - In this paper, we consider the eigenvalue PDE problem of the infinitesimal generators of metastable diffusion processes. We propose a numerical algorithm based on training artificial neural networks for solving the leading eigenvalues and eigenfunctions of such high-dimensional eigenvalue problem. The algorithm is useful in understanding the dynamical behaviors of metastable processes on large timescales. We demonstrate the capability of our algorithm on a high-dimensional model problem, and on the simple molecular system alanine dipeptide. Y1 - 2021 U6 - https://doi.org/10.1016/j.jcp.2022.111377 VL - 465 ER - TY - JOUR A1 - Zhang, Wei A1 - Klus, Stefan A1 - Conrad, Tim A1 - Schütte, Christof T1 - Learning chemical reaction networks from trajectory data JF - SIAM Journal on Applied Dynamical Systems (SIADS) N2 - We develop a data-driven method to learn chemical reaction networks from trajectory data. Modeling the reaction system as a continuous-time Markov chain and assuming the system is fully observed,our method learns the propensity functions of the system with predetermined basis functions by maximizing the likelihood function of the trajectory data under l^1 sparse regularization. We demonstrate our method with numerical examples using synthetic data and carry out an asymptotic analysis of the proposed learning procedure in the infinite-data limit. Y1 - 2019 U6 - https://doi.org/10.1137/19M1265880 VL - 18 IS - 4 SP - 2000 EP - 2046 ER - TY - JOUR A1 - Zhao, Yue A1 - Zhang, Wei A1 - Li, Tiejun T1 - EPR-Net: Constructing non-equilibrium potential landscape via a variational force projection formulation JF - National Science Review N2 - We present EPR-Net, a novel and effective deep learning approach that tackles a crucial challenge in biophysics: constructing potential landscapes for high-dimensional non-equilibrium steady-state (NESS) systems. EPR-Net leverages a nice mathematical fact that the desired negative potential gradient is simply the orthogonal projection of the driving force of the underlying dynamics in a weighted inner-product space. Remarkably, our loss function has an intimate connection with the steady entropy production rate (EPR), enabling simultaneous landscape construction and EPR estimation. We introduce an enhanced learning strategy for systems with small noise, and extend our framework to include dimensionality reduction and state-dependent diffusion coefficient case in a unified fashion. Comparative evaluations on benchmark problems demonstrate the superior accuracy, effectiveness, and robustness of EPR-Net compared to existing methods. We apply our approach to challenging biophysical problems, such as an 8D limit cycle and a 52D multi-stability problem, which provide accurate solutions and interesting insights on constructed landscapes. With its versatility and power, EPR-Net offers a promising solution for diverse landscape construction problems in biophysics. Y1 - 2024 U6 - https://doi.org/10.1093/nsr/nwae052 VL - 11 IS - 7 ER - TY - JOUR A1 - Akhyar, Fatima-Zahrae A1 - Zhang, Wei A1 - Stoltz, Gabriel A1 - Schütte, Christof T1 - Generative modeling of conditional probability distributions on the level-sets of collective variables N2 - Given a probability distribution $\mu$ in $\mathbb{R}^d$ represented by data, we study in this paper the generative modeling of its conditional probability distributions on the level-sets of a collective variable $\xi: \mathbb{R}^d \rightarrow \mathbb{R}^k$, where $1 \le k t0. We obtain promising results assuming Newtonian cooling for boxes filled with air, heaps of clothes, or books in two different experimental scenarios. Two different parameter estimators, (TA0^, t0^) based on four temperature measurements and (TA0*, t0*) for 4N measurements were tested. Our results in a climate chamber were partially appropriate for TTDE input. A decline at time t0 from TA0 = 22.5°C ↓ TA1 = 14°C was reconstructed at t = t0 + 95min with relative deviations ρt0^ = 27% and ρTA0^ = 19% relative to t - t0 and TA0 – TA1 respectively, for N = 1 with span Δt = 50min. For N = 200 in a time interval [t0 + 95min, t0 + 295min] we found ρt0^ = 5% and ρTA0^ = 11% with the same Δt. Further research is necessary to guarantee applicability in routine casework, in particular with respect to more elaborate cooling models, estimation algorithms, and evaluation localization. Y1 - 2026 ER - TY - JOUR A1 - Villani, Paolo A1 - Andrés Arcones, Daniel A1 - Unger, Jörg F. A1 - Weiser, Martin T1 - Gaussian mixture models for model improvement N2 - Modeling complex physical systems such as they arise in civil engineering applications requires finding a trade-off between physical fidelity and practicality. Consequently, deviations of simulation from measurements are ubiquitous even after model calibration due to the model discrepancy, which may result from deliberate modeling decisions, ignorance, or lack of knowledge. If the mismatch between simulation and measurements are deemed unacceptable, the model has to be improved. Targeted model improvement is challenging due to a non-local impact of model discrepancies on measurements and the dependence on sensor configurations. Many approaches to model improvement, such as Bayesian calibration with additive mismatch terms, gray-box models, symbolic regression, or stochastic model updating, often lack interpretability, generalizability, physical consistency, or practical applicability. This paper introduces a non-intrusive approach to model discrepancy analysis using mixture models. Instead of directly modifying the model structure, the method maps sensor readings to clusters of physically meaningful parameters, automatically assigning sensor readings to parameter vector clusters. This mapping can reveal systematic discrepancies and model biases, guiding targeted, physics-based refinements by the modeler. The approach is formulated within a Bayesian framework, enabling the identification of parameter clusters and their assignments via the Expectation-Maximization (EM) algorithm. The methodology is demonstrated through numerical experiments, including an illustrative example and a real-world case study of heat transfer in a concrete bridge. Y1 - 2026 ER -