24小时热门版块排行榜    

查看: 2428  |  回复: 9

行侠仗义

铁虫 (正式写手)

[求助] MP2做一个简单分子,结构优化收敛,频率计算不收敛 已有3人参与

只有两个原子,Li2+

输入文件如下:

%mem=12gb
%chk=123456.chk
%nprocs=12
# mp2/6-311++g(3df,3pd) opt(tight,maxcyc=500) freq scf=maxcyc=1000

****

1 2
Li 0. 0. 0.
Li 3.14 0. 0.


输出文件:

                      1
                     SGG
Frequencies --   -682.3791
Red. masses --      7.0160
Frc consts  --      1.9248
IR Inten    --      0.0000
  Atom  AN      X      Y      Z
     1   3     0.00   0.00   0.71
     2   3     0.00   0.00  -0.71

         Item               Value     Threshold  Converged?
Maximum Force            0.016439     0.000450     NO
RMS     Force            0.016439     0.000300     NO
Maximum Displacement     0.132967     0.001800     NO
RMS     Displacement     0.188043     0.001200     NO

有虚频率,并且频率不收敛。请问各位如何解决?

我用HF, B3LYP算这个分子的时候,是没有虚频率的!
回复此楼

» 猜你喜欢

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

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

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

tcclab

银虫 (正式写手)

基组里面为什么两个加号?
2楼2014-11-30 20:51:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
gmy1990: 金币+2, 感谢应助 2014-12-01 03:31:36
这个计算只有1个活性电子,因此MP2等同于HF计算。这个例子说明Gaussian没有考虑一些特殊情况。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2014-11-30 21:51:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

行侠仗义

铁虫 (正式写手)

引用回帖:
3楼: Originally posted by beefly at 2014-11-30 21:51:03
这个计算只有1个活性电子,因此MP2等同于HF计算。这个例子说明Gaussian没有考虑一些特殊情况。

可是优化的结构,频率与HF结果不同啊。HF能正常结束,没有虚频率。
4楼2014-12-01 14:20:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxing0827

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
这个问题嘛,首先MP2在目前来说都是相对精确的方法,比常用的DFT方法要精确的多,所以在计算时间上也要多的多。一般来说不收敛,采用的方法是增加迭代次数,也就是maxcyc这一项。也可以试试加上opt=noeigen,这个关键词的意思是不必每一步都计算。
自信,理解
5楼2014-12-01 14:48:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

行侠仗义

铁虫 (正式写手)

引用回帖:
5楼: Originally posted by yangxing0827 at 2014-12-01 14:48:15
这个问题嘛,首先MP2在目前来说都是相对精确的方法,比常用的DFT方法要精确的多,所以在计算时间上也要多的多。一般来说不收敛,采用的方法是增加迭代次数,也就是maxcyc这一项。也可以试试加上opt=noeigen,这个关 ...

计算很快啊。只有两个原子而已。优化部分是收敛的。但是频率部分不收敛。我所用的迭代次数500次。
6楼2014-12-01 15:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

引用回帖:
4楼: Originally posted by 行侠仗义 at 2014-12-01 14:20:15
可是优化的结构,频率与HF结果不同啊。HF能正常结束,没有虚频率。...

上面已经说了,这个计算就不应该用mp2。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
7楼2014-12-01 21:02:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MSLABNK

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
为什么不调节一下键长重新优化试一试呢?我觉得有时是开始设的键长不合适,建议稍微调整一下。
m/z234.3456
8楼2014-12-02 09:35:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kekexili_08

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
行侠仗义: 金币+5, ★★★★★最佳答案, 解决问题。 2014-12-03 16:12:47
用mp2(full) 做, 结果会正常。
9楼2014-12-03 13:16:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

行侠仗义

铁虫 (正式写手)

引用回帖:
9楼: Originally posted by kekexili_08 at 2014-12-03 13:16:48
用mp2(full) 做, 结果会正常。

谢谢!
10楼2014-12-03 16:12:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 行侠仗义 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:47 by zNcEhTcH7ihq
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 18:59 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-15 15:27 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见