计算化学公社

标题: 大分子体系结构优化后计算bsse出错 [打印本页]

作者
Author: liuhui1118 时间: 2023-8-29 08:39
标题: 大分子体系结构优化后计算bsse出错
本帖最后由 liuhui1118 于 2023-8-29 08:41 编辑

报错信息如下：
Petite list used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on          energy=1.00D-06.
No special actions if energy rises.
Density matrix is not changing but DIIS error= 2.25D-06 CofLast=-2.28D-02.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB3LYP) =  -150.310088705    A.U. after  129 cycles
         NFock=128  Conv=0.76D-06    -V/T= 2.0032
Convergence failure -- run terminated.
Error termination via Lnk1e in /gdata/lun5/hlliu/g09/l502.exe at Mon Aug 28 20:02:10 2023.
Job cpu time:    7 days  4 hours 28 minutes 23.4 seconds.
File lengths (MBytes):  RWF= 1222 Int=    0 D2E=    0 Chk=    30 Scr=    2

90个原子的大分子结构
之前进行频率优化时，也总是收敛失败，但是某一次突然就成功了。然后计算bsse是出现了这样的问题。

作者
Author: zjxitcc 时间: 2023-8-29 09:04
你这里面有氧气O2吗？如果有，那你这用的是限制性闭壳层计算（RB3LYP），计算没意义，至于报错与否就不用看了。单重态O2（或含单重态O2的）体系要用对称破缺UDFT来算，建议阅读
《用Gaussian做UHF计算》
《谈谈Gaussian软件中的guess=mix》
《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82
之前的各种计算也要检查、如果发现问题就纠正。

作者
Author: liuhui1118 时间: 2023-8-29 16:36

zjxitcc 发表于 2023-8-29 09:04
你这里面有氧气O2吗？如果有，那你这用的是限制性闭壳层计算（RB3LYP），计算没意义，至于报错与否就不用 ...

明白了。原来这样弄是错的

。我想简单了。感谢

作者
Author: liuhui1118 时间: 2023-8-29 17:06

liuhui1118 发表于 2023-8-29 16:36
明白了。原来这样弄是错的。我想简单了。感谢

看了一下，不知道我的理解对不对。大分子单体那边不能用RB3LYP计算？

用UB3LYP/6-311+g**方法优化氧气分子和单分子单体，再加其他的一些关键词。
然后再将单体和氧气分子放到一块，把两个片段的电荷和多重度分开写，进行优化。
这样理解对吗？

作者
Author: sobereva 时间: 2023-8-30 01:48

liuhui1118 发表于 2023-8-29 17:06
看了一下，不知道我的理解对不对。大分子单体那边不能用RB3LYP计算？

用UB3LYP/6-311+g**方法优化氧气 ...

氧气基态是三重态，如果那个大分子是单重态的，二者物理吸附的复合物也是三重态（除非你研究单重态氧气的情况，那是另一码事，即便这样也不是当闭壳层算）

作者
Author: liuhui1118 时间: 2023-8-30 08:29

sobereva 发表于 2023-8-30 01:48
氧气基态是三重态，如果那个大分子是单重态的，二者物理吸附的复合物也是三重态（除非你研究单重态氧气的 ...

明白了。谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)