| 5 | 1/1 | 返回列表 |
| 查看: 2190 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
孫1986金蟲 (正式寫手)
|
[交流]
高斯大俠留步,看看這個高斯輸入文件有何問題?
|
||
|
#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é)果,輸入這樣一個文件,有錯誤,請哪位大俠指導(dǎo)一二 |

新蟲 (初入文壇)
|
該例子有一定的典型性,有很多高斯初學(xué)者會遇到的問題,在此詳細討論一下: 首先,這個輸入文件可以運行,但是容易導(dǎo)致錯誤意想不到的結(jié)果 (1)當(dāng)給出opt=z-matrix時,高斯會去讀取輸入文件中的Z-matrix以及變量。可是樓主給出的是直角坐標(biāo)(cartesian),所以程序會按照opt=cartesian來進行優(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)收斂速度遠遠低于默認(rèn)的冗余內(nèi)坐標(biāo)(redundant internal coordinate) (2) fopt是配合opt=z-matrix使用的,要求對所有變量進行優(yōu)化。這個選項現(xiàn)在已經(jīng)很少人用了,通常直接使用opt即可。 (3)iop(5/13=1,1/11=1)。iop(5/13)=1要求哪怕SCF不收斂都繼續(xù)計算,該選項是用來調(diào)試用的,在實際應(yīng)用中不應(yīng)該出現(xiàn)。iop(1/11=1)即noeigentest,建議在結(jié)構(gòu)已經(jīng)很接近過渡態(tài)時使用,否則容易得到其他意想不到的虛頻,比如甲基旋轉(zhuǎn)。一般的建議是,能使用關(guān)鍵詞盡量避免使用IOP。首先IOP在多步計算中(比如opt freq)只能適用于第一步,而使用關(guān)鍵詞對所有步驟都適用。此外,在不同版本的高斯程序中,IOP參數(shù)可能會發(fā)生改變。 (4)optcyc=N, 即opt(maxcyc=N),還有scf(maxcyc=300)。如果結(jié)構(gòu)或者SCF不收斂,增大默認(rèn)的周期數(shù)通常不會有改善,反而耗費更多的計算時間。 (5) 如果optcyc寫成opt(maxcyc=N),同時又使用opt=(ts,z-matrix),同一個opt出現(xiàn)兩次,這種情況應(yīng)該避免,因為程序在生成非標(biāo)準(zhǔn)輸入(即輸出文件開始部分生成的一大堆IOP選項)時會對相同選項進行加和,可能導(dǎo)致意想不到的問題。 (6)總的來說,初學(xué)高斯,建議從最少的選項開始,完成一個基本計算,然后添加適當(dāng)?shù)倪x項以實現(xiàn)額外的功能。直接從其他人那里拷貝過來輸入文件,可能適用于特定體系,但不一定在任何情況都適用。尤其是比如IOP(5/13=1),在正式計算中絕對應(yīng)該避免出現(xiàn) |
木蟲 (正式寫手)
木蟲 (正式寫手)
至尊木蟲 (文壇精英)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +3 | 學(xué)而思兮知 2026-04-01 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 一志愿哈工大 材料與化工 +28 | 洛神哥哥 2026-03-31 | 28/1400 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +21 | 荷包蛋hyj 2026-03-26 | 21/1050 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +22 | 材料學(xué)257求調(diào)劑 2026-03-28 | 23/1150 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +15 | 柳溪邊 2026-03-26 | 16/800 |
|
|
[考研] 0710生物學(xué)考研調(diào)劑 +3 | 李多米lee. 2026-03-27 | 4/200 |
|
|
[考研] 求生物學(xué)調(diào)劑 +7 | 15172915737 2026-04-01 | 7/350 |
|
|
[考研] 08工科,295,接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-31 | 6/300 |
|
|
[考研] 生物學(xué)308分求調(diào)劑(一志愿華東師大) +3 | 相信必會光芒萬?/a> 2026-03-31 | 3/150 |
|
|
[考研] 277跪求調(diào)劑 +8 | 1915668 2026-03-27 | 13/650 |
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +12 | Li李魚 2026-03-26 | 12/600 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|