|
|
G16,b3lyp,def2svp,gd3,PCM溶剂化，125个原子，540个电子，过渡态一次优化非常顺利结束。
从过渡态虚频震动方向的结构作为初始结构计算中间体INT1，开始优化，脚本设置6小时杀掉。看结果，在第20+个点的时候两个力YES，RMS位移YES，Maximum位移NO，结构是想要的INT1结构，虽然后续的点YES全消失了，但是不影响开心的心情；
然而之后心态崩溃开始了：
1.取该三个YES的结构做一个输入文件A，加上maxstep=5,notrust,优化6小时，能量基本不变，四个参数的YES却全部消失；
2.继续把最后一个结构做成输入文件B，继续优化6小时，仍然没有一个YES；
3.输入文件A中加上int=ultrafine,做成输入文件C，优化12小时，能量稳定，没有YES；
4.输入文件A中加上opt=gdiis，calcfc做成输入文件D，优化6小时，能量稳定，没有YES；
5.取三个YES的结构并且不改变命令，继续优化24小时，消失的YES并未出现；
5.取已经优化完成的INT2的结构（该结构不到12小时就opt+freq结束），将H换成Cl，做成INT1的初始结构文件，优化12小时，两个力参数YES，但是不收敛，重复步骤1的maxstep=5,notrust，YES仍然消失；
6.取过渡态柔性扫描曲线最低点结构优化，得到两个力参数的YES后消失，再也没有出现YES；
7.作为有强迫症的人，努力说服自己，用自己向来都不喜欢的opt=loose，成功收敛，一看频率，崩溃，一个-700+的大虚频！
从在老师那里学计算开始也几年了，第一次遇到这么难优化的结构。现在已经走上最后一条路了，opt=tight，打算拼48小时的机时，如果这次再优化不好，就努力压制自己的强迫症，彻底放弃这个结构了
|
|
发表于 Post on 2024-10-24 10:43:03
|
只看该作者 Only view this author
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author