24小时热门版块排行榜    

查看: 2154  |  回复: 8

alanzxw

金虫 (初入文坛)

[求助] 【感激涕零求助】关于含有Mn原子体系的自旋多重度确认,而且体系计算不成功 已有2人参与

,本人菜鸟,第一次接触过渡金属的高斯优化。
如题,我的体系是一个有机分子+一个Mn原子的优化,体系有机分子和Mn电荷都是零,命令和坐标如下:

#P B3LYP/genecp scrf=(solvent=acetone) Opt freq=noraman
0 2
C                  0.00000000    0.00000000    0.00000000
C                  0.00000000    0.00000000    1.52892272
H                  1.00274739    0.00000000   -0.42788123
H                 -0.59787084    0.79368669   -0.44037257
H                  0.92369613    0.36628535    1.96930727
H                 -0.85609160    0.52211507    1.95680667
C                 -0.58185992   -2.07147763    0.76444215
O                 -0.90640325   -3.22689151    0.76442698
O                 -0.12906425   -1.40837558    1.84587390
O                 -0.62315255   -1.26958314   -0.31697957
Mn                -2.78061675   -2.20345148    0.91949219

C H O 0
6-311++G(d)
****
Mn 0
lanl2dz
****

Mn 0
lanl2dz

在论坛看到对于确定Mn和Fe的自旋多重度是比较难的,所以优化错误我一度怀疑是多重度设置错误,另外错误信息的最后几行如下:
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -446.282939914     A.U. after  129 cycles
             Convg  =    0.3981D-03             -V/T =  2.1194
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 1.8989 S= 0.9659
<L.S>= 0.000000000000E+00
KE= 3.986750007943D+02 PE=-1.786236867846D+03 EE= 5.696120174817D+02
Annihilation of the first spin contaminant:
S**2 before annihilation     1.8989,   after     1.2869
Convergence failure -- run terminated.
Error termination via Lnk1e in /opt/gauss/g09/l502.exe at Tue Apr 21 00:24:06 2015.
Job cpu time:  0 days  1 hours 11 minutes 18.1 seconds.
File lengths (MBytes):  RWF=     62 Int=      0 D2E=      0 Chk=      4 Scr=      1
求大神!!!在线等!
回复此楼

» 猜你喜欢

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

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

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

alanzxw

金虫 (初入文坛)

补充:这个计算一个点都没有跑,输出文件就一个点,GaussView打开已经Read Intermediate Geometries(Gaussian Optimizations Only)了。
2楼2015-04-20 22:12:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kong890109

金虫 (小有名气)


gmy1990: 金币+1, 谢谢交流 2015-04-21 17:52:09
你这是scf不收敛,可以增大scf的循环次数,scf=maxcycle=512,但是感觉还是你的初始结构有问题,可能还不管用。
3楼2015-04-21 09:45:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
gmy1990: 金币+1, 谢谢交流 2015-04-21 17:52:17
如果只跑了一步说明初始结构有问题,建议调整初始结构。
和多重度关系不大,多重度就是孤电子数+1,还是比较好确认的
4楼2015-04-21 12:54:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alanzxw

金虫 (初入文坛)

引用回帖:
4楼: Originally posted by 枪下游魂 at 2015-04-21 12:54:38
如果只跑了一步说明初始结构有问题,建议调整初始结构。
和多重度关系不大,多重度就是孤电子数+1,还是比较好确认的

我设置多重度为6算出一个优化结构了,我在纳闷,多重度=孤电子数+1,以Mn原子为例,这里孤电子是一下哪种算法呢:
1、25/2=12......1,取孤电子数=1;
2、按Mn的价电子排布,3d(5)4s(2),半满规则,这里孤电子数取5.
5楼2015-04-21 16:02:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kong890109

金虫 (小有名气)

自旋2,4,6都可以的,看来你要比较不同自旋的能量啦,看看哪种自旋情况稳定。
6楼2015-04-21 16:17:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alanzxw

金虫 (初入文坛)

引用回帖:
6楼: Originally posted by kong890109 at 2015-04-21 16:17:09
自旋2,4,6都可以的,看来你要比较不同自旋的能量啦,看看哪种自旋情况稳定。

自旋2,4,6我都试了,最后只有6算出来,而且这只是一个双分子体系,结果里面Mn的位置和我初始摆放的位置差别很大,所以我感觉应该不存在初始结构不合理的问题了。我也想比较自旋2,4,6的能量取最低,可是只有一个算得出来,唉。
7楼2015-04-21 17:05:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kong890109

金虫 (小有名气)

以前我也做过一些,我的体系是自旋小的能量高,并且非常不好收敛,可能跟基组和方法也有一定关系。
8楼2015-04-22 09:05:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

引用回帖:
5楼: Originally posted by alanzxw at 2015-04-21 16:02:31
我设置多重度为6算出一个优化结构了,我在纳闷,多重度=孤电子数+1,以Mn原子为例,这里孤电子是一下哪种算法呢:
1、25/2=12......1,取孤电子数=1;
2、按Mn的价电子排布,3d(5)4s(2),半满规则,这里孤电子数 ...

可以阅读下ligand field theory相关资料,可以帮助你理解过度金属的电子在d轨道占居
9楼2015-04-24 19:03:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 alanzxw 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 10/500 2026-03-18 20:41 by Wangjingyue
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +4 1孙悟空 2026-03-17 4/200 2026-03-18 17:59 by fivewind
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
请填处理意见