24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1337  |  回复: 4

zhficcas

银虫 (正式写手)

[求助] 混合基组中,如何读入前一基组的CHK文件?

在有些时候的计算中,对分子优化可以先用小基组优化得到合理的初始构型,然后通过Guess和Geom关键词读入小基组的CHK文件,接着进行更大基组的优化.但对混合基组如何读取呢?采用GUESS=READ和GEOM=ALLCHECK就不能用了.举个例子,首先采用下面的命令优化构型:
%mem=1200MB
%rwf=文件名
%chk=文件.chk

分子构型
原子1
lanl2dz
****
原子2,原子3 0
3-21G*
****

原子1 0
lanl2dz
优化完成后,打算对原子2和原子3 采用更高基组6-31g*来优化, 这时如何从上面的CHK文件读取原来的信息呢.采用:
geom=allcheck guess=read 就不行了.

[ Last edited by zhficcas on 2011-8-7 at 18:21 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lxzhq

铁杆木虫 (著名写手)

不知啊
2楼2011-08-07 18:40:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖


gmy1990(金币+1): 2011-08-08 11:03:36
zhficcas(金币+1): 多谢交流 2011-08-08 11:22:13
其实,我没看出来为什么guess=read 和geom=allcheck会不行...
3楼2011-08-08 10:36:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhficcas

银虫 (正式写手)

引用回帖:
3楼: Originally posted by fatpig8832 at 2011-08-08 10:36:27:
其实,我没看出来为什么guess=read 和geom=allcheck会不行...

举个例子,在命令行输入:
#p geom=allcheck guess=read chkgen b3lyp/genecp OPT 后在分子坐标栏输入:
空格
原子1 0
LANL2DZ
****
原子2,原子3 0
6-31G*****
空格
原子1 0
LANL2DZ
运行后会提示和原来的基组不匹配(原来用3-21g*的),没有混合基组时,我们可以直接把6-31g*基组写在命令栏里(b3lyp/6-31g*)的,这样就能运行的。现在混合基组后,基组要写在分子坐标栏里,这样就运行不了了。
4楼2011-08-08 11:35:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

zhficcas(金币+1): 愕! 我无语...., :o,多谢 2011-08-08 16:05:33
引用回帖:
4楼: Originally posted by zhficcas at 2011-08-08 11:35:57:
举个例子,在命令行输入:
#p geom=allcheck guess=read chkgen b3lyp/genecp OPT 后在分子坐标栏输入:
空格
原子1 0
LANL2DZ
****
原子2,原子3 0
6-31G*****
空格
原子1 0
LANL2DZ  ...

这...难道你没想过为什么会冲突...

因为chkgen是指从chk中读取gen,即3-21g*,你后面又指定了6-31g*,显然不匹配...
所以只要把chkgen去掉就好了...
5楼2011-08-08 15:00:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhficcas 的主题更新
信息提示
请填处理意见