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

幸运兔

请问高斯中对振动计算会出现三部分力常数，想知道他们之间都什么关系
力常数第一次出现的地方（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

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

小庄师弟

zjxitcc 发表于 2021-11-2 18:36
高斯输出文件内容多少，是会随原子数而改变的。你可以用fch文件里的力常数，不管原子数多少 一定是会输出 ...

谢谢老师，我去试试

zjxitcc

小庄师弟 发表于 2021-10-30 20:00
有这个可能，我测试的几个小分子（6、7个原子）是显示这部分内容的，但是算的其他几个比较大的分子（有40 ...

高斯输出文件内容多少，是会随原子数而改变的。你可以用fch文件里的力常数，不管原子数多少 一定是会输出的，注意方位朝向和单位。

小庄师弟

paramecium86 发表于 2021-10-30 16:59
我记得以前算的时候感觉log文件里的force constants in cartesian coordinates这部分 分子超过一定大小好 ...

有这个可能，我测试的几个小分子（6、7个原子）是显示这部分内容的，但是算的其他几个比较大的分子（有40个原子的，有100左右原子的）不显示这部分内容。

paramecium86

小庄师弟 发表于 2021-10-29 20:23
计算时的gjf文件

我记得以前算的时候感觉log文件里的force constants in cartesian coordinates这部分 分子超过一定大小好像就不输出了。
比如记得以前用同样的关键词 算苯环就输出这段儿，但是 算蒽和芘这段儿就不见了。

小庄师弟

wzkchem5 发表于 2021-10-29 20:37
可能是报错了。
上传完整输出文件

我算了好几个分子都是这种情况，也用GaussView画了几个分子测试了一下，是有“ Force constants in Cartesian coordinates: ...”这部分内容的。

wzkchem5

小庄师弟 发表于 2021-10-29 13:23
计算时的gjf文件

可能是报错了。
上传完整输出文件

小庄师弟

计算时的gjf文件

小庄师弟

wzkchem5 发表于 2021-10-29 18:59
因为由其他文件格式转换成的gjf文件不会自动写上opt freq，需要你手动添加。

老师，您说加上opt freq是指哪里的？在计算的时候gjf文件里写了opt freq

小庄师弟

zjxitcc 发表于 2021-10-29 19:07
保存成gjf文件，然后你不打开检查关键词的么？

关键词已经贴了#p opt freq B3LYP/6-311g(d,p)  scale=0.9679 em=GD3BJ，老师您好，我的表达有一点错误，用MS导出的时mol格式，然后用GaussView打开保存为gjf的

zjxitcc

小庄师弟 发表于 2021-10-29 17:36
计算了，头文件#p opt freq B3LYP/6-311g(d,p)  scale=0.9679 em=GD3BJ 。从GaussView中画的小分子然后再计 ...

保存成gjf文件，然后你不打开检查关键词的么？

wzkchem5

小庄师弟 发表于 2021-10-29 10:36
计算了，头文件#p opt freq B3LYP/6-311g(d,p)  scale=0.9679 em=GD3BJ 。从GaussView中画的小分子然后再计 ...

因为由其他文件格式转换成的gjf文件不会自动写上opt freq，需要你手动添加。

小庄师弟

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

zjxitcc

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

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

小庄师弟

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