大分子体系结构优化后计算bsse出错

liuhui1118 · 发表于 Post on 2023-8-29 08:39:58

本帖最后由 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是出现了这样的问题。

zjxitcc · 发表于 Post on 2023-8-29 09:04:05

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

liuhui1118 · 发表于 Post on 2023-8-29 16:36:41

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

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

。我想简单了。感谢

liuhui1118 · 发表于 Post on 2023-8-29 17:06:08

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

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

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

sobereva · 发表于 Post on 2023-8-30 01:48:37

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

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

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

liuhui1118 · 发表于 Post on 2023-8-30 08:29:08

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

明白了。谢谢老师

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

[Gaussian/gview] 大分子体系结构优化后计算bsse出错