| 23 | 1/1 | 返回列表 |
| 查看: 4370 | 回復: 22 | |||||
[交流]
【求助】鍵長鍵角等固定后的優(yōu)化..懸賞30...【已經(jīng)解決~感謝~】
|
|||||
|
新手問題........... 現(xiàn)在正在找一個過渡態(tài),沒有固定鍵長鍵角之前,構型偏離較大最后就不能收斂,或者分子散了,若構型偏離較小就會優(yōu)化直接到穩(wěn)定態(tài) 現(xiàn)在準備固定某個鍵長后繼續(xù)優(yōu)化,用gaussview建立模型之后,用自帶的redundant,然后手動再修改某幾個數(shù)值,輸入文件中有: D 4 1 2 3 -120.000000 F 這樣的描述,計算過程中也沒有出錯,但是最后鍵角仍舊偏離較大,鍵長也是如此 想請教一下這個如果固定鍵長之后優(yōu)化,gaussian是在這個輸入結構附近尋找收斂點么?若在這個點附近找不到最小點是不是就會破壞這個限制?加上了鍵長鍵角限制之后和沒有加鍵長鍵角有什么區(qū)別呢? 或者....悲情的就是固定方法不對,輸入文件不對? 求大神詳解...或者能麻煩提供詳解的鏈接地址么?.... [ Last edited by 亍寞 on 2011-3-12 at 21:17 ] |
量化計算有意義的討論 | Gaussian |
» 搶金幣啦!回帖就可以得到:
+2/124
+1/79
+1/65
+2/56
+1/40
+1/38
+1/37
+1/32
+1/17
+1/12
+1/11
+1/10
+1/9
+1/7
+1/5
+1/4
+1/3
+1/2
+1/2
+1/1
|
謝謝,現(xiàn)在您所說的主要的反應坐標已經(jīng)找到了,是鍵長,但是就是再固定的時候除了點問題,固定之后計算,但是計算結果中這個希望固定的鍵長你卻恢復了穩(wěn)定態(tài)的鍵長,因為以前沒有固定過鍵長進行計算,所以也不知道是不是固定方法有問題 這是部分輸入文件:(z-matrix) %mem=3500MB %nprocshared=2 #p opt=(modredundant,maxcyc=200) freq b3lyp/genecp scf=tight 鍵說明部分。。。略 B1 4.900000 基組說明部分。。。略 請問這個是一定要用popt命令么?那么popt(參數(shù))中的參數(shù)又應該怎么設置比較好?您所說的反應坐標的逐漸改變是用scan命令還是人為的固定不同鍵長? 十分感謝您的回復! |
|
謝謝回復,這個原理我明白,是輸入文件的問題,麻煩您看下,這是我剛才用乙烷分子做了嘗試,想固定兩個碳原子之間的距離之后優(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 原子就是兩個碳原子) 計算過程中沒有出錯,但是在最后看輸出文件的時候發(fā)現(xiàn)鍵長又回去了,不知道應該怎么輸入,像上面那個輸入文件,應該怎么改可以改對呢? |
|
乙烷你可以試試這個: #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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學材料與化工 264分各科過A區(qū)國家線 +10 | 哈哈157349 2026-03-21 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 348求調劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 321求調劑 +4 | 璞玉~~ 2026-03-25 | 4/200 |
|
|
[考研] 085600材料與化工306 +10 | z1z2z3879 2026-03-21 | 11/550 |
|
|
[考研] 求調劑 +3 | 張zz111 2026-03-27 | 4/200 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +3 | 辰沐5211314 2026-03-26 | 3/150 |
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調劑 +3 | 美味可樂雞翅 2026-03-26 | 3/150 |
|
|
[考研] 284求調劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 351求調劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 324求調劑 +4 | wysyjs25 2026-03-21 | 4/200 |
|
|
[考研] 化學工程085602 305分求調劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 環(huán)境專碩324分求調劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 291 求調劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
[考研] 材料與化工304求B區(qū)調劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 考研化學308分求調劑 +10 | 你好明天你好 2026-03-23 | 11/550 |
|
|
[考研] 調劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 求調劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
[考研] 315分,誠求調劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 考研調劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|