Force fields for macromolecular assemblies containing Diketopyrrolopyrrole and Thiophene

Jiang, Ling, Rogers, David M., Hirst, Jonathan D. and Do, Hainam (2020) Force fields for macromolecular assemblies containing Diketopyrrolopyrrole and Thiophene. Journal of Chemical Theory and Computation, 16 (8). pp. 5150-5162. ISSN 1549-9618

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
Available under Licence Creative Commons Attribution.
Download (13MB) | Preview


Utilizing a force-matching procedure, we parametrize new force fields systematically for large conjugated systems. We model both conjugated polymers and molecular crystals that contain diketopyrrolopyrrole, thiophene, and thieno[3,2-b]thiophene units. These systems have recently been found to have low band gaps, which exhibit high efficiency for photovoltaic devices. The equilibrium structures, forces, and energies of the building block chromophores, diketopyrrolopyrrole, thiophene, and thieno[3,2-b]thiophene computed using our parameters are comparable to those computed using the reference electronic structure method. We assess the suitability of this new force field for electronic property calculations by comparing the electronic excitation properties computed along classical and ab initio molecular dynamics trajectories. For both trajectories, we find similar distributions of TDDFT-calculated excitation energies and oscillator strengths for the building block chromophore diketopyrrolopyrrole-thieno[3,2-b]thiophene. The structural, dynamic, and electronic properties of the macromolecular assemblies built upon these chromophores are characterized. For both polymers and molecular crystals, pronounced peaks around 0° or 180° are observed for the torsions between chromophores under ambient conditions. The high planarity in these systems can promote local ordering and π–π stacking, thereby potentially facilitating charge transport across these materials. For the model conducting polymers, we found that the fluctuations in the density of states per chain per monomer is negligibly small and does not vary significantly with chains comprising 20–40 monomers. Analysis of the electron–hole distributions and the transition density matrices indicates that the delocalized length is approximately 4–6 monomers, which is in good agreement with other theoretical and experimental studies of different conducting polymers. For the molecular crystals, our investigation of the characteristic time scale of the fluctuation in the excitonic couplings shows that a low-frequency vibration below 100 cm–1 is observed for the nearest neighbors. These observations are in line with previous studies on other molecular crystals, in which low-frequency vibrations are believed to be responsible for the large modulation of the excitonic coupling. Thus, our approach and the new force fields provide a direct route for studying the structure–property relations and the molecular level origins of the high efficiency of these classes of materials.

Item Type: Article
Schools/Departments: University of Nottingham Ningbo China > Faculty of Science and Engineering > Department of Chemical and Environmental Engineering
Identification Number:
Depositing User: Yu, Tiffany
Date Deposited: 11 Sep 2020 07:15
Last Modified: 11 Sep 2020 07:15

Actions (Archive Staff Only)

Edit View Edit View