Beryllium subphthalocyanines have been recently shown to be suitable candidates for photochemical devices if combined with appropriate donor systems. The ability of beryllium subphthalocyanines to self-assemble is explored for the first time by means of density functional theory calculations. Free dimers of beryllium subphtalocyanine and their corresponding complexes with water and pyridine are computed at the wB97X-D/6-311+G(d,p) level of theory. In contrast with the behavior reported for beryllium phthalocyanines, for beryllium subphtalocyanines, beryllium–aza–nitrogen intermolecular interactions are observed, suggesting that these species are likely to self-assemble. Aggregates of related structures such as beryllium subporphyrazines with axial groups confirm the importance of hydrogen bonds in the stacking.