24小时热门版块排行榜

返回列表

当前主题已经存档。

九层楼

木虫 (小有名气)

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 ]

回复此楼

» 猜你喜欢

超导转变温度NaN 已经有4人回复
华南师范大学先进光电子研究院电子纸团队诚聘2026级博士研究生已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有200人回复
QE利用声子计算的一系列lambda和Tc如何选已经有8人回复
理论计算合作已经有2人回复
小木虫的论文辅导靠谱吗？有没有用过的同学？已经有1人回复
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有已经有0人回复
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后已经有81人回复
2026博士申请—化学能源方向，本人熟练分子动力学模拟的使用已经有4人回复
书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复

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

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

九层楼

木虫 (小有名气)

顶一下

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

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

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

返回列表