%nproc=4
%mem=10GB
%chk=1.chk
#p opt freq b3lyp/genecp <<<<<<<关键词部分 b3lyp泛函 全电子/赝势混合基组
<<<<<<<<<<一定要空一行
Title Card Required
<<<<<<<<<<一定要空一行
0 2 <<<<<<<电荷和多重度
O 1.23641237 1.02449231 0.48357982
O 3.79340942 -0.41974574 0.40393821
C 1.28546627 -0.30626188 0.36753494
N 2.59517275 -1.04941707 -0.02178549
H 2.57396972 -1.96339548 0.38342297
C 0.00854381 -1.12716700 0.62673335
H 0.27531773 -2.09151263 1.00590241
H -0.60292444 -0.61901871 1.34282824
Fe 2.53412498 0.45179272 -1.25580654
C -0.77227654 -1.29148902 -0.69042946
H -1.65948890 -1.86185816 -0.51033706
H -1.03905046 -0.32714339 -1.06959852
H -0.16080829 -1.79963731 -1.40652435
<<<<<<<<<<一定要空一行
Fe 0 <<<<<<<先定义电子基组部分 最后我一般加一个0
SDD
****
C H O N 0
6-311g(d,p)
**** <<<<<<<<<<每个不同电子基组之间要用****作为分隔符,写完最后一个电子基组也要写
<<<<<<<<<<一定要空一行
Fe 0 <<<<<<<<<<<这里开始是赝势部分 最后也要写个0
SDD
<<<<<<<<<这里一定要有个空行 ,
这个原因有好多吧,你的计算路径有没有中文字符或路径过长,你试试在桌面建文件夹再投任务看看能不能成功。
从你上传的输出文件可以看出来这个是SCF不收敛。
1 你首先需要修改你的 基组。你对所有原子都使用 Lanl2dz基组这本身极不合理。并且你这样写等于没有给过渡金属原子定义赝势。当然scf没法收敛。你需要修改成混合基组写法,对于那些小原子用 6-31g(d) , 对过渡金属用lanl2dz. 并且要正确定义赝势才行,要不然即使强行收敛的结果也是根本算出来就不对。
如果不想用赝势基组,也不想用混合基组,那也可以把lanl2dz 换成 def2svp 也可以解决基组问题。
2 其次可以加上scf=xqc 帮助收敛。但是前提是把基组改好以后。
给你个例子,例如想让Fe是SDD基组,其他是6-311g(d,p),则输入文件为
%nproc=4
%mem=10GB
%chk=1.chk
#p opt freq b3lyp/genecp <<<<<<<关键词部分 b3lyp泛函 全电子/赝势混合基组
<<<<<<<<<<一定要空一行
Title Card Required
<<<<<<<<<<一定要空一行
0 2 <<<<<<<电荷和多重度
O 1.23641237 1.02449231 0.48357982
O 3.79340942 -0.41974574 0.40393821
C 1.28546627 -0.30626188 0.36753494
N 2.59517275 -1.04941707 -0.02178549
H 2.57396972 -1.96339548 0.38342297
C 0.00854381 -1.12716700 0.62673335
H 0.27531773 -2.09151263 1.00590241
H -0.60292444 -0.61901871 1.34282824
Fe 2.53412498 0.45179272 -1.25580654
C -0.77227654 -1.29148902 -0.69042946
H -1.65948890 -1.86185816 -0.51033706
H -1.03905046 -0.32714339 -1.06959852
H -0.16080829 -1.79963731 -1.40652435
<<<<<<<<<<一定要空一行
Fe 0 <<<<<<<先定义电子基组部分 最后我一般加一个0
SDD
****
C H O N 0
6-311g(d,p)
**** <<<<<<<<<<每个不同电子基组之间要用****作为分隔符,写完最后一个电子基组也要写
<<<<<<<<<<一定要空一行
Fe 0 <<<<<<<<<<<这里开始是赝势部分 最后也要写个0
SDD
<<<<<<<<<这里一定要有个空行
,
不同基组算出来能垒会有些许不同的。至于基组是干什么用的,如果这个还都不清楚,建议去随便找本量化书把薛定谔方程大概了解一下。然后找本计算化学的书把HF-SCF过程了解一下,就知道基组本身是干啥的了。