| 23 | 1/1 | 返回列表 |
| 查看: 4368 | 回復(fù): 22 | |||||
[交流]
【求助】鍵長鍵角等固定后的優(yōu)化..懸賞30...【已經(jīng)解決~感謝~】
|
|||||
|
新手問題........... 現(xiàn)在正在找一個(gè)過渡態(tài),沒有固定鍵長鍵角之前,構(gòu)型偏離較大最后就不能收斂,或者分子散了,若構(gòu)型偏離較小就會(huì)優(yōu)化直接到穩(wěn)定態(tài) 現(xiàn)在準(zhǔn)備固定某個(gè)鍵長后繼續(xù)優(yōu)化,用gaussview建立模型之后,用自帶的redundant,然后手動(dòng)再修改某幾個(gè)數(shù)值,輸入文件中有: D 4 1 2 3 -120.000000 F 這樣的描述,計(jì)算過程中也沒有出錯(cuò),但是最后鍵角仍舊偏離較大,鍵長也是如此 想請(qǐng)教一下這個(gè)如果固定鍵長之后優(yōu)化,gaussian是在這個(gè)輸入結(jié)構(gòu)附近尋找收斂點(diǎn)么?若在這個(gè)點(diǎn)附近找不到最小點(diǎn)是不是就會(huì)破壞這個(gè)限制?加上了鍵長鍵角限制之后和沒有加鍵長鍵角有什么區(qū)別呢? 或者....悲情的就是固定方法不對(duì),輸入文件不對(duì)? 求大神詳解...或者能麻煩提供詳解的鏈接地址么?.... [ Last edited by 亍寞 on 2011-3-12 at 21:17 ] |
量化計(jì)算有意義的討論 | Gaussian |
» 搶金幣啦!回帖就可以得到:
+1/297
+1/270
+1/84
+2/78
+1/76
+1/46
+1/31
+1/16
+3/13
+1/13
+1/11
+1/9
+1/8
+1/7
+1/7
+1/7
+1/6
+1/5
+1/5
+1/1
|
謝謝,現(xiàn)在您所說的主要的反應(yīng)坐標(biāo)已經(jīng)找到了,是鍵長,但是就是再固定的時(shí)候除了點(diǎn)問題,固定之后計(jì)算,但是計(jì)算結(jié)果中這個(gè)希望固定的鍵長你卻恢復(fù)了穩(wěn)定態(tài)的鍵長,因?yàn)橐郧皼]有固定過鍵長進(jìn)行計(jì)算,所以也不知道是不是固定方法有問題 這是部分輸入文件:(z-matrix) %mem=3500MB %nprocshared=2 #p opt=(modredundant,maxcyc=200) freq b3lyp/genecp scf=tight 鍵說明部分。。。略 B1 4.900000 基組說明部分。。。略 請(qǐng)問這個(gè)是一定要用popt命令么?那么popt(參數(shù))中的參數(shù)又應(yīng)該怎么設(shè)置比較好?您所說的反應(yīng)坐標(biāo)的逐漸改變是用scan命令還是人為的固定不同鍵長? 十分感謝您的回復(fù)! |
|
謝謝回復(fù),這個(gè)原理我明白,是輸入文件的問題,麻煩您看下,這是我剛才用乙烷分子做了嘗試,想固定兩個(gè)碳原子之間的距離之后優(yōu)化,這是輸入文件 %mem=3200MB %nprocshared=2 #p opt(modredundant) freq b3lyp/6-31g scf=tight Title Card Required 0 1 C H 1 B1 H 1 B2 2 A1 H 1 B3 2 A2 3 D1 C 1 B4 2 A3 3 D2 H 5 B5 1 A4 2 D3 H 5 B6 1 A5 2 D4 H 5 B7 1 A6 2 D5 B1 1.07000000 B2 1.07000000 B3 1.07000000 B4 2.00000000 B5 1.07000000 B6 1.07000000 B7 1.07000000 A1 109.47120255 A2 109.47123134 A3 109.47120255 A4 109.47120255 A5 109.47120255 A6 109.47123134 D1 119.99999260 D2 -120.00001480 D3 179.88890000 D4 -60.11108520 D5 59.88890740 constant: B 4 5 2.00000000 (4 5 原子就是兩個(gè)碳原子) 計(jì)算過程中沒有出錯(cuò),但是在最后看輸出文件的時(shí)候發(fā)現(xiàn)鍵長又回去了,不知道應(yīng)該怎么輸入,像上面那個(gè)輸入文件,應(yīng)該怎么改可以改對(duì)呢? |
|
乙烷你可以試試這個(gè): #p opt(modredundant) b3lyp/6-31g scf=tight Title Card Required 0 1 C 0.000000 0.000000 0.000000 C 0.000000 0.000000 1.450000 H 1.026719 0.000000 1.813000 H -0.513360 0.889165 1.813000 H -0.513360 -0.889165 1.813000 H -1.026719 0.000000 -0.363000 H 0.513360 0.889165 -0.363000 H 0.513360 -0.889165 -0.363000 B 1 2 2.0 F 空行 |
| 23 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 7/350 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +4 | 北道巷 2026-03-24 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 299求調(diào)劑 +7 | shxchem 2026-03-20 | 9/450 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +3 | 鋼鐵大炮 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | 林之夕 2026-03-24 | 5/250 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|