Fast and accurate rare-event simulation by Markov chain Monte Carlo for general quantum codes
Abstract
The splitting method estimates logical error rates of quantum error-correcting codes deep in the rare-event regime, where direct Monte Carlo sampling is prohibitively expensive. It relies on a Markov chain Monte Carlo (MCMC) sampler of faulty configurations, usually driven by one-bit-flip proposals. However, one-bit flips have been reported to be trapped on some codes, which biases the estimate. To address this issue, we propose a new kernel that combines locally balanced one-bit-flip proposals with a cluster-exchange move. We benchmark both kernels and direct Monte Carlo on twelve codes under code-capacity bit-flip noise in two ways, by comparing the mean fault weight with an MCMC-independent reference and by measuring the computation time to reach a fixed accuracy. In the first, our kernel agrees with the reference on every code, including those where one-bit flips are trapped. In the second, it is faster than one-bit flips wherever they are trapped, by orders of magnitude on some codes. Furthermore, at the lowest error rate, it is faster than direct Monte Carlo on every code, by more than four orders of magnitude on some.
Related
- Quantum Simulation of Hadronization: From String Breaking to Event Generation
- Rare Event Simulation of Quantum Error-Correcting Circuits
- Fast and Accurate Equivariant Foundation Models for Atomistic Simulation
- Neutrino Direct Simulation Monte Carlo: Accurate modeling of out-of-equilibrium decaying cosmological relics
- Reducing simulation-related emissions in rare-event searches through optimised event biasing
BibTeX
@misc{haruna2026fast,
title = {{Fast and accurate rare-event simulation by Markov chain Monte Carlo for general quantum codes}},
author = {Junichi Haruna and Kotaro Murakami},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.35498},
eprint = {2610.35498},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.35498},
}
RIS
TY - UNPB
TI - Fast and accurate rare-event simulation by Markov chain Monte Carlo for general quantum codes
AU - Junichi Haruna
AU - Kotaro Murakami
PY - 2026
DA - 2026/10/26
PB - arXiv
DO - 10.48550/arxiv.2610.35498
UR - https://arxiv.org/abs/2610.35498
AB - The splitting method estimates logical error rates of quantum error-correcting codes deep in the rare-event regime, where direct Monte Carlo sampling is prohibitively expensive. It relies on a Markov chain Monte Carlo (MCMC) sampler of faulty configurations, usually driven by one-bit-flip proposals. However, one-bit flips have been reported to be trapped on some codes, which biases the estimate. To address this issue, we propose a new kernel that combines locally balanced one-bit-flip proposals with a cluster-exchange move. We benchmark both kernels and direct Monte Carlo on twelve codes under code-capacity bit-flip noise in two ways, by comparing the mean fault weight with an MCMC-independent reference and by measuring the computation time to reach a fixed accuracy. In the first, our kernel agrees with the reference on every code, including those where one-bit flips are trapped. In the second, it is faster than one-bit flips wherever they are trapped, by orders of magnitude on some codes. Furthermore, at the lowest error rate, it is faster than direct Monte Carlo on every code, by more than four orders of magnitude on some.
ER -