计算化学公社

标题: 请教一个审稿意见 [打印本页]

作者
Author: captain 时间: 2015-8-27 16:07
标题: 请教一个审稿意见
我使用ωB97X-D/6-31G*方法和基组对一个超分子体系进行了结构优化，并在优化结构基础上使用ωB97X-D/6-31++G**方法和基组对单点能进行了计算，
其中有个审稿人有这么两条意见，不清楚具体怎么修改，请大家帮忙看看
1、建议作者对结构优化中的6-31G*基组的收敛性进行确认并报道相关数据。6-31G*基组尺寸较小，且不包含H原子的高角动量基函数，考虑到本文体系涉及氢键，H原子的极化基组可能对结构有显著的影响。
2、同理建议作者确认并报道单点计算中6-31++G**基组的收敛性。考虑到本文需要比较较小（约1kcal/mol左右）的能量差值，double-zeta基组，尤其是Pople style的double-zeta基组是否足够是一个值得确认的问题。

作者
Author: sobereva 时间: 2015-8-27 16:48
对于和氢键有关的氢，建议改用6-31G**再优化，其余的氢的基组不变。

1kcal/mol用这样的基组太糙了。实际上氢的弥散函数几乎没必要还费时间，明显不如用6-311+G**，耗时也不会增加。建议你用ORCA，开RI估计上BLYP-D3/def2QZVP(-g,-f)都没问题。

作者
Author: captain 时间: 2015-8-27 20:21

sobereva 发表于 2015-8-27 16:48
对于和氢键有关的氢，建议改用6-31G**再优化，其余的氢的基组不变。

1kcal/mol用这样的基组太糙了。实际 ...

那请问收敛性具体怎么确认啊？
用您所说的方法重新优化将结果再与我的方法的结果做比较吗？

作者
Author: yflchx 时间: 2015-8-27 20:48
比如DZ、TZ、QZ、5Z、6Z的结果分别为5、4、3.9、3.89、3.889 kcal/mol，我想用TZ就差不多了。

选个体系做个测试。

作者
Author: sobereva 时间: 2015-8-27 22:24

captain 发表于 2015-8-27 20:21
那请问收敛性具体怎么确认啊？
用您所说的方法重新优化将结果再与我的方法的结果做比较吗？

其实都不用比较也不用收敛性测试，一次性用个像样的能让审稿人接受的基组就完事了。6-31++G**太烂了。wB97XD结合6-311+G**并且用counterpoise考虑BSSE校正就不错。如果审稿人还不满意，就用ORCA跑PWPB95-D3/def2-QZVPP，怎么都够了。

作者
Author: captain 时间: 2015-8-27 22:29

sobereva 发表于 2015-8-27 22:24
其实都不用比较也不用收敛性测试，一次性用个像样的能让审稿人接受的基组就完事了。6-31++G**太烂了。wB9 ...

明白了十分感谢！

作者
Author: captain 时间: 2015-8-27 22:30

yflchx 发表于 2015-8-27 20:48
比如DZ、TZ、QZ、5Z、6Z的结果分别为5、4、3.9、3.89、3.889 kcal/mol，我想用TZ就差不多了。

选个体系 ...

十分感谢！

作者
Author: 左撇子wy 时间: 2021-2-16 11:04
我是刚接触这方面的小白，想问大神一个问题ωb97xd和ωb97x-d有区别吗？

作者
Author: sobereva 时间: 2021-2-16 11:22

左撇子wy 发表于 2021-2-16 11:04
我是刚接触这方面的小白，想问大神一个问题ωb97xd和ωb97x-d有区别吗？

本来就是同一个东西

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