计算化学公社

标题: 请教计算激发态想得到吸收光谱，报错l801原因 [打印本页]

作者
Author: 运气加持 时间: 2021-7-19 15:51
标题: 请教计算激发态想得到吸收光谱，报错l801原因
新手计算激发态，分别在PBE0,B3LYP泛函下计算吸收光谱，同样的基态结构，在PBE0下计算出现错误【Excessive mixing of frozen core and valence orbitals. Error termination via Lnk1e in /home/share//g09/l801.exe 】，在B3LYP时正常结束，请问这个错误是什么造成的？谢谢。
输入文件关键词分别为#p pbe1pbe/lanl2dz em=gd3bj TD(nstates=800) ，以及#p B3LYP/lanl2dz em=gd3bj TD(nstates=800)

作者
Author: tiandikuoyuan 时间: 2021-7-19 20:17
计算的什么体系，需要计算800个激发态？

作者
Author: zjxitcc 时间: 2021-7-19 20:22
本帖最后由 zjxitcc 于 2021-7-19 20:58 编辑

哪有直接对非金属、金属全部使用LANL2DZ基组的。看sob老师博文
《谈谈量子化学中基组的选择》http://sobereva.com/336
《谈谈赝势基组的选用》http://sobereva.com/373

作者
Author: sobereva 时间: 2021-7-19 22:40
你这是瞎算

lanl2DZ对前四周期主族几乎没意义
该用什么仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

算800个态也不知道你想做什么。对于一般大小的体系，算800个态根本没意义，都涉及到里德堡态了，还拿lanl2DZ算根本就不可思议。
把此文以及里面提到的其它博文看了了解点常识性知识再计算
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

作者
Author: 运气加持 时间: 2021-7-20 14:27

sobereva 发表于 2021-7-19 22:40
你这是瞎算

lanl2DZ对前四周期主族几乎没意义

谢谢老师，提问的时候我没讲清楚，计算涉及的体系较大，有近100个原子，包括H,C,N,Pb,I,Cs等。算这个是因为我看到一篇文章体系与我的相似，是用PBE(B3LYP)-LANL2DZ来算的激发态继而得到吸收光谱（它涉及的元素有H,C,N,Pb,I），算了500个态，为了涵盖整个光谱范围【文章相关信息在下面的附件以及截图】。我的800确实设置的过大了，因为一开始觉得我的体系很大，现在发现500个态足够了。我现在的问题是，1、对于这样的体系要定义一个混合基组吗？2、如果不定义混合基组，只是想用同样的计算方法去比较和文章中相似的体系的结果，输出文件的报错信息是什么问题呢【现在发现相似体系的计算，除了PBE0泛函会报错，B3LYP有时也会有同样的报错】？

作者
Author: sobereva 时间: 2021-7-20 18:53

运气加持发表于 2021-7-20 14:27
谢谢老师，提问的时候我没讲清楚，计算涉及的体系较大，有近100个原子，包括H,C,N,Pb,I,Cs等。算这个是因 ...

近100个原子根本不算多大，现在较好的计算资源下TDDFT都能算几百个原子

你这种体系，图结果基本靠谱，又不希望昂贵的话，首选是def2-SV(P)，Gaussian里写个def2SV就完事，都不用自定义基组或混合基组。如果体系里有许多带显著负电荷的，那些原子可以改为带弥散的ma-def2-SVP

先用了得当的基组再说别的。

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