计算化学公社

标题: 针对甲醛使用不用方法计算偶极矩得到的结果一样？ [打印本页]

作者
Author: mhliu 时间: 2023-11-20 21:38
标题: 针对甲醛使用不用方法计算偶极矩得到的结果一样？
老师您好，请问我针对同一个甲醛构型，基于同一个基组分别换用不同的方法算单点，得到的偶极矩为什么是一样的？
我理解的是不同方法下波函数改变了，虽然构型不变但是偶极矩应该是变化的？
然后我还尝试用一样的方法，变换基组，偶极矩就发生变化了，具体过如下图

作者
Author: 啊不错的飞过海 时间: 2023-11-20 22:25
本帖最后由啊不错的飞过海于 2023-11-21 21:27 编辑

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

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

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

作者
Author: mhliu 时间: 2023-11-21 10:05

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

好的谢谢~加了density变化了

作者
Author: hebrewsnabla 时间: 2023-11-21 20:50

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

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

作者
Author: 啊不错的飞过海 时间: 2023-11-21 21:28

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

感谢提醒，修改了2楼

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