计算化学公社

标题: 求助chemshell调用gaussian时，qm区域混合基组的设定 [打印本页]

作者
Author: chengjunhui 时间: 2020-4-11 21:34
标题: 求助chemshell调用gaussian时，qm区域混合基组的设定
老师们好
我chemshell调用gaussian时，qm区域混合基组的设置总是有问题，查了一下手册，上面例子说basisspec={ { sto-3g all } {3-21g {c* h*} } }，但我按照这个格式设， basisspec={ { sto-3g all } {3-21g {c* h*} } }，总是提示错误get_basis : could not assign basis sets for atom symbols: x x。希望老师能告诉我哪里不对

作者
Author: chenjinfeng850 时间: 2020-4-12 08:02
根据你的输入文件来看，首先确定一下单一基组是否能正常运行。

然后根据ChemShell 对Gaussian 的调用基组关键词为：
sto3g, 321g, 431g, 631g, 631gd, 631gdp, lanl2dz, sddall, dz, and dzp; respectively set the basis sets: STO-3G, 3-21G, 4-31G, 6-31G, 6-31G(d), 6-31G(d,p), LANL2DZ, SDDALL, D95, and D95(d,p).
（http://www.cse.scitech.ac.uk/ccg ... anual/gaussian.html）
测试下是否因为关键词写的不标准，导致ChemShell 不认

作者
Author: chengjunhui 时间: 2020-4-12 17:17

chenjinfeng850 发表于 2020-4-12 08:02
根据你的输入文件来看，首先确定一下单一基组是否能正常运行。

然后根据ChemShell 对Gaussian 的调用 ...

老师您好，我试了一下发现，对于我的体系即使不用混合基组，从431g开始就开始出现基组问题，下面是gaussian输出文件所提示的错误，这是什么原因呢，谢谢老师
Standard basis: 4-31G (6D, 7F)
Atomic number out of range in LEEAO.
Error termination via Lnk1e in /opt/g03c02/g03/l301.exe at Sun Apr 12 07:26:53 2020.

作者
Author: hebrewsnabla 时间: 2020-4-12 17:42

chengjunhui 发表于 2020-4-12 17:17
老师您好，我试了一下发现，对于我的体系即使不用混合基组，从431g开始就开始出现基组问题，下面是gaussi ...

4-31G仅支持H-Ne的元素。超过这个范围当然报错

作者
Author: chengjunhui 时间: 2020-4-12 17:50

hebrewsnabla 发表于 2020-4-12 17:42
4-31G仅支持H-Ne的元素。超过这个范围当然报错

明白了，非常感谢！

作者
Author: chengjunhui 时间: 2020-4-13 10:24

chenjinfeng850 发表于 2020-4-12 08:02
根据你的输入文件来看，首先确定一下单一基组是否能正常运行。

然后根据ChemShell 对Gaussian 的调用 ...

老师您好，我现在已经可以进行普通的混合基组的计算了，但我现在要用混合赝势基组， basisspec= { { 6-31g(d) all } { sddall Fe } }，总是提示错误Fatal error : incomplete basis set，似乎是因为我的格式写的不对，它识别不了，老师您能告诉我格式应该是什么样的嘛，谢谢老师！

作者
Author: chenjinfeng850 时间: 2020-4-13 16:04

chengjunhui 发表于 2020-4-13 10:24
老师您好，我现在已经可以进行普通的混合基组的计算了，但我现在要用混合赝势基组， basisspec= { { 6-31 ...

具体的报错信息是什么？
1.你试试铁用LANL2DZ基组可不可
2.以把铁改成 fe 试试
再有就是ChemShell关于自定义机组有一句话：
A basis specification of "null" can be used to denote that no basis set is required (this is needed since it is otherwise a fatal error to assign an incomplete basis set).
（http://www.cse.scitech.ac.uk/ccg ... l/manual/basis.html）
不知道这句话能不能帮到你

作者
Author: chengjunhui 时间: 2020-4-13 22:21

chenjinfeng850 发表于 2020-4-13 16:04
具体的报错信息是什么？
1.你试试铁用LANL2DZ基组可不可
2.以把铁改成 fe 试试

太谢谢您了，是fe的问题

作者
Author: 找雷 时间: 2020-8-18 17:57
老师，能否分享一下chemshell程序？

作者
Author: zjxitcc 时间: 2020-8-19 11:01

找雷发表于 2020-8-18 17:57
老师，能否分享一下chemshell程序？

别把这儿当小mu虫呀

最新python版本的ChemShell是免费开源的，无需在这儿求，自己去官网下载即可https://www.chemshell.org/licence

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