| 5 | 1/1 | 返回列表 |
| 查看: 1407 | 回復(fù): 8 | |||
| 當(dāng)前主題已經(jīng)存檔。 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】如何正確輸入BSSE文件
|
|||
| 我的體系含有Ru, O , H ,C 原子,采用B3LYP/GENECP 基組 請(qǐng)問如何輸入正確的 計(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方法吧。但原理都是一樣的。 |

木蟲 (正式寫手)

|
我試驗(yàn)了一下 但是輸出的提示如下 %chk=AA4_bsse.chk %mem=1700mb %nproc=4 Will use up to 4 processors via shared memory. ----------------------------------------------------------------- # oniom(b3lyp/genecp:b3lyp/3-21g) nosymm scf=tight Counterpoise=2 ----------------------------------------------------------------- Cannot combine ONIOM and Counterpoise. 請(qǐng)問 為什么? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296材料專碩求調(diào)劑 +4 | 202451007219 2026-04-02 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +10 | 舟自梗 2026-04-01 | 11/550 |
|
|
[考研] 266求調(diào)劑 +10 | 陽陽哇塞 2026-04-01 | 10/500 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +7 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 生物學(xué)327,求調(diào)劑 +5 | 書上的梅子 2026-04-01 | 6/300 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +14 | 一定上岸aaa 2026-03-27 | 18/900 |
|
|
[考研] 生物與醫(yī)藥考研調(diào)劑 +5 | 鐵憨憨123425 2026-03-31 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
[考研] 調(diào)劑申請(qǐng) +8 | 張張張張zy 2026-03-31 | 9/450 |
|
|
[考研] 311(085601)求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +6 | Ong3 2026-03-27 | 6/300 |
|
|
[考研] 085404 22408 315分 +5 | zhuangyan123 2026-03-31 | 6/300 |
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|