| 查看: 1468 | 回复: 5 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
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 |
» 猜你喜欢
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
关于合成MIL-125的求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有260人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
» 本主题相关商家推荐: (我也要在这里推广)
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
立方区域大小问题
★ ★ ★ ★ ★ ★ ★ ★ ★
yjcmwgk(金币+5,VIP+0):极为感谢！ 5-16 16:37
lei0736(金币+4,VIP+0):谢谢 5-17 15:13
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 ] |
5楼2009-05-16 16:21:15
2楼2009-05-16 15:44:09
erylingjet
铁杆木虫 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 5.372
- 金币: 9933.4
- 散金: 752
- 红花: 2
- 帖子: 2248
- 在线: 241.5小时
- 虫号: 415302
- 注册: 2007-06-28
- 性别: GG
- 专业: 计算化学
3楼2009-05-16 15:44:57
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
4楼2009-05-16 16:04:05
28