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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Multiwfn » 【求助】关于cubegen的立方区域大小问题
61/1返回列表
查看: 1241  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

[交流] 【求助】关于cubegen的立方区域大小问题

各位Gaussian高手：
最近计算了一个分子轨道，是反键轨道。在linux下使用的命令如下：
cubegen 0 mo=29 xxx.fchk yyy.cube 100 h
但是发现分子轨道图的大小超过了立方区域的大小。如图（isovalue=0.02）
Gaussian说明文件中讲到可以通过手工建立输入文件的方法来得到自定义的立方区域，但是语焉不详，不知道如何建立输入文件，并将自己建立的立方格子输入文件加入到cubegen命令中去。此处请教高人指教，非常感谢！

又及：增大isovalue可以解决问题，但是我不想这样做。我希望的是不改变isovalue的前提下，如何通过增大立方格子的大小来得到完整的轨道图？谢谢各位高手！！

A value of -1 says to read the cube specification from the input stream, according to the following format:
IFlag, X0, Y0, Z0 Output unit number and initial point.
N1, X1, Y1, Z1 Number of points and step-size in the X-direction.
N2, X2, Y2, Z2 Number of points and step-size in the Y-direction.
N3, X3, Y3, Z3 Number of points and step-size in the Z-direction.

又又及：
下面的是高斯03使用手册中的话语。我能看懂，但是不会用：
IFlag is the output unit number. If IFlag is less than 0, then a formatted file will be produced; otherwise, an unformatted file will be written.
If N1<0 the input cube coordinates are assumed to be in Bohr, otherwise, they are interpreted as Angstroms. |N1| is used as the number of X-direction points in any case; N2 and N3 specify the number of points in the Y and Z directions, respectively. Note that the three axes are used exactly as specified; they are not orthogonalized, so the grid need not be rectangular.
The value -5 says to read in an arbitrary list of points from standard input. If you enter this input by hand, terminate the input with an end-of-file (i.e., ^D under Unix). Alternatively, you can redirect standard input to a file containing the list of points (do not place a blank line or ^D at the end of the file).

[ Last edited by yjcmwgk on 2009-11-26 at 13:45 ]
回复此楼

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

Gaussian

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-05-16 11:37:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)
★
lei0736(金币+1,VIP+0):谢谢 5-17 15:13
GaussView吗？不是可以右键拖动来缩放吗
赞一下(10人) 回复此楼
2楼2009-05-16 15:44:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)
★
lei0736(金币+1,VIP+0):谢谢 5-17 15:13
这么麻烦啊？


不知道这样做和打开CHK文件做分子轨道有什么区别啊？
赞一下(10人) 回复此楼
3楼2009-05-16 15:44:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
回复recoli，我说的是cube文件的立方格子的大小，不是GView中分子显示比例的大小
回复erylingjet，这样做与你在GView中打开fchk或chk文件，然后直接点击MO产生的分子轨道是一样的，区别在于这样可调的参数比较多
赞一下 回复此楼
4楼2009-05-16 16:04:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

立方区域大小问题

★ ★ ★ ★ ★ ★ ★ ★ ★
yjcmwgk(金币+5,VIP+0):极为感谢！ 5-16 16:37
lei0736(金币+4,VIP+0):谢谢 5-17 15:13
我没有在linux下使用的经验，我在WINDOWS下遇到这种情况是这样解决的。
   当选定了格点总数之后，cubegen会自动根据所计算的分子空间大小，确定起始点和步长，通常都是可以把MO都装得进去的。
   但对你算的反键轨道，MO会有超常的弥散，空间选定就不够用了，需要使用关键字cube=(Orbitals,cards)，来指定cube的格点数、初始点、步长。如：
# HF/6-311++G** POP=FULL IOP(6/7=3) cube=(Orbitals,cards)

H2-LUMO

0,1
1  0.0 0.0 0.0
1  0.735396 0.0 0.0

H2-MO2-CUB.cube
-51 -8.0 -8.0 -8.0
80 0.200 0.000 0.000
80 0.000 0.200 0.000
80 0.000 0.000 0.200
2

    这是一个H2分子反键轨道的例子。
   在成键时，起点取-2就可以了，现在取到了-8，才能将图装下。
   linux也应有相应的操作吧。

[ Last edited by zhou2009 on 2009-5-16 at 16:36 ]
赞一下(11人) 回复此楼
5楼2009-05-16 16:21:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qzhaosdu

金虫 (著名写手)
gwdavid(金币+0,VIP+0):此话怎么讲啊,兄弟 6-9 09:22
齐球球，鄙视你
回复此楼
Anewday,anewhour,anewminute,anewpeople.
6楼2009-06-08 23:04:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rainandstar

银虫 (小有名气)
帖子真精彩!
已经收录到淘贴专辑《Gaussian》
赞一下 回复此楼
7楼2011-09-06 10:06:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yjcmwgk 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭