 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2222 | 回復(fù): 14 | |||
| 【有獎(jiǎng)交流】積極回復(fù)本帖子，參與交流，就有機(jī)會(huì)分得作者 hcmzj668 的 2 個(gè)金幣 ，回帖就立即獲得 1 個(gè)金幣，每人有 1 次機(jī)會(huì) | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
怎么基于STEOM-DLPNO-CCSD計(jì)算rISC
|
|||
|
各位同行們，正常的td-dft計(jì)算risc (或isc)已按照orca手冊(cè)成功，現(xiàn)在想計(jì)算基于steom-dlpno-ccsd的risc，不知道步驟和輸入文件怎么編寫，遂來求助。 手冊(cè)中寫道："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è)中沒有給出具體的輸入例子，就有點(diǎn)無從下手了，于是產(chǎ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得到的能量，寫入如下輸入文件中來計(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ì)寫，嘗試著寫了如下輸入文件，不知道對(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方式寫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ù)此樓» 猜你喜歡
陳教授團(tuán)隊(duì)招收2026級(jí)學(xué)術(shù)/專業(yè)型碩士研究生（調(diào)劑）
已經(jīng)有0人回復(fù)
-
招聘研究助理
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有113人回復(fù)
-
國(guó)家級(jí)人才課題組招收2026年入學(xué)博士
已經(jīng)有0人回復(fù)
-
申博求助
已經(jīng)有3人回復(fù)
-
請(qǐng)問紅外光譜除了賽默飛以外國(guó)產(chǎn)哪家最強(qiáng)？
已經(jīng)有1人回復(fù)
-
國(guó)家級(jí)人才課題組招收2026年入學(xué)博士
已經(jīng)有0人回復(fù)
-
304求調(diào)劑
已經(jīng)有3人回復(fù)
-
國(guó)家級(jí)人才課題組招收2026年入學(xué)博士
已經(jīng)有1人回復(fù)
-
304求調(diào)劑
已經(jīng)有4人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
有機(jī)合成崗位招聘
+1/939
-
[求助]日本動(dòng)漫里用手表和三角板測(cè)經(jīng)緯度是怎么回事？
+1/188
-
海南大學(xué)徐月山老師招生第二批博士名額2~3個(gè)，2026年9月份入學(xué)（高端設(shè)備開發(fā)方向）
+1/95
-
江西南部某一本高校招收環(huán)境、生物、計(jì)算機(jī)、農(nóng)業(yè)等專業(yè)調(diào)劑生
+1/82
-
北京信息科技大學(xué)智能機(jī)器人技術(shù)方向，聯(lián)合清華培養(yǎng)碩生、博生，并長(zhǎng)期招博士后
+1/80
-
課題組宣傳
+1/35
-
【博士/研究型碩士/訪問學(xué)者】 澳大利亞國(guó)立大學(xué) (ANU) | 電化學(xué)儲(chǔ)能材料方向
+1/34
-
山東大學(xué)（青島校區(qū)）招副研究員及博士后（COF\MOF\催化\電池）
+1/17
-
河南大學(xué)張琳教授團(tuán)隊(duì)招收博士研究生
+1/15
-
綠色催化氧化團(tuán)隊(duì)接收調(diào)劑研究生
+1/15
-
【碩士調(diào)劑招生】武漢科技大學(xué)招收材料、化工、環(huán)境、冶金類2026級(jí)碩士生
+1/8
-
招收土木工程學(xué)科擬調(diào)劑碩士生1名，急！
+3/7
-
教育部長(zhǎng)江學(xué)者和創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃”入選團(tuán)隊(duì)招收 材料，化學(xué)與化工博士研究生
+1/6
-
306求調(diào)劑
+1/6
-
【博士招生】北方工業(yè)大學(xué)招收 機(jī)械工程方向 申請(qǐng)考核制博士
+1/6
-
-分子探針等化學(xué)/生物方向博士后和研究助理
+1/4
-
關(guān)于文章OE投稿大修求助
+1/3
-
大灣區(qū)大學(xué)姚瑤課題組2026屆博士生招生
+1/3
-
中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所博士后招聘啟事
+1/3
-
2026年浙江清華長(zhǎng)三角研究院智慧能源中心招聘博士后-電氣方向
+1/1
paramecium86
版主 (著名寫手)
- QC強(qiáng)帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26376.6
- 帖子: 2025
- 在線: 1441.9小時(shí)
- 蟲號(hào): 2466443
11樓2024-11-26 11:53:04
paramecium86
版主 (著名寫手)
- QC強(qiáng)帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26376.6
- 帖子: 2025
- 在線: 1441.9小時(shí)
- 蟲號(hào): 2466443
2樓2024-11-25 13:53:13
|
好嘞，謝謝版主大人的回復(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? |
3樓2024-11-25 14:27:07
12樓2024-11-26 11:57:07
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩調(diào)劑 +5 | 慕辰123 2026-03-05 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
|
[考博] 26申博自薦 +7 | 烷基八氮- 2026-03-06 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +4 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +5 | forgman95 2026-03-05 | 7/350 |
|
|
|
[考研] 一志愿211，0860總分286食品工程專業(yè)求調(diào)劑 +6 | 林林Winnie 2026-03-05 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +5 | wxz2 2026-03-03 | 8/400 |
|
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|
|
|
[考研] 材料類考研調(diào)劑 +7 | gemmgemm 2026-03-01 | 8/400 |
|
|
|
[考研] 一志愿清華深研院材料專碩294分，專業(yè)課111分，本科中南大學(xué)材料，有六級(jí)，有工作經(jīng)驗(yàn) +3 | H14528 2026-03-04 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|