 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 15 | 1/1 | 返回列表 |
| 查看: 2223 | 回復(fù): 14 | |||
| 【有獎(jiǎng)交流】積極回復(fù)本帖子，參與交流，就有機(jī)會(huì)分得作者 hcmzj668 的 2 個(gè)金幣 ，回帖就立即獲得 1 個(gè)金幣，每人有 1 次機(jī)會(huì) | |||
[交流]
怎么基于STEOM-DLPNO-CCSD計(jì)算rISC
|
|||
|
各位同行們，正常的td-dft計(jì)算risc (或isc)已按照orca手冊(cè)成功，現(xiàn)在想計(jì)算基于steom-dlpno-ccsd的risc，不知道步驟和輸入文件怎么編寫，遂來(lái)求助。 手冊(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è)中沒(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得到的能量，寫入如下輸入文件中來(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ì)寫，嘗試著寫了如下輸入文件，不知道對(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ù)此樓» 猜你喜歡
硅碳半電池循環(huán)異常失效
已經(jīng)有1人回復(fù)
-
26年儲(chǔ)能、電池方向博士申請(qǐng)
已經(jīng)有1人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有121人回復(fù)
-
吉林大學(xué)物理化學(xué)專業(yè)招生考核制博士
已經(jīng)有20人回復(fù)
-
吉林大學(xué)物理化學(xué)專業(yè)招生2026級(jí)考核制博士
已經(jīng)有20人回復(fù)
-
一志愿湖南大學(xué)0703，356分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
0703化學(xué)求調(diào)劑學(xué)校
已經(jīng)有2人回復(fù)
-
070300（化學(xué)）一志愿東師 總分320分求調(diào)劑
已經(jīng)有1人回復(fù)
-
澳大利亞皇家墨爾本理工大學(xué)（RMIT）理學(xué)院 院士團(tuán)隊(duì)2026招收CSC博士生
已經(jīng)有11人回復(fù)
-
大連工業(yè)大學(xué)杰青/長(zhǎng)江團(tuán)隊(duì)招收儲(chǔ)能電池方向博士
已經(jīng)有0人回復(fù)
-
267求調(diào)劑，一志愿廈大0856，有競(jìng)賽有科研經(jīng)歷
已經(jīng)有2人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
清華-西電集成電路方向招收聯(lián)培博士（普博）
+2/218
-
2026級(jí)博士招生｜蘇州國(guó)家實(shí)驗(yàn)室 —上海交通大學(xué)聯(lián)合培養(yǎng)（授上交學(xué)位）
+2/86
-
華中科技大學(xué)物理學(xué)院于濤教授課題組 2026年擬招收優(yōu)秀考研生1名
+2/86
-
arXiv 預(yù)印本求推薦
+2/84
-
北京信息科技大學(xué)聯(lián)合清華大學(xué)某國(guó)家重點(diǎn)實(shí)驗(yàn)室，聯(lián)合培養(yǎng)碩生、博生，并長(zhǎng)期招博士后
+1/83
-
長(zhǎng)江大學(xué)武漢校區(qū)新能源博士招聘-2026
+1/41
-
招博士，申請(qǐng)考核制，環(huán)境化學(xué)、理論計(jì)算、環(huán)境毒理/健康方向
+1/25
-
【已截止】
+1/19
-
青年長(zhǎng)江學(xué)者團(tuán)隊(duì) 2026年 碩士研究生招生（包括第一志愿和調(diào)劑考生）
+1/15
-
上海中醫(yī)藥大學(xué) 謝老師課題組招收2026年申請(qǐng)考核博士生1名
+1/13
-
【碳足跡與碳中和戰(zhàn)略】招收碩士、博士、博士后、科研助理
+1/11
-
中南林業(yè)科技大學(xué)《微生物學(xué)》2026年博士生招生（“申請(qǐng)-審核”錄取）
+1/8
-
溫州大學(xué)化材學(xué)院王娟課題組招生
+1/8
-
306求調(diào)劑
+1/6
-
電子科技大學(xué)巫江團(tuán)隊(duì)招收博士研究生
+1/4
-
關(guān)于文章OE投稿大修求助
+1/3
-
北京高校副校長(zhǎng)團(tuán)隊(duì)招收機(jī)械類，環(huán)境類學(xué)碩和專碩
+1/2
-
青島大學(xué)-馬麗春教授-材料學(xué)院2026年招收博士研究生
+1/2
-
清華大學(xué)深圳國(guó)際研究生院張旻課題組招收2027級(jí)博士生_機(jī)械、儀器方向
+1/1
-
青島科技大學(xué)高分子學(xué)院老師招收碩士生2名
+1/1
paramecium86
版主 (著名寫手)
- QC強(qiáng)帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26376.6
- 帖子: 2025
- 在線: 1441.9小時(shí)
- 蟲號(hào): 2466443
hcmzj66813樓2024-11-26 19:33:33
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
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
12樓2024-11-26 11:57:07
14樓2024-11-26 22:40:43
簡(jiǎn)單回復(fù)
nono20094樓
2024-11-25 15:43
回復(fù)
hcmzj668(金幣+1): 謝謝參與
·
2024-11-25 16:39
回復(fù)
hcmzj668(金幣+1): 謝謝參與
庚辰辛丑6樓
2024-11-25 21:22
回復(fù)
hcmzj668(金幣+1): 謝謝參與
psylhh7樓
2024-11-25 22:07
回復(fù)
hcmzj668(金幣+1): 謝謝參與
tzynew8樓
2024-11-25 23:10
回復(fù)
hcmzj668(金幣+1): 謝謝參與
lymjjh9樓
2024-11-26 08:34
回復(fù)
hcmzj668(金幣+1): 謝謝參與
lymjjh10樓
2024-11-26 08:35
回復(fù)
月夜風(fēng)波15樓
2024-11-29 15:13
回復(fù)
hcmzj668(金幣+1): 謝謝參與
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hcmzj668 的主題更新
| 15 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)（0805）338 求調(diào)劑 +3 | xiaokang3286 2026-03-07 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 第一志愿上海大學(xué)，專業(yè)化學(xué)工程與技術(shù)，總分288，求調(diào)劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
[考研] 一志愿南京大學(xué)資源與環(huán)境專業(yè)283分求調(diào)劑 +8 | 追云?? 2026-03-04 | 8/400 |
|
|
|
[考研] 一志愿211，0860總分286食品工程專業(yè)求調(diào)劑 +6 | 林林Winnie 2026-03-05 | 6/300 |
|
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +9 | llx0610 2026-03-02 | 10/500 |
|
|
|
[考研] 301求調(diào)劑 +7 | 朝天椒的雙馬尾 2026-03-05 | 7/350 |
|
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-05 | 5/250 |
|
|
|
[考研] 080500材料科學(xué)與工程 +13 | 202114020319 2026-03-03 | 13/650 |
|
|
|
[考博] 26申博-目前4篇SCI一作 +4 | chen_2024 2026-03-02 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 302材料工程求調(diào)劑 +7 | Doleres 2026-03-01 | 8/400 |
|
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
|
[考研] 0703化學(xué) 一志愿東北大學(xué) 279分求調(diào)劑 +4 | 聶樂(lè)平 2026-03-03 | 4/200 |
|
|
|
[考研]
 一志愿中南大學(xué)理學(xué)化學(xué)
+6
|
15779376950 2026-03-01 | 8/400 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|
|
|
[基金申請(qǐng)]
剛錄用，沒(méi)有期刊號(hào)，但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|