<div dir="ltr"><div>And some modification to the gradients calculation is needed.</div><div><br></div><div>    ! compute the gradients</div><div>    gt1 =  ( s32 * ism * ism )  * SIGN(1.0_dp , DOT_PRODUCT ( t12, tn)) * tm</div><div>    gt4 = -( s32 * isn * isn )  * SIGN(1.0_dp , DOT_PRODUCT ( t12, tn)) * tn</div><br>On Friday, August 14, 2015 at 11:28:30 AM UTC-4, 姚懿 wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">In the mol_force.F file. there is one line in force_imp_torsions subroutine<div><br><div>phi = SIGN( ACOS ( cosphi ), DOT_PRODUCT ( t12, tn ))<br></div></div><div><br></div><div>I suggest to change it to </div><div><br></div><div>phi = ACOS( cosphi )</div><div><br></div><div>since improper potential is not periodic. For a dihedral close to 180 degree it can varies between -180.0 and 180.0 that will influence the result a lot.</div><div><br></div><div>YY</div></div></blockquote></div>