 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2489 | 回復(fù): 5 | |||
| 【獎勵】 本帖被評價3次，作者ghcacj增加金幣 2.4 個 | |||
| 本帖產(chǎn)生 1 個 QC強帖 ，點擊這里進(jìn)行查看 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[資源]
【gchcaj個人文集】淺談Orca計算BSSE校正
|
|||
|
今天看到lihb734兄寫了個關(guān)于BSSE校正的帖子，頗有啟發(fā)，所以我也來寫一個關(guān)于使用Orca校正BSSE的，大家看到的Gaussian的BSSE校正比較多，可以說，Gaussian在計算BSSE的input文件上比較方便，給出兩個分子的坐標(biāo)以及Counterpoise的關(guān)鍵詞即可，上手容易，但是細(xì)究中間的原理過程，我個人orca在這里好些，廢話不多說，先貼出BSSE的定義是，這個是Boys和Bernardi二位搞出來的，如公式1： 初看起來比較復(fù)雜，其實不難，表達(dá)式中的下表表示計算的片段，AB表示二聚體，A或B表示單聚體；上標(biāo)表示分子所處的優(yōu)化幾何坐標(biāo)，AB表示處于二聚體的優(yōu)化幾何坐標(biāo)，A或B表示單聚體優(yōu)化幾何坐標(biāo)；括號里表示使用的基組。這樣就可以把這個分成七項四個大項，第一項是二聚體的能量，第二項是兩個單聚體能量之和，第三項是處于二聚體優(yōu)化幾何坐標(biāo)下的使用單聚體部分的基組計算兩個單聚體的能量之和，第四項是處于二聚體優(yōu)化幾何坐標(biāo)下的使用二聚體部分的基組計算兩個單聚體的能量之和。所以公式2可以改寫為： 這里dimer表示二聚體，monomer2表示原公式的第四項，monomer1表示原公式的第三項，monomer0表示原公式的第二項。如果只計算前三項，就是被校正后的能量。 現(xiàn)在用一個例子來說明下，也是兩個水分子的作用，在Orca手冊91-93頁： Orca的常規(guī)語法這里不再解釋，只是介紹下輸入文件中的冒號作用，這個冒號是專門為這個BSSE計算考慮的，把冒號跟在原子后面，該原子的電子和核電荷均被刪除，但是保留原子以被分配基組。 這個任務(wù)其實是通過建立多個任務(wù)，再用關(guān)鍵詞%id將每大項的兩個結(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，這時候就不難計算了，ΔE=-152.655623625 -(-152.648364970+152.647006948) +152.647062118=-0.007203485 Eh 計算到了這里，應(yīng)該說Orca不如Gaussian那么智能，因為Gaussian一個關(guān)鍵詞，Orca要寫這么多，但是做計算的經(jīng)驗告訴我，越是書寫麻煩，說明自定義的空間就越大，這里就體現(xiàn)如果我們要做一個Scan的話，Orca不用一個點一個點的計算，而且可以將BSSE校正中的每項拆開，分別做Scan或者SP，最后通過公式進(jìn)行處理即可。這點Gaussian我好像發(fā)現(xiàn)不行，因為Gaussian的Counterpoise包含了太多的東西，中間的東西無法單獨調(diào)用（本人用Gaussian較少，如果說錯了，請指出，不過別罵娘就好）。這樣在計算時方便些（個人感覺） 我也來個免責(zé)聲明： 本人主業(yè)是Monte Carlo，對于QC是外行，所以這個基本是本人的一家之言，如果有問題請指出，不承擔(dān)任何學(xué)術(shù)上的責(zé)任。 [ Last edited by 375642546 on 2011-5-6 at 17:25 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化軟件學(xué)習(xí) |
» 猜你喜歡
26年儲能、電池方向博士申請
已經(jīng)有2人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有13人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有175人回復(fù)
-
青島大學(xué)化學(xué)化工學(xué)院2026年招收調(diào)劑碩士研究生
已經(jīng)有0人回復(fù)
-
北方工業(yè)大學(xué)085808儲能技術(shù)-02方向（電池）接收調(diào)劑生
已經(jīng)有0人回復(fù)
-
272求調(diào)劑，接受跨專業(yè)調(diào)劑！
已經(jīng)有3人回復(fù)
-
有一區(qū)TOP論文，求調(diào)劑
已經(jīng)有22人回復(fù)
-
大連工業(yè)大學(xué)杰青/長江團(tuán)隊-生物質(zhì)材料-儲能電池方向招收2026級博士生
已經(jīng)有14人回復(fù)
-
齊齊哈爾大學(xué)材料科學(xué)與工程學(xué)院有調(diào)劑名額
已經(jīng)有0人回復(fù)
-
2026貴州大學(xué)博士招生
已經(jīng)有0人回復(fù)
-
青島大學(xué)化學(xué)化工學(xué)院2026年招收調(diào)劑碩士研究生
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 很好的計算繪圖軟件教程
已經(jīng)有14人回復(fù)
- origin計算曲線面積-表面輪廓數(shù)據(jù)處理
已經(jīng)有12人回復(fù)
- 計算旋轉(zhuǎn)能壘用勢能面掃描之后還用找過渡態(tài)嗎？我是新手，希望大家給與幫助。[/font]
已經(jīng)有3人回復(fù)
- Materials studio 計算單點能
已經(jīng)有9人回復(fù)
- 純理論計算投稿的問題
已經(jīng)有21人回復(fù)
- |求助| 統(tǒng)計學(xué)中在進(jìn)行相關(guān)性檢驗時，P值如何計算？
已經(jīng)有9人回復(fù)
- 有關(guān)opt+freq+bsse校正輸出文件的意思
已經(jīng)有3人回復(fù)
- BSSE校正中 片段2的變化如何改寫
已經(jīng)有8人回復(fù)
- 求助DFT計算中CuCl2在γ-Al2O3 表面的吸附結(jié)構(gòu)？
已經(jīng)有5人回復(fù)
- 【求助】用Orca計算一個片簇的單點能，不收斂。
已經(jīng)有7人回復(fù)
- 【求助】BSSE計算后如何進(jìn)行零點能校正？
已經(jīng)有8人回復(fù)
- 【求助】Orca 2.8的一個計算贗勢的例子報錯
已經(jīng)有12人回復(fù)
- 【求助】BSSE計算出錯
已經(jīng)有16人回復(fù)
- 【求助】用BSSE校正計算分子片之間的相互作用能出現(xiàn)錯誤
已經(jīng)有3人回復(fù)
- 【求助】請教使用orca計算出的這個自旋密度圖是怎么畫的？
已經(jīng)有5人回復(fù)
3樓2011-05-06 18:00:43
2樓2011-05-06 15:55:02
4樓2011-05-08 07:16:58
5樓2011-05-09 10:34:46
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程322分 +5 | 哈哈哈吼吼吼哈 2026-04-01 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-04-01 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +7 | 莫xiao 2026-04-01 | 7/350 |
|
|
|
[考研] 070300一志愿211，312分求調(diào)劑院校 +14 | 小黃鴨寶 2026-03-30 | 14/700 |
|
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +10 | 腳滑的守法公民 2026-03-27 | 10/500 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +6 | 李瀟可 2026-03-26 | 6/300 |
|
|
|
[考研] 336求調(diào)劑 +3 | kiyy 2026-04-01 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +11 | Amber00 2026-03-31 | 11/550 |
|
|
|
[考研] 土木304求調(diào)劑 +5 | 兔突突突， 2026-03-31 | 6/300 |
|
|
|
[考研] 求調(diào)劑：一志愿：南京大學(xué) 專業(yè)：0705 總分320 ，本科985，四六級已過 +3 | lfy760306 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿華東師范大學(xué)有機化學(xué)專業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +9 | 真名有冰 2026-03-29 | 10/500 |
|
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
|
[考研] 354求調(diào)劑 +3 | lxb598 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑 六級521 +4 | 愛吃大鴨梨 2026-03-31 | 4/200 |
|
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 求佛 +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|