24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1586  |  回复: 8

布布_

铜虫 (初入文坛)

[求助] 请大家帮助我解决下找过渡态中的燃眉之急 谢谢!卡好久了。。。已有3人参与

最近找过渡态都是TS先用HF/3-21G找到后再用大基组m062x/6-31g*跑,但是最近的几个过渡态都是用小基组找到了,但是读CHK时候 跑不出来,两个途径的问题分别是两种情况:1.一种是过渡态能量一直小范围震荡,都到0.000X了但就是不平,需要加什么关键词或者怎样?加gdiis管用么?或者加IOP=(1/8=8)? 2:一种是一个成环(C=C和一个N成环,成环后下一步是H转移到C=C另一个C上),小基组虚频都是400多500多,但是读CHK后成环两个原子反而跑远了,无法成环,体系中氢键作用比较多,是不是HF和M062X方法之间的误差问题?有什么好办法解决么?基函数一千一左右,不用小基组直接TS的话跑1002就得一天多。。。太费时间了。。请大家指点下我该咋整,一个多月一点进展都没有,很焦急。。先谢谢大家了,请大家不吝赐教,感激不尽!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2014-07-04 18:11:06
布布_: 金币+30, ★★★很有帮助 2014-07-08 08:50:18
1、可以加gdiis,如果再没用的话我一般是把振荡的那几个结构拿出来单独算频率看看。
2、6-31g*我觉得不算大机组。大机组怎么也得到6-311++g(d,p)吧。。。而且你氢键多,为何不加弥散?我觉得加到6-31++g(d)。用HF给个初始结构和初始猜我觉得没问题,但是用来算弱相互作用我个人认为是“无意义”的, 如果是高斯09d以上版本的话可以用DFT-D,解决弱相互作用比较不错。

http://hi.baidu.com/sobereva/item/3cbf8904defdc4803d42e2e9
大神的这篇文章可供参考。

个人见解,欢迎批评指正。
2楼2014-07-04 10:10:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

再加一句,算弱相互作用,还是很多的,就不要考虑小机组了,能上大的就上大 的吧。。。
3楼2014-07-04 10:12:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jyzhao1981

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2014-07-04 18:11:14
布布_: 金币+30, ★★★很有帮助 2014-07-08 08:50:30
对于你的第一个问题,加gdiis能解决。第二个反应直接用m062x/6-31g**做一个柔性扫描去寻找过渡态。你的体系中有氢键,泛函选的对,但基组至少要用6-31g**了。算1002花一天多的时间多正常的事情啊,你怎么那么急呢。对于含有非键作用的体系,HF很难给出正确的结果,它忽略了大部分电子相关啊。
人生需要多少等待?
4楼2014-07-04 12:32:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

布布_

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2014-07-04 10:10:53
1、可以加gdiis,如果再没用的话我一般是把振荡的那几个结构拿出来单独算频率看看。
2、6-31g*我觉得不算大机组。大机组怎么也得到6-311++g(d,p)吧。。。而且你氢键多,为何不加弥散?我觉得加到6-31++g(d)。用HF给 ...

谢谢~
5楼2014-07-08 08:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

布布_

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by jyzhao1981 at 2014-07-04 12:32:12
对于你的第一个问题,加gdiis能解决。第二个反应直接用m062x/6-31g**做一个柔性扫描去寻找过渡态。你的体系中有氢键,泛函选的对,但基组至少要用6-31g**了。算1002花一天多的时间多正常的事情啊,你怎么那么急呢。 ...

好吧 其实我觉得最好的么单独找也很费时间 因为一直不对。。
6楼2014-07-08 08:51:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

布布_

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by 枪下游魂 at 2014-07-04 10:10:53
1、可以加gdiis,如果再没用的话我一般是把振荡的那几个结构拿出来单独算频率看看。
2、6-31g*我觉得不算大机组。大机组怎么也得到6-311++g(d,p)吧。。。而且你氢键多,为何不加弥散?我觉得加到6-31++g(d)。用HF给 ...

那我做反应时候是要直接用m062x-6-31++g(d)比较好么、、那样就是会算的很慢。。 用ORCA可以么?
7楼2014-07-08 10:17:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crosschannel

金虫 (正式写手)

【答案】应助回帖

有环的话,优化时候出现断裂我也遇见过
我试过把坐标换一下,直角换为内坐标,然后重新算,就没事了
8楼2014-07-09 07:42:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

布布_

铜虫 (初入文坛)

引用回帖:
8楼: Originally posted by crosschannel at 2014-07-09 07:42:39
有环的话,优化时候出现断裂我也遇见过
我试过把坐标换一下,直角换为内坐标,然后重新算,就没事了

您好 我的体系是这步过渡态要成环 但是远离了 您说的是本来就是一个环还是要成环?谢谢
9楼2014-07-09 11:03:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 布布_ 的主题更新
信息提示
请填处理意见