Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法

wzkchem5 · 发表于 Post on 2024-3-15 23:33:21

Sui1999 发表于 2024-3-15 02:46
老师您好，我在做一个由CHONP元素组成的体系，使用orca501优化结构，关键词是! b3lyp/G cc-pVTZ d3 opt fre ...

可以试试关掉RIJCOSX。
对于杂化泛函，ORCA默认打开RIJCOSX，其中的RI部分需要辅助基组，如果用户没有给定辅助基组，ORCA会默认对非相对论基组用def2/J，对相对论基组用SARC/J。但是def2/J恐怕未必适合cc-pVTZ。高斯没有类似RIJCOSX的方法，因此不会引入这种近似，但与此同时也会算得更慢。
当然正如sob老师所说，cc-pVTZ不太适合给DFT，尤其是非双杂化泛函用，一定要用triple zeta基组优化的话，不如用def-TZVP、def2-TZVP(-f)或def2-TZVP。这样一来也就自动解决了辅助基组的问题，因为def2/J完全适合搭配这些基组使用

qingmo · 发表于 Post on 2024-6-1 20:53:43

本帖最后由 qingmo 于 2024-6-1 22:00 编辑

老师您好！我对一个分子结构做常规优化后，再做freq计算。没有虚频，但其中的最大位移是NO。因此我取最后结构又做了opt＝calcall计算，查看log文件结果是没有虚频，四个也都是YES。
之后我打开opt＝calcall计算后的log文件取出最后结构再做freq计算，虽然还是没虚频，但关于位移的两项却变成NO了。。。
请问老师这是什么原因，calcall不是和freq一样都是精确计算Hessian计算吗

sobereva · 发表于 Post on 2024-6-2 09:29:00

qingmo 发表于 2024-6-1 20:53
老师您好！我对一个分子结构做常规优化后，再做freq计算。没有虚频，但其中的最大位移是NO。因此我取最后结 ...

优化过程计算位移的算法不仅依赖于当前步的信息，还依赖于之间一些步。而且你若要100%精确基于opt最后的结构算，应该用geom=check从chk里读取全精度的坐标。

qingmo · 发表于 Post on 2024-6-2 16:26:57

本帖最后由 qingmo 于 2024-6-2 16:29 编辑

sobereva 发表于 2024-6-2 09:29
优化过程计算位移的算法不仅依赖于当前步的信息，还依赖于之间一些步。而且你若要100%精确基于opt最后的 ...

感谢老师您的答复！已改成geom＝cheak提取chk坐标计算。
老师您说的第一句话“依赖于之间一些步”，我不是太懂。。请问可以理解为此时体系结构和上一步优化的结构之间的关系吗？
另外，geom计算是读取了全精度坐标，如果把全精度坐标提取出来应怎么做？我在log和fchk文件里都没找到全精度的。。之前都是GV里打开out文件计算，确实菜鸟行为了

qingmo · 发表于 Post on 2024-6-3 11:41:58

本帖最后由 qingmo 于 2024-6-3 12:49 编辑

按照老师您所说的用geom=check从chk里读取全精度的坐标步骤后，freq得到的确实如opt=calcall步骤一样，没有出现虚频和NO了，看来应该就是坐标问题。
之后用GV观察到fchk文件坐标小数点后有效数字达八位，而原来out文件只有六位，因此复制该坐标做freq计算，此时结果又出现了NO。
不知道是因为GV观察到的fchk文件坐标其实也并不是精确坐标？还是因为即使是，但gaussian要想读取到精确的，只能根据geom=check从chk步骤才能读取到？

（然后分别用out文件坐标，fchk文件坐标，和geom=check读取chk文件计算激发能，激发能和振子强度基本一致。是不是只要没有虚频，就算位移相关项出现了NO，关于算跃迁偶极矩，激发能这类物理量，此时的结构也是可以接受的？）

lurensan · 发表于 Post on 2024-6-3 12:04:53

各位老师，请问这个怎样把虚频消掉，尝试了很多方法都不好使

irccw-4-5-o-quetong-6-31.7z (96.7 KB, 下载次数 Times of downloads: 1)

sobereva · 发表于 Post on 2024-6-3 23:51:19

lurensan 发表于 2024-6-3 12:04
各位老师，请问这个怎样把虚频消掉，尝试了很多方法都不好使

本帖里该说的都说了，没有可补充的

sobereva · 发表于 Post on 2024-6-3 23:52:51

qingmo 发表于 2024-6-3 11:41
按照老师您所说的用geom=check从chk里读取全精度的坐标步骤后，freq得到的确实如opt=calcall步骤一样，没有 ...

gview载入、保存时的有效位数问题。gjf里的数据记录精度比chk里的低

qingmo · 发表于 Post on 2024-6-4 13:32:14

sobereva 发表于 2024-6-3 23:52
gview载入、保存时的有效位数问题。gjf里的数据记录精度比chk里的低

好的，谢谢老师

blue11 · 发表于 Post on 2024-9-5 10:13:22

sobereva 发表于 2016-3-6 02:17
只要比较相互的能量、求差值，积分格点必须完全一致，int=ultrafine要加就都加。

其它分子倒不一定 ...

请问sob老师，如果A和B优化过程中，适用的坐标不一致，需要该成一致的坐标吗？
如果A的不同构型A1和A2使用的坐标也不一致，需要放在一样的坐标下进行比较吗？
还是只需要一样的积分格点就可以了？谢谢sob老师！

sobereva · 发表于 Post on 2024-9-5 10:49:17

blue11 发表于 2024-9-5 10:13
请问sob老师，如果A和B优化过程中，适用的坐标不一致，需要该成一致的坐标吗？
如果A的不同构型A1和A2使 ...

看不懂
A和B都不是一个物质，不可能坐标一致

blue11 · 发表于 Post on 2024-9-5 19:34:42

sobereva 发表于 2024-9-5 10:49
看不懂
A和B都不是一个物质，不可能坐标一致

不好意思，sob老师，没说明白我的意思。我想问一下，我的A在高斯默认的冗余内坐标下优化成功了一个构型，但是在笛卡尔坐标下优化成功了另一个构型。请问一下，我可以计算一下这两个构型的的能量，然后选出一个能量最稳定的吗？

sobereva · 发表于 Post on 2024-9-6 00:05:54

blue11 发表于 2024-9-5 19:34
不好意思，sob老师，没说明白我的意思。我想问一下，我的A在高斯默认的冗余内坐标下优化成功了一个构型， ...

如果你是指在两种坐标系下优化得到了不同构型，应当用能量更低的（确切来说是自由能更低的）。但注意未必这就是真正能量最低结构，有可能需要用molclus做构象/构型搜索来确定真正的能量最低结构。

虎王 · 发表于 Post on 2024-10-8 18:55:10

请问一下，用了opt(gdiis)后可以优化收敛，但log出现了Warning! Use of Opt=GDIIS is deprecated since it is seldom a good choice. 是否有什么影响？例如用了opt(gdiis)是否会导致优化构型不精确？

zjxitcc · 发表于 Post on 2024-10-8 18:58:24

本帖最后由 zjxitcc 于 2024-10-8 18:59 编辑

虎王发表于 2024-10-8 18:55
请问一下，用了opt(gdiis)后可以优化收敛，但log出现了Warning! Use of Opt=GDIIS is deprecated since it ...

如果正常结束，没有影响，不会不精确。可以google了解一下优化方法的基本原理，便不会有此困惑。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法