24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2023  |  回复: 12
当前主题已经存档。

xiaowandouer

木虫 (正式写手)

[交流] 【求助】高斯中的BOMD

最近想试着用高斯中的BOMD算点动力学,把手册中例子的基组换成了自定义基组(因为自己的体系会用到自定义基组),可是总报L118错
WANTED A FLOATING POINT NUMBER AS INPUT. (有时候是wanted an integer)
  FOUND A STRING AS INPUT.
c o h 0                                                                           
                                                                                 
                                                                                 
                 
?
Error termination via Lnk1e in /home/gauss/g03/l118.exe at Wed Dec 17 15:38:52 2008.


难道BOMD和自定义基组无法同时使用?

麻烦知道的虫友指点一下,谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

xiaowandouer

木虫 (正式写手)

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 8-10 12:10
to erylingjet: 这些是对应命令行设置的参数。

仔细考虑了下上面列出的问题,其实问题不在自定义基组上,只要把六项终止标准全部列上(还是期望能像手册中提到的那样,只列一两项,可是现在没搞懂),就可以跑。如果跑不出预想的构型,加个动能就可以了。morse参数似乎不是必须的。

另外上面说的神奇问题其实不神奇,因为对H来说3-21g和6-31g*不会引起基函数个数的变化,只是高斯函数会发生变化,这个自己动手算下就好。昨天自己有点犯傻,呵呵

还是希望能和做BOMD的虫友多多交流,上面只是自己一点肤浅的认识
9楼2008-12-18 16:12:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 8-10 12:10
现在看来如果想让输入文件跑起来,只要仔细阅读手册上关于BOMD的介绍就好,所得结果的好赖,后续的分析等等这方面的经验我也很缺乏,很希望有经验的虫友能来讲讲

补充:我前几天比较迷惑的终止标准问题,现在稍有进展。就像手册上说的那样你可以任意设置六项中的一项或者若干项,只是把不设置的项写上0(0表示关闭这些终止标准),见手册原话“a value of zero for any parameter turns off testing for the corresponding stopping criteria.”
11楼2008-12-19 09:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yyx19840628

木虫 (著名写手)

鹰的眼睛,狼的耳朵,豹的速 ...

我也想做呢?苦于没有头绪,楼主你会了教我。
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
2楼2008-12-17 18:08:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)

★ ★
lei0736(金币+2,VIP+0):谢谢
个人认为,输入文件有问题

我算其它的也出现过同样的问题

不知道你的输入文件是什么样的?
3楼2008-12-17 19:18:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quantumor

金虫 (著名写手)

快乐兔子


lei0736(金币+1,VIP+0):谢谢
估计是你定义的基组的格式不符号Gaussian的输入文件的要求。
愿好运与快乐伴随你!
4楼2008-12-17 19:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

崩溃,整理了半个多小时的问题,就因为网络一个错误给关了,全没了!!!!
只好明天再来写了
5楼2008-12-17 22:11:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 8-10 12:09
最近这两天刚开始接触BOMD,计算的经验是一点都没有,所以问题一大堆,恳请有相关计算经验的虫友们多多指点。,本人会用金币答谢(虽然热心的虫友并不是冲着金币来的,O(∩_∩)O…………)下面上问题

**************************************************************
手册上说,在大多数情况下,终止标准取一两项就足够了,利用手册上给的例子,我做了下测验,发现完全不是这么回事。
1.
# 输入文件与例子相同,只是按照从后往前的顺序,每次去掉一个标准或者n<6个标准,可以算;

# 如果是自定义基组则不能算(这儿可以确定不是自定义基组格式的问题,后面会提到);

# 如果任意去掉1-5中的任意一个标准,也不能算;

我的问题是:出现上述情况的原因是什么呢?什么样的情况(也就是例子说的大多数情况)下,终止标准只需要一两项就可以?什么样的情况下需要输入更多项?如果我只输入R1,R4两项标准,那么程序能否给识别出来?也就是说终止标准的输入次序打乱后,程序能否识别?
***********************************************************************************

2.关于自定义基组
#与例子用的基组相同,只是将其写成自定义的形式,可以算;
#把基组改变,C,O定义为3-21g基组;H定义为6-31g*基组,可以算,但是奇怪的问题产生了,检查两个输

出文件,竟然发现基函数的个数相同!!!很神奇!!!(可以排除自定义基组格式的问题,因为如果不

对,程序不会往下跑);
#用自定义基组跑成功的输入文件,去掉终止标准的任何一项都不能跑。

问题:自定义基组跟终止标准间有什么关系?

附上我的自定义基组的输入文件

# hf gen BOMD(Phase=(1,3),RTemp=300,NSample=1,ReadStop) Geom=Crowd

xxxx

0 1
C
O 1 r1
H 1 r2 2 a
H 1 r3 3 b 2 180.
     
r1 1.15275608
r2 1.74415774
r3 1.09413376
a 114.81897892
b 49.08562961
     
1
1 1 2 2
13.0 11.0 1.3 2.5 0.0000005
1 5.145
C O -112.09329898 1.12895435 0.49458169 2.24078955
H H -1.12295984 0.73482237 0.19500473 1.94603924

c o  h 0
3-21g
****

基本就是在这个文件上修修改改
**********************************************************************************

3. 初始动能如何确定?
4. Morse参数如何确定?

[ Last edited by xiaowandouer on 2008-12-18 at 12:05 ]
6楼2008-12-18 09:30:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

还有一个问题是:初始构型是否必须从过渡态开始?
7楼2008-12-18 09:43:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎再次光临计算模拟版!
1
1 1 2 2
13.0 11.0 1.3 2.5 0.0000005
1 5.145

以上这些是输入文件自动产生的吧?去掉试试。
8楼2008-12-18 13:35:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spkeey

金虫 (正式写手)

标记。我也很想做BOMD,苦于自己很愚钝阿,LZ学会了可以开个贴讲讲。
SPKEEY
10楼2008-12-18 21:25:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowandouer 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见