[交流] 【资源】求助一个关于单点能的问题

Originally posted by zhx2134 at 2010-12-30 09:26:38:
向各位诚恳的请教一个问题,我现在用CCSD(T)/6-311+G(d,p)在算一个分子的单点能,这个分子的结构是用B3LYP/6-31G(d)优化过的,因为CCSD(T)不能做freq,而我又想进行零点能的修正,请问我应该怎么做呢?请大家帮忙解决一下.

Originally posted by sculhf at 2010-12-30 10:01:53:

CCSD(T)不能做freq?
谁说的？
你试过吗

Originally posted by zhx2134 at 2010-12-30 10:15:02:

试过了,如果加上freq的话,交在服务器上的作业就不能算,自动停下来.如果不加freq只算sp的话就可以算的.我用的是gaussian 03

Originally posted by sculhf at 2010-12-30 10:19:42:

或许你体系太大了

Originally posted by zhx2134 at 2010-12-30 10:23:22:

谢谢回复，体系不大，就是SiCl4的一个结构

Originally posted by zhx2134 at 2010-12-30 09:26:38:
向各位诚恳的请教一个问题,我现在用CCSD(T)/6-311+G(d,p)在算一个分子的单点能,这个分子的结构是用B3LYP/6-31G(d)优化过的,因为CCSD(T)不能做freq,而我又想进行零点能的修正,请问我应该怎么做呢?请大家帮忙解决一下.

Originally posted by sculhf at 2010-12-30 10:34:57:

%chk=NaCl2-
%mem=40MW
# opt freq CCSD(T)/6-311+g* geom=connectivity

NaCl2

-1 1
Na
Cl                 1              B1
Cl                 1              B2    2              A1

  ...

Originally posted by zhx2134 at 2010-12-30 12:33:03:

非常感谢能这么耐心的回复我的帖子，我把我的输入文件帖上来，是这样的
%chk=SiH4.chk
%mem=2GB
%nproc=4
#p CCSD(T)/6-311++G(d,p) sp freq test

Title Card Required

0 1
Si                 0.0 ...

Originally posted by lingznpmg at 2010-12-30 22:09:40:
有好多文献上都是先用低基组（如B3LYP/6-31G(d)）优化并作频率分析得到零点能校正数值，然后用高基组（如CCSD(T)/6-311+G(d,p)）算一个分子的单点能,然后加上用低基组算的零点能来校正的。楼主可以找类似的文献看 ...