计算化学公社

标题: quantum espresso计算MoS2晶体的焓值求助 [打印本页]

作者
Author: laifeng20151107 时间: 2021-5-13 10:58
标题: quantum espresso计算MoS2晶体的焓值求助
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指代不同（数值对应不上）吗？

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

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

作者
Author: laifeng20151107 时间: 2021-5-13 16:54

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

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

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

作者
Author: 卡开发发 时间: 2021-5-14 00:04

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)