24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1448  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 297,工科调剂?河南农业大学本科 +15 河南农业大学-能 2026-04-14 15/750 2026-04-20 12:39 by Equinoxhua
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
信息提示
请填处理意见