计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: linima
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 混合溶剂该选用哪种溶剂模型?比如乙腈:H2O=9:1的比例

[复制链接 Copy URL]

7

帖子

0

威望

25

eV
积分
32

Level 2 能力者

16#
发表于 Post on 2020-11-19 16:13:40 | 只看该作者 Only view this author
本帖最后由 ZnPMg 于 2020-11-19 16:34 编辑
sobereva 发表于 2020-11-15 08:25
显然是0.9和0.1的混合比

各种溶剂的epsinf都在3以下,几十的显然是eps。没说是epsinf的情况一律是eps

Sob先生,
Solvent              : Acetonitrile, Eps=  35.688000 Eps(inf)=   1.806874
Solvent              : Water, Eps=  78.355300 Eps(inf)=   1.777849
V(Acetonitrile:Water)=9:1  
混合溶剂的Eps=35.688*0.9+78.355300*0.1=39.95473;
Eps(inf)=0.9*1.806874+0.1* 1.777849=1.8039715

算出混合溶剂的Eps和Eps(inf)值,如下编辑输入文件就可以了吗?
#p M062X/6-31+G(d,p) Opt  Freq  SCRF=(smd,solvent=generic)  

a

0 1
分子结构

Eps=39.95473
Eps(inf)=1.8039715


运算结果出错了,如下: Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned on.
Using the following non-standard input for PCM:
Eps=39.95473
Epsinf=1.8039715
End of line while reading PCM input.
Error termination via Lnk1e in /home/slgao/CCsoft/g16/l301.exe at Thu Nov 19 01:46:18 2020.
Job cpu time:       0 days  0 hours  0 minutes  7.0 seconds.
Elapsed time:       0 days  0 hours  0 minutes  1.1 seconds.
File lengths (MBytes):  RWF=      6 Int=      0 D2E=      0 Chk=      1 Scr=      1

去掉后solvent=generic后可以计算了。



5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

17#
发表于 Post on 2020-11-20 04:34:45 | 只看该作者 Only view this author
ZnPMg 发表于 2020-11-19 16:13
Sob先生,
Solvent              : Acetonitrile, Eps=  35.688000 Eps(inf)=   1.806874
Solvent      ...

显然要去掉generic,你又没提供非极性部分的参数
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

7

帖子

0

威望

25

eV
积分
32

Level 2 能力者

18#
发表于 Post on 2020-11-20 16:58:10 | 只看该作者 Only view this author
本帖最后由 ZnPMg 于 2020-11-23 14:09 编辑

这个帖子(http://bbs.keinsci.com/thread-12720-1-1.html)说的很详细,就是把它漏掉了,打扰各位老师了,不好意思。

2

帖子

0

威望

37

eV
积分
39

Level 2 能力者

19#
发表于 Post on 2022-5-21 16:23:59 | 只看该作者 Only view this author
ZnPMg 发表于 2020-11-19 14:56
先只用一种溶剂,例如:#p M062X/6-31G(d,p) SCRF=(smd,solvent=Acetonitrile),简单算一下,.log文件中会 ...

我的想计算一个过渡态和反应物的单点能,反应 溶剂是 乙腈:水= 1:1,可以直接 用这个数值吗?

1万

帖子

0

威望

7396

eV
积分
18149

Level 6 (一方通行)

20#
发表于 Post on 2022-5-21 16:36:31 | 只看该作者 Only view this author
yiwind 发表于 2022-5-21 09:23
我的想计算一个过渡态和反应物的单点能,反应 溶剂是 乙腈:水= 1:1,可以直接 用这个数值吗?

要用0.5和0.5代替那个公式里的0.9和0.1,自己重新计算
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

1

帖子

0

威望

11

eV
积分
12

Level 1 能力者

21#
发表于 Post on 2022-9-30 11:27:28 | 只看该作者 Only view this author
本帖最后由 Dreamlovesu 于 2022-9-30 11:29 编辑
sobereva 发表于 2019-11-28 02:26
是明尼苏达
看原文就知道水在SMD模型里是一个特例

老师,那既然水是特例,那请问与水的混合溶剂的其它参数应该取什么呢?我看完整个帖子,还是没有找到应该怎样处理。我也查不到水的参数数据。我也是想用与水的混合溶剂做结构优化和单点能计算。谢谢老师!

5万

帖子

99

威望

5万

eV
积分
112374

管理员

公社社长

22#
发表于 Post on 2022-10-2 04:41:03 | 只看该作者 Only view this author
Dreamlovesu 发表于 2022-9-30 11:27
老师,那既然水是特例,那请问与水的混合溶剂的其它参数应该取什么呢?我看完整个帖子,还是没有找到应该 ...

结构优化本来也不推荐用SMD,下文明确说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327http://bbs.keinsci.com/thread-3345-1-1.html

用默认的IEFPCM,自定义混合溶剂的eps就完了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 22:24 , Processed in 0.325621 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list