| 5 | 1/1 | 返回列表 |
| 查看: 2340 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
溶劑化效應(yīng)輸入
|
||
|
請問溶劑化效應(yīng)的輸入命令scrf=(pcm,solvent=CH2Cl2,read)和scrf=(pcm,solvent=CH2Cl2)有什么區(qū)別? 在優(yōu)化好的構(gòu)型的基礎(chǔ)上用大基組進(jìn)行計(jì)算,分別用了上面兩個命令,發(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)時(shí)輸入報(bào)錯,把read去掉就能正常運(yùn)行。請問這是為什么? |
金蟲 (小有名氣)
至尊木蟲 (著名寫手)
我叫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報(bào)錯,是因?yàn)槟闱懊娑嗔艘粋pop=nboread而隨后你給了$nbo bndidx $end ,這兩個是相對應(yīng)的,然而如果加了scrf(read),Gaussian會先試圖去讀你指定修改的esp espinf以及rsolv等這些信息,此時(shí)讀到$nbo bndidx $end 讓Gaussian迷惑而報(bào)錯。至于前面不出錯,那只是Gaussian的容錯功能,雖然你給了scrf(read),但你沒有給任何修正的esp espinf以及rsolv等,那么Gaussian認(rèn)為你開了個玩笑,它會容忍你的,它默認(rèn)采用自己保存的CH2Cl2的相應(yīng)數(shù)據(jù)而正常忽略scrf(read)的命令。 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
Gymno 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
[考研] 一志愿廈門大學(xué)材料工程專碩354找調(diào)劑。。 +5 | 貝唄鋇鋇 2026-03-30 | 5/250 |
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|