请问QE中怎样计算反应焓变和吉布斯自由能变？

风飞

老师，您好，
          使用QE软件通过CI-NEB方法可以计算H2O在MgO(100）面的解离能垒（电子能垒），

          请问怎样计算H2O在MgO(100）面解离的的焓变（H）和吉布斯自由能变（G）呢？

卡开发发

用QE配合phonopy进行处理会比较方便，有限差分和线性响应的结果现在基本上都可以处理。需要小心的是，对于气体分子的平动、转动以及体积功等需要自行修正。

风飞

卡开发发 发表于 2021-5-5 10:07
用QE配合phonopy进行处理会比较方便，有限差分和线性响应的结果现在基本上都可以处理。需要小心的是，对于 ...

嗯，好的 谢谢老师，

但是有几个疑问，
1）之前看到他们输算声子特别的费时，并且还需要扩胞，我的体系有80-90个原子甚至120多个原子，如果在经过扩胞，那么体系将会变得特别的大，算声子会不会算不动呢？（电脑配置：双路8171M， 52核）


2）请问怎样对气体分子的平动、转动以及体积功等进行修正呢？

卡开发发

风飞 发表于 2021-5-5 19:49
嗯，好的 谢谢老师，

但是有几个疑问，

1、表面体系可能不需要扩，因为按道理分子应该吸附在足够大的表面。另外QE所使用的DFPT的方法也不必使用扩超晶格的方法来做，对于小的表面或者晶体理论上可以设置多个q点进行处理。

2、Phonopy对气体分子还是当成固体做的，所以缺少了平动转动和体积功，平动和转动的修正可以查统计力学的公式，体积功用pV=nRT处理。你可以用ASE或者Shermo处理，也可以自己写程序处理，解决方法不是唯一的。

3、也有文献固定部分原子（振动时不考虑这些原子的贡献），但这种做法的可靠性有待讨论，这个问题你需要自行考虑。

风飞

卡开发发 发表于 2021-5-6 00:58
1、表面体系可能不需要扩，因为按道理分子应该吸附在足够大的表面。另外QE所使用的DFPT的方法也不必使用 ...

嗯，好的 谢谢老师