Self-assembling of zinc phthalocyanines on ZnO (1010) surface through multiple time scales

We adopt a hierarchic combination of theoretical methods to study the assembling of zinc phthalocyanines (ZnPcs) on a ZnO (1010) surface through multiple time scales. Atomistic simulations, such as model potential molecular dynamics and metadynamics, are used to study the energetics and short time evolution (up to ∼100 ns) of small ZnPc aggregates. The stability and the lifetime of large clusters is then studied by means of an atomistically informed coarse-grained model using classical molecular dynamics. Finally, the macroscopic time scale clustering phenomenon is studied by Metropolis Monte Carlo algorithms as a function of temperature and surface coverage. We provide evidence that at room temperature the aggregation is likely to occur at sufficiently high coverage, and we characterize the nature, morphology, and lifetime of ZnPc's clusters. We identify the molecular stripes oriented along [010] crystallographic directions as the most energetically stable aggregates.