高斯给出的力常数的计算原理求助

幸运兔

请问高斯中对振动计算会出现三部分力常数，想知道他们之间都什么关系
力常数第一次出现的地方（NO1）
，                      1                      2                      3                     A'                     A'                     A'
Frequencies --    410.4464               620.1512              2469.4234
Red. masses --     34.8764                 1.0531                 1.0371
Frc consts  --      3.4617                 0.2386                 3.7262
这个地方给出这三个原子的力（NO2）
Center     Atomic                   Forces (Hartrees/Bohr)
Number     Number              X              Y              Z
-------------------------------------------------------------------
      1       29           0.000015174    0.000000000    0.000005543
      2       16           0.000083506    0.000000000   -0.000007000
      3        1          -0.000098680    0.000000000    0.000001456
-------------------------------------------------------------------

接着又分别给出笛卡尔坐标和内坐标下的力常数，（NO3）

Force constants in Cartesian coordinates:
                1             2             3             4             5
      1  0.754426D-02
      2  0.000000D+00 -0.283154D-03
      3  0.139549D-01  0.000000D+00  0.138514D+00
      4 -0.464779D-02  0.000000D+00 -0.124626D-01  0.214795D+00
      5  0.000000D+00  0.282784D-03  0.000000D+00  0.000000D+00 -0.815863D-03
      6 -0.286608D-02  0.000000D+00 -0.130844D+00  0.518039D-01  0.000000D+00
      7 -0.289647D-02  0.000000D+00 -0.149236D-02 -0.210147D+00  0.000000D+00
      8  0.000000D+00  0.000000D+00  0.000000D+00  0.000000D+00  0.533080D-03
      9 -0.110888D-01  0.000000D+00 -0.766951D-02 -0.393413D-01  0.000000D+00
                6             7             8             9
      6  0.149949D+00
      7 -0.489378D-01  0.213044D+00
      8  0.000000D+00  0.000000D+00 -0.533449D-03
      9 -0.191049D-01  0.504301D-01  0.000000D+00  0.267744D-01
FormGI is forming the generalized inverse of G from B-inverse, IUseBI=4.
Force constants in internal coordinates:
                1             2             3
      1  0.139774D+00
      2  0.345149D-02  0.225477D+00
      3  0.204044D-01  0.190794D-01  0.969444D-01

就想知道，笛卡尔坐标和内坐标下的力常数和给出的力那个表之间有什么样的关系，如何计算的，
非常感谢指导

granvia

力常数是能量关于位置坐标的二阶导数，因此取决于位置坐标的定义。所以，在简正坐标（每个坐标往往对应为一群原子），笛卡尔坐标（每个只对应一个原子），内坐标下的力常数是不一样的。之间的关系通过二阶微分的坐标变换应该能够推出吧，比如雅克比行列式。不过没做过这方面的东西

让你变成回忆

granvia 发表于 2019-11-28 13:31
力常数是能量关于位置坐标的二阶导数，因此取决于位置坐标的定义。所以，在简正坐标（每个坐标往往对应为一 ...

sob老师的这篇博文你仔细看了就知道基本的原理了。
http://sobereva.com/328

幸运兔

granvia 发表于 2019-11-28 13:31
力常数是能量关于位置坐标的二阶导数，因此取决于位置坐标的定义。所以，在简正坐标（每个坐标往往对应为一 ...

非常感谢您的解答。力常数是能量对坐标的二阶导数，那我可以这样理解么：力对坐标一阶导数就可以得到力常数，也就是将NO2中力的数据对坐标进行一阶求导，也就是NO2中的力除以坐标X,Y,Z就应该得到NO3中笛卡尔坐标下的数据对不？但是我找到坐标                        Standard orientation:
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
      1         29           0        0.028667   -0.772854    0.000000
      2         16           0        0.028667    1.296145    0.000000
      3          1           0       -1.290034    1.674452    0.000000
进行相除，得到的数据和NO3中不一致。也就是用 0.000015174 / 0.028667期望得到0.754426D-02，但是不是，特别苦恼，非常感谢解答

幸运兔

让你变成回忆 发表于 2019-11-28 13:39
sob老师的这篇博文你仔细看了就知道基本的原理了。
http://sobereva.com/328

非常感谢您的回答，您给的链接也正是我想去学习的内容，谢谢推荐，但是我理解的是Hess2freq这个程序的源代码解释了最基本的问题，但是我在链接里下载的确实两个计算振动的高斯文件，请问你下载到了这个程序了吗

让你变成回忆

幸运兔 发表于 2019-11-28 14:59
非常感谢您的回答，您给的链接也正是我想去学习的内容，谢谢推荐，但是我理解的是Hess2freq这个程序的源 ...

链接里面已经给了源程序以及intro的PDF文件啊。

幸运兔

让你变成回忆 发表于 2019-11-28 15:41
链接里面已经给了源程序以及intro的PDF文件啊。

原谅我，仔细看了，还是没有找到PDF，我发现是2016写的，现在找不到了呢

sobereva

幸运兔 发表于 2019-11-28 15:46
原谅我，仔细看了，还是没有找到PDF，我发现是2016写的，现在找不到了呢

intro.pdf明明就在博文里的程序包里

sobereva

幸运兔 发表于 2019-11-28 14:54
非常感谢您的解答。力常数是能量对坐标的二阶导数，那我可以这样理解么：力对坐标一阶导数就可以得到力常 ...

有限差分求数值导数和除以坐标完全是两码事
去看有限差分的基础知识

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容（仔细看http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的标题“高斯力常数的计算原理”改了，以后务必注意，下次将扣分处理。

幸运兔

sobereva 发表于 2019-11-29 07:56
有限差分求数值导数和除以坐标完全是两码事
去看有限差分的基础知识

多谢老师

幸运兔

sobereva 发表于 2019-11-29 07:55
intro.pdf明明就在博文里的程序包里

谢谢老师，从新下载，得到了程序和intro.pdf，谢谢您辛苦的付出。

幸运兔

让你变成回忆 发表于 2019-11-28 15:41
链接里面已经给了源程序以及intro的PDF文件啊。

谢谢老师，再次下载，找到相应的文件，感谢告知

小庄师弟

请问一下，我的.out文件中没有NO3力常数矩阵，请问你是在计算的时候加了什么关键字吗？

zjxitcc

小庄师弟 发表于 2021-10-29 16:47
请问一下，我的.out文件中没有NO3力常数矩阵，请问你是在计算的时候加了什么关键字吗？

你不加freq、只算个单点，那当然就没有力常数了。另外，加freq之前要确保你的结构已经优化过了。

小庄师弟

计算了，头文件#p opt freq B3LYP/6-311g(d,p)  scale=0.9679 em=GD3BJ 。从GaussView中画的小分子然后再计算就有力常数矩阵，从CCDC下载到的cif晶体文件，经MS-Build-Symmetry-Nonperiodic Superstructure转化，再用GaussView保存成gjf文件进行计算，然后输出文件里就没有力常数矩阵。