计算化学公社

 找回密码 Forget password
 注册 Register
Views: 567|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 求助gaussian优化,因力小于收敛限制100倍造成opt和freq都有no

[复制链接 Copy URL]

6

帖子

0

威望

97

eV
积分
103

Level 2 能力者

本帖最后由 shashengfo 于 2025-5-23 20:32 编辑

求助各位老师,gaussian16c02优化此分子由于力小于收敛限制100倍就直接收敛了,虽然没有虚频,但造成opt和freq都有no。Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法 - 思想家公社的门口:量子化学·分子模拟·二次元就是sob老师这篇博文里3.1说的那种情况,我不知道后面应该怎么处理了,有看到Gaussian做结构优化时遇到受力小于收敛标准100倍而收敛但位移还未达到收敛标准的问题 - 量子化学 (Quantum Chemistry) - 计算化学公社这篇帖子,但是不知道我这个分子能不能按柔性大分子那样不考虑位移收敛。希望能得到各位老师的指点,谢谢!

首先用#P PBE1PBE/def2SVP EmpiricalDispersion=GD3BJ freq opt=recalc=8优化后得到结果是opt四个yes,freq两个no。
Maximum Force            0.000034     0.000450     YES
RMS     Force            0.000004     0.000300     YES
Maximum Displacement     0.000490     0.001800     YES
RMS     Displacement     0.000099     0.001200     YES


Maximum Force            0.000034     0.000450     YES
RMS     Force            0.000004     0.000300     YES
Maximum Displacement     0.017973     0.001800     NO
RMS     Displacement     0.004174     0.001200     NO

再此结构上使用#p opt=(GDIIS,Recalc=5) freq def2svp empiricaldispersion=gd3bj int=superfine pbe1pbe 继续计算得到的结果是opt和freq显示最大位移是no且超出收敛限两倍,数值一样。(也试了calcall结果差不多)

Maximum Force            0.000001     0.000450     YES
RMS     Force            0.000000     0.000300     YES
Maximum Displacement     0.005438     0.001800     NO
RMS     Displacement     0.000542     0.001200     YES


202505232032388406..png (150.69 KB, 下载次数 Times of downloads: 15)

202505232032388406..png

opt.gjf

4.87 KB, 下载次数 Times of downloads: 1

reopt.gjf

4 KB, 下载次数 Times of downloads: 3

6

帖子

0

威望

97

eV
积分
103

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-5-23 21:54:50 | 只看该作者 Only view this author
chands 发表于 2025-5-23 21:24
你的分子柔性很大,势能面平缓,做结构优化未必能做到四个YES。现在做的已经很好了,注意点柔性分子的势能 ...

感谢老师指教!

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

2#
发表于 Post on 2025-5-23 21:24:44 | 只看该作者 Only view this author
你的分子柔性很大,势能面平缓,做结构优化未必能做到四个YES。现在做的已经很好了,注意点柔性分子的势能面可能有多个极小值,可能有多个稳定的结构。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 21:48 , Processed in 0.166896 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list