Thank you for your response. It did not occur to me that by scaling B by a factor of Chi would solve the problem. (That would have been too easy ). Actually I modified LAMMPS to include ‘Chi’ as an input parameter. Now as I read your mail my efforts look unnecessary. I see your point.