| 5 | 1/1 | 返回列表 |
| 查看: 2339 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
溶劑化效應(yīng)輸入
|
||
|
請(qǐng)問(wèn)溶劑化效應(yīng)的輸入命令scrf=(pcm,solvent=CH2Cl2,read)和scrf=(pcm,solvent=CH2Cl2)有什么區(qū)別? 在優(yōu)化好的構(gòu)型的基礎(chǔ)上用大基組進(jìn)行計(jì)算,分別用了上面兩個(gè)命令,發(fā)現(xiàn)得出的能量,以及相關(guān)的自旋密度等數(shù)據(jù)不一樣,請(qǐng)問(wèn)這是怎么回事? 輸入文件如下: # 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. 另一個(gè)輸入文件: # 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.,在這個(gè)文件中,當(dāng)用scrf=(pcm,solvent=CH2Cl2,read)時(shí)輸入報(bào)錯(cuò),把read去掉就能正常運(yùn)行。請(qǐng)問(wèn)這是為什么? |
至尊木蟲(chóng) (著名寫手)
我叫MT

至尊木蟲(chóng) (著名寫手)
我叫MT
|
1.你沒(méi)有說(shuō)你的Gaussian版本。我假設(shè)是Gaussian03 e01 windows版 在scrf中加read和不加read,只要你已經(jīng)指定溶劑而不是generic方式自定義溶劑,那么,只要你后面沒(méi)有修改esp espinf 以及rsolv,那么得到的結(jié)果能量就應(yīng)該是一致的,你說(shuō)得到能量不一致,有待于分析你的log文件來(lái)尋找原因,你確定分子結(jié)構(gòu)一致? 2.至于第二個(gè)加read報(bào)錯(cuò),是因?yàn)槟闱懊娑嗔艘粋(gè)pop=nboread而隨后你給了$nbo bndidx $end ,這兩個(gè)是相對(duì)應(yīng)的,然而如果加了scrf(read),Gaussian會(huì)先試圖去讀你指定修改的esp espinf以及rsolv等這些信息,此時(shí)讀到$nbo bndidx $end 讓Gaussian迷惑而報(bào)錯(cuò)。至于前面不出錯(cuò),那只是Gaussian的容錯(cuò)功能,雖然你給了scrf(read),但你沒(méi)有給任何修正的esp espinf以及rsolv等,那么Gaussian認(rèn)為你開(kāi)了個(gè)玩笑,它會(huì)容忍你的,它默認(rèn)采用自己保存的CH2Cl2的相應(yīng)數(shù)據(jù)而正常忽略scrf(read)的命令。 |

|
謝謝您的回復(fù)。也就是說(shuō),既然我已經(jīng)指定了溶劑solvent=CH2Cl2,我加和不加read是一樣的(反正我沒(méi)在最后一行自定義溶劑的一些參數(shù))。 那看來(lái)導(dǎo)致結(jié)果不一樣的原因是因?yàn)榈谝粋(gè)計(jì)算是用linux版的G03,而第二個(gè)計(jì)算是用linux版的G09。高斯手冊(cè)上也說(shuō)“PCM 形式的某些細(xì)節(jié)與執(zhí)行方式已不同于 Gaussian 03”。再請(qǐng)教您,這個(gè)應(yīng)該怎么理解呢?我要選擇能量的話,G03的合理呢?還是G09的合理呢? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[有機(jī)交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-03-29 | 8/400 |
|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 6/300 |
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué),化工,生物,生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入! +3 | llss0711 2026-03-28 | 6/300 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
[考研] 一志愿上海理工能源動(dòng)力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛(ài)吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語(yǔ)已過(guò)四六級(jí),六級(jí)521,求調(diào)劑 +4 | 美味可樂(lè)雞翅 2026-03-26 | 4/200 |
|
|
[考研] 打過(guò)很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|