应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】用cube文件画HOMO和LUMO轨道
81/1返回列表
查看: 1103  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zzl7337

金虫 (正式写手)

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

我用生成.cube文件的方式画HOMO和LUMO轨道，计算正常结束，也生成了.cube文件，但是为什么导入.cube文件的时候出错，.cube文件很小，才10几K
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-05-11 12:48:26
已阅   回复此楼   关注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人) 回复此楼
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人) 回复此楼
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没问题。
赞一下(2人) 回复此楼
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
赞一下(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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭