GrandCanonical:适用于SW和LJ流体的腔偏大规范Monte Carlo模拟代码 源码

腔偏向马尔可夫链蒙特卡罗大正则乐团 根据Mezei的基于空腔的算法,在μVT集成体中进行流体模拟的半交互代码,可改进在高密度系统中的插入。 该环境包括一个以原点为中心的固定大小的模拟框,其中填充了基于系统密度的球形分子。 它们的相互作用是根据精简单元系统中的Lennard-Jones电位或Square-Well电位(附加电位的实现很简单)进行的。 根据即将发生的方案中的概率,通过插入,移除或置换分子来生成新状态: 与官方报告的结果相比,设计的代码具有成功的参数。 这包括一个20σ20σ20σ模拟盒,其中包含2个0.5 V V的10E4弛豫步长和2.5个0.5 V V 10E5的平衡步长,总共25,000个样本。 通过解析参数,执行很简单: julia GrandCanonical.jl -h usage: GrandCanonical.jl [-h] μ T positional