首页 > 技术文章 > MD笔记

sysu 2019-06-19 21:16 原文

1.力场中的例子电荷是有效电荷(clayff),有别于化学式中的电荷。

2.游离状态的阳离子(如层间阳离子)的电荷不能变动;而Al-O八面体、Si-O四面体中的离子(Al、Si等)电荷可以微调。

3.minimize是在绝对0度、动能为0的条件下进行结构优化(分子静力学),只对从数据库中得来的模型进行minimize,加水后的模型不需要优化。

先对最小超晶胞(如Na-mmt的超小晶胞为4a*b*c)进行minimize,然后再进行下一步建模(如加水分子或者继续supercell),最后在进行优化并计算。

优化:

nve:保持能量不变;

nvt:提到常温来优化(298或者300k);

npt:提到常压来优化(1个标准大气压)。

4.警告可以忽略,错误需要修改?

(待确认,如模型总电荷也显示warning,但minimize时,总势能会很大(正值1e+16))

如:以下的warning能忽略吗?因为在minimize和优化的过程当中都会出现这个警告,这个是data文件的问题(之前请教过一个老师)。

5.元胞盒子大小从数据库中得到的a、b及3个方向的角度不需要修改,根据自己需要来改c(层间距)。

6.能量最小化minimize如何判断最小化已经成功?

(1)log文件上的minimize的收敛情况;

(2)结构是否稳定(VMD中查看),目前的问题是minimize后原子会跑出盒子外边去?

 

7.优化过程中,timestep、run总步数,这个没有很严格的规定,看能量收敛和结构稳定与否来判断即可。

8.MS——data——lammps中minimize后所得的data——VMD——pdb——MS(主要是坐标信息)——建模(加水或者supercell)。。。

问题1:转换过程中只有坐标信息,没有盒子和键、角信息,不能够识别各个原子类型等??

问题2:pdb导入MS后,由于不能识别,故无法继续建模?

(1)lammps所得的data

     

 

(2)VMD生成的pdb

(3)MS(主要是坐标信息)

 

 

9.模型的总电荷一定要为0,影响很大。

10.Al-oh-ho这里的oh-ho指的是体系中的结合水,层间水是自由水。

11.NPT中的P指的是静水压力,随着P增大(深度加大),c层间距减小;

12.timestep:从1.0fs——0.5fs——0.1fs;

13.原子丢失的原因(lost atoms:.......),Temp、Epair等增大

(1)timestep有可能太大了,导致数值失稳,改小一点;

(2)输入的结构不好

14.计算中Dangerous builds不等于0:

(1)说明粒子在跑的时候“Rcut skin ”的范围内已经跑出去了,或者外边的粒子跑进来了,因此,neighbor list要建得频繁一点才可以;

(2)但是不影响数据收集;

(3)如果是小概率事件没有关系,比如,neighbor list运行1万次,而Dangerous只有几十或更小;

(4)措施:

a、把时间步长timestep调小(推荐做法);

b、把skin变大(不推荐,会影响计算效率)。

 

推荐阅读