24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2147  |  回复: 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 的主题更新
信息提示
请填处理意见