{"title":"利用集合变量和Jarzynski-Crooks路径对亚稳态系统进行采样","authors":"C. Schönle , M. Gabrié , T. Lelièvre , G. Stoltz","doi":"10.1016/j.jcp.2025.113806","DOIUrl":null,"url":null,"abstract":"<div><div>We consider the problem of sampling a high dimensional multimodal target probability measure. We assume that a good proposal kernel to move only a subset of the degrees of freedoms (also known as collective variables) is known a priori. This proposal kernel can for example be built using normalizing flows <span><span>[32]</span></span>, <span><span>[27]</span></span>, <span><span>[16]</span></span>. We show how to extend the move from the collective variable space to the full space and how to implement an accept-reject step in order to get a reversible chain with respect to a target probability measure. The accept-reject step does not require to know the marginal of the original measure in the collective variable (namely to know the free energy). The obtained algorithm admits several variants, some of them being very close to methods which have been proposed previously in the literature, in particular in <span><span>[3]</span></span>, <span><span>[29]</span></span>, <span><span>[6]</span></span>, <span><span>[7]</span></span>, <span><span>[31]</span></span>. We show how the obtained acceptance ratio can be expressed in terms of the work which appears in the Jarzynski–Crooks equality, at least for some variants. Numerical illustrations demonstrate the efficiency of the approach on various simple test cases and allow us to compare the variants of the algorithm.</div></div>","PeriodicalId":352,"journal":{"name":"Journal of Computational Physics","volume":"527 ","pages":"Article 113806"},"PeriodicalIF":3.8000,"publicationDate":"2025-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Sampling metastable systems using collective variables and Jarzynski–Crooks paths\",\"authors\":\"C. Schönle , M. Gabrié , T. Lelièvre , G. Stoltz\",\"doi\":\"10.1016/j.jcp.2025.113806\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>We consider the problem of sampling a high dimensional multimodal target probability measure. We assume that a good proposal kernel to move only a subset of the degrees of freedoms (also known as collective variables) is known a priori. This proposal kernel can for example be built using normalizing flows <span><span>[32]</span></span>, <span><span>[27]</span></span>, <span><span>[16]</span></span>. We show how to extend the move from the collective variable space to the full space and how to implement an accept-reject step in order to get a reversible chain with respect to a target probability measure. The accept-reject step does not require to know the marginal of the original measure in the collective variable (namely to know the free energy). The obtained algorithm admits several variants, some of them being very close to methods which have been proposed previously in the literature, in particular in <span><span>[3]</span></span>, <span><span>[29]</span></span>, <span><span>[6]</span></span>, <span><span>[7]</span></span>, <span><span>[31]</span></span>. We show how the obtained acceptance ratio can be expressed in terms of the work which appears in the Jarzynski–Crooks equality, at least for some variants. Numerical illustrations demonstrate the efficiency of the approach on various simple test cases and allow us to compare the variants of the algorithm.</div></div>\",\"PeriodicalId\":352,\"journal\":{\"name\":\"Journal of Computational Physics\",\"volume\":\"527 \",\"pages\":\"Article 113806\"},\"PeriodicalIF\":3.8000,\"publicationDate\":\"2025-02-03\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Journal of Computational Physics\",\"FirstCategoryId\":\"101\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0021999125000890\",\"RegionNum\":2,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational Physics","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0021999125000890","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
Sampling metastable systems using collective variables and Jarzynski–Crooks paths
We consider the problem of sampling a high dimensional multimodal target probability measure. We assume that a good proposal kernel to move only a subset of the degrees of freedoms (also known as collective variables) is known a priori. This proposal kernel can for example be built using normalizing flows [32], [27], [16]. We show how to extend the move from the collective variable space to the full space and how to implement an accept-reject step in order to get a reversible chain with respect to a target probability measure. The accept-reject step does not require to know the marginal of the original measure in the collective variable (namely to know the free energy). The obtained algorithm admits several variants, some of them being very close to methods which have been proposed previously in the literature, in particular in [3], [29], [6], [7], [31]. We show how the obtained acceptance ratio can be expressed in terms of the work which appears in the Jarzynski–Crooks equality, at least for some variants. Numerical illustrations demonstrate the efficiency of the approach on various simple test cases and allow us to compare the variants of the algorithm.
期刊介绍:
Journal of Computational Physics thoroughly treats the computational aspects of physical problems, presenting techniques for the numerical solution of mathematical equations arising in all areas of physics. The journal seeks to emphasize methods that cross disciplinary boundaries.
The Journal of Computational Physics also publishes short notes of 4 pages or less (including figures, tables, and references but excluding title pages). Letters to the Editor commenting on articles already published in this Journal will also be considered. Neither notes nor letters should have an abstract.