24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1439  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 296求调剂 +10 汪!?! 2026-04-09 11/550 2026-04-15 09:17 by fenglj492
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +6 626776879 2026-04-08 6/300 2026-04-14 20:09 by gwjxiaolang
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 化工学硕294分,求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 071000生物学调剂求助 +18 zzzzwww 2026-04-09 21/1050 2026-04-14 15:39 by zs92450
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 复试调剂,一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
信息提示
请填处理意见