| 查看: 2347 | 回复: 6 | ||
[求助]
溶剂化效应输入
|
|
请问溶剂化效应的输入命令scrf=(pcm,solvent=CH2Cl2,read)和scrf=(pcm,solvent=CH2Cl2)有什么区别？ 在优化好的构型的基础上用大基组进行计算，分别用了上面两个命令，发现得出的能量，以及相关的自旋密度等数据不一样，请问这是怎么回事？ 输入文件如下： # ub3lyp/gen scrf=(pcm,solvent=CH2Cl2,read) pop=full gfoldprint Title Card Required 0 5 坐标数据 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 坐标数据 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.，在这个文件中，当用scrf=(pcm,solvent=CH2Cl2,read)时输入报错，把read去掉就能正常运行。请问这是为什么？ |
回复此楼» 猜你喜欢
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有17人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有289人回复
-
070300总分298求调剂
已经有2人回复
-
一志愿厦门大学化学学硕307求调剂
已经有13人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有1人回复
-
调剂招生
已经有7人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
硕士研究生调剂
已经有0人回复
-
投稿期刊求助
已经有1人回复
-
江西师大化学工程学院调剂今晚0点开放！ 缺额100个，速冲！
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】溶剂化效应自由能的输出问题
已经有10人回复
- 溶剂效应和优化过渡态的问题
已经有10人回复
- 研究溶剂效应出现虚频问题，请教各位大虾指点
已经有9人回复
- Gaussian对分子构型优化加溶剂化与否对最后结果有影响吗？
已经有8人回复
- Gaussian 溶剂化效应 优化不收敛，大家常用什么办法解决？
已经有7人回复
- 【讨论】溶剂对醇氧化反应的影响
已经有6人回复
- 【求助】如何考察环氧丙烷和氧化苯乙烯的溶剂化效应
已经有13人回复
- 【求助】过渡态能计算溶剂化效应吗
已经有4人回复
- 【求助】tBuOH作为溶剂，在gaussian中如何输入？
已经有4人回复
- 【求助】吸收光谱的溶剂化效应算法
已经有8人回复
- 【求助】溶剂化效应 出错
已经有4人回复
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2012-04-06 15:37:07
雨萍: 金币+20, ★★★★★最佳答案 2012-04-06 17:00:27
感谢参与，应助指数 +1
gmy1990: 金币+3 2012-04-06 15:37:07
雨萍: 金币+20, ★★★★★最佳答案 2012-04-06 17:00:27
|
1.你没有说你的Gaussian版本。我假设是Gaussian03 e01 windows版 在scrf中加read和不加read，只要你已经指定溶剂而不是generic方式自定义溶剂，那么，只要你后面没有修改esp espinf 以及rsolv，那么得到的结果能量就应该是一致的，你说得到能量不一致，有待于分析你的log文件来寻找原因，你确定分子结构一致？ 2.至于第二个加read报错，是因为你前面多了一个pop=nboread而随后你给了$nbo bndidx $end ，这两个是相对应的，然而如果加了scrf(read)，Gaussian会先试图去读你指定修改的esp espinf以及rsolv等这些信息，此时读到$nbo bndidx $end 让Gaussian迷惑而报错。至于前面不出错，那只是Gaussian的容错功能，虽然你给了scrf(read),但你没有给任何修正的esp espinf以及rsolv等，那么Gaussian认为你开了个玩笑，它会容忍你的，它默认采用自己保存的CH2Cl2的相应数据而正常忽略scrf(read)的命令。 |
2楼2012-04-06 15:16:33
3楼2012-04-06 17:10:42
4楼2012-04-07 16:17:13
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
5楼2012-04-07 22:52:07
6楼2012-04-08 11:58:11
wenxianliu
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 938.5
- 散金: 25
- 红花: 3
- 帖子: 179
- 在线: 41.5小时
- 虫号: 1127140
- 注册: 2010-10-20
- 性别: GG
- 专业: 理论和计算化学
7楼2012-09-22 15:41:40
