分子团簇计算偶极矩的泛函基组

大源

老师，我想问一下，计算大约20多个原子的团簇偶极矩，B3lyp aug-cc-pvdz 进行优化得到的偶极矩结果和 ccsd(t) aug-cc-pvdz 计算单点得到的偶极矩的结果有点差距，是否应该采用计算单点的结果作为分子偶极矩结果

zjxitcc

（1）CCSD(T)计算偶极矩需要计算CCSD(T)弛豫密度，只有极少的程序支持（Molpro, CFOUR, PySCF）计算这个，并且三个程序的结果尚未实现严格一致。因此我怀疑你搞错了，用的是Hartree-Fock (HF)偶极矩在做对比。例如Gaussian输入文件写CCSD(T)，做单点计算，这样获得的就是HF偶极矩。

（2）如果你没获得CCSD(T)偶极矩，请不用看这一点；如果你确认拿到了正确的CCSD(T)偶极矩，那你需要讲明“有点差距”是什么意思，两个结果分别是多少，相差多少，这样大家可以帮你判断。

（3）“20多个原子的团簇偶极矩”信息过少，例如二肽的计算难度与Fe_20金属团簇的计算难度是天壤之别，CCSD(T)完全无法描述Fe_20金属团簇，因此讨论偶极矩也就没什么意义。

大源

zjxitcc 发表于 2026-8-24 09:46
（1）CCSD(T)计算偶极矩需要计算CCSD(T)弛豫密度，只有极少的程序支持（Molpro, CFOUR, PySCF）计算这个， ...

谢谢老师指出这个问题。我之前确实是直接在 Gaussian 中使用 ccsd(t) aug-cc-pvtz 做单点，然后读取输出文件中的 dipole moment，并把这个值理解成了 CCSD(T) 水平的偶极矩。看了您的解释之后，我意识到这里可能存在误解：Gaussian 中虽然总能量已经计算到 CCSD(T) 水平，但默认输出的偶极矩不一定包含 CCSD(T) 相关密度 。我计算的分子团簇都是前三周期的金属团簇，B3LYP/aug-cc-pVDZ 优化计算中 该结构的偶极矩大约是0.17debye, 在ccsd(t) aug-cc-pvtz 计算单点得到的dipole moment 是0.6Debye 左右, 计算结果有点差距。

zjxitcc

大源 发表于 2026-8-24 10:00
谢谢老师指出这个问题。我之前确实是直接在 Gaussian 中使用 ccsd(t) aug-cc-pvtz 做单点，然后读取输出 ...

（1）你用的是HF偶极矩，不宜再说“在ccsd(t) aug-cc-pvtz 计算单点得到的dipole moment”。既然你没有CCSD(T)偶极矩，也就无从比较，原问题不成立了。
（2）对于金属团簇，你用B3LYP做opt freq，一定要读取优化好的几何结构和波函数，检验波函数稳定性，以免误用不稳定波函数 计算频率、计算偶极矩。

大源

zjxitcc 发表于 2026-8-24 10:07
（1）你用的是HF偶极矩，不宜再说“在ccsd(t) aug-cc-pvtz 计算单点得到的dipole moment”。既然你没有CC ...

好的，谢谢老师指教。