I am using the HOOMD-blue package, and I was advised to use force-shifted Leanard Jones potential, with a tree neighbor list. I do not know which parameters I should write for those charged molecules (hard spheres) "'N, charge'" The system would also have some counterions, and they have an impact on the system.
classhoomd.md.pair.ForceShiftedLJ(nlist, r_cut=None, r_on=0.0, mode='none')
nlist = hoomd.md.nlist.tree(
r_buff=0.4, check_period=1, d_max=None, dist_check=True, name=None)
Thanks in advance!