 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2462 | 回復(fù): 5 | |||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)3次，作者ghcacj增加金幣 2.4 個(gè) | |||
| 本帖產(chǎn)生 1 個(gè) QC強(qiáng)帖 ，點(diǎn)擊這里進(jìn)行查看 | |||
ghcacj榮譽(yù)版主 (著名寫手)
|
[資源]
【gchcaj個(gè)人文集】淺談Orca計(jì)算BSSE校正
|
||
|
今天看到lihb734兄寫了個(gè)關(guān)于BSSE校正的帖子，頗有啟發(fā)，所以我也來寫一個(gè)關(guān)于使用Orca校正BSSE的，大家看到的Gaussian的BSSE校正比較多，可以說，Gaussian在計(jì)算BSSE的input文件上比較方便，給出兩個(gè)分子的坐標(biāo)以及Counterpoise的關(guān)鍵詞即可，上手容易，但是細(xì)究中間的原理過程，我個(gè)人orca在這里好些，廢話不多說，先貼出BSSE的定義是，這個(gè)是Boys和Bernardi二位搞出來的，如公式1： 初看起來比較復(fù)雜，其實(shí)不難，表達(dá)式中的下表表示計(jì)算的片段，AB表示二聚體，A或B表示單聚體；上標(biāo)表示分子所處的優(yōu)化幾何坐標(biāo)，AB表示處于二聚體的優(yōu)化幾何坐標(biāo)，A或B表示單聚體優(yōu)化幾何坐標(biāo)；括號(hào)里表示使用的基組。這樣就可以把這個(gè)分成七項(xiàng)四個(gè)大項(xiàng)，第一項(xiàng)是二聚體的能量，第二項(xiàng)是兩個(gè)單聚體能量之和，第三項(xiàng)是處于二聚體優(yōu)化幾何坐標(biāo)下的使用單聚體部分的基組計(jì)算兩個(gè)單聚體的能量之和，第四項(xiàng)是處于二聚體優(yōu)化幾何坐標(biāo)下的使用二聚體部分的基組計(jì)算兩個(gè)單聚體的能量之和。所以公式2可以改寫為： 這里dimer表示二聚體，monomer2表示原公式的第四項(xiàng)，monomer1表示原公式的第三項(xiàng)，monomer0表示原公式的第二項(xiàng)。如果只計(jì)算前三項(xiàng)，就是被校正后的能量。 現(xiàn)在用一個(gè)例子來說明下，也是兩個(gè)水分子的作用，在Orca手冊(cè)91-93頁： Orca的常規(guī)語法這里不再解釋，只是介紹下輸入文件中的冒號(hào)作用，這個(gè)冒號(hào)是專門為這個(gè)BSSE計(jì)算考慮的，把冒號(hào)跟在原子后面，該原子的電子和核電荷均被刪除，但是保留原子以被分配基組。 這個(gè)任務(wù)其實(shí)是通過建立多個(gè)任務(wù)，再用關(guān)鍵詞%id將每大項(xiàng)的兩個(gè)結(jié)果能量聯(lián)系起來，部分輸出文件如下： Monomer : -152.647062118 Eh Dimer : -152.655623625 Eh -5.372 kcal/mol Monomer at dimer geometry: -152.647006948 Eh 0.035 kcal/mol Same with AB Basis set : -152.648364970 Eh -0.818 kcal/mol 按照公式2，這時(shí)候就不難計(jì)算了，ΔE=-152.655623625 -(-152.648364970+152.647006948) +152.647062118=-0.007203485 Eh 計(jì)算到了這里，應(yīng)該說Orca不如Gaussian那么智能，因?yàn)镚aussian一個(gè)關(guān)鍵詞，Orca要寫這么多，但是做計(jì)算的經(jīng)驗(yàn)告訴我，越是書寫麻煩，說明自定義的空間就越大，這里就體現(xiàn)如果我們要做一個(gè)Scan的話，Orca不用一個(gè)點(diǎn)一個(gè)點(diǎn)的計(jì)算，而且可以將BSSE校正中的每項(xiàng)拆開，分別做Scan或者SP，最后通過公式進(jìn)行處理即可。這點(diǎn)Gaussian我好像發(fā)現(xiàn)不行，因?yàn)镚aussian的Counterpoise包含了太多的東西，中間的東西無法單獨(dú)調(diào)用（本人用Gaussian較少，如果說錯(cuò)了，請(qǐng)指出，不過別罵娘就好）。這樣在計(jì)算時(shí)方便些（個(gè)人感覺） 我也來個(gè)免責(zé)聲明： 本人主業(yè)是Monte Carlo，對(duì)于QC是外行，所以這個(gè)基本是本人的一家之言，如果有問題請(qǐng)指出，不承擔(dān)任何學(xué)術(shù)上的責(zé)任。 [ Last edited by 375642546 on 2011-5-6 at 17:25 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化軟件學(xué)習(xí) |
» 猜你喜歡
CSC訪學(xué)博后項(xiàng)目獲批，外方學(xué)校暫停合作該怎么辦？
已經(jīng)有51人回復(fù)
-
求助各位大佬，球球了，在這科薩頭
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有193人回復(fù)
-
[電子教材]天津大學(xué)物理化學(xué)（第七版）上、下冊(cè)
已經(jīng)有3人回復(fù)
-
陰離子交換膜電解二氧化碳還原
已經(jīng)有0人回復(fù)
-
分子動(dòng)力學(xué)模擬合作需求
已經(jīng)有2人回復(fù)
-
九江學(xué)院2026年最新高層次人才招聘公告
已經(jīng)有0人回復(fù)
-
深圳信息職業(yè)技術(shù)學(xué)院2026年教師招聘公告（最新）
已經(jīng)有0人回復(fù)
-
研究生招生
已經(jīng)有0人回復(fù)
-
山西大學(xué)資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有2人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 很好的計(jì)算繪圖軟件教程
已經(jīng)有14人回復(fù)
- origin計(jì)算曲線面積-表面輪廓數(shù)據(jù)處理
已經(jīng)有12人回復(fù)
- 計(jì)算旋轉(zhuǎn)能壘用勢能面掃描之后還用找過渡態(tài)嗎？我是新手，希望大家給與幫助。[/font]
已經(jīng)有3人回復(fù)
- Materials studio 計(jì)算單點(diǎn)能
已經(jīng)有9人回復(fù)
- 純理論計(jì)算投稿的問題
已經(jīng)有21人回復(fù)
- |求助| 統(tǒng)計(jì)學(xué)中在進(jìn)行相關(guān)性檢驗(yàn)時(shí)，P值如何計(jì)算？
已經(jīng)有9人回復(fù)
- 有關(guān)opt+freq+bsse校正輸出文件的意思
已經(jīng)有3人回復(fù)
- BSSE校正中 片段2的變化如何改寫
已經(jīng)有8人回復(fù)
- 求助DFT計(jì)算中CuCl2在γ-Al2O3 表面的吸附結(jié)構(gòu)？
已經(jīng)有5人回復(fù)
- 【求助】用Orca計(jì)算一個(gè)片簇的單點(diǎn)能，不收斂。
已經(jīng)有7人回復(fù)
- 【求助】BSSE計(jì)算后如何進(jìn)行零點(diǎn)能校正？
已經(jīng)有8人回復(fù)
- 【求助】Orca 2.8的一個(gè)計(jì)算贗勢的例子報(bào)錯(cuò)
已經(jīng)有12人回復(fù)
- 【求助】BSSE計(jì)算出錯(cuò)
已經(jīng)有16人回復(fù)
- 【求助】用BSSE校正計(jì)算分子片之間的相互作用能出現(xiàn)錯(cuò)誤
已經(jīng)有3人回復(fù)
- 【求助】請(qǐng)教使用orca計(jì)算出的這個(gè)自旋密度圖是怎么畫的？
已經(jīng)有5人回復(fù)
lihb734
鐵桿木蟲 (職業(yè)作家)
- QC強(qiáng)帖: 2
- 應(yīng)助: 132 (高中生)
- 金幣: 7495.4
- 帖子: 3107
- 在線: 1177.3小時(shí)
- 蟲號(hào): 448403
2樓2011-05-06 15:55:02
ghcacj
榮譽(yù)版主 (著名寫手)
- QC強(qiáng)帖: 2
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 2.764
- 金幣: 2543.5
- 帖子: 1376
- 在線: 1199小時(shí)
- 蟲號(hào): 513281
3樓2011-05-06 18:00:43
4樓2011-05-08 07:16:58
ghcacj
榮譽(yù)版主 (著名寫手)
- QC強(qiáng)帖: 2
- 應(yīng)助: 26 (小學(xué)生)
- 貴賓: 2.764
- 金幣: 2543.5
- 帖子: 1376
- 在線: 1199小時(shí)
- 蟲號(hào): 513281
5樓2011-05-09 10:34:46
6樓2011-05-09 10:56:50
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ghcacj 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|