Computational Study of Intramolecular Heterocyclic Ring Formation with Cyclic Phosphazenes

Computational Study of Intramolecular Heterocyclic Ring Formation with Cyclic Phosphazenes

Whelton A. Miller III

Department of Bioengineering School of Engineering and Applied Science

University of Pennsylvania Philadelphia, PA 19104 USA

Preston B. Moore Department of Chemistry and Biochemistry University of the Sciences in Philadelphia

Philadelphia, PA 19104 USA

AbstractPolyphosphazenes, 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.

KeywordsPhosphazenes; Cyclic Phosphazenes; Ring Formation; Quantum Mechanics (QM); Density Functional Theory (DFT); Phosphorus Nuclear Magnetic Resonance (31P NMR); Becke Lee, Yang, and Parr hybride Density Functional Theory (B3LYP)

1. INTRODUCTION

   This Polyphosphazenes have been used in a variety of application, including several areas of biomedical research, such as drug delivery systems [1,2], and in material science, such as fire-resistant polymers [3]. Due to their substituition patterns, phosphazenes potenially have more variations than benzene analogs. Here we present a computational study of the chemical modifications to a group of cyclic phosphazenes (figure 1). We report the general conformation of the compounds, as well as relative energies, and potential applications.

   In the past, phosphazenes have been studied for their properties both as linear polymers [4,5,6], and cyclic structures [7,8,9]. Their general electronic structure has also been studied [10,11,12], focusing on the interactions between the nitrogen and phosphorus atoms linear or cyclic geometries. Intramolecular reactions of difunctional nucleophiles, amino-alcohols, diols, and diamines with non-

   substituted cyclic phosphazenes such as hexachlorophosphazene (NPCl2)3, can yield three types of products [13-16]. First, the dinucleophile may replace two chlorine atoms at the same phosphorus atom to form geminal or spirocyclic structures [17-21]. Second, chlorine units on adjacent phosphorus atoms may be replaced to produce a vicinal substituted transannular-bridged cis-structure [20,21]. Third, substitution can produce a trans-vicinal bridge (figure 1) [13,20].

   Fig 1. Substitution patterns for cyclized phosphazenes. X and Y are heteroatoms, i.e., N or O.

   In this paper we focus on the relative energies of hexachlorophosphazenes with heteroatom chains, specifically the energy difference between geminal and vicinal structures. By using Quantum Mechanics (QM) to calculate the energy of formation of the resulting geminal and vicinal substituted structures, we strive to understand the relationship between heteroatom chain length, the resulting cyclized compound, and the likelihood of compound formation.

   Phosphorus Nuclear Magnetic Resonance (31P NMR) spectroscopy is widely used to investigate the structure of

   phosphorus compounds. We use QM calculations to predict the 31P NMR shifts of our structures, and therefore correlate structural changes to changes in shielding. We analyzed our results, including the 31P NMR predictions, to further understand the relative energies of dinucleophiles with cyclic phosphazenes. Our analysis allows the prediction of reactivity and the rational design of future reactions and materials based on phosphazene chemistry.

2. METHODOLOGY

   We examined three distinct conformations, consisting of geminal (spirocyclic) and vicinal (transannular) bridged compounds (cis or trans) (figure 1). We attached heteroatom chains 1-9 atoms in length, including heteroatom caps, but excluding hydrogen (HaX-(CH2)n-YHb), where X and Y can consist of any combination of O or N. Conformational analysis on the expected products was performed with Molecular Operating Environment (MOE) [22]. The geometry was optimized using Hartree Fock (HF), and Becke Lee, Yang, and Parr hybrid Density Functional Theory (DFT) B3LYP with basis set 6-31g(d,p), and 6-31+g(d,p) [23-25]. Comparative calculations were done under tight optimization parameters, with diffusion functions added to the calculation to account for the possible interaction of the lone pair electrons on N, O, P, and Cl atoms. Calculations were carried out using Gaussian03 [26]. The lowest energy structures after optimization are reported and used for analysis. We then used equation 1 to determine the reaction energy (E), allowing for discrimination between favorable and unfavorable reactions (scheme 1). 31P NMR calculations were performed with the GIAO [26,27] method from the Gaussian03 package, using DFT with basis set 6-31g(d,p).

   1. E(final)=E(PhosphazeneFP) + E(2HCl) (1)

   2. E(initial)=E(Phosphazene ring (P3N3Cl6)+Amino Alcohol) + E(HaX-(CH2)n-YHb)

   3. E(reaction)=E(final) – E(initial)

   Scheme 1. Reaction schemes for cyclic phosphazene with Heteroatom chains. X and Y are heteroatom caps, i.e., N or O for geminal substitution, vicinal (cis) and vicinal (trans) substitution.

3. RESULTS AND DISCUSSION

   When compared to the cyclic phosphazenes, carbon analogs such as benzene are limited because they permit only single substitution at a given position on the ring. Due to the substitution pattern of phosphazenes, branching can be significantly different than in carbon analogs. One such substitution pattern is geminal substitution. The geminal spirocyclic structures form a bicyclic ring at the phosphorous. Vicinal structures consist of attachments at two adjacent phosphorous atoms that can be either cis where the attachment is on the same side of the ring, or trans where the attachment is on different sides of the ring. In our investigations we will use a carbon chain capped with heteroatoms, where X and Y can be either N or O in any substitution pattern to react with the phosphazene (scheme 1). The result of a secondary intra-atomic reaction made after the initial heteroatom bond can give either the geminal or vicinal (cis or trans) substitution pattern.

   1. Geminal Cyclic Phosphazene Complexes

      Geminal structures (figure 2) can be precursors to larger phosphazene macromolecules. These compounds in turn can be incorporated as dendrimers and linear polymers for numerous uses. Understanding the energetics of these systems, therefore, is crucial to understanding the structural properties, preferred products of reactions, and by extension, the properties and reactivity of larger phosphazene polymers.

      Fig 2. Geminal Substituted Cyclized Phosphazene. X and Y can equal N or O.

      Fig 3. Plot of energy (E) verses chain length for geminal cyclic phosphazenes. X-ais represents HaX-(CH2)n-YHb, where X and Y can be O or N. Each product has been calculated using three different levels of theory.

      Figure 3 shows the energy of formation of geminal substituted phosphazenes. For each product we used three different levels of theory to assess basis set dependence, and observed that the trend is independent of the theory used. Specifically, the energy of formation is highly dependent on the ability to form cyclic rings and is lowest for the 5-atom (3- carbon chain capped with heteroatoms) chain that forms a six- membered ring, completed by the phosphorus atom. This preference parallels the well-known strain of ring formation in carbon analogs such as cyclohexane. The dip in energy at the 5-atom chain length is similar to the dip in energy seen in the formation of cyclohexane, when compared to cyclobutane, cyclopentane, cycloheptane, etc [27]. This can be seen in the energy difference between the 5-atom, 6-atom, and 7-atom chain lengths and the rest of the complexes. We find that for structures resulting in geminal substitution, the energy is almost solely dependent on the size and geometry of the resulting geminal ring with little influence from the phosphazene ring.

   2. Vicinal Bridged Cyclic Phosphazene Complexes Reactions of difunctional nucleophiles with

      cyclophosphazenes to yield vicinal substitution patterns have

      been of considerable interest, even though they are found to be thermodynamically less stable compared to geminal or spiro-substituted cyclophosphazenes [29]. Vicinally substituted structures create strain in the phosphazene ring for example, as the ring tries to bend the P-N-P angle to accommodate a bond between the heteroatom chain and the adjacent phosphorus atom. The bend in the P-N-P angle is less pronounced in cis, which depends on the heteroatom chains ability to reach across the face of the cyclophosphazene, than in trans where the heteroatom chain crosses the back face of the cyclophosphazene to reach the adjacent phosphorus.

      Vicinal Bridged (Cis)

      Fig 4. Cis Vicinal Substituted Cyclized Phosphazene, X and Y can equal N, or O. View of optimized geometries for Vicinal (cis) phosphazene complexes.

      Fig. 5 Plot of energy (E) verses chain length for cis Vicinal substituted phosphazene. X-axis represents HaX-(CH2)n-YHb, where X and Y can be O or N. Each product was calculated using three different levels of theory.

      As mentioned earlier, vicinal structures substituted in a cis (Figure 4) manner do not have a large additional ring strain arising from the carbon chain attempting to traverse cyclic phosphazenes. Without this bending at the P-N-P angle, the structures formed relate to chain length as expected. Favorable reactions begin at chain lengths of 4 atoms, depending on the heteroatom caps, with oxygen being more reactive (Figure 5). Again, we have used three different levels of theory to assess the basis set dependence on E, and as before, the trends are independent of the basis set used. Chain lengths of 4 atoms or more are favored for any substitution pattern, diamines, diols, or amino-alcohols (i.e., E)is < 0). There is about a ~6-8 kcal/mol energy difference between the 7-atom chain length and the 8-atom chain length. As shown in Figure 5 the optimum chain length for vicinal (cis) substitution is 8 atoms long, which also happens to be the crossover point for vicinal (cis) and geminal substitution; both patterns have the same E at this chain length (Figures 3 & 5). Although Figure 5 shows a consistent decrease in energy with an increase in chain length for 3-8 atoms, the 9- atom chain is slightly less favorable than the 8-atom chain. Initially, one might have expected the 9-atom chain to be more favorable because its greater number of degrees of freedom can accommodate ring constraints, but chain lengths larger than 8-atoms cause steric crowding at ring closure points, which accounts for the increase in energy for these lengths. We conclude that there are no geometries that will accommodate the constraints of the ring with an energy lower than that observed for 8-atoms chains.

      Vicinal Bridged (Trans)

      Fig. 6 Trans Vicinal Substituted Cyclized Phosphazene. X and Y can equal N, or O.

      With respect to the geminal and vicinal (cis) compounds, vicinal (trans) substitution (Figure 6) is less dependent on the chain length or the ability to cyclize, but more dependent on the ring strain of the phosphazene (Figure 7). Ring strain breaks the planarity of the cyclic phosphazene, and decreases the pseudo-resonance, which is distinct from full- resonance in benzene [11]. The distortion of the phosphazene ring with respect to planarity adds a significant amount of energy to the formation energy of the resulting vicinal (trans) substituted phosphazene [30], which dominates the reaction energy of the heteroatom chain with the phosphazene. As expected, increased chain length reduces distortion, which favors reaction (Figure 7).

      Fig. 7 Plot of energy (E) verses chain length for vicinal (trans) substituted phosphazene. X-axis represents HaX-(CH2)n-YHb, where X and Y

      can be O or N. Each product was calculated using three different levels of

      theory.

      Fig 8 Plot of energy (E) versus chain length for geminal (8a), vicinal (cis)(8b), vicinal (trans)(8c), and vicinal (trans, 6-9 atom chain length). X- axis represents HaX-(CH2)n-YHb, where X and Y can be O or N. Calculated using (HF/6-31+G(dp)).

      A comparison of Figure 5 to Figure 7 shows that cis substitution is favored over trans substitution in most cases. The energy difference between trans and cis vicinal substituted phosphazenes includes contributions from ring strain, and the ability of the chain to reach around the ring. We did find, however, that oxygen capped chains, i.e. diols and amino-alcohols, were more favored than diamines in vicinal (cis) and vicinal (trans) substitution.

   3. Comparison of Substitution Patterns

   As expected, in all cases oxygen serves as a stronger reactant than Nitrogen. Figure 8 shows that diols are more favorable, followed by the mixed heteroatom chains amino- alcohols, then diamines. Depending upon chain length, the second displacement will either add geminally, or trans or cis vicinally, with the geminal reactions being far more favorable.

   Geminal addition is highly dependent upon the ability to form an intra-molecular ring, wth the stability of this ring being dependent upon its size. A conventional pattern of stability such as that seen in classic organic chemistry with respect to cycloalkane rings may be expected. For chain lengths 7 in all cases, diols, amino-alcohols, and diamines favor geminal substitution (Figure 9), with 5- and 6- membered rings being most stable. These rings are further stabilized by chair or boat conformation, with chair being prefered for six-membered geminal structures. Reactivity follows the previously observed trend (Figure 8) of diols first, followed by amino-alcohols and then diamines (Figures 9,10). Intra-atom reactions, which produce geminal substitution, are more favorable when compared to all vicinal substitutions at lower chain lengths.

   Once the chain length increases beyond 5- and 6- membered rings, a cross-over of preference occurs due to steric interactions (Figure 9). These steric interactions are relieved by reacting with the adjacent (vicinal) phosphorous. Competition between cis and trans vicinal substitution is also dependent on steric interactions, although the distortion of the

   Fig. 9 Plot of energy (E) verses chain length for comparison of geminal versus vicinal (cis) substituted phosphazene, amino-alcohols (9a), diols (9b), and diamines (9c). X-axis represents HaX-(CH2)n-YHb, where X and Y can be O or N. Calculated using (HF/6-31+G(dp)).

   phosphazene ring plays a more important role (Figure 10). As the chain length increases and the ring strain decreases, the vicinal (trans) product becomes comparable in reactivity to the geminal and vicinal (is) substituted products.

   Fig 10. Plot of energy (E) verses chain length and reaction energy for geminal versus vicinal (cis and trans) substituted phosphazenes, amino- alcohols (10a), diols (10b), and diamines (10c). X-axis represents HaX- (CH2)n-YHb, where X and Y can be O or N. Calculated using (HF/6- 31+G(dp)).

   31P NMR of cyclic Phosphazene structures

   31P NMR has been widely used for the analysis of cyclophosphazenes and their analogs. Electronegative atoms

   can influence the phosphazene skeleton due to charge density at the neighboring phosphorus atoms [30-32]. Therefore, 31P NMR data can provide structural information while probing the nature of the electronic interaction within the phosphazene ring and any substituents. Electron withdrawal

   from the substituent-phosphorus bond can lead to either shielding or deshielding of the phosphorus atom. Substituent effects on phosphorus chemical shifts will depend on electronegativity and the bond angle of the substituent attached [31]. QM calculations determine the difference in shielding, which distinguishes between compounds substituted geminally or vicinally. Shielding can then be adjusted using further QM calculations with respect to a reference (330ppm, hexachlorophosphazene), or further adjusted from experimental data (21.28ppm, hexachlorophosphazene in CDCl3) for more accurate calculated spectra. Using this method we determined the degeneracy, relative chemical shifts, and neglect splitting patterns of these molecules. For geminal substitution, we expected that non-substituted phosphorus atoms (PCl2) would have similar chemical shifts, and this hypothesis was confirmed by the calculated 31P NMR. Because substitution only takes place on a single phosphorus atom, the chemical shift should only be reflected at that atom. Table 1 shows changes in shielding only in the PR2 phosphorus atom. We found similar chemical shifts in the two remaining PCl2 phosphorus atoms for geminally substituted diols, diamines, and amino-alcohols.

   Tab. 1 Calculated isotropic 31P NMR data using hexachlorophosphazene as a reference (351.28ppm), for geminal substituted amino alcohols, diamines, and diols by chain length. Calculated using DFT (B3LYP/6-31G(dp)).

   Geminal 31P data

   Fig 11. Calculated Isotropic 31P NMR spectra using hexachlorophosphazene as a reference (351.28ppm), for geminal substituted amino alcohol with 5- atom heteroatom chain. Calculated using DFT (B3LYP/6-31G(dp)).

   For vicinal (cis) molecules, the calculated 31P NMR shows similar chemical shifts for phosphorus atoms of like substitution (PClR). Becuase there is little or no disturbance to the planarity of the phosphazene ring system, distinction between phosphorus of like substitution is minimal. As expected, the two PClR phosphorus atoms showed approximately the same chemical shifts, while PCl phosphorus atoms had a significant difference in shielding.

   Vicinal (cis) 31P data

   Fig. 12 Calculated Isotropic 31P NMR spectra using hexachlorophosphazene as a reference (351.28ppm) for vicinal (cis) substituted amino alcohol with 5- atom heteroatom chain. Calculated using DFT (B3LYP/6-31G(dp)).

   Tab. 2 Calculated isotropic 31P NMR data using hexachlorophosphazene in CDCl3 as a reference (351.28ppm), for vicinal (cis) Substituted amino alcohols, diamines, and diols by chain length. Calculated using DFT (B3LYP/6-31G(dp)).

   For vicinal (trans) substitution, the expected result is similar chemical shifts at two of the three phosphorus atoms, similar to vicinal (cis). However, for compounds bridged with 3-5 atoms, the chemical shifts are different for all three phosphorus nuclei, which can be attributed to the ring distortion mentioned above. This distortion breaks the planarity of the cyclophosphorus ring and therefore disrupts the rings pseudo-resonance. With this resonance disturbed, charge can become more localized, and therefore will affect the shielding of the local phosphorus atom. The vicinal (trans) phosphazene analogs affected by this distortion are in shown in italics (Table 3). This distortion effect occurs in products of short heteroatom chains; as chain length increases, distortion decreases similarly to the trend seen in Figure 7. Again, these observations indicate that ring distortion/loss of planarity of the phosphazene ring becomes less apparent as the chain becomes longer. Ring strain is no longer significant in chain lengths of 6 or more atoms.

   Vicinal (trans) 31P data

   Fig. 13 Calculated Isotropic 31P NMR spectra using hexachlorophosphazene as a reference (351.28ppm), for vicinal (trans) substituted amino alcohol with 5-atom heteroatom chain. Calculated using DFT (B3LYP/6-31G(dp)).

   Tab. 3 Calculated isotropic 31P NMR data using hexachlorophosphazene in CDCl3 as a reference (351.28ppm), for vicinal (trans) substituted amino alcohols, diamines, and diols by chain length. Calculated using DFT (B3LYP/6-31G(dp)). Those phosphorus atoms affected by ring puckering are shown in italics and bold type.

4. CONCLUSION

The factors that influence the reactivity pattern of cyclic phosphazenes include the reactivity of a strong nucleophile,

i.e. the basicity of the attaching reactant, and the steric influence of reactants and substituents already attached. The main influence highlighted by this work is steric hindrance by a bulky substituent. Based on the size and nature of the substituent, we can distinguish the different types of products expected.

Reactivity = O-(CH2)n-O > N-(CH2)n-O > N-(CH2)n-N (2)

Geminal > Vicinal Cis > Vicinal Trans up to n=7. At n= 8 vicinal cis

> geminal (3)

For diols, diamines, and amino alcohols, reactivity favors the oxygen-containing nucleophiles, diols, and amino alcohols over diamines (Equation 2, and Figure 8). This would be expected, because oxygen is a more electronegative nucleophile. Geminal and vicinal (cis) substitution is favored until chain length falls below 7

atoms. At chain lengths greater than 7, vicinal (cis) outcompetes geminal substitution (Equation 3).

Geminal substitution is dependent on its ability to form 5- or 6-membered rings, similar to carbon analogs such as cyclohexane. However, vicinal substitution is dependent on steric factors and the overall geometry of the ring. If the planarity of the phosphazene ring is distorted as seen in vicinal (trans) substitution, the ring distortion becomes the primary energy parameter. To avoid this ring distortion, the dinucleophillic chain must be able to reach the adjacent phosphorus atom. We found that chain lengths greater than 7 atoms have the best chance of avoiding this ring puckering phenomenon.

The reliance on longer chain lengths and the larger amount of steric interactions involved in vicinal substitution is why geminal substitution is favored for chain lengths greater than 7 atoms. Cis or trans substitution in vicinal compounds is, again, dictated by sterics. Steric hindrance from ring distortion versus the ability to reach across the face of the phosphazene distinguishes between vicinal (cis) and vicinal (trans) substitution. At chain lengths greater than 7, this hindrance becomes less of a factor, and trans substitution becomes competitive.

We have shown that 31P NMR can be used to track the substitution patterns and the differences between geminal, vicinal (cis), and vicinal (trans) substitution. Geminal substitution distinguishes a change only on a single phosphorus atom, which corresponds to substitution on a single phosphorus atom. For vicinal (trans), we found similar chemical shifts between the substituted phosphorus atoms at chain lengths n 4. At chain lengths n 3, we observed differences in shielding of all phosphorus nuclei due to a loss of resonance resulting from ring distortion. Vicinal (cis) shows similar chemical shifts for similarly substituted phosphorus atoms, as expected, and the ring distortion effect is negligible in these compounds. Changing the basis set does show a change in relative energy however, the pattern of relative energy remains the same. In summary, we found that different properties can be created by varying heteroatom chain length. Using this information, we hope to contribute to future design of larger phosphazene polymers and other biomedical compounds. We believe that better understanding of the substitution patterns is critical to better design of these compounds.

ACKNOWLEDGEMENT

This work has been partially supported by the National Science Foundation under Grant No NSF #CHE 1229564, the University of the Sciences in Philadelphia, the H. O. West Foundation, and the West Center for Computational Chemistry and Drug Design. The authors would like to thank the researchers from the West Center for their continued helpful discussions, in particular Dr. Zhiwei Liu for her insightful and helpful comments.

REFERENCES

1. S. A. Ibim, A. A. Ambrosio, D. Larrier, H. R. Allcock, C. T. Laurencin, Controlled macromolecule release from poly(phosphazene) matrices J. Controlled Release. 1996, 40, pp. 31-39.

2. S. Y. Cho, H. R. Allcock, Dendrimers Derived from PolyphosphazenePoly(propyleneimine) Systems: Encapsulation and Triggered Release of Hydrophobic Guest Molecules Macromolecules 2007, 40, pp. 3115-3121.

3. J. Y. Allen, H. R. Allcock, The Use of Phosphazenes as Fire Resistant Materials J. Fire Sciences 1993, 11, pp 320-328.

4. H. R Allcock, S. Kwon, An ionically crosslinkable polyphosphazene: poly[bis(carboxylatophenoxy)phosphazene] and its hydrogels and membranes Macromolecules 1989, 22, pp. 75-79.

5. H. R. Allcock, M. L. Turner, Ring expansion and polymerization of transannular bridged cyclotriphosphazenes and their spirocyclic analogs Macromolecules 1993, 26, pp. 3-10.

6. H. R. Allcock, G. K. Dudley, Lower Critical Solubility Temperature Study of Alkyl Ether Based Polyphosphazenes Macromolecules 1996, 29, pp. 1313-1319.

7. D. B. Davies, T. A. Clayton, R. E. Eaton, R.A. Shaw, A. Egan, M. B. Hursthouse, G. D. Sykara, I. Porwolik-Czomperlik, M. Siwy, K. Brandt, Chiral Configurations of Cyclophosphazenes J. Am. Chem. Soc. 2000, 122, pp. 12447-12457.

8. P. I. Richards, , A. Steiner Cyclophosphazenes as Nodal Ligands in Coordination Polymers Inorg. Chem. 2004, 43, pp. 2810-2817.

9. Allcock, H. R.; Kellam III, E. C. Incorporation of Cyclic Phosphazene Trimers into Saturated and Unsaturated Ethylene-like Polymer Backbones Macromolecules 2002, 35, pp. 40-47.

10. M. Breza, The electronic structure of planar phosphazene rings

    Polyhedron. 2000, 19, 389-397.

11. A. B. Chaplin, J. A. Harrison, P. J. Dyson, Revisiting the Electronic Structure of Phosphazenes Inorg. Chem. 2005, 44, pp. 8407-8417.

12. R. C. Haddon Chem. Phys. Letters. 1985, pp. 120, 372-374.

13. H. R. Allcock, M. L. Turner, K.B. Visscher, Synthesis of transannular- and spiro-substituted cyclotriphosphazenes: x-ray crystal structures of 1,1-[N3P3(OCH2CF3)4{O2C12H8}],1,3-[N3P3(OCH2CF3)4{O2C12H8}], 1,1- [N3P3(OCH2CF3)4{O2C10H6}], and 1,3-[N3P3(OCH2CF3)4}O2C10H6}] Inorg. Chem. 1992, 31, pp. 4354-4364.

14. H. R. Allcock, Recent advances in phosphazene (phosphonitrilic) chemistry Chem. Rev. 1972, 72(4), 315-356.

15. H. R. Allcock, J. S. Rutt, M. Parvez, Synthesis of cyclic phosphazenes with isothiocyanato, thiourethane, and thiourea side groups: x-ray crystal structure of N3P3(NMe2)3(NCS)3 Inorg. Chem. 1991, 30, pp. 1776-1782.

16. C. W. Allen, Regio- and stereochemical control in substitution reactions of cyclophosphazenes Chem. Rev. 1991, 91, 119-135.

17. Iter, et. al. PhosphorusNitrogen Compounds. 14. Synthesis, Stereogenism, and Structural Investigations of Novel N/O Spirocyclic Phosphazene Derivatives Inorg. Chem. 2007, 46, 9931-9944.

18. K. Muralidharan, P.; Venugopalan, A. J. Elias, Ansa versus Spiro Substitution of Cyclophosphazenes: Is Fluorination Essential for Ansa

    to Spiro Transformation of Cyclophosphazenes? Inorg. Chem. 2003,

    42, 3176-3182.

19. G. A. Gabino, Synthesis of New Phosphazene High Molecular Weight Polymers Containing Functionalized and Optically Active Spirocyclic Groups Macromolecules 1998, 31, pp. 3189-3196.

20. A. J. Elias, Twamley, B.; Shreeve, J. M. Syntheses and Experimental Studies on the Relative Stabilities of Spiro, Ansa, and Bridged Derivatives of Cyclic Tetrameric Fluorophosphazene Inorg. Chem. 2001, 40, pp. 2120-2126.

21. Molecular Operating Environment (MOE), 2011.10; Chemical Computing Group Inc., 1010 Sherbooke St. West, Suite #910, Montreal, QC, Canada, H3A 2R7, 2011.

22. A. D. Becke, A new mixing of HartreeFock and local densityfunctional theories J. Chem. Phys. 1993, 98, 1372.

23. C. T. Lee, W. T. Yang, R. G. Parr, Development of the Colle-Salvetti conelation energy formula into a functional of the electron density Physical Review B 1988, 37, pp. 785-789.

24. T. H. Dunning, J.; Huzinaga, P. J. In Modern Theoretical Chemistry; Plenum: New York, 1976;

25. Frisch, et al.; Gaussian03; Gaussian, Inc.: Wallingford CT, 2004.

26. E. D. Jemmis, Kiran, B. Aromaticity in X Y H (X = B, Al, Ga; Y = N, P, As), X Z H (Z = O, S, Se), and Phosp3ha3ze6nes. Theoretical Study of the Struc3tur3es3, Energetics, and Magnetic Properties Inorg. Chem 1998, 37, pp. 2110-2116.

27. McMurray, J. Organic Chemsitry; Brooks/Cole: Pacific Grove CA, 2000.

28. K. Muralidharan, A. Elias, Preparation of the First Examples of AnsaSpiro Substituted Fluorophosphazenes and Their Structural Studies: Analysis of CH···FP Weak Interactions in Substituted Fluorophosphazenes J. Inorg. Chem. 2003, 42, pp. 7535-7543.

29. S. S. Krishnamurthy, M. Woods, Nuclear Magnetic Resonance of Cyclophosphazenes; Academic Press Inc.: London, 1987;

30. H. R. Allcock, P. J. Harris, R. A. Nissan, Organometallic phosphazenes: synthesis and rearrangement of propynyl- and propadienylcyclotriphosphazenes J. Am. Chem. Soc. 1981, 103, pp. 2256-2261.

31. M. A. Benson, S. Zacchini, R. Boomishanker, A. Steiner, Alkylation and acylation of cyclotriphosphazenes Inorg. Chem. 2007, 46, pp. 7097-7108.

SYNOPSIS

Phosphazenes have potentially more variations than benzene analogs because of their unique substitution patterns. This study focuses on the properties of phosphazenes, in particular hexachlorophosphazene, to understand their geometry and the relationship of preferred structure to reactivity. To compare the interactions, we tethered diols, amino alcohols, and diamines with a carbon linker of 1-7 atoms. These were attached to form ring structures of geminal, vicinal-cis, and vicinal-trans moities.