【有獎(jiǎng)交流】積極回復(fù)本帖子，參與交流，就有機(jī)會(huì)分得作者 hcmzj668 的 2 個(gè)金幣，回帖就立即獲得 1 個(gè)金幣，每人有 1 次機(jī)會(huì)

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

hcmzj668

應(yīng)助: 0 (幼兒園)
金幣: 737.3
帖子: 38
在線: 101.6小時(shí)
蟲(chóng)號(hào): 2565211

[交流] 怎么基于STEOM-DLPNO-CCSD計(jì)算rISC

各位同行們，正常的td-dft計(jì)算risc (或isc)已按照orca手冊(cè)成功，現(xiàn)在想計(jì)算基于steom-dlpno-ccsd的risc，不知道步驟和輸入文件怎么編寫(xiě)，遂來(lái)求助。

手冊(cè)中寫(xiě)道："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。
但是，手冊(cè)中沒(méi)有給出具體的輸入例子，就有點(diǎn)無(wú)從下手了，于是產(chǎn)生了如下問(wèn)題和想法，麻煩各位老師幫忙建議。

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

2. 基于上述得到的gs結(jié)構(gòu)，做steom-dlpno-ccsd計(jì)算得到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. 最后，最重要的步驟計(jì)算risc:
a)是用手冊(cè)中提供的例子作為輸入文件（td-dft），將上述td-dft得到的s1-t1結(jié)構(gòu)和hessian， steom-dlpno-ccsd得到的能量，寫(xiě)入如下輸入文件中來(lái)計(jì)算？
%tddft
nroots 5
sroot 1
troot 1
trootssl 0 1 -1 #t1的三個(gè)亞能級(jí)
dosoc
%esd
iscishess "t1.hess"
iscfshess "s1.hess"
dele xxx # 用steom-dlpno-ccsd計(jì)算得到的
（應(yīng)該是3個(gè)job，此處省略了）

b)還是直接用td-dft得到的s1-t1結(jié)構(gòu)和hessian，按照steom-dlpno-ccsd的輸入文件計(jì)算？（這里不會(huì)寫(xiě)，嘗試著寫(xiě)了如下輸入文件，不知道對(duì)還是錯(cuò)）
!steom-dlpno-ccsddef2-tzvp def2-tzvp/c def2/j tightscf tightpno esd(isc)
%maxcore 5000
%scf
maxiter 7000
end
%mdci
nroots 3 #如何控制t1的3個(gè)亞能級(jí)呢？（直接按照手冊(cè)中的td-dft方式寫(xiě)3個(gè)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 # 這個(gè) keyword會(huì)直接去自動(dòng)獲取dele不呢？
usej true
doht false
temp 298.15
%pal nprocs 24
end
* xyzfile s0-opt.xyz 0 1
最后，麻煩各位行家了，祝大家生活愉快，科研順利。

回復(fù)此樓

» 猜你喜歡

光電子賦能·半導(dǎo)體創(chuàng)芯！已經(jīng)有1人回復(fù)
光電子賦能·半導(dǎo)體創(chuàng)芯！已經(jīng)有0人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有136人回復(fù)
QE計(jì)算電聲耦合的時(shí)候報(bào)錯(cuò)Error in routine lambda (100)wrong or too many modes 已經(jīng)有2人回復(fù)
中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復(fù)
cif求助下載已經(jīng)有1人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有5人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有0人回復(fù)
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有2人回復(fù)
化學(xué)化工學(xué)院招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）已經(jīng)有0人回復(fù)
348求調(diào)劑已經(jīng)有4人回復(fù)

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

查看全部散金貼

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

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

hcmzj668

金蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 737.3
帖子: 38
在線: 101.6小時(shí)
蟲(chóng)號(hào): 2565211

引用回帖:

11樓: Originally posted by paramecium86 at 2024-11-26 11:53:04
這個(gè)沒(méi)啥問(wèn)題的。...

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

贊一下

回復(fù)此樓

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

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

查看全部 15 個(gè)回答

paramecium86

版主 (著名寫(xiě)手)

QC強(qiáng)帖: 5
應(yīng)助: 587 (博士)
貴賓: 0.01
金幣: 26376.6
帖子: 2025
在線: 1441.9小時(shí)
蟲(chóng)號(hào): 2466443

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

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

贊一下

回復(fù)此樓

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

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

hcmzj668

金蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 737.3
帖子: 38
在線: 101.6小時(shí)
蟲(chóng)號(hào): 2565211

引用回帖:

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

好嘞，謝謝版主大人的回復(fù)。

版主大人，關(guān)于能量，這里找得對(duì)嗎？
基于上述得到的gs結(jié)構(gòu)，做steom-dlpno-ccsd計(jì)算得到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?

贊一下

回復(fù)此樓

3樓2024-11-25 14:27:07

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

paramecium86

版主 (著名寫(xiě)手)

QC強(qiáng)帖: 5
應(yīng)助: 587 (博士)
貴賓: 0.01
金幣: 26376.6
帖子: 2025
在線: 1441.9小時(shí)
蟲(chóng)號(hào): 2466443

引用回帖:

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

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

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

這個(gè)沒(méi)啥問(wèn)題的。

回復(fù)此樓

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

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

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

[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國(guó)金及主持省級(jí)大創(chuàng)經(jīng)歷 +3	熊境喆 2026-03-05	3/150	2026-03-07 14:32 by lishahe
[考研] 材料調(diào)劑 +10	ounce. 2026-03-04	12/600	2026-03-07 09:14 by Mornach1988
[考研] 材料調(diào)劑 +5	ounce. 2026-03-05	13/650	2026-03-06 20:34 by Eember.
[考研] 一志愿211，0860總分286食品工程專業(yè)求調(diào)劑 +6	林林Winnie 2026-03-05	6/300	2026-03-06 20:06 by 我想上研究生.
[考研] 332材料求調(diào)劑 +6	zjy101327 2026-03-05	7/350	2026-03-06 00:18 by wutongshun
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4	啊歐歐歐 2026-03-04	5/250	2026-03-05 19:27 by Leding1356
[考研] 材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/ 5+3	蘇堯幺幺 2026-03-03	8/400	2026-03-05 16:06 by njzyff
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +5	Seele_v5 2026-03-04	5/250	2026-03-05 13:26 by Iveryant
[考研] 中國(guó)科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑，有科研和獲獎(jiǎng)經(jīng)歷 +11	wsxw 2026-03-02	12/600	2026-03-05 09:49 by oxidpl
[考研] 化學(xué) 0703求調(diào)劑總分293 一志愿211 +9	土土小蟲(chóng) 2026-03-03	11/550	2026-03-05 09:26 by vfdbghn
[考研] 0856材料與化工，270求調(diào)劑 +17	YXCT 2026-03-01	20/1000	2026-03-05 09:22 by 一切OK
[考研] 0856材料專碩274能調(diào)劑去哪里？ +3	22735 2026-03-04	4/200	2026-03-05 09:06 by 斬魂滴兔子！
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3	Hz795795 2026-03-04	3/150	2026-03-04 20:54 by wutongshun
[考研] 一志愿中南大學(xué)理學(xué)化學(xué) +6	15779376950 2026-03-01	8/400	2026-03-04 16:19 by zhukairuo
[考研] 0703化學(xué)求調(diào)劑 +3	將未有期 2026-03-04	3/150	2026-03-04 16:19 by sslc1985
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3	zju51 2026-03-04	3/150	2026-03-04 13:01 by zhukairuo
[考研] 298求調(diào)劑一志愿中海洋 +3	lour. 2026-03-03	3/150	2026-03-03 20:41 by gxg2025
[考研] 材料270求調(diào)劑 6+6	Eiiiio 2026-03-01	11/550	2026-03-03 15:04 by Phoebelala
[考研] 295求調(diào)劑。一志愿報(bào)考鄭州大學(xué)化學(xué)工藝學(xué)碩，總分295分 +8	yl1 2026-03-02	9/450	2026-03-02 21:21 by sunny81
[考博] 誠(chéng)招農(nóng)業(yè)博士 +3	心欣向榮 2026-02-28	3/150	2026-03-02 13:33 by 時(shí)間不狗