24小时热门版块排行榜    

查看: 1530  |  回复: 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 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见