lym563557413 发表于 2022-8-24 09:09 提醒一下，量子化学软件处理基函数线性相关最常用的阈值是sthresh 1e-6，这也是Gaussian和GAMESS的默认阈值，不建议使用sthresh 1e-7。如果想100%确认是否存在基函数线性相关，可以用一样的结构和基组在Gaussian中算几秒，看轨道数和基函数数目即可下定论。 |
lym563557413 发表于 2022-8-29 17:59 下结论为时过早。先确保两个计算的坐标一样，然后确保基组一样，然后确保ORCA中使用了TightSCF或VeryTightSCF，并且使用了TightPNO，再然后检查两个程序的HF能量是否一致（相差多少），这些因素都排除之后，来这里说一下你的检查结论，然后再提问。 |
本帖最后由 zjxitcc 于 2024-2-21 15:12 编辑 yuyangzap 发表于 2022-8-26 10:43 不建议。推荐不要尝试各种关键词（很多时候看运气收敛，人工成本过高），而是直接采用别的程序（如Gaussian）算好了，甚至检验好了波函数稳定性，然后用fch2mkl小程序传轨道给ORCA。如果程序花500圈收敛了，它敢收敛我还不敢信。 |
lym563557413 发表于 2022-8-29 17:59 这年头别用MP2，要不用DLPNO-CCSD(T)，要不就用SAPT。 |
yuyangzap 发表于 2022-8-26 10:43 老师我还有个问题请问一下，我用orca计算DLPNO-MP2/aTZ基组下的分子间相互作用能是-18.44kcal/mol，跟GAUSSIAN计算的MP2/aTZ的-11.46相比，相差好多，请问关于这种问题，可能是方法不同造成的吗，我的体系是一个堆积体系，我把结构图放下面老师，请老师指点 |
lym563557413 发表于 2022-8-24 09:09 老师，我这里现在不是读取，是正常计算DLmp2/atz基组下的单点能，可是算了好几天了还没好，现在是这样子的，我把截图放下面，请老师看看，这个趋势还能继续等吗 |
| 算半天搞不定，就别读取，直接加大迭代，用slowconv算法 |
wzkchem5 发表于 2022-8-24 01:45 好的老师，多谢解答，我先按照您说的降低sthresh参数试一下 |
lym563557413 发表于 2022-8-23 15:04 这也是典型的基组线性相关问题。凡是能量显著低于合理值（比如比其他基组下正常计算的能量低100Hartree以上），基本都是基组线性相关，要靠提高sthresh以及降低计算数值误差（把积分阈值设严、增加格点）的方法解决 |
lym563557413 发表于 2022-8-23 15:01 %scf sthresh 1e-7 end |
lym563557413 发表于 2022-8-23 22:01 顺便说一下，我这个本来是想直接算atz基组下的单点能，也是因为直接算能量趋势不对，所以才想着通过读取tz的波函数来计算atz的能量，当时直接计算atz基组下的问题我也把图贴下面一下老师 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-1-25 23:25 , Processed in 0.176996 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites