24小时热门版块排行榜    

查看: 1805  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 大佬帮着分析今年的评审意见?该课题连续申请了三年了,明年要不要换个方向。 100+3 cai-11138 2026-08-31 11/550 2026-09-06 20:49 by archvillain
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见