| 查看: 1983 | 回复: 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去掉就能正常运行。请问这是为什么？ |
回复此楼» 猜你喜欢
切片机滑丝了怎么办
已经有1人回复
-
上海交大招收材料/化学方向“申请考核” 学术型博士（请勿回复帖子或站内投条）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有150人回复
-
镁金属电解液不溶，怎么办
已经有1人回复
-
新西兰 奥塔哥大学 化学系 分析传感方向 招收2025年入学的硕博研究生
已经有0人回复
-
新西兰 奥塔哥大学 化学系 分析传感方向 招收CSC公派博士生
已经有0人回复
-
2025申博
已经有4人回复
-
薄膜测TEM怎么制样
已经有0人回复
-
水系锌电测试LSV如何确定参数设置
已经有0人回复
-
反应釜维修
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】溶剂化效应自由能的输出问题
已经有10人回复
- 溶剂效应和优化过渡态的问题
已经有10人回复
- 研究溶剂效应出现虚频问题，请教各位大虾指点
已经有9人回复
- Gaussian对分子构型优化加溶剂化与否对最后结果有影响吗？
已经有8人回复
- Gaussian 溶剂化效应 优化不收敛，大家常用什么办法解决？
已经有7人回复
- 【讨论】溶剂对醇氧化反应的影响
已经有6人回复
- 【求助】如何考察环氧丙烷和氧化苯乙烯的溶剂化效应
已经有13人回复
- 【求助】算溶剂化效应的时候老是出错
已经有3人回复
- 【求助】过渡态能计算溶剂化效应吗
已经有4人回复
- 【求助】tBuOH作为溶剂，在gaussian中如何输入？
已经有4人回复
- 【求助】吸收光谱的溶剂化效应算法
已经有8人回复
- 【求助】溶剂化效应 出错
已经有4人回复
4楼2012-04-07 16:17:13
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
zhangmt
至尊木虫 (著名写手)
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
5楼2012-04-07 22:52:07
5
