24小时热门版块排行榜    

查看: 2719  |  回复: 7

fxh0402

铁虫 (小有名气)

[求助] dmol3算PtRu结构优化失败

我在MS里面导入了Pt晶胞,然后把顶角的Pt原子替换成了Ru原子,然后进行切面111面,2*2的晶胞,4层,固定住了下面两层,参数:use symmetry没有选,spin unrestricted 勾不勾选都是失败,迭代次数设置为500,也挺大了吧?smearing设置为0.009或者0.01也够大了吧?不勾选optimize cell只优化原子,不能再简单了吧?还有没有别的参数修改办法了?怎么办??急!!
Message: SCF not converging. Choose larger smearing value in DMol3 SCF panel or modify/delete ''Occupation Thermal" in the input file. You may also need to change spin or use symmetry.
Resubmit DMol3.
Message : DMol3 job failed
Error : DMol3 exiting
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

smearing可加速收敛,但是也有一定的代价,并且不解决所有问题。
既然Pt和Ru都是50%占据,建议对晶格进行优化,Pt原来的结构是面心立方的结构,替换原子后如果还是这样的结构,我建议拟合三阶B-M方程确定晶格常数,Dmol3本身的优化晶格方法是不准确的。

另外以Ru这样的电子结构截断半径最好大一些,多极展开最好也是大一些(Hex)。至于是否有磁性我不确定,如果有的话最好指定初始自旋态。Pt的表面Dmol3计算确实稍微难收敛一些。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-10-21 11:25:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxh0402

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by 卡开发发 at 2014-10-21 11:25:36
smearing可加速收敛,但是也有一定的代价,并且不解决所有问题。
既然Pt和Ru都是50%占据,建议对晶格进行优化,Pt原来的结构是面心立方的结构,替换原子后如果还是这样的结构,我建议拟合三阶B-M方程确定晶格常数, ...

把smearing设大,有什么代价呢?是结构不可靠吗?但是我发现我的周期平板结构一直没有变化啊。
Pt原来的结构是面心立方的,我把顶角的原子替换为Ru原子,Pt:Ru应该是3:1的
3楼2014-10-22 08:25:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by fxh0402 at 2014-10-22 08:25:11
把smearing设大,有什么代价呢?是结构不可靠吗?但是我发现我的周期平板结构一直没有变化啊。
Pt原来的结构是面心立方的,我把顶角的原子替换为Ru原子,Pt:Ru应该是3:1的...

能量计算不可靠,不信可以看看总能以及外推到0K的能量,这个问题论坛应该有不少人论述过。不管比例吧(上年纪了智商捉鸡了一下),最好考虑晶格结构的弛豫,否则后面优化过程搞不好力难收敛。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-10-22 10:11:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxh0402

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by 卡开发发 at 2014-10-22 10:11:09
能量计算不可靠,不信可以看看总能以及外推到0K的能量,这个问题论坛应该有不少人论述过。不管比例吧(上年纪了智商捉鸡了一下),最好考虑晶格结构的弛豫,否则后面优化过程搞不好力难收敛。...

我是初学者,不是很懂您的意思,您可以说的详细点儿吗?总能是在优化好的.outmol ?OK的能量是单点能energy吗?
还有您说的晶格结构的弛豫怎么考虑呢?
5楼2014-10-23 09:15:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
5楼: Originally posted by fxh0402 at 2014-10-23 09:15:41
我是初学者,不是很懂您的意思,您可以说的详细点儿吗?总能是在优化好的.outmol ?OK的能量是单点能energy吗?
还有您说的晶格结构的弛豫怎么考虑呢?...

如果是采用了合适的smearing,我觉得采用的能量应该取
Df  binding energy extrapolated to T=0K     xxxx Ha      xxxx eV
不一定挂在论坛,计算问题问题欢迎留言。
6楼2014-10-23 13:04:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
5楼: Originally posted by fxh0402 at 2014-10-23 09:15:41
我是初学者,不是很懂您的意思,您可以说的详细点儿吗?总能是在优化好的.outmol ?OK的能量是单点能energy吗?
还有您说的晶格结构的弛豫怎么考虑呢?...

一般而言smearing的值不宜过大,输出文件中有
Fermi Energy:          -0.190638 Ha      -5.188 eV   -kTS_e= -0.0049465Ha
其中 -kTS_e项最好不要超过1meV/atoms
晶格弛豫如果满足立方晶格,直接做不同晶格长度与能量的曲线来拟合三阶BM方程,具体做法论坛有人讨论。
不一定挂在论坛,计算问题问题欢迎留言。
7楼2014-10-23 13:06:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浴雪飘零_

新虫 (小有名气)

这种情况不优化晶包不好吧

发自小木虫Android客户端
8楼2016-12-10 01:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fxh0402 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 7/350 2026-11-02 12:46 by n0P4fnsXKbBK
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 12:25 by n0P4fnsXKbBK
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 06:29 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:13 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 05:47 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:37 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
信息提示
请填处理意见