VASP计算差分电荷和时候怎样改INCAR里的数据

Question

VASP计算差分电荷和时候怎样改INCAR里的数据,计算到vi incar然后就进入了另一个界面也返回不了

Answer 1

引自免费微信小程序：皆我百晓生

在VASP中进行差分电荷密度（Charge Density Difference）的计算时，你需要对INCAR文件中的几个关键参数进行调整。差分电荷密度是两个系统（通常是同一结构但电子态不同的两个系统）电荷密度的差异，通常用于分析化学键合、反应机制等。

首先，确保你的输入文件（CONTCAR或WAVECAR）包含了正确的结构信息。接下来，按照以下步骤修改INCAR文件：

1. 取消注释 ICHARG 参数并设置为合适的值。这将告诉VASP从何处读取或创建波函数。对于差分电荷密度计算，通常使用一个系统的自洽场波函数（SCF）作为起点，因此可以设置为1或2。例如，如果你想要从上一任务的WAVECAR文件继续计算，可以设置为1：

   ICHARG = 1
   

2. 取消注释 NELECT 参数并设置正确的总电子数。对于差分电荷密度计算，这个数值应该与你要比较的两个体系的总电子数之和减去1相等。例如，如果你的两个体系分别有n和m个电子，则总的电子数为n+m-1。

   NELECT = n + m - 1
   

3. 保留 或者适当调整其他参数，如ENCUT、SIGMA、NELM、EDIFF等，以保证电子步的收敛性。

4. 取消注释 LCHARG 和 ADDGRID 参数，它们分别控制是否写入CHGCAR文件（包含电荷密度信息）以及是否增加网格密度以提高电荷密度和势能的计算精度。

   LCHARG = .TRUE.
   ADDGRID = .TRUE.
   

5. 取消注释 LVTOT 和 LVHAR 参数，并设置为.TRUE.，以便同时输出总电势（包括静电势和交换相关势）和离子+Hartree电势到LOCPOT文件。这样你就可以得到差分电势了。

   LVTOT = .TRUE.
   LVHAR = .TRUE.
   

6. 如果你的结构有对称性，确保设置适当的ISYM值。对于差分电荷密度计算，通常建议关闭对称性处理，因为对称可能会干扰结果的直观理解。

   ISYM = 0
   

7. 完成上述更改后，运行VASP。当计算完成后，你会在输出目录下找到CHGCAR和LOCPOT文件。你可以使用一些后处理工具，如XCrySDen、VESTA或其他可视化软件来查看和比较这两个文件，从而得到差分电荷密度和差分势能图。

这是一个例子，展示了如何修改INCAR文件：

ISTART = 1
ISPIN = 1
ICHARG = 1
LREAL = .FALSE.
PREC = Normal
LWAVE = .TRUE.
LCHARG = .TRUE.
ADDGRID = .TRUE.
LVTOT = .TRUE.
LVHAR = .TRUE.
NELECT = n + m - 1
ISYM = 0

# 其他参数保持不变...

请确保替换n和m为你实际的电子数。完成这些设置后，你可以重新运行VASP计算，然后使用后处理工具分析差分电荷密度和势能。