如何用VASP计算分子的键能

1、POSCAR文件为O dimer in a box 1.0 ! universal scaling parameters 8.0 0.0 0.0 ! l锾攒揉敫attice vector a(1)0.0 8.0 0.0 ! lattice vector a(2)0.0 0.0 8.0 ! lattice vector a(3)2 ! number of atomscart ! positions in cartesian coordinates0 0 0 ! first atom0 0 1.22 ! second atom将氧分子放在一个立方体中

2、INCAR文件为SYSTEM = O2 dimer 坡纠课柩in a boxISMEAR = 0! Gaussian sme锾攒揉敫aringISPIN = 2! spin polarized calculationNSW = 5 ! 5 ionic stepsIBRION = 2! use the conjugate gradient algorithm（使用共轭梯度算法）

3、KPOINTS文件为Gamma-point only0Monkhorst Pack1 1 10 0 0

4、我们应该在INCAR文件中选择几何优化。本例中最多算5个离子步。优化中使用了共轭梯度(CG)算法(IBRION=2)。CG算法需要沿着搜索方向最小化优化。这是布伦特算法的一种变体。

5、计算过程如下，其中quantitiy trial-energy change是试验步骤中能量的改变

6、实验步数可以由POTIM参数控制，优化后的模型会存于CONTCAR中，如果要继续优化，可以将CONTCAR复制为POSCAR。