We performed nonadiabatic, surface hopping molecular dynamics simulations for H-type azobenzene tetramers to investigate exciton localization and transfer in these aggregates upon ππ * excitation. Free and constrained (SAM-like) models were considered, with short (a = 3.5 Å) and longer (a = 5.5 Å) intermolecular, nearest-neighbour distances. The electronic structure of the tetramers were treated at the semiempirical rAM1/FOMO-CIS level, allowing us to afford surface hopping dynamics simulations for the relatively large QM part including four azobenzene molecules. The exciton dynamics were explored using transition density matrix analysis.