insert rattle or shake algorithm to fix bond length during MD simulation.