24小时热门版块排行榜

>论坛更新日志 (3977)
>基金申请 (125)
>导师招生 (111)
>文献求助 (110)
>虫友互识 (80)
>硕博家园 (78)
>博后之家 (54)
>考博 (41)
>找工作 (35)
>论文投稿 (32)
>休闲灌水 (30)
>论文道贺祈福 (27)
>教师之家 (25)
>公派出国 (23)
>考研 (21)
>有机交流 (15)

返回列表

当前主题已经存档。

九层楼

木虫 (小有名气)

QC强帖: 2
应助: 0 (幼儿园)
金币: 3273.2
帖子: 260
在线: 204.5小时
虫号: 86618
注册: 2005-08-15
性别: GG
专业: 理论和计算化学

[交流] 【求助】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 ]

回复此楼

» 猜你喜欢

1楼 2010-01-20 19:37:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

九层楼

木虫 (小有名气)

顶一下

2楼2010-01-22 17:13:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主九层楼的主题更新

返回列表