计算化学公社

标题: 求助castep计算的SnO2(110)晶面的体系能量过大 [打印本页]

作者
Author: 573994810 时间: 2020-3-6 10:54
标题: 求助castep计算的SnO2(110)晶面的体系能量过大
请教各位大神：
我建立的SnO2(110)晶面，算气体的吸附能，用CASTEP计算的SnO2(110)晶面的体系能量特别大呢，17000＋eV，2×2的超晶胞，看其他文献里都是1100左右eV的能量，怎么差这么多呢？
另外按照教程重复了经典安全Pd(110)晶面的计算，CASTEP算出来也是11000＋eV，网上教程里Pd(110)晶面的能量才1000eV左右啊。
新手小白，实在困惑，请各位大神指点迷津！

作者
Author: zyj19831206 时间: 2020-3-6 11:39
你的参数选取的是否与文献一样呢？各种参数不同，计算结果不同就很正常。

作者
Author: 卡开发发 时间: 2020-3-6 13:28
可能是赝势方面的问题，绝对能量没什么参考价值，应该比较键能或能级差这些信息。

作者
Author: 573994810 时间: 2020-3-6 15:15

卡开发发发表于 2020-3-6 13:28
可能是赝势方面的问题，绝对能量没什么参考价值，应该比较键能或能级差这些信息。

CASTEP help里的经典案例CO吸附在Pd(110)面，函数用的都是GGA PBE
1.优化bulk Pd quality选medium，勾选optimize cell，优化后　final enthalpy为－3493eV
2.构造优化CO　按help里来做的，优化后能量为－596　eV，这个还基本正确
3.切面，优化Pd(110) cleave plane (110),　top fractional 为0，thickness为1.5，真空层厚度为8埃，然后固定下层原子，只弛豫表面，castep优化，不勾选optimize cell，截断能为300eV, k points 3*4*1, 优化后能量为-10477.5eV
4. co吸附到Pd(110)表面，然后几何优化，优化后能量为－11075.3eV
5.另外还算了CO吸附到（1×2×1）超晶胞Pd(110)晶面的几何优化，优化后能量为－11075　eV, 跟吸附在单胞上的一样呢？
按照我的计算结果　吸附能＝CO吸附在单胞Pd(110)-Pd(110)-CO=2eV，跟网络上相关教程的计算结果（1.8eV）近似，这样以单胞吸附的来算是不是也正确？　网络教程的计算结果（CO -598.18eV; Pd(110) -1594.53eV; (1*1)CO on Pd(110) -2185.38eV; (2*1)CO on Pd(110) -4370.99eV）.
但是超晶胞计算的结果就不对了，仔细对照网上的教程看了看，应该是因为超晶胞吸附的键长没设置好。超晶胞吸附不是把吸附分子直接拖进真空层就可以吗？　我的计算结果差将近10倍呢，这样也是对的吗？
谢谢！

作者
Author: 卡开发发 时间: 2020-3-6 18:29

573994810 发表于 2020-3-6 15:15
CASTEP help里的经典案例CO吸附在Pd(110)面，函数用的都是GGA PBE
1.优化bulk Pd quality选medium，勾 ...

计算参数没有，也不知道具体是何种结构，光凭盲猜是没办法给你诊断的。

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