查看: 1400  |  回复: 28
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 流浪猫小琪 的 296 个金币 ,回帖就立即获得 2 个金币,每人有 1 次机会
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

流浪猫小琪

铁虫 (正式写手)


[交流] 如何使用CBS-QB3寻找过渡态

已在其他水平下找到了过渡态,打算在此基础上使用CBS-QB3水平下通过IRCmax寻找过渡态,不知道如何写输入文件。想请教一个CBS-QB3做IRCmax的输入文件。
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (正式写手)


这个是因为你用zmatrix来写坐标 确没有给出要优化的变量

发自小木虫IOS客户端
28楼2019-05-24 01:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 29 个回答

paramecium86

版主 (正式写手)



流浪猫小琪(金币+2): 谢谢参与
这个没法一步完成。  
一般来说IRCMax写法如下:
IRCMax(方法2:方法1) 含义是是在方法1的IRC途径上 用方法2在每一个点都做单点能计算,找出其中能量最高的点。  但是CBS-QB3是一个多步骤的方法 其中还包含了第一步一个opt步骤 所以这个不能直接写的。
高斯说明书的原话也是这样说的  : Analytic gradients are required for the IRC portion of the calculation (model1 above). Any non-compound energy method and basis set may be used for model2.  CBS-QB3就是属于 compound energy method。 所以没法直接写。如果需要就只能分步骤人工计算
第一步 先用 方法1 算IRC ,
第二步 自己提取IRC上每一个点的结构 逐一用CBS-QB3做一次单点能 (如果是Gaussian16 就支持CBS(NoOpt)能自动跳过其中的opt步骤 ,给出能量)。但是G09不支持这个写法,那么就要人工来做整个CBS步骤 人为跳过opt步骤。

CBS-QB3除去优化的话需要进行四步骤
1. #p B3LYP/CBSB7 Freq
2. #p  CCSD(T)/6-31+G(d')
3. #p Geom=AllCheck Guess=Read MP4SDQ/CBSB4
4. #P Geom=AllCheck Guess=Read MP2/CBSB3 CBSExtrap=(NMin=10,MinPop)

你可以把这四部分的计算都用--link1--来隔开 写在同一个文件里。

http://www.chem.ucla.edu/~mccarren/CBSQB3_worksheet/  这个网站可以帮你按照CBS-qb3的公式分别取这几个计算中有用的部分加加减减 得到最后的能量。

然后对每一个点都要做一次这些运算。整个运算步骤还是挺多的。

» 本帖已获得的红花(最新10朵)

7楼2019-05-15 14:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪猫小琪

铁虫 (正式写手)


送红花一朵
引用回帖:
7楼: Originally posted by paramecium86 at 2019-05-15 14:27:13
这个没法一步完成。  
一般来说IRCMax写法如下:
IRCMax(方法2:方法1) 含义是是在方法1的IRC途径上 用方法2在每一个点都做单点能计算,找出其中能量最高的点。  但是CBS-QB3是一个多步骤的方法 其中还包含了第一步 ...

那还有没有其他的高精度的方法可以直接计算过渡态呢?我这个过渡态还是蛮多的。
13楼2019-05-16 11:42:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
rocisky73楼
2019-05-15 12:08   回复  
流浪猫小琪(金币+2): 谢谢参与
发自小木虫IOS客户端
yundaclgc5楼
2019-05-15 12:18   回复  
流浪猫小琪(金币+2): 谢谢参与
发自小木虫Android客户端
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 审稿累死了,都喜欢过年前投稿 +31 lc1984 2020-01-19 44/2200 2020-01-20 00:11 by shen2159
[考博] 拟录取的小纠结 +28 zl4738 2020-01-13 31/1550 2020-01-19 23:59 by 从新开始1
[考博] 生活不易,考非全日制博士也不易! +3 dyp8848 2020-01-14 3/150 2020-01-19 23:42 by enao313
[硕博家园] 谈谈人生谈谈心情,小年文章被拒 +16 133456 2020-01-17 26/1300 2020-01-19 22:47 by 执明
[论文投稿] 爱思唯尔发表的论文不给页码了吗 +7 小弟无知 2020-01-19 12/600 2020-01-19 22:35 by 小弟无知
[基金申请] 有没有已经知道2020年教育部人文社科立项结果的虫友? +14 laomama 2020-01-17 20/1000 2020-01-19 17:15 by Dlliu5858
[版块工场] 天然产物化学谱图解析 +3 不放糖的汤圆 2020-01-17 5/250 2020-01-19 15:48 by 不放糖的汤圆
[第一性原理] CASTEP 源代碼對全球學術界免費公開(官網正式上線) +5 mhslee 2020-01-18 5/250 2020-01-19 11:46 by hxxywwp
[育儿交流] 胚胎停育检查结果——22号染色体三体,以及一些给大龄备孕家庭的建议 (金币+5) +22 spacr 2020-01-17 24/1200 2020-01-19 11:41 by syj525123
[论文投稿] required reviews completed。该状态维持半月了,忐忑啊? +7 tzjs731 2020-01-16 13/650 2020-01-19 09:34 by 天天看书123
[论文投稿] 论文投稿 +4 清新的你 2020-01-18 4/200 2020-01-18 21:32 by 逝水无痕1990
[论文投稿] JACS Au:终于连浓眉大眼的JACS也受不了金钱的诱惑了 +9 Tiriosh 2020-01-17 13/650 2020-01-18 19:47 by Tiriosh
[考博] 金属材料考博 +9 77497749 2020-01-16 17/850 2020-01-18 00:12 by 学习君cb
[硕博家园] 我想请教一个的问题 +9 左亚雪 2020-01-17 19/950 2020-01-17 21:11 by 向宽处行
[基金申请] 谁有教育部社科司的电话? +22 royinwest 2020-01-16 41/2050 2020-01-17 20:49 by royinwest
[论文投稿] 有关论文小修求助 +6 矢车菊2012 2020-01-15 11/550 2020-01-17 06:40 by 矢车菊2012
[硕博家园] 工程博士和一般博士的区别 +4 对夏天过敏 2020-01-13 7/350 2020-01-16 18:11 by 对夏天过敏
[考研] 学科数学求调剂 +4 刘乐乐啊 2020-01-15 6/300 2020-01-15 11:54 by 机车飞驰
[硕博家园] 材料的请看过来 +6 睿黑 2020-01-14 6/300 2020-01-14 12:22 by 天煞将军
[版块工场] 【早起签到贴】2020年1月14日(q q 48626260 进群申请:小木虫——您个人ID) +47 8475 2020-01-14 47/2350 2020-01-14 11:28 by waitingforu
信息提示
请填处理意见