科学上网

 在计算二维材料单点能计算的时候发现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