计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › CCSD(T)方法计算单点读出的偶极矩与用5点数值差分公式计 ...

CCSD(T)方法计算单点读出的偶极矩与用5点数值差分公式计算得出的偶极矩区别在哪里？

查看数: 1295 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

liurq

发布时间: 2024-7-8 19:10

正文摘要:

请问各位老师，CCSD(T)方法计算单点读出的偶极矩与用5点数值差分公式计算得出的偶极矩区别在哪里？为什么数值上会有较大的差异？我用Gaussian先进行了单点计算，计算级别为CCSD(T)/def2TZVPP，在log文件里读出偶极 ...

回复 Reply

sobereva 发表于 Post on 2024-7-8 23:07:30

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “偶极矩问题” 改了，以后务必注意！

liurq 发表于 Post on 2024-7-8 20:08:29

本帖最后由 liurq 于 2024-7-8 20:11 编辑

zjxitcc 发表于 2024-7-8 20:05
建议：使用外加电场时最好加上nosymm关键词做计算，不然分子朝向有点迷惑。
回答：算个CCSD(T)单点能然 ...

好的，谢谢老师，很感谢您的解惑，我会加以修改

zjxitcc 发表于 Post on 2024-7-8 20:05:32

本帖最后由 zjxitcc 于 2024-7-8 20:07 编辑

liurq 发表于 2024-7-8 19:40
老师，请问我在log文件的这里读出的CCSD(T)的能量，再基于这个公式dE/dFx = [E(-2Dx) -8E(-Dx)+8E(+Dx)-E ...

建议：使用外加电场时最好加上nosymm关键词做计算，不然分子朝向有点迷惑。
回答：算个CCSD(T)单点能然后输出文件读偶极矩，这肯定是错的。你的另一种做法（能量对外电场的导数）思路是对的，但我不知道有没有意义，因为有两点我不确定：（1）没加nosymm，不知道是什么朝向，有没被旋转；（2）Gaussian的CCSD(T)是否严格考虑了外加电场的响应。

liurq 发表于 Post on 2024-7-8 19:35:09

本帖最后由 liurq 于 2024-7-8 19:43 编辑

zjxitcc 发表于 2024-7-8 19:24
高斯CCSD(T)没有非弛豫密度，无法解析计算偶极矩，你在log文件中看到的是HF方法的低精度偶极矩。

谢谢您答疑，请问这两种方法做出来的数值都没有意义吗？可是它们又为什么会差别这么明显呢？

wjc404 发表于 Post on 2024-7-8 19:24:39

你从高斯的log文件里直接读出的偶极矩，是在"population analysis using the SCF density"下面吧，这个是基于HF的密度做的。高斯做不了CCSD(T)的密度，CFOUR可以。

zjxitcc 发表于 Post on 2024-7-8 19:24:28

高斯CCSD(T)没有非弛豫密度，无法解析计算偶极矩，你在log文件中看到的是HF方法的低精度偶极矩。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register