24小时热门版块排行榜    

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

流浪猫小琪

银虫 (正式写手)


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

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

» 猜你喜欢

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

查看全部散金贴

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

流浪猫小琪

银虫 (正式写手)


引用回帖:
28楼: Originally posted by paramecium86 at 2019-05-24 01:37:11
这个是因为你用zmatrix来写坐标 确没有给出要优化的变量

谢谢前辈指点。CCSD(T)/aug-cc-pVTZ做IRCmax怎么耗这么多时间,我在超算开60核,已经运行了40个小时了还没有结果,体系只有S1O2H2C六个原子。
29楼2019-05-26 01:21:09
已阅   回复此楼   关注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分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 找工作千万不能来差学校 +9 YN完蛋 2024-05-02 9/450 2024-05-05 13:24 by 溪山游客
[论文投稿] 想问下电池发在nano letters含金量高不高 +6 未来富婆蛙 2024-04-30 11/550 2024-05-05 10:53 by eemdl
[教师之家] 博士毕业找工作选择! +9 gsq嘿嘿 2024-04-30 12/600 2024-05-05 09:21 by Nano-Firefly
[论文投稿] 投稿某一SCI期刊,论文被拒绝,是否可以再次投稿到这个期刊? 8+3 ilyzm 2024-05-03 7/350 2024-05-05 09:20 by bobvan
[考研] 要不要二战? +5 ajjznskk 2024-05-05 5/250 2024-05-05 09:02 by 花满楼Gsir
[基金申请] 活页超过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 钠与水
[论文投稿] 想请问下IEEE Trans最后一轮,你们等了多久? +4 DD4Dr 2024-05-03 5/250 2024-05-04 20:52 by topedit
[考博] 2024申博 化学 一作sci两篇 +8 zzzjjjttt 2024-04-29 8/400 2024-05-04 19:14 by lqmkl
[基金申请] 关于公布2024年度国家自然科学基金项目申请初审结果的通告 (金币+2) +12 hdzw9071 2024-04-30 18/900 2024-05-04 17:49 by felicity6056
[考博] 北京师范大学人工智能学院2024年第二批博士研究生招生报名 +4 taoto 2024-05-02 5/250 2024-05-04 15:46 by LYS1200
[考研] 考研日记 +5 九三hs 2024-05-02 5/250 2024-05-04 14:49 by 鱼翔浅底1
[硕博家园] 小硕盲审没过 +17 xmc2007 2024-04-28 17/850 2024-05-04 09:45 by dong5391
[考博] 食品专业寻博士导师 +3 里氏什么棱 2024-04-29 6/300 2024-05-04 08:00 by 努力搬砖中呀
[硕博家园] 去美国的华人导师底下读博后会很难出站吗? +6 Aries-FZU 2024-05-03 6/300 2024-05-04 07:32 by 晓目崇
[基金申请] 今年面上预计又打酱油了 +15 speedxyl 2024-04-30 17/850 2024-05-04 02:07 by jbb0523
[考博] 博士申请 +6 Farssh 2024-04-29 11/550 2024-05-03 18:05 by 头秃食品狗
[无机/物化] 过氧化氢氧化 10(金币+5) +3 15132328891 2024-04-29 6/300 2024-05-03 13:39 by tan151646
[论文投稿] 第一篇论文终于进去审稿状态了 +5 lizhengke06 2024-05-01 8/400 2024-05-03 09:14 by nono2009
[基金申请] 第二次申请青基,祝自己好运! 中! 中! 中! +13 阿豪v5 2024-04-29 15/750 2024-05-01 14:44 by Howard28
信息提示
请填处理意见