静态计算INCAR
时间: 2023-10-10 19:06:30 浏览: 137
INCAR 是 VASP 软件的一个输入文件,用于控制计算的参数和条件。静态计算通常指的是计算晶体结构的能量、电子密度和电子态密度等物理量。在进行静态计算之前,需要先准备好 INCAR 文件,其中包括定义计算类型、精度、计算算法等参数。下面是一个简单的 INCAR 文件示例:
```
SYSTEM = crystal ! 定义计算体系名称
ISTART = 0 ! 从头开始计算
ICHARG = 2 ! 从 CHGCAR 文件读入电荷密度
ENCUT = 500 ! 截断能
ISMEAR = 0 ! 无平滑
SIGMA = 0.1 ! 能级展宽
EDIFF = 1E-6 ! 能量收敛标准
NELM = 100 ! 最大电子步数
LREAL = .FALSE. ! 关闭实空间投影
```
以上是一个基本的静态计算 INCAR 文件,其中定义了计算体系名称为 "crystal",从头开始计算,从 CHGCAR 文件读入电荷密度,使用截断能为 500 eV,无平滑,能级展宽为 0.1 eV,能量收敛标准为 1E-6 eV,最大电子步数为 100 步,关闭实空间投影。这只是一个示例,具体参数需要根据具体的计算需求进行调整。
相关问题
力常数计算INCAR设置
力常数(force constant)在VASP计算中是非常重要的参数,影响到计算的精度和收敛速度。在INCAR文件中设置力常数可以通过以下步骤完成:
1. 执行第一次计算,其中包含几何优化和静态计算,记录POSCAR文件。
2. 使用VASP软件中的phonopy工具计算力常数。这可以通过在命令行中输入以下命令完成:
```bash
phonopy -f vasprun.xml
```
其中,vasprun.xml是第一步计算得到的VASP输出文件。该命令将生成一个FORCE_SETS文件,其中包含原子振动的力常数信息。
3. 在INCAR文件中添加如下行:
```bash
ISTART = 0
IBRION = -1
NSW = 0
NFREE = 2
FCONST = 0.02
```
其中,FCONST是力常数的值,单位为eV/Å^2。根据第二步计算得到的FORCE_SETS文件中的值来确定FCONST的具体数值。
4. 将第一步计算得到的POSCAR文件和FORCE_SETS文件与INCAR文件一起放入计算的工作目录中,运行VASP计算即可。
需要注意的是,力常数的值应该足够小,以保证计算的精度和收敛速度。如果FCONST值过大,可能会导致计算不稳定或无法收敛。因此,建议在计算之前先进行一些测试,以确定最适合的FCONST值。
MoS2计算弹性常数INCAR设置
在进行MoS2弹性常数计算时,可以在VASP的INCAR文件中设置以下参数:
1. 禁用自旋极化
```
ISPIN=1
```
2. 设定K点网格
```
KPOINTS
0
Monkhorst-Pack
N1 N2 1
0 0 0
```
其中N1和N2是K点网格的数量,可以根据计算精度的要求进行设置。
3. 设定计算类型为静态计算
```
IBRION=-1
ISIF=3
```
4. 设定静态应变
```
ISYM=0
EPSILON
1.01 0.0 0.0
0.0 1.01 0.0
0.0 0.0 1.01
```
其中EPSILON是应变矩阵,这里将每个方向的应变设置为1.01,可以根据需求进行调整。在计算中,需要依次将应变矩阵的每一个元素赋值为1.01,然后将其余元素赋值为0,进行单点计算,得到应变下的能量。
5. 设定计算参数
在计算过程中,需要根据MoS2的大小和对计算精度的要求进行计算参数的设置。以下是一些常用的计算参数:
```
ENCUT: 截断能
EDIFF: 能量收敛精度
EDIFFG: 力收敛精度
ISMEAR: 选择平滑函数类型
SIGMA: 平滑函数宽度
```
以上是MoS2弹性常数计算中常用的一些INCAR参数设置,根据具体情况进行调整。需要注意的是,在计算中应该考虑MoS2的对称性,以减少计算量和提高计算精度。
阅读全文