The intermolecular interactions underlying self-assembly are important for designing supramolecular architectures with defined function and composition over various length scales. Here, an integrated quantum mechanical/molecular dynamics (QM/MD) approach is introduced to explore the intermolecular interactions governing the two-dimensional selfassembly of a shape-persistent macrocycle, that is, tricarbazolo triazolophane (tricarb), on highly ordered pyrolytic graphite (HOPG). Three structural motifs (row, zig-zag, and rosette) for the self-association and growth of tricarb oligomers are analyzed using dispersion-corrected density functional theory to reveal their relative stability. QM analysis shows that the major factor in the self-assembly, via H-bonding contributions to the interaction energies, originates from favorable local dipole orientations between triazole units in the zig-zag and rosette motifs. Closing of the macrocyclic loop leads to the formation of a tertiary structure, which induces stability for the rosette formation relative to other motifs. MD simulations of these motifs on HOPG demonstrate the stability of the rosette, corroborating the QM results. In addition, the effect of long alkyl chains on the selfassembly dynamics is explored by simulating preexisting honeycomb domains of tricarbs with various alkyl chain lengths, without and with the presence of an explicit solvent. MD simulations reveal that alkyl chains mediate the formation of a honeycomb pattern, as observed experimentally.