求助lammps中如何将MEAM势函数与ZBL势结合(短程作用时使用ZBL，长程时使用MEAM)？

墨茶茶

各位老师好！我最近在做MEAM势函数描述的材料的级联碰撞，我在阅读文献过程中发现有几篇文章提到将MEAM势函数与ZBL势结合，从而能够更好的描述级联碰撞。我在翻阅内外网之后发现在国外论坛有人提到过MEAM势函数本身就自带ZBL，只需要在势函数文件中将带有ZBL字段的关键词enable就行。但我在具体实践过程中发现，当我enab ZBL势之后原先较为正常的级联碰撞过程变得不太正常，PKA原子在模拟过程中几乎不与周围原子发生作用产生缺陷，径直穿透模拟盒子。我自我纠错过程中没有发现问题，因此在这希望能有做过MEAM与ZBL结合进行级联碰撞的老师能够解答一下我的疑惑？

势函数文件中的ZBL关键词
MEAM, zbl=1 -blend the MEAM with the ZBL potential for small atom separations - LAMMPS / LAMMPS Beginners - Materials Science Community Discourse (matsci.org)
国外论文lammps开发团队人员的回答

修改ZBL字段前后级联碰撞产生缺陷对比(左有ZBL，右无。初级级联碰撞原子从右下角向左上角赋予初速度)

墨茶茶

wuke 发表于 2024-6-13 15:49
你好，请问你的解决了嘛，我的在用meam势模拟级联碰撞时也会出现pka射穿现象，该怎么解决这个问题呢？

您好，不好意思，没有看到。这个问题并没有解决。我看其他论文都采用了ZBL势，所以我最后也采用了

墨茶茶

dxdxdx 发表于 2024-2-24 14:29
同学你好 请问你pka射穿的问题解决了吗，是否可以交流讨论下

您好，不好意思，没有看到。这个问题并没有解决。我看其他论文都采用了ZBL势，所以我最后也采用了

wuke

你好，请问你的解决了嘛，我的在用meam势模拟级联碰撞时也会出现pka射穿现象，该怎么解决这个问题呢？

dxdxdx

同学你好 请问你pka射穿的问题解决了吗，是否可以交流讨论下