"monomer" unit of which is shown in Fig. , forms alternating layers that may or may not have an offset, as the bulk materials are amorphous or polycrystalline and in some cases the micro-domains are sensitive to the presence or absence of solvent molecules occupying the pore. For applications in catalysis it is desirable that the pores be accessible, but while it has been suggested that most COF materials exhibit eclipsed-cofacial stacking (sometimes called "AA" stacking, borrowing the nomenclature used for graphene), detailed analysis of powder diffraction along with theoretical calculations suggests that slipped-or offset-stacking (sometimes called "inclined stacking") is common. The preference for eclipsed-cofacial ("AA") versus fully-staggered ("AB") stacking may be controllable via functionalization. In the present work, we examine stacking interactions in homo-dimers of the four monomers shown in Fig. , with the goal of elucidating the intermolecular forces that drive offset-stacking. The computational tool that we will employ is extended symmetryadapted perturbation theory with many-body dispersion (XSAPT+MBD), which was also used in some of the aforementioned analyses of offset-stacking in polycyclic aromatic hydrocarbons. For small noncovalent dimers where high-level calculations are feasible, XSAPT+MBD is comparable in accuracy (to within 1 kcal/mol) with the best ab initio benchmarks, ,34 yet for large systems its monomer-based nature makes it more affordable even than density functional theory (DFT). Crucially, XSAPT+MBD comes equipped with a physically-motivated energy decomposition analysis, which allows us to investigate the forces responsible for offset-stacking, and its energy decomposition is faithful component-by-component to highlevel ab initio benchmarks, again at the level of 1 kcal/mol.