计算分子的圆锥交叉点的时候出现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

snljty

感觉卢老师讲义上那个乱用Gaussian关键词的例子库可以更新了。。。

杨小狗

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

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

snljty

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

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

杨小狗

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

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

杨小狗

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

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

sobereva

把scrf去了

杨小狗

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

谢谢社长搭救

pwzhou

支持在溶剂中优化交叉点构型的软件据我所知，目前只有Firefly支持，Gaussian只支持在气相中优化，所以必须去掉scrf才可以。

杨小狗

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

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

陈小北

想问一下楼主现在已经成功做出来了吗？

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

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

biogon

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

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

biogon

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

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

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

陈小北

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

好的，了解了

谢谢！

杨小狗

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

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