计算化学公社

标题: 计算二钠离子盐体系偶极矩，输出结果没有dipole moment [打印本页]

作者
Author: 吕安琪 时间: 2019-7-7 15:18
标题: 计算二钠离子盐体系偶极矩，输出结果没有dipole moment
本帖最后由吕安琪于 2019-7-13 22:18 编辑

各位老师好，我在计算二钠离子盐体系偶极矩时，输入文件如下，但是输出结果没有dipole moment，请问各位老师应该如何正确计算该体系的偶极矩？

%chk=dandian.chk
%mem=30GB
%nproc=12
# m062x/cc-pvdz charge

test

-2 1
笛卡尔坐标

钠离子坐标 1.0
钠离子坐标 1.0

作者
Author: sobereva 时间: 2019-7-7 19:22
不可能没有输出偶极矩。说不清楚就把输出文件压缩后上传

当前这种基组根本不可能给出定性合理的偶极矩，更何况当前还带负电荷，不加弥散函数算的结果毫无意义。
算偶极矩一定要用带弥散函数的基组，至少def2-SVPD，有条件的用def2-TZVPD或aug-cc-pVTZ

作者
Author: 吕安琪 时间: 2019-7-13 22:15

sobereva 发表于 2019-7-7 19:22
不可能没有输出偶极矩。说不清楚就把输出文件压缩后上传

当前这种基组根本不可能给出定性合理的偶极矩， ...

谢谢老师，按照您的建议已经更改了基组，所得.log为上传文件，但是输出的偶极矩信息在gaussview中观察发现过短，请问老师这种情况属于正常吗

作者
Author: sobereva 时间: 2019-7-14 01:30

吕安琪发表于 2019-7-13 22:15
谢谢老师，按照您的建议已经更改了基组，所得.log为上传文件，但是输出的偶极矩信息在gaussview中观察发 ...

你的体系本来就是中心对称的，当然偶极矩理应为0（不精确为0是因为你没对体系做对称化）

作者
Author: 吕安琪 时间: 2019-7-15 10:18

sobereva 发表于 2019-7-14 01:30
你的体系本来就是中心对称的，当然偶极矩理应为0（不精确为0是因为你没对体系做对称化）

知道了谢谢老师

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