ex12o2单点计算和优化结果分析.pdf

Ex12 O2 单点计算和优化结果分析 前面几节，我们讲解了 O2 单点计算的输入文件的一些细节部分。这一节我们主要讲解一下 O2 的分子结构分析以及如何初步进行构型优化计算。首先提交 O2 静态计算任务，运行， 等待结束。如果你的任务出错，请跟大师兄的输入文件进行对比、改正，直至计算正常结束。 如下图： 任务结束后查看 OSZICAR 从 OSZICAR 最后，得到体系的磁矩为 2μB，你应该知道这个磁矩是怎么回事，由哪两个电 子贡献。如果不知道，看下图： O2 的分子结构： 看完该图，相比大家心里对于 O2 的成键方式有了一个更加深刻的印象。我们对比下VASP 的输出结果。首先： 分析 α 电子的排布情况： 在这里，你会发现能带 3 和 4 是简并的，应该是 π （p2x）和 π(p2y）轨道中的 α 电子。能 带 5 对应的应该是 σ(2pz) 的电子。 在 O2 分子的电子构型中，两个 O 原子的 2pz 轨道以头碰头的方式形成一个σ键，其能量要 比2px （2py）以肩并肩方式形成的π键能量要低。但是，能带5 的能量（-13.3126）比3 和 4 的（-13.3870）要高些，这与O2 的电子构型不一致，表明VASP 的单点计算结果是不可 靠的。 再看一下 β 电子的排布情况： 能带4 和 5 对应的是应该是 π （2px）和 π （2py ）轨道中的 β 电子。且总的轨道能量与前 面图中一致，这说明 VASP 对β电子的描述是合理的。 为什么出现这样的情况呢？难道跟前面 O 原子的计算一样，VASP 又算不准啦？ 不是的，VASP 怎么着也是个老牌的，响当当的计算程序，总不能让人每天指着鼻子骂算不 准！ 这里的主要原因是：来自于实验的键长值未必就是计算程序所认可的。 也就是说，实验值和理论值之间存在偏差，实验的结构不能直接用来计算其性质，只可以作 为一个理想的初始值。所以，O2 的分子结构需要优化一下。 VASP 优化的时候，需要用到一个参数：IBRION 引用官网的话：IBRIONdetermines how the ions are updated and moved. 也就是说 IBRION 这个参数决定了结构的优化过程。 当你去官网查看的时候（google 搜索 VASP IBRION 这两个关键词），会发现 IBRION 有 很多值。 想要正确进行计算，你就需要去硬着头皮去了解各个值的含义了，这个过程必须自己去做， 只听别人的建议去设置参数，而不自己去主动学习的，你的能力永远不会得到提升！！！ 一般来说，优化结构的时候有 3 个选择： IBRION=3：你的初始结构很差的时候； IBRION=2：共轭梯度算法，很可靠的一个选择，一般来说都用它。 IBRION=1：用于小范围内稳定结构的搜索。 如果你的体系遇到结构不收敛的时候，可以尝试下换下 IBRION 的参数。 下面，我们在 INCAR 中加上 IBRION 参数（IBRION=2），其他输入文件保持不变，重新进 行计算： 你会发现任务很快就算完了。而且只有一步，难道输入的结构就是 VASP 计算出来的的稳 定结构吗？ 有这种可能，但几率极低。 在这个例子里，如果你看下OUTCAR 中的电子构型，发现它的信息和前面的单点计算一样。 这说明，vasp 并没有优化，而是又运行了一次单点计算。 为什么呢？ 这是因为另一个参数：NSW。 NSW 控制几何结构优化的步数。也就是VASP 进行多少离子步。 官网查看下 NSW 选项，发现默认值是 0，也就是没有进行优化。 （默认值，也叫缺省值，英文里面是Default。 意思是，如果你不输入这个参数，程序将 默认使用 XXX 的数值） 现在原因找到了，继续优化吧： 那么 NSW 怎么设置呢？ 首先，它必须是大于等于 0 的整数。 其次，一般来说，简单的体系 200 步内就可以正常结束。不知道什么时候收敛，初始结构 很差，或者设置了很严格的收敛标准，那么你就要增大一下 NSW 了，比如 NSW=500 或者 更大。 因为我们的这个例子很简单，设置了NSW=10 （你也可以设置为100,200 或者 500，不会影 响计算结果的。） 计算完