24小时热门版块排行榜    

CyRhmU.jpeg
查看: 348  |  回复: 1
当前主题已经存档。

九层楼

木虫 (小有名气)

[交流] 【求助】Gaussian输出文件中,Cartesian Force和Redundant Variable的Force的关系

Gaussian输出文件中、每个SCF中、
(1)先输出坐标,例如:
CODE:
   
     56         8           0        1.029963   -0.957091    0.017046
    129          8           0        0.504390   -1.405935    1.710342

(2)SCF Done后、会先输出Cartesian Forces,例如:
CODE:
   
56        8          -0.001336908   -0.001073193    0.002309701
129        8           0.000304343    0.000710809   -0.002832933

(3)后面又按照Redundant Variable输出包含Force等的信息
CODE:
   
Variable       Old X    -DE/DX   Delta X   Delta X   Delta X     New X
                                 (Linear)    (Quad)   (Total)
R72          3.45615  -0.00226   0.00000  -0.00081  -0.00081   3.45534

这儿R72是 R(56,129) 、由上面那2个原子组成

请问从Cartesian Force如何得到这个R72的值啊?
之所以想知道这个、是因为、虽然Gaussian同时输出Cartesian Force和Redundant Variable的Force、但见过的其他软件还是输出Cartesian Force的多,如果能了解了如何转化、就可以任意转化为原子间的力、感觉可能对理解也许会有些帮助?

我尝试过
把O56的力(作为矢量)、和O129的力、分别投影在、按O56->O129的坐标得到的矢量上,但是、对得到的2个力、加了求Length是9.9E-5、或减了求是0.00567、结果都不是0.00226(带符号-0.00226),很困惑

谢谢大家的帮助!

[ Last edited by 九层楼 on 2010-1-20 at 19:47 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

九层楼

木虫 (小有名气)


顶一下
2楼2010-01-22 17:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 九层楼 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见