Polyphosphazenes, because of their unique properties, have generated many opportunities to explore a variety of applications. These applications include areas such as biomedical research (e.g. drug delivery) and material science (e.g. fire-resistant polymers). Phosphazenes potentially have more variations then benzene analogues because of different substitution patterns. Here we present A computational study of the chemical modifications to a group of cyclic phosphazenes mainly hexachlorophosphazene (PNCl2)3. This study focuses on the relative energies of reactivity of hexachlorophosphazene to understand their geometry and the complexes they likely form. We compare diols, amino alcohols, and diamines with a carbon linker of 1-7 atoms. These heteroatom chains are attached to a single phosphorus atom or adjoining phosphorus atoms to form ring structures of geminal, vicinal (cis), and vicinal (trans) moieties. We find that the reactivities of “heteroatom caps” are predicted to be O,O (diol) > N,O (amino alcohol) > N,N (diamine). These results can be used to predict energetics and thus the stability of new compounds for biomedical and industrial applications.