计算化学公社

标题: 溶剂化效应的计算请教 [打印本页]

作者
Author: sun0615 时间: 2021-2-11 13:21
标题: 溶剂化效应的计算请教

我参考http://sobereva.com/327?replyto=159里面的计算例子，但是结构改用SMD模型下的优化结构，详细能量计算如下

1. 使用M06-2X/6-31G*优化结构与计算热校正能量，得到溶剂化下的电子能量(a)与自由能热校正量(b)

2. 使用(1)之结构，以M06-2X (D3) / 6-31G*方法计算真空下的单点能量，得到真空下的电子能量(c)

3. 使用(1)之结构，以PWPB95 (D3) / def2-QZVPP方法计算真空下的单点能量，得到真空下的高精确电子能量(d)

最后的总能为(d) + (b) + (c-a) +1.89/627.5095

不知道这样是否正确?

作者
Author: wzkchem5 时间: 2021-2-11 13:29
(c)不应该加D3，而且(c-a)应该是(a-c)。其他没什么问题

作者
Author: sun0615 时间: 2021-2-11 13:48
本帖最后由 sun0615 于 2021-2-11 13:53 编辑

wzkchem5 发表于 2021-2-11 13:29
(c)不应该加D3，而且(c-a)应该是(a-c)。其他没什么问题

谢谢您的回复，是(a-c)，上面笔误

至于计算溶解自由能时是否应该加D3是否有比较详细的解释，或是例子可以参考?

上面有写错的地方，本来「1.使用M06-2X/6-31G*优化结构与计算热校正能量」是要写M06-2X (D3)/6-31G*

因此如果(c)不应该加D3，那应该是结构优化加D3，之后另外计算不加D3的溶剂化单点能量-真空单点能量，或是结构优化就不加D3?

作者
Author: wzkchem5 时间: 2021-2-11 16:46

sun0615 发表于 2021-2-11 13:48
谢谢您的回复，是(a-c)，上面笔误

至于计算溶解自由能时是否应该加D3是否有比较详细的解释，或是例子 ...

同一个结构下的D3色散校正和是否加溶剂无关，所以a和c的D3贡献刚好抵消，加了等于白加

作者
Author: sun0615 时间: 2021-2-11 23:24

wzkchem5 发表于 2021-2-11 16:46
同一个结构下的D3色散校正和是否加溶剂无关，所以a和c的D3贡献刚好抵消，加了等于白加

谢谢老师，所以a-c的结果与色散校正无关，那么为了让结构较准确，在结构优化时使用D3，则溶剂化下的电子能量(a)为加上D3后的能量，并在在计算真空能量(c)时同样加上D3让他抵销，这样做是否是正确的呢?

作者
Author: wzkchem5 时间: 2021-2-12 10:24

sun0615 发表于 2021-2-11 23:24
谢谢老师，所以a-c的结果与色散校正无关，那么为了让结构较准确，在结构优化时使用D3，则溶剂化下的电子 ...

可以

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