MECP程序搜索交叉点遇到的问题咨询，谢谢指点

小范范1989 · 发表于 Post on 2016-6-23 19:35:27

自己学习了sob老师的帖子和MECP程序：http://bbs.keinsci.com/forum.php?mod=viewthread&tid=865
我的问题如下：
1：我想寻找S1-T2的势能面交叉点。所以，我的Input——Header A和B的文件为：
%mem=2GB
%nproc=12
%chk=S1-gas.chk
#n pbe1pbe/6-311G* force guess(read) td

First State

0 1

%mem=GB
%nproc=12
%chk=T2-gas.chk
#n pbe1pbe/6-311G* force guess(read) td(triplet,root=2)

Second State

0 1

通过把prepare.sh里开头的natm改成实际原子数。然后运行./prepare.sh。运行./runfirst.sh。出现的Job0_A.log和Job0_B.log没有正常结束，末尾为：
6                   5.03728  -3.64617  -0.35
1                   6.5061 -2.5512 -1.48056
1                   3.40645  -4.46848 0.78964
1                   5.55526  -4.59437  -0.44384
End of file in ZSymb.
Error termination via Lnk1e in /home/software/g09/l101.exe at Thu Jun 23 12:01:10 2016.
Job cpu time:  0 days  0 hours  0 minutes  6.6 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=    1

我猜测是不是我输入文件的问题？因为我用例子能正常结束。
谢谢指点。

lyfchem · 发表于 Post on 2016-6-23 20:04:10

为什么要加上td呢

小范范1989 · 发表于 Post on 2016-6-23 20:27:25

lyfchem 发表于 2016-6-23 20:04
为什么要加上td呢

你好，我是想找S1和T2的，所以，S1我就写了td。
不知道对不对，谢谢你

lyfchem · 发表于 Post on 2016-6-23 21:14:04

小范范1989 发表于 2016-6-23 20:27
你好，我是想找S1和T2的，所以，S1我就写了td。
不知道对不对，谢谢你

肯定不对撒

liyuanhe211 · 发表于 Post on 2016-6-24 01:22:05

末尾加几个空行

小范范1989 · 发表于 Post on 2016-6-24 07:29:47

liyuanhe211 发表于 2016-6-24 01:22
末尾加几个空行

谢谢。我在geom末尾空了一个空行，同时在Input_Header_A和B 中加了空行。为什么还是这个错误？谢谢指点。附件是我的文件，是不是我的表头写错了？这样找S1-T2的表头有错误吗？谢谢

sobereva · 发表于 Post on 2016-6-24 16:28:30

小范范1989 发表于 2016-6-24 07:29
谢谢。我在geom末尾空了一个空行，同时在Input_Header_A和B 中加了空行。为什么还是这个错误？谢谢指点。 ...

geom后面加俩空行

高斯Linux输入文件末尾至少有俩空行是常识

遗忘爱神 · 发表于 Post on 2016-12-21 22:54:07

所以你最后找到了这个圆锥交叉点了吗？为什么我把geom里的结构提出来算单点两个态还差了挺多

杨小狗 · 发表于 Post on 2017-2-18 20:49:56

lyfchem 发表于 2016-6-23 21:14
肯定不对撒

这个到底对不对呢加不加td呢我感觉现在改版的sobMECP的extract_energy可以对理论方法进行选择应该不必了

小范范1989 · 发表于 Post on 2017-4-12 08:25:54

再次咨询sob老师一下，@sobereva
就是我用这个sob-MECP模拟test中例子，C6H5+能够正常结束，但是在模拟这个C6H5+_A2-3B1的时候，拷贝的初始的文件，计算不出例子的结果。
geom没有更新，我计算我自己的分子找S1-T1交叉点，也是这个问题。reprotfile没有内容，可能是我运行runfirst和runMECP是通过fan-one.pbs脚本运行，
输出都在fan-one.o55884和55885文件中了。
文件已经上传（我把chk文件删除了，其他都保留了），望sob老师指点一二，谢谢sob老师。

小范范1989 · 发表于 Post on 2017-4-12 09:51:27

刚刚有测试了TEST中最后一个分子，还是这样的问题呢、reportfile是0字节，可能我是脚本提交，输出在了.o文件中了。
谢谢老师指点

sobereva · 发表于 Post on 2017-4-12 10:47:52

小范范1989 发表于 2017-4-12 09:51
刚刚有测试了TEST中最后一个分子，还是这样的问题呢、reportfile是0字节，可能我是脚本提交，输出在了.o文 ...

我并不保证程序和pbs之类脚本提交方式兼容
先用单机，把自带的例子里的输入文件跑一下

小范范1989 · 发表于 Post on 2017-4-12 23:30:37

sobereva 发表于 2017-4-12 10:47
我并不保证程序和pbs之类脚本提交方式兼容
先用单机，把自带的例子里的输入文件跑一下

谢谢sob老师指点。执行老师的例子，找s0-T1的交叉点能找到。脚本应该是可以执行的。
今天在群里咨询，想用背景电荷的方式，runfirst没问题，runMECP的时候，采用这个方式遇到如下的错误，（用脚本执行的，脚本文件我也上传了）ProgFile Exists - OK
First Input OK
At line 230 of file MECP.f (unit = 8, file = 'ProgFile')
Fortran runtime error: End of file
Problem with Fortran program
Step Number 0 -- MECP not yet converged
An error has occurred, possibly in the Gaussian Job

高斯任务有错误吗？还是需要修改MECP.f的内容呢？谢谢sob老师

sobereva · 发表于 Post on 2017-4-12 23:44:30

小范范1989 发表于 2017-4-12 23:30
谢谢sob老师指点。执行老师的例子，找s0-T1的交叉点能找到。脚本应该是可以执行的。
今天在群里咨询，想 ...

原理上和用背景电荷并不矛盾
检查下高斯输出文件，看怎么报错的
input tail最后也多加一些空行试试

小强 · 发表于 Post on 2017-5-14 11:38:43

小范范1989 发表于 2017-4-12 23:30
谢谢sob老师指点。执行老师的例子，找s0-T1的交叉点能找到。脚本应该是可以执行的。
今天在群里咨询，想 ...

请问楼主，你这问题解决了吗？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] MECP程序搜索交叉点遇到的问题咨询，谢谢指点

评分 Rate

浏览过的版块