I am trying to develop phase coexistence curves for generic polymers. Of course, I can use standard Flory theory free energy functionals to do this, but I want to define my nearest neighbor polymer-polymer and polymer-solvent interaction energies explicitly, and run numerical experiments.
I am thinking about running Monte Carlo simulations of polymers on a discrete lattice. My question is, does such software exist? I want to plant a polymer on a lattice along with solvent molecules and examine the behavior of the system after $N$ moves.