计算化学公社

标题: molpro如何计算uccsd(t)-f12方法的偶极矩 [打印本页]

作者
Author: dlldwtz 时间: 2023-5-24 10:10
标题: molpro如何计算uccsd(t)-f12方法的偶极矩
molpro在使用uccsd(t)-f12方法计算OH自由基时，无法使用expec,dm这样的关键词计算偶极矩，有什么方法可以通过输出文件推导得到偶极矩吗？

作者
Author: dlldwtz 时间: 2023-5-24 10:50
已经解决，利用有限场关键词dip可以计算

作者
Author: zyy1997 时间: 2023-5-24 15:31
你好，可以上传一下输入文件吗？谢谢！

作者
Author: dlldwtz 时间: 2023-5-24 15:39

zyy1997 发表于 2023-5-24 15:31
你好，可以上传一下输入文件吗？谢谢！

***,oh
memory,100,m

nosym
geometry={
H
O H R
}
R = 0.7525 Ang

basis=avtz

dip,,,0.00001
rhf
uccsd(t)-f12a
ep=energy

dip,,,-0.00001
rhf
uccsd(t)-f12a
en=energy

de=ep-en
dm=de/2/0.00001

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