A Replica Exchange Markov Chain Monte Carlo Method for Disconnected Implicit Manifolds via Tubular Relaxation
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2026
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866918464922320896 |
|---|---|
| author | Wang, Xuyuan Han, Donglin |
| author_facet | Wang, Xuyuan Han, Donglin |
| contents | Markov chain Monte Carlo (MCMC) methods provide powerful framework for sampling unknown probability measures across a wide range of scientific applications. In some settings, the target distribution is supported on a lower-dimensional submanifold of Euclidean space defined by nonlinear constraints, motivating the development of constrained Hamiltonian Monte Carlo (CHMC) methods. Most existing CHMC algorithms rely on the assumption that the implicit manifold is connected, allowing local constrained integrators such as RATTLE to explore the posterior ergodically. In practice, this assumption is occasionally violated due to complex geometric structures induced by nonlinear constraints of a model. We propose a replica exchange MCMC framework that couples a constrained chain evolving on the implicit manifold with a relaxed auxiliary chain defined in a tubular neighborhood of the constraint. The relaxed chain enables transitions between disconnected components. We show that the resulting algorithm enables sampling from a broader class of implicit manifolds, including those with disconnected components. We prove that the proposed sampler satisfies detailed balance, irreducibility, ergodicity, and convergence. We also demonstrate its effectiveness on examples from molecular and biological dynamical systems. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2604_22055 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | A Replica Exchange Markov Chain Monte Carlo Method for Disconnected Implicit Manifolds via Tubular Relaxation Wang, Xuyuan Han, Donglin Numerical Analysis 65C05 Markov chain Monte Carlo (MCMC) methods provide powerful framework for sampling unknown probability measures across a wide range of scientific applications. In some settings, the target distribution is supported on a lower-dimensional submanifold of Euclidean space defined by nonlinear constraints, motivating the development of constrained Hamiltonian Monte Carlo (CHMC) methods. Most existing CHMC algorithms rely on the assumption that the implicit manifold is connected, allowing local constrained integrators such as RATTLE to explore the posterior ergodically. In practice, this assumption is occasionally violated due to complex geometric structures induced by nonlinear constraints of a model. We propose a replica exchange MCMC framework that couples a constrained chain evolving on the implicit manifold with a relaxed auxiliary chain defined in a tubular neighborhood of the constraint. The relaxed chain enables transitions between disconnected components. We show that the resulting algorithm enables sampling from a broader class of implicit manifolds, including those with disconnected components. We prove that the proposed sampler satisfies detailed balance, irreducibility, ergodicity, and convergence. We also demonstrate its effectiveness on examples from molecular and biological dynamical systems. |
| title | A Replica Exchange Markov Chain Monte Carlo Method for Disconnected Implicit Manifolds via Tubular Relaxation |
| topic | Numerical Analysis 65C05 |
| url | https://arxiv.org/abs/2604.22055 |