| 29 | 1/1 | 返回列表 |
| 查看: 2731 | 回復(fù): 28 | |||
[交流]
如何使用CBS-QB3尋找過渡態(tài)
|
|||
| 已在其他水平下找到了過渡態(tài),打算在此基礎(chǔ)上使用CBS-QB3水平下通過IRCmax尋找過渡態(tài),不知道如何寫輸入文件。想請教一個CBS-QB3做IRCmax的輸入文件。 |
» 搶金幣啦!回帖就可以得到:
+1/905
+2/412
+5/90
+1/84
+1/83
+1/80
+2/58
+1/33
+1/30
+1/18
+1/18
+2/16
+1/14
+1/11
+2/10
+1/9
+1/9
+1/6
+1/3
+1/1
版主 (著名寫手)
|
這個沒法一步完成。 一般來說IRCMax寫法如下: IRCMax(方法2:方法1) 含義是是在方法1的IRC途徑上 用方法2在每一個點(diǎn)都做單點(diǎn)能計算,找出其中能量最高的點(diǎn)。 但是CBS-QB3是一個多步驟的方法 其中還包含了第一步一個opt步驟 所以這個不能直接寫的。 高斯說明書的原話也是這樣說的 : Analytic gradients are required for the IRC portion of the calculation (model1 above). Any non-compound energy method and basis set may be used for model2. CBS-QB3就是屬于 compound energy method。 所以沒法直接寫。如果需要就只能分步驟人工計算 第一步 先用 方法1 算IRC , 第二步 自己提取IRC上每一個點(diǎn)的結(jié)構(gòu) 逐一用CBS-QB3做一次單點(diǎn)能 (如果是Gaussian16 就支持CBS(NoOpt)能自動跳過其中的opt步驟 ,給出能量)。但是G09不支持這個寫法,那么就要人工來做整個CBS步驟 人為跳過opt步驟。 CBS-QB3除去優(yōu)化的話需要進(jìn)行四步驟 1. #p B3LYP/CBSB7 Freq 2. #p CCSD(T)/6-31+G(d') 3. #p Geom=AllCheck Guess=Read MP4SDQ/CBSB4 4. #P Geom=AllCheck Guess=Read MP2/CBSB3 CBSExtrap=(NMin=10,MinPop) 你可以把這四部分的計算都用--link1--來隔開 寫在同一個文件里。 http://www.chem.ucla.edu/~mccarren/CBSQB3_worksheet/ 這個網(wǎng)站可以幫你按照CBS-qb3的公式分別取這幾個計算中有用的部分加加減減 得到最后的能量。 然后對每一個點(diǎn)都要做一次這些運(yùn)算。整個運(yùn)算步驟還是挺多的。 |
版主 (著名寫手)
版主 (著名寫手)
版主 (著名寫手)
送紅花一朵 |
前輩,這是今天的IRCmax的輸入輸出文件。我B2LYP/6-311++(d,p)做的測試。 input文件 %chk=RE-W12IRCmax.chk SCF=(MaxCycle=100,xqc) #p IRCmax(B3LYP/6-311++G(d,p):B3LYP/cc-pVTZ) NoSymm RE-W8-W8 TS Search 0 1 O H 1 0.96474388 H 1 2.52324319 2 96.16865489 C 1 1.37453100 2 110.36391138 3 -68.07987935 0 O 4 1.18920427 1 124.02515701 2 -96.13403058 0 S 4 1.79329724 1 114.57558758 5 -177.11988964 0 1 2 1.0 4 1.0 2 3 6 1.0 4 5 2.0 6 1.0 5 6 output 文件 %NProcShared=12 Will use up to 12 processors via shared memory. %LindaWorker=ys0716,yg0926,ys1047,ys1147,ys1438,ys0428,ys1521,ys1112, %chk=RE-W12IRCmax.chk SetLPE: input flags="" SetLPE: new flags=" -nodelist 'ys0716 yg0926 ys1047 ys1147 ys1438 ys0428 ys1521 ys1112'" Will use up to 8 processors via Linda. --------------------------------------------------- #p IRCmax(B3LYP/6-311++G(d,p):B3LYP/cc-pVTZ) NoSymm --------------------------------------------------- You must use ReadCartesianFC, FCCards, or CalcFC with IRC. Error termination via Lnk1e in /home-yw/Soft/Gaussian.D01.Linda/g09/l1.exe at Thu May 23 01:19:37 2019. Job cpu time: 0 days 0 hours 0 minutes 1.5 seconds. File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1 以前在超算上面計算IRC時候,我都是直接從過渡態(tài)結(jié)構(gòu)直接進(jìn)行計算,沒有用過chk文件。 |
版主 (著名寫手)
|
還是有問題。 %chk=RE-W12IRCmax.chk SCF=(MaxCycle=100,xqc) #p IRCmax(calcfc,B3LYP/6-311++G(d,p):B3LYP/cc-pVTZ) NoSymm RE-W8-W8 TS Search 0 1 O H 1 0.96474388 H 1 2.52324319 2 96.16865489 C 1 1.37453100 2 110.36391138 3 -68.07987935 0 O 4 1.18920427 1 124.02515701 2 -96.13403058 0 S 4 1.79329724 1 114.57558758 5 -177.11988964 0 1 2 1.0 4 1.0 2 3 6 1.0 4 5 2.0 6 1.0 5 6 輸出文件 %NProcShared=12 Will use up to 12 processors via shared memory. %LindaWorker=yg0923,yg1048,ys1030,ys0736,yg0935,ys0730,ys0934,yg0946, %chk=RE-W12IRCmax.chk SetLPE: input flags="" SetLPE: new flags=" -nodelist 'yg0923 yg1048 ys1030 ys0736 yg0935 ys0730 ys0934 yg0946'" Will use up to 8 processors via Linda. ---------------------------------------------------------- #p IRCmax(calcfc,B3LYP/6-311++G(d,p):B3LYP/cc-pVTZ) NoSymm ---------------------------------------------------------- 1/10=4,14=-1,18=40,26=3,29=120000,38=1,44=3,58=1/1,15; 2/15=1,29=3/2; 3/5=4,6=6,7=1111,11=2,16=1,25=1,30=1,74=-5/1,2,3; 4//1; 5/5=2,38=5/2; 6/7=2,8=2,9=2,10=2,28=1/1; 1/10=4,14=-1,26=3,44=3,48=2,58=1/15; 3/5=16,6=1,11=2,16=1,25=1,30=1,71=2,74=-5,140=1/1,2,3; 4//1; 5/5=2,38=5/2; 8/6=4,10=90,11=11/1; 11/6=1,8=1,9=11,15=111,16=1,31=1/1,2,10; 10/6=1,31=1/2; 6/7=2,8=2,9=2,10=2,28=1/1; 7/8=1,10=1,18=20,25=1,30=1/1,2,3,16; 1/10=4,14=-1,26=3,44=3,58=1/15(1); 99/5=20/99; 2/15=1,29=3/2; 3/5=16,6=1,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3; 4/5=1/1; 5/5=2/2; 7/30=1/1,2,3,16; 1/14=-1,26=3,44=3,48=3,58=1/15(-5); 2/15=1,29=3/2; 3/5=4,6=6,7=1111,11=2,16=1,25=1,30=1,74=-5/1,2,3; 4/5=1/1; 5/5=2/2; 1/14=-1,26=3,44=3,48=2,58=1/15; 3/5=16,6=1,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3; 4/5=1/1; 5/5=2/2; 7/8=1,30=1,45=1/1,2,3,16; 1/14=-1,26=3,44=3,48=1,58=1/15(-15); 2/15=1,29=3/2; 99/5=20,9=1/99; Leave Link 1 at Thu May 23 17:20:46 2019, MaxMem= 0 cpu: 1.2 (Enter /home-yw/Soft/Gaussian.D01.Linda/g09/l101.exe) ------------------ RE-W8-W8 TS Search ------------------ Symbolic Z-matrix: Charge = 0 Multiplicity = 1 O H 1 0.96474 H 1 2.52324 2 96.16865 C 1 1.37453 2 110.36391 3 -68.07988 0 O 4 1.1892 1 124.02516 2 -96.13403 0 S 4 1.7933 1 114.57559 5 -177.11989 0 NAtoms= 6 NQM= 6 NQMF= 0 NMMI= 0 NMMIF= 0 NMic= 0 NMicF= 0. Isotopes and Nuclear Properties: (Nuclear quadrupole moments (NQMom) in fm**2, nuclear magnetic moments (NMagM) in nuclear magnetons) Atom 1 2 3 4 5 6 IAtWgt= 16 1 1 12 16 32 AtmWgt= 15.9949146 1.0078250 1.0078250 12.0000000 15.9949146 31.9720718 NucSpn= 0 1 1 0 0 0 AtZEff= 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 NQMom= 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 NMagM= 0.0000000 2.7928460 2.7928460 0.0000000 0.0000000 0.0000000 AtZNuc= 8.0000000 1.0000000 1.0000000 6.0000000 8.0000000 16.0000000 Leave Link 101 at Thu May 23 17:20:46 2019, MaxMem= 33554432 cpu: 1.5 (Enter /home-yw/Soft/Gaussian.D01.Linda/g09/l115.exe) IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC Berny reaction coordinate optimization. There are no variables to optimize! Error termination via Lnk1e in /home-yw/Soft/Gaussian.D01.Linda/g09/l115.exe at Thu May 23 17:20:46 2019. Job cpu time: 0 days 0 hours 0 minutes 2.8 seconds. File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1 |
| 29 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請] 面上項(xiàng)目還需要AI說明嗎? +3 | liyundong 2026-03-08 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑 +4 | 調(diào)劑的考研學(xué)生 2026-03-09 | 4/200 |
|
|
[考研] 考研一志愿長安大學(xué)材料與化工309分請求調(diào)劑 +5 | dtdxzxx 2026-03-06 | 7/350 |
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
5+3
|
蘇堯幺幺 2026-03-03 | 9/450 |
|
|
[考研] 材料與化工求調(diào)劑 +8 | 與冬清寧 2026-03-07 | 10/500 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實(shí)人 2026-03-09 | 4/200 |
|
|
[考博] 申博 +5 | 添菜了哈 2026-03-04 | 8/400 |
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
[考研] 083000,總分284,求調(diào)劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|
[考研] 化學(xué)工程求調(diào)劑 +12 | 化工人999 2026-03-04 | 12/600 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +9 | 土土小蟲 2026-03-03 | 11/550 |
|
|
[考研] 0856材料專碩274能調(diào)劑去哪里? +3 | 22735 2026-03-04 | 4/200 |
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +8 | 拾玖壹 2026-03-04 | 8/400 |
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
[考研] 0856材料工程,初試313調(diào)劑 +7 | 賣個關(guān)子吧 2026-03-03 | 7/350 |
|