计算方法和基组求组

yanmu

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

反应温度50呢

我本是个娃娃

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

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

yanmu

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

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

我本是个娃娃

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

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

yanmu

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

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

yanmu

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)
我在气相优化时就加了温度，计算单点能也加了温度，这样得到的结构和能量会不会有问题？谢谢老师

liyuanhe211

输入文件里破坏掉所有对称性。

liyuanhe211

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

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

sobereva

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时候写这个才有意义

yanmu

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

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

yanmu

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)
这样写命令对吗？

liyuanhe211

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重算。

liyuanhe211

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)做差

yanmu

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

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

好的，谢谢老师

yanmu

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)，精度不够吗