Chin-Yi Lin, Kuo-Chin Chen, Philippe Lemey, Marc A Suchard, Andrew J Holbrook, Min-Hsiu Hsieh
{"title":"多提案MCMC的量子加速。","authors":"Chin-Yi Lin, Kuo-Chin Chen, Philippe Lemey, Marc A Suchard, Andrew J Holbrook, Min-Hsiu Hsieh","doi":"10.1214/25-ba1546","DOIUrl":null,"url":null,"abstract":"<p><p>Multiproposal Markov chain Monte Carlo (MCMC) algorithms choose from multiple proposals to generate their next chain step in order to sample from challenging target distributions more efficiently. However, on classical machines, these algorithms require <math><mi>𝒪</mi> <mo>(</mo> <mi>P</mi> <mo>)</mo></math> target evaluations for each Markov chain step when choosing from <math><mi>P</mi></math> proposals. Recent work demonstrates the possibility of quadratic quantum speedups for one such multiproposal MCMC algorithm. After generating <math><mi>P</mi></math> proposals, this quantum parallel MCMC (QPMCMC) algorithm requires only <math><mi>𝒪</mi> <mo>(</mo> <msqrt><mi>P</mi></msqrt> <mo>)</mo></math> target evaluations at each step, outperforming its classical counterpart. However, generating <math><mi>P</mi></math> proposals using classical computers still requires <math><mi>𝒪</mi> <mo>(</mo> <mi>P</mi> <mo>)</mo></math> time complexity, resulting in the overall complexity of QPMCMC remaining <math><mi>𝒪</mi> <mo>(</mo> <mi>P</mi> <mo>)</mo></math> . Here, we present a new, faster quantum multiproposal MCMC strategy, QPMCMC2. With a specially designed Tjelmeland distribution that generates proposals close to the input state, QPMCMC2 requires only <math><mi>𝒪</mi> <mo>(</mo> <mn>1</mn> <mo>)</mo></math> target evaluations and <math><mi>𝒪</mi> <mo>(</mo> <mtext>log</mtext> <mspace></mspace> <mi>P</mi> <mo>)</mo></math> qubits when computing over a large number of proposals <math><mi>P</mi></math> . Unlike its slower predecessor, the QPMCMC2 Markov kernel (1) maintains detailed balance exactly and (2) is fully explicit for a large class of graphical models. We demonstrate this flexibility by applying QPMCMC2 to novel Ising-type models built on bacterial evolutionary networks and obtain significant speedups for Bayesian ancestral trait reconstruction for 248 observed salmonella bacteria.</p>","PeriodicalId":55398,"journal":{"name":"Bayesian Analysis","volume":" ","pages":""},"PeriodicalIF":2.5000,"publicationDate":"2025-08-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12456418/pdf/","citationCount":"0","resultStr":"{\"title\":\"Quantum Speedups for Multiproposal MCMC.\",\"authors\":\"Chin-Yi Lin, Kuo-Chin Chen, Philippe Lemey, Marc A Suchard, Andrew J Holbrook, Min-Hsiu Hsieh\",\"doi\":\"10.1214/25-ba1546\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<p><p>Multiproposal Markov chain Monte Carlo (MCMC) algorithms choose from multiple proposals to generate their next chain step in order to sample from challenging target distributions more efficiently. However, on classical machines, these algorithms require <math><mi>𝒪</mi> <mo>(</mo> <mi>P</mi> <mo>)</mo></math> target evaluations for each Markov chain step when choosing from <math><mi>P</mi></math> proposals. Recent work demonstrates the possibility of quadratic quantum speedups for one such multiproposal MCMC algorithm. After generating <math><mi>P</mi></math> proposals, this quantum parallel MCMC (QPMCMC) algorithm requires only <math><mi>𝒪</mi> <mo>(</mo> <msqrt><mi>P</mi></msqrt> <mo>)</mo></math> target evaluations at each step, outperforming its classical counterpart. However, generating <math><mi>P</mi></math> proposals using classical computers still requires <math><mi>𝒪</mi> <mo>(</mo> <mi>P</mi> <mo>)</mo></math> time complexity, resulting in the overall complexity of QPMCMC remaining <math><mi>𝒪</mi> <mo>(</mo> <mi>P</mi> <mo>)</mo></math> . Here, we present a new, faster quantum multiproposal MCMC strategy, QPMCMC2. With a specially designed Tjelmeland distribution that generates proposals close to the input state, QPMCMC2 requires only <math><mi>𝒪</mi> <mo>(</mo> <mn>1</mn> <mo>)</mo></math> target evaluations and <math><mi>𝒪</mi> <mo>(</mo> <mtext>log</mtext> <mspace></mspace> <mi>P</mi> <mo>)</mo></math> qubits when computing over a large number of proposals <math><mi>P</mi></math> . Unlike its slower predecessor, the QPMCMC2 Markov kernel (1) maintains detailed balance exactly and (2) is fully explicit for a large class of graphical models. We demonstrate this flexibility by applying QPMCMC2 to novel Ising-type models built on bacterial evolutionary networks and obtain significant speedups for Bayesian ancestral trait reconstruction for 248 observed salmonella bacteria.</p>\",\"PeriodicalId\":55398,\"journal\":{\"name\":\"Bayesian Analysis\",\"volume\":\" \",\"pages\":\"\"},\"PeriodicalIF\":2.5000,\"publicationDate\":\"2025-08-18\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12456418/pdf/\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Bayesian Analysis\",\"FirstCategoryId\":\"100\",\"ListUrlMain\":\"https://doi.org/10.1214/25-ba1546\",\"RegionNum\":2,\"RegionCategory\":\"数学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"MATHEMATICS, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Bayesian Analysis","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1214/25-ba1546","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
Multiproposal Markov chain Monte Carlo (MCMC) algorithms choose from multiple proposals to generate their next chain step in order to sample from challenging target distributions more efficiently. However, on classical machines, these algorithms require target evaluations for each Markov chain step when choosing from proposals. Recent work demonstrates the possibility of quadratic quantum speedups for one such multiproposal MCMC algorithm. After generating proposals, this quantum parallel MCMC (QPMCMC) algorithm requires only target evaluations at each step, outperforming its classical counterpart. However, generating proposals using classical computers still requires time complexity, resulting in the overall complexity of QPMCMC remaining . Here, we present a new, faster quantum multiproposal MCMC strategy, QPMCMC2. With a specially designed Tjelmeland distribution that generates proposals close to the input state, QPMCMC2 requires only target evaluations and qubits when computing over a large number of proposals . Unlike its slower predecessor, the QPMCMC2 Markov kernel (1) maintains detailed balance exactly and (2) is fully explicit for a large class of graphical models. We demonstrate this flexibility by applying QPMCMC2 to novel Ising-type models built on bacterial evolutionary networks and obtain significant speedups for Bayesian ancestral trait reconstruction for 248 observed salmonella bacteria.
期刊介绍:
Bayesian Analysis is an electronic journal of the International Society for Bayesian Analysis. It seeks to publish a wide range of articles that demonstrate or discuss Bayesian methods in some theoretical or applied context. The journal welcomes submissions involving presentation of new computational and statistical methods; critical reviews and discussions of existing approaches; historical perspectives; description of important scientific or policy application areas; case studies; and methods for experimental design, data collection, data sharing, or data mining.
Evaluation of submissions is based on importance of content and effectiveness of communication. Discussion papers are typically chosen by the Editor in Chief, or suggested by an Editor, among the regular submissions. In addition, the Journal encourages individual authors to submit manuscripts for consideration as discussion papers.