The heats of formation for the molecules BH(3)PH(3), BH(2)PH(2), HBPH, AlH(3)NH(3), AlH(2)NH(2), HAlNH, AlH(3)PH(3), AlH(2)PH(2), HAlPH, AlH(4)(-), PH(3), PH(4), and PH(4)(+), as well as the diatomics BP, AlN, and AlP, have been calculated by using ab initio molecular orbital theory. The coupled cluster with single and double excitations and perturbative triples method (CCSD(T)) was employed for the total valence electronic energies. Correlation consistent basis sets were used, up through the augmented quadruple-zeta, to extrapolate to the complete basis set limit. Additional d core functions were used for Al and P. Core/valence, scalar relativistic, and spin-orbit corrections were included in an additive fashion to predict the atomization energies. Geometries were calculated at the CCSD(T) level up through at least aug-cc-pVTZ and frequencies were calculated at the CCSD(T)/aug-cc-pVDZ level. The heats of formation of the salts [BH(4)(-)][PH(4)(+)](s), [AlH(4)(-)][NH(4)(+)](s), and [AlH(4)(-)][PH(4)(+)](s) have been estimated by using an empirical expression for the lattice energy and the calculated heats of formation of the two component ions. The calculations show that both AlH(3)NH(3)(g) and [AlH(4)(-)][NH(4)(+)](s) can serve as good hydrogen storage systems that release H(2) in a slightly exothermic process. In addition, AlH(3)PH(3) and the salts [AlH(4)(-)][PH(4)(+)] and [BH(4)(-)][PH(4)(+)] have the potential to serve as H(2) storage systems. The hydride affinity of AlH(3) is calculated to be -70.4 kcal/mol at 298 K. The proton affinity of PH(3) is calculated to be 187.8 kcal/mol at 298 K in excellent agreement with the experimental value of 188 kcal/mol. PH(4) is calculated to be barely stable with respect to loss of a hydrogen to form PH(3).