计算化学公社

标题: 计算分子的圆锥交叉点的时候出现L1003报错 [打印本页]

作者
Author: 杨小狗 时间: 2020-5-25 21:16
标题: 计算分子的圆锥交叉点的时候出现L1003报错
PCM not available in this version.
请问各位老师，这是为什么呢？是因为 NOPCMCSCF option不支持SCRF吗？

我的关键词是：
#p opt=conical casscf(10,10,nroot=2,nocpmcscf)/sto-3g
scrf=(smd,solvent=acetonitrile) nosymm guess=read pop=full iop(5/7=128)

报错消息的尾部信息：
MCSCF converged.
Leave Link  510 at Mon May 25 19:59:06 2020, MaxMem= 536870912 cpu:          28463.6 elap:          7116.4
(Enter /apps/chem/Gaussian/G16_A03/AVX2/g16/l801.exe)
Range of M.O.s used for correlation:    1 111
NBasis= 111 NAE= 67 NBE= 67 NFC=    0 NFV=    0
NROrb= 111 NOA= 67 NOB= 67 NVA= 44 NVB= 44
Leave Link  801 at Mon May 25 19:59:11 2020, MaxMem= 536870912 cpu:             0.2 elap:             0.7
(Enter /apps/chem/Gaussian/G16_A03/AVX2/g16/l1101.exe)
Not using compressed storage, NAtomX= 35.
Will process    35 centers per pass.
Leave Link 1101 at Mon May 25 19:59:15 2020, MaxMem= 536870912 cpu:             2.0 elap:             0.5
(Enter /apps/chem/Gaussian/G16_A03/AVX2/g16/l1003.exe)
Gradient Difference/Derivative Coupling Calculation
CP-MCSCF equations will not be solved
NO. OF ORBITALS       =111    NO. OF CORE-ORBITALS = 62
NO. OF VALENCE-ORBITALS = 10    NO. OF VIRTUAL-ORBITALS = 39
PCM not available in this version.
Error termination via Lnk1e in /apps/chem/Gaussian/G16_A03/AVX2/g16/l1003.exe at Mon May 25 19:59:20 2020.
Job cpu time:    0 days  7 hours 54 minutes 37.4 seconds.
Elapsed time:    0 days  1 hours 58 minutes 40.7 seconds.
File lengths (MBytes):  RWF= 230 Int=    0 D2E=    0 Chk=    11 Scr=    1

作者
Author: snljty 时间: 2020-5-26 01:03
感觉卢老师讲义上那个乱用Gaussian关键词的例子库可以更新了。。。

作者
Author: 杨小狗 时间: 2020-5-26 01:10

snljty 发表于 2020-5-26 01:03
感觉卢老师讲义上那个乱用Gaussian关键词的例子库可以更新了。。。

iop（5/7=128）≡ scf(maxcycle=128) nopcmcscf是为了做太平均节省计算时间和所需空间。其他的没问题啊除非scrf中做casscf是错误的。

作者
Author: snljty 时间: 2020-5-26 01:13

杨小狗发表于 2020-5-26 01:10
iop（5/7=128）≡ scf(maxcycle=128) nopcmcscf是为了做太平均节省计算时间和所需空间。其他的没问题啊 ...

http://bbs.keinsci.com/thread-3460-1-1.html

作者
Author: 杨小狗 时间: 2020-5-26 01:16

snljty 发表于 2020-5-26 01:13
http://bbs.keinsci.com/thread-3460-1-1.html

casscf有的时候不太好做，默认的scf圈数是64 我加到128就是想看看ITN是不是有减小的趋势。是不是我贴的错误和scf圈数设置有关啊，谢谢

作者
Author: 杨小狗 时间: 2020-5-26 01:22

snljty 发表于 2020-5-26 01:13
http://bbs.keinsci.com/thread-3460-1-1.html

啊谢谢 pop=full这个我其实真实最开始的时候就用了几次。后来真没用到它。到时候我把这个去掉

作者
Author: sobereva 时间: 2020-5-26 07:22
把scrf去了

作者
Author: 杨小狗 时间: 2020-5-26 09:19

sobereva 发表于 2020-5-26 07:22
把scrf去了

谢谢社长搭救

作者
Author: pwzhou 时间: 2020-5-26 11:10
支持在溶剂中优化交叉点构型的软件据我所知，目前只有Firefly支持，Gaussian只支持在气相中优化，所以必须去掉scrf才可以。

作者
Author: 杨小狗 时间: 2020-5-26 12:03

pwzhou 发表于 2020-5-26 11:10
支持在溶剂中优化交叉点构型的软件据我所知，目前只有Firefly支持，Gaussian只支持在气相中优化，所以必须 ...

感谢老师回复和帮助谢谢！

作者
Author: 陈小北 时间: 2020-12-15 09:56
想问一下楼主现在已经成功做出来了吗？

不知道楼主的体系是不是激发态呢？

如果是的话，想问一下如果想找激发态S1态和S2态的圆锥交叉点应该以哪一个结构为输入文件呢？是基态优化下的还是激发态优化的？

作者
Author: biogon 时间: 2020-12-15 10:44

杨小狗发表于 2020-5-26 01:16
casscf有的时候不太好做，默认的scf圈数是64 我加到128就是想看看ITN是不是有减小的趋势。是不是我贴的错 ...

高斯算这么大活性空间的cas收敛性非常差速度还很慢

作者
Author: biogon 时间: 2020-12-15 10:51

陈小北发表于 2020-12-15 09:56
想问一下楼主现在已经成功做出来了吗？

不知道楼主的体系是不是激发态呢？

建议用S2的极小点结构出发，找到比较接近CI处的结构再优化CI结构

作者
Author: 陈小北 时间: 2020-12-15 16:03

biogon 发表于 2020-12-15 10:51
建议用S2的极小点结构出发，找到比较接近CI处的结构再优化CI结构

好的，了解了

谢谢！

作者
Author: 杨小狗 时间: 2020-12-17 16:44

biogon 发表于 2020-12-15 10:44
高斯算这么大活性空间的cas收敛性非常差速度还很慢

是的，用orca能强一些嘛
谢谢啦

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