计算化学公社

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

[Gaussian/gview] gaussian自定义基组溶剂化计算报错

[复制链接 Copy URL]

10

帖子

0

威望

69

eV
积分
79

Level 2 能力者

做溶剂化计算练习 input文件如图
SCRF=SMD 能算 但是改为SCRF=(SMD,Read)  末尾加上RADII=KLAMT时  报错图片里的内容 请问是什么原因呢

我参考的是《显著改进SMD溶剂模型描述Br和I元素精度的溶剂模型SMD18的介绍》http://sobereva.com/613 这个文章 文章末尾说“SMD18模型体现出隐式溶剂模型计算能量的准确度对原子半径是多么的敏感。顺带一提,J. Phys. Chem. A, 123, 9498 (2019)提出的SMD-B模型将SMD定义的原子半径都改为了Bondi范德华半径,发现对离子体系(特别是含硫的)的溶解自由能得到一定改进。既可以如上通过modifysph来对各个元素一一修改半径,也可以在scrf里写read的同时在输入文件末尾后面空一行直接加上Radii=Bondi来实现,当然后者省事得多。”所以我想学习一下

微信图片_20220618213537.png (28.25 KB, 下载次数 Times of downloads: 20)

input文件

input文件

微信图片_20220618213842.png (23.55 KB, 下载次数 Times of downloads: 28)

报错信息

报错信息

928

帖子

1

威望

8271

eV
积分
9219

Level 6 (一方通行)

2#
发表于 Post on 2022-6-18 21:57:45 | 只看该作者 Only view this author
显然,自定义基组信息和自定义scrf信息的顺序并不是任意的,试试调整顺序就知道了

10

帖子

0

威望

69

eV
积分
79

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-6-18 22:05:19 | 只看该作者 Only view this author
hebrewsnabla 发表于 2022-6-18 21:57
显然,自定义基组信息和自定义scrf信息的顺序并不是任意的,试试调整顺序就知道了

我刚调整了两者关键词的顺序,报错结果不变,然后调整了坐标后二者信息的顺序,报错信息显示:Unrecognized atomic symbolRADII

10

帖子

0

威望

69

eV
积分
79

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2022-6-18 22:08:56 | 只看该作者 Only view this author
自定义基组是在https://www.basissetexchange.org上下载的def2-tzvpd gaussian版本 没有经过修改

928

帖子

1

威望

8271

eV
积分
9219

Level 6 (一方通行)

5#
发表于 Post on 2022-6-18 22:35:48 | 只看该作者 Only view this author
本帖最后由 hebrewsnabla 于 2022-6-18 22:42 编辑
zhangyu 发表于 2022-6-18 22:05
我刚调整了两者关键词的顺序,报错结果不变,然后调整了坐标后二者信息的顺序,报错信息显示:Unrecogniz ...

确实是基组信息应该在前面。不过你的gbs文件里面可能有空行。导致高斯把空行后的内容认为是scrf信息了

10

帖子

0

威望

69

eV
积分
79

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2022-6-19 14:02:59 | 只看该作者 Only view this author
hebrewsnabla 发表于 2022-6-18 22:35
确实是基组信息应该在前面。不过你的gbs文件里面可能有空行。导致高斯把空行后的内容认为是scrf信息了

C:\Users\11730\Desktop\微信图片_20220619135554
感谢您的建议 我翻阅gbs文件 发现确实在这个地方有个空行 报错原因应该是这个地方引起的
请问空行后面部分是赝势信息吗?如果我没用到赝势的话可以把后面这部分删除掉吗?如果说我需要使用赝势的话,空行后面这部分我该如何处理,才能避免因为没用到基组中的某些原子而报错“No RB atoms found in this molecule.”这类信息?在原子名称前加“-”可以吗?

10

帖子

0

威望

69

eV
积分
79

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2022-6-19 14:57:16 | 只看该作者 Only view this author
hebrewsnabla 发表于 2022-6-18 22:35
确实是基组信息应该在前面。不过你的gbs文件里面可能有空行。导致高斯把空行后的内容认为是scrf信息了

谢谢老师 问题找到并且已经解决了 非常非常感谢

本版积分规则 Credits rule

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

GMT+8, 2026-2-28 12:41 , Processed in 0.167632 second(s), 23 queries , Gzip On.

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