Growth of two-dimensional covalent organic frameworks on substrates: insight from microsecond atomistic simulations

While growing two-dimensional covalent organic frameworks (2D COFs) on substrates holds promise for producing functional monolayers, the presence of many defects in the resulting crystals often hinders their practical applications. Achieving structural order while suppressing defect formation necessitates a detailed atomic-level understanding. The key lies in understanding the polymerization process with high nano-scale accuracy, which presents significant challenges. Here, we perform microsecond atomistic molecular dynamics simulations to describe the deposition and polymerization of cyclohexa-m-phenylene on metal substrates, closely mimicking experimental conditions. Our improved approach highlights that 2D polymerization occurs through monomer addition and island coalescence, with a pre-bonding stage allowing monomers/oligomers to dynamically adjust their configurations to the expanding island structures. Our results elucidate the mechanisms underlying the formation of vacancy and dislocation defects during 2D polymerization as well as their healing processes. Overall, our findings underscore the significant roles that high surface mobility, effective monomer-substrate anchoring, high framework rigidity, moderate monomer coordination, and low bonding rate play in forming large, extended 2D crystals while suppressing vacancy and dislocation defects. We demonstrate how these factors can be tuned through substrate selection, deposition rate modulation, and temperature control, thereby offering valuable insight for strategically optimizing on-surface 2D polymerizations. Microsecond atomistic simulations reveal the complex interplay among monomer mobility and coordination, monomer-substrate anchoring, framework rigidity, and bonding rates in polymerization on metal surfaces.