请问高斯里面的Forces怎么算的？怎么正负号对不上呢？

413

加了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
-------------------------------------------------------------------

sobereva

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

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

413

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

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