VASP遇到一些问题总结Word格式.docx
《VASP遇到一些问题总结Word格式.docx》由会员分享，可在线阅读，更多相关《VASP遇到一些问题总结Word格式.docx（15页珍藏版）》请在冰豆网上搜索。
（6）关于几何优化
包括很多种了，比如晶格常数和原子位置同时优化，只优化原子位置，只优化晶格常数，还有晶格常数和原子位置分开优化等等。
在PRL一篇文章中见到过只优化原子位置，晶格常数用实验值的例子（PRL100,186402（2008））；
也见到过晶格常数先优化，之后固定晶格常数优化原子位置的情况；
更多的情况则是Fullgeometryoptimization。
一般情况下，也有不优化几何结构直接计算电子结构的，但是对于缺陷形成能的计算则往往要优化。
（7）关于软件
软件大致分为基于平面波的软件，如CASTEP、PWSCF和ABINIT等等，计算量大概和体系原子数目的三次方相关；
还有基于原子轨道线性组合的软件（LCAO），比如openmx，siesta，dmol等，计算量和体系原子数目相关，一般可模拟较多原子数目的体系。
VASP是使用赝势和平面波基组，进行从头量子力学分子动力学计算的软件包，它基于CASTEP1989版开发。
VAMP/VASP中的方法基于有限温度下的局域密度近似（用自由能作为变量）以及对每一MD步骤用有效矩阵对角方案和有效Pulay混合求解瞬时电子基态。
这些技术可以避免原始的Car-Parrinello方法存在的一切问题，而后者是基于电子、离子运动方程同时积分的方法。
离子和电子的相互作用超缓Vanderbilt赝势（US-PP）或投影扩充波（PAW）方法描述。
两种技术都可以相当程度地减少过渡金属或第一行元素的每个原子所必需的平面波数量。
力与张量可以用VAMP/VASP很容易地计算，用于把原子衰减到其瞬时基态中。
02、VASP程序的亮点：
1.VASP使用PAW方法或超软赝势，因此基组尺寸非常小，描述体材料一般需要每原子不超过100个平面波，大多数情况下甚至每原子50个平面波就能得到可靠结果。
2.在平面波程序中，某些部分代码的执行是三次标度。
在VASP中，三次标度部分的前因子足可忽略，导致关于体系尺寸的高效标度。
因此可以在实空间求解势的非局域贡献，并使正交化的次数最少。
当体系具有大约2000个电子能带时，三次标度部分与其它部分可比，因此VASP可用于直到4000个价电子的体系。
3.VASP使用传统的自洽场循环计算电子基态。
这一方案与数值方法组合会实现有效、稳定、快速的Kohn-Sham方程自洽求解方案。
程序使用的迭代矩阵对角化方案（RMM-DISS和分块Davidson）可能是目前最快的方案。
4.VASP包含全功能的对称性代码，可以自动确定任意构型的对称性。
5.对称性代码还用于设定Monkhorst-Pack特殊点，可以有效计算体材料和对称的团簇。
Brillouin区的积分使用模糊方法或四面体方法。
四面体方法可以用Bl&
ouml;
chl校正去掉线性四面体方法的二次误差，实现更快的k点收敛速度。
03、VASP的新功能：
1.大规模并行计算需要较少的内存。
2.加入新的梯度校正泛函AM05和PBEsol；
用标准PBEPOTCAR文件提供新泛函；
改善了单中心处理。
3.离子位置和格矢中加入有限差分，从而得到二阶导，用于计算原子间力常数和声子（需要超晶胞近似），和弹性常数。
计算中自动考虑对称性。
4.离子位置和静电场中加入线性响应，从而得到二阶导，用于计算原子间力常数和声子（需要超晶胞近似），Born有效电荷张量，静态介电张量（电子和离子贡献），内应变张量，压电张量（电子和离子贡献）。
线性响应只能用于局域和半局域泛函。
5.精确的非局域交换和杂化泛函：
Hartree-Fock方法；
杂化泛函，特别是PBE0和HSE06；
屏蔽交换；
（实验性的）简单模型势GW-COHSEX，用于经验的屏蔽交换内核；
（实验性的）杂化泛函B3LYP。
6.通过本征态求和计算含频介电张量：
使用粒子无关近似，或通过GW的随机相近似。
可用于局域，半局域，杂化泛函，屏蔽交换，和Hartree-Fock。
7.完全含频GW，速度达到等离子极点模型：
单发G0W0；
在G和W中迭代本征矢直至自洽；
（实验性的）迭代G（也可以选W）本征矢的自洽GW；
（实验性的）对相关能使用RPA近似的GW总能量；
用LDA计算G和W的顶点校正（局域场效应），仅能用于非自旋极化的情况；
（实验性的）W的多体顶点校正，仅能用于非自旋极化的情况。
8.实验性的功能：
用TD-HF和TD-杂化泛函求解Cassida方程（仅能用于非自旋极化的Tamm-Dancoff近似）；
GW顶点的Bethe-Salpeter（仅能用于非自旋极化的Tamm-Dancoff近似）。
1、VASP能够进行哪些过程的计算怎样设置
我们平时最常用的研究方法是做单点能计算，结构优化、从头计算的分子动力学和电子结构相关性质的计算。
一般我们的研究可以按照这样的过程来进行
如果要研究一个体系的最优化构型问题可以首先进行结构弛豫优化，然后对优化后的结构进行性质计算或者单点能计算。
如果要研究一个体系的热力学变化过程可以首先进行分子动力学过程模拟，然后在某个温度或压强下进行性质计算或者单点能计算。
如果要研究一个体系的热力学结构变化可以首先在初始温度下进行NVT计算，然后进行分子动力学退火，然后在结束温度下进行性质计算研究。
2、什么是单点能计算（singlepointenergy）如何计算
跟其它软件类似，VASP具有单点能计算的功能。
也就是说，对一个给定的固定不变的结构（包括原子、分子、表面或体材料）能够计算其总能，即静态计算功能。
单点能计算需要的参数最少，最多只要在KPOINTS文件中设置一下合适的K点或者在INCAR文件中给定一个截断能ENCUT就可以了。
还有一个参数就是电子步的收敛标准的设置EDIFF，默认值为EDIFF=1E-4，一般不需要修改这个值。
具体来说要计算单点能，只要在INCAR中设置IBRION=-1也就是让离子不移动就可以了。
3、什么是结构优化（structureoptimization）如何计算
结构优化又叫结构弛豫（structurerelax），是指通过对体系的坐标进行调整，使得其能量或内力达到最小的过程，与动力学退火不同，它是一种在０Ｋ下用原子间静力进行优化的方法。
可以认为结构优化后的结构是相对稳定的基态结构，能够在实验之中获得的几率要大些（当然这只是理论计算的结果，必须由实验来验证）。
一般要做弛豫计算，需要设置弛豫收敛标准，也就是告诉系统收敛达成的判据（convergencebreakcondition），当系统检测到能量变化减小到一个确定值时例如EDIFFG=1E-3时视为收敛中断计算，移动离子位置尝试进行下一步计算。
EDIFFG这个值可以为负，例如EDIFFG=，这时的收敛标准是当系统发现所有离子间作用力都小于给定的数值，如０.０２eV/A时视为收敛而中断。
弛豫计算主要有两种方式：
准牛顿方法（quasi-NewtonRMM-DIIS）和共轭梯度法（CG）两种。
准牛顿方法计算速度较快，适合于初始结构与平衡结构（势能面上全局最小值）比较接近的情况，而CG方法慢一些，找到全局最小的可能性也要大一些。
选择方法为IBRION=1时为准牛顿方法而IBRION=2时为CG方法。
具体来说要做弛豫计算，设置IBRION=1或者2就可以了，其它参数根据需要来设置。
NSW是进行弛豫的最大步数，例如设置NSW=100，当计算在100步之内达到收敛时计算自动中断，而100步内没有达到收敛的话系统将在第100步后强制中止（平常计算步数不会超过100步，超过100步可能是计算的体系出了问题）。
参数通常可以从文献中发现，例如收敛标准EDIFFG等。
有的时候我们需要一些带限制条件的弛豫计算，例如冻结部分原子、限制自旋的计算等等。
冻结部分原子可以在POSCAR文件中设置selectivedynamic来实现。
自旋多重度限制可以在INCAR中以NUPDOWN选项来设置。
另外ISIF选项可以控制弛豫时的晶胞变化情况，例如晶胞的形状和体积等。
费米面附近能级电子分布的smearing是一种促进收敛的有效方法，可能产生物理意义不明确的分数占据态情况，不过问题不大。
在INCAR文件中以ISMEAR来设置。
一般来说K点只有一两个的时候采用ISMEAR=0，金属体材料用ISMEAR=1或2,半导体材料用ISMEAR=-5等等。
不过有时电子步收敛速度依然很慢，还需要设置一些算法控制选项，例如设置ALGO=Very_Fast，减小真空层厚度，减少K点数目等。
弛豫是一种非常有效的分析计算手段，虽然是静力学计算但是往往获得一些动力学得不到的结果。
4、vasp的分子动力学模拟
vasp做分子动力学的好处，由于vasp是近些年开发的比较成熟的软件，在做电子scf速度方面有较好的优势。
缺点：
可选系综太少。
尽管如此，对于大多数有关分子动力学的任务还是可以胜任的。
主要使用的系综是NVT和NVE。
一般做分子动力学的时候都需要较多原子，一般都超过100个。
当原子数多的时候，k点实际就需要较少了。
有的时候用一个k点就行，不过这都需要严格的测试。
通常超过200个原子的时候，用一个k点，即Gamma点就可以了。
INCAR：
EDIFF
一般来说，用1E-4或者1E-5都可以，这个参数只是对第一个离子步的自洽影响大一些，对于长时间的分子动力学的模拟，精度小一点也无所谓，但不能太小。
IBRION=0分子动力学模拟
IALGO=48一般用48，对于原子数较多，这个优化方式较好。
NSW=1000
多少个时间步长。
POTIM=3
时间步长,单位fs,通常1到3.
ISIF=2
计算外界的压力.
NBLOCK=1
多少个时间步长，写一次CONTCAR，CHG和CHGCAR，PCDAT.
KBLOCK=50NBLOCK*KBLOCK个步长写一次XDATCAR.（个离子步写一次PCDAT.）
ISMEAR=-1
费米迪拉克分布.
SIGMA=单位:
电子伏
NELMIN=8
一般用6到8,最小的电子scf数.太少的话,收敛的不好.
LREAL=A
APACO=10径向分布函数距离,单位是埃.
NPACO=200
径向分布函数插的点数.
LCHARG=F尽量不写电荷密度,否则CHG文件太大.
TEBEG=300
初始温度.
TEEND=300终态温度。
不设的话，等于TEBEG.
SMASS=-3
NVEensemble;
-1用来做模拟退火。
大于0NVT系综。
正确：
SMASS=1,2,3是没有区别的。
都是NVTensemble。
SMASS只要是大于0就是NVT系综。
CONTCAR是每个离子步之后都会写出来的，但是会用新的把老的覆盖
CHG是在每10个离子步写一次，不会覆盖
CHGCAR是在任务正常结束之后才写的。
5、收敛判据的选择
结构弛豫的判据一般有两中选择：
能量和力。
这两者是相关的，理想情况下，能量收敛到基态，力也应该是收敛到平衡态的。
但是数值计算过程上的差异导致以二者为判据的收敛速度差异很大，力收敛速度绝大部分情况下都慢于能量收敛速度。
这是因为力的计算是在能量的基础上进行的，能量对坐标的一阶导数得到力。
计算量的增大和误差的传递导致力收敛慢。
到底是以能量为收敛判据，还是以力为收敛判据呢关心能量的人，觉得以能量为判据就够了；
关心力相关量的人，没有选择，只能用力作为收敛标准。
对于超胞体系的结构优化，文献大部分采用Gamma点做单点优化。
这个时候即使采用力为判据（EDIFFG=），在做静态自洽计算能量的时候，会发现，原本已经收敛得好好的力在不少敏感位置还是超过了结构优化时设置的标准。
这个时候，是不是该怀疑对超胞仅做Gamma点结构优化的合理性呢是不是要提高K点密度再做结构优化呢。
在我看来，这取决于所研究的问题的复杂程度。
我们的计算从原胞开始，到超胞，到掺杂结构，到吸附结构，到反应和解离。
每一步都在增加复杂程度。
结构优化终点与初始结构是有关的，如果遇到对初始结构敏感的优化，那就头疼了。
而且，还要注意到，催化反应不仅与原子本身及其化学环境有关，还会与几何构型有关。
气固催化反应过程是电子的传递过程，也是分子拆分与重新组合的过程。
如果优化终点的构型不同，可能会导致化学反应的途径上的差异。
仅从这一点来看，第一性原理计算的复杂性，结果上的合理性判断都不是手册上写的那么简单。
对于涉及构型敏感性的结构优化过程，我觉得，以力作为收敛判据更合适。
而且需要在Gamma点优化的基础上再提高K点密度继续优化，直到静态自洽计算时力达到收敛标准的。
6、结构优化参数设置
结构优化，或者叫弛豫，是后续计算的基础。
其收敛性受两个主要因素影响：
初始结构的合理性和弛豫参数的设置
初始结构
初始结构包括原子堆积方式，和自旋、磁性、电荷、偶极等具有明确物理意义的模型相关参数。
比如掺杂，表面吸附，空位等结构，初始原子的距离，角度等的设置需要有一定的经验积累。
DFT计算短程强相互作用（相对于范德华力），如果初始距离设置过远（如超过4埃），则明显导致收敛很慢甚至得到不合理的结果。
比较好的设置方法可以参照键长。
比如CO在O顶位的吸附，可以参照CO2中C-O键长来设置（如增长20%）。
也可以参照文献。
记住一些常见键长，典型晶体中原子间距离等参数，有助于提高初始结构设置的合理性。
实在不行，可以先在小体系上测试，然后再放到大体系中算。
弛豫参数
弛豫参数对收敛速度影响很大，这一点在计算工作没有全部铺开时可能不会觉察到有什么不妥，反正就给NSW设置个“无穷大”的数，最后总会有结果的。
但是，时间是宝贵的，恰当的设置3小时就收敛的结果，不恰当的设置可能要一个白天加一个黑夜。
如果你赶文章或者赶着毕业，你就知道这意味这什么。
结构优化分电子迭代和离子弛豫两个嵌套的过程。
电子迭代自洽的速度，有四个响很大的因素：
初始结构的合理性，k点密度，是否考虑自旋和高斯展宽（SIGMA）；
离子弛豫的收敛速度，有三个很大的影响因素：
弛豫方法（IBRION）,步长（POTIM）和收敛判据（EDIFFG）.
一般来说，针对理论催化的计算，初始结构都是不太合理的。
因此一开始采用很粗糙的优化（EDIFF=，EDIFFG=），很低的K点密度（Gamma），不考虑自旋就可以了，这样NSW<
60的设置就比较好。
其它参数可以默认。
经过第一轮优化，就可以进入下一步细致的优化了。
就我的经验，EDIFF=1E-4,EDIFFG=，不考虑自旋，IBRION=2，其它默认，NSW=100;
跑完后可以设置IBRION=1，减小OPTIM（默认为，可以设置）继续优化。
优化的时候让它自己闷头跑是不对的，经常看看中间过程，根据情况调节优化参数是可以很好的提高优化速度。
这个时候，提交两个以上的任务排队是好的方式，一个在调整的时候，下一个可以接着运行，不会因为停下当前任务导致机器空闲。
无论结构优化还是静态自洽，电子步的收敛也常常让新手头痛。
如果电子步不能在40步内收敛，要么是参数设置的问题，要么是初始模型太糟糕（糟糕的不是一点点）。
静态自洽过程电子步不收敛一般是参数设置有问题。
这个时候，改变迭代算法（ALGO），提高高斯展宽（SIGMA增加）,设置自洽延迟（NELMDL）都是不错的方法。
对于大体系比较难收敛的话，可以先调节AMIN,BMIX跑十多步，得到电荷密度和波函数，再重新计算。
实在没办法了，可以先放任它跑40步，没有收敛的迹象的话，停下来，得到电荷密度和波函数后重新计算。
一般都能在40步内收敛。
对于离子弛豫过程，不调节关系也不大。
开始两个离子步可能要跑满60步（默认的），后面就会越来越快了。
总的说来，一般入门者，多看手册，多想多理解，多上机实践总结，比较容易提高到一个熟练操作工的水平。
如果要想做到“精确打击”，做到能在问题始发的时候就立刻采取有效措施来解决，就需要回归基础理论和计算方法上来了。
7、优化结果对初始结构和“优化路径”的依赖
原子吸附问题不大，但是小分子吸附，存在初始构型上的差异。
slab上水平放置，还是垂直放置，可能导致收敛结果上的差异。
根据H-K理论，理想情况下，优化得到的应该是全局最小，但在数值计算的时候可能经常碰到不是全局最小的情况。
实际操作中发现，多个不同初始结构优化收敛后在能量和结构上存在一定差异。
为了加快收敛速度，特别是对于表面-分子吸附结构，初始放松约束，比如EDIFF=1E-3，EDIFFG=,NSW=30可能是很好的设置。
但是下面的情况应当慎重：
EDIFF=1E-3;
EDIFFG=;
！
或者更小
NSW=500；
或者更大
电子步收敛约束较小，而离子步约束偏大，离子步数又很多，这种情况下，可能导致的结果是结构弛豫到严重未知的区间。
再在这个基础上提高约束来优化，可能就是徒劳的了——结果不可逆转的偏向不正常的区间。
好的做法，是对初始结构做比较松弛的约束，弛豫离子步NSW应该限制在一个较小的数值内。
EDIFF=1E-3的话，EDIFFG也最好是偏大一些，如而不是.这样可以在较少的步数内达到初步收敛。
对于远离基态的初始结构，一开始在非常松弛的约束下跑若干离子步，时间上带来的好处是很大的。
对于100个原子的体系用vasp做Gamma点优化，如果一开始就是正常优化（EDIFF=1E-4,EDIFFG=）设置，开始十个离子步可能都要花上几个小时。
如果这个时候才发现输入文件有错误，那下午的时间就白费了，顺便带上晚上机器空转。
所以，我习惯的做法，是在初始几步优化后，会用xcrysden检查一下XDATCAR中的数据，用生成，然后看看弛豫过程是不是按照设想的那样。
后续过程跑完一个收敛过程，就再检查一下。
如此这般，才放心的展开后续计算。
8、目的导向的结构优化
结构优化到这个阶段，是高级的了。
为了得到特定结构，或者为了验证某些猜想，需要设计合理的初始结构，然后在这个基础上小心优化，比如POTIM=跑几步看看，然后修改优化参数。
我遇到过的一件跟结构优化关系很大的算例是CeO2氧空位结构电子局域的问题bbs/tid=2954558）。
按照一般方式（从优化好的bulk建slab模型，然后优化）得到一个O空位留下的两个电子均匀局域到O次外层三个Ce原子上，得到空位形成能.经高人指点后，调节空位附近O原子位置，打破对称性后重新优化，两个电子完美的局域到两个Ce原子上了。
并且空位形成能降低到eV。
从这个例子可以看到，结构优化存在不少技巧的，这些技巧建立在研究者对模拟对象的物理意义的理解上。
对物理图像的直观深入理解，才能做好模型预设，在此引导下才可能有目的的优化出不比寻常的结果。
目前第一性原理理论中的交换关联泛函部分包含经验参数。
考虑这一点对优化结果的影响也很有意思。
比如有专家提到，DFT+U参数对某些结构的收敛终态构型有影响。
构型的变化可能影响表面反应过程。
基于这一点，一个好的计算研究可能就出来了。
真实过程总是复杂多变的。
无论何种模拟，估计都可以找到一些试验现象来验证。
但是到底应该如何评判模拟结果，如何从第一性原理研究中得出有意义的结论需要很好的洞察力。
这样的模拟不见得就必须建立的试验的基础上，完全凭空设计的模型有可能更能优美的解释本质。
9、在单机上计算态密度好像不会出问题。
我先谈一下我的看法：
第一个WARNING，可以在INCAR文件中设置NGX,NGY和NGZ的值，设置的值要足够大，就可以消除这个warning。
设置多大合适呢这就要用到编译vasp时，同时也编译得到的makeparam小程序，makeparam可以帮助你预先检查你设置的文件是否正确，以及某些参数的值是否合适。
要得到合适的NGX,NGY,NGZ以及NBANDS,先在INCAR中不设置这些参数的值，然后运行makeparam>
，其中是包含了输出结果的文件，在文件中你可以看到这些参数的值，以及计算大概需要多少的内存。
然后把文件中的NGX,NGY,NGZ和NBANDS的值拷贝到INCAR文件中。
第二个是计算态密度时，我个人的做法是，一般把KPOINTS文件中的k点增多，然后把INCAR文件中的ISTART=1,ICHARG=11，当然还设置RWIGS。
最后把静止自洽计算得到的CHG和CHGCAR文件拷贝到当前目录下。
从我在单机上的计算来看，没有WAVECAR文件也是可以计算态密度的。
我想你出现的这个问题，可能是你cluster上计算时，每个节点上的CHGCAR和WAVECAR文件不一致造成的。
第三个是当k点数增加了，会出现一个WARING，要把此WARNING消失掉，在INCAR文件中设置NELMDL，它的值小于等于默认值（默认值好像是-5,你可以设为-6）。
没有cluster的系统用来计算，也没有这样的经历，我仅从在单机上的计算经验来谈，有错还请包涵
10、如何用VASP计算铁磁、反铁磁和顺磁
顺磁，意味进行non-spinpolarized的计算，也就是ISPIN=1。
铁磁，意味进行spin-polarized的计算，ISPIN=2
升级会员 