[lammps-users] lj12-6 with spherical harmonics

Dear list:
I'm attempting to simulate a highly anisotropic pes and have found a paper where the effective interactions have been fitted to a lennard jones 12-6 potential with spherical harmonics to handle the angular dependence. Before I went ahead and coded this into LAMMPS myself, I was wondering (hoping!) that someone may already have this or similar code they would be willing to share. If not, any other pointers would be greatly appreciated

Thanks in advanced

I haven't heard of anyone doing this.