计算化学公社

标题: 虚频影响ZPE的结果吗？ [打印本页]

作者
Author: wxhwbh 时间: 2020-2-12 21:25
标题: 虚频影响ZPE的结果吗？
我使用B3LYP/aug-cc-pvtz 计算得到某个构型，想计算它的热力学量。计划使用B3LYP/6-31G(d)结合ZPE校正因子得到ZPE。可是算出来的结果有6个虚频（原来的计算级别是没有虚频的），想问问虚频的存在对ZPE的计算有影响吗？

作者
Author: zjxitcc 时间: 2020-2-12 22:19
结构优化与freq必须在同一级别下进行。就算没有虚频，你也是不能用的，因为你切换了基组，力不为零，算热力学量没有意义。

作者
Author: wxhwbh 时间: 2020-2-16 16:29

zjxitcc 发表于 2020-2-12 22:19
结构优化与freq必须在同一级别下进行。就算没有虚频，你也是不能用的，因为你切换了基组，力不为零，算热力 ...

谢谢您。我主要是想把频率校正因子考虑进去。
比如说我现在已经在较大的基组下得到了体系的最优结构，想计算它的热力学量。而这个基组的频率校正因子查不到，就只能使用已有校正因子的级别（比如BELYP/6-31G（d))计算。如果这样的话，是不是要做BELYP/6-31G（d)下的opt freq呢？但这样得到的构型必定不是原来的构型，请问这样的计算结果有意义吗？

作者
Author: niobium 时间: 2020-2-16 18:13

wxhwbh 发表于 2020-2-16 16:29
谢谢您。我主要是想把频率校正因子考虑进去。
比如说我现在已经在较大的基组下得到了体系的最优结构，想 ...

不妨把你的计算水平写出来，你找不到校正因子，说不定别人有

作者
Author: niobium 时间: 2020-2-16 18:15

wxhwbh 发表于 2020-2-16 16:29
谢谢您。我主要是想把频率校正因子考虑进去。
比如说我现在已经在较大的基组下得到了体系的最优结构，想 ...

哦，主贴里面有

作者
Author: niobium 时间: 2020-2-16 18:17

wxhwbh 发表于 2020-2-16 16:29
谢谢您。我主要是想把频率校正因子考虑进去。
比如说我现在已经在较大的基组下得到了体系的最优结构，想 ...

http://bbs.keinsci.com/thread-3805-1-1.html

作者
Author: wxhwbh 时间: 2020-2-16 19:47

niobium 发表于 2020-2-16 18:17
http://bbs.keinsci.com/thread-3805-1-1.html

我的计算级别是 B3LYP-D3/aug-cc-pVTZ (aug-cc-pV(T+d)Z for sulfur)。这种混合基组不可能有吧[捂脸]

作者
Author: sobereva 时间: 2020-2-18 12:35

wxhwbh 发表于 2020-2-16 16:29
谢谢您。我主要是想把频率校正因子考虑进去。
比如说我现在已经在较大的基组下得到了体系的最优结构，想 ...

可以自己拟合，看
自行拟合基频校正因子与ZPE校正因子的简单方法
http://sobereva.com/391（http://bbs.keinsci.com/thread-6907-1-1.html）

PS：几个月内我会发表一篇几乎所有主流方法结合几乎所有常见基组的频率校正因子数据库

作者
Author: wxhwbh 时间: 2020-2-18 13:03

sobereva 发表于 2020-2-18 12:35
可以自己拟合，看
自行拟合基频校正因子与ZPE校正因子的简单方法
http://sobereva.com/391（http://bbs ...

谢谢社长，顺便问一下，像这种拟合的数据集在哪里可以找到呢？因为我这个基组有S原子，而帖子里的数据集都没有S。

作者
Author: sobereva 时间: 2020-2-19 13:03

wxhwbh 发表于 2020-2-18 13:03
谢谢社长，顺便问一下，像这种拟合的数据集在哪里可以找到呢？因为我这个基组有S原子，而帖子里的数据集 ...

不需要专门考虑某种元素。频率校正因子的思想不考虑元素的依赖性

作者
Author: wxhwbh 时间: 2020-2-19 19:27

sobereva 发表于 2020-2-19 13:03
不需要专门考虑某种元素。频率校正因子的思想不考虑元素的依赖性

可是我的基组定义的对S原子用aug-cc-pV(T+d)Z ，其他元素用aug-cc-pVTZ。这样必须得用有S的分子做测试吧。
或者我采用另一种思路：直接取大部分元素用的基组aug-cc-pVTZ。不过我的方法是B3LYP-D3，而只找得到不加D3校正的频率校正因子。想问问B3LYP/aug-cc-pVTZ的校正因子用在B3LYP-D3/aug-cc-pVTZ的体系上合适吗？

作者
Author: sobereva 时间: 2020-2-21 19:40

wxhwbh 发表于 2020-2-19 19:27
可是我的基组定义的对S原子用aug-cc-pV(T+d)Z ，其他元素用aug-cc-pVTZ。这样必须得用有S的分子做测试吧 ...

没必要那么讲究，是否+d的影响完全可忽略不计
合适

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