| 查看: 1494 | 回复: 9 | |||
[交流]
超胞计算声子的问题
|
|||
|
超胞计算声子。(1)把原胞转化为2*2*1的超胞计算,出现虚频;(2)原胞转化为2*2*1的超胞,并且截断半径设置为4倍,还是出现虚频(3)采用原胞,截断半径设为4倍,计算的结构显示稳定。1和2为什么会出现虚频呢? (1) (2) (3) |
» 猜你喜欢
售SCI文章,我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急
已经有6人回复
各位大神,目前国内有哪些比较好用的逆合成软件?
已经有8人回复
售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
2楼2012-07-20 10:58:16
★ ★ ★ ★ ★
zhuzhudgyg(金币+1): 谢谢参与
sunyang1988: 金币+2, 谢谢交流 2012-07-20 18:21:33
sunyang1988: 专家考核 2012-07-20 18:23:39
zhuzhudgyg: 金币+2 2012-07-21 10:17:05
zhuzhudgyg(金币+1): 谢谢参与
sunyang1988: 金币+2, 谢谢交流 2012-07-20 18:21:33
sunyang1988: 专家考核 2012-07-20 18:23:39
zhuzhudgyg: 金币+2 2012-07-21 10:17:05
|
赞一下计算资源~ 超胞会引起相互作用的复杂化,对收敛的要求有时会更高,以前见过单胞不稳定,搞成超胞后由于限制作用,却声子稳定了的…… 首先尽量提高K点和截断能做好收敛(几何优化时需要的K点和截断要高于收敛条件才比较可靠),然后声子的计算同时也要高的精度,高的收敛标准,因为声子是要求能量的微分的,这说明了声子要远比能量本身敏感,没有好的收敛做出来的声子可靠性很差~ 可以先做一系列的计算让单胞的声子谱收敛了,比如同位置的频率差别在10(cm^-1)之内这样的(这要求不算太高)。然后再考虑跳到超胞~ |
3楼2012-07-20 11:09:54
4楼2012-07-21 11:01:35
★ ★ ★
zhuzhudgyg(金币+1): 谢谢参与
sunyang1988: 金币+1, 谢谢交流 2012-07-21 13:01:07
zhuzhudgyg: 金币+1 2012-07-22 09:52:39
zhuzhudgyg(金币+1): 谢谢参与
sunyang1988: 金币+1, 谢谢交流 2012-07-21 13:01:07
zhuzhudgyg: 金币+1 2012-07-22 09:52:39
|
MS提供的精度基本没什么科学意义,所有参数的收敛要自己测试,更重要的,所有的参数都可以自己手动设置,根本就没有所谓的“精度调到最高”一说。MS这货很大的一个毛病就是提供了几个精度选项,让使用者容易发懒,甚至达不到需要的计算精度。 我不知道你是用的超软还是模守赝势,根据我自己做的收敛测试,在CASTEP中,含C的超软收敛截断能至少是400eV以上,而模守更是要至少800eV以上。其他所有的含原子序数10以内元素的截断能也都只高不低。 并且我说过,声子计算对结构等各种参数非常敏感,所以尽量都在收敛的基础上取得更精确一些才可以。LZ向来不喜欢说计算的体系和结构,但是10*10*2的K点,总让人感觉像是在计算表面而不是晶体。 超胞怎么设置显然没有任何规定,自己多计算多分析结果的异同就是了。最重要的,还是要参数收敛。 |
5楼2012-07-21 11:12:16
6楼2012-07-22 09:51:42
7楼2012-12-11 11:03:38
10楼2012-12-13 09:48:19
简单回复
su-b088楼
2012-12-11 14:15
回复
su-b089楼
2012-12-11 14:15
回复










回复此楼
楼主电脑配置咋样?超胞算声子谱要比原包多很久吧