科学上网

 在计算二维材料单点能计算的时候发现120个电子步后仍无法手链，一番请教老师过后，发现可以尝试修改NELMDL这个参数 有关NELMDL参数的解释 一开始尝试NELMDL=-12和-18，并没帮助收敛（其中-18在108个电子步才收敛），这里说在vasp6以上版本使用ALGO=all和NELMDL=1会帮助很快收敛，然后我尝试修改这两个参数，发现确实可以帮助收敛，40个电子步就可以快速收敛了 （而且只花费了2600s，不到一小时） ，但是我仔细查阅ALGO=all这个参数发现 这个算法似乎对于绝缘体比较好用，但是不知道能否用于半导体，持着谨慎的态度又讲ALGO改回了Normal，发现也能够收敛，不过需要50步左右的电子步，所以我们还是选择这个参数会比较保险一点。

.meam/c设置方法 把library.meam和ZrCu.meam文件保存到in文件同一文件夹。 假设in文件中只有两种原子，原子类型1为Zr，2为Cu，势函数设置为： pair_style meam / c pair_coeff * * library . meam Zr Cu ZrCu . meam Zr Cu 假设原子类型1为Cu，2为Zr，势函数设置为： pair_style meam / c pair_coeff * * library . meam Zr Cu ZrCu . meam Cu Zr 假设有四种原子类型，1和2为Cu，3和4为Zr： pair_style meam / c pair_coeff * * library . meam Zr Cu ZrCu . meam Cu Cu Zr Zr 小结：通过上面例子可以看出，library.meam文件名后紧跟的参数Zr Cu功能是从library.meam中提取Zr和Cu的参数，并将Zr和Cu的索引设为1和2. ZrCu.meam文件名的参数是将Zr和Cu分别映射到in文件中的原子类型。 在多体势函数使用的时候  需要把library中所有的元素都写上 pair_style   meam/c   #势函数设置 pair_coeff   * * library.meam AlS SiS MgS CuS FeS SiCu.meam SiS CuS

 首先按照百度上的解决办法将默认的AMIX=0.4改为0.2，但尝试过后仍然报同样的错误 同时也尝试过修改POTIM，减少POTIM也是报这个错误。 vasp的wiki中给出了几种可能： 1、因为给出了不合理的输入几何结构 解决办法：修改输入结构 2、 初始化并行化 ISIF=3 松弛时遇到错误 解决方案：在某些情况下，将 NPAR=4 替换为 NCORE=16 可以解决问题 3、读取 CHGCAR 文件时出现问题（例如抱怨在 k 点网格更改后尺寸已更改） 解决方案：删除 CHGCAR 文件并从头开始密度

 第一种办法  改EDIFF = 1E-8 另一个可能是，POTCAR里面的原子没有按照正确的顺序排列。ENMAX大的应该放在前面，小的应该放在后面。

 如果deepmd版本过低（2.0.3）会出现  what():  incompatable model: version 1.1 in graph, but version 1.0  supported  的报错 如果deepmd版本过高（2.1.4）会出现 ERROR: Unrecognized pair style 'deepmd' (src/force.cpp:279) 的报错 只有当deepmd安装版本为（2.1.0）时可以运行成功（势函数也是通过2.1.0版本获得的） phonolammps in.lammps --dim 2 2 2 --logshow 使用--logshow可以看到报错 之前一直看不到报错也没有输出 导致不知道从哪里下手解决

 和收敛有关的参数主要有能量的CUTOFF，SCF的次数，SMEAR，以及K点选取等，可以尝试一下。 SCF最大仍不收敛 一般来说，MS默认的SCF次数100对于小系统是够用的，但是，原子数一多，就不一定了。 a) 首先结构的合理性，如果自建的结构偏离最低能量太大（或失配度太大），可能会难以收敛。 b) 取消一切对称性，充分驰豫晶格。 c) 有些结构本身就是亚稳态，scf的各项收敛指标如果设的太高，对于亚稳态就可能达不到如此的精度。 d) 改变收敛的条件，降低精度(能量的截断值)。 e) 增加循环次数。 f) 改变赝势。 优化不收敛 增加设Max.Interations的大小 interations 是定义积分精度的，相当于gaussian里的int选项；根据gauss的经验，对重原子如果不用细的积分网格，结果就不准确，特别是频率计算，可能会把正的频率算成负的频率。提高 interations的确可以提高精度，特别是对于过渡态和频率计算。 总能和文献不一致，而且多种方法都不同 很正常，只要相对值近似就可以。 在倒空间进行电子结构计算的方法中,都存在一个能量零点的取法的问题,不同的程序,取的不同。根本原因是V(G)在G=0是发散的(可以去R.T. Martin的电子结构那本书)。 castep是利用周期性,在倒空间进行电子结构计算的.我没有记错的话,Dmol3是在实空间来进行电子结构计算的.能量零点的取法就跟不同前者了。 因此在进行电子结构计算中，不要比较任何绝对能量值的大小，这样是毫无意义的，即使是同一个程序中也是这样，更不用说是两个不同的程序。 能量的相对值才具有物理意义。

 cat /public/software/vasp/POTCAR/PBE/{As,B,H,Ga,N}/POTCAR >|POTCAR 单个元素不加括号