%nproc=16
%mem=10GB
%chk=1.chk
#p opt freq b3lyp/genecp <<<<<<<关键词部分 b3lyp泛函 全电子/赝势混合基组
<<<<<<<<<<一定要空一行
Title Card Required
<<<<<<<<<<一定要空一行
0 2 <<<<<<<电荷和多重度,之下是分子坐标 这个坐标可以把你的分子用gview画出来 然后让gview给出个坐标直接复制过来用就行。
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
<<<<<<<<<这里一定要有个空行 ,
你要知道的具体是什么呢?是要问算重金属和轻原子的不同么?那基本上最大的不同就是混合基组设置 还有多重度不太容易确定的问题。其他的在input写法上和算一般的轻原子区别不太大的
给你个例子,例如算一个带Fe的分子 的结构优化频率计算 如果想让Fe是SDD基组,其他是6-311g(d,p),则输入文件为
%nproc=16
%mem=10GB
%chk=1.chk
#p opt freq b3lyp/genecp <<<<<<<关键词部分 b3lyp泛函 全电子/赝势混合基组
<<<<<<<<<<一定要空一行
Title Card Required
<<<<<<<<<<一定要空一行
0 2 <<<<<<<电荷和多重度,之下是分子坐标 这个坐标可以把你的分子用gview画出来 然后让gview给出个坐标直接复制过来用就行。
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
<<<<<<<<<这里一定要有个空行
,
谢谢大佬