I found an error in my topology, I was also changing one bond length from 0.135 to 0.128 nm. Forget about my previous messages, now I obtain exactly 0 as it should be. Stefano