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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】静电势计算
101/1返回列表
查看: 1345  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

spring965

金虫 (正式写手)

[交流] 【求助】静电势计算

哪位高手知道静电势计算需要哪些关键词，不知道我的计算过程错在哪了，


下面是输入和输出文件，哪位高手给分析分析。

chk=Gfcube                              
#p  b3lyp/3-21g  5d cube=(potential,read)  
                                          
opt                                       
                                          
-1 1  

输出
Writing cube to file .                        
End of file reading cube parameters.           
Error termination via Lnk1e in C:\G03W\l602.exe
回复此楼

» 猜你喜欢
1楼 2008-06-20 16:36:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)
★ ★ ★
小红豆(金币+3,VIP+0):辛苦了。
看样子好像是:
%chk=Gfcube.chk                              
#p  b3lyp/3-21g  5d cube=(potential) Guess=Read Geom=Check  
                                          
opt                                       
                                          
-1 1
不知你的Gfcube是不是chk文件？
赞一下(10人) 回复此楼
悠悠天地，奈何独立苍茫？
2楼2008-06-20 17:31:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spring965

金虫 (正式写手)
%chk=Gfcube.chk                              
#p  b3lyp/3-21g  5d cube=(potential)
                                          
cube                                       
                                          
-1 1
6             0        0.000000    0.682382    0.000000
6             0        0.702701   -0.724862    0.000000
7             0       -0.204764   -1.895182    0.000000
8             0        0.782573    1.663192    0.000000
8             0       -1.261098    0.650934    0.000000
1             0        1.363720   -0.790337    0.876813
1             0        1.363720   -0.790337   -0.876813
1             0       -0.841049   -1.705586   -0.780125
1             0       -0.841049   -1.705586    0.780125

但是还是报错

您看有时间在帮忙看看
赞一下 回复此楼
3楼2008-06-20 22:20:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ferlich

禁虫 (小有名气)
★
小红豆(金币+1,VIP+0):辛苦了。
本帖内容被屏蔽
4楼2008-06-20 22:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spring965

金虫 (正式写手)
★ ★ ★
小红豆(金币+3,VIP+0):辛苦了。
先使用formchk工具生成格式化的fchk

formchk jobname.chk jobname.fchk

这一步如何做到？

再用cubegen生成cub文件

cubegen 0 potential=scf jobname.fchk test.cube 80

这一步如何做到？

cubegen 0 potential=scf jobname.fchk test.cube 80 }是什么意思？

请不吝赐教
赞一下(10人) 回复此楼
5楼2008-06-20 22:52:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spring965

金虫 (正式写手)
是不是无论是opt还是freq还是其它的计算得到的*.chk文件，只要scratch目录下有*.chk文件，都可以在view中得到分子的静电势呢？
赞一下 回复此楼
6楼2008-06-20 22:59:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

北京一叶567

木虫 (小有名气)
★
lei0736(金币+1,VIP+0):谢谢 十大“耐思”版主评选欢迎参与投票 http://emuch.net/bbs/viewthread.php?tid=845527&fpage=1
不是的。只有与静电势相关的计算，最后才能通过chk生成cub文件。如mk等。
赞一下(10人) 回复此楼
7楼2008-06-21 00:03:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spring965

金虫 (正式写手)
%chk=Gfcube.chk                              
#p  b3lyp/3-21g  5d cube=(potential)
                                          
cube                                       
                                          
-1 1
6             0        0.000000    0.682382    0.000000
6             0        0.702701   -0.724862    0.000000
7             0       -0.204764   -1.895182    0.000000
8             0        0.782573    1.663192    0.000000
8             0       -1.261098    0.650934    0.000000
1             0        1.363720   -0.790337    0.876813
1             0        1.363720   -0.790337   -0.876813
1             0       -0.841049   -1.705586   -0.780125
1             0       -0.841049   -1.705586    0.780125

这个输入文件到底错在哪里了，正确的应该如何输入？
赞一下 回复此楼
8楼2008-06-21 08:13:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spring965

金虫 (正式写手)
这个体系很小，2分钟就出错
602死掉，报错信息

                Electrostatic Grid Calculations

--------------------------------------------------------------
Grid  1 is a regular grid
Map a grid with origin  -3.891498 -4.796782 -3.446400
    Take  80 steps of     0.100723  0.000000  0.000000
    Take  90 steps of     0.000000  0.100723  0.000000
    Take  68 steps of     0.000000  0.000000  0.100723
Compute potential derivative range 0 to 0 and Write to file
赞一下 回复此楼
9楼2008-06-21 09:25:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxjbuilder

木虫 (知名作家)
★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢
你的那个chk我以为是已经算出来的呢，我是想通过调用chk文件直接生成out文件，
不过好像你不该取过来坐标，应该是Guess=Read，Geom=Check试试
当然，你也可以试试用二楼的方法试试，那个肯定对。在windows下的时候，要在MSDOS中运行cubugen操作
当然也可以这样算

%chk=name.chk                              
#p  b3lyp/3-21g   cube=potential test
                                          
cube                                       
                                          
-1 1
6             0        0.000000    0.682382    0.000000
6             0        0.702701   -0.724862    0.000000
7             0       -0.204764   -1.895182    0.000000
8             0        0.782573    1.663192    0.000000
8             0       -1.261098    0.650934    0.000000
1             0        1.363720   -0.790337    0.876813
1             0        1.363720   -0.790337   -0.876813
1             0       -0.841049   -1.705586   -0.780125
1             0       -0.841049   -1.705586    0.780125

eof
我换了一个例子试了一下，可以算的
赞一下(10人) 回复此楼
悠悠天地，奈何独立苍茫？
10楼2008-06-21 17:52:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 spring965 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭