通过gaussian进行限制性优化后没有虚频，但是拿输出结构进行优化，结构改变

wuxzd

限制性优化的输入文件在下面


#p opt=modredundant freq 6-311g(d,p) pop=nbo empiricaldispersion=gd3 pbe1pbe


123


0 1
C                  0.10145500   -0.73671300   -0.20876200
C                 -0.72966700   -1.61989200    0.49329100
C                 -0.52069500   -2.98769400    0.43617700
C                  0.52069500   -3.51718800   -0.32699000
C                  1.34552300   -2.64934500   -1.04091600
C                  1.13432700   -1.28290000   -0.96825100
C                 -0.52069500    3.51718800   -0.32699000
C                 -1.34552300    2.64934500   -1.04091600
C                 -1.13432700    1.28290000   -0.96825100
C                  0.72966700    1.61989200    0.49329100
C                  0.52069500    2.98769400    0.43617700
C                 -0.10145500    0.73671300   -0.20876200
N                 -0.67483200    4.90146300   -0.33974100
N                  0.67483200   -4.90146300   -0.33974100
C                  1.72955500   -5.59322700   -0.76690300
C                  1.79443700   -6.97611000   -0.77239500
C                  0.65980200   -7.77930800   -0.29637400
O                 -0.39203900   -7.26409200    0.13195600
C                 -1.72955500    5.59322700   -0.76690300
C                 -1.79443700    6.97611000   -0.77239500
C                 -0.65980200    7.77930800   -0.29637400
O                  0.39203900    7.26409200    0.13195600
C                  2.98373800   -7.61339900   -1.24112100
C                  3.09112400   -8.96735100   -1.26590200
C                  1.99598400   -9.76311200   -0.81365200
C                  0.84176000   -9.21015000   -0.35362500
C                 -2.98373800    7.61339900   -1.24112100
C                 -3.09112400    8.96735100   -1.26590200
C                 -1.99598400    9.76311200   -0.81365200
C                 -0.84176000    9.21015000   -0.35362500
C                 -1.86590000   -1.10299500    1.33594000
F                 -2.40416400   -2.06574800    2.09520500
F                 -2.86009800   -0.61247300    0.57123700
F                 -1.47942400   -0.12041200    2.15611200
C                  1.86590000    1.10299500    1.33594000
F                  2.40416400    2.06574800    2.09520500
F                  2.86009800    0.61247300    0.57123700
F                  1.47942400    0.12041200    2.15611200
H                 -1.16030700   -3.65391600    1.00312300
H                  2.14186400   -3.02729900   -1.67095300
H                  1.78323400   -0.61467200   -1.52189800
H                 -2.14186400    3.02729900   -1.67095300
H                 -1.78323400    0.61467200   -1.52189800
H                  1.16030700    3.65391600    1.00312300
H                  0.08249100    5.50502700    0.01376500
H                 -0.08249100   -5.50502700    0.01376500
H                  2.58666000   -5.02223100   -1.11469400
H                 -2.58666000    5.02223100   -1.11469400
H                  3.80626600   -6.98882400   -1.58137800
H                  3.99465400   -9.44667700   -1.62369500
H                  2.09349700  -10.84488000   -0.83914100
H                  0.01178800   -9.81650900   -0.00868000
H                 -3.80626600    6.98882400   -1.58137800
H                 -3.99465400    9.44667700   -1.62369500
H                 -2.09349700   10.84488000   -0.83914100
H                 -0.01178800    9.81650900   -0.00868000


B 18 46 F
B 22 45 F




得到的结果没有虚频，然后在输出文件的基础上进行了先限制性优化，得到的结果却出现了虚频，这该怎么办呢，我想要得到的是限制性优化的结构，因为和势能面扫描的结果吻合



这是第二次优化的输入文件信息：
#p opt freq 6-311g(d,p) empiricaldispersion=gd3 pbe1pbe

456

0 1
C                  0.10145500   -0.73671300   -0.20876200
C                 -0.72966700   -1.61989200    0.49329100
C                 -0.52069500   -2.98769400    0.43617700
C                  0.52069500   -3.51718800   -0.32699000
C                  1.34552300   -2.64934500   -1.04091600
C                  1.13432700   -1.28290000   -0.96825100
C                 -0.52069500    3.51718800   -0.32699000
C                 -1.34552300    2.64934500   -1.04091600
C                 -1.13432700    1.28290000   -0.96825100
C                  0.72966700    1.61989200    0.49329100
C                  0.52069500    2.98769400    0.43617700
C                 -0.10145500    0.73671300   -0.20876200
N                 -0.67483200    4.90146300   -0.33974100
N                  0.67483200   -4.90146300   -0.33974100
C                  1.72955500   -5.59322700   -0.76690300
C                  1.79443700   -6.97611000   -0.77239500
C                  0.65980200   -7.77930800   -0.29637400
O                 -0.39203900   -7.26409200    0.13195600
C                 -1.72955500    5.59322700   -0.76690300
C                 -1.79443700    6.97611000   -0.77239500
C                 -0.65980200    7.77930800   -0.29637400
O                  0.39203900    7.26409200    0.13195600
C                  2.98373800   -7.61339900   -1.24112100
C                  3.09112400   -8.96735100   -1.26590200
C                  1.99598400   -9.76311200   -0.81365200
C                  0.84176000   -9.21015000   -0.35362500
C                 -2.98373800    7.61339900   -1.24112100
C                 -3.09112400    8.96735100   -1.26590200
C                 -1.99598400    9.76311200   -0.81365200
C                 -0.84176000    9.21015000   -0.35362500
C                 -1.86590000   -1.10299500    1.33594000
F                 -2.40416400   -2.06574800    2.09520500
F                 -2.86009800   -0.61247300    0.57123700
F                 -1.47942400   -0.12041200    2.15611200
C                  1.86590000    1.10299500    1.33594000
F                  2.40416400    2.06574800    2.09520500
F                  2.86009800    0.61247300    0.57123700
F                  1.47942400    0.12041200    2.15611200
H                 -1.16030700   -3.65391600    1.00312300
H                  2.14186400   -3.02729900   -1.67095300
H                  1.78323400   -0.61467200   -1.52189800
H                 -2.14186400    3.02729900   -1.67095300
H                 -1.78323400    0.61467200   -1.52189800
H                  1.16030700    3.65391600    1.00312300
H                  0.08249100    5.50502700    0.01376500
H                 -0.08249100   -5.50502700    0.01376500
H                  2.58666000   -5.02223100   -1.11469400
H                 -2.58666000    5.02223100   -1.11469400
H                  3.80626600   -6.98882400   -1.58137800
H                  3.99465400   -9.44667700   -1.62369500
H                  2.09349700  -10.84488000   -0.83914100
H                  0.01178800   -9.81650900   -0.00868000
H                 -3.80626600    6.98882400   -1.58137800
H                 -3.99465400    9.44667700   -1.62369500
H                 -2.09349700   10.84488000   -0.83914100
H                 -0.01178800    9.81650900   -0.00868000

zjxitcc

你第一个任务得到的不是真正的极小点，是坐标限制下的极小点，“没有虚频”也是假象，是程序把坐标限制考虑进去了。你放开坐标限制之后，首先结构的力不为零，你直接算freq看有无虚频意义不大，此时热力学量也没有意义。放开限制后，它结构就是应该会变化，不会变化才是奇怪的（除非运气特别好，做了一个没有用的坐标限制，可以落到真正极小点）。

“和势能面扫描的结果吻合”这话比较模糊，不知道想讲什么。

wuxzd

zjxitcc 发表于 2023-11-23 15:03
你第一个任务得到的不是真正的极小点，是坐标限制下的极小点，“没有虚频”也是假象，是程序把坐标限制考虑 ...

哦哦哦好的理解了，势能面扫描是我对结构里的两个氢键距离进行了刚性扫描，但是结果显示能量最低的点O-H的距离是1.79左右，但是我实际上优化的结果是1.59左右，我不知道怎么解释这种现象，所以试了下限制性优化