Development of hybrid coarse-grained atomistic models for rapid assessment of local structuring of polymeric semiconductors

Decades of work in the field of computational study of semiconducting polymers using atomistic models illustrate the challenges of generating equilibrated models for this class of materials. While adopting a coarse-grained model can be helpful, the process of developing a suitable model is particularly non-trivial and time-consuming for semiconducting polymers due to a large number of different interactions with some having an anisotropic nature. This work introduces a procedure for the rapid generation of a hybrid model for semiconducting polymers where atoms of secondary importance (those in the alkyl side chains) are transformed into coarse-grained beads to reduce the computational cost of generating an equilibrated structure. The parameters are determined from easy-to-equilibrate simulations of very short oligomers and the model is constructed to enable a very simple back-mapping procedure to reconstruct geometries with atomistic resolution. The model is illustrated for three related polymers containing DPP (diketopyrrolopyrrole) to evaluate the transferability of the potential across different families of polymers. The accuracy of the model, determined by comparison with the results of fully equilibrated simulations of the same material before and after back-mapping, is fully satisfactory for two out of the three cases considered. We noticed that accuracy can be determined very early in the workflow so that it is easy to assess when the deployment of this method is advantageous. The hybrid representation can be used to evaluate directly the electronic properties of structures sampled by the simulations.

Automated summary

Decades of work have shown the challenges of generating equilibrated models for semiconducting polymers, especially due to the anisotropic nature of their interactions. A rapid hybrid model generation procedure is introduced in this study, which proves to be accurate for two out of three cases considered. The hybrid representation can be used to directly evaluate the electronic properties of structures sampled by simulations.