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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【感激涕零求助】关于含有Mn原子体系的自旋多重度确认，而且体系计算不成功
91/1返回列表
查看: 1965  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
求大神！！！在线等！
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2015-04-20 21:52:34
已阅   回复此楼   关注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轨道占居
赞一下(1人) 回复此楼
9楼2015-04-24 19:03:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 alanzxw 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭