| 查看: 561 | 回复: 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),如何设置。请高手指教！ |
回复此楼» 猜你喜欢
分子动力学 Plumed教程求助
已经有0人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有141人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
2楼2017-04-06 19:06:55
zxy964777869
铜虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 9.5
- 散金: 2255
- 红花: 8
- 帖子: 501
- 在线: 900.4小时
- 虫号: 1355217
- 注册: 2011-07-25
- 性别: GG
- 专业: 理论和计算化学
|
可以去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
40