Molecular modeling of crystalline alkylthiophene oligomers and polymers.

We present the results of a thorough molecular modeling study of several alkylthiophene-based oligomers and polymers. In particular, we consider two polymers whose limit-ordered crystal structures have been recently reported by our group, on the basis of powder X-ray data analysis: poly(3-(S)-2-methylbutylthiophene) (P3MBT) and form I' of poly(3-butylthiophene) (P3BT). We first describe the development of a series general purpose force fields for the simulation of these and related systems. The force fields incorporate the results of ab initio calculations of the bond torsion energies of selected oligomers and differ in the set of atomic charges used to represent the electrostatic interactions. We then present the results of an extensive validation of these force fields, by means of molecular mechanics (MM) energy minimizations and molecular dynamics (MD) simulations of the crystal structures of these oligomers and polymers. While our "best" force field does not outperform the others on each of the investigated systems, it provides a balanced description of their overall structure and energetics. Finally, our MM minimizations and MD simulations confirm that the reported crystal structures of P3MBT and P3BT are stable and correspond to well-defined energetic minima. The room-temperature MD simulations reveal a certain degree of side-chain disorder, even in our virtually defect-free polymer crystal models.