计算化学公社

标题: 算弱相互作用和算单点能时所用泛函基组要保持一致吗 [打印本页]

作者
Author: yuna 时间: 2023-3-30 20:41
标题: 算弱相互作用和算单点能时所用泛函基组要保持一致吗
1.各位老师好，我算的是硅烷还原二氧化碳的过程，体系中一个中间体是通过弱相互作用结合的，所有的结构优化频率分析是在B3LYP-D3(BJ)/6-311G(d,p)下算的，A和B的弱相互作用能是在M06-2X-D3/6-311+G(d,p)下计算的。
原本单点能计算想用的泛函基组是M062x/def2TZVP，但是一开始就涉及到弱相互作用的问题了，我要不要用M06-2X-D3/6-311+G(d,p)算体系中其他结构的单点能呢？

2.弱相互作用能的计算是只用电子能量求差值对吗？
谢谢

作者
Author: wzkchem5 时间: 2023-3-30 21:43
要保持一致，因为形成弱相互作用复合物这一步要和后续反应画在同一张图里，至少也要在同一篇文章里，而且可能涉及彼此之间的加减，如果理论级别不统一不能互相加减。
用电子能求差值得到的结果适合用来做能量分解分析等，如果要和实验的能量对比，需要加零点能校正以及对内能的热校正；如果要和实验的平衡常数等对比，需要加自由能校正

作者
Author: yuna 时间: 2023-3-30 22:37

wzkchem5 发表于 2023-3-30 21:43
要保持一致，因为形成弱相互作用复合物这一步要和后续反应画在同一张图里，至少也要在同一篇文章里，而且可 ...

明白了，谢谢老师的解答

作者
Author: yuna 时间: 2023-3-31 11:47

wzkchem5 发表于 2023-3-30 21:43
要保持一致，因为形成弱相互作用复合物这一步要和后续反应画在同一张图里，至少也要在同一篇文章里，而且可 ...

再请问老师下，我体系中会有一些阳离子或者碱金属离子，单点能的计算也用M062X-D3/6-311+G(d,p)这样合适吗？我看sob老师说阳离子不能加弥散函数。

作者
Author: ChemG 时间: 2023-3-31 16:10

yuna 发表于 2023-3-31 11:47
再请问老师下，我体系中会有一些阳离子或者碱金属离子，单点能的计算也用M062X-D3/6-311+G(d,p)这样合适 ...

阴离子因为电子在空间分布更广所以要加弥散。阳离子不是不能加，应该说相比阴离子，作用不是那么大，并且对计算量、SCF收敛等都有影响

作者
Author: yuna 时间: 2023-3-31 16:59

ChemG 发表于 2023-3-31 16:10
阴离子因为电子在空间分布更广所以要加弥散。阳离子不是不能加，应该说相比阴离子，作用不是那么大，并且 ...

好的，谢谢

作者
Author: wzkchem5 时间: 2023-3-31 19:22

yuna 发表于 2023-3-31 04:47
再请问老师下，我体系中会有一些阳离子或者碱金属离子，单点能的计算也用M062X-D3/6-311+G(d,p)这样合适 ...

阳离子是通常不需加弥散函数，而不是不能加。“阳离子不加弥散函数”的优先级低于“阴离子加弥散函数”的优先级，而两者的优先级都低于“加不加弥散函数必须统一”的优先级

作者
Author: sobereva 时间: 2023-4-1 09:45

yuna 发表于 2023-3-31 11:47
再请问老师下，我体系中会有一些阳离子或者碱金属离子，单点能的计算也用M062X-D3/6-311+G(d,p)这样合适 ...

我从来没说过什么“阳离子不能加弥散函数”
不能加和不需要加完全是两码事

作者
Author: yuna 时间: 2023-4-2 13:21
本帖最后由 yuna 于 2024-3-15 17:13 编辑

sobereva 发表于 2023-4-1 09:45
我从来没说过什么“阳离子不能加弥散函数”
不能加和不需要加完全是两码事

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