求教EOM-CCSD在高斯里的计算方法

静静子

各位老师好，本人最近学习sob老师的帖子【新提醒】乱谈激发态的计算方法 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)，“EOM-CCSD：使用很普遍，Gaussian就有。单参考态体系结果不错，对于单电子激发误差<0.3eV，强于TDDFT。对多参考态体系不好，双电子激发特征明显的体系误差也较大。略昂贵，能算的体系大小远小于TDDFT能处理的尺寸。”本人想要学习EOM-CCSD的计算方法，但不知道如何在高斯里去具体操作，请各位老师指教~

wzkchem5

参见http://gaussian.com/eom/
另外，如果要算的体系稍微大一些（比如20~50个原子），建议用ORCA的DLPNO-STEOM-CCSD，精度和EOM-CCSD差不多，但是算得快得多

静静子

wzkchem5 发表于 2022-10-13 17:29
参见http://gaussian.com/eom/
另外，如果要算的体系稍微大一些（比如20~50个原子），建议用ORCA的DLPNO-S ...

想请问下老师，这个EOM-CCSD是需要花钱买的是嘛？

wzkchem5

静静子 发表于 2022-10-13 11:20
想请问下老师，这个EOM-CCSD是需要花钱买的是嘛？

你是想问orca吗？不需要，是免费的，在https://orcaforum.kofo.mpg.de/index.php注册账号以后免费下载

zjxitcc

静静子 发表于 2022-10-13 18:20
想请问下老师，这个EOM-CCSD是需要花钱买的是嘛？

EOM-CCSD是方法，ORCA是免费的程序。没有叫EOM-CCSD的程序

静静子

刚接触EOM-CCSD确实很多地方不了解，想请教各位老师可以先用DFT理论获得分子的基态优化结构，然后基于此结构采用EOM-CCSD方法去计算分子的激发能嘛？此外，想问下这里该如何调用关键词去计算singlet和triplet的激发能？
%Chk=my_eom
# EOMCCSD(NStates=10,EnergyOnly)/Aug-CC-PVDZ

chrinide

wzkchem5 发表于 2022-10-13 19:08
你是想问orca吗？不需要，是免费的，在https://orcaforum.kofo.mpg.de/index.php注册账号以后免费下载

话说 现在不挂梯子见不到ORCA了，就算挂了梯子最近也不稳定

sobereva

静静子 发表于 2022-10-13 19:20
刚接触EOM-CCSD确实很多地方不了解，想请教各位老师可以先用DFT理论获得分子的基态优化结构，然后基于此结 ...

这个流程是标准操作

EOMCCSD里写singlet（默认）、triplet控制要求解的激发态的自旋多重度

静静子

sobereva 发表于 2022-10-14 21:32
这个流程是标准操作

EOMCCSD里写singlet（默认）、triplet控制要求解的激发态的自旋多重度

计算S0→S1激发能：
%Chk=my_eom
# EOMCCSD(singlet,NStates=10,EnergyOnly)/Aug-CC-PVDZ

[No Title]

0 1
计算S0→T1激发能：
%Chk=my_eom
# EOMCCSD(triplet,NStates=10,EnergyOnly)/Aug-CC-PVDZ

[No Title]

0 1
麻烦sob老师帮忙看一下这里关键词编辑是否合适？如果考虑溶剂环境，在Aug-CC-PVDZ后面加上溶剂条件可以嘛？感谢老师~

ionexchangeC

当前基组太差，中性体系弥散函数不如加一级zeta数，换成cc-pVTZ更好

sobereva

静静子 发表于 2022-10-16 13:57
计算S0→S1激发能：
%Chk=my_eom
# EOMCCSD(singlet,NStates=10,EnergyOnly)/Aug-CC-PVDZ

只算S1、T1而已，一般没必要算10个态，会造成大量额外不必要的耗时
通常只有里德堡激发、阴离子体系才需要带弥散函数，如10L所说，此时远不如用cc-pVTZ。DZ档次的基组对于EOM-CCSD太low了
可以带溶剂模型

静静子

sobereva 发表于 2022-10-16 18:51
只算S1、T1而已，一般没必要算10个态，会造成大量额外不必要的耗时
通常只有里德堡激发、阴离子体系才需 ...

好的，感谢老师的回复

静静子

输入文件：
%chk=/scratch/QAO-PhDPA-s.chk
# EOMCCSD(singlet,NStates=5,EnergyOnly)/cc-pVTZ scrf=(solvent=toluene)

[No Title]

0 1

输出报错文件：
**** Warning!!: The largest alpha MO coefficient is  0.23429492D+02

MOs do not have abelian symmetry and doing EOM, so symmetry turned off.
Semi-Direct transformation.
ModeAB=           4 MOrb=            99 LenV=   32158300557
LASXX= 178959448206 LTotXX=178959448206 LenRXX=358267399281
LTotAB=179307951075 MaxLAS=132631098138 LenRXY=           0
NonZer=537226847487 LenScr=************ LnRSAI=132631098138
LnScr1=266302442496 LExtra=           0 Total= ************
MaxDsk=          -1 SrtSym=           T ITran=            4
JobTyp=0 Pass  1:  I=   1 to  99.
Erroneous write. Write -1 instead of 2994280.
fd = 4
orig len = 8384512 left = 2994280
g_write
想请各位老师帮忙指点一下，感激不尽~
还想问下输入这个关键词“EnergyOnly”的含义是什么？

静静子

输入文件：
%chk=/scratch/QAO-PhDPA-s.chk
# EOMCCSD(singlet,NStates=5,EnergyOnly)/cc-pVTZ scrf=(solvent=toluene)

[No Title]

0 1

输出报错文件：
**** Warning!!: The largest alpha MO coefficient is  0.23429492D+02

MOs do not have abelian symmetry and doing EOM, so symmetry turned off.
Semi-Direct transformation.
ModeAB=           4 MOrb=            99 LenV=   32158300557
LASXX= 178959448206 LTotXX=178959448206 LenRXX=358267399281
LTotAB=179307951075 MaxLAS=132631098138 LenRXY=           0
NonZer=537226847487 LenScr=************ LnRSAI=132631098138
LnScr1=266302442496 LExtra=           0 Total= ************
MaxDsk=          -1 SrtSym=           T ITran=            4
JobTyp=0 Pass  1:  I=   1 to  99.
Erroneous write. Write -1 instead of 2994280.
fd = 4
orig len = 8384512 left = 2994280
g_write
想请各位老师帮忙指点一下，不知该如何处理？还有个小问题：输入文件的关键词“EnergyOnly”代表的含义是什么？本人是从官网的学习实例中看到的，但对此并不了解，特来请教各位老师，下面压缩文件中给出了我的输入及输出报错文件。

函数与激情

sobereva 发表于 2022-10-16 18:51
只算S1、T1而已，一般没必要算10个态，会造成大量额外不必要的耗时
通常只有里德堡激发、阴离子体系才需 ...

请教下卢老师，我的体系是一个小分子+O2的复合物，基态是三重态，我想做基于基态的垂直激发过程，输入文件是按照下面吗？谢谢！
%chk=t1.chk
%nprocshared=16
%mem=8GB
#p EOMCCSD(nstates=5,triplet)/cc-pvtz

Title Card Required

0 3