quantum espresso计算MoS2晶体的焓值求助

laifeng20151107

    quantum espresso中计算类型设置的是vc-relax，输出结果中会显示total energy的值，同时下面会说明the total energy is the sum of the following terms (one-electron contribution, hartree contribution xc contribution, ewald contribution, smearing contrib)，下面同时会显enthalpy old和enthalpy new的值，请问这个焓值可以直接用吗？在.save文件夹中会有一个data-file-schema.xml文件，这个文件中有output栏，里头也会有total energy，但这个total energy中的组成不同，有etot, eband, ehart, vtxc, etxc, ewald, demet，所以请问这个部分的total energy和运行过程中输出显示中的total energy指代不同（数值对应不上）吗？

卡开发发

不能，这个是用来估计变晶格优化的参数。

要得到焓的话要么用MD，要么通过谐振近似来使用声子计算模块计算焓。

laifeng20151107

卡开发发 发表于 2021-5-13 11:12
不能，这个是用来估计变晶格优化的参数。

要得到焓的话要么用MD，要么通过谐振近似来使用声子计算模块计 ...

老师，第二种方法能不能具体一点？是不是先结构优化然后计算声子，跟gaussian里面做振动分析类似？

卡开发发

laifeng20151107 发表于 2021-5-13 16:54
老师，第二种方法能不能具体一点？是不是先结构优化然后计算声子，跟gaussian里面做振动分析类似？

对，基本原理完全相同，建议使用phonopy来做数据的前后处理，目前来看会比较方便。