24小时热门版块排行榜    

查看: 1593  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 18/900 2026-08-01 13:56 by gltch
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 8/400 2026-08-01 13:53 by y4XuShGyllvT
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 13:49 by y4XuShGyllvT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:49 by y4XuShGyllvT
[基金申请] 小木虫看见有人已经知道结果了 +16 1234567wang 2026-07-25 20/1000 2026-08-01 08:57 by type2type
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见