| 查看: 2421 | 回复: 21 | |||||
| 当前主题已经存档。 | |||||
yjcmwgk禁虫 (文坛精英)
密度泛函·小卒
|
[交流]
【求助】关于自定义cube文件的一个疑问，20金币，谢谢诸位高手
|
||||
|
请教各位大侠： 我想计算静电势（自定义格点文件） 我所知道的方法是： ==gjf文件开始== %chk=G:\sundries\card test\c.chk # B3LYP/6-31+G* cube=(potential,cards) pot 0,1 H 0.0 0.0 0.0 H 0.735396 0.0 0.0 G:\sundries\card test\pot.cube -1 -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 ==gjf文件结束== 我现在想用cubegen命令 如何自定义格点呢？ 高斯03的说明书中关于cubegen的解释我看了，它对自定义格点的指定方法我也看懂了，跟cube一样。但是对格点的说明如何读取呢？我看貌似说是写一个格点说明文件，然后读取这个文件？ 这一点说明书中貌似没给解释啊。也许是我太笨了，呵呵 谢谢诸位大侠 [ Last edited by yjcmwgk on 2009-11-26 at 13:44 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 Gaussian | 作图技巧 |
» 猜你喜欢
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
关于合成MIL-125的求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有235人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
lihb734
铁杆木虫 (职业作家)
站在计算化学入门的门槛上
- QC强帖: 2
- 应助: 132 (高中生)
- 金币: 7495.4
- 散金: 9378
- 红花: 74
- 帖子: 3107
- 在线: 1177.3小时
- 虫号: 448403
- 注册: 2007-11-01
- 性别: GG
- 专业: 理论和计算化学
2楼2009-10-19 17:52:25
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
★ ★ ★ ★ ★ ★
yjcmwgk(金币+3,VIP+0): 10-19 20:24
lei0736(金币+3,VIP+0):谢谢 10-27 20:10
yjcmwgk(金币+3,VIP+0): 10-19 20:24
lei0736(金币+3,VIP+0):谢谢 10-27 20:10
|
你用的gjf文件自定义cube是没有错的。 小小的缺点是：定义的cube起始点取-8.0，按现在的格点数和步长，要到+8.0，前面格点数是正值80，就是18Å边长。二你算的H2所展开的potential空间的有效分布大约只需边长4～6Å已足够。构建了一个过大的空间，装了一个小图形，实际上降低了作图的精度。 再，既用cube关键词，它直接会生成cube文件，可以省去%行。 如果不用cube关键词而用cubegen命令，手册上如下说，我没有这样做过。 cubegen命令的格式如下， Cubegen memory kind fchkfile cubefile npts format 如果npts取0，按默认格点数等运行。 如果npts取-1，则在下面读入自定义格点文件，格式和你在cube关键词时一样，即： -1 -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 |
3楼2009-10-19 18:02:32
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
4楼2009-10-19 19:12:45
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
5楼2009-10-19 19:13:41
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
6楼2009-10-19 19:18:32
zhou2009
版主 (著名写手)
- QC强帖: 15
- 应助: 337 (大学生)
- 贵宾: 4.053
- 金币: 42654.6
- 红花: 153
- 帖子: 1729
- 在线: 1274.7小时
- 虫号: 676800
- 注册: 2008-12-18
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
7楼2009-10-19 19:30:36
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
8楼2009-10-19 19:42:55
yjcmwgk
禁虫 (文坛精英)
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 555.4
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
我自己已经把问题解决了。
|
问题已经解决 用cubegen，我习惯于敲一长串字符进去，所有的信息都在一行中，这样不行。 必须用以下方式才行： D:\ProgramFiles\gaussian03w>cubegen Property [Density]? potential Formatted Checkpoint file? G:\sundries\a.fchk Cube file [ ]? test2.cube Approximate points per side [0]? -1 -1 -8.0 -8.0 -8.0 80 0.2 0.2 0.2 80 0.2 0.2 0.2 80 0.2 0.2 0.2 Header in cube file [H]? h [ Last edited by yjcmwgk on 2009-10-19 at 19:50 ] |
9楼2009-10-19 19:45:08
10楼2009-10-19 20:18:27
28