应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Turbomole/ORCA » 怎么基于STEOM-DLPNO-CCSD计算rISC
151/1返回列表
查看: 4066  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hcmzj668 的 2 个金币 ，回帖就立即获得 1 个金币，每人有 1 次机会

hcmzj668

金虫 (初入文坛)

[交流] 怎么基于STEOM-DLPNO-CCSD计算rISC

各位同行们，正常的td-dft计算risc (或isc)已按照orca手册成功，现在想计算基于steom-dlpno-ccsd的risc，不知道步骤和输入文件怎么编写，遂来求助。

手册中写道："theinput is similar to those discussed above. here esd(isc) should be used on themain input to indicate an intersystem crossing calculation, and the hessians should be provided byiscishessian and iscfshessian for the initial and final states, respectively. please note that the geometryused in the input file should correspond to that of the final state, specified through the iscfshessian flag. the relevantmatrix elements can be computed using any method available in orca.as a simple example, one could compute theexcited singlet and ground triplet geometries and hessians for anthracene usingtd-dft. then, compute the spin-orbit coupling (soc) matrix elements for aspecific triplet spin-sublevel using the same method (see the details below),potentially employing methods like casscf, mrci,steom-ccsd"

即：先得到gs，es和三重态的结构和hessians，然后基于gs结构做steom-dlpno-ccsd下的单线态和三重态的能量以用于dele。
但是，手册中没有给出具体的输入例子，就有点无从下手了，于是产生了如下问题和想法，麻烦各位老师帮忙建议。

1. “gs，es和三重态的结构和hessians”可以用td-dft方法得到，后直接使用？

2. 基于上述得到的gs结构，做steom-dlpno-ccsd计算得到s1和t1的能量，是如下结果不呢？
steom-ccsdresults (singlets)
iroot= 1: 0.109392 au 2.977 ev 24008.9 cm**-1
steom-ccsd results (triplets)
iroot= 1: 0.089216 au 2.428 ev 19580.7 cm**-1
所以dele for risc be : 19580.7 - 24008.9 = -4,428.2 cm-1?


3. 最后，最重要的步骤计算risc:
a)是用手册中提供的例子作为输入文件（td-dft），将上述td-dft得到的s1-t1结构和hessian， steom-dlpno-ccsd得到的能量，写入如下输入文件中来计算？
%tddft
nroots 5
sroot 1
troot 1
trootssl 0 1 -1 #t1的三个亚能级
dosoc
%esd
iscishess "t1.hess"
iscfshess "s1.hess"
dele xxx # 用steom-dlpno-ccsd计算得到的
（应该是3个job，此处省略了）

b)还是直接用td-dft得到的s1-t1结构和hessian，按照steom-dlpno-ccsd的输入文件计算？（这里不会写，尝试着写了如下输入文件，不知道对还是错）
!steom-dlpno-ccsddef2-tzvp def2-tzvp/c def2/j tightscf tightpno esd(isc)
%maxcore 5000
%scf
maxiter 7000
end
%mdci
nroots 3 #如何控制t1的3个亚能级呢？（直接按照手册中的td-dft方式写3个job？）
dorootwise true
othresh 0.001
vthresh 0.001
tcutpnosingles 1e-11
ndav 400
dostoresteom true
dosimpledens false
addl2term true
dtol 1e-5
dotriplet true
steomsoc true
end
%esd
iscishess "t1-opt.hess"
iscfshess "s1-opt.hess"
approxaden true # 这个 keyword会直接去自动获取dele不呢？
usej true
doht false
temp 298.15
%pal nprocs 24
end
* xyzfile s0-opt.xyz 0 1
最后，麻烦各位行家了，祝大家生活愉快，科研顺利。
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2024-11-25 12:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

paramecium86

版主 (著名写手)
引用回帖:
12楼: Originally posted by hcmzj668 at 2024-11-25 22:57:07
版主 这是激发能，不是绝热能，对吧？...

哦 你的重点是问绝热还是垂直的能量差。我理解得用绝热的。也就是你分别优化出s1和t1结构, 分别算steom-cc得到的s1 和 t1能量。然后再作差。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

hcmzj668
13楼2024-11-26 19:33:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

paramecium86

版主 (著名写手)
★
hcmzj668(金币+1): 谢谢参与
ORCA用ECD方式算RISC。这个计算需要方法支持解析梯度（Hessian）那基本就限定住了TD方式。然后如你第一种方式所说那就自己算了STEOM下的DELE和SOC 填到输入文件里就可以了。近似当做是在STEOM-DLPNO-CC方式下的RISC了。
赞一下 回复此楼
2楼2024-11-25 13:53:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hcmzj668

金虫 (初入文坛)
引用回帖:
2楼: originally posted by paramecium86 at 2024-11-25 13:53:13
orca用ecd方式算risc。这个计算需要方法支持解析梯度（hessian）那基本就限定住了td方式。然后如你第一种方式所说那就自己算了steom下的dele和soc 填到输入文件里就可以了。近似当做是在steom-dlpno-cc方式下的risc ...

好嘞，谢谢版主大人的回复。

版主大人，关于能量，这里找得对吗？
基于上述得到的gs结构，做steom-dlpno-ccsd计算得到s1和t1的能量，是如下结果不呢？

steom-ccsdresults (singlets)
iroot= 1: 0.109392 au 2.977 ev 24008.9 cm**-1
steom-ccsd results (triplets)
iroot= 1: 0.089216 au 2.428 ev 19580.7 cm**-1

所以DELE for risc be : 19580.7 - 24008.9 = -4,428.2 cm-1?
赞一下 回复此楼
3楼2024-11-25 14:27:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
引用回帖:
3楼: Originally posted by hcmzj668 at 2024-11-25 01:27:07
好嘞，谢谢版主大人的回复。

版主大人，关于能量，这里找得对吗？
基于上述得到的gs结构，做steom-dlpno-ccsd计算得到s1和t1的能量，是如下结果不呢？

steom-ccsdresults (singlets)
iroot= 1: 0.109 ...

这个没啥问题的。
回复此楼
11楼2024-11-26 11:53:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hcmzj668

金虫 (初入文坛)
引用回帖:
11楼: Originally posted by paramecium86 at 2024-11-26 11:53:04
这个没啥问题的。...

版主 这是激发能，不是绝热能，对吧？
赞一下 回复此楼
12楼2024-11-26 11:57:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hcmzj668

金虫 (初入文坛)
送红花一朵
引用回帖:
13楼: Originally posted by paramecium86 at 2024-11-26 19:33:33
哦 你的重点是问绝热还是垂直的能量差。我理解得用绝热的。也就是你分别优化出s1和t1结构, 分别算steom-cc得到的s1 和 t1能量。然后再作差。...

好嘞 谢谢您

[ 发自手机版 https://muchong.com/3g ]
回复此楼
14楼2024-11-26 22:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
nono20094楼
2024-11-25 15:43   回复  
hcmzj668(金币+1): 谢谢参与
·
miaojiabing5楼
2024-11-25 16:39   回复  
hcmzj668(金币+1): 谢谢参与
庚辰辛丑6楼
2024-11-25 21:22   回复  
hcmzj668(金币+1): 谢谢参与
psylhh7楼
2024-11-25 22:07   回复  
hcmzj668(金币+1): 谢谢参与
tzynew8楼
2024-11-25 23:10   回复  
hcmzj668(金币+1): 谢谢参与
lymjjh9楼
2024-11-26 08:34   回复  
hcmzj668(金币+1): 谢谢参与
lymjjh10楼
2024-11-26 08:35   回复  
月夜风波15楼
2024-11-29 15:13   回复  
hcmzj668(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 hcmzj668 的主题更新
151/1返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
信息提示
关闭
请填处理意见
关闭