运行前需要准备工作，就是四个文本文件。 这四个文件是告诉软件，你要算什么？你要怎么算？你想得到什么结果？ POSCAR,POTCAR,KPOINTS,INCAR 1） 第一：POSCAR POSCAR（position + car）：结构车厢。 里面放着就是初始结构的位置。 第一行就是名字，无所谓，可以随便取。 第二行是1.0，晶胞的缩放系数，别管那么多，反正1.0就对了。 第三到第五行，晶胞参数，不懂的自己去翻材基课本。 第六行：晶胞内元素种类，这里是Li，示意锂元素。 第七行：原子个数，3，示意三个锂原子。 第八行：不用管，direct就对了。 第九行到第十一行：三个原子的具体坐标。 展示以下具体的结构示意图，可用vesta查看。 2） 第二：POTCAR POTCAR是用来描述原子的物理性质的，比如晶胞内有O原子和Li原子两种原子，他们的性质肯定是不一样的，你总不能只给出一个是O一个是Li后就什么也不管了吧。 还得给出原子的具体信息，如质量，电荷量，电荷分布什么的，这些信息很复杂，你不用管，也没法管，直接去网上下载就行，下载后直接拿来用。 这个文件，也被称为赝（yan）势，赝就假的的意思，势指势场，赝势就是指假的势场，在这里的意思是模拟的势场，它和真实势场存在一定误差。 3） 第三：KPOINTS K点文件，具体涉及布里渊区什么的，入门阶段，你只需要知道你取的点越密，计算精度越高，耗时越长。 入门阶段，你唯一需要改动的就是第四行，就是K点，具体取多少，没有硬性规定，打个比方，你可以取10 10 10，也可以取 5 5 5。但是 10 10 10 肯定比 5 5 5 算的精度要高，但是速度也慢！ 4） 第四：INCAR INCAR里放的是计算参数，比如精度，计算步数。 里面的参数很多，对于初学者，你只需要知道我给出的下面8个，然后再慢慢拓展。 第一行 SYSTEM 就是计算体系的名词，随便取。 第二行 ISIF 就是设置优化的参数，一般取2或者取3。2代表结构优化过程，晶格参数不会变；3代表晶格参数会变。两者有什么区别呢？如果你是在优化某一晶体，你给的初始晶格参数不是最稳定的，你就必须得用3，让他自己找到最稳定的经过参数。但如果你是在计算某一气体分子，那你就可以取2。 第三行 NSW 最大优化步数，在结构优化中，每移动一次原子我们称为一步。我这边设置NSW=500，意味着如果原子移动了500步，结构还没有达到我给的稳定条件，那么也停止计算。一般500步还不收敛，那结构可能就有问题了，也没必要继续计算。 第四行 IBRION 就是原子移动的算法，一般取2就行，如果是分子动力学，则需要取别的值。 第五行 POTIM 原子移动步长。尽管我们不知道具体的算法，但是我们应该意识到，移动的dx和受力F应该是有个递增关系的：dx=k*F，k是个系数。如果k越小，每步原子移动的越慢。这个POTIM就是和K有关，POTIM越大K越大。 第六行 ENCUT 截断能，具体意思自己百度，对于初学者，你只需要知道截断能越大，计算越准确，耗时越长。一般520满足绝大多数体系的计算。 第七行 EDIFF 能量收敛标准，当相邻两步的能量差小于设定值，我们就认为能量稳定，这里取得是0.1E-04 eV。 第八行 EDIFFG 力的收敛标准，当体系中任意两个原子的相互作用力的最大值小于设定值，我们就认为结构稳定，这里取得是 0.02 eV/A。 除了这些参数，还有ISPIN，ISTART，ISYM，ISMEAR等参数需要初学者知道，不过你可以先用之前的八个参数，完成简单的计算。 |
说点什么...