应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » dmol3算PtRu结构优化失败
81/1返回列表
查看: 2092  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-10-20 20:23:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

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

另外以Ru这样的电子结构截断半径最好大一些，多极展开最好也是大一些(Hex)。至于是否有磁性我不确定，如果有的话最好指定初始自旋态。Pt的表面Dmol3计算确实稍微难收敛一些。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭