| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2183 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
孫1986金蟲(chóng) (正式寫(xiě)手)
|
[交流]
高斯大俠留步,看看這個(gè)高斯輸入文件有何問(wèn)題?
|
||
|
#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 自己為高斯菜鳥(niǎo),現(xiàn)在初學(xué)階段,模擬別人的結(jié)果,輸入這樣一個(gè)文件,有錯(cuò)誤,請(qǐng)哪位大俠指導(dǎo)一二 |

至尊木蟲(chóng) (文壇精英)

木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
新蟲(chóng) (初入文壇)
|
該例子有一定的典型性,有很多高斯初學(xué)者會(huì)遇到的問(wèn)題,在此詳細(xì)討論一下: 首先,這個(gè)輸入文件可以運(yùn)行,但是容易導(dǎo)致錯(cuò)誤意想不到的結(jié)果 (1)當(dāng)給出opt=z-matrix時(shí),高斯會(huì)去讀取輸入文件中的Z-matrix以及變量?墒菢侵鹘o出的是直角坐標(biāo)(cartesian),所以程序會(huì)按照opt=cartesian來(lái)進(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)是用來(lái)調(diào)試用的,在實(shí)際應(yīng)用中不應(yīng)該出現(xiàn)。iop(1/11=1)即noeigentest,建議在結(jié)構(gòu)已經(jīng)很接近過(guò)渡態(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寫(xiě)成opt(maxcyc=N),同時(shí)又使用opt=(ts,z-matrix),同一個(gè)opt出現(xiàn)兩次,這種情況應(yīng)該避免,因?yàn)槌绦蛟谏煞菢?biāo)準(zhǔn)輸入(即輸出文件開(kāi)始部分生成的一大堆IOP選項(xiàng))時(shí)會(huì)對(duì)相同選項(xiàng)進(jìn)行加和,可能導(dǎo)致意想不到的問(wèn)題。 (6)總的來(lái)說(shuō),初學(xué)高斯,建議從最少的選項(xiàng)開(kāi)始,完成一個(gè)基本計(jì)算,然后添加適當(dāng)?shù)倪x項(xiàng)以實(shí)現(xiàn)額外的功能。直接從其他人那里拷貝過(guò)來(lái)輸入文件,可能適用于特定體系,但不一定在任何情況都適用。尤其是比如IOP(5/13=1),在正式計(jì)算中絕對(duì)應(yīng)該避免出現(xiàn) |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +10 | jhybd 2026-03-23 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 281求調(diào)劑 +4 | 亞克西good 2026-03-26 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +5 | 傾____蕭 2026-03-27 | 6/300 |
|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +4 | 篤行致遠(yuǎn). 2026-03-28 | 4/200 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專(zhuān)業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|