如何在ORCA中用ADC2方法计算S1和T1的激发能？

暖空 · 发表于 Post on 2021-8-29 10:42:15

本帖最后由暖空于 2021-8-29 10:47 编辑

大家好！
请问如何在ORCA中用ADC2方法计算S1和T1态的激发能？

交代下背景，这是实验组的一篇文章，里面的理论是用dft计算的，但是在审稿意见中，审稿人提到有篇理论文献中的分子和他们的类似（文章见附件），那篇理论文章是在TURBOMOLE中用SCS-CC2方法计算的，遂抨击了一番DFT，然后让他们用TURBOMOLE计算。但是他们和审稿人说明了没有TURBOMOLE的版权，审稿人又建议他们在ORCA中用ADC2方法计算：

The authors acknowledge the value of coupled cluster calculations but note that they do not have the Turbomole software. An alternative are ADC(2) calculations, which could be run using ORCA which is free in lieu of SCS-CC2. Note that DFT does a very poor job, not only in predicting excited states but even the order of the HOMO levels is not reproduced. This puts into doubt all of the calculations. The authors must address this weakness in the manuscript.
实验组就找我帮忙。我搜索了一番，发现ADC2只在今年新出的ORCA5.0以后的版本才实现，网上相关的内容也比较少，手册中关于ADC2的描述也不多。虽然检索和咨询后发现ADC2计算激发能的精确度甚至还不如双杂化泛函，并且较为耗时，一般只能计算50个原子以内的分子。和实验组的老师沟通以后对方还是想先按照审稿人的意见在ORCA中先用ADC2计算一下S1态和T1态的激发能，如果数据不行的话再试试双杂化泛函，到时候也好给审稿人回复。

ORCA 5.0.1的手册上关于ADC2的描述不多（8.5.5 Excited States with ADC2），计算只给了一个算例：

我的理解是这里只计算了9个单重激发态。

请问我的理解是否正确？如果正确的话如何用ADC2计算T1态的激发能呢？

谢谢大家！

Pershin-2019-Highly emissive excitons with red.pdf (814.14 KB, 下载次数 Times of downloads: 35)

beefly · 发表于 Post on 2024-12-6 18:14:52

wal 发表于 2024-12-3 14:32
诶？不是只能算单重态嘛？我看到老师加了dotriplet true
另外做eom-ccsd，为什么关键词用的还是steom呢

也可以算三重态。

steom考虑了一些三激发，对于电荷转移态有改进，比如我之前研究过的有机分子体系，eom-ccsd的激发能是6 eV以上，steom-ccsd是2 eV，接近实验值。但是steom容易出现活性空间激发振幅小于98%的问题，导致结果不可靠，需要反复调试参数。如果不是特殊体系，改用eom-ccsd就可以了，输入差不多

wal · 发表于 Post on 2024-12-3 14:32:25

beefly 发表于 2024-12-3 12:01

诶？不是只能算单重态嘛？我看到老师加了dotriplet true
另外做eom-ccsd，为什么关键词用的还是steom呢

beefly · 发表于 Post on 2024-12-3 12:01:08

wal 发表于 2024-11-27 17:53
求问一下老师，EOM-DLPNO-CCSD的输入是怎么写的，手册上没有查到
直接把EOM-CCSD的输入抄过来加上DLPNO ...

! rhf steom-dlpno-ccsd def2-tzvp(-f) def2-tzvp/c def2/j tightscf cpcm
%maxcore 256000
%pal nprocs 4 end
%cpcm
smd true
smdsolvent "dichloromethane"
end
%mdci
nroots 1
dotriplet true
dosolv true
othresh 0.0001
vthresh 0.0001
doeommp2 true
dlpnolinear true
newdomains true
docosxeom true
doaox3e true
kcopt kc_aoblas
maxiter 500
dorootwise true
end

复制代码

wal · 发表于 Post on 2024-11-27 17:53:18

beefly 发表于 2021-9-1 19:22
orc可以做eom-dlpno-ccsd。我用def2-tzvpp(-f)基组配200gb内存算50-60原子没问题，不过目前能算单态

求问一下老师，EOM-DLPNO-CCSD的输入是怎么写的，手册上没有查到
直接把EOM-CCSD的输入抄过来加上DLPNO相关的辅助基组之类的能行吗

biogon · 发表于 Post on 2021-9-28 10:07:32

beefly 发表于 2021-9-1 19:22
orc可以做eom-dlpno-ccsd。我用def2-tzvpp(-f)基组配200gb内存算50-60原子没问题，不过目前能算单态

我又仔细看了下手册，那个是IP and EA-EOM-DLPNO-CCSD，这个也能用来算一般性的激发态么？

biogon · 发表于 Post on 2021-9-2 18:13:27

beefly 发表于 2021-9-2 13:21
花两万块钱买个组装机，再扩一万块钱的内存（给orca的ice-ci预备的，eom用不了这么多），只能一个任务独 ...

这个不吃硬盘？我跑STEOM都要用8T SSD 1T内存的机器跑

beefly · 发表于 Post on 2021-9-2 13:21:24

biogon 发表于 2021-9-2 08:28
这个我看只能算单重态就还没用过，这个和dlpno-steom一样挺费硬盘的吧

花两万块钱买个组装机，再扩一万块钱的内存（给orca的ice-ci预备的，eom用不了这么多），只能一个任务独占，可以随便糟蹋

biogon · 发表于 Post on 2021-9-2 08:28:31

beefly 发表于 2021-9-1 19:22
orc可以做eom-dlpno-ccsd。我用def2-tzvpp(-f)基组配200gb内存算50-60原子没问题，不过目前能算单态

这个我看只能算单重态就还没用过，这个和dlpno-steom一样挺费硬盘的吧

beefly · 发表于 Post on 2021-9-1 19:22:39

biogon 发表于 2021-9-1 16:41
eom就是太昂贵了

orc可以做eom-dlpno-ccsd。我用def2-tzvpp(-f)基组配200gb内存算50-60原子没问题，不过目前能算单态

biogon · 发表于 Post on 2021-9-1 16:41:42

beefly 发表于 2021-9-1 15:38
对，但是在eom计算中不会涉及，所以更方便使用

eom就是太昂贵了

beefly · 发表于 Post on 2021-9-1 15:38:40

biogon 发表于 2021-9-1 14:53
中性分子单参考但是离子是多参考的情况是存在的吧

对，但是在eom计算中不会涉及，所以更方便使用

biogon · 发表于 Post on 2021-9-1 14:53:42

beefly 发表于 2021-9-1 14:19
steom假设中性分子、阴阳离子都是单参考，因此比一般的eom受限制更多

中性分子单参考但是离子是多参考的情况是存在的吧

beefly · 发表于 Post on 2021-9-1 14:19:30

biogon 发表于 2021-8-29 20:49
STEOM在算的时候为什么容易在IP-EOM或者EA-EOM的时候不收敛

steom假设中性分子、阴阳离子都是单参考，因此比一般的eom受限制更多

beefly · 发表于 Post on 2021-9-1 14:17:44

暖空发表于 2021-8-29 18:36
感谢回复！我现在用的服务器的内存和磁盘容量差了点儿，在用ADC2测试60个原子左右的分子时磁盘容量已经爆 ...

没有可比性。你的adc2没有用到密度拟合，只能算小分子

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] 如何在ORCA中用ADC2方法计算S1和T1的激发能？

浏览过的版块