24小时热门版块排行榜    

查看: 1441  |  回复: 7

zzl7337

金虫 (正式写手)

[交流] 【求助】用cube文件画HOMO和LUMO轨道 已有3人参与

我用生成.cube文件的方式画HOMO和LUMO轨道,计算正常结束,也生成了.cube文件,但是为什么导入.cube文件的时候出错,.cube文件很小,才10几K
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+2):鼓励应助!! 2010-05-11 15:38:31
.cube文件很小,才10几K ,这肯定错了!如果是80*80*80的格点,应有6676kB。

你的计算cube的操作是怎样进行的?贴出来。
2楼2010-05-11 14:39:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+2): 2010-05-11 15:38:23
heyo_123(金币+2):鼓励应助!! 2010-05-11 15:38:46
现在作cube,其实可以这样作:

1、在#加关键词output=wfn;

2、将得到的波函数文件wfn读入Multiwfn程序,选择你要做的HOMO的MO序号,选择做格点(grid),选择做波函数,输入格点数,如80,80,80。然后对问格点起点和步长回答n,即由程序自动给出,就可得到HOMO的cube文件了。同样再次操作,可以得到LUMO的。

3、剩下就是作图了。
3楼2010-05-11 15:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzl7337

金虫 (正式写手)

我的.cube文件是:
%mem=200MW
%nproc=4
#n sp rb3lyp/6-31g* 5d cube=(orbitals) test

HOMO and LUMO in default cube

0 1
   6             0       -0.608687   -3.665801   -0.635413
   6             0       -0.662189   -1.281138   -0.585523
   6             0        0.687175   -1.265334   -0.064312
   6             0        1.408704   -2.427413    0.165055
………………

S-SBTD.cube

483
484
486

有错误吗?大家帮帮忙,就差一个分子了
4楼2010-05-11 15:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heyo_123

至尊木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yzcluster(金币+1): 2010-05-12 17:56:21
按照gaussview的方法,做HOMO.LUMO没问题。
5楼2010-05-11 15:39:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzl7337

金虫 (正式写手)

那为什么我的.cube文件有错误啊?
6楼2010-05-11 15:45:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+2):鼓励应助!! 2010-05-11 17:30:41
我模拟了你的计算,发现:

你的输入:

S-SBTD.cube

483
484

S-SBTD.cube之后要紧跟MO序号,你空了一格,它读了一个“空”,便只能输出它默认的一套参数,没有MO数据,正常结束,只有10几K。

正确例如:

#  POP=FULL  cube=(Orbitals)

H2

0,1
1  0.0 0.0 0.0
1  0.735396 0.0 0.0

H2-H-L.cube
1
2
7楼2010-05-11 17:25:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzl7337

金虫 (正式写手)

非常感谢,我算一下看看
8楼2010-05-11 17:33:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zzl7337 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 01:37 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:16 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:41 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zICmwzsBXjbN 2026-08-29 5/250 2026-08-30 20:08 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 G6APbkg8SA6w 2026-08-29 6/300 2026-08-30 19:02 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:51 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:40 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-30 17:56 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见