前面我们介绍完了算表面吸附，以及过渡态的一些基本的操作，和注意事项。当我们面对一个新的课题时，就需要运用所学到的这些技能来完成所必需的计算，来验证我们的想法，思路等。后面几节，主要参考本人2014年发表的一篇关于甲醇分解反应计算的文章来介绍一下怎么运用所学习到的这些本领。（https://pubs.acs.org/doi/abs/10.1021/cs501698w） 搬砖这个词用来形容我们完成一个 ...

伴随着Python2的落幕，Python3的兴起，很多基于2的软件也面临着着要么被遗忘，要么主动拥抱Python3的选择。其中就有我们最喜欢的p4vasp。https://github.com/orest-d/p4vasp 或者www.p4vasp.at。本人不认识作者，也不知道啥时候出新的版本。使用p4vasp需要有python2，但其他软件又需要python3。为此我们不得不做出一些选择， 最 ...

前面一节，我们讲解了对于扩散这一类反应中，nebmake.pl这个命令生成IMAGES时候的一个坑。通过学习，要了解到一些脚本或者程序本身存在一些缺陷，我们在使用的时候，要避免盲目相信，直接把脚本的结构直接拿来用。总之，原则就是尽可能获取一些具有明确物理化学意义的，比较理想的初始构型。本节，讲解一下旋转过程中，nebmake.pl的坑。这里与其说是坑，不如说是nebmake.pl不适用的情况。因为 ...

汇集本人日常学习中的一些数据库，个人感觉不错的网站。也会不断地更新完善。如果你有好的数据库，链接，推荐的书籍，也欢迎发邮件至lqcata@gmail.com留言帮助我们更新。 结构查找The Materials Project https://www.materialsproject.org/ 使用方法： 1） 点击右上方：Login 2） 点击： Sign in with your ema ...

<<<<<<< HEAD
最近整理材料，发现遗落在角落里的一个python脚本（cssm.py），可以批量切常见金属稳定表面的slab模型。 cssm.py 是：cleave_stable_surfaces_from_metals的缩写。废话不多说，下面的是在天河II上运行的实例。

被批评公众号长草了，简单介绍一下ASE转换文件的小技巧，顺便锄下草。顺利学会下面操作的前提: 1) 已经在自己电脑上安装好了ASE。 2) 电脑里面有POSCARtoolkit.py脚本。见（https://www.bigbrosci.com/2018/11/11/ex73/） 终端操作下面是终端里面转换格式的懒人命令： 1ase gui file_from -o POSCAR 1 ase g ...

算过渡态，光知道检查虚频是不对的。如果算出来虚频多，很大程度是因为你的NEB初始结构有问题，也就是说你的过渡态路径背后的物理化学意义不是那么地理想。而对于NEB的初始结构，大部分人都是通过VTST的nebmake.pl脚本来实现的。 1nebmake.pl IS FS 8 分析上面这个命令：有4个关键的信息: 脚本：nebmake.pl 初始结构: IS 末态结构: FS 插的点数：8 本节默 ...

前面几节，我们从gamma点到3x3x1计算过渡态。首先复习一下计算的流程： 1） 我们计算的是金属表面上H原子的扩散； 2） 使用gamma点计算的时候，slab固定住了； 3）gamma点计算结束后，检查结果： 看到有一个漂亮的NEB图； 能量变化也是稳稳妥妥滴； 结构查看也没啥异常情况； 4） 在确认第3)步检查OK之后，将gamma点的计算备份；增大K点至3x3x1继续算。 5）计算结 ...

什么是ASE？ASE 是Atomic Simulation Environment的简称。详细的介绍大家可以浏览官网：https://wiki.fysik.dtu.dk/ase/about.html。ASE是一个基于python的程序库，可以做的事情很多，不同程序任务的设置，结构的搭建，结构分析都可以做。 为什么用ASE？ ASE里面包含了很多的功能，其中个人认为最为方便地就是计算模型不同文件类型 ...

前面一节我们通过粗糙的方法演示了一遍怎么在超算上跑过渡态的NEB计算。一般来说，通过这一步我们可以大体上确定这个基元反应中的路径是不是靠谱的了。不论你得到一个粗糙的NEB能量曲线： 亦或者一个比较精细的： 这都表明，你快要得到这个过渡态结构了。在进行后面操作讲解之前，我们先做2个事情： 1） 由于本人也是老司机了，算过渡态的时候遇到的那些零零碎碎的问题一时也很难记起来，如果你正在算过渡态，并且遇 ...

通过job-ID快速进入计算目录，可以有效地避免狂输cd命令，提高我们的计算效率。本节我们通过在国科智算超算中心的演示，介绍一个Slurm任务管理系统，保存任务ID，计算目录，并快速根据ID进入计算目录的方法。 核心思想：有很多种方法可以实现：通过Job-ID进入计算目录的办法。这里只介绍本节的思路。 通过squeue 命令获取用户所有的任务的ID； 通过scontrol 命令获取所有ID对应的 ...

国科智算提交VASP任务今天QQ群里，有人问怎么在国科智算的超算中心提交VASP的任务。本着一言不合就写教程的态度，这一节我们就看下VASP的任务是怎么提交的。提交任务的脚本见群文件：vasp_qiangli.sh。下载后重命名成vasp.sh即可。想试用的，购买机时的可以加QQ群：608307988咨询一下。 Slurm国科智算上系统采用的是slurm 任务调度工具。Simple Linux U ...

前面的准备工作完成了，INCAR，KPOINTS，POSCAR（IMAGES）， POTCAR也检查完了。剩下的就是准备脚本提交文件了。本节主要是在天河II号超算中心上给大家简单示范一下： 1）准备脚本，提交任务； 2）过渡态任务运行时候的查看； 3） VTST脚本nebresults.pl的安装和使用。 提交任务首先通过命令操作，熟悉下天河II号超算中心提交NEB计算的流程。 1234567 ...

前面一节，我们提到说，在Ubuntu或者其他Linux系统下面，可以使用p4vasp或者ASE将结构批量打开，查看我们初步设置的NEB路径是否合理。但是在Windows系统下，我们不方便使用命令进行查看。这里介绍一下在Windows下面通过VASPkit结合VMD查看NEB路径结构的方法。 软件的获取： VASPkit程序的下载 链接：https://sourceforge.net/project ...

前面我们介绍了怎么准备过渡态计算的INCAR，POSCAR，POTCAR和POTCAR，这一节理所当然，就是要提交任务了。过渡态计算的时候，和大多数的计算是一样的：检查输入文件，准备好脚本，敲命令提交任务。所以，我们就在这三个方面展开，具体讨论一下里面的细节部分，避免不必要的出错，提高计算的成功率和效率。 心态算过渡态的时候，要把自己想象成一个整装待命的狙击手，要熟悉自己的枪（VASP），周围的环 ...

过渡态计算INCAR的设置vasp的过渡态任务所需要的五个文件：INCAR, KPOINTS, POSCAR, POTCAR以及脚本。我们已经学习了三个，还差INCAR和脚本，这一节，我们主要介绍INCAR中参数的设置。很多读者都没有接触过渡态的计算，总是幻想着过渡态是如何地高大上，令人可望而不可即。其实很简单，记住前面我们说的那句话，过渡态的计算无非就是优化一个特殊的结构而已。既然是优化结构，那 ...

前面2节我们讲了如何通过vtst的脚本进行插点,得到过渡态计算中特定数目的images，以及一个非常实用的小脚本来固定slab模型中底层的原子。本节我们就准备一下过渡态计算中：POSCAR、KPOINTS和POTCAR文件的准备工作。 关于Images前面我们插了8个点，得到10个文件夹，分别为00， 01， 02，….09 细心的你会发现，每个文件夹中都会有一个POSCAR文件。00 中的POS ...

在学习下面neb的计算前，思考了很久，感觉有必要先给大家介绍一个非常实用的python小脚本，经历过反复的修改，最终在群友小塔玛斯的努力下，完成了一个稳定的过渡态计算中非常实用小脚本：POSCARtoolkit.py 本节我们介绍一下这个脚本的具体使用方法。 版本和版权问题：版本： 这个脚本有2个版本，分别适用于python3和python2的2.6及以上版本。 大家根据自己使用的python版本 ...

不好意思，让大家久等了。最近一直很忙，后面还是会很忙。为了不让大家失望，抽时间把简单的过渡态计算讲解一下，这一节我们需要学习的是通过CI-NEB（Climbing Image Nudged Elastic Band）计算H原子在Ru(0001)表面上的扩散过程。 NEB计算过渡态的准备工作准备工作1： 说到NEB或者CINEB，不得不提的就是Henkelman课题组的VTST：http://he ...

前面一节我们讲了直接优化NH3翻转的过渡态。这一节我们继续直接优化过渡态的介绍。计算乙烷分子（CH3CH3）绕C—C键旋转的过渡态。我们知道：乙烷有交错式和重叠式构型。而重叠式能量较高。想象乙烷分子绕着C—C旋转，从交错式到重叠式，再到交错式这一个过程。很容易就得到这个结论：重叠式的结构就是过渡态。所以，如果计算这个过渡态的工作就转化成了优化重叠式的构型。这也是本节的一个主要思想。 手动搭建初始 ...