Investigating the Stacking Interactions Responsible for Driving 3D Self-Association of Tricarb Macrocycles

We have investigated the noncovalent forces that play a crucial role in the three-dimensional (3D) self-association of the tricarb macrocycle (composed of alternating triazoles and carbazoles) to understand the multilayer stacks observed through electron microscopy. To explore this idea quantitatively, we have investigated a stacked dimer model of tricarb, where we consider homochiral as well as heterochiral forms of the dimer. We have computed the rotational potential energy surface of the dimer by conducting an angle-dependent scan between the two macrocycles using different levels of theory including the RI-MP2 ab initio method. We observe that dimers oriented at an angle of 60° exhibit the highest stability, while a secondary minimum is observed at an angle of 30°. While density functional theory (DFT) describes the behavior of both minima very close to that obtained with RI-MP2, semiempirical and MM models appear to obtain only a shoulder instead of the second minimum. To further understand the underlying interactions responsible for stabilizing the self-assembly of the macrocycles, we employed energy decomposition analysis (EDA) using SAPT0. This quantitative assessment allowed us to identify the major contributing noncovalent interactions, including electrostatic, exchange-repulsion, dispersion, and induction. Finally, we expanded our study to evaluate the accuracy of the MIM (molecules-in-molecules) fragmentation methodology in capturing the crucial π-stacking interactions. Our benchmarking results using the MIM method accurately replicated the angle-dependent PES results. This shows the efficacy of MIM in predicting the noncovalent interactions responsible for the construction of 3D and other higher-order nanoarchitectures for tricarb and related compounds.