TY - JOUR A1 - Cvetković, Nada A1 - Conrad, Tim A1 - Lie, Han Cheng T1 - A convergent discretization method for transition path theory for diffusion processes JF - Multiscale modeling & simulation : a SIAM interdisciplinary journal N2 - Transition path theory (TPT) for diffusion processes is a framework for analyzing the transitions of multiscale ergodic diffusion processes between disjoint metastable subsets of state space. Most methods for applying TPT involve the construction of a Markov state model on a discretization of state space that approximates the underlying diffusion process. However, the assumption of Markovianity is difficult to verify in practice, and there are to date no known error bounds or convergence results for these methods. We propose a Monte Carlo method for approximating the forward committor, probability current, and streamlines from TPT for diffusion processes. Our method uses only sample trajectory data and partitions of state space based on Voronoi tessellations. It does not require the construction of a Markovian approximating process. We rigorously prove error bounds for the approximate TPT objects and use these bounds to show convergence to their exact counterparts in the limit of arbitrarily fine discretization. We illustrate some features of our method by application to a process that solves the Smoluchowski equation on a triple-well potential. KW - ergodic diffusion processes KW - transition paths KW - rare events KW - Monte Carlo KW - methods Y1 - 2021 U6 - https://doi.org/10.1137/20M1329354 SN - 1540-3459 SN - 1540-3467 VL - 19 IS - 1 SP - 242 EP - 266 PB - Society for Industrial and Applied Mathematics CY - Philadelphia ER - TY - JOUR A1 - Koltai, Peter A1 - Lie, Han Cheng A1 - Plonka, Martin T1 - Frechet differentiable drift dependence of Perron-Frobenius and Koopman operators for non-deterministic dynamics JF - Nonlinearity N2 - We prove the Fréchet differentiability with respect to the drift of Perron–Frobenius and Koopman operators associated to time-inhomogeneous ordinary stochastic differential equations. This result relies on a similar differentiability result for pathwise expectations of path functionals of the solution of the stochastic differential equation, which we establish using Girsanov's formula. We demonstrate the significance of our result in the context of dynamical systems and operator theory, by proving continuously differentiable drift dependence of the simple eigen- and singular values and the corresponding eigen- and singular functions of the stochastic Perron–Frobenius and Koopman operators. KW - stochastic differential equations KW - transfer operator KW - Koopman operator KW - Perron-Frobenius operator KW - smooth drift dependence KW - linear response KW - pathwise expectations Y1 - 2019 U6 - https://doi.org/10.1088/1361-6544/ab1f2a SN - 0951-7715 SN - 1361-6544 VL - 32 IS - 11 SP - 4232 EP - 4257 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Ayanbayev, Birzhan A1 - Klebanov, Ilja A1 - Lie, Han Cheng A1 - Sullivan, Tim J. T1 - Gamma-convergence of Onsager-Machlup functionals BT - II. Infinite product measures on Banach spaces JF - Inverse problems : an international journal of inverse problems, inverse methods and computerised inversion of data N2 - We derive Onsager-Machlup functionals for countable product measures on weighted l(p) subspaces of the sequence space R-N. Each measure in the product is a shifted and scaled copy of a reference probability measure on R that admits a sufficiently regular Lebesgue density. We study the equicoercivity and Gamma-convergence of sequences of Onsager-Machlup functionals associated to convergent sequences of measures within this class. We use these results to establish analogous results for probability measures on separable Banach or Hilbert spaces, including Gaussian, Cauchy, and Besov measures with summability parameter 1 <= p <= 2. Together with part I of this paper, this provides a basis for analysis of the convergence of maximum a posteriori estimators in Bayesian inverse problems and most likely paths in transition path theory. KW - Bayesian inverse problems KW - Gamma-convergence KW - maximum a posteriori KW - estimation KW - Onsager-Machlup functional KW - small ball probabilities KW - transition path theory Y1 - 2021 U6 - https://doi.org/10.1088/1361-6420/ac3f82 SN - 0266-5611 SN - 1361-6420 VL - 38 IS - 2 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Lie, Han Cheng A1 - Stahn, Martin A1 - Sullivan, Tim J. T1 - Randomised one-step time integration methods for deterministic operator differential equations JF - Calcolo N2 - Uncertainty quantification plays an important role in problems that involve inferring a parameter of an initial value problem from observations of the solution. Conrad et al. (Stat Comput 27(4):1065-1082, 2017) proposed randomisation of deterministic time integration methods as a strategy for quantifying uncertainty due to the unknown time discretisation error. We consider this strategy for systems that are described by deterministic, possibly time-dependent operator differential equations defined on a Banach space or a Gelfand triple. Our main results are strong error bounds on the random trajectories measured in Orlicz norms, proven under a weaker assumption on the local truncation error of the underlying deterministic time integration method. Our analysis establishes the theoretical validity of randomised time integration for differential equations in infinite-dimensional settings. KW - Time integration KW - Operator differential equations KW - Randomisation KW - Uncertainty quantification Y1 - 2022 U6 - https://doi.org/10.1007/s10092-022-00457-6 SN - 0008-0624 SN - 1126-5434 VL - 59 IS - 1 PB - Springer CY - Milano ER - TY - JOUR A1 - Lie, Han Cheng A1 - Stuart, A. M. A1 - Sullivan, Tim J. T1 - Strong convergence rates of probabilistic integrators for ordinary differential equations JF - Statistics and Computing N2 - Probabilistic integration of a continuous dynamical system is a way of systematically introducing discretisation error, at scales no larger than errors introduced by standard numerical discretisation, in order to enable thorough exploration of possible responses of the system to inputs. It is thus a potentially useful approach in a number of applications such as forward uncertainty quantification, inverse problems, and data assimilation. We extend the convergence analysis of probabilistic integrators for deterministic ordinary differential equations, as proposed by Conrad et al. (Stat Comput 27(4):1065-1082, 2017. ), to establish mean-square convergence in the uniform norm on discrete- or continuous-time solutions under relaxed regularity assumptions on the driving vector fields and their induced flows. Specifically, we show that randomised high-order integrators for globally Lipschitz flows and randomised Euler integrators for dissipative vector fields with polynomially bounded local Lipschitz constants all have the same mean-square convergence rate as their deterministic counterparts, provided that the variance of the integration noise is not of higher order than the corresponding deterministic integrator. These and similar results are proven for probabilistic integrators where the random perturbations may be state-dependent, non-Gaussian, or non-centred random variables. KW - Probabilistic numerical methods KW - Ordinary differential equations KW - Convergence rates KW - Uncertainty quantification Y1 - 2019 U6 - https://doi.org/10.1007/s11222-019-09898-6 SN - 0960-3174 SN - 1573-1375 VL - 29 IS - 6 SP - 1265 EP - 1283 PB - Springer CY - Dordrecht ER -