In this paper we present molecular dynamics simulations of silicon nitride, both in bulk and as an interface to crystalline silicon. We investigate, in particular, the bonding structure of the silicon nitride and analyze the simulations to search for defective geometries which have been identified as potential charge carrier traps when silicon nitride forms an interface with silicon semiconductors. The simulations reveal how the bonding patterns in silicon nitride are dependent upon the stoichiometry of the system. Furthermore we demonstrate how having an "interphase", where the nitrogen content in silicon gradually reduces towards pure silicon across a boundary region, as opposed to an interface where there is an abrupt drop in nitrogen concentration at the boundary, can result in significantly different numbers of certain important carrier trap