Computations on all the possible positional isomers of the closo-azaboranes NBn-1Hn (n = 5-12) reveal substantial differences in the relative energies. Data at the B3LYP/6-311+G** level of density functional theory (DFT) agree well with expectations based on the topological charge stabilization, with the qualitative connectivity preferences of Williams, and with the Jemmis-Schleyer six interstitial electron rules. The energetic relationship involving each of the most stable positional isomers, 1-NB4H5, NB5H6, 2-NB6H7, 1-NB7H8, 4-NB8H9, 1-NB9H10, 2-NB11H11, NB11H12, was based on the energies (DeltaH) of the model reaction: NBH2 + (n-1)BHincrement --> NBnHn+1 (n = 4-11). This evaluation shows that the stabilities of closo-azaboranes NBn-1Hn (n = 5-12) increase with increasing cluster size from 5 to 12 vertexes. The "three-dimensional aromaticity" of these closo-azaboranes NBn-1Hn (n = 5-12) is demonstrated by their the nucleus-independent chemical shifts (NICS) and their magnetic susceptibilities (chi), which match one another well. However, there is no direct relationship between these magnetic properties and the relative stabilities of the positional isomers of each cluster. As expected, other energy contributions such as topological charge stabilization and connectivity can be equally important.