24小时热门版块排行榜    

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

流浪猫小琪

银虫 (正式写手)


[交流] 如何使用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客户端
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 2024还有没有补录的博士 +6 crazy peng 2024-04-30 7/350 2024-05-05 10:17 by 宋mlk
[硕博家园] 大家能不能给点意见啊 +8 是小蚂蚁啦 2024-05-01 15/750 2024-05-05 10:14 by Elio-
[考博] 西安交通大学转化医学研究院蔡文课题组招收2025年入学的申请-考核制博士生 +3 caiwen2016 2024-05-03 3/150 2024-05-05 09:39 by 宋mlk
[教师之家] 博士毕业找工作选择! +9 gsq嘿嘿 2024-04-30 12/600 2024-05-05 09:21 by Nano-Firefly
[教师之家] 七级岗,到手不到8000,带本科论文12教学分/生,额定教学分280 +20 河西夜郎 2024-05-02 23/1150 2024-05-05 08:57 by Quakerbird
[教师之家] 校级优秀本科毕业论文,比例多少?看指导教师资历吗? +3 河西夜郎 2024-05-04 3/150 2024-05-05 08:35 by Xiaolin81
[基金申请] 活页超过7000字,有影响吗? +4 青竹82 2024-05-04 5/250 2024-05-04 23:28 by 华氏911.
[高分子] 水性聚氨酯 +5 钠与水 2024-04-28 10/500 2024-05-04 22:48 by 钠与水
[论文投稿] 审稿人评审完成状态还是under review 5+3 ZH0826 2024-04-30 10/500 2024-05-04 21:05 by topedit
[基金申请] 关于公布2024年度国家自然科学基金项目申请初审结果的通告 (金币+2) +12 hdzw9071 2024-04-30 18/900 2024-05-04 17:49 by felicity6056
[硕博家园] 等盲审 +6 暮色恋伊人 2024-05-03 7/350 2024-05-04 16:51 by shl2112501
[论文投稿] 想问下电池发在nano letters含金量高不高 +5 未来富婆蛙 2024-04-30 10/500 2024-05-04 16:02 by jujuhenin
[考博] 北京师范大学人工智能学院2024年第二批博士研究生招生报名 +4 taoto 2024-05-02 5/250 2024-05-04 15:46 by LYS1200
[硕博家园] 博士后和科研助理好难找啊 +13 找工作啊 2024-04-28 13/650 2024-05-04 14:53 by 白菜1992
[硕博家园] 找工作千万不能来差学校 +8 YN完蛋 2024-05-02 8/400 2024-05-04 12:46 by kanyechris
[考博] 食品专业寻博士导师 +3 里氏什么棱 2024-04-29 6/300 2024-05-04 08:00 by 努力搬砖中呀
[考博] 博士申请 +6 Farssh 2024-04-29 11/550 2024-05-03 18:05 by 头秃食品狗
[论文投稿] 第一篇论文终于进去审稿状态了 +5 lizhengke06 2024-05-01 8/400 2024-05-03 09:14 by nono2009
[基金申请] 多大算大龄? +4 wyjecho666 2024-05-01 5/250 2024-05-02 13:47 by AD_com
[教师之家] 一般发什么文章才能中青基啊 +7 保持肌肤 2024-04-29 7/350 2024-05-01 17:03 by zju2000
信息提示
请填处理意见