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

流浪猫小琪

银虫 (正式写手)


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

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

» 猜你喜欢

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
这个是因为你用zmatrix来写坐标 确没有给出要优化的变量

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

流浪猫小琪(金币+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分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 不知道自己为啥读博,只想快点可以回家! +15 5253709 2021-01-10 20/1000 2021-01-16 18:21 by kinghuatan
[基金申请] 2020年博后中青基,去新单位还能申面上吗? +6 guoguo215 2021-01-16 15/750 2021-01-16 17:37 by guoguo215
[有机交流] Suzuki偶联反应的小疑问 +3 化学小白ii 2021-01-14 3/150 2021-01-16 17:30 by 浓缩铀
[硕博家园] 博士负能量 +7 山之巅322 2021-01-15 11/550 2021-01-16 17:28 by balapili
[考研] 辜负了阳光 +3 考研女孩人 2021-01-14 3/150 2021-01-16 17:18 by 无人的
[教师之家] 文科offer求助:厦大和暨南大学助理教授,怎么选 +9 伐克奥芙 2021-01-14 11/550 2021-01-16 16:51 by wht3198
[教师之家] 迟到的2020年总结 +12 zhousl 2021-01-13 12/600 2021-01-16 16:14 by balapili
[基金申请] 2021指南发布 +5 forever411 2021-01-15 7/350 2021-01-16 12:37 by zhuzg0628
[硕博家园] 读博申请 +3 文文清清 2021-01-16 3/150 2021-01-16 10:58 by 傻猪咪
[找工作] 知道怎么选了 +15 wmyislord 2021-01-14 26/1300 2021-01-16 10:55 by terminato2
[找工作] 合同制教师离职求助!!! +6 cclz 2021-01-15 8/400 2021-01-16 09:25 by xianxin81
[论文投稿] 中科院分区升级版 +9 utaz3582 2021-01-10 13/650 2021-01-16 07:24 by aaa45688123
[分子生物] 核酸检测实验中如何防范气溶胶污染 +3 小小书虫飞 2021-01-14 3/150 2021-01-16 04:20 by junjieshao
[文学芳草园] 偶尔脆弱 +7 海阔天空skn 2021-01-13 7/350 2021-01-16 00:24 by Fulfill
[基金申请] 周五了,北自然会不会公布呢? +3 chhq0102 2021-01-15 4/200 2021-01-15 18:00 by chhq0102
[论文投稿] 这个帖子已不存在 +8 huyingjie8 2021-01-13 14/700 2021-01-14 09:11 by 海宝心声
[论文投稿] 文章大修后一个月状态从Decision in process变成Reviewer Assignment Pending? 5+3 佩奇拱白菜 2021-01-13 4/200 2021-01-14 08:05 by 国际科学编辑
[博后之家] 国家纳米科学中心高薪诚聘人工光合作用特别研究助理 +3 pwyqe 2021-01-09 3/150 2021-01-13 22:41 by asdfsh
[考博] 教育类 思政类 读博 +3 JssCoral 2021-01-13 4/200 2021-01-13 19:03 by berry1984
[电化学] 为什么锂离子电池集流体负极要用铜箔,正极要用铝箔? +4 MJtongshao 2021-01-09 6/300 2021-01-11 16:58 by xmcssj
信息提示
请填处理意见