计算化学公社

标题: 关于混合基组的问题 [打印本页]

作者
Author: panger 时间: 2015-4-29 01:25
标题: 关于混合基组的问题
做一个体系涉及到两种元素需要添加赝势，按照下面的格式输入:
%chk=name
# 方法/genecp Opt

Title Card Required

0 1
分子系统说明部分

基组
****
M1
基组
****
M2
基组
****

M1
赝势

M2
赝势

各种空行，星号基本都试过了，无论如何第二种元素的赝势gaussian也读不出来？
大家有试过这种情况吗？该怎么办呢？谢谢。

作者
Author: sobereva 时间: 2015-4-29 04:22
M1
赝势

M2
赝势

二者之间不要有空行。

作者
Author: 小范范1989 时间: 2015-4-29 08:24
在sob的这个帖子中有下面一行：http://sobereva.com/60
对于不同元素可以使用不同赝势，在赝势定义部分写多种即可，例如对镧采用MWB46，对镥采用MWB60，应写为
C H O N 0
6-31G**
****
La
MWB46
****
Lu
MWB60
****
[空行]
La 0
MWB46
Lu 0
MWB60
[空行]
[空行]

可见，两个赝势应该没有空格。
顺带问一句，那三个四个五个都不空行？
谢谢

作者
Author: blueyangliu 时间: 2015-4-29 08:43

小范范1989 发表于 2015-4-29 08:24
在sob的这个帖子中有下面一行：http://sobereva.com/60

有出错信息吗？不会是基组不支持这种元素吧

作者
Author: panger 时间: 2015-4-29 11:23
谢谢大家，实话说上面的方法（格式）我都试过了，无论如何就是报错。不知道是不是跟我昨天用的几台机器有关，今天在实验室换了另外的机器以后任务就能交了。另外我试过了
M1
赝势

M2
赝势
两者之间有空行和没空行的情况均能提交，也能正常结束。
谢谢大家。

作者
Author: sobereva 时间: 2015-4-29 17:15

panger 发表于 2015-4-29 11:23
谢谢大家，实话说上面的方法（格式）我都试过了，无论如何就是报错。不知道是不是跟我昨天用的几台机器有关 ...

M2前面一定别留空行。留了空行没准就不读取空行后面的赝势了，即便正常计算了，结果也是离谱的。
高斯对空行非常敏感

作者
Author: panger 时间: 2015-4-29 23:33

sobereva 发表于 2015-4-29 17:15
M2前面一定别留空行。留了空行没准就不读取空行后面的赝势了，即便正常计算了，结果也是离谱的。
高斯 ...

谢谢sob老师。

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