科学上网

 首先按照百度上的解决办法将默认的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 单个元素不加括号

  通过在终端中运行，在终端窗口中找到 conda 环境的目录 。 echo   $CONDA_PREFIX 输入该目录并创建这些子目录和文件： cd $CONDA_PREFIX mkdir -p ./etc/conda/activate.d mkdir -p ./etc/conda/deactivate.d touch ./etc/conda/activate.d/env_vars.sh touch ./etc/conda/deactivate.d/env_vars.sh 编辑 ./etc/conda/activate.d/env_vars.sh 如下： #!/bin/sh export LD_LIBRARY_PATH=/public/software/anaconda3/lib:$LD_LIBRARY_PATH 编辑 ./etc/conda/deactivate.d/env_vars.sh 如下： #!/bin/sh unset LD_LIBRARY_PATH

首先我们百度发现这个库是在conda的lib文件夹下 然后我们打开集群公共软件中conda下的lib文件夹有这个so库 然后在环境变量中加入这个路径

 第一种 VASPkit结构微扰 使用4，structure editor  然后使用407Apply Random-Displacement on Selected Atoms  接着使用功能1 原子选择all，最大偏移选择0.03（推荐） 第二种dpdata微扰 from multiprocessing . spawn import import_main_path import dpdata n0 = 25 #微扰结构数量 c0 = 0.05 #单元收到5%的扰动 d0 = 0.5 #原子位置受到0.5埃的扰动 perturbed_system = dpdata .System( './POSCAR' ).perturb( pert_num = n0 ,     cell_pert_fraction = c0 ,     atom_pert_distance = d0 ,     atom_pert_style = 'normal' ) #print(perturbed_system.data) for i in range ( len ( perturbed_system )):     perturbed_system .to( 'vasp/poscar' ,     'POSCAR' + str ( c0 )+ '.' + str ( d0 )+ '.' + str ( i + 1 ), frame_idx = i )    

  关于热导率在在 0 附近震荡的情况，在 github 上已经有了讨论， https://github.com/deepmodeling/deepmd-kit/pull/1093 。 主要原因是 lammps 根据樊老师 15 的 PRB ，更新了计算热流的方法， EMD 的输入文件 compute heat/flux 命令需要做个小的修改。大家可以下载最新版的 deepmd ，然后改下 EMD 输入文件就可以了。应该 EMD 比以前更 physical 了

 

 https://blog.csdn.net/qq_43689832/article/details/124164624?utm_source=app&app_version=5.3.0&utm_source=app

 

  https://github.com/deepmodeling/dpgen/pull/567

 lmp_gpu改成lmp

  解决办法 删除文件中的 # 中文解释

  Lammps 手册 https://docs.lammps.org/fix_ave_time.html lammps 教程 https://docs.lammps.org/fix_ave_time.html lammps python 接口 https://blog.sciencenet.cn/home.php?mod=space&uid=3366501&do=blog&id=1190121 寇享（很多教程 包括分子动力学， deepmd 等） https://www.koushare.com/home/searchRes?key=deepmd b 站 up 天玑算 - 科研服务

  环境变量改法：      # 2022.3.4  An huang export PATH=/public/bin:$PATH export PATH=/public/software/anaconda3/envs/deepmd/bin:$PATH 同时卸载了我自己安装的 conda 同时 python 调用的位置也发生了改变

   以后每次使用 deepmd 前都需要激活这个环境，只使用 vasp 和 lammps 不用激活这个环境 之后只需要 conda activate deepmd 就可以了       PBE54 是赝势 赝势 （ pseudopotential ），或 有效势 （ effective potential ），是指在对 能带结构 进行 数值计算 时所引入的一个虚拟的 势 。引入赝势有助于实现一个复杂的系统的近似计算。事实上，赝势近似法是正交平面波方法（ Orthogonalized Plane Wave method ， OPW method ）的延伸，其应用范围包括 原子物理学 和 中子散射 。 “ 赝势 ” 这个概念是由汉斯 · 赫尔曼于 1934 年首先发表的。 PBS命令与使用 https://blog.csdn.net/educast/article/details/7210453 conda 的安装与使用（ 2021-04-27 更新） - 简书 (jianshu.com) export   PATH=/home/xiehanpeng22/bin/anaconda3/bin:$PATH 进入 conda 环境 chmod 777 activate . ./activate conda activate /public/home/xiehanpeng22/bin/anaconda3 退出环境 . ./deactivate # 或者用 conda deactivate 扩胞的原则是，扩胞之后 3 个方向的晶格差不多一样，介于 10-20 之间