| 查看: 1463 | 回复: 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 |
» 猜你喜欢
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
厦门大学侯旭团队招聘青年教师、博士后及博硕士研究生
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有200人回复
-
哈工大深圳 表界面物理化学方向（胶体与界面，光谱探测研究） 博士招生
已经有0人回复
-
元素硫的化合价和氧化态
已经有1人回复
-
江苏科技大学2026年招聘启事
已经有1人回复
-
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
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
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
20