Hello users,
I have a question concerning the pair style gran/hertz/history. By comparing the normal force defined in LAMMPS with the cited paper by Zhang,Makse, I get difference of sqrt(2) coming essentially from the polydisperse factor. I believe it should be given by sqrt(2R_iR_j/(R_i+R_j)) instead of sqrt(R_iR_j/(R_i+R_j)) as written in the LAMMPS doc. This of course is just an overall factor which can be easily accounted for by adjusting the values of all the parameters. But then, looking at page 477, it says that for monodisperse systems, the parameters should be multiplied by 2. I believe they should rather be multiplied by sqrt(2). Is that right?
François Fillion-Gourdeau, PhD
Scientifique de la défense | Defence Scientist
Section Protection et Effets d'armes | Weapons Effects and Protection Section
Recherche et développement pour la Défense Canada - Valcartier | Defence Research and Development Canada - Valcartier
2459, Boul. Pie-XI Nord | 2459, Pie-XI Blvd. North
Québec, Quebec Canada G3J 1X5
Téléphone | Telephone (418) 844-4000 x4395
Télécopieur | Facsimile (418) 844-4876
Gouvernement du Canada | Government of Canada
www.valcartier.drdc-rddc.gc.ca <http://www.valcartier.drdc-rddc.gc.ca/>