| 查看: 1529 | 回复: 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 |
回复此楼» 猜你喜欢
2026年中南民族大学化学与材料科学学院调剂公告
已经有40人回复
-
武汉纺织大学招生物材料、化学、高分子相关专业硕士生
已经有0人回复
-
有机高分子材料论文润色/翻译怎么收费?
已经有257人回复
-
武汉纺织大学招生物材料、化学、高分子相关专业硕士生-另有少量士兵计划可咨询
已经有0人回复
-
[武汉纺织大学]材料学院[姜伟团队名称]课题组-2026年硕士研究生调剂信息
已经有0人回复
-
中国科学院·过程工程研究所·招聘海内外高层次人才---过程所“百人计划”
已经有1人回复
-
中国科学院·过程工程研究所·招聘海内外高层次人才---过程所“百人计划”
已经有5人回复
-
2026考博自荐
已经有4人回复
-
温敏水凝胶收缩
已经有0人回复
-
PDMS膜在溶液中的离子透过率
已经有0人回复
-
攒人品 论文祈福
已经有139人回复
» 本主题相关商家推荐: (我也要在这里推广)
40