24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1402  |  回复: 14
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 weiliu6762 的 2 个金币

weiliu6762

铜虫 (小有名气)

[交流] 【讨论】DMol3计算Fe(110)表面不收敛

不知道大家有没有遇到这种情况。我用dmol3计算bcc Fe(110)面,不知道为什么,怎么算都是不收敛。我用的是GGA-PW91,选上了自旋,DSPP,DNP,smearing设成0.004。计算后迭代了1000步仍然不收敛。请大家帮我分析分析。多谢!!
回复此楼

» 猜你喜欢

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

徐亮8090

金虫 (正式写手)

2楼2009-06-09 13:00:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

徐亮8090

金虫 (正式写手)


zxzj05(金币+1,VIP+0):^_^ 3Q!鼓励答疑! 6-9 22:01
你算的是Fe团簇吗
3楼2009-06-09 13:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiliu6762

铜虫 (小有名气)

不是啊,就是平的金属表面。
4楼2009-06-09 13:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archhector

★ ★
zxzj05(金币+2,VIP+0):^_^ 3Q!鼓励答疑! 6-9 22:01
引用回帖:
Originally posted by weiliu6762 at 2009-6-9 11:29:
不知道大家有没有遇到这种情况。我用dmol3计算bcc Fe(110)面,不知道为什么,怎么算都是不收敛。我用的是GGA-PW91,选上了自旋,DSPP,DNP,smearing设成0.004。计算后迭代了1000步仍然不收敛。请大家帮我分析分 ...

你选了自旋,设置了Fe原子的初始自旋值么?你的体系初始自旋值是不是显示的0啊。一般磁性体系不收敛的话必须调整磁性原子的spin值
5楼2009-06-09 17:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)


wuchenwf(金币+1,VIP+0):谢谢 6-10 21:33
smear调大些
储氢家族欢迎储氢研究者!
6楼2009-06-09 22:02:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiliu6762

铜虫 (小有名气)

我选的就是“use formal spin as initial”,看来这样是不行的。请问你知道对于Fe到底该怎么选spin数值吗?我的是(2*2)的Fe(110)表面,4层,一共16个原子。谢谢!
引用回帖:
Originally posted by archhector at 2009-6-9 17:41:

你选了自旋,设置了Fe原子的初始自旋值么?你的体系初始自旋值是不是显示的0啊。一般磁性体系不收敛的话必须调整磁性原子的spin值

7楼2009-06-09 22:14:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

archhector

★ ★
wuchenwf(金币+2,VIP+0):谢谢 6-10 21:33
引用回帖:
Originally posted by weiliu6762 at 2009-6-9 22:14:
我选的就是“use formal spin as initial”,看来这样是不行的。请问你知道对于Fe到底该怎么选spin数值吗?我的是(2*2)的Fe(110)表面,4层,一共16个原子。谢谢!

做表面的话你取的层数应该少了点吧。选中铁原子,在modify里面的spin修改。
一般来说原子的high spin和low spin它是给了你的,在最下面一行改成它们两个值的差就行了。Fe的话一般是4,每个原子都要改。然后再用use formal spin as initial。我也有磁性体系无法收敛的时候,需要修改spin,尝试把个别原子spin值改小,后面就很繁琐了。总之磁性体系的调试很麻烦。
8楼2009-06-10 15:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水月慧心

金虫 (正式写手)

引用回帖:
Originally posted by archhector at 2009-6-10 15:59:

做表面的话你取的层数应该少了点吧。选中铁原子,在modify里面的spin修改。
一般来说原子的high spin和low spin它是给了你的,在最下面一行改成它们两个值的差就行了。Fe的话一般是4,每个原子都要改。然后再用 ...

同感啊 磁性体系真的好难收敛,我做一个氧化铁的表面结构,优化了不下几十遍了 还是没法收敛,我用的是默认的4,看来改小点试试看,我现在都没辙了 只能改smearing的值,而且要改的 很大才可能收敛........你说的后面的调式很繁琐 能不能具体的讨论一下呢?
9楼2009-12-15 10:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goodtimega

铁杆木虫 (著名写手)


zxzj05(金币+1,VIP+0):谢谢回帖交流 12-15 20:36
精度调低点
10楼2009-12-15 10:29:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 weiliu6762 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学 +29 妮妮ninicgb 2026-04-04 32/1600 2026-04-08 22:50 by 猪会飞
[考研] 311求调剂 +5 surte 2026-04-08 11/550 2026-04-08 22:10 by 土木硕士招生
[考研] 307分材料专业求调剂 +12 Hll胡 2026-04-05 12/600 2026-04-08 16:33 by luoyongfeng
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +3 momo(上岸版) 2026-04-06 6/300 2026-04-08 14:53 by 化学521
[考研] 085404,285分求调剂 +4 薇薇考研 2026-04-07 4/200 2026-04-08 14:32 by zhq0425
[考研] 312求调剂 +4 Say Never 2026-04-04 4/200 2026-04-08 08:41 by barlinike
[考研] 0854电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-05 6/300 2026-04-07 22:16 by hemengdong
[考研] 315求调剂 +17 小羊小羊_ 2026-04-02 18/900 2026-04-07 22:01 by lijunpoly
[考研] 081700,311,求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 一志愿武理车辆专硕总分 281 求调剂 +4 上岸研究生. 2026-04-02 4/200 2026-04-07 09:52 by 加油向未来啊
[考研] 华南理工0703化学,总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 一志愿武汉理工大学-085601材料工程(专硕)-总分353求调剂 +3 2626262626li 2026-04-02 3/150 2026-04-06 09:08 by 无际的草原
[考研] 求调剂 +6 朔朔话 2026-04-02 7/350 2026-04-04 19:16 by 蓝云思雨
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 土木水利328分求调剂 +6 疾风知劲草666 2026-04-02 6/300 2026-04-03 11:38 by znian
[考研] 338求调剂,一志愿能源动力,外语是日语203 +5 zzz,,r 2026-04-02 5/250 2026-04-03 09:45 by 蓝云思雨
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
信息提示
请填处理意见