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

木蟲 (正式寫手)
|
用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的基組計算HBr 舉個例子,比如計算AB體系的BSSE矯正: 需要得出的五個能量分別是: 第一個是整體的能量E(AB) 第二個是單獨計算A的能量E(A) 第三個是單獨計算B的能量E(B) 第四個是計算在AB基組下,A的能量E(A/AB) 第五個是計算在AB基組下,B的能量E(B/AB) 那么,BSSE energy =E(A)+E(B)-E(A/AB)-E(B/AB)。計算E(A/AB)的基函數(shù)是計算AB的基函數(shù),而計算E(A)的基函數(shù)只是計算AB基函數(shù)的一部分,基函數(shù)越多,能量計算的越準(zhǔn)確,能量計算越準(zhǔn)確,就越低,所以E(A/AB) 輸出文件里面5個scf第一個是整體的能量(記為E(AB)),后面4個是單體能量,分別是2個部分在整體基組和各自基組下的能量。 我當(dāng)時是用-Bq原子算的BSSE矯正, 估計你是用CP方法吧。但原理都是一樣的。 |

至尊木蟲 (著名寫手)

|
我試驗了一下 但是輸出的提示如下 %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. 請問 為什么? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工調(diào)劑,339 +13 | 10213207 2026-03-31 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑11408 +4 | 劉小鹿lu 2026-03-27 | 4/200 |
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +19 | Youth_ 2026-03-30 | 19/950 |
|
|
[考研] 262求調(diào)劑 +4 | 勵志一定發(fā)文章 2026-04-02 | 5/250 |
|
|
[考研] 材料調(diào)劑 +11 | 一樣YWY 2026-03-31 | 11/550 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-31 | 7/350 |
|
|
[考研] 275求調(diào)劑 +12 | waltzh 2026-04-01 | 12/600 |
|
|
[考研] 求調(diào)劑 +5 | 零八# 2026-03-27 | 5/250 |
|
|
[考研] 330分求調(diào)劑 +11 | qzenlc 2026-03-29 | 11/550 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
[碩博家園] 博一被送出聯(lián)培感覺不適應(yīng)怎么辦 +3 | 全村的狗 2026-03-31 | 3/150 |
|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +5 | 雙馬尾痞老板2 2026-03-31 | 5/250 |
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +7 | 鐵憨憨123425 2026-03-31 | 7/350 |
|
|
[考研] 309分085801求調(diào)劑 +7 | 學(xué)員Gtwj7W 2026-03-31 | 7/350 |
|
|
[考研] 生物考研337分求調(diào)劑 +4 | cgxin 2026-03-30 | 6/300 |
|
|
[考研] 085600,專業(yè)課化工原理,320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進合金團隊招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|