“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register

求助优化和单点算出的选择性相反问题解决

查看数: 132 | 评论数: 9 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2026-4-22 11:46

正文摘要:

本帖最后由 Jiang127127 于 2026-4-22 20:50 编辑 请问各位老师,在计算反应的中自由优化算出来的选择性是合理的,选用b3lyp-d3计算方法,6-31G**基组 但是在进行单点计算后,选择性反转了,还高出6Kcal/mol,选 ...

回复 Reply

sobereva 发表于 Post on yesterday 18:58
Jiang127127 发表于 2026-4-23 16:12
感谢Sob老师的悉心指导和无私帮助!因为我用oniom分层计算,我尝试用您给出的方法优化,发现两个过渡态结 ...

没具体文件没法说
如果高层原子太多,肯定非常昂贵。即便作为低层的半经验级别很便宜,如果低层特别大也照样跑不动。并且隐式溶剂模型会使得原本耗时很低的半经验方法耗时不再那么低。
Jiang127127 发表于 Post on yesterday 16:12
sobereva 发表于 2026-4-22 14:45
信息不充分,没法确切回答
一个带溶剂模型一个没带溶剂模型、一个大基组一个小基组、并且泛函特征截然不同 ...

感谢Sob老师的悉心指导和无私帮助!因为我用oniom分层计算,我尝试用您给出的方法优化,发现两个过渡态结构仅跑了6步就已经占据约50G内存,耗时为6小时,我的命令如下:#p opt=(modredundant,noeigen) freq oniom(wB97XD/6-311G*:pm6d3) scrf=(IEFPCM, Solvent=H2O);请问老师对于这种情况是不是要换种方法呀?非常抱歉之前没有指出oniom分层计算,给老师添麻烦了,恳请老师再帮忙看一下问题如何解决,谢谢老师!
sobereva 发表于 Post on yesterday 01:14
建议wB97XD/6-311G*或def-TZVP结合IEFPCM优化和振动分析,wB97XD/def2-TZVP结合SMD算高精度单点

没事甭用M06
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html


Jiang127127 发表于 Post on 前天 20:55
sobereva 发表于 2026-4-22 14:45
信息不充分,没法确切回答
一个带溶剂模型一个没带溶剂模型、一个大基组一个小基组、并且泛函特征截然不同 ...

谢谢sob老师的悉心指导,反应的详细信息补充上啦,辛苦老师再帮忙指导一下问题解决,谢谢老师的耐心纠错!我以后一定以更加严谨的态度对待科研!不滥用关键词!
Jiang127127 发表于 Post on 前天 20:52
阮文升-河内国大 发表于 2026-4-22 12:14
你应该说清楚这是什么反应,给出反应物结构,结构保密的话也要尽量详细地把能告诉我们的信息都告诉我们,否 ...

谢谢老师的耐心指导!已经给出反应相关信息,辛苦老师再帮助看一下问题出在哪里,谢谢老师!
sobereva 发表于 Post on 前天 14:45
信息不充分,没法确切回答
一个带溶剂模型一个没带溶剂模型、一个大基组一个小基组、并且泛函特征截然不同,哪怕其中一个因素都可能对能量相对高低产生定性影响。在不给出体系非常具体信息的情况下,只能让你每次换一个因素对照结果判断情况。并且没有体系的具体信息,也没法直接告诉你该用什么级别做优化和单点最合适从而免得你乱试。
并且算单点时如果没收敛到几何优化最后一帧对应的电子结构状态,也使得结果没有可比性

没事甭写毫无意义的noraman,认真看下文
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331http://bbs.keinsci.com/thread-3460-1-1.html
Jiang127127 发表于 Post on 前天 14:37
Luciferly 发表于 2026-4-22 12:35
你算自由能的时候加校正量了吗

你好,我的命令是#p opt=noeigen freq=noraman b3lyp/genecp empiricaldispersion=gd3bj
请问您说的校正量是指温度或者熵的校正吗?
Luciferly 发表于 Post on 前天 12:35
你算自由能的时候加校正量了吗
阮文升-河内国大 发表于 Post on 前天 12:14
你应该说清楚这是什么反应,给出反应物结构,结构保密的话也要尽量详细地把能告诉我们的信息都告诉我们,否则根本没法回答

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-4-24 05:43 , Processed in 0.167056 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list