LAMMPS处理overlap atoms的几个命

如果我们给随机的初始条件,那么原子重叠是不可避免的。模拟的开始必须把重叠原子push开,否则模拟就无法进行了。这和能量最小化有区别,因为推开不一定能量就最小了。

如何把原子推开呢,LAMMPS里大概有如下几种方面:

pair_style soft command 加上 small timestep, soft style 势使系统即使原子坐标重叠力也不会无穷大,如果加上小的时间步,这样原子一次移动的距离就不会太大【最好<0.1个单位】,从而保证原子的 blowing up.

fix nve/limit command 是通过限制最大移动距离来实现原子的分离,这个对势能和时间步长没有要求。

minimize command 通过一定的算法使系统能量最小化,再一定时候可以用来push overlap的原子,一般用于周期性系统,而且系统中没有需要固定不动的原子,而上面两个命令没有这个限制。

目前来说,最方便的是第二个命令了,对系统没有特别的要求

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值