计算ECD时出现ORCA finished by error termination in CIS问题

ljc050512

wzkchem5 发表于 2021-11-27 16:32
对，可以试一下。因为第0步迭代用的内存不会完全释放，有一部分内存是即使第0步迭代结束了也还在继续占用 ...

感谢老师的耐心解答。跟您学会了不少知识。

ljc050512

wzkchem5 发表于 2021-11-27 16:32
对，可以试一下。因为第0步迭代用的内存不会完全释放，有一部分内存是即使第0步迭代结束了也还在继续占用 ...

老师您好，我用了之前的泛函和基组，算完之后与文献的紫外光谱差距很大。我想用文献给出泛函和基组再试一下。涉及到写混合基组的事情，我不太清楚，想请教一下。文献中用泛函的是CAM-B3LYP，对Au用的基组是LANL2dz，对C、P、O、Cl和H用的是6-31G(*)。我查阅了咱这个论坛的一些贴子，写了两种设置，想麻烦您给看看。
（一）：
! CAM-B3LYP 6-31G(d) ECP{LanL2=[Au],LANL2DZ} RIJCOSX tightSCF noautostart miniprint nopop
%maxcore  100000
%pal nprocs   2 end
%cpcm
smd true
SMDsolvent "DICHLOROMETHANE"
end
%tddft
nroots 100
TDA false
end
（二）：
! CAM-B3LYP 6-31G(d) RIJCOSX tightSCF noautostart miniprint nopop
%maxcore  100000
%pal nprocs   2 end
%cpcm
smd true
SMDsolvent "DICHLOROMETHANE"
end
%basis
NewGTO Au "LANL2DZ" end # end前面需要写"DelECP Au"吗？
%tddft
nroots 100
TDA false
end
哪一种比较合适？还是二者皆可？
另外， 6-31G(d) 是不是后面加上RIJCOSX 不合适？直接去掉RIJCOSX吗？直接去掉的话，是不是就没有“加速”效果了？还是改基组为def2-SV(P)，然后写成“def2-SV(P) def2/J RIJCOSX”比较合适？

wzkchem5

ljc050512 发表于 2021-12-3 13:03
老师您好，我用了之前的泛函和基组，算完之后与文献的紫外光谱差距很大。我想用文献给出泛函和基组再试一 ...

LANL2DZ和6-31G*都属于比较过时的基组了，而且对于Au这么重的元素，用赝势考虑相对论效应未必靠谱。建议对金属至少用ZORA-def2-TZVP(-f)，并且用ZORA方法显式考虑相对论效应，非金属至少用ZORA-def2-SV(P)

ljc050512

wzkchem5 发表于 2021-12-3 20:13
LANL2DZ和6-31G*都属于比较过时的基组了，而且对于Au这么重的元素，用赝势考虑相对论效应未必靠谱。建议 ...

好的老师，我试试看您推荐的混合基组。我想请教一下，我写的那两个混合机组的格式可以用吗？第一次写混合机组，有点糊涂呢......

wzkchem5

ljc050512 发表于 2021-12-3 14:17
好的老师，我试试看您推荐的混合基组。我想请教一下，我写的那两个混合机组的格式可以用吗？第一次写混合 ...

第一种写法我不确定行不行，第二种需要在%tddft前面再加一个end，因为NewGTO那一行的end对应的是NewGTO，%basis还需要另一个end，此外需要用NewECP（不是DelECP）指定ECP。
类似“这样写行不行”的问题，最简单的办法是自己写完了跑一下，如果输出文件一上来就报错，就证明这样写不对，也省得你发帖问了。如果输出文件没报错，那么仔细地、一行行地读输出文件，验证输出文件的内容是否和你设想的相符，如果还是不确定这样做对不对，再来问。

ljc050512

wzkchem5 发表于 2021-12-3 21:47
第一种写法我不确定行不行，第二种需要在%tddft前面再加一个end，因为NewGTO那一行的end对应的是NewGTO， ...

是的老师，我都去试一下。谢谢老师指点