计算化学公社

标题: 计算方法和基组求组 [打印本页]

作者
Author: yanmu 时间: 2017-5-30 11:05
标题: 计算方法和基组求组
最近在做硅烷做还原剂的反应计算，里面还有负离子存在，温度对反应产物影响较大，因此对计算精度要求很高，用什么基组和方法更好一点？谢谢大家

作者
Author: sobereva 时间: 2017-5-30 11:25
谈谈量子化学中基组的选择
http://sobereva.com/336
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272
各种后HF方法精度简单横测
http://sobereva.com/378

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

作者
Author: yanmu 时间: 2017-5-30 11:59

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

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

作者
Author: sobereva 时间: 2017-5-30 14:13

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

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

作者
Author: yanmu 时间: 2017-5-30 15:48

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)计算，老师说先用这个，可以再找个备用的以防最后结果不合适，所以一直想问问有没有更合适的基组和方法。

作者
Author: sobereva 时间: 2017-5-30 16:01

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你说的精度很高，要更精确结果必须提升算单点的级别，提升基组、上双杂化，必然明显增加计算量

作者
Author: yanmu 时间: 2017-5-30 19:51

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

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

作者
Author: yanmu 时间: 2017-5-30 20:48
标题: 几何优化和频率分析出错
用高斯优化醋酸根离子时，优化显示正确，但是频率分析有虚频，调整了好几次结构还是有虚频，怎么解决呢？谢谢

作者
Author: 我本是个娃娃 时间: 2017-5-30 21:00
http://sobereva.com/278

作者
Author: yanmu 时间: 2017-5-30 21:12

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

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

作者
Author: 我本是个娃娃 时间: 2017-5-30 21:18

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

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

作者
Author: yanmu 时间: 2017-5-30 21:35

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

我一直是10的-6

作者
Author: 我本是个娃娃 时间: 2017-5-30 21:35

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

贴出来关键词

作者
Author: yanmu 时间: 2017-5-30 21:36

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

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

作者
Author: 我本是个娃娃 时间: 2017-5-30 21:45
#p opt freq m062x/6-31+g(d)即可

作者
Author: yanmu 时间: 2017-5-30 21:48

我本是个娃娃发表于 2017-5-30 21:45
#p opt freq m062x/6-31+g(d)即可

反应温度50呢

作者
Author: 我本是个娃娃 时间: 2017-5-30 21:50

yanmu 发表于 2017-5-30 21:48
反应温度50呢

优化的时候不必写，算热力学数据的时候才需要

作者
Author: yanmu 时间: 2017-5-30 21:55

我本是个娃娃发表于 2017-5-30 21:50
优化的时候不必写，算热力学数据的时候才需要

你的意思是如果有温度和压强条件改变，计算溶剂化单点能的时候加温度和压强吗？那如果我j现在几何优化时加上了，会有什么错误的结果吗?

作者
Author: 我本是个娃娃 时间: 2017-5-30 21:57

yanmu 发表于 2017-5-30 21:55
你的意思是如果有温度和压强条件改变，计算溶剂化单点能的时候加温度和压强吗？那如果我j现在几何优化时 ...

不会，温度和压强并不会改变势能面的形状，所以除非是热力学数据，否则没必要添加

作者
Author: yanmu 时间: 2017-5-30 22:10

我本是个娃娃发表于 2017-5-30 21:57
不会，温度和压强并不会改变势能面的形状，所以除非是热力学数据，否则没必要添加

%nprocshared=8
%mem=13000MB
# scrf=(solvent=CH3CN,read,smd) temperature=323.15 m062x/6-311++g(d,p)
这是计算溶剂化单点能时的命令。那我直接还是用之前的优化结构再做单点能计算没有错吧？因为时第一次独立做一个反应机理，总觉得挺担心的，毕竟优化结构时就加了温度在里面。谢谢你的耐心解释

作者
Author: yanmu 时间: 2017-5-30 22:32

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

%nprocshared=8
%mem=13000MB
# opt freq pop=nbo temperature=323.15 m062x/6-31+g(d)
我在气相优化时就加了温度，计算单点能也加了温度，这样得到的结构和能量会不会有问题？谢谢老师

作者
Author: liyuanhe211 时间: 2017-5-30 22:35
输入文件里破坏掉所有对称性。

作者
Author: liyuanhe211 时间: 2017-5-31 00:36

yanmu 发表于 2017-5-30 22:32
%nprocshared=8
%mem=13000MB
# opt freq pop=nbo temperature=323.15 m062x/6-31+g(d)

不影响opt的结果，但影响freq的热力学部分的结果。

作者
Author: sobereva 时间: 2017-5-31 05:41

yanmu 发表于 2017-5-30 22:10
%nprocshared=8
%mem=13000MB
# scrf=(solvent=CH3CN,read,smd) temperature=323.15 m062x/6-311++g(d, ...

算单点、振动分析时设temperature、pressure毫无意义，写了也白写，哪怕自己对理论一无所知，自己一对比测试就知道根本不影响结果
只有freq时候写这个才有意义

作者
Author: yanmu 时间: 2017-5-31 08:31

liyuanhe211 发表于 2017-5-31 00:36
不影响opt的结果，但影响freq的热力学部分的结果。

如果我做一个反应温度是50度的机理时，计算其能量，能垒，在
# opt freq=noraman pop=nbo m062x/6-31+g(d) 命令中要加temperature=323.15这个吗

作者
Author: yanmu 时间: 2017-5-31 08:39

sobereva 发表于 2017-5-31 05:41
算单点、振动分析时设temperature、pressure毫无意义，写了也白写，哪怕自己对理论一无所知，自己一对 ...

气相优化求校正因子 # opt freq pop=nbo temperature=323.15 m062x/6-31+g(d)
过渡态 #p opt=(calcfc,noeigen,ts) freq pop=nbo temperature=323.15 m062x/6-31+g(d)
溶剂化单点能 # scrf=(solvent=CH3CN,read,smd) m062x/6-311++g(d,p)
这样写命令对吗？

作者
Author: liyuanhe211 时间: 2017-5-31 09:04

yanmu 发表于 2017-5-31 08:31
如果我做一个反应温度是50度的机理时，计算其能量，能垒，在
# opt freq=noraman pop=nbo m062x/6-31+g( ...

#p opt freq m062x/6-31+g(d) temperature=323

另外即使freq计算时不加temperature命令，只要保留了chk文件，可以快速转换得到任意温度下的热力学量数值，不必因为没写temperature重算。

作者
Author: liyuanhe211 时间: 2017-5-31 09:07
本帖最后由 liyuanhe211 于 2017-5-31 09:13 编辑

yanmu 发表于 2017-5-31 08:39
气相优化求校正因子 # opt freq pop=nbo temperature=323.15 m062x/6-31+g(d)
过渡态 #p opt=(calcfc, ...

热力学校正因子是查出来的，不是求出来的，也不是对某一个物质计算出来的。不明白的话读：
谈谈谐振频率校正因子
http://sobereva.com/221

不用写那“点15”，不论是实验还是计算根本精确不到 0.15°C 带来的差异，非要写个“点15”就像“去买1斤6两零15.28378g白菜一样”令人发笑

气相优化和频率计算 #p opt freq temperature=323 m062x/6-31+g(d)
过渡态优化和频率计算 #p opt=(calcfc,noeigen,ts) freq temperature=323 m062x/6-31+g(d)
溶剂化【自由】能 #p scrf=(solvent=acetonitrile,smd) m052x/6-31G(d) 的单点能与#p m052x/6-31G(d)做差

作者
Author: yanmu 时间: 2017-5-31 15:05

liyuanhe211 发表于 2017-5-31 09:04
#p opt freq m062x/6-31+g(d) temperature=323

另外即使freq计算时不加temperature命令，只要保留了ch ...

好的，谢谢老师

作者
Author: yanmu 时间: 2017-5-31 15:12

yanmu 发表于 2017-5-31 15:05
好的，谢谢老师

谢谢老师的耐心解答。那最后一个求溶剂化单点能时#p scrf=(solvent=acetonitrile,smd) m052x/6-31G(d)为什么变成了m052x/6-31G(d)，我原先用的m062x/6-311++g(d,p)，精度不够吗

作者
Author: yanmu 时间: 2017-5-31 15:13

liyuanhe211 发表于 2017-5-31 09:04
#p opt freq m062x/6-31+g(d) temperature=323

另外即使freq计算时不加temperature命令，只要保留了ch ...

作者
Author: sobereva 时间: 2017-5-31 15:25

yanmu 发表于 2017-5-31 15:13
谢谢老师的耐心解答。那最后一个求溶剂化单点能时#p scrf=(solvent=acetonitrile,smd) m052x/6-31G(d)为 ...

看此文
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327
SMD结合M052X/6-31G*时溶解自由能计算最准确

作者
Author: yanmu 时间: 2017-5-31 15:38

sobereva 发表于 2017-5-31 15:25
看此文
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

谢谢老师，明白了好多

作者
Author: liyuanhe211 时间: 2017-5-31 20:38

yanmu 发表于 2017-5-31 15:12
谢谢老师的耐心解答。那最后一个求溶剂化单点能时#p scrf=(solvent=acetonitrile,smd) m052x/6-31G(d)为 ...

谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

作者
Author: yanmu 时间: 2017-6-1 15:08

liyuanhe211 发表于 2017-5-31 20:38
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

谢谢，现在又懂了好多

作者
Author: yanmu 时间: 2017-6-8 21:34

sobereva 发表于 2017-5-31 15:25
看此文
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

sob老师，我找过度态时出现 Error in internal coordinate system，然后我就加了opt=cartesian,能量一直下降，就是不收敛，需要加什么关键词

作者
Author: sobereva 时间: 2017-6-8 23:23

yanmu 发表于 2017-6-8 21:34
sob老师，我找过度态时出现 Error in internal coordinate system，然后我就加了opt=cartesian,能量一直 ...

量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

作者
Author: yanmu 时间: 2017-6-9 08:19

sobereva 发表于 2017-6-8 23:23
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

谢谢

作者
Author: yanmu 时间: 2017-6-29 09:31

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

老师，我目前用到的溶剂化模型smd结合m062x/6-311++g(d,p)或m052x/6-31g（d），能垒算的特别高，有点不符合实验现象，想问一下有没有更好的计算溶剂化单点能的基组能将能垒降低一些?

作者
Author: sobereva 时间: 2017-6-29 09:43

yanmu 发表于 2017-6-29 09:31
老师，我目前用到的溶剂化模型smd结合m062x/6-311++g(d,p)或m052x/6-31g（d），能垒算的特别高，有点不符 ...

溶解自由能就用m052x/6-31g*，不要考虑别的
单点可以尝试wB97XD、MN15、B3LYP。或者双杂化DSD-PBEP86-D3(BJ)，此时基组起码用def2-TZVP
给氢加弥散无意义，见
谈谈量子化学中基组的选择
http://sobereva.com/336
如果不是非得用弥散的体系，用def2-TZVP比6-311++G**强得多

作者
Author: yanmu 时间: 2017-6-29 11:32

sobereva 发表于 2017-6-29 09:43
溶解自由能就用m052x/6-31g*，不要考虑别的
单点可以尝试wB97XD、MN15、B3LYP。或者双杂化DSD-PBEP86- ...

老师，我几何优化用的就是m062x/6-31g(d)了，所以计算单点能级别更好的基组是哪个呢?

作者
Author: sobereva 时间: 2017-6-29 21:23

yanmu 发表于 2017-6-29 11:32
老师，我几何优化用的就是m062x/6-31g(d)了，所以计算单点能级别更好的基组是哪个呢?

谈谈量子化学中基组的选择
http://sobereva.com/336
def2-TZVP

作者
Author: yanmu 时间: 2017-6-30 08:24

sobereva 发表于 2017-6-29 21:23
谈谈量子化学中基组的选择
http://sobereva.com/336
def2-TZVP

#p scrf=(solvent=acetonitrile,smd) m062x/def2-tzvp 是这样写吗

作者
Author: 冰释之川 时间: 2017-6-30 09:40

yanmu 发表于 2017-6-30 08:24
#p scrf=(solvent=acetonitrile,smd) m062x/def2-tzvp 是这样写吗

RTFM！！ http://sobereva.com/g09/m_basis_sets.htm

#p scrf=(solvent=acetonitrile,smd) m062x/def2tzvp

作者
Author: yanmu 时间: 2017-6-30 10:21

冰释之川发表于 2017-6-30 09:40
RTFM！！ http://sobereva.com/g09/m_basis_sets.htm

#p scrf=(solvent=acetonitrile,smd) m062x/def ...

谢谢

作者
Author: yanmu 时间: 2017-7-4 22:19

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

在帖子中看到这样一段话，“对于Gaussian、Crystal这样基于高斯基函数的程序，算孤立体系常用的6-31G**基组也适合有机周期性体系计算。但很多专门给分子计算的基组不宜直接用在周期性计算上，如6-311G*计算Si。因为这些基组往往有指数较小的呈一定弥散特征的基函数，对于描述原子轨道的尾巴特征有益，但在周期性体系中并不需要描述这点，还可能造成数值问题，故需要将这些基函数的指数人为调大或去掉。”上面这一段的意思是用6-311G*计算Si反而会出问题，不合适吗?因为目前计算反应体系主要就是硅，所以想确定一下。谢谢老师

作者
Author: sobereva 时间: 2017-7-5 06:24

yanmu 发表于 2017-7-4 22:19
在帖子中看到这样一段话，“对于Gaussian、Crystal这样基于高斯基函数的程序，算孤立体系常用的6-31G**基 ...

不用PBC方式来计算的话，6-311G*算Si毫无问题。
对于PBC计算含Si的晶体，也不一定出问题，但是可能遇到一些数值稳定性问题出现诡异结果。

作者
Author: yanmu 时间: 2017-7-5 09:00

sobereva 发表于 2017-7-5 06:24
不用PBC方式来计算的话，6-311G*算Si毫无问题。
对于PBC计算含Si的晶体，也不一定出问题，但是可能遇 ...

好的，谢谢老师

作者
Author: yanmu 时间: 2017-7-12 11:37

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

老师，优化一些对称性高的分子有虚频，于是我用了opt=tight,int=ultrafine,nosymm,虚频消除掉了。但是反应体系中有其他反应底物中间体等并没有用这个，我统计热力学数据，计算势垒时会不会不准确。

作者
Author: sobereva 时间: 2017-7-12 17:43

yanmu 发表于 2017-7-12 11:37
老师，优化一些对称性高的分子有虚频，于是我用了opt=tight,int=ultrafine,nosymm,虚频消除掉了。但是反 ...

不会不准确。
但如果横向对比不同体系能量，一个用了ultrafine则其余都得用

作者
Author: yanmu 时间: 2017-7-12 22:15

sobereva 发表于 2017-7-12 17:43
不会不准确。
但如果横向对比不同体系能量，一个用了ultrafine则其余都得用

反应中涉及到三个反应底物，5个中间体，目的是把整个路径势能图计算出来，但其中一个中间体因为对称性高，用到了ultrafine，消了虚频，其他的都没用这个关键词，怕影响能量计算结果，是不是应该破坏掉其所有对称性，重新优化那个中间体计算更准确呢？谢谢老师

作者
Author: sobereva 时间: 2017-7-13 10:53

yanmu 发表于 2017-7-12 22:15
反应中涉及到三个反应底物，5个中间体，目的是把整个路径势能图计算出来，但其中一个中间体因为对称性高 ...

你可以用优化出没虚频结构，再用默认的fine算一次单点，这样其它结构也就不需要再在ultrafine下计算一遍能量了

作者
Author: yanmu 时间: 2017-7-14 09:17

sobereva 发表于 2017-7-13 10:53
你可以用优化出没虚频结构，再用默认的fine算一次单点，这样其它结构也就不需要再在ultrafine下计算一 ...

用int=ultrafine优化出的没虚频的结构，再计算单点能时也加上int=ultrafine计算一下单点能，这样理解对吗?

作者
Author: sobereva 时间: 2017-7-14 10:12

yanmu 发表于 2017-7-14 09:17
用int=ultrafine优化出的没虚频的结构，再计算单点能时也加上int=ultrafine计算一下单点能，这样理解对吗 ...

优化用不用ultrafine和单点用不用ultrafine是两回事，可以都ultrafine，也可以只有一个用ultrafine

作者
Author: yanmu 时间: 2017-7-14 11:03

sobereva 发表于 2017-7-14 10:12
优化用不用ultrafine和单点用不用ultrafine是两回事，可以都ultrafine，也可以只有一个用ultrafine

那我想解决的由ultrafine优化消虚频带来的问题怎么处理呢？我可能没有理解到老师的意思

作者
Author: sobereva 时间: 2017-7-14 12:05

yanmu 发表于 2017-7-14 11:03
那我想解决的由ultrafine优化消虚频带来的问题怎么处理呢？我可能没有理解到老师的意思

我不知道你具体想解决什么
反正你就记住，算能量的时候，如果要横向对比，一个带ultrafine则所有其它的必须也得带，否则没可比性，误差没法抵消

作者
Author: yanmu 时间: 2017-10-23 12:08

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

看到了一篇咪唑类化合物在乙腈溶液中负氢解离焓的文献，气相结构与能量的计算采用 MP2/6-311 + +G( d，p) / / B3LYP /6-31 + G( d) 方法，溶剂化能的计算采用 HF/6-31 + G( d，p) / /IEFPCM/UA0( 1. 20) 方法，这个是2012年的。现在我想计算一下硅烷的负氢解离焓，用什么方法比较合适呢?x谢谢

作者
Author: sobereva 时间: 2017-10-23 14:18

yanmu 发表于 2017-10-23 12:08
看到了一篇咪唑类化合物在乙腈溶液中负氢解离焓的文献，气相结构与能量的计算采用 MP2/6-311 + +G( d，p) ...

溶剂化能按这里面用M052X/6-31G*结合SMD
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

优化用M062X/6-311G*，单点用M062X/def2-TZVP（要求更高可以用双杂化泛函）
我不晓得你说的负氢解离是个什么过程，如果计算时牵扯到阴离子，那么单点计算时建议用ma-TZVP

作者
Author: yanmu 时间: 2017-10-23 20:58

sobereva 发表于 2017-10-23 14:18
溶剂化能按这里面用M052X/6-31G*结合SMD
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobe ...

这个体系主要就是一个供氢給体分子，解离成一个阳离子和一个H负离子，需要用到气相中的焓和溶剂化能数据。一个物质的溶剂化能等于溶剂模型化单点能减去气相单点能，那么按照老师所说的，结构优化用M062X/6-311G*，溶剂模型化单点能用M052X/6-31G*结合SMD，气相单点能用M062X/def2-TZVP，我这样理解对吗?
溶剂型单点能#p scrf=(solvent=acetonitrile,read,smd) ，气相单点能以前没算过，关键词直接写#p scrf 吗？谢谢老师

作者
Author: sobereva 时间: 2017-10-24 00:04

yanmu 发表于 2017-10-23 20:58
这个体系主要就是一个供氢給体分子，解离成一个阳离子和一个H负离子，需要用到气相中的焓和溶剂化能数据 ...

对
用于算溶解自由能的气相单点能直接写# M052X/6-31G*

作者
Author: yanmu 时间: 2017-10-24 12:25

sobereva 发表于 2017-10-24 00:04
对
用于算溶解自由能的气相单点能直接写# M052X/6-31G*

溶剂模型化单点能#p scrf=(solvent=acetonitrile,read,smd) M052X/6-31G*，气相单点能# M052X/6-31G*。
那溶剂相和气相单点能的计算都是用M052X/6-31G*，溶剂化能就等于两者相减。那我看的那篇文献中原文描述是---气相结构与能量的计算采用 MP2/6-311 + +G( d，p) / / B3LYP /6-31 + G( d) 方法，溶剂化能的计算采用 HF/6-31 + G( d，p) / /IEFPCM/UA0( 1. 20) 方法，这两种单点能的计算方法是不同的。我都用M052X/6-31G*是可以的吗

作者
Author: sobereva 时间: 2017-10-24 15:41

yanmu 发表于 2017-10-24 12:25
溶剂模型化单点能#p scrf=(solvent=acetonitrile,read,smd) M052X/6-31G*，气相单点能# M052X/6-31G*。
...

甭管你之前看的文章了，结果肯定不如SMD好。那个计算气相能量的级别在现在来看性价比也极低

作者
Author: yanmu 时间: 2017-10-24 21:39

sobereva 发表于 2017-10-24 15:41
甭管你之前看的文章了，结果肯定不如SMD好。那个计算气相能量的级别在现在来看性价比也极低

主要是因为他的气相和溶剂型单点用的不同基组，所以我犹豫了一下。谢谢老师的耐心解答

欢迎光临计算化学公社 (http://bbs.keinsci.com/)