Monte Carlo morphological modelling of a P3HT:PCBM bulk heterojunction organic solar cell

To deepen the understanding of morphology evolution in bulk heterojunction P3HT:PCBM organic photovoltaics system by thermal treatment, domain‐size‐dependent interfacial energies were first determined by coarse‐grained molecular dynamics modelling and then used in Monte Carlo simulations of the morphology evolution. Thereby initial conditions associated with optimal interfacial surface area, continuous volume, as well as domain sizes, and spatial distributions of the phase separated domains were identified. In line with earlier studies, a 1:1 P3HT:PCBM blend ratio is found to exhibit the most efficient morphology for exciton dissociation and charge transport. Our simulations reveal that preseeding of P3HT crystal at the anode side prior to the annealing process will be instrumental to pin the formation of P3HT at the favorable electrode especially when seeding exceeds a threshold of 10% surface coverage, whereas denser seeding patterns beyond the threshold did not improve the active layer morphology further. The observed trilayer depth profile (in the absence of preseeded P3HT crystals) implies that the commonly used thickness 100 nm of the active layer is not ideal for ensuring that donor and acceptor phases dominate at opposite ends of the active layer. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 270–279