计算方法和基组求组

yanmu · 发表于 Post on 2017-5-30 11:05:48

最近在做硅烷做还原剂的反应计算，里面还有负离子存在，温度对反应产物影响较大，因此对计算精度要求很高，用什么基组和方法更好一点？谢谢大家

sobereva · 发表于 Post on 2017-5-30 11:25:24

谈谈量子化学中基组的选择
http://sobereva.com/336
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272
各种后HF方法精度简单横测
http://sobereva.com/378

该怎么选决定因素很多。要考虑哪些因素第一个博文里说了
总不能一说高精度就让你用CCSD(T)/aug-cc-pVQZ

yanmu · 发表于 Post on 2017-5-30 11:59:04

sobereva 发表于 2017-5-30 11:25
谈谈量子化学中基组的选择
http://sobereva.com/336
简谈量子化学计算中DFT泛函的选择

反应物是二苯硅烷和二氧化碳和苯基胺，温度50，极性溶剂。看了一篇计算文献用的是M052X/6-31g(d)||m052x/6-311++g(d,p),我目前要计算的和他差别在于使用催化剂醋酸根离子亲核进攻硅烷，促进还原。那我继续用他的方法吗？还想问一下，m062x和m052x有什么差别，用哪一个更好一点？真心感谢

sobereva · 发表于 Post on 2017-5-30 14:13:46

yanmu 发表于 2017-5-30 11:59
反应物是二苯硅烷和二氧化碳和苯基胺，温度50，极性溶剂。看了一篇计算文献用的是M052X/6-31g(d)||m052x/ ...

m062x，整体精度比M052X好得多
不要效仿文献，很多都是菜鸟写的，你看的那篇用的级别就不理想。而且根本没有||这种写法，写正确了那也是m052x/6-311++g(d,p)//M052X/6-31g(d)

yanmu · 发表于 Post on 2017-5-30 15:48:55

sobereva 发表于 2017-5-30 14:13
m062x，整体精度比M052X好得多
不要效仿文献，很多都是菜鸟写的，你看的那篇用的级别就不理想。而且根 ...

我刚开始想用b3lyp/6-31+g(d)气相优化，m062x/6-311++g(d,p)溶剂化单点能，但b3lyp优化出的醋酸根离子和硅烷不成键，有弱相互作用；后来用m062x/6-31+g(d)气相优化结果成键，想用m062x/6-311+g(d,p)//m062x/6-31+g(d)计算，老师说先用这个，可以再找个备用的以防最后结果不合适，所以一直想问问有没有更合适的基组和方法。

sobereva · 发表于 Post on 2017-5-30 16:01:14

yanmu 发表于 2017-5-30 15:48
我刚开始想用b3lyp/6-31+g(d)气相优化，m062x/6-311++g(d,p)溶剂化单点能，但b3lyp优化出的醋酸根离子和 ...

可以用m062x/6-311+g(d,p)//m062x/6-31+g(d)这个计算
结果不会定性错误，但算不上1L你说的精度很高，要更精确结果必须提升算单点的级别，提升基组、上双杂化，必然明显增加计算量

yanmu · 发表于 Post on 2017-5-30 19:51:06

sobereva 发表于 2017-5-30 16:01
可以用m062x/6-311+g(d,p)//m062x/6-31+g(d)这个计算
结果不会定性错误，但算不上1L你说的精度很高， ...

更高级别的算单点能的基组可以用哪个呢？对于一些带好几个苯环的结构，优化结构特别慢，计算时间长，显示收敛了但一直不结束。

yanmu · 发表于 Post on 2017-5-30 20:48:35

用高斯优化醋酸根离子时，优化显示正确，但是频率分析有虚频，调整了好几次结构还是有虚频，怎么解决呢？谢谢

我本是个娃娃 · 发表于 Post on 2017-5-30 21:00:11

http://sobereva.com/278

yanmu · 发表于 Post on 2017-5-30 21:12:04

我本是个娃娃发表于 2017-5-30 21:00
http://sobereva.com/278

谢谢。还想问一个问题，就是我几何优化时显示收敛了，到10的-6次方了，但是一直不结束下来，是分子对称性太高了吗？

我本是个娃娃 · 发表于 Post on 2017-5-30 21:18:37

yanmu 发表于 2017-5-30 21:12
谢谢。还想问一个问题，就是我几何优化时显示收敛了，到10的-6次方了，但是一直不结束下来，是分子对称性 ...

10E-6?不应该是10E-8吗？

yanmu · 发表于 Post on 2017-5-30 21:35:17

我本是个娃娃发表于 2017-5-30 21:18
10E-6?不应该是10E-8吗？

我一直是10的-6

我本是个娃娃 · 发表于 Post on 2017-5-30 21:35:41

yanmu 发表于 2017-5-30 21:35
我一直是10的-6

贴出来关键词

yanmu · 发表于 Post on 2017-5-30 21:36:32

我本是个娃娃发表于 2017-5-30 21:35
贴出来关键词

%nprocshared=8
%mem=13000MB
# opt freq=noraman pop=nbo temperature=323.15 m062x/6-31+g(d)

我本是个娃娃 · 发表于 Post on 2017-5-30 21:45:22

#p opt freq m062x/6-31+g(d)即可

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 计算方法和基组求组

几何优化和频率分析出错

浏览过的版块