【有獎交流】積極回復本帖子，參與交流，就有機會分得作者 hcmzj668 的 2 個金幣，回帖就立即獲得 1 個金幣，每人有 1 次機會

hcmzj668

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 737.3
帖子: 38
在線: 101.6小時
蟲號: 2565211

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

各位同行們，正常的td-dft計算risc (或isc)已按照orca手冊成功，現(xiàn)在想計算基于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和三重態(tài)的結(jié)構(gòu)和hessians，然后基于gs結(jié)構(gòu)做steom-dlpno-ccsd下的單線態(tài)和三重態(tài)的能量以用于dele。
但是，手冊中沒有給出具體的輸入例子，就有點無從下手了，于是產(chǎn)生了如下問題和想法，麻煩各位老師幫忙建議。

1. “gs，es和三重態(tài)的結(jié)構(gòu)和hessians”可以用td-dft方法得到，后直接使用？

2. 基于上述得到的gs結(jié)構(gòu)，做steom-dlpno-ccsd計算得到s1和t1的能量，是如下結(jié)果不呢？
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結(jié)構(gòu)和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結(jié)構(gòu)和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
最后，麻煩各位行家了，祝大家生活愉快，科研順利。

回復此樓

» 猜你喜歡

PbS量子點紫外吸收已經(jīng)有0人回復
比利時魯汶大學與國家留學基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship 已經(jīng)有0人回復
物理化學論文潤色/翻譯怎么收費? 已經(jīng)有265人回復
請問四氫呋喃溶解的聚合物用甲醇沉淀時，如何使沉淀過程加速？已經(jīng)有2人回復
七嗪類物質(zhì)合成求助已經(jīng)有0人回復
倫敦大學學院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生已經(jīng)有11人回復
Chemical Bonding at Surfaces and Interfaces,最經(jīng)典的一本表面上化學相互作用教材已經(jīng)有0人回復
中國科學院成都山地災害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復
誠邀有志之士加入江蘇大學環(huán)境學院合成生物學課題組！(長期接收學碩專碩、碩博連讀) 已經(jīng)有16人回復
一志愿南大理論與計算化學方向（070300）342分，求調(diào)劑已經(jīng)有0人回復
歡迎化學、化工、物理、材料、或計算機背景的考研咨詢已經(jīng)有0人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2024-11-25 12:51:30

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

paramecium86

版主 (著名寫手)

引用回帖:

12樓: Originally posted by hcmzj668 at 2024-11-25 22:57:07
版主這是激發(fā)能，不是絕熱能，對吧？...

哦你的重點是問絕熱還是垂直的能量差。我理解得用絕熱的。也就是你分別優(yōu)化出s1和t1結(jié)構(gòu), 分別算steom-cc得到的s1 和 t1能量。然后再作差。

贊一下(1人)

回復此樓

» 本帖已獲得的紅花（最新10朵）

hcmzj668

13樓2024-11-26 19:33:33

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

paramecium86

版主 (著名寫手)

★
hcmzj668(金幣+1): 謝謝參與

ORCA用ECD方式算RISC。這個計算需要方法支持解析梯度（Hessian）那基本就限定住了TD方式。然后如你第一種方式所說那就自己算了STEOM下的DELE和SOC 填到輸入文件里就可以了。近似當做是在STEOM-DLPNO-CC方式下的RISC了。

贊一下

回復此樓

2樓2024-11-25 13:53:13

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hcmzj668

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 737.3
帖子: 38
在線: 101.6小時
蟲號: 2565211

引用回帖:

2樓: originally posted by paramecium86 at 2024-11-25 13:53:13
orca用ecd方式算risc。這個計算需要方法支持解析梯度（hessian）那基本就限定住了td方式。然后如你第一種方式所說那就自己算了steom下的dele和soc 填到輸入文件里就可以了。近似當做是在steom-dlpno-cc方式下的risc ...

好嘞，謝謝版主大人的回復。

版主大人，關(guān)于能量，這里找得對嗎？
基于上述得到的gs結(jié)構(gòu)，做steom-dlpno-ccsd計算得到s1和t1的能量，是如下結(jié)果不呢？

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

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

paramecium86

版主 (著名寫手)

引用回帖:

3樓: Originally posted by hcmzj668 at 2024-11-25 01:27:07
好嘞，謝謝版主大人的回復。

版主大人，關(guān)于能量，這里找得對嗎？
基于上述得到的gs結(jié)構(gòu)，做steom-dlpno-ccsd計算得到s1和t1的能量，是如下結(jié)果不呢？

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

這個沒啥問題的。

回復此樓

11樓2024-11-26 11:53:04

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hcmzj668

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 737.3
帖子: 38
在線: 101.6小時
蟲號: 2565211

引用回帖:

11樓: Originally posted by paramecium86 at 2024-11-26 11:53:04
這個沒啥問題的。...

版主這是激發(fā)能，不是絕熱能，對吧？

贊一下

回復此樓

12樓2024-11-26 11:57:07

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hcmzj668

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 737.3
帖子: 38
在線: 101.6小時
蟲號: 2565211

送紅花一朵

引用回帖:

13樓: Originally posted by paramecium86 at 2024-11-26 19:33:33
哦你的重點是問絕熱還是垂直的能量差。我理解得用絕熱的。也就是你分別優(yōu)化出s1和t1結(jié)構(gòu), 分別算steom-cc得到的s1 和 t1能量。然后再作差。...

好嘞謝謝您

[ 發(fā)自手機版 http://m.gaoyang168.com/3g ]

回復此樓

14樓2024-11-26 22:40:43

已閱回復此樓關(guān)注TA 給TA發(fā)消息送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): 謝謝參與

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 274求調(diào)劑0856材料化工 +9	z2839474511 2026-03-11	10/500	2026-03-11 23:46 by BruceLiu320
[考研] 0703求調(diào)劑 +6	jtyq001 2026-03-10	6/300	2026-03-11 22:18 by peike
[考研] 293求調(diào)劑，一志愿陜師大生物學 +3	??????.?.??? 2026-03-09	3/150	2026-03-11 10:02 by 學員8dgXkO
[考研] 279求調(diào)劑 +3	莫xiao 2026-03-10	4/200	2026-03-11 08:06 by 斬魂滴兔子！
[考研] 化工0817調(diào)劑 +8	燦若星晨 2026-03-10	8/400	2026-03-10 22:44 by 星空星月
[考研] 085600 材料與化工 295 求調(diào)劑 +9	dream…… 2026-03-10	9/450	2026-03-10 22:24 by 劍詩杜康
[考研] 調(diào)劑 +5	調(diào)劑的考研學生 2026-03-09	5/250	2026-03-10 22:13 by peike
[考研] 086000生物與醫(yī)藥319分求調(diào)劑 +4	Tolkien 2026-03-07	8/400	2026-03-10 21:34 by Tolkien
[考研] 0856材料與化工309分求調(diào)劑 +4	ZyZy…… 2026-03-10	4/200	2026-03-10 15:20 by houyaoxu
[基金申請] PDF在線壓縮軟件 10+3	haxia 2026-03-08	4/200	2026-03-10 13:44 by lfy8008
[考研] 復試調(diào)劑 +6	呼呼？~+123456 2026-03-08	8/400	2026-03-10 11:36 by 呼呼？~+123456
[基金申請] 面上項目還需要AI說明嗎？ +3	liyundong 2026-03-08	3/150	2026-03-09 22:30 by kingkocxr
[考研] 求調(diào)劑，數(shù)一英一274分 +4	小菲會努力 2026-03-08	4/200	2026-03-09 12:40 by 一定上岸喲_
[考研] 293一志愿華東理工 0817化學工程與技術(shù) 調(diào)劑 +5	fjj0912 2026-03-07	5/250	2026-03-09 09:13 by 30660438
[考博] 26申博自薦 +9	烷基八氮- 2026-03-06	11/550	2026-03-08 23:20 by gaobiao
[考研] 考研求調(diào)劑，工科，總分304 +4	矮子不想長高了 2026-03-08	7/350	2026-03-08 22:21 by 小木蟲tim
[考研] 288求調(diào)劑（0703）一志愿東北大學 +5	好好- 2026-03-07	5/250	2026-03-08 18:22 by zzpnuli111
[考研] 求調(diào)劑，一志愿華中科大0702，數(shù)一英一，293 +4	小羅露一二 2026-03-07	4/200	2026-03-08 16:36 by 星空星月
[考研] 070300化學求調(diào)劑292分 +3	打烊eee 2026-03-07	3/150	2026-03-07 20:04 by houyaoxu
[考研] 復試調(diào)劑 +7	呼呼？~+123456 2026-03-05	10/500	2026-03-06 16:07 by @颯颯颯颯