计算化学公社

标题: 求助-算单态自由基的反应过程方法基组比较 [打印本页]

作者
Author: 风飞 时间: 2019-11-29 16:36
标题: 求助-算单态自由基的反应过程方法基组比较
本帖最后由风飞于 2019-12-19 16:40 编辑

各位老师，大家好，
我利用def2TZVPP empiricaldispersion=gd3 m062x 计算单态自由基的反应过程中（过渡态，中间体，反应能垒和速率常数）：审稿人发现在另外类似体系文章中用的方法和基组是：M06–2X/augcc-
pVTZ计算进行优化，并且使用CCSD（T）计算单点，审稿人问我这种方法和之前的方法相比有什么优势，请问这个该怎么回答合适呢？是否可以说，在精度满足的情况下，计算速度快呢？，对于类似体系，我也有看到用B3LYP-D3/6-311+G(2df,2p)进行优化，然后用MP2计算单点

谢谢！

作者
Author: liyuanhe211 时间: 2019-11-29 17:36
精度上没啥优势，M06-2X/aug-cc-pVTZ的唯一缺点是太浪费了。只能说虽然精度没有它高，但是也够用。如果愿意补的话，可以CCSD(T)或DLPNO-CCSD(T)补个单点，这样就可以说精度差不多，你的级别还更匹配、不浪费。

CCSD(T)还是CCSDT写清楚。
不要用MP2做这个目的。

作者
Author: sobereva 时间: 2019-11-30 00:43
又是审稿人不懂还瞎找茬

你当前的级别完全没问题，提升到aug-cc-pVTZ优化纯粹是瞎浪费时间，几乎没有任何可查觉的改进。本来几何优化就对基组很不敏感，甚至def-TZVP都已经完全够了。

看此文
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387（http://bbs.keinsci.com/thread-6600-1-1.html）

如果审稿人就是说的CCSDT，而不是你把CCSD(T)写成了CCSDT，说明审稿人连基本理论计算常识都没有。CCSDT多数情况下结果反倒还不如CCSD(T)，而耗时多一个数量级以上，Gaussian、ORCA等大多数量化程序还不支持。看
各种后HF方法精度简单横测
http://sobereva.com/378（http://bbs.keinsci.com/thread-6005-1-1.html）

6-311+G(2df,2p)用于优化是非常烂的选择，本身这个档次的基组就没有存在意义（远不如用ma-TZVP、def2-TZVP），普通泛函做优化根本没必要f极化。
用MP2算单点简直是有病，干出这种事说明作者连MP2和DFT的相对精度关系好坏都不懂，仔细看此文的相关内容了解常识
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

作者
Author: 风飞 时间: 2019-12-1 19:30

sobereva 发表于 2019-11-30 00:43
又是审稿人不懂还瞎找茬

你当前的级别完全没问题，提升到aug-cc-pVTZ优化纯粹是瞎浪费时间，几乎没有任 ...

好的，谢谢老师

作者
Author: 风飞 时间: 2019-12-1 19:33

liyuanhe211 发表于 2019-11-29 17:36
精度上没啥优势，M06-2X/aug-cc-pVTZ的唯一缺点是太浪费了。只能说虽然精度没有它高，但是也够用。如果愿意 ...

谢谢老师，有两个审稿人，有个说是CCSDT，有个说是CCSDT, 那我这边在用CCSD（T）做个单点，

作者
Author: 风飞 时间: 2019-12-19 16:29

liyuanhe211 发表于 2019-11-29 17:36
精度上没啥优势，M06-2X/aug-cc-pVTZ的唯一缺点是太浪费了。只能说虽然精度没有它高，但是也够用。如果愿意 ...

老师，您好，我的结构优化时在m062x-D3/def2TZVPP 下进行优化的，根据您说的用CCSD(T) 计算单点，那么采用的基组是def2TZVPP，还是aug-cc-pVTZ？我尝试了直接写CCSD(T) 不写基组情况下得到的能量和CCSD(T) /def2TZVPP得到的能量以及CCSD(T) /aug-cc-pVTZ得到的能量都不一样。

作者
Author: liyuanhe211 时间: 2019-12-19 17:30

风飞发表于 2019-12-19 16:29
老师，您好，我的结构优化时在m062x-D3/def2TZVPP 下进行优化的，根据您说的用CCSD(T) 计算单点，那么采 ...

采用你能算的、适合CCSD(T)的基组，不必追求和优化一致。
当然得写基组，不写基组默认是STO-3G

作者
Author: liyuanhe211 时间: 2019-12-19 17:31

风飞发表于 2019-12-1 19:33
谢谢老师，有两个审稿人，有个说是CCSDT，有个说是CCSDT, 那我这边在用CCSD（T）做个单点，

"有个说是CCSDT，有个说是CCSDT"没看懂

作者
Author: 风飞 时间: 2019-12-19 17:52

liyuanhe211 发表于 2019-12-19 17:30
采用你能算的、适合CCSD(T)的基组，不必追求和优化一致。
当然得写基组，不写基组默认是STO-3G

嗯，谢谢老师，我没有采用CCSD(T，T1diag)/aug-cc-pVTZ 计算单点，而采用的是CCSD(T，T1diag)/def2TZVPP 计算的单点，这样可以吗？因为我之前在采用CCSD(T，T1diag)/def2TZVPP 计算时，每个基本都要花上>16h,用CCSD(T，T1diag)/aug-cc-pVTZ 计算单点时间会更长,

作者
Author: sobereva 时间: 2019-12-19 18:44

风飞发表于 2019-12-19 17:52
嗯，谢谢老师，我没有采用CCSD(T，T1diag)/aug-cc-pVTZ 计算单点，而采用的是CCSD(T，T1diag)/def2TZVPP ...

def2-TZVPP尚且要花那么长，aug-cc-pVTZ你根本不可能算得动
def2-TZVPP近似等于cc-pVTZ，如果弥散函数不是必需的情况，结果当然可以用

作者
Author: 风飞 时间: 2019-12-19 20:48

liyuanhe211 发表于 2019-12-19 17:30
采用你能算的、适合CCSD(T)的基组，不必追求和优化一致。
当然得写基组，不写基组默认是STO-3G

好的谢谢老师

作者
Author: 风飞 时间: 2019-12-19 20:48

sobereva 发表于 2019-12-19 18:44
def2-TZVPP尚且要花那么长，aug-cc-pVTZ你根本不可能算得动
def2-TZVPP近似等于cc-pVTZ，如果弥散函数不 ...

好的谢谢老师

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