用EOMCCSD计算激发能的多重度设置和报错问题

雨飞小雁塔

想用Eom-CCSD计算激发能，Fe2+(H2O)6-Fe3+(H2O)6   六水亚铁和六水铁的复合物的多重度应该怎么设置？输出报下面的错是怎么回事？谢谢解答
输入
%chk=ferric-ferrous-eom
%mem=20gb
%nproc=28
# EOMCCSD(energyonly)/gen pseudo=read

ferric-ferrous-eom

5 10

输出
Iteration Nr.  50
**********************
DD1Dir will call FoFMem   1 times, MxPair=      8514
NAB=  2842 NAA=  1653 NBB=  1176.
DE(Corr)=  -2.6256660     E(CORR)=     -1158.5773215     Delta=-6.05D-11
NORM(A)=   0.13247837D+01
*************
*MAX. CYCLES*
*************
Largest amplitude= 7.29D-02
Error termination via Lnk1e in /home/jjgu/g09E01/g09/l913.exe at Thu Mar  7 01:53:36 2019.

sobereva

注意求助帖必须在帖子标题明确体现出“提问”、“求助”及具体内容（见http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的标题“用EOMCCSD计算激发能”改了，以后请务必注意

尽量直接把输入文件上传，特别是你的关键词里包含了从其它地方读入信息的情况

假设两个位点的Fe的自旋方向都相同，设10重态

那是耦合簇迭代求解不收敛。先确保基态结构已经充分优化了，然后尝试用不同的基组

雨飞小雁塔

sobereva 发表于 2019-3-11 03:30
注意求助帖必须在帖子标题明确体现出“提问”、“求助”及具体内容（见http://bbs.keinsci.com/thread-9348 ...

好的，谢谢老师。基态结构我是用cam-B3lyp去分别优化六水亚铁和六水铁，然后把它们放在一起。不知老师说的基态结构充分优化是指用CCSD去优化复合物吗？

sobereva

雨飞小雁塔 发表于 2019-3-11 09:31
好的，谢谢老师。基态结构我是用cam-B3lyp去分别优化六水亚铁和六水铁，然后把它们放在一起。不知老师说 ...

CCSD优化基态一般来说是极其不合适的做法，又贵又没什么好处，对你当前的情况根本算不动
分别优化然后人工摆在一起得到的结构并没有物理意义，要算激发能也得是对复合物先优化。
CAM-B3LYP绝对不适合优化基态，尤其是配合物，要弄清楚CAM-B3LYP到底适合干什么，仔细看
量子化学研究中切换泛函应当注意的问题
http://sobereva.com/415（http://bbs.keinsci.com/thread-9878-1-1.html）