针对甲醛使用不用方法计算偶极矩得到的结果一样？

mhliu · 发表于 Post on 2023-11-20 21:38:57

老师您好，请问我针对同一个甲醛构型，基于同一个基组分别换用不同的方法算单点，得到的偶极矩为什么是一样的？
我理解的是不同方法下波函数改变了，虽然构型不变但是偶极矩应该是变化的？
然后我还尝试用一样的方法，变换基组，偶极矩就发生变化了，具体过如下图

啊不错的飞过海 · 发表于 Post on 2023-11-20 22:25:26

本帖最后由啊不错的飞过海于 2023-11-21 21:27 编辑

你用来计算偶极矩的密度矩阵驰豫过吗？或者说，你是否加了density的关键词？
一般来讲偶极矩这样相应电场的性质使用relaxed的密度算出来物理意义更明确，更多信息可以参考http://bbs.keinsci.com/thread-5738-1-1.html。

ionexchangeC · 发表于 Post on 2023-11-20 22:44:26

需要加density关键词，否则还是HF级别下的偶极矩
另外顺便说一句，高斯其实也能计算CCSD(T)级别的偶极矩，用数值方法就行了（在X,Y,Z方向各加一个很小的电场）

sobereva · 发表于 Post on 2023-11-21 02:14:07

你若想把偶极矩算得很准，用这样的基组根本不合适。仔细看此文
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

mhliu · 发表于 Post on 2023-11-21 10:05:08

ionexchangeC 发表于 2023-11-20 22:44
需要加density关键词，否则还是HF级别下的偶极矩
另外顺便说一句，高斯其实也能计算CCSD(T)级别的偶极矩， ...

好的谢谢~加了density变化了

hebrewsnabla · 发表于 Post on 2023-11-21 20:50:14

啊不错的飞过海发表于 2023-11-20 22:25
你用来计算偶极矩的密度矩阵驰豫过吗？或者说，你是否加了density=relaxed/unrelaxed的关键词？
一般来讲 ...

高斯应该没有 density=relaxed/unrelaxed 关键词

啊不错的飞过海 · 发表于 Post on 2023-11-21 21:28:10

hebrewsnabla 发表于 2023-11-21 20:50
高斯应该没有 density=relaxed/unrelaxed 关键词

感谢提醒，修改了2楼

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

[Gaussian/gview] 针对甲醛使用不用方法计算偶极矩得到的结果一样？