简单量子化学问题答疑专帖（旧版，已关帖）

mfdsrax2 · 发表于 Post on 2021-11-27 21:40:29

sobereva 发表于 2021-11-27 17:30
“连续20次失败”语义不明
IRC第一个点SCF没收敛的话根本不会往后跑

就是校正步不收敛，但是分开算，或者用calcfc就能正常收敛，rcfc就不行

archer · 发表于 Post on 2021-11-27 22:32:35

计算Cl-的吸附可以不考虑弥散吗？弥散计算是否需要很大内存？算几步就出现：

Error termination in NtrErr:
NtrErr Called from FileIO.

sobereva · 发表于 Post on 2021-11-28 14:15:07

archer 发表于 2021-11-27 22:32
计算Cl-的吸附可以不考虑弥散吗？弥散计算是否需要很大内存？算几步就出现：

Error termination in NtrE ...

需要考虑（至少对Cl这个原子而言），毕竟负电荷那么密集

和内存需求量没有什么直接关系，看具体体系、算什么任务

sobereva · 发表于 Post on 2021-11-28 14:16:32

mfdsrax2 发表于 2021-11-27 21:40
就是校正步不收敛，但是分开算，或者用calcfc就能正常收敛，rcfc就不行

甭用rcfc就完了
rcfc读取之前最后的Hessian，倘若你之前的任务不是精确计算TS的Hessian的任务，那读的Hessian根本就不是像calcfc得到的那样严格的

mfdsrax2 · 发表于 Post on 2021-11-28 17:15:50

sobereva 发表于 2021-11-28 14:16
甭用rcfc就完了
rcfc读取之前最后的Hessian，倘若你之前的任务不是精确计算TS的Hessian的任务，那读的He ...

我发现之前的优化用calcall，然后IRC用rcfc就会出问题，如果优化是calcfc+freq，IRC用rcfc就是正确的，calcall之后也要用freq才能在IRC用rcfc吗？

mfdsrax2 · 发表于 Post on 2021-11-28 17:23:21

本帖最后由 mfdsrax2 于 2021-11-28 17:25 编辑

我在用CCSD(T)算能量，体系稍微一大就会卡在最后一步，从任务开始到这里用了2天，但是现在又过去快10天了，进度还是在这里一点没动，RWF文件倒是一直在更新，每分钟都要更新好几次，就是LOG和CHK不动，这是什么问题，还应该继续等下去吗

wzkchem5 · 发表于 Post on 2021-11-28 17:41:54

mfdsrax2 发表于 2021-11-28 10:23
我在用CCSD(T)算能量，体系稍微一大就会卡在最后一步，从任务开始到这里用了2天，但是现在又过去快10天了， ...

换用orca，做DLPNO-CCSD(T)计算。高斯算小体系的CCSD(T)还可以，稍微大一点的体系不要用高斯算CCSD(T)，高斯不是用来干这个的

329678001 · 发表于 Post on 2021-11-28 18:52:59

请教一下各位老师，我的结构有个-10的虚频，我每次想要去掉这个虚频的时候，结构就会出现很大变化。想请教一下是不是加一个去掉对称化的关键词会去掉虚频的同时还不影响几何结构呢？

落雁 · 发表于 Post on 2021-11-28 20:51:02

各位老师好，我在拜读了sob老师《Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法》http://sobereva.com/314 一文吼，我想计算几个分子的自发辐射速率，但是计算出来的S1的振子强度是0，S2和S3不是，请问这种情况是由于有非常强烈的内转换吗？

yiwei · 发表于 Post on 2021-11-29 10:43:05

请问Gaussian中Orca的Int(Grid=99590)对应Orca的什么关键词呀?
我使用ωB97M-V泛函，用Grid5显示如下：
UNRECOGNIZED OR DUPLICATED KEYWORD(S) IN SIMPLE INPUT LINE
GRID5

请问可以怎么解决吗？（使用的是Orca 5.0.1 )

rendong · 发表于 Post on 2021-11-29 10:48:18

本帖最后由 rendong 于 2021-11-29 15:44 编辑

老师好，我在用CAM-B3LYP(full)/gen （混合基组C,O,H用的def2SVP，Cd用的def2TZVP）计算Cd配合物的CT激发态时，错误提示关键词不对，如果不加(full)，又会出现l801冻芯错误，想请问老师我需要如何修改合适。（另外在用CAM-B3LYP/def2SVP 单一机组时没有出现冻芯，我想对Cd提升一下精度）谢谢老师

rendong · 发表于 Post on 2021-11-29 10:48:25

本帖最后由 rendong 于 2021-11-29 15:41 编辑

另外想请教一下老师对于Gaussian程序，计算配合物电荷转移（CT）激发除了CAM-B3LYP还有其他合适的替代泛函么

rendong · 发表于 Post on 2021-11-29 10:50:46

rendong 发表于 2021-11-29 10:48
老师好，我在用cam-b3lyp(full)/gen （混合基组C,O,H用的defSVP，Cd用的def2TZVP）计算Cd配合物的CT激发态 ...

非常抱歉，卡了一下，点了两次提交，

Yaqi · 发表于 Post on 2021-11-29 12:15:01

卢老师您好，最近做了一个主客体分子识别的相关工作，计算体系类似下图等一系列分子：

计算过程大致如下：
1. 使用Moclus进行构象搜索，找到最优复合物构象；（半经验方法：GFN1-xTB 和 GFN2-xTB）
2. 进行单点计算和结合能计算；（优化：B3LYP/6-31G(d)；单点能：B3LYP/6-311G (d, p)）
3. 通过Multiwfn进行了IGM的弱相互作用分析。所有计算加了Grimme’s D3BJ 色散校正和水的隐式溶剂模型。

审稿意见中有两条我有点儿迷惑，不知道从什么方面去回答，请卢老师指点。
相关审稿意见如下：
1. The Gaussian calculations are also ill-advised, as solvent displacement is a strong driver of the selectivity and appears not to be included in the calculations；
2. However, the computational calculations do not have description of the water treatment in any detail, which is essential for binding studies of CBs. “Here we took the solvent effects of water (ε = 78.3553) into account by using the SMD-flavor of self-consistent reaction field (SCRF) theory” is not sufficient, and a more explicit description of entropic effects is essential.

谢谢卢老师。

wzkchem5 · 发表于 Post on 2021-11-29 16:29:23

Yaqi 发表于 2021-11-29 05:15
卢老师您好，最近做了一个主客体分子识别的相关工作，计算体系类似下图等一系列分子：

计算过程大致如下 ...

就引SMD原文献，说SMD模型已经包括了溶剂排阻效应和熵效应就行了，然后把SMD原文献的相关内容paraphrase一下，说具体是怎么考虑的

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 简单量子化学问题答疑专帖（旧版，已关帖）