Aristotelis P. Sgouros*, and , Doros N. Theodorou,
{"title":"为粒子场布朗动力学模拟开发基于无网格核的方案","authors":"Aristotelis P. Sgouros*, and , Doros N. Theodorou, ","doi":"10.1021/acs.jpcb.4c01441","DOIUrl":null,"url":null,"abstract":"<p >We develop a meshless discretization scheme for particle-field Brownian dynamics simulations. The density is assigned on the particle level using a weighting kernel with finite support. The system’s free energy density is derived from an equation of state (EoS) and includes a square gradient term. The numerical stability of the scheme is evaluated in terms of reproducing the thermodynamics (equilibrium density and compressibility) and dynamics (diffusion coefficient) of homogeneous samples. Using a reduced description to simplify our analysis, we find that numerical stability depends strictly on reduced reference compressibility, kernel range, time step in relation to the friction factor, and reduced external pressure, the latter being relevant under isobaric conditions. Appropriate parametrization yields precise thermodynamics, further improved through a simple renormalization protocol. The dynamics can be restored exactly through a trivial manipulation of the time step and friction coefficient. A semiempirical formula for the upper bound on the time step is derived, which takes into account variations in compressibility, friction factor, and kernel range. We test the scheme on realistic mesoscopic models of fluids, involving both simple (Helfand) and more sophisticated (Sanchez–Lacombe) equations of state.</p>","PeriodicalId":60,"journal":{"name":"The Journal of Physical Chemistry B","volume":null,"pages":null},"PeriodicalIF":2.8000,"publicationDate":"2024-07-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/epdf/10.1021/acs.jpcb.4c01441","citationCount":"0","resultStr":"{\"title\":\"Development of a Meshless Kernel-Based Scheme for Particle-Field Brownian Dynamics Simulations\",\"authors\":\"Aristotelis P. Sgouros*, and , Doros N. Theodorou, \",\"doi\":\"10.1021/acs.jpcb.4c01441\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<p >We develop a meshless discretization scheme for particle-field Brownian dynamics simulations. The density is assigned on the particle level using a weighting kernel with finite support. The system’s free energy density is derived from an equation of state (EoS) and includes a square gradient term. The numerical stability of the scheme is evaluated in terms of reproducing the thermodynamics (equilibrium density and compressibility) and dynamics (diffusion coefficient) of homogeneous samples. Using a reduced description to simplify our analysis, we find that numerical stability depends strictly on reduced reference compressibility, kernel range, time step in relation to the friction factor, and reduced external pressure, the latter being relevant under isobaric conditions. Appropriate parametrization yields precise thermodynamics, further improved through a simple renormalization protocol. The dynamics can be restored exactly through a trivial manipulation of the time step and friction coefficient. A semiempirical formula for the upper bound on the time step is derived, which takes into account variations in compressibility, friction factor, and kernel range. We test the scheme on realistic mesoscopic models of fluids, involving both simple (Helfand) and more sophisticated (Sanchez–Lacombe) equations of state.</p>\",\"PeriodicalId\":60,\"journal\":{\"name\":\"The Journal of Physical Chemistry B\",\"volume\":null,\"pages\":null},\"PeriodicalIF\":2.8000,\"publicationDate\":\"2024-07-10\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"https://pubs.acs.org/doi/epdf/10.1021/acs.jpcb.4c01441\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"The Journal of Physical Chemistry B\",\"FirstCategoryId\":\"1\",\"ListUrlMain\":\"https://pubs.acs.org/doi/10.1021/acs.jpcb.4c01441\",\"RegionNum\":2,\"RegionCategory\":\"化学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q3\",\"JCRName\":\"CHEMISTRY, PHYSICAL\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"The Journal of Physical Chemistry B","FirstCategoryId":"1","ListUrlMain":"https://pubs.acs.org/doi/10.1021/acs.jpcb.4c01441","RegionNum":2,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"CHEMISTRY, PHYSICAL","Score":null,"Total":0}
Development of a Meshless Kernel-Based Scheme for Particle-Field Brownian Dynamics Simulations
We develop a meshless discretization scheme for particle-field Brownian dynamics simulations. The density is assigned on the particle level using a weighting kernel with finite support. The system’s free energy density is derived from an equation of state (EoS) and includes a square gradient term. The numerical stability of the scheme is evaluated in terms of reproducing the thermodynamics (equilibrium density and compressibility) and dynamics (diffusion coefficient) of homogeneous samples. Using a reduced description to simplify our analysis, we find that numerical stability depends strictly on reduced reference compressibility, kernel range, time step in relation to the friction factor, and reduced external pressure, the latter being relevant under isobaric conditions. Appropriate parametrization yields precise thermodynamics, further improved through a simple renormalization protocol. The dynamics can be restored exactly through a trivial manipulation of the time step and friction coefficient. A semiempirical formula for the upper bound on the time step is derived, which takes into account variations in compressibility, friction factor, and kernel range. We test the scheme on realistic mesoscopic models of fluids, involving both simple (Helfand) and more sophisticated (Sanchez–Lacombe) equations of state.
期刊介绍:
An essential criterion for acceptance of research articles in the journal is that they provide new physical insight. Please refer to the New Physical Insights virtual issue on what constitutes new physical insight. Manuscripts that are essentially reporting data or applications of data are, in general, not suitable for publication in JPC B.