TY - THES A1 - Lie, Han Cheng T1 - A Meshless Method for Conformational Analysis Y1 - 2012 ER - TY - JOUR A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Cameron--Martin theorems for sequences of Cauchy-distributed random variables JF - arXiv N2 - Given a sequence of Cauchy-distributed random variables defined by a sequence of location parameters and a sequence of scale parameters, we consider another sequence of random variables that is obtained by perturbing the location or scale parameter sequences. Using a result of Kakutani on equivalence of infinite product measures, we provide sufficient conditions for the equivalence of laws of the two sequences. Y1 - 2016 SP - 1608.03784 ER - TY - GEN A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Cameron--Martin theorems for sequences of Cauchy-distributed random variables N2 - Given a sequence of Cauchy-distributed random variables defined by a sequence of location parameters and a sequence of scale parameters, we consider another sequence of random variables that is obtained by perturbing the location or scale parameter sequences. Using a result of Kakutani on equivalence of infinite product measures, we provide sufficient conditions for the equivalence of laws of the two sequences. T3 - ZIB-Report - 16-40 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60230 SN - 1438-0064 ER - TY - GEN A1 - Quer, Jannes A1 - Lie, Han Cheng T1 - Some connections between importance sampling and enhanced sampling methods in molecular dynamics N2 - Enhanced sampling methods play an important role in molecular dynamics, because they enable the collection of better statistics of rare events that are important in many physical phenomena. We show that many enhanced sampling methods can be viewed as methods for performing importance sampling, by identifying important correspondences between the language of molecular dynamics and the language of probability theory. We illustrate these connections by highlighting the similarities between the rare event simulation method of Hartmann and Schütte (J. Stat. Mech. Theor. Exp., 2012), and the enhanced sampling method of Valsson and Parrinello (Phys. Rev. Lett. 113, 090601). We show that the idea of changing a probability measure is fundamental to both enhanced sampling and importance sampling. T3 - ZIB-Report - 17-30 Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-64289 SN - 1438-0064 ER - TY - THES A1 - Lie, Han Cheng T1 - On a strongly convex approximation of a stochastic optimal control problem for importance sampling of metastable diffusions Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:188-fudissthesis000000101680-4 ER - TY - GEN A1 - Lie, Han Cheng A1 - Sullivan, T. J. A1 - Teckentrup, Aretha T1 - Random forward models and log-likelihoods in Bayesian inverse problems T2 - SIAM/ASA Journal on Uncertainty Quantification N2 - We consider the use of randomised forward models and log-likelihoods within the Bayesian approach to inverse problems. Such random approximations to the exact forward model or log-likelihood arise naturally when a computationally expensive model is approximated using a cheaper stochastic surrogate, as in Gaussian process emulation (kriging), or in the field of probabilistic numerical methods. We show that the Hellinger distance between the exact and approximate Bayesian posteriors is bounded by moments of the difference between the true and approximate log-likelihoods. Example applications of these stability results are given for randomised misfit models in large data applications and the probabilistic solution of ordinary differential equations. T3 - ZIB-Report - 18-03 KW - Bayesian inverse problem KW - random likelihood KW - surrogate model KW - posterior consistency KW - probabilistic numerics KW - uncertainty quantification KW - randomised misfit Y1 - 2018 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-66324 SN - 1438-0064 VL - 6 IS - 4 SP - 1600 EP - 1629 ER - TY - JOUR A1 - Lie, Han Cheng A1 - Sullivan, T. J. A1 - Stuart, Andrew 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 model error, at scales no larger than errors inroduced 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.\ (\textit{Stat.\ Comput.}, 2016), 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. Y1 - 2019 U6 - https://doi.org/10.1007/s11222-019-09898-6 VL - 29 IS - 6 SP - 1265 EP - 1283 ER - TY - JOUR A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Equivalence of weak and strong modes of measures on topological vector spaces JF - Inverse Problems N2 - Modes of a probability measure on an infinite-dimensional Banach space X are often defined by maximising the small-radius limit of the ratio of measures of norm balls. Helin and Burger weakened the definition of such modes by considering only balls with centres in proper subspaces of X, and posed the question of when this restricted notion coincides with the unrestricted one. We generalise these definitions to modes of arbitrary measures on topological vector spaces, defined by arbitrary bounded, convex, neighbourhoods of the origin. We show that a coincident limiting ratios condition is a necessary and sufficient condition for the equivalence of these two types of modes, and show that the coincident limiting ratios condition is satisfied in a wide range of topological vector spaces. Y1 - 2018 U6 - https://doi.org/10.1088/1361-6420/aadef2 VL - 34 IS - 11 SP - 115013 ER - TY - JOUR A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Quasi-invariance of countable products of Cauchy measures under non-unitary dilations JF - Electronic Communications in Probability Y1 - 2018 U6 - https://doi.org/10.1214/18-ECP113 VL - 23 IS - 8 SP - 1 EP - 6 ER - TY - GEN A1 - Bojarovski, Stefan A1 - Hege, Hans-Christian A1 - Lie, Han Cheng A1 - Weber, Marcus T1 - Topological analysis and visualization of scalar functions characterizing conformational transitions of molecules on multiple time-scales T2 - Shape Up 2015 - Exercises in Materials Geometry and Topology, 14-18 Sept. 2015, Berlin, Germany N2 - Molecular processes such as protein folding or ligand-receptor-binding can be understood by analyzing the free energy landscape. Those processes are often metastable, i.e. the molecular systems remain in basins around local minima of the free energy landscape, and in rare cases undergo gauche transitions between metastable states by passing saddle-points of this landscape. By discretizing the configuration space, this can be modeled as a discrete Markov process. One way to compute the transition rates between conformations of a molecular system is by utilizing Transition Path Theory and the concept of committor functions. A fundamental problem from the computational point of view is that many time-scales are involved, ranging from 10^(-14) sec for the fastest motion to 10^(-6) sec or more for conformation changes that cause biological effects. The goal of our work is to provide a better understanding of such transitions in configuration space on various time-scales by analyzing characteristic scalar functions topologically and geometrically. We are developing suitable visualization and interaction techniques to support our analysis. For example, we are analyzing a transition rate indicator function by computing and visualizing its Reeb graph together with the sets of molecular states corresponding to maxima of the transition rate indicator function. A particular challenge is the high dimensionality of the domain which does not allow for a straightforward visualization of the function. The computational topology approach to the analysis of the transition rate indicator functions for a molecular system allows to explore different time scales of the system by utilizing coarser or finer topological partitioning of the function. A specific goal is the development of tools for analyzing the hierarchy of these partitionings. This approach tackles the analysis of a complex and sparse dataset from a different angle than the well-known spectral analysis of Markov State Models. Y1 - 2015 ER - TY - CHAP A1 - Teymur, Onur A1 - Lie, Han Cheng A1 - Sullivan, T. J. A1 - Calderhead, Ben T1 - Implicit probabilistic integrators for ODEs T2 - Advances in Neural Information Processing Systems 31 (NIPS 2018) Y1 - 2018 UR - http://papers.nips.cc/paper/7955-implicit-probabilistic-integrators-for-odes ER - TY - JOUR A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Erratum: Equivalence of weak and strong modes of measures on topological vector spaces (2018 Inverse Problems 34 115013) JF - Inverse Problems Y1 - 2018 U6 - https://doi.org/10.1088/1361-6420/aae55b VL - 34 IS - 12 SP - 129601 ER - TY - JOUR A1 - Lie, Han Cheng A1 - Sullivan, T. J. A1 - Teckentrup, Aretha T1 - Error bounds for some approximate posterior measures in Bayesian inference JF - Numerical Mathematics and Advanced Applications ENUMATH 2019 Y1 - 2020 U6 - https://doi.org/10.1007/978-3-030-55874-1_26 SP - 275 EP - 283 PB - Springer ER - TY - JOUR A1 - Ayanbayev, Birzhan A1 - Klebanov, Ilja A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Γ-convergence of Onsager–Machlup functionals: I. With applications to maximum a posteriori estimation in Bayesian inverse problems JF - Inverse Problems Y1 - 2022 U6 - https://doi.org/10.1088/1361-6420/ac3f81 VL - 38 IS - 2 ER - TY - JOUR A1 - Ayanbayev, Birzhan A1 - Klebanov, Ilja A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Γ-convergence of Onsager–Machlup functionals: II. Infinite product measures on Banach spaces JF - Inverse Problems Y1 - 2022 U6 - https://doi.org/10.1088/1361-6420/ac3f82 VL - 38 IS - 2 ER - TY - JOUR A1 - Koltai, Péter A1 - Lie, Han Cheng A1 - Plonka, Martin T1 - Fréchet 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. Y1 - 2019 U6 - https://doi.org/10.1088/1361-6544/ab1f2a VL - 32 IS - 11 SP - 4232 EP - 4257 ER - TY - JOUR A1 - Lie, Han Cheng A1 - Fackeldey, Konstantin A1 - Weber, Marcus T1 - A Square Root Approximation of Transition Rates for a Markov State Model JF - SIAM. J. Matrix Anal. Appl. Y1 - 2013 U6 - https://doi.org/10.1137/120899959 VL - 34 IS - 2 SP - 738 EP - 756 ER - 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 - GEN A1 - Lie, Han Cheng A1 - Fackeldey, Konstantin A1 - Weber, Marcus T1 - A square root approximation of transition rates for a Markov State Model N2 - Trajectory- or mesh-based methods for analyzing the dynamical behavior of large molecules tend to be impractical due to the curse of dimensionality - their computational cost increases exponentially with the size of the molecule. We propose a method to break the curse by a novel square root approximation of transition rates, Monte Carlo quadrature and a discretization approach based on solving linear programs. With randomly sampled points on the molecular energy landscape and randomly generated discretizations of the molecular configuration space as our initial data, we construct a matrix describing the transition rates between adjacent discretization regions. This transition rate matrix yields a Markov State Model of the molecular dynamics. We use Perron cluster analysis and coarse-graining techniques in order to identify metastable sets in configuration space and approximate the transition rates between the metastable sets. Application of our method to a simple energy landscape on a two-dimensional configuration space provides proof of concept and an example for which we compare the performance of different discretizations. We show that the computational cost of our method grows only polynomially with the size of the molecule. However, finding discretizations of higher-dimensional configuration spaces in which metastable sets can be identified remains a challenge. T3 - ZIB-Report - 13-43 KW - Markov State Models KW - Markov chains KW - meshfree methods KW - metastability KW - Voronoi KW - linear programming Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-42195 SN - 1438-0064 ER - TY - JOUR A1 - Quer, Jannes A1 - Lie, Han Cheng T1 - Some connections between importance sampling and enhanced sampling methods in molecular dynamics JF - Journal of Chemical Physics N2 - Enhanced sampling methods play an important role in molecular dynamics, because they enable the collection of better statistics of rare events that are important in many physical phenomena. We show that many enhanced sampling methods can be viewed as methods for performing importance sampling, by identifying important correspondences between the language of molecular dynamics and the language of probability theory. We illustrate these connections by highlighting the similarities between the rare event simulation method of Hartmann and Schütte (J. Stat. Mech. Theor. Exp., 2012), and the enhanced sampling method of Valsson and Parrinello (Phys. Rev. Lett. 113, 090601). We show that the idea of changing a probability measure is fundamental to both enhanced sampling and importance sampling. Y1 - 2017 ER -