计算化学公社

标题: gaussian自定义基组溶剂化计算报错 [打印本页]

作者
Author: zhangyu 时间: 2022-6-18 21:43
标题: gaussian自定义基组溶剂化计算报错
做溶剂化计算练习 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来实现，当然后者省事得多。”所以我想学习一下

作者
Author: hebrewsnabla 时间: 2022-6-18 21:57
显然，自定义基组信息和自定义scrf信息的顺序并不是任意的，试试调整顺序就知道了

作者
Author: zhangyu 时间: 2022-6-18 22:05

hebrewsnabla 发表于 2022-6-18 21:57
显然，自定义基组信息和自定义scrf信息的顺序并不是任意的，试试调整顺序就知道了

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

作者
Author: zhangyu 时间: 2022-6-18 22:08
自定义基组是在https://www.basissetexchange.org上下载的def2-tzvpd gaussian版本没有经过修改

作者
Author: hebrewsnabla 时间: 2022-6-18 22:35
本帖最后由 hebrewsnabla 于 2022-6-18 22:42 编辑

zhangyu 发表于 2022-6-18 22:05
我刚调整了两者关键词的顺序，报错结果不变，然后调整了坐标后二者信息的顺序，报错信息显示：Unrecogniz ...

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

作者
Author: zhangyu 时间: 2022-6-19 14:02

hebrewsnabla 发表于 2022-6-18 22:35
确实是基组信息应该在前面。不过你的gbs文件里面可能有空行。导致高斯把空行后的内容认为是scrf信息了

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

作者
Author: zhangyu 时间: 2022-6-19 14:57

hebrewsnabla 发表于 2022-6-18 22:35
确实是基组信息应该在前面。不过你的gbs文件里面可能有空行。导致高斯把空行后的内容认为是scrf信息了

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)