| 查看: 1566 | 回复: 0 | ||
[求助]
用高斯搜索过渡态反应物和产物自旋多重度不一样怎么办？
|
|
如题，用高斯搜索过渡态，反应物自旋多重度为1，产物自旋多重度为5，过渡态自旋多重度应该怎么设？我设自旋多重度为1搜到过渡态后，跑IRC很容易搜到反应物，但是怎么跑到自旋多重度为5的产物很迷惑，我可以用过渡态的结构把自旋多重度改为5来进行计算吗？这是否意味着我的结构在过渡态瞬间电子结构改变？这样做对吗？ 附上计算脚本 1.先扫描键长 #p opt=modredundant b3lyp/gen geom=connectivity pseudo=read scfcyc=200 symm=veryloose Fe(H2B(pz)2)2phen_2 6-31g 0 1 Fe 0.00000000 0.00000000 2.62521400 S -1.24849800 1.31269500 -6.46666200 S 1.24849800 -1.31269500 -6.46666200 N -1.25380300 0.39111200 1.08503100 分子坐标好长，不写了 2.找到可能的过渡态，计算频率 #p freq ub3lyp/gen pseudo=read scfcyc=200 symm=veryloose Fe_phen_2freq 0 1 Fe 2.54767100 0.00002800 0.00002100 S -6.28746100 -2.17776300 0.42803200 S -6.28746900 2.17775800 -0.42805800 这里还是用自旋多重度=1计算的，结果发现只有1个-200的虚频，很开心 3.人为把自旋多重度改为5进行优化计算，跑出来的结果跟产物很相像，但是有一个-8.71的虚频，震动方向如fig1红色箭头所示，我是否需要把它重新调关键词优化？要怎么调关键词呢？可以直接沿震动方向调坐标优化吗？ #p opt freq=noraman ub3lyp/gen iop(3/24=1,5/13=1) pseudo=read symm=veryloose ts1 6-31+g* optfreq 0 5 Fe 0.00000000 0.00000000 2.53767100 S 2.27317600 -0.66823400 -6.22746100 S -2.27317600 0.66823400 -6.22746100 4.分别跑IRC找产物和反应物，问题来了 a）自旋多重度设为1找反应物 #p irc=(CalcFC,forward,maxpoints=10,stepsize=20) b3lyp/gen pseudo=read ts1_irc forward 0 1 Fe 2.54767100 0.00002800 0.00002100 S -6.28746100 -2.17776300 0.42803200 S -6.28746900 2.17775800 -0.42805800 很顺利结束 b）自旋多重度设为5找产物 #p irc=(CalcFC,reverse,maxpoints=30,stepsize=20,maxcycle=400) IOp(5/13=1) b3lyp/gen pseudo=read ts1_irc reverse 0 5 Fe 2.54767100 0.00002800 0.00002100 S -6.28746100 -2.17776300 0.42803200 S -6.28746900 2.17775800 -0.42805800 各种报错，各种调参数,报错信息如下： Cartesian Forces: Max 0.027937985 RMS 0.005666714 Z-matrix is all fixed cartesians, so copy forces. Leave Link 716 at Wed Dec 9 11:03:35 2015, MaxMem= 4026531840 cpu: 1.0 (Enter /public/software/g09/l123.exe) IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC SCF convergence failure! Path Direction Number 1 Point Number 2 Unrecoverable SCF failure Error termination via Lnk1e in /public/software/g09/l123.exe at Wed Dec 9 11:03:35 2015. Job cpu time: 24 days 21 hours 43 minutes 58.5 seconds. File lengths (MBytes): RWF= 17086 Int= 0 D2E= 0 Chk= 106 Scr= 2 求助各位大神   fig1.JPG |
回复此楼» 猜你喜欢
science advances的AAAS系统上返修问题
已经有4人回复
-
LDPE黑十字为什么那么大
已经有0人回复
-
有机高分子材料论文润色/翻译怎么收费?
已经有296人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有159人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有159人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有157人回复
-
西南交通大学医学院2027年接收推荐免试攻读研究生（含直博生）
已经有152人回复
-
博士申请
已经有15人回复
-
【中科院杭州医学研究所27级博士招生】化学分析/生物材料与器件方向
已经有15人回复
-
静电纺丝膜厚度问题
已经有94人回复
-
常规静电纺丝制备
已经有106人回复
» 本主题相关商家推荐: (我也要在这里推广)
