Patchy particle simulation to prepare coarse-grain simulations of microgels

I have read this Paper about preparing a coarse-grain simulation of microgel particles from patchy particles and I want to replicate this procedure for other microgel systems.

They were using patchy particles (bi- and tetravalent) to set up the structure with custom interactions between the patches. When the structure was self-assembled they replace the patchy particles with a bead-spring model.

I have limited knowledge about LAMMPS, but I worked with Gromacs before. How can I implement patchy particles with custom forces and patch-patch interactions (Eq. 1-4 from the Paper) in LAMMPS?

This is the potential they have used, already implemented in LAMMPS:

pair_style cosine/squared command

For the same patch-patch interactions:

bond_style fene command

1 Like