计算化学公社

标题: 关于范德华系数C6计算问题 [打印本页]

作者
Author: hiramqdl 时间: 2023-5-29 14:38
标题: 关于范德华系数C6计算问题

想请教一下，论坛里面有人知道范德华常数C6是如何计算的吗？查阅了很久依然没有收获。唯一相关的是下图中所示的算法，但受限于已知辐射体的电离能。想请教一下，是否有别的计算公式？
另外，计算过程中涉及到光谱相关知识，为什么这里的电离能的单位是cm-1，电离能不应该是J/mol吗？

作者
Author: wzkchem5 时间: 2023-5-29 16:08
这个估算方法极其粗糙，不建议使用。一个比这个算得还快、但是更准的方法是做一个DFT-D3（或DFT-D4）计算，并把print level调得很大，这样程序会在DFT-D级别下计算这个分子的C6。另一个更准但是耗时较多的方法是用TDDFT算，参见https://doi.org/10.1063/1.3382344里TDDFT级别下的C6的计算方法。
能量如果表示成cm-1的单位，指的是具有这么多能量的光子的波数。波数固然严格讲不是能量，但是因为两者之间有简单的正比关系，因此习惯上常常把能量表成以cm-1为单位，却不说它是“具有这么多能量的光子的波数”，而直接说这个波数就是能量。

作者
Author: hiramqdl 时间: 2023-5-29 16:28
感谢老师的回答，也就是说C6是可以通过Gaussian计算求解的，明白了。
另外还有一个问题，计算共振态展宽的时候，需要共振态振子强度，在本站看到了sob老师提供的相关公式，但是依然不清楚应该如何计算，尤其是对于跃迁电偶极矩的计算，用Gaussian可以实现吗？

作者
Author: wzkchem5 时间: 2023-5-29 16:46

hiramqdl 发表于 2023-5-29 09:28
感谢老师的回答，也就是说C6是可以通过Gaussian计算求解的，明白了。
另外还有一个问题，计算共振态展宽的 ...

跃迁电偶极矩当然可以计算，高斯（以及所有其他主流量化软件）TDDFT计算默认就会计算出来。
但是高斯可能不是特别适合算共振态，如果你指的是continuous ionization background里的共振峰的话，可能需要用支持complex absorbing potential的程序计算

作者
Author: hiramqdl 时间: 2023-5-29 21:45
振子强度和共振态振子强度，有什么区别呢，或者这属于哪方面的知识？打算计算NEC辐射模型，在研究的过程中发现很多难点。

作者
Author: hiramqdl 时间: 2023-5-29 21:48
我只用高斯算电偶极矩，然后计算振子强度，可以理解此时的振子强度=共振态振子强度吗？

作者
Author: sobereva 时间: 2023-5-29 21:59

hiramqdl 发表于 2023-5-29 16:28
感谢老师的回答，也就是说C6是可以通过Gaussian计算求解的，明白了。
另外还有一个问题，计算共振态展宽的 ...

引用我的培训班ppt的时候不得把顶端的版权信息截掉，这里也明确说了：
北京科音办的培训班FAQ
http://bbs.keinsci.com/thread-5098-1-1.html

已经把你的图删了，下次注意！

作者
Author: sobereva 时间: 2023-5-29 22:10
这是严格的计算C6的方法，利用Dalton可以直接算

(, 下载次数 Times of downloads: 18)

DFT-D3、D4只不过是经验近似得到C6。具体来说，DFT-D4是通过配位数、原子电荷对事先算好的参考体系的原子动态极化率进行scaling得到的实际体系中原子动态极化率，然后数值积分这个式子计算C6,8系数

“C6是可以通过Gaussian计算求解的”说法不确切。Grimme的独立的DFT-D程序能给你输出某些中间信息不代表Gaussian也能给你输出

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