| 5 | 1/1 | 返回列表 |
| 查看: 1410 | 回復(fù): 8 | |||
| 當(dāng)前主題已經(jīng)存檔。 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
bents鐵蟲 (小有名氣)
|
[交流]
【求助】如何正確輸入BSSE文件
|
||
| 我的體系含有Ru, O , H ,C 原子,采用B3LYP/GENECP 基組 請(qǐng)問(wèn)如何輸入正確的 計(jì)算相互作用的文件 謝謝大家的幫助 |
木蟲 (正式寫手)

木蟲 (正式寫手)
|
用counterpoise關(guān)鍵字: #uB3LYP/6-31G Counterpoise=2 SCF=(tight) C2 + NO --> C2NO 0 2 0 1 0 2 C 0.000000 1.600184 0.000000 1 C 0.000000 0.960265 0.000000 1 N 0.208401 3.031919 0.000000 2 O 1.164068 3.666273 0.000000 2 0 2 0 1 0 2分別代表總體系電荷和多重度,片段一的電荷和多重度,片段二的電荷和多重度。 分子坐標(biāo)后的1和2是定義的片段 用Bq原子: # b3lyp/3-21G** nosymm scf=tight test HBr + H2O manual counterpoise calculation, H2O removed 0 1 H 0.685176 -0.004924 -0.026973 Br -0.771917 0.000050 0.001967 O-Bq 2.536864 -0.000136 -0.051401 H-Bq 3.015128 0.789231 0.184042 H-Bq 3.021888 -0.784986 0.185282 這就相當(dāng)于用HBr + H2O的基組計(jì)算HBr 舉個(gè)例子,比如計(jì)算AB體系的BSSE矯正: 需要得出的五個(gè)能量分別是: 第一個(gè)是整體的能量E(AB) 第二個(gè)是單獨(dú)計(jì)算A的能量E(A) 第三個(gè)是單獨(dú)計(jì)算B的能量E(B) 第四個(gè)是計(jì)算在AB基組下,A的能量E(A/AB) 第五個(gè)是計(jì)算在AB基組下,B的能量E(B/AB) 那么,BSSE energy =E(A)+E(B)-E(A/AB)-E(B/AB)。計(jì)算E(A/AB)的基函數(shù)是計(jì)算AB的基函數(shù),而計(jì)算E(A)的基函數(shù)只是計(jì)算AB基函數(shù)的一部分,基函數(shù)越多,能量計(jì)算的越準(zhǔn)確,能量計(jì)算越準(zhǔn)確,就越低,所以E(A/AB) 輸出文件里面5個(gè)scf第一個(gè)是整體的能量(記為E(AB)),后面4個(gè)是單體能量,分別是2個(gè)部分在整體基組和各自基組下的能量。 我當(dāng)時(shí)是用-Bq原子算的BSSE矯正, 估計(jì)你是用CP方法吧。但原理都是一樣的。 |

木蟲 (正式寫手)

至尊木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑22408 288分 +4 | new382 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料專碩322分 +9 | 哈哈哈吼吼吼哈 2026-04-02 | 9/450 |
|
|
[考博] 材料工程專業(yè)碩士申博 +3 | 麟正宇 2026-03-30 | 3/150 |
|
|
[考研] 266分,一志愿電氣工程,本科材料,求材料專業(yè)調(diào)劑 +4 | 哇呼哼呼哼 2026-04-02 | 4/200 |
|
|
[考研] 08開頭看過(guò)來(lái)。! +4 | wwwwffffff 2026-03-31 | 6/300 |
|
|
[考研] 385分 生物學(xué)(071000)求調(diào)劑 +6 | qf626 2026-04-01 | 6/300 |
|
|
[考研] 生物學(xué)308分求調(diào)劑(一志愿華東師大) +5 | 相信必會(huì)光芒萬(wàn)?/a> 2026-03-31 | 6/300 |
|
|
[考研] 322求調(diào)劑 +5 | 熹僖XX 2026-03-31 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 12138。、m?x?b 2026-03-26 | 4/200 |
|
|
[考研] 274求調(diào)劑 +8 | 薛定諤的虎。 2026-04-01 | 8/400 |
|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +14 | 小黃鴨寶 2026-03-30 | 14/700 |
|
|
[考研] 調(diào)劑 +5 | 好好讀書。 2026-03-28 | 7/350 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +15 | minht 2026-03-29 | 16/800 |
|
|
[考研] 一志愿華南師范361分,化學(xué)求調(diào)劑 +4 | Nicole88888 2026-04-01 | 4/200 |
|
|
[考研] 0856 335分 +9 | cccchenso 2026-03-29 | 9/450 |
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 13/650 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|