计算化学公社

标题: 高斯尾部输出PolarDeriv含义 [打印本页]

作者
Author: nuke 时间: 2017-6-1 08:10
标题: 高斯尾部输出PolarDeriv含义
我想问一下 gaussian 尾部输出的问题
hf1.gjf

%mem=60MB
%nproc=1
# opt=(calcall,tight) freq=raman rb3lyp/6-311++g(d,p)  nosymm

Title Card Required

0 1
F  0  0. 0. 0.
H  0  0. 0. 0.92223200

hf1.out 结尾有个
PolarDeriv=0.,0.,0.,-1.3075358,0.,0.,0.,0.,0.,0.,-1.3075357,0.,-0.3696203,0.,
-0.3696203,0.,0.,-5.7465674,0.,0.,0.,1.3075358,0.,0.,0.,0.,0.,0.,1.307
5357,0.,0.3696203,0.,0.3696203,0.,0.,5.7465674

如果我要看氢原子的PolarDeriv应该是从第19个到第36个那些数表示。
他们的顺序应该是 (\partial \alpha_{xx})/(\partial x), (\partial \alpha_{xy})/(\partial x), (\partial \alpha_{yy})/{\partial x), (\partial \alpha_{xz})/(\partial x)
。。。。类似的这种顺序吧？这里的\alpha指polarizability。  谢谢

作者
Author: sobereva 时间: 2017-6-1 08:59
做有限差分验证一下就知道了

作者
Author: nuke 时间: 2017-6-1 09:42

sobereva 发表于 2017-6-1 08:59
做有限差分验证一下就知道了

你的意思是改变一下氢原子的z轴坐标，再看一下结果吧

作者
Author: sobereva 时间: 2017-6-1 17:53

nuke 发表于 2017-6-1 09:42
你的意思是改变一下氢原子的z轴坐标，再看一下结果吧

对

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