24小时热门版块排行榜    

查看: 2750  |  回复: 28
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

流浪猫小琪

银虫 (正式写手)


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

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

» 猜你喜欢

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
25楼: Originally posted by 流浪猫小琪 at 2019-05-22 12:28:35
前辈,这是今天的IRCmax的输入输出文件。我B2LYP/6-311++(d,p)做的测试。

input文件
%chk=RE-W12IRCmax.chk
SCF=(MaxCycle=100,xqc)
#p IRCmax(B3LYP/6-311++G(d,p):B3LYP/cc-pVTZ)
NoSymm

RE-W8-W ...

嗯 可以直接从结构开始 。那是提示的意思不是说要你一定得读取chk  但是和IRC一样 是必须要规定如何计算FC 也就是 如果从chk读取FC 就要计算FC,用CalcFC命令。 或者是每一个irc点都重新计算一次FC,那就是用calcall关键词。 总之因为要计算Force constant 所以必须加入一个命令去调控去计算FC。所以你可以写成 #p IRCmax(calcfc, B3LYP/6-311++G(d,p):B3LYP/cc-pVTZ) NoSymm
26楼2019-05-23 07:52:40
已阅   回复此楼   关注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客户端
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +7 调剂的考研学生 2026-03-09 7/350 2026-03-15 00:58 by 秋有木北
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 308求调剂 +5 是Lupa啊 2026-03-08 8/400 2026-03-14 04:25 by JourneyLucky
[考研] 一志愿中国科学技术大学,086000,生物与医药,求调剂 +3 12341234789 2026-03-08 4/200 2026-03-14 03:14 by JourneyLucky
[考研] 理学化学调剂 +3 040517FYfy? 2026-03-09 3/150 2026-03-14 03:09 by JourneyLucky
[考研] 296求调剂 +5 Xinyu Wu311 2026-03-09 5/250 2026-03-14 03:05 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
信息提示
请填处理意见