<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>