{"title":"Concepts in Monte Carlo sampling","authors":"Gabriele Tartero, Werner Krauth","doi":"10.1119/5.0176853","DOIUrl":null,"url":null,"abstract":"We discuss contemporary ideas in Monte Carlo algorithms in the simplified setting of the one-dimensional anharmonic oscillator. After reviewing the connection between molecular dynamics and Monte Carlo, we introduce the Metropolis and the factorized Metropolis algorithms and lifted non-reversible Markov chains. We, furthermore, illustrate the concept of thinning, where moves are accepted by simple bounding potentials rather than the harmonic and quartic contributions to the anharmonic oscillator. We point out the multiple connections of our example algorithms with real-world sampling problems. This paper is self-contained, and Python implementations are provided.","PeriodicalId":7589,"journal":{"name":"American Journal of Physics","volume":"127 1","pages":""},"PeriodicalIF":0.8000,"publicationDate":"2023-09-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"American Journal of Physics","FirstCategoryId":"101","ListUrlMain":"https://doi.org/10.1119/5.0176853","RegionNum":4,"RegionCategory":"教育学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"EDUCATION, SCIENTIFIC DISCIPLINES","Score":null,"Total":0}
引用次数: 0
Abstract
We discuss contemporary ideas in Monte Carlo algorithms in the simplified setting of the one-dimensional anharmonic oscillator. After reviewing the connection between molecular dynamics and Monte Carlo, we introduce the Metropolis and the factorized Metropolis algorithms and lifted non-reversible Markov chains. We, furthermore, illustrate the concept of thinning, where moves are accepted by simple bounding potentials rather than the harmonic and quartic contributions to the anharmonic oscillator. We point out the multiple connections of our example algorithms with real-world sampling problems. This paper is self-contained, and Python implementations are provided.
期刊介绍:
The mission of the American Journal of Physics (AJP) is to publish articles on the educational and cultural aspects of physics that are useful, interesting, and accessible to a diverse audience of physics students, educators, and researchers. Our audience generally reads outside their specialties to broaden their understanding of physics and to expand and enhance their pedagogical toolkits at the undergraduate and graduate levels.