24小时热门版块排行榜    

查看: 3169  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

2328980933

新虫 (著名写手)

[求助] dmol优化铁簇的问题

我最近在研究铁簇与高分子吸附,但是铁簇优化比较难,参数设置经常会导致铁簇优化失败,通过与文献铁簇键长2.35埃米对比,文献是用高斯算的。我调整了参数,找到了一个最接近文献的键长2.37埃米,但是其他的参数都在2.41左右。我按照最好的那个参数优化高分子是成功的,但是把他们放在一块就不能成功,而我又不能改变参数,因为目前只找到了一个参数是铁簇的键长与文献比较接近。求助大家畅所欲言,看看有没有什么办法,有没有遇到此问题的或者算过这个的,谢谢大家
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★ ★ ★
2328980933(ben_ladeng代发): 金币+5, 代楼主发放金币 2014-06-03 06:35:38
引用回帖:
8楼: Originally posted by 2328980933 at 2014-05-26 15:44:02
参数为# Task parameters
Calculate                     optimize
Opt_energy_convergence        2.0000e-005
Opt_gradient_convergence      4.0000e-003 A
Opt_displacement_convergence  5.0000e-003 A
O ...

信息仍然不够多,实在不行把结构的截图发来看看。

这里暂时只能知道SCF没有收敛,对此给出几个建议:
1、调整铁原子的自旋状态
2、增加计算精度
3、适当减小charge mix
4、加大DIIS size

另外,题外话,个人观点是包含纯金属的体系不建议用DFT-D,一般得到的金属键长偏短。
不一定挂在论坛,计算问题问题欢迎留言。
9楼2014-05-26 17:55:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

Fe要注意调整自旋,这个一般是造成scf不收敛比较重要的原因。另外Fe这种带有磁性的体系Gaussian处理可能不会很理想,对于文献的结果考证后在去模仿也不迟。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-05-26 13:25:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2328980933

新虫 (著名写手)

引用回帖:
2楼: Originally posted by 卡开发发 at 2014-05-26 13:25:31
Fe要注意调整自旋,这个一般是造成scf不收敛比较重要的原因。另外Fe这种带有磁性的体系Gaussian处理可能不会很理想,对于文献的结果考证后在去模仿也不迟。

自旋我也都试过,自动的是最小的,就是和高分子在一块的时候不能优化,很奇怪,而且根本不出收敛曲线
3楼2014-05-26 14:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by 2328980933 at 2014-05-26 14:10:35
自旋我也都试过,自动的是最小的,就是和高分子在一块的时候不能优化,很奇怪,而且根本不出收敛曲线...

那是因为scf 比较困难,Fe自旋推荐先设置到高自旋(+4),不设定的话我试过,有时候会计算成反铁磁性。
不如把outmol文件上传过来,以便于诊断。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-05-26 14:12:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:23 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 17:01 by Xc2FxHsHkBNc
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:06 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
请填处理意见