24小时热门版块排行榜    

查看: 1293  |  回复: 9

zhuzhudgyg

金虫 (正式写手)


[交流] 超胞计算声子的问题

超胞计算声子。(1)把原胞转化为2*2*1的超胞计算,出现虚频;(2)原胞转化为2*2*1的超胞,并且截断半径设置为4倍,还是出现虚频(3)采用原胞,截断半径设为4倍,计算的结构显示稳定。1和2为什么会出现虚频呢?

(1)



(2)



(3)
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

乐石||

木虫 (正式写手)



zhuzhudgyg(金币+1): 谢谢参与
关注一下
2楼2012-07-20 10:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzhudgyg

金虫 (正式写手)


引用回帖:
3楼: Originally posted by souledge at 2012-07-20 11:09:54
赞一下计算资源~
超胞会引起相互作用的复杂化,对收敛的要求有时会更高,以前见过单胞不稳定,搞成超胞后由于限制作用,却声子稳定了的……
首先尽量提高K点和截断能做好收敛(几何优化时需要的K点和截断要高于收 ...

谢谢夸奖!
这些都是几何优化后的计算,几何优化的截断能350eV,K点设置为10*10*2,收敛还可以,我用的是MS4.3,这个精度已经最高了,声子的精度也是调到最高。另外想问一句,超胞计算声子,怎样设置超胞比较好?我一直用以上的第三种。
4楼2012-07-21 11:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
zhuzhudgyg(金币+1): 谢谢参与
sunyang1988: 金币+1, 谢谢交流 2012-07-21 13:01:07
zhuzhudgyg: 金币+1 2012-07-22 09:52:39
引用回帖:
4楼: Originally posted by zhuzhudgyg at 2012-07-21 11:01:35
谢谢夸奖!
这些都是几何优化后的计算,几何优化的截断能350eV,K点设置为10*10*2,收敛还可以,我用的是MS4.3,这个精度已经最高了,声子的精度也是调到最高。另外想问一句,超胞计算声子,怎样设置超胞比较好? ...

MS提供的精度基本没什么科学意义,所有参数的收敛要自己测试,更重要的,所有的参数都可以自己手动设置,根本就没有所谓的“精度调到最高”一说。MS这货很大的一个毛病就是提供了几个精度选项,让使用者容易发懒,甚至达不到需要的计算精度。
我不知道你是用的超软还是模守赝势,根据我自己做的收敛测试,在CASTEP中,含C的超软收敛截断能至少是400eV以上,而模守更是要至少800eV以上。其他所有的含原子序数10以内元素的截断能也都只高不低。
并且我说过,声子计算对结构等各种参数非常敏感,所以尽量都在收敛的基础上取得更精确一些才可以。LZ向来不喜欢说计算的体系和结构,但是10*10*2的K点,总让人感觉像是在计算表面而不是晶体。
超胞怎么设置显然没有任何规定,自己多计算多分析结果的异同就是了。最重要的,还是要参数收敛。
5楼2012-07-21 11:12:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzhudgyg

金虫 (正式写手)


引用回帖:
5楼: Originally posted by souledge at 2012-07-21 11:12:16
MS提供的精度基本没什么科学意义,所有参数的收敛要自己测试,更重要的,所有的参数都可以自己手动设置,根本就没有所谓的“精度调到最高”一说。MS这货很大的一个毛病就是提供了几个精度选项,让使用者容易发懒, ...

首先表示感谢!
我用的MS4.3的castep,超软默认最高截断能是350,5.0就能达到450了,模守恒最高660,我以为这是MS本身最高的设置,虽然知道可以改,但还以为即使设置了,4.3也达不到这种计算精度呢。我做的是六角密排结构的晶体,含有过渡金属。另外计算表面和晶体有什么区别?怎样才算是计算晶体呢?说得我有点迷茫了。
6楼2012-07-22 09:51:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youqin5912

金虫 (正式写手)



zhuzhudgyg(金币+1): 谢谢参与
楼主电脑配置咋样?超胞算声子谱要比原包多很久吧
7楼2012-12-11 11:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzhudgyg

金虫 (正式写手)


引用回帖:
7楼: Originally posted by youqin5912 at 2012-12-11 11:03:38
楼主电脑配置咋样?超胞算声子谱要比原包多很久吧

我现在用服务器,但是也就四个或者八个核,超胞肯定比原胞要久一些,但是也要准确一些
10楼2012-12-13 09:48:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
su-b088楼
2012-12-11 14:15   回复  
su-b089楼
2012-12-11 14:15   回复  
相关版块跳转 我要订阅楼主 zhuzhudgyg 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见