24小时热门版块排行榜    

查看: 2297  |  回复: 24

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教! 2011-11-02 09:54:36
引用回帖:
9楼: Originally posted by huazhorg at 2011-10-27 10:52:50:
ecutwfc=45已经不小了,这里的单位是Ry, 45Ry就相当于600多eV 的截断,已经相当的高了

这个其实并不一定的,要因元素以及赝势的类型不同来说。比如说,对于rrkj类型的USPP,要到42 Ry以上才能收敛,如果要做高压的,就需要更高了。而且lz的计算里面包含了d电子,这也是需要很大的cutoff的。所以不能一概而论。
11楼2011-10-27 13:31:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★
贺仪(金币+3): 多谢指教! 2011-11-02 09:54:49
引用回帖:
10楼: Originally posted by huazhorg at 2011-10-27 10:59:09:
以他这么大的单胞,我认为8X8X8的K点就够了,振动出现虛频,有可能是结构没有优化好,也有可能是势能过软,换用较硬的赝势试试。另外,In原子较重,可以试试加U,再看看结果如何。

这个单胞一点都不大,声子的计算需要非常密的K点计算波函,才能给出比较准确的结果。对于一个6方的单胞来说,8x8x8的一般来说肯定不够。
至于优化,这个真不能这么说。没有经过收敛测试的优化,意义并不大,并且一般不支持建立个模型就直接优化。不要被优化带个“优”字就给唬住觉得有它就成了。做第一性计算是这样,用操作系统也是这样。一般我们都管操作系统中的各种优化软件叫做“忧”化软件的。
至于加U,个人至少不确定声子计算支持不支持+U后的自洽结果,因为这个需要在代码级别上修改对E的二次求导~
思想重于技巧,内涵重于表象
12楼2011-10-27 14:15:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by goldenfisher at 2011-10-26 22:36:37:
测试很简单,就是把截断从小到达取不同数值进行自洽场迭代计算,看不同截断的时候的能量收敛精度。如果基本前后两步的数值基本一致,则可认为收敛。能量或者压强的收敛与截断、K点选择等有关系,你应该不同的参数 ...

您好,我现在想做收敛测试,我想问一下,在对截断收敛测试时,那网格怎么取呢?难道随便取网格比如4*4*4的,然后把截断从小到大取不同的数值吗?还有,在取不同的数值时,步长一般取多大呢?
什么时候能算完!!!!
13楼2011-10-27 15:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

送鲜花一朵
引用回帖:
6楼: Originally posted by goldenfisher at 2011-10-26 22:36:37:
测试很简单,就是把截断从小到达取不同数值进行自洽场迭代计算,看不同截断的时候的能量收敛精度。如果基本前后两步的数值基本一致,则可认为收敛。能量或者压强的收敛与截断、K点选择等有关系,你应该不同的参数 ...

您好,我已经做了测试,贗势pbe超软,测试后截断70Ry,kmesh 8*8*8,现在已经算了已经算了五天了,感觉应该明天出结果,但是现在出来的结果第一个动力学矩阵还有-84cm-1频支,总感觉做时空间矩阵后还是负值。我晕!您能再给我点建议吗?现在我老板抓住这个问题不放,而且必须用现在含4d电子的贗势还做,我晕,拜托了!!!
什么时候能算完!!!!
14楼2011-11-01 10:59:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
7楼: Originally posted by souledge at 2011-10-26 23:26:21:
你这么说也不管用啊,毕竟测试也是个比较耗时间精力的过程,并且还可能证明了之前看起来跟实验非常“符合”的结果有问题了~
事实上去年我在实验室里做组会报告时,做的就是如何测试收敛的问题,被组里几个也做计 ...

笑了……
低截断能低K-mesh发文章不会被人批死啊?
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
15楼2011-11-01 12:54:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

引用回帖:
15楼: Originally posted by iamikaruk at 2011-11-01 12:54:14:
笑了……
低截断能低K-mesh发文章不会被人批死啊?

所以那哥们的文章,尽管老板觉得非常有创意,但是投了一年多了,愣是没中……我觉得中了才不正常,因为老板从来不肯投1.5以下的期刊~
思想重于技巧,内涵重于表象
16楼2011-11-01 14:34:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
16楼: Originally posted by souledge at 2011-11-01 14:34:27:
所以那哥们的文章,尽管老板觉得非常有创意,但是投了一年多了,愣是没中……我觉得中了才不正常,因为老板从来不肯投1.5以下的期刊~

基本问题,这个很要命啊。
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
17楼2011-11-01 16:39:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

引用回帖:
10楼: Originally posted by huazhorg at 2011-10-27 10:59:09:
以他这么大的单胞,我认为8X8X8的K点就够了,振动出现虛频,有可能是结构没有优化好,也有可能是势能过软,换用较硬的赝势试试。另外,In原子较重,可以试试加U,再看看结果如何。

请问你在这里所说的贗势过软和较硬,在物理上怎么解释呢?是和原子的价电子对贗势的贡献有关系吗?
什么时候能算完!!!!
18楼2011-11-01 20:50:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢指教! 2011-11-02 21:04:04
引用回帖:
18楼: Originally posted by hanyanli0475 at 2011-11-01 20:50:36:
请问你在这里所说的贗势过软和较硬,在物理上怎么解释呢?是和原子的价电子对贗势的贡献有关系吗?

大致来说,较硬的赝势,如NC赝势,赝波函数在定义的核心区域的截止半径以上是符合全电子波函数的。它要求改造后的波函数其在截止半径Rc之内的总电荷量仍要等于未改造前Rc之内总量的大小,这样能保证赝势具有较高的精度,但这也牺牲了计算效率,计算量增加。较软的赝势,如超软赝势,为了减少计算量而不损失太多的精度,波函数被处理得更平滑,所需的平面波基底函数更少,因而计算量也相对较少。

» 本帖已获得的红花(最新10朵)

19楼2011-11-02 09:21:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)

送鲜花一朵
引用回帖:
19楼: Originally posted by huazhorg at 2011-11-02 09:21:00:
大致来说,较硬的赝势,如NC赝势,赝波函数在定义的核心区域的截止半径以上是符合全电子波函数的。它要求改造后的波函数其在截止半径Rc之内的总电荷量仍要等于未改造前Rc之内总量的大小,这样能保证赝势具有较高 ...

我在计算时发现,用NC做计算速度很快,算声子也很快,8*8*8的20个小时就能算完,但是用超软算声子8*8*8的,用了5天。按照您的解释,超软应该更快一些呀,我不太理解。对了本人用的8核处理器但是没有装并行。呵呵!请您再进一步解释一下吧,我对这个问题真的迷惑,多谢指教!
什么时候能算完!!!!
20楼2011-11-02 09:44:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hanyanli0475 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料专业298分有科研经历 硕士研究生调剂自荐信 +5 zyf上岸 2026-03-01 5/250 2026-03-01 11:50 by ms629
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +3 虫虫虫虫虫7 2026-03-01 6/300 2026-03-01 11:50 by gaoxiaoniuma
[考研] 290求调剂 +7 材料专硕调剂; 2026-02-28 8/400 2026-03-01 11:45 by 王伟要上岸啊
[考研] 材料类求调剂 +8 wana_kiko 2026-02-28 8/400 2026-03-01 11:44 by 王伟要上岸啊
[考研] 材料化工调剂 +6 今夏不夏 2026-03-01 6/300 2026-03-01 11:28 by gaoxiaoniuma
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考研] 高分子化学与物理调剂 +5 好好好1233 2026-02-28 9/450 2026-03-01 10:59 by fengyu211
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
信息提示
请填处理意见