lammps教程:nve/nvt/npt系综设置方法

系综是指大量性质和结构完全相同的、处于各种运动状态的、各自独立的系统的集合。


在lammps模拟中,常用的系综有nve、nvt、npt三种系综,下面说说这三种系综的使用方法。


1. nve(微正则系综)

nve系综保证系统中原子数量n、体系总体积v和体系总能量e保持不变。


nve系综没有控温的功能,初始条件确定后,在力场的作用下,原子速度发生变化,相应的体系温度发生变化。我们知道,体系总能量e=势能+动能,温度发生变化,动能就会变化,势能和动能相互转换,总能量保持不变。


nve使用方法:

fix 1 all nve

nve语法比较简单,使用fix语句确定原子组就可以了,nve系综对边界条件没有要求。


2. nvt(等温等压)


nvt系综保证体系的原子数量n、体积v和温度t保持不变。


nvt系综下,模拟盒子box的尺寸不会发生变化,lammps通过改变原子的速度对体系的温度进行调节。


nvt语法:

fix ID group-ID nvt temp Tstart Tstop Tdamp

其中:
ID:fix语句ID标识符
group-ID:原子组标识
temp:温度标识
Tstart:初始温度
Tstop:结束温度
Tdamp:阻尼系数
Tstart、Tstop主要设置初始和结束温度,Tdamp称为阻尼系数,主要是控制温度调节的幅度,在nvt系综中,Tdamp一般设置为100*timestep。

#nvt系综下恒温:
fix 1 all nvt temp 300.0 300.0 100
#nvt系综下升温(比如融化):
fix 1 all nvt temp 300.0 2000.0 100
#nvt系综下降温(比如结晶):
fix 1 all nvt temp 2000.0 300.0 100
nvt系综对边界条件没有要求。


3. npt(等压等温)


npt系综保证体系的原子数量n、压强p和温度t保持不变。


npt系综不仅进行控温,还进行控压。和nvt一样,npt系综通过调节原子速度调控温度,不同的是,npt系综下box的尺寸可以发生变化。


npt系综通过改变box的尺寸调节压力,比如,当体系压力超过设定值时,扩大box尺寸降低压力。


npt语法:

fix ID group-ID npt temp Tstart Tstop Tdamp Pstart Pstop Pdamp

npt语法和nvt差不多,只不过多了Pstart Pstop Pdamp三个参数进行压力控制。


Pstart、Pstop分别为初始压力、结束压力,Pdamp压力调节阻尼系数,通常为1000*timestep。


需要注意的是,npt在控压方向上必须保证边界条件为p,否则会出错。

#温度300k,iso表示在各个方向上进行控压,压力为0
fix 1 all npt temp 300.0 300.0 100.0 iso 0.0 0.0 1000.0
#温度300k,对x、y方向进行控压,z方向不控压
fix 1 all npt temp 300.0 300.0 100.0 x 0.0 0.0 1000 y 0.0 0.0 1000

4. 选择依据

在模拟的时候到底应该选择哪一个系综,可能是初学者最头疼的事了。


选择系综主要与模拟的体系、模拟的目的有关。


nve系综是一个封闭的系统,与外界没有能量交换,一般用于不需要控温、体系内能量相互转换的模拟,比如石墨烯卷曲到纳米棒。


在大多数模拟中,用nvt或者npt都可以。


以拉伸为例,使用deform进行拉伸,系综可以为nvt,也可以为npt,不过在npt下更有容易模拟泊松效应。


使用velocity方式进行拉伸,在拉伸方向只能为s边界,使用nvt系综更方便。


nvt和npt系综也可以配合使用,比如,对于比较复杂的模型,可以先进行nvt系综下弛豫,然后再进行npt系综下的弛豫。


系综的选择没有统一的规则,最好的方法是多读文献,相同的模拟环境下,看看文献中用的是哪种系综。

--------------------------

更多lammps教程请关注微信公众号:lammps加油站

lammps 交流咨询,请加QQ群:

 

 

  • 39
    点赞
  • 121
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

lammps加油站

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值