24小时热门版块排行榜    

查看: 494  |  回复: 2

lixiaona158

木虫 (著名写手)

自由人

[求助] BOMD里面phase如何设置计算的是reverse

# hf/lanl1dz bomd=(maxpoints=1000,rtemp=300,phase=(8,10)) iop(5/13=1)

过度态是在8号和10号原子之间震动,如上所设是使8和10间距离拉远(forward),但是现在想让他们靠近(reverse),如何设置。请高手指教!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whoisme

铁虫 (正式写手)

请问楼主,这个问题解决了吗?我现在也遇到这个问题了,但是基本没见过有人讨论,是因为BOMD不靠谱?还是怎么回事
2楼2017-04-06 19:06:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxy964777869

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by whoisme at 2017-04-06 19:06:55
请问楼主,这个问题解决了吗?我现在也遇到这个问题了,但是基本没见过有人讨论,是因为BOMD不靠谱?还是怎么回事

可以去log文件中找MWVelocity或这Velocity速度,所有的乘以-1,再次构建输入文件,加上ReadMWVelocity或者readVelocity就变成reverse了.
下面是一个DA反应的reverse例子,供参考(这个里面是读check里的构型,但无关紧要)
# geom=check bomd(stepsize=10000, maxpoint=300,ReadMWVelocity,Rtemp=298) iop(1/8=10000) iop(5/13=1) scf(maxcyc=100)  6-31g(d) m062x

Title Card Required

0 1

0

  3.132557010225E+13   -6.911830025237E+13   -3.619079660215E+13
  2.184843855621E+13  2.346101919885E+13   -7.192318300385E+12
  7.236554150938E+12  4.518401467820E+13   -1.033604779883E+14
-5.901432369073E+13   -2.873581935557E+13  2.634628443915E+13
-2.462377574252E+13  1.962558701538E+13  1.294434156125E+13
  -2.402221576176E+13   -6.294839639573E+13   -1.076133684978E+13
  -2.764690649314E+13   -7.641353769502E+13   -4.006480513063E+12
  -4.330065394631E+13  1.030069652573E+14  2.967441957725E+13
6.055399605237E+13  2.376695319430E+12  7.528180549885E+13
  -4.461406762243E+12  4.393605446418E+13   -1.099335646164E+13
  -1.908798493371E+13  1.110358503041E+14  8.110635738800E+13
4.173341984029E+13   -1.001721394708E+14  5.520888573981E+13
4.523747331444E+12  3.391164704303E+13   -7.099431624973E+12
9.252794058984E+13   -5.214144150282E+12   -8.523698667319E+12
4.519872839516E+13  1.172208942949E+14  9.379250217085E+13
1.300040448118E+13  2.059225057387E+13  1.034724982639E+14
3楼2017-04-06 22:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lixiaona158 的主题更新
信息提示
请填处理意见