计算化学公社

标题: 求教EOM-CCSD在高斯里的计算方法 [打印本页]

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

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

作者
Author: 静静子 时间: 2022-10-13 18:20

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

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

作者
Author: wzkchem5 时间: 2022-10-13 19:08

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

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

作者
Author: zjxitcc 时间: 2022-10-13 19:13

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

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

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

作者
Author: chrinide 时间: 2022-10-13 20:45

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

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

作者
Author: sobereva 时间: 2022-10-14 21:32

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

这个流程是标准操作

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

作者
Author: 静静子 时间: 2022-10-16 13:57

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后面加上溶剂条件可以嘛？感谢老师~

作者
Author: ionexchangeC 时间: 2022-10-16 15:11
当前基组太差，中性体系弥散函数不如加一级zeta数，换成cc-pVTZ更好

作者
Author: sobereva 时间: 2022-10-16 18:51

静静子发表于 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了
可以带溶剂模型

作者
Author: 静静子 时间: 2022-10-17 00:31

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

好的，感谢老师的回复

作者
Author: 静静子 时间: 2022-10-18 10:31
输入文件：
%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”的含义是什么？

作者
Author: 静静子 时间: 2022-10-18 10:42
输入文件：
%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”代表的含义是什么？本人是从官网的学习实例中看到的，但对此并不了解，特来请教各位老师，下面压缩文件中给出了我的输入及输出报错文件。

作者
Author: 函数与激情 时间: 2024-5-15 11:19

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

作者
Author: zjxitcc 时间: 2024-5-15 15:06
本帖最后由 zjxitcc 于 2024-5-15 15:08 编辑

函数与激情发表于 2024-5-15 11:19
请教下卢老师，我的体系是一个小分子+O2的复合物，基态是三重态，我想做基于基态的垂直激发过程，输入文 ...

这么算没啥意义。EOM-CCSD的计算步骤是先算HF，然后算基态CCSD，最后算EOM。你这基态是三重态，说白了基态是用UHF算，它存在自旋污染，然后UCCSD存在自旋污染，接着EOM得到的每个态都存在自旋污染，到时候S1/T1/S2一个都辨认不出来，不如不算。

TD(triplet)和EOM(triplet)关键词只有在基态是闭壳层单重态的时候有效。你这个例子里写的triplet是不会生效的。你这种体系适合用MRSF-TDDFT做计算，可以参考《Spin-flip方法中RODFT难收敛解决办法》http://bbs.keinsci.com/thread-45194-1-1.html

另外，CC计算耗内存很大，你的内存（单位GB）应该远大于你的核数。

作者
Author: 函数与激情 时间: 2024-5-15 20:30

zjxitcc 发表于 2024-5-15 15:06
这么算没啥意义。EOM-CCSD的计算步骤是先算HF，然后算基态CCSD，最后算EOM。你这基态是三重态，说白了基 ...

好的谢谢！

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