| 5 | 1/1 | 返回列表 |
| 查看: 2342 | 回復(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認為你開了個玩笑,它會容忍你的,它默認采用自己保存的CH2Cl2的相應(yīng)數(shù)據(jù)而正常忽略scrf(read)的命令。 |

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 一志愿哈工大 材料與化工 +9 | 洛神哥哥 2026-03-31 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +27 | RichLi_ 2026-03-25 | 31/1550 |
|
|
[考研] 一志愿浙江大學(xué)工科動力工程370,數(shù)一121,專業(yè)課135,現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
[考研] 吉大生物學(xué)326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[考研] 求收留 +6 | 1943443204 2026-03-28 | 6/300 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +5 | 朵朵- 2026-03-26 | 7/350 |
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個好覺_16 2026-03-24 | 4/200 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|