计算化学公社

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

[Gaussian/gview] 关于G16中G4方法bug的一个解决办法

[复制链接 Copy URL]

2479

帖子

11

威望

6864

eV
积分
9563

Level 6 (一方通行)

关于G16中G4方法bug的一个解决办法
文/我本是个娃娃

摘要:在G16(A02, B01)中,使用G4方法计算存在不收敛的情况,使得整体计算结果错误。而G09中,又无法实现G16独有的G4(SP)和G4(noopt)功能。据此,本文提出一个解决这个bug的方案。经过测试表明,方案不仅可以解决这个问题,并使得最终的计算结果符合标准,并且可以灵活修改,完成Gaussian程序内置G4方法所无法完成的计算任务。

G4.md (7.58 KB, 下载次数 Times of downloads: 37)

G4.zip (68 KB, 下载次数 Times of downloads: 37)


评分 Rate

参与人数
Participants 2
eV +11 收起 理由
Reason
sobereva + 8
anson + 3 好物!

查看全部评分 View all ratings

356

帖子

0

威望

1919

eV
积分
2275

Level 5 (御坂)

2#
发表于 Post on 2022-11-14 09:06:50 | 只看该作者 Only view this author
为什么一篇文章只有标题作者和摘要,没有正文

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

3#
发表于 Post on 2022-11-14 09:57:20 | 只看该作者 Only view this author
mfdsrax2 发表于 2022-11-14 09:06
为什么一篇文章只有标题作者和摘要,没有正文

正文当然是最关键的脚本与安装包啦
Stand on the shoulders of giants

477

帖子

0

威望

887

eV
积分
1364

Level 4 (黑子)

4#
发表于 Post on 2024-4-11 09:23:59 | 只看该作者 Only view this author
请教一下老师,我利用“ #p ub3lyp 6-31G(2df,p) guess=mix stable=opt nosymm”和“ #p G4=SP guess=read geom=allcheck nosymm”计算双自由基的单点能,out文件末尾如图,能量结果是异常的。out文件中只有最后一个SCF的<S**2>= 0.0000,其余都是接近1的,这是不是您提到的“G4方法计算存在不收敛的情况,使得整体计算结果错误”的原因呢?或者说明G4不适合计算这个双自由基的单点能呢?谢谢老师

202404110923583912..png (57.87 KB, 下载次数 Times of downloads: 10)

202404110923583912..png

2479

帖子

11

威望

6864

eV
积分
9563

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2024-4-15 10:56:09 | 只看该作者 Only view this author
糖糖糖豆9988 发表于 2024-4-11 09:23
请教一下老师,我利用“ #p ub3lyp 6-31G(2df,p) guess=mix stable=opt nosymm”和“ #p G4=SP guess=read  ...

软件是用的G09还是G16?
如果是G16,建议改成G09试试。

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 22:41 , Processed in 0.198808 second(s), 24 queries , Gzip On.

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