| 5 | 1/1 | 返回列表 |
| 查看: 2225 | 回復(fù): 14 | |||
| 【有獎交流】積極回復(fù)本帖子,參與交流,就有機會分得作者 hcmzj668 的 2 個金幣 ,回帖就立即獲得 1 個金幣,每人有 1 次機會 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
怎么基于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計算得到的 (應(yīng)該是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 最后,麻煩各位行家了,祝大家生活愉快,科研順利。 |
» 搶金幣啦!回帖就可以得到:
+1/90
+1/81
+1/38
+1/33
+1/29
+1/26
+1/21
+1/19
+1/17
+1/15
+1/11
+1/8
+1/7
+1/5
+1/5
+1/5
+1/4
+1/3
+1/1
+1/1
版主 (著名寫手)
版主 (著名寫手)
|
好嘞,謝謝版主大人的回復(fù)。 版主大人,關(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? |
版主 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] A區(qū)一本交叉課題組,低分調(diào)劑,招收機械電子信息通信等交叉方向 +39 | lisimayy 2026-03-04 | 52/2600 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑085600材料與化工 +9 | Daunrin 2026-03-07 | 10/500 |
|
|
[考研] 081700學(xué)碩一志愿北京化工大學(xué)數(shù)二英一過六級有競賽求調(diào)劑 +4 | galaxary 2026-03-07 | 6/300 |
|
|
[考博] 申博 +3 | 添菜了哈 2026-03-04 | 6/300 |
|
|
[考研]
|
Sixuan wang 2026-03-06 | 5/250 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
[考研] 085602高分子方向求調(diào)劑 +7 | tlgudy 2026-03-04 | 7/350 |
|
|
[考博] 26申博-目前4篇SCI一作 +4 | chen_2024 2026-03-02 | 4/200 |
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
[考研] 283求調(diào)劑 +10 | 鹿沫笙 2026-03-02 | 11/550 |
|
|
[考研] 281電子信息求調(diào)劑 +5 | jhtfeybgj 2026-03-02 | 9/450 |
|
|
[考研] 0856材料專碩274能調(diào)劑去哪里? +3 | 22735 2026-03-04 | 4/200 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 0856材料工程,初試313調(diào)劑 +7 | 賣個關(guān)子吧 2026-03-03 | 7/350 |
|
|
[基金申請] 請問大家,研究風(fēng)險與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|