24小时热门版块排行榜    

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

zhuzhudgyg

金虫 (正式写手)


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

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

(1)



(2)



(3)
回复此楼

» 猜你喜欢

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

查看全部散金贴

已阅   回复此楼   关注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的回帖
查看全部 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

金虫 (正式写手)


引用回帖:
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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学求调剂 +9 小黄鸭宝 2026-03-30 9/450 2026-03-31 13:49 by anny19840123
[考研] 354求调剂 +3 lxb598 2026-03-31 4/200 2026-03-31 13:42 by sophie2180
[考研] 070300化学354求调剂 +13 101次希望 2026-03-28 13/650 2026-03-31 13:38 by anny19840123
[考研] 322求调剂 +10 旧吢 2026-03-24 10/500 2026-03-31 13:31 by 无际的草原
[考研] 311求调剂一志愿合肥工业大学 +9 秋二十二 2026-03-30 9/450 2026-03-31 13:15 by Dyhoer
[考研] 085602 307分 求调剂 +8 不知道叫什么! 2026-03-26 8/400 2026-03-31 12:19 by 1018329917
[考研] 266分,求材料相关专业调剂 +10 哇呼哼呼哼 2026-03-30 12/600 2026-03-31 11:00 by 熊一刀
[考研] 276求调剂 +3 赵久华 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 291求调剂 +12 Y-cap 2026-03-29 15/750 2026-03-31 09:25 by Huaxue_Wang
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 286求调剂 +5 Faune 2026-03-30 5/250 2026-03-30 17:15 by 无际的草原
[考研] 332求调剂 +6 @MZB382400 2026-03-28 6/300 2026-03-30 16:57 by 无际的草原
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 330分求调剂 +5 qzenlc 2026-03-29 5/250 2026-03-29 07:37 by 无际的草原
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
[考研] 材料277求调剂 +7 min3 2026-03-24 7/350 2026-03-28 11:39 by xuxiang
[考研] 0856,材料与化工321分求调剂 +12 大馋小子 2026-03-27 13/650 2026-03-28 10:56 by self2008
[考研] 化学调剂 +4 爱吃番茄的旭 2026-03-24 5/250 2026-03-27 17:50 by kiokin
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-03-26 4/200 2026-03-27 08:04 by chemisry
[考研] 打过很多竞赛,085406控制工程300分,求调剂 +3 askeladz 2026-03-26 3/150 2026-03-26 09:08 by 给你你注意休息
信息提示
请填处理意见