| 查看: 1654 | 回复: 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就得一天多。。。太费时间了。。请大家指点下我该咋整，一个多月一点进展都没有，很焦急。。 先谢谢大家了，请大家不吝赐教，感激不尽！ |
回复此楼» 猜你喜欢
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！
已经有1人回复
-
光电子赋能·半导体创芯！
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有67人回复
-
光电子赋能·半导体创芯！
已经有0人回复
-
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes
已经有2人回复
-
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
cif求助下载
已经有1人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有2人回复
-
高效微纳米储能材料与产业化推进实验室 苏州 招收硕士研究生（含调剂）
已经有0人回复
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2014-07-04 18:11:06
布布_: 金币+30, ★★★很有帮助 2014-07-08 08:50:18
感谢参与，应助指数 +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
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
3楼2014-07-04 10:12:02
jyzhao1981
木虫 (正式写手)
- 应助: 19 (小学生)
- 金币: 3345.5
- 红花: 3
- 帖子: 558
- 在线: 90.7小时
- 虫号: 152951
- 注册: 2006-01-02
- 专业: 物理有机化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2014-07-04 18:11:14
布布_: 金币+30, ★★★很有帮助 2014-07-08 08:50:30
感谢参与，应助指数 +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
5楼2014-07-08 08:49:55
6楼2014-07-08 08:51:00
7楼2014-07-08 10:17:35
crosschannel
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 998.5
- 散金: 28
- 红花: 13
- 帖子: 998
- 在线: 236.3小时
- 虫号: 2615554
- 注册: 2013-08-25
- 专业: 质谱分析
8楼2014-07-09 07:42:39
9楼2014-07-09 11:03:12