| 查看: 1470 | 回复: 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 |
» 猜你喜欢
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有0人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有290人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有63人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
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
qzhaosdu
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1401.6
- 散金: 2052
- 红花: 7
- 帖子: 1969
- 在线: 263.5小时
- 虫号: 493351
- 注册: 2008-01-10
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
6楼2009-06-08 23:04:20
rainandstar
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 315.5
- 红花: 2
- 帖子: 91
- 在线: 69.8小时
- 虫号: 1385353
- 注册: 2011-08-31
- 专业: 质谱分析
7楼2011-09-06 10:06:54
