计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9785|回复 Reply: 3

[Gaussian/gview] 【分享】监督高斯结构优化进度的脚本

[复制链接 Copy URL]

375

帖子

1

威望

2378

eV
积分
2773

Level 5 (御坂)

发表于 Post on 2021-1-2 19:29:56 | 显示全部楼层 Show all |阅读模式 Reading model
本帖最后由 北大-陶豫 于 2021-5-4 20:30 编辑

https://github.com/Tao-Yu-Chem/GaussianScripts/ 里的 Converge_Reader.cpp(或最后的附件)

(2021年3月6日更新:添加了对双杂化泛函的支持,为说明添加了英文版。)(Update on March 6, 2021: Added support for double-hybrid functionals, and added an English version for the description.)

这是本人写的监督高斯结构优化进度的脚本,能显示出每步的能量、受力、位移等信息,如果受力、位移等达到了收敛标准会加上下划线,十分直观。
This is a script I wrote to monitor the progress of Gaussian structure optimization. It can display the energy, force, displacement and other information of each step. If the force and displacement reach the convergence standard, they will be underlined, which is very intuitive.

下面是一个 opt freq job 的使用效果:
Below is the use effect of an opt freq job:
基本示意图.png

(2021年3月6日新增)对双杂化泛函也支持得很好,能给出加入了 MP2 关联能的正确能量(而 GView 会给出错误的能量,即 SCF Done 的能量,详见 http://bbs.keinsci.com/thread-13450-1-1.html):(Added on March 6, 2021) It also supports double hybrid functionals very well (While GView will give the wrong energy, i. e. , the energy of "SCF Done", see http://bbs.keinsci.com/thread-13450-1-1.html for details):
opt_freq_B2PLYP.png

会自动找出受力、位移最小的步,因此可以帮助判断结构优化是否发生了振荡。比如下面就是一个发生振荡的例子:
It will automatically find the step with the smallest force and displacement, so it can help determine whether the structure optimization had an oscillation. Below is an example of oscillation:
结构优化不收敛示意图.png
当前已经走了 184 步,但能量最低的是第 147 步,受力最小的是第 85 步,位移最小的是第 36 步,因此显然发生了振荡。
At present, 184 steps have been taken, but the lowest energy is on step 147, the lowest force is on step 85, and the lowest displacement is on step 36, so it is obvious that there is an oscillation.

会自动检测 scan job,并给出 scan summary:
It can automatically detect a scan job and give a scan summary:
scan summary.png

有时 scan job 未完成,或结构优化走步太多(发生振荡、未能收敛),会给对应的能量显示警告信息:
Sometimes the scan job is not completed, or the structure optimization takes too many steps (oscillation, failure to converge), a warning message will be displayed for the corresponding energy:
scan_B2PLYP_InProgress.png
(注意,对于 scan job 中结构优化发生振荡的情况,我的脚本会给出最低能量,而不是像 GView 一样给出最后一个结构的能量。)
(Note that for the oscillation of structure optimization in scan job, my script will give the lowest energy instead of giving the energy of the last structure like GView does.)

能自动检测错误并显示出来:(但是还不太支持 IRC job)
It can also automatically detect and display errors: (but IRC jobs are poorly supported)
检测错误.png

现在已知的缺点:
Currently known disadvantages:
(1)基本只支持 sp,opt,freq 等 job type,尚不支持 IRC 等 job type。
(1) Basically, only job types such as sp, opt, freq are supported, while job types such as IRC are not yet supported.
(2)能量读取是从有 "SCF Done" 的一行中读取,因此对涉及 MPn 的任务,如 MPn、CCSD 等任务无法得到正确的能量(但已经支持双杂化泛函得到正确能量了)。
(2) The energy read is read from the line with "SCF Done", so for tasks involving MPn, such as MPn, CCSD and other tasks, the correct energy cannot be obtained (But double hybrid functional has been supported to obtain the correct energy).

Converge_Reader.cpp

6.73 KB, 下载次数 Times of downloads: 26

Converge_Reader_For_Windows.cpp

7.34 KB, 下载次数 Times of downloads: 16

评分 Rate

参与人数
Participants 7
eV +32 收起 理由
Reason
zmjsce + 4 好物!
biogon + 5
Sway + 4 とてもいい!
IreneAlder + 4 谢谢
kay + 5 赞!
函数与激情 + 5 GJ!
冰释之川 + 5 好物!

查看全部评分 View all ratings

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

1061

帖子

16

威望

5780

eV
积分
7161

Level 6 (一方通行)

計算化学の社畜

发表于 Post on 2021-1-2 19:39:45 | 显示全部楼层 Show all
本帖最后由 冰释之川 于 2021-1-2 19:52 编辑

IRC 过渡态啥的可以参考钟成老师的程序:《(更新3.0版)Gaussian优化过渡态的监控/诊断小脚本-tsmonitor 》http://bbs.keinsci.com/thread-14081-1-1.html

P.s. (1)其实用当前指标与指标的收敛限之比来跟踪几何优化收敛情况更加直观,比值小于1代表达到收敛限,(2)如果能显示最小的振动本征值就更好了。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
北大-陶豫 + 2 赞!

查看全部评分 View all ratings

Stand on the shoulders of giants

1852

帖子

1

威望

4098

eV
积分
5970

Level 6 (一方通行)

发表于 Post on 2021-1-4 13:53:16 | 显示全部楼层 Show all
molden 可视化监控
免费,可用。
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
App,  https://gitee.com/hpc4you/linux
hpc4you toolkit,  https://gitee.com/hpc4you/hpc
电邮hpc4you@163.com VX: hpc4you

375

帖子

1

威望

2378

eV
积分
2773

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2021-5-4 17:08:51 | 显示全部楼层 Show all
本帖最后由 北大-陶豫 于 2021-5-4 20:34 编辑

5月4日更新:加了个Windows版,并改了2个bug:
1 计算单点能的时候受力会显示乱码,因为没有读入力导致力的值是未定义值;

2 处理了运行时忘记传入log文件名参数的问题。

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

本版积分规则 Credits rule

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

GMT+8, 2023-2-2 23:24 , Processed in 0.203827 second(s), 25 queries .

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