24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2796  |  回复: 23

萱月琪琪

金虫 (正式写手)

[求助] 在gaussian计算中,先用小基组计算,在改成大基组计算出现的问题

在gaussian计算中,先用小基组b3lyp/6-31+g(d,p)计算opt freq,选出最小稳定构型之后,再用b3lyp/6-311++g(d,p) 计算该构型下的opt freq,输入的命令如下,我怎么感觉我这是从头开始计算了呢。怎样输入命令是利用b3lyp/6-31+g(d,p)计算的结果进行的b3lyp/6-311++g(d,p) 呢?感谢您的帮助
(1)基组b3lyp/6-31+g(d,p)计算opt freq输入文件是:
%chk=GDC4ofq4-1.chk
%mem=600MW
%nproc=6
#p opt freq b3lyp/6-31+g(d,p) geom=connectivity

GDC4ofq4-1

1 1
C
C                  1              B1


(2)基组b3lyp/6-311++g(d,p)计算opt freq输入文件是:
%chk= GDC4ofq4-1.chk
%mem=600MW
%nproc=6
#p opt freq b3lyp/6-311++g(d,p) geom=connectivity

GDC4ofq4-1

1 1
C
C                  1              B1
怎样输入是在6-31基组的基础之上进行的计算呢,我这样输入是从新开始计算的,计算的时间特别久。输入文件我采用的是win下的gaussview软件,通过ssh上传到linux下计算的
先表示感谢了

[ Last edited by 萱月琪琪 on 2013-3-31 at 11:27 ]
回复此楼

» 收录本帖的淘帖专辑推荐

Gaussian

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

世上无难事只怕有心人。。。。信心。。。毅力。。。努力。。。成功
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xyz7819227

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
等小基组算完了直接读取chk文件,再采用大基组算不就行了吗
海纳百川,有容乃大
8楼2013-03-31 16:59:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lxhlxh052c

木虫 (著名写手)

小木虫之无虫

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
萱月琪琪: 金币+5, 有帮助, 谢谢哈 还是没有计算成功 2013-04-01 00:12:28
把geom=connectivity 改成geom=allcheck,然后把title 及以下的部分删去就可以了吧。
这就是我想要的
2楼2013-03-31 11:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stalart

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
建议楼主把郭勇编辑的那本入门好好看看,自己就能解决很多问题。
3楼2013-03-31 13:39:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

萱月琪琪

金虫 (正式写手)

引用回帖:
3楼: Originally posted by stalart at 2013-03-31 13:39:15
建议楼主把郭勇编辑的那本入门好好看看,自己就能解决很多问题。

我看过了呀,也没看出怎么回事!说明书也看来这,就是没弄明白呢
世上无难事只怕有心人。。。。信心。。。毅力。。。努力。。。成功
4楼2013-03-31 16:05:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

萱月琪琪

金虫 (正式写手)

引用回帖:
3楼: Originally posted by stalart at 2013-03-31 13:39:15
建议楼主把郭勇编辑的那本入门好好看看,自己就能解决很多问题。

我之前看的是66页的那版,刚找到了459页的那版了
世上无难事只怕有心人。。。。信心。。。毅力。。。努力。。。成功
5楼2013-03-31 16:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

萱月琪琪

金虫 (正式写手)

引用回帖:
2楼: Originally posted by lxhlxh052c at 2013-03-31 11:56:19
把geom=connectivity 改成geom=allcheck,然后把title 及以下的部分删去就可以了吧。

是这个意思不?
使用gv输入文件
%chk=GDC4ofq5.chk
%mem=60MW
%nproc=6
#p opt freq b3lyp/6-311++g(d,p) geom=allcheck

31时 %chk=GDC4of5.chk  ,311时%chk=GDC4ofq5.chk   
手动删除下面的题目及以下  gaussian计算的时候不到两秒钟就出现的是die link
世上无难事只怕有心人。。。。信心。。。毅力。。。努力。。。成功
6楼2013-03-31 16:35:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maynard

金虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
萱月琪琪: 金币+15, ★★★很有帮助, 你的意思是额外输入关键词RCFC ,说实在的我还是不怎么明白。chk文件的功能是啥我清楚,就是不知道具体怎么操作。感谢了 2013-04-01 00:15:58
如果你想得到b3lyp/6-311++g(d,p)下优化的构型,那应该没办法减少计算量了,geom=allcheck也只是读取结构,手动提取差不多
geom=connectivity是确定原子键结关系的,就是那些原子间成键,哪些不成键,如果没有特殊需要就没必要写这个
当然是有个方法,就是在opt里加上RCFC,从chk文件里提取小基组计算的力矩阵,再在大基组条件下优化,至于计算量能下降多少就不好说了
7楼2013-03-31 16:39:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

萱月琪琪

金虫 (正式写手)

引用回帖:
8楼: Originally posted by xyz7819227 at 2013-03-31 16:59:40
等小基组算完了直接读取chk文件,再采用大基组算不就行了吗

我这每次使用gv打开.chk的时候都显示的是 Can't find Gaussian formchk utility.  那我要怎么办呢?
世上无难事只怕有心人。。。。信心。。。毅力。。。努力。。。成功
9楼2013-03-31 17:06:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xyz7819227

木虫 (正式写手)

【答案】应助回帖

引用回帖:
9楼: Originally posted by 萱月琪琪 at 2013-03-31 17:06:55
我这每次使用gv打开.chk的时候都显示的是 Can't find Gaussian formchk utility.  那我要怎么办呢?...

你把chk文件剪切出来放在桌面上试试
海纳百川,有容乃大
10楼2013-03-31 17:25:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 萱月琪琪 的主题更新
信息提示
请填处理意见