I'm thinking of running MD simulation of poling process (heating and cooling in external electric field) of host-guest polymer system: organic chromophore in polymer matrix (PMMA or polycarbonate).
So, I'm trying to find how to build the system. I've found that Packmol can be used to randomly place prebuilt polymer chains but I don't understand how can I introduce my chromophore (I also need to control polymer density and chromophore weight content).
What tools can you recommend for building such systems and feeding them to some MD code (at the moment I'm considering GROMACS, LAMMPS or CP2k, but I'm not sure)?