3。弹性常数似乎也没办法计算
弹性常数可以用侯老师的提供的脚本defvector.f手动计算,是基于能量更应变的关系。也可以用 vasp自带的基于应力和应变,基于线性响应的方法来计算的。
第一种基于能量更应变的关系,就是算出来一系列能量,然后拟合就可以知道弹性常数,木虫上有侯老师的教程
第二种基于线性响应的方法,直接使用vasp自带的就可以了,OUTCAR中会直接给出,可以参考链接https://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html
Finally, IBRION=6 and ISIF3 allows to calculate the elastic constants. The elastic tensor is determined by performing six finite distortions of the lattice and deriving the elastic constants from the strain-stress relationship [4]. The elastic tensor is calculated both, for rigid ions, as well, as allowing for relaxation of the ions. The elastic moduli for rigid ions are written after the line ,
加了真空层的。
可以的 但是我只算过能带和 太密度。
能带的k点就在垂直于c方向的平面上取就好
1.K点该怎么选择呢,只能选择x和y方向的吗?
只用选择平面内的就可以,异质结和多层材料是一样的,可以参考附件中的能带
2。用emc计算有效质量总是报错。
这个要贴出来具体错误
3。弹性常数似乎也没办法计算
弹性常数可以用侯老师的提供的脚本defvector.f手动计算,是基于能量更应变的关系。也可以用 vasp自带的基于应力和应变,基于线性响应的方法来计算的。
第一种基于能量更应变的关系,就是算出来一系列能量,然后拟合就可以知道弹性常数,木虫上有侯老师的教程
第二种基于线性响应的方法,直接使用vasp自带的就可以了,OUTCAR中会直接给出,可以参考链接https://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html
Finally, IBRION=6 and ISIF3 allows to calculate the elastic constants. The elastic tensor is determined by performing six finite distortions of the lattice and deriving the elastic constants from the strain-stress relationship [4]. The elastic tensor is calculated both, for rigid ions, as well, as allowing for relaxation of the ions. The elastic moduli for rigid ions are written after the line ,