计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3375|回复 Reply: 5

[其它量化程序] 【BAGEL使用心得】Gaussian、BAGEL联用做CASPT2几何优化的脚本

[复制链接 Copy URL]

197

帖子

5

威望

1689

eV
积分
1986

Level 5 (御坂)

发表于 Post on 2021-12-23 15:25:59 | 显示全部楼层 Show all |阅读模式 Reading model
本帖最后由 Freeman 于 2022-7-31 17:00 编辑

【2022.07.31更新】


原来的脚本功能单一,代码不规范,应对不同的场景还要用户自己改代码,十分不便。现在gau_bag.sh升级重写了。细节如下。 gau_bag.sh (6.73 KB, 下载次数 Times of downloads: 22)

评分 Rate

参与人数
Participants 6
威望 +1 eV +25 收起 理由
Reason
chrinide + 5 GJ!
joeson + 5 好物!
hebrewsnabla + 5 GJ!
sobereva + 1
doublezhang + 5
biogon + 5 とてもいい!

查看全部评分 View all ratings

422

帖子

7

威望

4764

eV
积分
5326

Level 6 (一方通行)

BSJ Institute

发表于 Post on 2022-2-15 16:46:20 | 显示全部楼层 Show all
感谢楼主的工作。

楼主是否有测试过用该脚本优化激发态?我在尝试使用这段脚本进行构型优化时,优化基态非常正常,但优化激发态则会在算完初始构型后退出,而生成的各输出文件内容都是正常的。以下为.sh文件中的相应内容:
  1. cat>>$base/${base}.json<<EOF
  2. },
  3. {
  4. "title":"force",
  5. "target":1,
  6. "export":"true",
  7. "method":[{
  8. "title":"casscf",
  9. "nact":6,
  10. "nclosed":27,
  11. "nstate":3,
  12. "maxiter":100
  13. }]
  14. }
复制代码

197

帖子

5

威望

1689

eV
积分
1986

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2022-2-15 17:14:19 | 显示全部楼层 Show all
Accelerator 发表于 2022-2-15 16:46
感谢楼主的工作。

楼主是否有测试过用该脚本优化激发态?我在尝试使用这段脚本进行构型优化时,优化基态 ...

95行改为awk 'NR==2{printf "%20.12e",$1}' $base/ENERGY.out > $outputfile
95行以下的所有FORCE_0改成FORCE_1
就好了

597

帖子

20

威望

4151

eV
积分
5148

Level 6 (一方通行)

发表于 Post on 2022-5-5 21:31:22 | 显示全部楼层 Show all
好象一些简单的优化算法遇到碳环就容易出问题。比如cfour,经常会在快要收敛的时候,突然把碳环上的一个配体掰到相反的方向

26

帖子

0

威望

320

eV
积分
346

Level 3 能力者

发表于 Post on 2022-5-11 01:52:48 | 显示全部楼层 Show all
写BAGEL输入文件有个技巧,如果熟悉Python和json的话,就当Python脚本来写,用IDE也可以,这样括号,逗号不容易出错。
机器学习非绝热动力学PyRAI2MD开发者
https://github.com/mlcclab/PyRAI2MD-hiam
课题组网站 https://www.x-mol.com/groups/hiam_mlcclab

197

帖子

5

威望

1689

eV
积分
1986

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2022-7-31 17:09:39 | 显示全部楼层 Show all
重写更新了一下。冒个泡。

本版积分规则 Credits rule

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

GMT+8, 2023-2-2 22:38 , Processed in 0.370310 second(s), 25 queries .

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