wbn 发表于 2016-12-19 10:47
谢谢SOB老师。那我就从现在的轨迹力提取一个或几个相邻的分子，用DFT优化以后再来调试一下力场。另外我想 ...

动力学程序里管优化叫minimization

sobereva 发表于 2016-12-19 09:14
比起用rdf峰位置和DFT比，应该用力场优化的结构和DFT优化的结构对比。如果把氢的范德华参数去掉可以得到明 ...

谢谢SOB老师。那我就从现在的轨迹力提取一个或几个相邻的分子，用DFT优化以后再来调试一下力场。另外我想问下在gromacs里怎么优化结构啊...和energy minimization一样吗？我用来energy minimization的mdp文件是这样的行不行？

cpp                 =  /usr/bin/cpp
define              =
constraints         =  none
integrator          =  steep
nsteps              =  100000
;
;       Energy minimizing stuff
;
emtol               =  10
emstep              =  0.01

nstcomm             =  1
ns_type             =  grid
rlist               =  1.5
rcoulomb            =  1.5
rvdw                =  1.5
Tcoupl              =  no
Pcoupl              =  no
gen_vel             =  no

jackie 发表于 2016-12-19 09:07
个人感觉哈，需要看下OPLS-AA力场文献中当初对于氢键类型的拟合。力场当初是只考虑类似N-H...O=C型氢键还是 ...

谢谢。OPLS-AA在创建的时候好像确实没有对C-H ...O 氢键进行专门的测试

比起用rdf峰位置和DFT比，应该用力场优化的结构和DFT优化的结构对比。如果把氢的范德华参数去掉可以得到明显更合理的结果就去掉

关于O----H-C，可以那么判断

个人感觉哈，需要看下OPLS-AA力场文献中当初对于氢键类型的拟合。力场当初是只考虑类似N-H...O=C型氢键还是也包含了C-H...O=C型氢键？是只针对中性型体系来拟合还是也考虑了带电体系？而且你的咪唑本身的结构也有改变，可能拟合的参数没有考虑你这样的体系？而对于C-H是一个弱极化键，形成的C-H...O=C型中性体系氢键属于较弱型的氢键，经验从头算计算出的氢键键长一般在2.5 A左右变化。而你算出的2.4 A 正好也接近这个的值，所以还应需要搞清楚opls力场当初对于氢键的细节是怎么处理的。