24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1352  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhuzhudgyg

金虫 (正式写手)


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

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

(1)



(2)



(3)
回复此楼

» 猜你喜欢

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

查看全部散金贴

已阅   回复此楼   关注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的回帖
查看全部 10 个回答
★ ★ ★ ★ ★
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(金币+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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑大材料工程290求调剂 +12 Youth_ 2026-03-30 12/600 2026-03-31 03:34 by 蒙奇奇521
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +9 @taotao 2026-03-29 9/450 2026-03-30 22:29 by 我是小康
[考研] 抱歉 +4 田洪有 2026-03-30 4/200 2026-03-30 21:26 by mumin1990
[考研] 085602化工求调剂(331分) +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 一志愿华东师范大学有机化学专业,初试351分,复试被刷求调剂! +7 真名有冰 2026-03-29 8/400 2026-03-30 21:08 by maaj02
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[硕博家园] 求调剂 有机化学考研356分 +10 Nadiums 2026-03-25 11/550 2026-03-30 11:56 by yjolah
[考研] 296求调剂 +10 彼岸t 2026-03-29 10/500 2026-03-30 10:50 by 探123
[考研] 375求调剂 +6 雨夏整夜 2026-03-29 6/300 2026-03-30 10:21 by herarysara
[考研] 求调剂 +4 QiMing7 2026-03-25 5/250 2026-03-29 21:10 by 唐沐儿
[考研] 数一英一271专硕(085401)求调剂,可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
[考研] 266分,求材料冶金能源化工等调剂 +7 哇呼哼呼哼 2026-03-27 9/450 2026-03-28 12:22 by zllcz
[考研] 283求调剂 +7 A child 2026-03-28 7/350 2026-03-28 12:05 by zllcz
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 0703化学338求调剂! +6 Zuhui0306 2026-03-26 7/350 2026-03-27 10:35 by shangxh
[考研] 324求调剂 +5 hanamiko 2026-03-26 5/250 2026-03-27 10:33 by wangjy2002
[硕博家园] 北京林业大学硕导招生广告 +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
信息提示
请填处理意见