计算化学公社

标题: 请问高斯里面的Forces怎么算的？怎么正负号对不上呢？ [打印本页]

作者
Author: 413 时间: 2024-3-19 22:28
标题: 请问高斯里面的Forces怎么算的？怎么正负号对不上呢？
本帖最后由 413 于 2024-3-19 22:30 编辑

加了nosymm关键词
1. 力是能量的梯度，以X方向为例，所有的梯度都应该是负数的呀，因为能量是负的，所有的X坐标都是正的。
为什么第一个X方向的力是负的，其余是正的呢？

2. 如果不加nosymm，分子会自动转一下。这个时候会出现Standard orientation，那么受力是按照Standard orientation算的对吗？

                        Input orientation:
---------------------------------------------------------------------
Center    Atomic    Atomic          Coordinates (Angstroms)
Number    Number    Type          X          Y          Z
---------------------------------------------------------------------
   1       8          0       0.846875 0.237500 0.000000
   2       1          0       1.806875 0.237500 0.000000
   3       1          0       0.526420 1.142436 0.000000
---------------------------------------------------------------------
。。。。。。
SCF循环
SCF Done:  E(RPBE-PBE) =  -76.3190708107
。。。。。。
-------------------------------------------------------------------
Center    Atomic                Forces (Hartrees/Bohr)
Number    Number             X             Y             Z
-------------------------------------------------------------------
   1       8       -0.017025505 -0.024104840 -0.000000000
   2       1          0.012450069 0.009274071 0.000000000
   3       1          0.004575435 0.014830769 -0.000000000
-------------------------------------------------------------------

作者
Author: sobereva 时间: 2024-3-20 09:21
先搞清楚梯度是怎么定义的、什么叫导数，试图自己通过能量对核坐标做有限差分得到受力，自然就能理解怎么回事。显然跟坐标的正负毫无关系

不用nosymm时受力给的是标准朝向的

作者
Author: 413 时间: 2024-3-20 10:39

sobereva 发表于 2024-3-20 09:21
先搞清楚梯度是怎么定义的、什么叫导数，试图自己通过能量对核坐标做有限差分得到受力，自然就能理解怎么回 ...

我搞糊涂了，是每个原子的受力。

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