VASP计算弹性常数

记录一下最近计算的不同压力下Cr3C2的弹性常数

1. vasp直接计算

vasp5.1版本后支持设置IBRION=6来直接计算材料的弹性常数，最直接的做法就是在INCAR中设置以下三个参数

IBRION = 6
ISIF = 3
NFREE = 4

具体的INCAR设置参见

ISTART=0
ICHARG=2
PREC=Accurate
GGA = PE
ADDGRID =.TRUE.# Electronic relaxation
ENCUT=550
EDIFF=1E-8ISMEAR=0
SIGMA=0.05
NELMIN=5# Ionic relaxation
NSW=1
EDIFFG=-1e-3
ISIF=3IBRION=6   #计算弹性常数
POTIM=0.015
NFREE = 4# Other Tags
PSTRESS=50   #计算在5GPa压强下的弹性常数#Write flags
LWAVE=.FALSE.
LCHARG=.FALSE.KSPACING = 0.25   #设置K点密度

计算结束以后就会在OUTCAR中输出弹性常数（注意单位是kBar，1GPa=10kBar）

Cr3C2属于正交晶系，总共有9个独立矩阵元(C11,C12, C13,C22，C23，C33，C44，C55，C66)

2. vaspkit计算

要计算不同压强下的弹性常数需要现在不同压强下进行结构优化，再分别计算弹性常数。具体计算细节后续参见博客。

计算结果如下

除了弹性常数外，还可以得到柔顺系数，并且杨氏模量，体积模量等都不需要手动去算，都会直接输出。

VASP计算弹性常数相关推荐

VASP计算笔记_弹性常数计算
二阶弹性常数计算主要采用能量–应变方法 AELAS软件计算必要的文件 elas_strain.sh elas_cal.sh POTCAR INPOS 优化后的POSCAR.可在脚本设置或者单独给出 ...
VASP计算光学性质
PS:请见文末的打赏选项. VASP计算光学性质的INCAR是根据HSE06计算的自洽INCAR来改的 SCF INCAR: SYSTEM = System ISTART = 1 ICHARG = 1 ...
VASP计算新手礼包：基于CentOS 7.6 系统安装VASP 5.4.1
2020-9-24更新:更新VASP5.4.1网盘链接: 链接:https://pan.baidu.com/s/1n9C9CmF0j7iMilo1cPe1cA 提取码:69dj 2019-9-9 更新 ...
VASP计算笔记_声子谱计算
声子谱计算主要使用phonopy软件计算声子谱.版本:phonopy-2.1.3 密度泛函微软理论/线性响应方法 (DFPT) 必要的输入文件: INCAR KPOINTS POSCAR-unitc ...
使用VASP计算材料能带和态密度
前言一.前期准备二.计算流程 1.静态自洽(scf) 2.能带计算 3.态密度计算总结前言能带结构是目前采用第一性原理计算所得到的常用信息,可用来结合解释金属.半导体和绝缘体的区别.能带用来 ...
使用vasp计算α-SiO2弹性模量(个人记录)
使用vasp计算α-SiO2弹性模量文件下载由于vasp wiki的example里有alpha-SiO2的样例,故INCAR.KPOINTS.POSCAR.POTCAR文件可以从vasp wik ...
VASP计算AgGaS2能带及态密度及光学性质
第一步:构型优化 1．准备四个输入文件 POSCAR INCAR POTCAR KPOINTS POSCAR: 从ms中导入AgGaS2结构,选择CASTEP,file,save,并保存成原包.这样, ...
vasp计算之输入输出文件
一.vasp文件: INCAR in ** STOPCAR in stout out POTCAR in ** KPOINTS in ** IBZKPT out POSCAR in ** CONTCA ...
vasp计算---出现各种错误情况
1.项目场景: 项目场景(1):VASP运行时,出现[forrtl : severe (174) : SIGSEGV, segmentation fault occurred] 出现问题原因原因:出 ...

VASP计算弹性常数

1. vasp直接计算

2. vaspkit计算

VASP计算弹性常数相关推荐

最新文章

热门文章