简单量子化学问题答疑专帖

hehe321123

各位老师好，高斯过渡态的初猜结构可以不优化吗？

limen

sobereva 发表于 2023-7-13 12:17
都没什么比较的意义
B3LYP算热力学量早过时了http://bbs.keinsci.com/thread-12773-1-1.html

社长，能量用这个算“# DSDPBEP86/基组 IOp(3/125=0079905785,3/78=0429604296,3/76=0310006900,3/74=1004) em=gd3bj IOp(3/174=0437700,3/175=-1,3/176=0,3/177=-1,3/178=5500000)”，还是差数量级。

wzkchem5

limen 发表于 2023-7-13 10:14
社长，能量用这个算“# DSDPBEP86/基组 IOp(3/125=0079905785,3/78=0429604296,3/76=0310006900,3/74=100 ...

差多少？别只说差数量级，0.1kcal/mol和1kcal/mol是差一个数量级，10kcal/mol和100kcal/mol也是差一个数量级，但是前者没问题，后者有问题。

wzkchem5

hehe321123 发表于 2023-7-13 10:09
各位老师好，高斯过渡态的初猜结构可以不优化吗？

为什么要不优化？

hehe321123

各位老师好，我想问一下自定义溶剂时，什么情况下考虑非极性部分？什么情况下不考虑

limen

wzkchem5 发表于 2023-7-13 18:14
差多少？别只说差数量级，0.1kcal/mol和1kcal/mol是差一个数量级，10kcal/mol和100kcal/mol也是差一个数 ...

算乙醇，查阅资料乙醇生成焓为278kJ/mol
方案一：-1374kJ/mol
方案二：-260kJ/mol

sobereva

hehe321123 发表于 2023-7-13 21:58
各位老师好，我想问一下自定义溶剂时，什么情况下考虑非极性部分？什么情况下不考虑

只能说考虑非极性部分的必要性，没有什么情况刻意不应当考虑非极性部分。
非极性部分对几何优化、振动分析的影响相对于对算能量问题的影响整体小得多

sobereva

limen 发表于 2023-7-14 08:33
算乙醇，查阅资料乙醇生成焓为278kJ/mol
方案一：-1374kJ/mol
方案二：-260kJ/mol

通常差得非常离谱，基本都是计算存在硬伤，诸如流程不对、数据读错了、计算方式没可比性，以及收敛到不稳定波函数等等

sobereva

hehe321123 发表于 2023-7-13 15:30
各位老师，我想问一下如果过渡态中断以后续算我可以直接用rcfc读取终端的那个过渡态文件的hessin矩阵吗

直接在opt里写个restart文件完了

limen

sobereva 发表于 2023-7-14 11:07
通常差得非常离谱，基本都是计算存在硬伤，诸如流程不对、数据读错了、计算方式没可比性，以及收敛到不稳 ...

谢谢社长

hehe321123

sobereva 发表于 2023-7-14 11:07
直接在opt里写个restart文件完了

好的，谢谢老师

hehe321123

sobereva 发表于 2023-7-14 11:05
只能说考虑非极性部分的必要性，没有什么情况刻意不应当考虑非极性部分。
非极性部分对几何优化、振动分 ...

谢谢老师

糖糖糖豆9988

请教各位老师，自由能垒低的反应用MESS计算出来的高压极限速率常数并不是最快的，这如何解释呢？谢谢

Roonry

各位老师好，我使用M06-2X/6-311G**  计算溴化双酚A的分解 C-Br键会首先断开，Br自由基可能和甲基自由基结合生成溴甲烷，也可能和氢自由基生成溴化氢，但为什么生成溴甲烷路径最终的能量是正的，生成溴化氢是负的？不应该都是负的吗
是自由基结合需要校正吗？
谢谢老师

naive

我理解的使用geom=check guess=read chk文件，那读取的是不是最后一个结果，但是我想要读取的是之前的构象，那要怎么读取我想要的构象呢。
还是选择我需要的那个结构保存位gif文件，再去计算，那保存的gif文件是否会有原来的信息呢