TY - JOUR A1 - Cvetkovic, Nada A1 - Conrad, Tim A1 - Lie, Han Cheng T1 - A Convergent Discretisation Method for Transition Path Theory for Diffusion Processes JF - Multiscale Modeling & Simulation Y1 - 2021 U6 - https://doi.org/10.1137/20M1329354 VL - 19 IS - 1 SP - 242 EP - 266 PB - Society for Industrial and Applied Mathematics ER - TY - JOUR A1 - Grewe, Carl Martin A1 - Liu, Tuo A1 - Kahl, Christoph A1 - Andrea, Hildebrandt A1 - Zachow, Stefan T1 - Statistical Learning of Facial Expressions Improves Realism of Animated Avatar Faces JF - Frontiers in Virtual Reality Y1 - 2021 U6 - https://doi.org/10.3389/frvir.2021.619811 VL - 2 SP - 1 EP - 13 PB - Frontiers ER - TY - JOUR A1 - Nava-Yazdani, Esfandiar A1 - Hege, Hans-Christian A1 - von Tycowicz, Christoph T1 - A Hierarchical Geodesic Model for Longitudinal Analysis on Manifolds JF - Journal of Mathematical Imaging and Vision N2 - In many applications, geodesic hierarchical models are adequate for the study of temporal observations. We employ such a model derived for manifold-valued data to Kendall's shape space. In particular, instead of the Sasaki metric, we adapt a functional-based metric, which increases the computational efficiency and does not require the implementation of the curvature tensor. We propose the corresponding variational time discretization of geodesics and employ the approach for longitudinal analysis of 2D rat skulls shapes as well as 3D shapes derived from an imaging study on osteoarthritis. Particularly, we perform hypothesis test and estimate the mean trends. Y1 - 2022 U6 - https://doi.org/10.1007/s10851-022-01079-x VL - 64 IS - 4 SP - 395 EP - 407 ER - TY - JOUR A1 - Udvary, Daniel A1 - Harth, Philipp A1 - Macke, Jakob H. A1 - Hege, Hans-Christian A1 - de Kock, Christiaan P. J. A1 - Sakmann, Bert A1 - Oberlaender, Marcel T1 - The Impact of Neuron Morphology on Cortical Network Architecture JF - Cell Reports N2 - The neurons in the cerebral cortex are not randomly interconnected. This specificity in wiring can result from synapse formation mechanisms that connect neurons depending on their electrical activity and genetically defined identity. Here, we report that the morphological properties of the neurons provide an additional prominent source by which wiring specificity emerges in cortical networks. This morphologically determined wiring specificity reflects similarities between the neurons’ axo-dendritic projections patterns, the packing density and cellular diversity of the neuropil. The higher these three factors are the more recurrent is the topology of the network. Conversely, the lower these factors are the more feedforward is the network’s topology. These principles predict the empirically observed occurrences of clusters of synapses, cell type-specific connectivity patterns, and nonrandom network motifs. Thus, we demonstrate that wiring specificity emerges in the cerebral cortex at subcellular, cellular and network scales from the specific morphological properties of its neuronal constituents. Y1 - 2022 U6 - https://doi.org/10.1016/j.celrep.2022.110677 VL - 39 IS - 2 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 - 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 - 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 - Zhang, Wei A1 - Hartmann, Carsten A1 - von Kleist, Max T1 - Optimal control of Markov jump processes: Asymptotic analysis, algorithms and applications to the modeling of chemical reaction systems JF - Communications in Mathematical Sciences N2 - Markov jump processes are widely used to model natural and engineered processes. In the context of biological or chemical applications one typically refers to the chemical master equation (CME), which models the evolution of the probability mass of any copy-number combination of the interacting particles. When many interacting particles (“species”) are considered, the complexity of the CME quickly increases, making direct numerical simulations impossible. This is even more problematic when one aims at controlling the Markov jump processes defined by the CME. In this work, we study both open loop and feedback optimal control problems of the Markov jump processes in the case that the controls can only be switched at fixed control stages. Based on Kurtz’s limit theorems, we prove the convergence of the respective control value functions of the underlying Markov decision problem as the copy numbers of the species go to infinity. In the case of the optimal control problem on a finite time-horizon, we propose a hybrid control policy algorithm to overcome the difficulties due to the curse of dimensionality when the copy number of the involved species is large. Two numerical examples demonstrate the suitability of both the analysis and the proposed algorithms. Y1 - 2018 U6 - https://doi.org/10.4310/CMS.2018.v16.n2.a1 SP - 293 EP - 331 ER - TY - JOUR A1 - Sharma, Upanshu A1 - Zhang, Wei T1 - Non-reversible sampling schemes on submanifolds JF - SIAM Journal on Numerical Analysis N2 - Calculating averages with respect to probability measures on submanifolds is often necessary in various application areas such as molecular dynamics, computational statistical mechanics and Bayesian statistics. In recent years, various numerical schemes have been proposed in the literature to study this problem based on appropriate reversible constrained stochastic dynamics. In this paper we present and analyse a non-reversible generalisation of the projection-based scheme developed by one of the authors [ESAIM: M2AN, 54 (2020), pp. 391-430]. This scheme consists of two steps - starting from a state on the submanifold, we first update the state using a non-reversible stochastic differential equation which takes the state away from the submanifold, and in the second step we project the state back onto the manifold using the long-time limit of a ordinary differential equation. We prove the consistency of this numerical scheme and provide quantitative error estimates for estimators based on finite-time running averages. Furthermore, we present theoretical analysis which shows that this scheme outperforms its reversible counterpart in terms of asymptotic variance. We demonstrate our findings on an illustrative test example. Y1 - 2020 U6 - https://doi.org/10.1137/20M1378752 VL - 59 IS - 6 SP - 2989 EP - 3031 ER - TY - JOUR A1 - Zhang, Wei T1 - Ergodic SDEs on submanifolds and related numerical sampling schemes JF - ESAIM: Mathematical Modelling and Numerical Analysis N2 - In many applications, it is often necessary to sample the mean value of certain quantity with respect to a probability measure $\mu$ on the level set of a smooth function ξ:R^d→R^k, 1≤k