 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2464 | 回復: 5 | |||
| 【獎勵】 本帖被評價3次，作者ghcacj增加金幣 2.4 個 | |||
| 本帖產(chǎn)生 1 個 QC強帖 ，點擊這里進行查看 | |||
[資源]
【gchcaj個人文集】淺談Orca計算BSSE校正
|
|||
|
今天看到lihb734兄寫了個關(guān)于BSSE校正的帖子，頗有啟發(fā)，所以我也來寫一個關(guān)于使用Orca校正BSSE的，大家看到的Gaussian的BSSE校正比較多，可以說，Gaussian在計算BSSE的input文件上比較方便，給出兩個分子的坐標以及Counterpoise的關(guān)鍵詞即可，上手容易，但是細究中間的原理過程，我個人orca在這里好些，廢話不多說，先貼出BSSE的定義是，這個是Boys和Bernardi二位搞出來的，如公式1： 初看起來比較復雜，其實不難，表達式中的下表表示計算的片段，AB表示二聚體，A或B表示單聚體；上標表示分子所處的優(yōu)化幾何坐標，AB表示處于二聚體的優(yōu)化幾何坐標，A或B表示單聚體優(yōu)化幾何坐標；括號里表示使用的基組。這樣就可以把這個分成七項四個大項，第一項是二聚體的能量，第二項是兩個單聚體能量之和，第三項是處于二聚體優(yōu)化幾何坐標下的使用單聚體部分的基組計算兩個單聚體的能量之和，第四項是處于二聚體優(yōu)化幾何坐標下的使用二聚體部分的基組計算兩個單聚體的能量之和。所以公式2可以改寫為： 這里dimer表示二聚體，monomer2表示原公式的第四項，monomer1表示原公式的第三項，monomer0表示原公式的第二項。如果只計算前三項，就是被校正后的能量。 現(xiàn)在用一個例子來說明下，也是兩個水分子的作用，在Orca手冊91-93頁： Orca的常規(guī)語法這里不再解釋，只是介紹下輸入文件中的冒號作用，這個冒號是專門為這個BSSE計算考慮的，把冒號跟在原子后面，該原子的電子和核電荷均被刪除，但是保留原子以被分配基組。 這個任務其實是通過建立多個任務，再用關(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 計算到了這里，應該說Orca不如Gaussian那么智能，因為Gaussian一個關(guān)鍵詞，Orca要寫這么多，但是做計算的經(jīng)驗告訴我，越是書寫麻煩，說明自定義的空間就越大，這里就體現(xiàn)如果我們要做一個Scan的話，Orca不用一個點一個點的計算，而且可以將BSSE校正中的每項拆開，分別做Scan或者SP，最后通過公式進行處理即可。這點Gaussian我好像發(fā)現(xiàn)不行，因為Gaussian的Counterpoise包含了太多的東西，中間的東西無法單獨調(diào)用（本人用Gaussian較少，如果說錯了，請指出，不過別罵娘就好）。這樣在計算時方便些（個人感覺） 我也來個免責聲明： 本人主業(yè)是Monte Carlo，對于QC是外行，所以這個基本是本人的一家之言，如果有問題請指出，不承擔任何學術(shù)上的責任。 [ Last edited by 375642546 on 2011-5-6 at 17:25 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 量化軟件學習 |
» 猜你喜歡
材料/化學相關(guān)專業(yè)2026級學術(shù)/專業(yè)型碩士研究生
已經(jīng)有19人回復
-
求催化作用基礎(chǔ)第三版 （甄開吉) 課件
已經(jīng)有3人回復
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有162人回復
-
26年材料/化學相關(guān)專業(yè)博士研究生招生 溫州大學化學與材料工程學院(26年秋季入學)
已經(jīng)有0人回復
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經(jīng)有0人回復
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經(jīng)有0人回復
-
304求調(diào)劑
已經(jīng)有3人回復
-
上海華東理工大學機械與動力工程學院先進材料與制造研究組招聘公告
已經(jīng)有0人回復
-
山西大學資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有2人回復
-
304求調(diào)劑
已經(jīng)有3人回復
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 很好的計算繪圖軟件教程
已經(jīng)有14人回復
- origin計算曲線面積-表面輪廓數(shù)據(jù)處理
已經(jīng)有12人回復
- 計算旋轉(zhuǎn)能壘用勢能面掃描之后還用找過渡態(tài)嗎？我是新手，希望大家給與幫助。[/font]
已經(jīng)有3人回復
- Materials studio 計算單點能
已經(jīng)有9人回復
- 純理論計算投稿的問題
已經(jīng)有21人回復
- |求助| 統(tǒng)計學中在進行相關(guān)性檢驗時，P值如何計算？
已經(jīng)有9人回復
- 有關(guān)opt+freq+bsse校正輸出文件的意思
已經(jīng)有3人回復
- BSSE校正中 片段2的變化如何改寫
已經(jīng)有8人回復
- 求助DFT計算中CuCl2在γ-Al2O3 表面的吸附結(jié)構(gòu)？
已經(jīng)有5人回復
- 【求助】用Orca計算一個片簇的單點能，不收斂。
已經(jīng)有7人回復
- 【求助】BSSE計算后如何進行零點能校正？
已經(jīng)有8人回復
- 【求助】Orca 2.8的一個計算贗勢的例子報錯
已經(jīng)有12人回復
- 【求助】BSSE計算出錯
已經(jīng)有16人回復
- 【求助】用BSSE校正計算分子片之間的相互作用能出現(xiàn)錯誤
已經(jīng)有3人回復
- 【求助】請教使用orca計算出的這個自旋密度圖是怎么畫的？
已經(jīng)有5人回復
2樓2011-05-06 15:55:02
3樓2011-05-06 18:00:43
4樓2011-05-08 07:16:58
5樓2011-05-09 10:34:46
6樓2011-05-09 10:56:50
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ghcacj 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +17 | rare12345 2026-03-18 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 0856調(diào)劑，是學校就去 +6 | sllhht 2026-03-19 | 7/350 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +13 | 陌の森林 2026-03-18 | 13/650 |
|
|
|
[考研] 307求調(diào)劑 +8 | 冷笙123 2026-03-17 | 8/400 |
|
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 12/600 |
|
|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 生物學071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 0703化學336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|