| 5 | 1/1 | 返回列表 |
| 查看: 2186 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
孫1986金蟲 (正式寫手)
|
[交流]
高斯大俠留步,看看這個(gè)高斯輸入文件有何問題?
|
||
|
#p fopt=(ts z-matrix,calcfc) iop(5/13=1,1/11=1) freq optcyc=500 scf(maxcyc=300) ub3lyp/6-31g(d,p) geom=connectivity NO+CH2 reaction 0 2 N 0.00000000 0.00000000 -0.61797712 O 0.00000000 0.00000000 0.54072998 C 1.31900152 -1.04941686 -1.22943436 H 1.31900152 -0.05250786 -1.64471336 H 1.31900152 -2.04632586 -1.64471336 1 2 2.0 2 3 4 1.0 5 1.0 4 5 自己為高斯菜鳥,現(xiàn)在初學(xué)階段,模擬別人的結(jié)果,輸入這樣一個(gè)文件,有錯(cuò)誤,請(qǐng)哪位大俠指導(dǎo)一二 |

木蟲 (正式寫手)
木蟲 (正式寫手)
至尊木蟲 (文壇精英)

新蟲 (初入文壇)
|
該例子有一定的典型性,有很多高斯初學(xué)者會(huì)遇到的問題,在此詳細(xì)討論一下: 首先,這個(gè)輸入文件可以運(yùn)行,但是容易導(dǎo)致錯(cuò)誤意想不到的結(jié)果 (1)當(dāng)給出opt=z-matrix時(shí),高斯會(huì)去讀取輸入文件中的Z-matrix以及變量。可是樓主給出的是直角坐標(biāo)(cartesian),所以程序會(huì)按照opt=cartesian來進(jìn)行優(yōu)化。參考http://www.gaussian.com/g_tech/g_ur/k_opt.htm (When a Z-matrix without any variables is used for the molecule specification, and Opt=Z-matrix is specified, then the optimization will actually be performed in Cartesian coordinates.) 使用該方法(opt=cartesian)收斂速度遠(yuǎn)遠(yuǎn)低于默認(rèn)的冗余內(nèi)坐標(biāo)(redundant internal coordinate) (2) fopt是配合opt=z-matrix使用的,要求對(duì)所有變量進(jìn)行優(yōu)化。這個(gè)選項(xiàng)現(xiàn)在已經(jīng)很少人用了,通常直接使用opt即可。 (3)iop(5/13=1,1/11=1)。iop(5/13)=1要求哪怕SCF不收斂都繼續(xù)計(jì)算,該選項(xiàng)是用來調(diào)試用的,在實(shí)際應(yīng)用中不應(yīng)該出現(xiàn)。iop(1/11=1)即noeigentest,建議在結(jié)構(gòu)已經(jīng)很接近過渡態(tài)時(shí)使用,否則容易得到其他意想不到的虛頻,比如甲基旋轉(zhuǎn)。一般的建議是,能使用關(guān)鍵詞盡量避免使用IOP。首先IOP在多步計(jì)算中(比如opt freq)只能適用于第一步,而使用關(guān)鍵詞對(duì)所有步驟都適用。此外,在不同版本的高斯程序中,IOP參數(shù)可能會(huì)發(fā)生改變。 (4)optcyc=N, 即opt(maxcyc=N),還有scf(maxcyc=300)。如果結(jié)構(gòu)或者SCF不收斂,增大默認(rèn)的周期數(shù)通常不會(huì)有改善,反而耗費(fèi)更多的計(jì)算時(shí)間。 (5) 如果optcyc寫成opt(maxcyc=N),同時(shí)又使用opt=(ts,z-matrix),同一個(gè)opt出現(xiàn)兩次,這種情況應(yīng)該避免,因?yàn)槌绦蛟谏煞菢?biāo)準(zhǔn)輸入(即輸出文件開始部分生成的一大堆IOP選項(xiàng))時(shí)會(huì)對(duì)相同選項(xiàng)進(jìn)行加和,可能導(dǎo)致意想不到的問題。 (6)總的來說,初學(xué)高斯,建議從最少的選項(xiàng)開始,完成一個(gè)基本計(jì)算,然后添加適當(dāng)?shù)倪x項(xiàng)以實(shí)現(xiàn)額外的功能。直接從其他人那里拷貝過來輸入文件,可能適用于特定體系,但不一定在任何情況都適用。尤其是比如IOP(5/13=1),在正式計(jì)算中絕對(duì)應(yīng)該避免出現(xiàn) |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +13 | 田洪有 2026-03-26 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +15 | leaflight 2026-03-24 | 15/750 |
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +6 | minht 2026-03-29 | 6/300 |
|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 一志愿華北電力大學(xué)能動(dòng)專碩,293,求調(diào)劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研]
|
13659058978 2026-03-24 | 4/200 |
|