| 9 | 1/1 | 返回列表 |
| 查看: 1406 | 回復(fù): 8 | |||
| 當(dāng)前主題已經(jī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. 請問 為什么? |
木蟲 (正式寫手)

木蟲 (正式寫手)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 318求調(diào)劑,計算材料方向 +7 | 吸喵有害笙命 2026-04-01 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +9 | 一樣YWY 2026-04-01 | 9/450 |
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +12 | 睿08 2026-03-27 | 12/600 |
|
|
[考研] 279求調(diào)劑 +7 | 莫xiao 2026-04-01 | 7/350 |
|
|
[考研] 292求調(diào)劑 +17 | 木蟲er12138 2026-04-01 | 17/850 |
|
|
[考研] 270求調(diào)劑 +7 | 小杰pp 2026-03-31 | 8/400 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +7 | 月月鳥木 2026-04-01 | 7/350 |
|
|
[考研] 318求調(diào)劑 +8 | 七憶77 2026-04-01 | 8/400 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +15 | minht 2026-03-29 | 16/800 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +17 | minht 2026-03-26 | 17/850 |
|
|
[考研] 0805求調(diào)劑 +6 | 是水分 2026-03-31 | 6/300 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 4/200 |
|
|
[考研] 318求調(diào)劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
[考研] 263求調(diào)劑 +3 | DDDDuu 2026-03-27 | 3/150 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|