Modeling oxidised polypyrrole in the condensed phase with a novel force field

A novel model potential is developed for simulating oxidised oligopyrroles in condensed phases. The force field is a coarse grained model that represents the pyrrole monomers as planar rigid bodies with fixed charge and dipole moment and the chlorine dopants as point atomic charges. The analytic function contains 17 adjustable parameters that are initially fitted on a database of small structures calculated within all-electron density functional theory. A subsequent potential function refinement is pursued with a battery of condensed phase isothermal-isobaric Metropolis Monte Carlo in-silico simulations at ambient conditions with the goal of implementing a hybrid parametrization protocol enabling agreement with experimentally known thermodynamic properties of oxidised polypyrrole. The condensed system is composed of oligomers containing 12 monomers with a 1:3 dopant-to-monomer concentration. The final set of force field optimised parameters yields an equilibrium density of the condensed system at ambient conditions in excellent agreement with oxidised polypyrrole samples synthesised in wet-laboratories.