| 5 | 1/1 | 返回列表 |
| 查看: 2341 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
溶劑化效應(yīng)輸入
|
||
|
請問溶劑化效應(yīng)的輸入命令scrf=(pcm,solvent=CH2Cl2,read)和scrf=(pcm,solvent=CH2Cl2)有什么區(qū)別? 在優(yōu)化好的構(gòu)型的基礎(chǔ)上用大基組進行計算,分別用了上面兩個命令,發(fā)現(xiàn)得出的能量,以及相關(guān)的自旋密度等數(shù)據(jù)不一樣,請問這是怎么回事? 輸入文件如下: # ub3lyp/gen scrf=(pcm,solvent=CH2Cl2,read) pop=full gfoldprint Title Card Required 0 5 坐標(biāo)數(shù)據(jù) C F N H O Se S 0 6-311+g(d,p) **** Fe 0 S 6 1.00 257539.0000000 0.0002900 38636.9000000 0.0022600 8891.4400000 0.0115200 2544.0100000 0.0456600 844.7770000 0.1403500 312.5270000 0.3142000 S 2 1.00 125.5930000 0.4087800 53.4987000 0.2116300 S 1 1.00 17.7151000 1.0000000 S 1 1.00 7.3767700 1.0000000 S 1 1.00 2.0184700 1.0000000 S 1 1.00 0.7799350 1.0000000 S 1 1.00 0.1142200 1.0000000 S 1 1.00 0.0418890 1.0000000 P 3 1.00 1678.4000000 0.0024900 396.3920000 0.0201500 128.5880000 0.0919900 P 3 1.00 49.1158000 0.2599100 20.5035000 0.4288700 8.9871200 0.3269100 P 1 1.00 3.6824900 1.0000000 P 2 1.00 1.5217500 0.5692500 0.5926840 0.2686300 P 1 1.00 0.2023700 1.0000000 P 1 1.00 0.0627600 1.0000000 D 3 1.00 41.4526000 0.0251100 11.5403000 0.1362600 3.8854300 0.3532300 D 1 1.00 1.3238000 1.0000000 D 1 1.00 0.4166800 1.0000000 D 1 1.00 0.1133000 1.0000000 F 3 1.00 4.5150397 0.1737857 1.6191972 0.5973381 0.6933595 0.3929396 **** 得出的能量是-5988.86472748 a.u. 另一個輸入文件: # ub3lyp/Gen pop=nboread scrf=(pcm,solvent=CH2Cl2) Title Card Required 0 5 坐標(biāo)數(shù)據(jù) C F N H O Se S 0 6-311+g(d,p) **** Fe 0 S 6 1.00 257539.0000000 0.0002900 38636.9000000 0.0022600 8891.4400000 0.0115200 2544.0100000 0.0456600 844.7770000 0.1403500 312.5270000 0.3142000 S 2 1.00 125.5930000 0.4087800 53.4987000 0.2116300 S 1 1.00 17.7151000 1.0000000 S 1 1.00 7.3767700 1.0000000 S 1 1.00 2.0184700 1.0000000 S 1 1.00 0.7799350 1.0000000 S 1 1.00 0.1142200 1.0000000 S 1 1.00 0.0418890 1.0000000 P 3 1.00 1678.4000000 0.0024900 396.3920000 0.0201500 128.5880000 0.0919900 P 3 1.00 49.1158000 0.2599100 20.5035000 0.4288700 8.9871200 0.3269100 P 1 1.00 3.6824900 1.0000000 P 2 1.00 1.5217500 0.5692500 0.5926840 0.2686300 P 1 1.00 0.2023700 1.0000000 P 1 1.00 0.0627600 1.0000000 D 3 1.00 41.4526000 0.0251100 11.5403000 0.1362600 3.8854300 0.3532300 D 1 1.00 1.3238000 1.0000000 D 1 1.00 0.4166800 1.0000000 D 1 1.00 0.1133000 1.0000000 F 3 1.00 4.5150397 0.1737857 1.6191972 0.5973381 0.6933595 0.3929396 **** $nbo bndidx $end 得出的能量是-5988.86327012 a.u.,在這個文件中,當(dāng)用scrf=(pcm,solvent=CH2Cl2,read)時輸入報錯,把read去掉就能正常運行。請問這是為什么? |
至尊木蟲 (著名寫手)
我叫MT
|
1.你沒有說你的Gaussian版本。我假設(shè)是Gaussian03 e01 windows版 在scrf中加read和不加read,只要你已經(jīng)指定溶劑而不是generic方式自定義溶劑,那么,只要你后面沒有修改esp espinf 以及rsolv,那么得到的結(jié)果能量就應(yīng)該是一致的,你說得到能量不一致,有待于分析你的log文件來尋找原因,你確定分子結(jié)構(gòu)一致? 2.至于第二個加read報錯,是因為你前面多了一個pop=nboread而隨后你給了$nbo bndidx $end ,這兩個是相對應(yīng)的,然而如果加了scrf(read),Gaussian會先試圖去讀你指定修改的esp espinf以及rsolv等這些信息,此時讀到$nbo bndidx $end 讓Gaussian迷惑而報錯。至于前面不出錯,那只是Gaussian的容錯功能,雖然你給了scrf(read),但你沒有給任何修正的esp espinf以及rsolv等,那么Gaussian認(rèn)為你開了個玩笑,它會容忍你的,它默認(rèn)采用自己保存的CH2Cl2的相應(yīng)數(shù)據(jù)而正常忽略scrf(read)的命令。 |

至尊木蟲 (著名寫手)
我叫MT

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 085602化工求調(diào)劑(331分) +8 | 111@127 2026-03-30 | 8/400 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +9 | hanamiko 2026-03-27 | 9/450 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|